Changeset 21369

Timestamp:

11/12/16 19:33:34 (8 years ago)

Author:

Eric.Larour

Message:

CHG: shpread: new ability to invert profile. shpwrite: added line support.

Location:

issm/branches/trunk-larour-NatGeoScience2016/src/m/shp

Files:

• shpread.m (modified) (3 diffs)
• shpwrite.m (modified) (1 diff)

issm/branches/trunk-larour-NatGeoScience2016/src/m/shp/shpread.m

@@ -1 @@
-function Struct=shpread(filename)
+function Struct=shpread(filename,varargin)
 %SHPREAD - read a shape file and build a Structure
 %
@@ -18 +18 @@
 %
 %   See also EXPDOC, EXPWRITEASVERTICES
+
+%recover options
+options=pairoptions(varargin{:});
 
 %some checks
@@ -102 +105 @@
 end
 
+invert=getfieldvalue(options,'invert',0);
+if invert,
+        for i=1:length(Struct),
+                Struct(i).x=flipud(Struct(i).x);
+                Struct(i).y=flipud(Struct(i).y);
+        end
 end
+
+end

issm/branches/trunk-larour-NatGeoScience2016/src/m/shp/shpwrite.m

@@ -18 +18 @@
         if strcmpi(shp(i).Geometry,'Point'),
                 contours(i).Geometry='Point';
-        else strcmpi(shp(i).Geometry,'Polygon'),
+        elseif strcmpi(shp(i).Geometry,'Polygon'),
                 contours(i).Geometry='Polygon';
+        elseif strcmpi(shp(i).Geometry,'Line'),
+                contours(i).Geometry='Line';
         end
         contours(i).id=i;