Fix WQM global attribute instrument_burst_duration value, was wrong b…

…y 1 sample. Add global attributes instrument_burst_duration/interval to NXIC parser.
aodn · Oct 4, 2016 · fb223a2 · fb223a2
1 parent 891a219
commit fb223a2
Show file tree

Hide file tree

Showing 3 changed files with 24 additions and 6 deletions.
diff --git a/Parser/NXICBinaryParse.m b/Parser/NXICBinaryParse.m
@@ -104,12 +104,30 @@
   sample_data.meta.instrument_model     = header.instrument_model;
   sample_data.meta.instrument_serial_no = header.instrument_serial_no;
   sample_data.meta.featureType          = mode;
-  if header.sampleRate > 0
-      sample_data.meta.instrument_sample_interval = 1/header.sampleRate;
-  else
-      sample_data.meta.instrument_sample_interval = median(diff(samples.time*24*3600));
+
+  % Let's find each start of bursts
+  dt = [0; diff(samples.time)];
+  iBurst = [1; find(dt>(1/24/60)); length(samples.time)+1];
+
+  % let's read data burst by burst
+  nBurst = length(iBurst)-1;
+  firstTimeBurst = zeros(nBurst, 1);
+  sampleIntervalInBurst = zeros(nBurst, 1);
+  durationBurst = zeros(nBurst, 1);
+  for i=1:nBurst
+      timeBurst = samples.time(iBurst(i):iBurst(i+1)-1);
+      if numel(timeBurst)>1 % deals with the case of a file with a single sample in a single burst
+          sampleIntervalInBurst(i) = median(diff(timeBurst*24*3600));
+          firstTimeBurst(i) = timeBurst(1);
+          durationBurst(i) = (timeBurst(end) - timeBurst(1))*24*3600 + sampleIntervalInBurst(i);
+      end
   end
 
+  sample_data.meta.instrument_sample_interval   = round(median(sampleIntervalInBurst));
+  sample_data.meta.instrument_burst_interval    = round(median(diff(firstTimeBurst*24*3600)));
+  sample_data.meta.instrument_burst_duration    = round(median(durationBurst));
+
+
   sample_data.dimensions = {};
   sample_data.variables  = {};
 

diff --git a/Parser/readWQMdat.m b/Parser/readWQMdat.m
@@ -249,7 +249,7 @@
       if numel(timeBurst)>1 % deals with the case of a file with a single sample in a single burst
           sampleIntervalInBurst(i) = median(diff(timeBurst*24*3600));
           firstTimeBurst(i) = timeBurst(1);
-          durationBurst(i) = (timeBurst(end) - timeBurst(1))*24*3600;
+          durationBurst(i) = (timeBurst(end) - timeBurst(1))*24*3600 + sampleIntervalInBurst(i);
       end
   end
 

diff --git a/Parser/readWQMraw.m b/Parser/readWQMraw.m
@@ -159,7 +159,7 @@
       timeBurst = time(iBurst(i):iBurst(i+1)-1);
       sampleIntervalInBurst(i) = median(diff(timeBurst*24*3600));
       firstTimeBurst(i) = timeBurst(1);
-      durationBurst(i) = (timeBurst(end) - timeBurst(1))*24*3600;
+      durationBurst(i) = (timeBurst(end) - timeBurst(1))*24*3600 + sampleIntervalInBurst(i);
   end
 
   sample_data.meta.instrument_sample_interval   = round(median(sampleIntervalInBurst));