Changeset 20070

Timestamp:

02/03/16 15:21:37 (9 years ago)

Author:

Mathieu Morlighem

Message:

CHG: more general lat/long inputs

Location:

issm/trunk-jpl/src/m/plot

Files:

• googlemaps.m (modified) (1 diff)
• plot_googlemaps.m (modified) (1 diff)

issm/trunk-jpl/src/m/plot/googlemaps.m

@@ -40 +40 @@
                         [linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
                         -1,0,71);
+        elseif md.mesh.epsg==26906, %UTM 6V Columbia Glacier Alaska
+                [latlist lonlist]= utm2ll(...
+                        [linspace(xlim(1),xlim(2),100) linspace(xlim(2),xlim(2),100) linspace(xlim(2),xlim(1),100) linspace(xlim(1),xlim(1),100)],...
+                        [linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
+                        6);
+        elseif numel(md.mesh.lat)==numel(md.mesh.x),
+                latlist = md.mesh.lat; %That might work?
+                lonlist = md.mesh.long;
         else
-                error('EPSG code not supported yet');
+                error('EPSG code not supported yet, and no lat long found in md.mesh');
         end
 

issm/trunk-jpl/src/m/plot/plot_googlemaps.m

@@ -25 +25 @@
         xlim=getfieldvalue(options,'xlim',[min(x) max(x)]);
         ylim=getfieldvalue(options,'ylim',[min(y) max(y)]);
-        if md.mesh.epsg==3413,
+        if md.mesh.epsg==3413, %UPS Greenland
                 [latlist lonlist]= xy2ll(...
                         [linspace(xlim(1),xlim(2),100) linspace(xlim(2),xlim(2),100) linspace(xlim(2),xlim(1),100) linspace(xlim(1),xlim(1),100)],...
                         [linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
                         +1,45,70);
-        elseif md.mesh.epsg==3031,
+        elseif md.mesh.epsg==3031, %UPS Antarctica
                 [latlist lonlist]= xy2ll(...
                         [linspace(xlim(1),xlim(2),100) linspace(xlim(2),xlim(2),100) linspace(xlim(2),xlim(1),100) linspace(xlim(1),xlim(1),100)],...
                         [linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
                         -1,0,71);
-        else
+        elseif md.mesh.epsg==26906, %UTM 6V Columbia Glacier Alaska
+                [latlist lonlist]= utm2ll(...
+                        [linspace(xlim(1),xlim(2),100) linspace(xlim(2),xlim(2),100) linspace(xlim(2),xlim(1),100) linspace(xlim(1),xlim(1),100)],...
+                        [linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
+                        6);
+        elseif numel(md.mesh.lat)==numel(md.mesh.x),
                 latlist = md.mesh.lat; %That might work?
                 lonlist = md.mesh.long;
-                error('EPSG code not supported yet');
+        else
+                error('EPSG code not supported yet, and no lat long found in md.mesh');
         end