Changeset 21783

Timestamp:

07/13/17 14:17:58 (8 years ago)

Author:

Mathieu Morlighem

Message:

CHG: more flexible plot_overlay.m for RGB images

File:

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

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

@@ -41 +41 @@
 end
 contrast = getfieldvalue(options,'contrast',1);  
-radar = (md.radaroverlay.pwr).^(contrast);
-radar = radar./max(radar(:));
+radar    = md.radaroverlay.pwr;
 if size(radar,3)>1,
         disp('WARNING: color image converted to greyscale intensity image');
-        radar=sum(radar,3)/3;
+        if strcmp(class(radar),'uint8'),
+                radar=double(sum(radar,3))/(255*3);
+        else
+                radar=sum(radar,3)/3;
+        end
 end
+radar = radar.^(contrast);
+radar = radar./max(radar(:));
+
 if getfieldvalue(options,'backgroundbtw',0)
         radar(find(radar==0))=1; %Change background from black to white
@@ -54 +60 @@
 xmin=min(md.radaroverlay.x);
 ymax=max(md.radaroverlay.y);
-xspacing=(max(md.radaroverlay.x)-min(md.radaroverlay.x))/(length(md.radaroverlay.x));
-yspacing=(max(md.radaroverlay.y)-min(md.radaroverlay.y))/(length(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);