Skip to content
Snippets Groups Projects
Commit 70e65c79 authored by Paul Best's avatar Paul Best
Browse files

detection 1min

parent 19384e04
Branches
No related tags found
No related merge requests found
...@@ -6,8 +6,8 @@ parser = argparse.ArgumentParser() ...@@ -6,8 +6,8 @@ parser = argparse.ArgumentParser()
parser.add_argument('filename') parser.add_argument('filename')
args = parser.parse_args() args = parser.parse_args()
RORQUAL_HOPSIZE, RORQUAL_WINSIZE, RORQUAL_FE, RORQUAL_SAMPLEPERSAMPLE = 256, 4096, 4000, 25600 RORQUAL_HOPSIZE, RORQUAL_WINSIZE, RORQUAL_FE, RORQUAL_SAMPLEPERSAMPLE, RORQUAL_LENSIG = 256, 4096, 4000, 25600, 60
CACHA_HOPSIZE, CACHA_WINSIZE, CACHA_FE, CACHA_SAMPLEPERSAMPLE = 256, 512, 64000, 12800 CACHA_HOPSIZE, CACHA_WINSIZE, CACHA_FE, CACHA_SAMPLEPERSAMPLE, CACHA_LENSIG = 256, 512, 64000, 12800, 60
CACHA_RAPPORT_FE, RORQUAL_RAPPORT_FE, RAPPORT_N_CHANNELS = 128_000, 128_00, 5 CACHA_RAPPORT_FE, RORQUAL_RAPPORT_FE, RAPPORT_N_CHANNELS = 128_000, 128_00, 5
CACHA_NUM_DETEC, RORQUAL_NUM_DETEC = 10, 3 CACHA_NUM_DETEC, RORQUAL_NUM_DETEC = 10, 3
rorqual_predPos = lambda p: int(((p * RORQUAL_HOPSIZE + (RORQUAL_WINSIZE/2 + RORQUAL_HOPSIZE * 13) / 2) / RORQUAL_FE - 1) * RORQUAL_RAPPORT_FE) rorqual_predPos = lambda p: int(((p * RORQUAL_HOPSIZE + (RORQUAL_WINSIZE/2 + RORQUAL_HOPSIZE * 13) / 2) / RORQUAL_FE - 1) * RORQUAL_RAPPORT_FE)
...@@ -20,7 +20,7 @@ for l in lines: ...@@ -20,7 +20,7 @@ for l in lines:
if l == ' rorqual predPeaks\n': if l == ' rorqual predPeaks\n':
l = next(lines) l = next(lines)
if ',' in l: if ',' in l:
sig = np.zeros((30 * RORQUAL_RAPPORT_FE, RAPPORT_N_CHANNELS), dtype='int16') sig = np.zeros((RORQUAL_LENSIG * RORQUAL_RAPPORT_FE, RAPPORT_N_CHANNELS), dtype='int16')
peaks = np.array(l[:-2].split(',')).astype(int) peaks = np.array(l[:-2].split(',')).astype(int)
l = next(lines) l = next(lines)
for p in peaks[:RORQUAL_NUM_DETEC]: for p in peaks[:RORQUAL_NUM_DETEC]:
...@@ -31,7 +31,7 @@ for l in lines: ...@@ -31,7 +31,7 @@ for l in lines:
if l == ' cacha predPeaks\n': if l == ' cacha predPeaks\n':
l = next(lines) l = next(lines)
if ',' in l: if ',' in l:
sig = np.zeros((10 * CACHA_RAPPORT_FE, RAPPORT_N_CHANNELS), dtype='int16') sig = np.zeros((CACHA_LENSIG * CACHA_RAPPORT_FE, RAPPORT_N_CHANNELS), dtype='int16')
peaks = np.array(l[:-2].split(',')).astype(int) peaks = np.array(l[:-2].split(',')).astype(int)
l = next(lines) l = next(lines)
for p in peaks[:CACHA_NUM_DETEC]: for p in peaks[:CACHA_NUM_DETEC]:
......
No preview for this file type
...@@ -8,7 +8,7 @@ ...@@ -8,7 +8,7 @@
#define NTOAS_MAX 200 #define NTOAS_MAX 200
#define RORQUAL_SAMPLE_RATE 4000 #define RORQUAL_SAMPLE_RATE 4000
#define RORQUAL_LENSIG RORQUAL_SAMPLE_RATE*30 // load 60sec #define RORQUAL_LENSIG RORQUAL_SAMPLE_RATE*60 // load 60sec
#define RORQUAL_WINSIZE 4096 #define RORQUAL_WINSIZE 4096
#define RORQUAL_LENSPEC (RORQUAL_LENSIG - RORQUAL_WINSIZE)/RORQUAL_HOPSIZE #define RORQUAL_LENSPEC (RORQUAL_LENSIG - RORQUAL_WINSIZE)/RORQUAL_HOPSIZE
#define RORQUAL_HOPSIZE 256 #define RORQUAL_HOPSIZE 256
...@@ -18,7 +18,7 @@ ...@@ -18,7 +18,7 @@
#define RORQUAL_RAPPORT_SAMPLESPERSAMPLE 25600 // see cacha #define RORQUAL_RAPPORT_SAMPLESPERSAMPLE 25600 // see cacha
#define CACHA_SAMPLE_RATE 64000 #define CACHA_SAMPLE_RATE 64000
#define CACHA_LENSIG CACHA_SAMPLE_RATE*10 // load 10sec #define CACHA_LENSIG CACHA_SAMPLE_RATE*60 // load 10sec
#define CACHA_WINSIZE 512 #define CACHA_WINSIZE 512
#define CACHA_LENSPEC (CACHA_LENSIG - CACHA_WINSIZE)/CACHA_HOPSIZE #define CACHA_LENSPEC (CACHA_LENSIG - CACHA_WINSIZE)/CACHA_HOPSIZE
#define CACHA_HOPSIZE 256 #define CACHA_HOPSIZE 256
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment