Changeset 6660 for issm/trunk/src/m/model/plot/plot_overlay.m

Timestamp:

11/24/10 12:06:14 (14 years ago)

Author:

Eric.Larour

Message:

radarpower now can accept different image path

File:

• issm/trunk/src/m/model/plot/plot_overlay.m (modified) (4 diffs)

issm/trunk/src/m/model/plot/plot_overlay.m

@@ -9 +9 @@
 %process mesh and data
 [x y z elements is2d]=processmesh(md,[],options);
-[data datatype]=processdata(md,data,options);
+if strcmpi(data,'none'),
+        radaronly=1;
+        data=NaN*ones(md.numberofgrids,1);
+        datatype=1;
+else
+        radaronly=0;
+        [data datatype]=processdata(md,data,options);
+end
 
 %check is2d
@@ -30 +37 @@
 else
         disp('Extracting radar image...');
-        md=radarpower(md,xlim2,ylim2,getfieldvalue(options,'highres',0));
+        options=addfielddefault(options,'xlim',xlim2);
+        options=addfielddefault(options,'ylim',ylim2);
+        md=radarpower(md,options);
 end
 
@@ -64 +73 @@
         %hue (H)
         h_data=(data_grid-data_min)/(data_max-data_min);
+        if radaronly,
+                h_data(:)=0;
+        end
         %saturation (S)
         s_data=max(min((0.1+h_data).^(1/transparency),1),0);
@@ -110 +122 @@
 %Apply options, without colorbar and without grid
 options=changefieldvalue(options,'colormap',colorm);           %We used an HSV colorbar
-options=changefieldvalue(options,'caxis',[data_min data_max]); %force caxis so that the colorbar is ready
+if ~isnan(data_min),
+        options=changefieldvalue(options,'caxis',[data_min data_max]); %force caxis so that the colorbar is ready
+end
 options=addfielddefault(options,'axis','equal off');           %default axis
 applyoptions(md,data,options);