Changeset 17204

Timestamp:

01/31/14 09:38:31 (11 years ago)

Author:

bdef

Message:

NEW: adding point ability for shapefile read and plot

Location:

issm/trunk-jpl/src/m/shp

Files:

• shpdisp.m (modified) (1 diff)
• shpread.m (modified) (2 diffs)

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

@@ -45 +45 @@
 figure(figurenumber),hold on
 for i=1:length(domain),
-        if (isnumeric(linestyle))
-                plot(domain(i).x*unitmultiplier,domain(i).y*unitmultiplier,'Color',linestyle,'linewidth',linewidth);
+        if(isfield(domain,'nods'))
+                if (isnumeric(linestyle))
+                        plot(domain(i).x*unitmultiplier,domain(i).y*unitmultiplier,'Color',linestyle,'linewidth',linewidth);
+                else
+                        plot(domain(i).x*unitmultiplier,domain(i).y*unitmultiplier,linestyle,'linewidth',linewidth);
+          end
         else
-                plot(domain(i).x*unitmultiplier,domain(i).y*unitmultiplier,linestyle,'linewidth',linewidth);
+                plot(domain(i).x*unitmultiplier,domain(i).y*unitmultiplier,'ro','MarkerSize',5);
         end
+        
 end

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

@@ -6 +6 @@
 %   the shp file, for the density, for the nodes, and a field closed to 
 %   indicate if the domain is closed. 
+%   If this initial shapefile is point only, the fields closed and
+%   points are ommited
 %   The first argument is the .shp file to be read and the second one (optional) 
 %   indicates if the last point shall be read (1 to read it, 0 not to).
@@ -53 +55 @@
                 end
         end
+
+        if strcmpi(shp(i).Geometry,'Point'),
+                x=shp(i).X'; y=shp(i).Y';
+                ids=find(isnan(x));
+                x(ids)=[]; y(ids)=[];
+
+                Struct(end+1).x=x;
+                Struct(end).y=y;
+                Struct(end).density=1;
+                if isfield(shp,'id'),
+                        Struct(end).name=num2str(shp(i).id);
+                else
+                        Struct(end).name='';
+                end
+                for j=1:length(fields),
+                        field=fields{j};
+                        if ~(strcmpi(field,'X') | strcmpi(field,'Y') | strcmpi(field,'id')),
+                                Struct(end).(field)=shp(i).(field);
+                        end
+                end
+        end
 end