mlosch/interp_llc/pcol.m

% Similar to pcolor except that pcol() doesn't drop the last column
% or row of data (ie. doesn't interpolated). It uses the shading flat
% method by default.
%
% See also  PCOLOR, IMAGESC
function [hh] = pcol(varargin);

%cmap=colormap;
%if cmap(size(cmap,1),:)==[0 0 0]
%else
% if sum(sum(isnan(data)))~=0
%  colormap( [colormap' [0 0 0]']')
% end
%end

% Scale data to fit colormap
%clim=[min(min(data)) max(max(data))];
%data=(data-clim(1))/(clim(2)-clim(1))*(size(colormap,1)-1)+1;
%data(find(isnan(data)==1))=size(colormap,1)+1;

if nargin == 1
 data=varargin{1};
 h = pcolor(data([1:end 1],[1:end 1]));
else
 x=varargin{1};
 y=varargin{2};
 data=varargin{3};
 if isvector(x) & isvector(y)
   x=x(:);
   y=y(:);
   h = pcolor([x' 2*x(end)-x(end-1)],...
              [y' 2*y(end)-y(end-1)],...
              data([1:end 1],[1:end 1]));
 else
   x=[x 2*x(:,end)-x(:,end-1)];
   x=[x; 2*x(end,:)-x(end-1,:)];
   y=[y 2*y(:,end)-y(:,end-1)];
   y=[y; 2*y(end,:)-y(end-1,:)];
   h = pcolor(x, ...
              y, ...
              data([1:end 1],[1:end 1]));
 end
end
set(gca,'Layer','top')
shading flat;

if nargout == 1;
  hh = h;
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function answer = isvector(v)
  
  if prod(size(v)) == length(v)
    answer = 1;
  else
    answer = 0;
  end
  
  return
1	mlosch	1.1	% Similar to pcolor except that pcol() doesn't drop the last column
2			% or row of data (ie. doesn't interpolated). It uses the shading flat
3			% method by default.
4			%
5			% See also PCOLOR, IMAGESC
6			function [hh] = pcol(varargin);
7
8			%cmap=colormap;
9			%if cmap(size(cmap,1),:)==[0 0 0]
10			%else
11			% if sum(sum(isnan(data)))~=0
12			% colormap( [colormap' [0 0 0]']')
13			% end
14			%end
15
16			% Scale data to fit colormap
17			%clim=[min(min(data)) max(max(data))];
18			%data=(data-clim(1))/(clim(2)-clim(1))*(size(colormap,1)-1)+1;
19			%data(find(isnan(data)==1))=size(colormap,1)+1;
20
21			if nargin == 1
22			data=varargin{1};
23			h = pcolor(data([1:end 1],[1:end 1]));
24			else
25			x=varargin{1};
26			y=varargin{2};
27			data=varargin{3};
28			if isvector(x) & isvector(y)
29			x=x(:);
30			y=y(:);
31			h = pcolor([x' 2*x(end)-x(end-1)],...
32			[y' 2*y(end)-y(end-1)],...
33			data([1:end 1],[1:end 1]));
34			else
35			x=[x 2*x(:,end)-x(:,end-1)];
36			x=[x; 2*x(end,:)-x(end-1,:)];
37			y=[y 2*y(:,end)-y(:,end-1)];
38			y=[y; 2*y(end,:)-y(end-1,:)];
39			h = pcolor(x, ...
40			y, ...
41			data([1:end 1],[1:end 1]));
42			end
43			end
44			set(gca,'Layer','top')
45			shading flat;
46
47			if nargout == 1;
48			hh = h;
49			end
50
51			%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
52			function answer = isvector(v)
53
54			if prod(size(v)) == length(v)
55			answer = 1;
56			else
57			answer = 0;
58			end
59
60			return