Changeset 15155

Timestamp:

05/29/13 16:28:09 (12 years ago)

Author:

Mathieu Morlighem

Message:

CHG: also works for Antarctica now (rotation needed though...)

Location:

issm/trunk-jpl/src/m/plot

Files:

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

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

@@ -26 +26 @@
         xlim=getfieldvalue(options,'xlim',[min(md.mesh.x) max(md.mesh.x)]);
         ylim=getfieldvalue(options,'ylim',[min(md.mesh.y) max(md.mesh.y)]);
-        [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);
+        if strcmpi(md.mesh.hemisphere,'n'),
+                [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 strcmpi(md.mesh.hemisphere,'s'),
+                [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
+                error('field hemisphere should either be ''n'' or ''s''');
+        end
 
         %Image corners in lat/long
@@ -76 +85 @@
                 position = [num2str(latn) ',' num2str(lonn)];
                 disp(['Google Earth tile: ' num2str(x) '/' num2str(cols-1) ' ' num2str(y) '/' num2str(rows-1) ' (center: ' position ')']);
+                %Google maps API
+                %http://developers.google.com/maps/documentation/staticmaps/
                 params = [...
                         'center=' position ...
@@ -97 +108 @@
 [gX gY]=meshgrid(ulx:ulx+size(final,2)-1,uly:-1:uly-size(final,1)+1);
 [LAT LON]=pixelstolatlon(gX,gY, zoom);
-[X Y]=ll2xy(LAT,LON,+1,45,70);
+if strcmpi(md.mesh.hemisphere,'n'),
+        [X Y]=ll2xy(LAT,LON,+1,45,70);
+elseif strcmpi(md.mesh.hemisphere,'s'),
+        [X Y]=ll2xy(LAT,LON,-1,0,71);
+else
+        error('field hemisphere should either be ''n'' or ''s''');
+end
 
 md.radaroverlay.pwr=final;

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)]);
-        [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);
+        if strcmpi(md.mesh.hemisphere,'n'),
+                [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 strcmpi(md.mesh.hemisphere,'s'),
+                [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
+                error('field hemisphere should either be ''n'' or ''s''');
+        end
 
         %Image corners in lat/long