Changeset 24884

Timestamp:

05/22/20 09:29:34 (5 years ago)

Author:

Eric.Larour

Message:

CHG: added export routine to output shapefiles.

File:

• issm/trunk-jpl/src/m/classes/mesh2d.m (modified) (1 diff)

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

@@ -189 +189 @@
 
                 end % }}}
+                function export(self,varargin) % {{{
+
+                        options=pairoptions(varargin{:});
+                        filename=getfieldvalue(options,'filename');
+                        format=getfieldvalue(options,'format','shp');
+                        geometry=getfieldvalue(options,'geometry','line');
+                        index=getfieldvalue(options,'index',[]);
+                        proj=getfieldvalue(options,'projection','');
+
+                        %prepare contours: 
+                        contours= struct([]);
+                        if strcmpi(geometry,'point'),
+                                for i=1:self.numberofvertices,
+                                        contours(i).x = self.x(i);
+                                        contours(i).y = self.y(i);
+                                        contours(i).id = i;
+                                        contours(i).Geometry = 'Point';
+                                end
+                        elseif strcmpi(geometry,'line'),
+                                counter=1;
+                                for i=1:self.numberofelements,
+                                        el=self.elements(i,:);
+                                        %first line:
+                                        contours(counter).x = [self.x(el(1)) self.x(el(2))];
+                                        contours(counter).y = [self.y(el(1)) self.y(el(2))];
+                                        contours(counter).Geometry = 'Line';
+
+                                        %second line:
+                                        contours(counter+1).x = [self.x(el(2)) self.x(el(3))];
+                                        contours(counter+1).y = [self.y(el(2)) self.y(el(3))];
+                                        contours(counter+1).Geometry = 'Line';
+
+                                        %third line:
+                                        contours(counter+2).x = [self.x(el(3)) self.x(el(1))];
+                                        contours(counter+2).y = [self.y(el(3)) self.y(el(1))];
+                                        contours(counter+2).Geometry = 'Line';
+                                        
+                                        %increase counter: 
+                                        counter=counter+3;
+                                end
+                        elseif strcmpi(geometry,'polygon'),
+                                if isempty(index),
+                                        counter=1;
+                                        for i=1:self.numberofelements,
+                                                el=self.elements(i,:);
+                                                contours(i).x=[self.x(el(1)) self.x(el(2)) self.x(el(3)) self.x(el(1))];
+                                                contours(i).y=[self.y(el(1)) self.y(el(2)) self.y(el(3)) self.y(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.x(el(1)) self.x(el(2)) self.x(el(3)) self.x(el(1))];
+                                                contours(i).y=[self.y(el(1)) self.y(el(2)) self.y(el(3)) self.y(el(1))];
+                                                contours(i).id = index(i);
+                                                contours(i).Geometry = 'Polygon';
+                                        end
+                                end
+                        else
+                                error(sprintf('mesh3dsurface ''export'' error message: geometry %s not supported yet (should be ''point'' or ''line''',geometry));
+                        end
+
+                        %write file: 
+                        if strcmpi(format,'shp'),
+                                shpwrite(contours,filename);
+                        elseif strcmpi(format,'exp'),
+                                expwrite(contours,filename);
+                        else
+                                error(sprintf('mesh3dsurface ''export'' error message: file format %s not supported yet',format));
+                        end
+
+                        %write projection file: 
+                        if ~isempty(proj),
+                                proj2shpprj(filename,proj);
+                        end
+                        %write style file: 
+                        applyqgisstyle(filename,'mesh');
+
+
+                end % }}}
         end
 end