Changeset 15169

Timestamp:

05/30/13 10:47:53 (12 years ago)

Author:

Mathieu Morlighem

Message:

CHG: faster interpolation since we have a grid

File:

• issm/trunk-jpl/src/m/plot/plot_googlemaps.m (modified) (4 diffs)

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

@@ -46 +46 @@
 end
 
-%Retrieve image from md
-[X Y] = meshgrid(md.radaroverlay.x,md.radaroverlay.y);
+%Process image from model
 final = double(md.radaroverlay.pwr)/double(max(md.radaroverlay.pwr(:))); %rescale between 0 and 1
 
@@ -54 +53 @@
 
 %Prepare grid
-data_grid=InterpFromMeshToMesh2d(md.mesh.elements,md.mesh.x,md.mesh.y,data,X(:),Y(:),'default',NaN); data_grid=reshape(data_grid,size(X));
+if size(md.radaroverlay.x,1)==1 | size(md.radaroverlay.x,2)==1,
+        xmin=min(md.radaroverlay.x);
+        ymax=max(md.radaroverlay.y);
+        xspacing=md.radaroverlay.x(2)-md.radaroverlay.x(1);
+        yspacing=md.radaroverlay.y(2)-md.radaroverlay.y(1);
+        nlines=length(md.radaroverlay.y);
+        ncols =length(md.radaroverlay.x);
+                [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);
+else
+        [X Y] = meshgrid(md.radaroverlay.x,md.radaroverlay.y);
+        data_grid=InterpFromMeshToMesh2d(md.mesh.elements,md.mesh.x,md.mesh.y,data,X(:),Y(:),'default',NaN); data_grid=reshape(data_grid,size(X));
+        x_m=X(1,:); y_m=Y(:,1);
+end
 data_nan=isnan(data_grid);
 if exist(options,'caxis'),
@@ -69 +81 @@
 image_rgb = ind2rgb(uint16((data_grid - data_min)*(length(colorm)/(data_max-data_min))),colorm);
 
-alpha=ones(size(X));
+alpha=ones(size(data_grid));
 alpha(find(~data_nan))=transparency;
 alpha=repmat(alpha,[1 1 3]);
@@ -78 +90 @@
 subplotmodel(plotlines,plotcols,i,options);
 
-h=imagesc(X(1,:),Y(:,1),final);
+h=imagesc(x_m,y_m,final);
 
 %last step: mesh gridded?