Changeset 6625
- Timestamp:
- 11/19/10 13:51:59 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/model/plot/plot_overlay.m
r6624 r6625 25 25 ylim2=ylim/getfieldvalue(options,'unit',1); 26 26 27 %radar power {{{127 %radar power 28 28 if ~any(isnan(md.sarxm)) & ~any(isnan(md.sarym)) & ~any(isnan(md.sarpwr)), 29 29 disp('plot_overlay info: the radar image held by the model is being used'); … … 31 31 disp('Extracting radar image...'); 32 32 md=radarpower(md,xlim2,ylim2,getfieldvalue(options,'highres',0)); 33 end %}}}34 %InterpFromMeshToGrid{{{1 35 cornereast =min(xlim); 36 corner north=max(ylim);33 end 34 35 %InterpFromMeshToGrid 36 cornereast =min(xlim); cornernorth=max(ylim); 37 37 xspacing=(max(xlim)-min(xlim))/(length(md.sarxm)); 38 38 yspacing=(max(ylim)-min(ylim))/(length(md.sarym)); … … 55 55 data_nan=find(isnan(data_grid)); 56 56 57 %}}} 58 %Generate HSV image{{{1 59 60 %options 57 %Generate HSV image 61 58 contrast=getfieldvalue(options,'contrast',1); 62 59 transparency=getfieldvalue(options,'alpha',1); 63 60 data_grid(data_nan)=data_min; 64 65 61 %hue (H) 66 62 h_data=(data_grid-data_min)/(data_max-data_min); … … 70 66 %intensity (V) 71 67 radar=(md.sarpwr).^(contrast); 72 v_data=radar/max(radar(:)); % For the principla image,use radar power as intensity68 v_data=radar/max(radar(:)); %use radar power as intensity 73 69 74 %}}} 75 %Transform HSV to RGB {{{1 70 %Transform HSV to RGB 76 71 image_hsv=zeros(size(data_grid,1),size(data_grid,2),3); 77 72 image_hsv(:,:,1)=h_data; … … 79 74 image_hsv(:,:,3)=v_data; 80 75 image_rgb=hsv2rgb(image_hsv); 81 %}}}82 76 83 77 %Select plot area
Note:
See TracChangeset
for help on using the changeset viewer.