Changeset 15111

Timestamp:

05/26/13 17:33:07 (12 years ago)

Author:

Mathieu Morlighem

Message:

CHG: better bounding box

File:

• issm/trunk-jpl/src/m/plot/plot_googleearth.m (modified) (1 diff)

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

@@ -19 +19 @@
 xlim=getfieldvalue(options,'xlim',[min(x) max(x)]);
 ylim=getfieldvalue(options,'ylim',[min(y) max(y)]);
-post=getfieldvalue(options,'posting',diff(xlim)/1000);
-
-%Interpolating data on grid
-[x_m y_m data_grid]=InterpFromMeshToGrid(elements,x,y,data,xlim(1),ylim(2),post,post,round(diff(ylim)/post),round(diff(xlim)/post),NaN);
-if size(data_grid,1)<3 | size(data_grid,2)<3,
-        error('data_grid size too small in plot_gridded, check posting and units');
-end
+[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);
 
 %limits in lat/long
-[ullat,ullon] = xy2ll(xlim(1),ylim(2),+1,45,70);
-[lrlat,lrlon] = xy2ll(xlim(2),ylim(1),+1,45,70);
+ullat = max(latlist); ullon = min(lonlist);
+lrlat = min(latlist); lrlon = max(lonlist);
 
 %Find optimal zoom