Changeset 21828 for issm/trunk-jpl/src/m/plot/plot_overlay.m

Timestamp:

07/20/17 04:16:12 (8 years ago)

Author:

Mathieu Morlighem

Message:

CHG: simplified the arguments of InterpFromMeshToGrid

File:

• issm/trunk-jpl/src/m/plot/plot_overlay.m (modified) (2 diffs)

issm/trunk-jpl/src/m/plot/plot_overlay.m ¶

@@ -58 +58 @@
 
 %InterpFromMeshToGrid
-xmin=min(md.radaroverlay.x);
-ymax=max(md.radaroverlay.y);
-xspacing=(max(md.radaroverlay.x)-min(md.radaroverlay.x))/(length(md.radaroverlay.x) -1);
-yspacing=(max(md.radaroverlay.y)-min(md.radaroverlay.y))/(length(md.radaroverlay.y) -1);
-if(md.radaroverlay.x(end)-md.radaroverlay.x(1)<0)
-        xspacing= -xspacing; 
-end
-if(md.radaroverlay.y(end)-md.radaroverlay.y(1)<0)
-        yspacing= -yspacing; 
-end
-nlines=length(md.radaroverlay.y);
-ncols =length(md.radaroverlay.x);
 disp('Interpolating data on grid...');
+x_m = md.radaroverlay.x;
+y_m = md.radaroverlay.y;
 if radaronly,
-        x_m=xmin:xspacing:xmin+ncols*xspacing;
-        y_m=ymax-nlines*yspacing:yspacing:ymax;
-        data_grid=NaN*ones(nlines,ncols);
+        data_grid=NaN(size(radar));
 else
-        [x_m y_m data_grid]=InterpFromMeshToGrid(elements,x/getfieldvalue(options,'unit',1),y/getfieldvalue(options,'unit',1),...
-                data,xmin,ymax,xspacing,yspacing,nlines,ncols,NaN);
+        data_grid=InterpFromMeshToGrid(elements,x/getfieldvalue(options,'unit',1),y/getfieldvalue(options,'unit',1),data,x_m,y_m,NaN);
 end
 
@@ -152 +139 @@
 
 %Plot: 
-imagesc(md.radaroverlay.x*getfieldvalue(options,'unit',1),md.radaroverlay.y*getfieldvalue(options,'unit',1),image_rgb);set(gca,'YDir','normal');
+imagesc(x_m*getfieldvalue(options,'unit',1),y_m*getfieldvalue(options,'unit',1),image_rgb);set(gca,'YDir','normal');
 
 %last step: mesh overlay?