Changeset 6461 for issm/trunk/src/m/kml/kml_part_elems.m

Timestamp:

11/01/10 15:54:12 (14 years ago)

Author:

jschierm

Message:

Implementation of kml class hierarchy in kml plotting drivers.

File:

• issm/trunk/src/m/kml/kml_part_elems.m (modified) (5 diffs)

issm/trunk/src/m/kml/kml_part_elems.m ¶

@@ -2 +2 @@
 %  create kml polygons for the partition elements.
 %
-%  []=kml_part_elems(fid,md,params)
+%  [kfold]=kml_part_elems(md,params)
 %
 %  where the required input is:
-%    fid           (char, file ID of .kml file)
 %    md            (model, model class object)
 %
@@ -19 +18 @@
 %    prtplt        (char, 'off'/'no' for partition segment plot)
 %
-function []=kml_part_elems(varargin)
+%  and the required output is:
+%    kfold         (kml_folder, folder of polygon placemarks)
+%
+function [kfold]=kml_part_elems(varargin)
 
 if ~nargin
@@ -30 +32 @@
 iarg=1;
 if (nargin >= 1)
-    fid=varargin{1};
-end
-if ~exist('fid','var') || isempty(fid) || (fid < 0)
-    error('File ID ''%d'' must be open.',fid);
-end
-
-iarg=iarg+1;
-if (nargin >= 2)
-    md=varargin{2};
+    md=varargin{1};
 end
 if ~exist('md','var') || isempty(md) || ~isa(md,'model')
@@ -109 +103 @@
 if (~exist('prtplt','var') || strncmpi(prtplt,'on' ,2) || strncmpi(prtplt,'y',1)) && ...
     md.npart
-    fprintf(fid,'    <Folder>\n');
-    fprintf(fid,'      <name>Partition Elements</name>\n');
-    fprintf(fid,'      <visibility>1</visibility>\n');
-    fprintf(fid,'      <description>Partitions=%d, Nodes=%d</description>\n',...
+    kfold=kml_folder();
+    kfold.name      ='Partition Elements';
+    kfold.visibility=1;
+    kfold.descript  =sprintf('Partitions=%d, Nodes=%d\n',...
         md.npart,md.numberofgrids);
+    kfold.feature   ={repmat(kml_placemark(),0,0)};
 
-%  write each partition loop as a polygon
+%  write each partition loop as a polygon placemark
 
     disp(['Writing ' num2str(md.npart) ' partitions as KML polygons.']);
@@ -144 +139 @@
         end
         for i=1:length(iloop)-1
-            fprintf(fid,'      <Placemark>\n');
+            kplace=kml_placemark();
             if (length(iloop)-1 > 1)
-                fprintf(fid,'        <name>Partition %d, Loop %d</name>\n',k,i);
+                kplace.name      =sprintf('Partition %d, Loop %d',k,i);
             else
-                fprintf(fid,'        <name>Partition %d</name>\n',k);
+                kplace.name      =sprintf('Partition %d',k);
             end
-            fprintf(fid,'        <visibility>1</visibility>\n');
+            kplace.visibility=1;
             if exist('pdata','var')
-                fprintf(fid,'        <description>Partition data: %g</description>\n',pdata(k));
+                kplace.descript  =sprintf('Partition data: %g',pdata(k));
                 imap = fix((pdata(k)-cmin)/(cmax-cmin)*size(cmap,1))+1;
                 if     (imap >= 1) && (imap <= size(cmap,1))
-                    fprintf(fid,'        <styleUrl>#MatlabColor%d</styleUrl>\n',imap);
+                    kplace.styleurl  =sprintf('#MatlabColor%d',imap);
                 elseif (pdata(k) == cmax)
-                    fprintf(fid,'        <styleUrl>#MatlabColor%d</styleUrl>\n',size(cmap,1));
+                    kplace.styleurl  =sprintf('#MatlabColor%d',size(cmap,1));
                 else
-                    fprintf(fid,'        <styleUrl>#BlackLineEmptyPoly</styleUrl>\n');
+                    kplace.styleurl  =sprintf('#BlackLineEmptyPoly');
                 end
             else
-                fprintf(fid,'        <styleUrl>#BlackLineRandomPoly</styleUrl>\n');
+                kplace.styleurl  =sprintf('#BlackLineRandomPoly');
             end
-            fprintf(fid,'        <Polygon>\n');
-            fprintf(fid,'          <extrude>1</extrude>\n');
-            fprintf(fid,'          <altitudeMode>relativeToGround</altitudeMode>\n');
-            fprintf(fid,'          <outerBoundaryIs>\n');
-            fprintf(fid,'            <LinearRing>\n');
-            fprintf(fid,'              <coordinates>\n');
+
+            kpoly=kml_polygon();
+            kpoly.extrude   =1;
+            kpoly.altmode   ='relativeToGround';
+
+            kring=kml_linearring();
+            kring.coords    =zeros(iloop(i+1)-iloop(i)+1,3);
+
             for j=iloop(i):iloop(i+1)-1
-                [lat(j),long(j)]=mapxy(md.x(edgeper(j,1)),md.y(edgeper(j,1)),'s');
-                fprintf(fid,'                %0.16g,%0.16g,%0.16g\n',long(j),lat(j),alt);
+                [lat,long]=mapxy(md.x(edgeper(j,1)),md.y(edgeper(j,1)),'s');
+                kring.coords(j-iloop(i)+1,:)=[long lat alt];
             end
-            fprintf(fid,'                %0.16g,%0.16g,%0.16g\n',long(iloop(i)),lat(iloop(i)),alt);
+            kring.coords(end,:)=kring.coords(1,:);
 
-            fprintf(fid,'              </coordinates>\n');
-            fprintf(fid,'            </LinearRing>\n');
-            fprintf(fid,'          </outerBoundaryIs>\n');
-            fprintf(fid,'        </Polygon>\n');
-            fprintf(fid,'      </Placemark>\n');
+            kpoly.outer=kring;
+            kplace.geometry=kpoly;
+            kfold.feature{1}(end+1)=kplace;
+            clear kring kpoly kplace
         end
     end
-    fprintf(fid,'    </Folder>\n');
 end