Changeset 17224

Timestamp:

02/06/14 16:18:43 (11 years ago)

Author:

Mathieu Morlighem

Message:

NEW: enableing landsat image overlay, and commented out imadjust, because no license is ever available and creates an error

Location:

issm/trunk-jpl/src/m/plot

Files:

• plot_overlay.m (modified) (1 diff)
• radarpower.m (modified) (2 diffs)

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

@@ -43 +43 @@
 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
 

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

@@ -90 +90 @@
                 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);
@@ -153 +160 @@
         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));