view_interp_s

status = view_interp(color_map,clims,contour_array,contour_array2)

view interpolation binary with header

color_map = what it says
clims = colormap limits (optional)
contour_array = first contour array (gets a solid blue contour line) 
contour_array2 = second contour array (gets a dashed red line).

0001 function status = view_interp_s(color_map,clims,contour_array,contour_array2,s)
0002 %status = view_interp(color_map,clims,contour_array,contour_array2)
0003 %
0004 %view interpolation binary with header
0005 %
0006 %color_map = what it says
0007 %clims = colormap limits (optional)
0008 %contour_array = first contour array (gets a solid blue contour line)
0009 %contour_array2 = second contour array (gets a dashed red line).
0010 fid = get_fid('r');
0011 nrow = fread(fid,1,'float');
0012 ncolumn = fread(fid,1,'float');
0013 plot_res = fread(fid,1,'float');
0014 interp = fread(fid,[ncolumn,nrow],'float');
0015 interp = interp';
0016 if nargin < 5
0017 s = 1;
0018 end;
0019 interp = interp/s;
0020 clims = clims/s;
0021 contour_array = contour_array/s;
0022 contour_array2 = contour_array2/s;
0023 if nargin == 1
0024 clims(2) = 0.9*max(max(interp));
0025 clims(1) = 0.9*min(min(interp));
0026 end;
0027 imagesc(interp,clims), colorbar, colormap(color_map), shading interp
0028 hold on
0029 if (nargin ==2|nargin == 1)
0030 contour(interp)
0031 elseif nargin == 3
0032 contour(interp,contour_array,'r-')
0033 elseif (nargin == 4|nargin == 5)
0034 contour(interp,contour_array,'b-');
0035 contour(interp,contour_array2,'r--');
0036 end;
0037 axis('off')
0038 fclose(fid);
0039 status = 1;
0040 
0041 
0042