Changeset 6660

Timestamp:

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

Author:

Eric.Larour

Message:

radarpower now can accept different image path

Location:

issm/trunk/src/m/model/plot

Files:

• plot_none.m (modified) (1 diff)
• plot_overlay.m (modified) (4 diffs)
• plotdoc.m (modified) (1 diff)

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

@@ -10 +10 @@
 options=addfielddefault(options,'axis','equal');
 
+if exist(options,'overlay'),
+        plot_overlay(md,'none',options,nlines,ncols,i);
+        return;
+end
+
 %apply options
 applyoptions(md,[],options);

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);

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

@@ -96 +96 @@
 disp('       ''yticklabel'': specifiy yticklabel');
 disp('       ''overlay'': yes or no. This will overlay a radar amplitude image behind');
+disp('       ''overlay_image'': path to overlay image. provide overlay_xlim, overlay_ylim, overlay_xposting and overlay_yposting options also');
 disp('       ''contrast'': (default 1) coefficient to add contrast to the radar amplitude image used in overlays');
 disp('       ''highres'': resolution of overlayed radar amplitude image (default is 0, high resolution is 1).');