Changeset 3174

Timestamp:

03/03/10 16:44:35 (15 years ago)

Author:

Eric.Larour

Message:

Fixed some problems in overlay, that was not working in km mode.

Location:

issm/trunk/src/m/classes/public/plot

Files:

• applyoptions.m (modified) (2 diffs)
• plot_overlay.m (modified) (3 diffs)

issm/trunk/src/m/classes/public/plot/applyoptions.m

@@ -264 +264 @@
         axes('pos',PosInset);
         axis equal off
-        box off
+        %box off
         if strcmpi(md.hemisphere,'n') | strcmpi(md.hemisphere,'north'),
                 A=expread([ ISSM_DIR '/../models/Greenland/Exp_Par/DomainOutline.exp']);
         elseif strcmpi(md.hemisphere,'s') | strcmpi(md.hemisphere,'south'),
                 A=expread([ ISSM_DIR '/../models/Antarctica/Exp_Par/DomainOutline.exp']);
+        else
+                error('applyoptions error message: hemisphere not defined');
         end
         Ax=[min(A.x) max(A.x)];
@@ -274 +276 @@
         mdx=[min(md.x) max(md.x)];
         mdy=[min(md.y) max(md.y)];
-        patch([Ax(1)  Ax(2)  Ax(2)  Ax(1)],[Ay(1)  Ay(1)  Ay(2)  Ay(2)],[1 1 1],'edgecolor','none')
         line(A.x,A.y,'color','b');
+        patch([Ax(1)  Ax(2)  Ax(2)  Ax(1) Ax(1)],[Ay(1)  Ay(1)  Ay(2)  Ay(2) Ay(1)],[1 1 1],'EdgeColor',[0 0 0],'LineWidth',1,'FaceLighting','none')
         patch( [mdx(1) mdx(2) mdx(2) mdx(1)],[mdy(1) mdy(1) mdy(2) mdy(2)],ones(4,1),'EdgeColor',[0 0 0],'FaceColor',[0.5 0.5 0.5],'FaceAlpha',0.5)
         colorbar('off');

issm/trunk/src/m/classes/public/plot/plot_overlay.m

@@ -20 +20 @@
 
 %get xlim and ylim
-xlim=getfieldvalue(options,'xlim',[min(md.x) max(md.x)]);
-ylim=getfieldvalue(options,'ylim',[min(md.y) max(md.y)]);
+xlim=getfieldvalue(options,'xlim',[min(x) max(x)])/getfieldvalue(options,'unit',1);
+ylim=getfieldvalue(options,'ylim',[min(y) max(y)])/getfieldvalue(options,'unit',1);
 
 %radar power
@@ -31 +31 @@
         t2=clock;fprintf('%s\n',[' done (' num2str(etime(t2,t1)) ' seconds)']);
 end
+
+
+xlim=getfieldvalue(options,'xlim',[min(x) max(x)]);
+ylim=getfieldvalue(options,'ylim',[min(y) max(y)]);
 
 %mesh2grid
@@ -101 +105 @@
 
 %Plot: 
-imagesc(md.sarxm,md.sarym,image_rgb);set(gca,'YDir','normal');
+imagesc(md.sarxm*getfieldvalue(options,'unit',1),md.sarym*getfieldvalue(options,'unit',1),image_rgb);set(gca,'YDir','normal');
 
 %last step: mesh overlay?