Changeset 26631

Timestamp:

11/16/21 13:28:18 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: better way of handling NaNs

File:

• issm/trunk-jpl/src/m/shp/shp2exp.m (modified) (1 diff)

issm/trunk-jpl/src/m/shp/shp2exp.m

@@ -48 +48 @@
         elseif strcmpi(shp(i).Geometry,'Point'),
                 x=shp(i).X; y=shp(i).Y;
-                expstruct(end+1).x=x;
-                expstruct(end).y=y;
-                expstruct(end).nods=length(x);
-                expstruct(end).density=1;
-                expstruct(end).closed=0;
-                %exp(end).name=num2str(shp(i).id);
+                if ~isnan(x) && ~isnan(y)
+                        expstruct(end+1).x=x;
+                        expstruct(end).y=y;
+                        expstruct(end).nods=length(x);
+                        expstruct(end).density=1;
+                        expstruct(end).closed=0;
+                        %exp(end).name=num2str(shp(i).id);
+                end
+        elseif strcmpi(shp(i).Geometry,'PointZ'),
+                shp(i)
         elseif strcmpi(shp(i).Geometry,'Line'),
-                x=shp(i).X; y=shp(i).Y;
-                expstruct(end+1).x=x;
-                expstruct(end).y=y;
-                expstruct(end).nods=length(x);
-                expstruct(end).density=1;
-                expstruct(end).closed=0;
+                x=shp(i).X(:); y=shp(i).Y(:);
+                pos = find(~isnan(x) & ~isnan(y));
+                idx=find(diff(pos)~=1);
+                A=[idx(1);diff(idx);numel(pos)-idx(end)];
+                Cx=mat2cell(x(pos),A,1);
+                Cy=mat2cell(y(pos),A,1);
+                for i=1:numel(Cx)
+                        expstruct(end+1).x=Cx{i};
+                        expstruct(end).y=Cy{i};
+                        expstruct(end).nods=length(Cx{i});
+                        expstruct(end).density=1;
+                        expstruct(end).closed=0;
+                end
         end
 end