Changeset 6625


Ignore:
Timestamp:
11/19/10 13:51:59 (14 years ago)
Author:
Mathieu Morlighem
Message:

minor

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/model/plot/plot_overlay.m

    r6624 r6625  
    2525ylim2=ylim/getfieldvalue(options,'unit',1);
    2626
    27 %radar power {{{1
     27%radar power
    2828if ~any(isnan(md.sarxm)) & ~any(isnan(md.sarym)) & ~any(isnan(md.sarpwr)),
    2929        disp('plot_overlay info: the radar image held by the model is being used');
     
    3131        disp('Extracting radar image...');
    3232        md=radarpower(md,xlim2,ylim2,getfieldvalue(options,'highres',0));
    33 end%}}}
    34 %InterpFromMeshToGrid{{{1
    35 cornereast =min(xlim);
    36 cornernorth=max(ylim);
     33end
     34
     35%InterpFromMeshToGrid
     36cornereast =min(xlim); cornernorth=max(ylim);
    3737xspacing=(max(xlim)-min(xlim))/(length(md.sarxm));
    3838yspacing=(max(ylim)-min(ylim))/(length(md.sarym));
     
    5555data_nan=find(isnan(data_grid));
    5656
    57 %}}}
    58 %Generate HSV image{{{1
    59 
    60 %options
     57%Generate HSV image
    6158contrast=getfieldvalue(options,'contrast',1); 
    6259transparency=getfieldvalue(options,'alpha',1);
    6360data_grid(data_nan)=data_min;
    64 
    6561%hue (H)
    6662h_data=(data_grid-data_min)/(data_max-data_min);
     
    7066%intensity (V)
    7167radar=(md.sarpwr).^(contrast);
    72 v_data=radar/max(radar(:)); %For the principla image, use radar power as intensity
     68v_data=radar/max(radar(:)); %use radar power as intensity
    7369
    74 %}}}
    75 %Transform HSV to RGB {{{1
     70%Transform HSV to RGB
    7671image_hsv=zeros(size(data_grid,1),size(data_grid,2),3);
    7772image_hsv(:,:,1)=h_data;
     
    7974image_hsv(:,:,3)=v_data;
    8075image_rgb=hsv2rgb(image_hsv);
    81 %}}}
    8276
    8377%Select plot area
Note: See TracChangeset for help on using the changeset viewer.