Changeset 16437

Timestamp:

10/16/13 17:01:19 (11 years ago)

Author:

cborstad

Message:

NEW: overlay_adjust parameters added for better tuning of brightness and contrast of overlay images.
Parameters a,b,c,d must be in [0,1]. In general:
To decrease contrast: increase a
To increase contrast: decrease b
To increase brightness: increase c
To decrease brightness: decrease d
Some trial-and-error is needed, and it helps to look at a histogram of the intensity first.

File:

• issm/trunk-jpl/src/m/plot/radarpower.m (modified) (4 diffs)

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

@@ -25 +25 @@
 ylim    = getfieldvalue(options,'ylim',[min(md.mesh.y) max(md.mesh.y)]);
 posting = getfieldvalue(options,'posting',0); % 0 -> image posting default
+a = getfieldvalue(options,'overlay_adjust_a',0);
+b = getfieldvalue(options,'overlay_adjust_b',1);
+c = getfieldvalue(options,'overlay_adjust_c',0);
+d = getfieldvalue(options,'overlay_adjust_d',1);
 
 %find gdal coordinates
@@ -85 +89 @@
                 %Read in temp.tif:
                 im=imread('temp.tif','TIFF');
+                %adjust contrast and brightness
+                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)));
@@ -114 +120 @@
                 %Read in temp.tif:
                 im=imread('temp.tif','TIFF');
+                %adjust contrast and brightness
+                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)));
@@ -144 +152 @@
         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=im(toplefty:toplefty+sizey,topleftx:topleftx+sizex);
         md.radaroverlay.pwr=double(flipud(im));