Changeset 9226

Timestamp:

08/09/11 14:38:01 (14 years ago)

Author:

Mathieu Morlighem

Message:

Added posting option for BIG overlays

Location:

issm/trunk/src/m/model

Files:

• plot/plot_overlay.m (modified) (1 diff)
• radarpower.m (modified) (5 diffs)

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

@@ -115 +115 @@
 %Select plot area 
 subplot(plotlines,plotcols,i);
-%uncomment following lines to have more space
-%P=get(gca,'pos');
-%P(3)=P(3)+0.07;
-%P(4)=P(4)+0.07;
-%%P(2)=P(2)+0.08;
-%P(1)=P(1)-0.03;
-%set(gca,'pos',P);
 
 %Plot: 

issm/trunk/src/m/model/radarpower.m

@@ -11 +11 @@
 
 %If gdal does not work, uncomment the following line
-%setenv('LD_LIBRARY_PATH','');
+setenv('LD_LIBRARY_PATH','');
 
 %Parse inputs
@@ -21 +21 @@
 xlim=getfieldvalue(options,'xlim',[min(md.x) max(md.x)]);
 ylim=getfieldvalue(options,'ylim',[min(md.y) max(md.y)]);
+posting=getfieldvalue(options,'posting',0); % 0 -> image posting default
 
 %find gdal coordinates
@@ -45 +46 @@
                 numcols=floor((xmax-xmin)/geom(1)); % x posting
                 numrows=floor((ymax-ymin)/geom(4)); % y posting
+                pixelskip=max(1,ceil(posting/geom(1)));
 
                 %Read and crop file
@@ -50 +52 @@
                 im=imread(jpgim);
                 im=im(firstrow:firstrow+numrows-1,firstcol:firstcol+numcols-1);
-                md.sarpwr=double(flipud(im));
+                md.sarpwr=double(flipud(im(1:pixelskip:end,1:pixelskip:end)));
                 md.sarxm=(xmin:(xmax-xmin)/(size(md.sarpwr,2)-1):xmax);
                 md.sarym=(ymin:(ymax-ymin)/(size(md.sarpwr,1)-1):ymax);
@@ -72 +74 @@
 
                 %Read in temp.tif:
-                md.sarpwr=double(flipud(imread('temp.tif','TIFF')));
+                im=imread('temp.tif','TIFF');
+                pixelskip=max(1,ceil(posting/((x1-x0)/(size(im,2)))));
+                md.sarpwr=double(flipud(im(1:pixelskip:end,1:pixelskip:end)));
                 md.sarxm=(x0:(x1-x0)/(size(md.sarpwr,2)-1):x1);
                 md.sarym=(y0:(y1-y0)/(size(md.sarpwr,1)-1):y1);