source: issm/oecreview/Archive/16554-17801/ISSM-17223-17224.diff@ 17802

Index: ../trunk-jpl/src/m/plot/radarpower.m
===================================================================
--- ../trunk-jpl/src/m/plot/radarpower.m        (revision 17223)
+++ ../trunk-jpl/src/m/plot/radarpower.m        (revision 17224)
@@ -89,9 +89,16 @@
                %Read in temp.tif:
                im=imread('temp.tif','TIFF');
                %adjust contrast and brightness
-               im=imadjust(im,[a b],[c d]);
+               %im=imadjust(im,[a b],[c d]);
                pixelskip=max(1,ceil(posting/((x1-x0)/(size(im,2)))));
-               md.radaroverlay.pwr=double(flipud(im(1:pixelskip:end,1:pixelskip:end)));
+               if size(im,3)==1,
+                       md.radaroverlay.pwr=double(flipud(im(1:pixelskip:end,1:pixelskip:end)));
+               else
+                       md.radaroverlay.pwr=double(im(1:pixelskip:end,1:pixelskip:end,1:pixelskip:end));
+                       md.radaroverlay.pwr(:,:,1)=flipud(md.radaroverlay.pwr(:,:,1));
+                       md.radaroverlay.pwr(:,:,2)=flipud(md.radaroverlay.pwr(:,:,2));
+                       md.radaroverlay.pwr(:,:,3)=flipud(md.radaroverlay.pwr(:,:,3));
+               end
                md.radaroverlay.x=(x0:(x1-x0)/(size(md.radaroverlay.pwr,2)-1):x1);
                md.radaroverlay.y=(y0:(y1-y0)/(size(md.radaroverlay.pwr,1)-1):y1);
 
@@ -152,7 +159,7 @@
        disp('Warning: expecting coordinates in polar stereographic (Std Latitude: 70ºN Meridian: 45º)');
        im=imread(overlay_image);
        %adjust contrast and brightness
-       im=imadjust(im,[a b],[c d]);
+       %im=imadjust(im,[a b],[c d]);
        im=im(toplefty:toplefty+sizey,topleftx:topleftx+sizex);
        md.radaroverlay.pwr=double(flipud(im));
        md.radaroverlay.x=(x0:(x1-x0)/(size(md.radaroverlay.pwr,2)-1):x1);
Index: ../trunk-jpl/src/m/plot/plot_overlay.m
===================================================================
--- ../trunk-jpl/src/m/plot/plot_overlay.m      (revision 17223)
+++ ../trunk-jpl/src/m/plot/plot_overlay.m      (revision 17224)
@@ -42,6 +42,10 @@
 contrast = getfieldvalue(options,'contrast',1);  
 radar = (md.radaroverlay.pwr).^(contrast);
 radar = radar./max(radar(:));
+if size(radar,3)>1,
+       disp('WARNING: color image converted to greyscale intensity image');
+       radar=sum(radar,3)/3;
+end
 %radar(find(radar==0))=1; %Change background from black to white
 
 %InterpFromMeshToGrid