Replace NaN in the SOG variable with ship velocity computed using the positions

tsg_util/shipVelocity.m

@@ -95,6 +95,14 @@ tsg.report.badvelocity = nBadVelocity;
 % ------------------------------------------------------
 if isempty(tsg.SPDC)
   tsg.SPDC = velocity(1:tsgLength);
+else
+  % Replace NaN with computed ship velocity
+  % ---------------------------------------
+  noSPDC = isnan( tsg.SPDC );
+  if ~isempty(noSPDC)
+    ind = find(noSPDC == 1);
+    tsg.SPDC(ind) = velocity(ind);
+  end 
 end
 
 % Save tsg structure