Changeset 17026

Timestamp:

12/17/13 14:46:22 (11 years ago)

Author:

cborstad

Message:

CHG: added support for line shapefiles

File:

• issm/trunk-jpl/src/m/shp/Shp2Exp.m (modified) (2 diffs)

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

@@ -15 +15 @@
         shp=shaperead(shapefilename);
 
-        exp=struct([]);
+        expstruct=struct([]);
         for i=1:length(shp),
                 if strcmpi(shp(i).Geometry,'Polygon'),
@@ -21 +21 @@
                         ids=find(isnan(x));
                         x(ids)=[]; y(ids)=[];
-
-                        exp(end+1).x=x;
-                        exp(end).y=y;
-                        exp(end).nods=length(x);
-                        exp(end).density=1;
-                        exp(end).closed=1;
-                        exp(end).name=num2str(shp(i).id);
+                        expstruct(end+1).x=x;
+                        expstruct(end).y=y;
+                        expstruct(end).nods=length(x);
+                        expstruct(end).density=1;
+                        expstruct(end).closed=1;
+                        expstruct(end).name=num2str(shp(i).id);
                 elseif strcmpi(shp(i).Geometry,'Point'),
-                                x=shp(i).X; y=shp(i).Y;
-
-                                exp(end+1).x=x;
-                                exp(end).y=y;
-                                exp(end).nods=length(x);
-                                exp(end).density=1;
-                                exp(end).closed=1;
-                                %exp(end).name=num2str(shp(i).id);
+                        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=1;
+                        %exp(end).name=num2str(shp(i).id);
+                elseif strcmpi(shp(i).Geometry,'Line'),
+                        x=shp(i).X; y=shp(i).Y;
+                        expstruct(end+1).x=x;
+                        expstruct(end).y=y;
+                        expstruct.x(end)=x(1);
+                        expstruct.y(end)=y(1);
+                        expstruct(end).nods=length(x);
+                        expstruct(end).density=1;
+                        expstruct(end).closed=1;
                 end
         end
 
-        expwrite(exp,expfilename);
+        expwrite(expstruct,expfilename);