Changeset 24885

Timestamp:

05/22/20 09:30:15 (5 years ago)

Author:

Eric.Larour

Message:

NEW: added new export routine to export 3D shape files in lat,long coordinates.

File:

• issm/trunk-jpl/src/m/classes/mesh3dsurface.m (modified) (3 diffs)

issm/trunk-jpl/src/m/classes/mesh3dsurface.m

@@ -200 +200 @@
                         format=getfieldvalue(options,'format','shp');
                         geometry=getfieldvalue(options,'geometry','line');
+                        index=getfieldvalue(options,'index',[]);
+                        proj=getfieldvalue(options,'projection','');
 
                         %prepare contours: 
@@ -233 +235 @@
                                 end
                         elseif strcmpi(geometry,'polygon'),
-                                counter=1;
-                                for i=1:self.numberofelements,
-                                        el=self.elements(i,:);
-                                        contours(i).x=[self.long(el(1)) self.long(el(2)) self.long(el(3)) self.long(el(1))];
-                                        contours(i).y=[self.lat(el(1)) self.lat(el(2)) self.lat(el(3)) self.lat(el(1))];
-                                        contours(i).Geometry = 'Polygon';
+                                if isempty(index),
+                                        counter=1;
+                                        for i=1:self.numberofelements,
+                                                el=self.elements(i,:);
+                                                contours(i).x=[self.long(el(1)) self.long(el(2)) self.long(el(3)) self.long(el(1))];
+                                                contours(i).y=[self.lat(el(1)) self.lat(el(2)) self.lat(el(3)) self.lat(el(1))];
+                                                contours(i).Geometry = 'Polygon';
+                                                contours(i).Id = i;
+                                        end
+                                else
+                                        counter=1;
+                                        for i=1:length(index),
+                                                el=self.elements(index(i),:);
+                                                contours(i).x=[self.long(el(1)) self.long(el(2)) self.long(el(3)) self.long(el(1))];
+                                                contours(i).y=[self.lat(el(1)) self.lat(el(2)) self.lat(el(3)) self.lat(el(1))];
+                                                contours(i).id = index(i);
+                                                contours(i).Geometry = 'Polygon';
+                                        end
                                 end
                         else
@@ -253 +267 @@
                         end
 
+                        %write projection file: 
+                        if ~isempty(proj),
+                                proj2shpprj(filename,proj);
+                        end
+
+                        %write style file: 
+                        applyqgisstyle(filename,'mesh');
+
+