warp.m

有关matlab的电子书籍有一定的帮助希望有用

function handle = warp(varargin)
%WARP Display image as texture-mapped surface.
%   WARP(X,MAP) displays the indexed image X with colormap MAP as
%   a texture map on a simple rectangular surface.
%
%   WARP(I,N) displays the intensity image I with gray scale
%   colormap of length N as a texture map on a simple rectangular
%   surface.
%
%   WARP(BW) displays the binary image BW as a texture map on a
%   simple rectangular surface.
%
%   WARP(RGB) displays the RGB image in the array RGB as a
%   texture map on a simple rectangular surface.
%
%   WARP(z,...) displays the image on the surface z.
%
%   WARP(x,y,z,...) displays the image on the surface (x,y,z).
%
%   H = WARP(...) returns a handle to the texture mapped
%   surface.
%
%   Class Support
%   -------------
%   The input image can be of class uint8 or double.
%
%   Example
%   -------
%       [x,y,z] = cylinder;
%       I = imread('testpat1.tif');
%       warp(x,y,z,I);
%
%   See also IMSHOW, IMAGE, IMAGESC, SURF.

%   Clay M. Thompson 3-5-93
%   Revised for MATLAB 5 by Steven L. Eddins, September 1996
%   Copyright 1993-1998 The MathWorks, Inc.  All Rights Reserved.
%   $Revision: 5.10 $  $Date: 1997/11/24 15:36:35 $

[x,y,z,cdata,cdatamapping,clim,map,likeimage,msg] = parse_inputs(varargin{:});
if (~isempty(msg))
    error(msg);
end

axHandle = newplot;
set(axHandle, 'YDir', 'reverse');
h = surface(x,y,z,cdata,'EdgeColor','none','FaceColor','texturemap', ...
        'CDataMapping',cdatamapping);
if (~isempty(clim))
    set(axHandle, 'CLim', clim);
end
if (~isempty(map))
    set(get(axHandle,'Parent'), 'Colormap', map);
end

if likeimage & ~ishold,
    view(2)
    axis([min(x(:)) max(x(:)) min(y(:)) max(y(:))])
else
    view(3)
end

if nargout, handle = h; end


%-----------------------------------------------------------
% Subfunction PARSE_INPUTS
%-----------------------------------------------------------
function [x,y,z,cdata,cdatamapping,clim,map,likeimage,msg] = ...
        parse_inputs(varargin)

msg = '';  % Empty string if no error encountered
x = [];
y = [];
z = [];
cdata = [];
map = [];
cdatamapping = 'direct';
clim = [];
likeimage = 0;
if (get(0,'ScreenDepth') > 16)
    defGrayMapLength = 256;
else
    defGrayMapLength = 64;
end

switch nargin
case 0
    msg = 'Not enough input arguments';
    
case 1
    % warp(I)
    % warp(RGB)
    likeimage = 1;

    if ((ndims(varargin{1}) == 3) & (size(varargin{1},3) == 3))
        % warp(RGB)
        cdata = varargin{1};
        if (isa(cdata,'uint8'))
            cdata = double(cdata)/255.0;
        end
        
    else
        % warp(I)
        cdata = varargin{1};
        cdatamapping = 'scaled';
        clim = [0 1];
        if (isa(cdata,'uint8'))
            if (islogical(cdata))
                cdata = double(cdata);
            else
                cdata = double(cdata)/255.0;
            end
        end
        map = gray(defGrayMapLength);
    end
    
case 2
    % warp(X,map)
    % warp(I,N)
    % warp(z,I)
    % warp(z,RGB)
    % warp(I,[a b])
    
    if ((ndims(varargin{2}) == 3) & (size(varargin{2},3) == 3))
        % warp(z,RGB)
        z = varargin{1};
        cdata = varargin{2};
        if (isa(cdata,'uint8'))
            cdata = double(cdata)/255.0;
        end
        
    elseif (prod(size(varargin{2})) == 1)
        % warp(I,N)
        cdata = varargin{1};
        map = gray(varargin{2});
        cdatamapping = 'scaled';
        clim = [0 1];
        if (isa(cdata,'uint8'))
            if (islogical(cdata))
                cdata = double(cdata);
            else
                cdata = double(cdata)/255.0;
            end
        end
        likeimage = 1;
        
    elseif (isequal(size(varargin{2}), [1 2]))
        % warp(I,[a b])
        cdata = varargin{1};
        cdatamapping = 'scaled';
        clim = varargin{2};
        map = gray(defGrayMapLength);
        if (isa(cdata,'uint8'))
            if (islogical(cdata))
                cdata = double(cdata);
            else
                cdata = double(cdata)/255.0;
                clim = clim/255.0;
            end
        end
        likeimage = 1;
        
    elseif (size(varargin{2},2) == 3)
        % warp(X,map)
        cdata = varargin{1};
        map = varargin{2};
        cdatamapping = 'direct';
        if (isa(cdata,'uint8'))
            cdata = double(cdata) + 1;
        end
        likeimage = 1;
        
    else
        % warp(z,I)
        z = varargin{1};
        cdata = varargin{2};
        cdatamapping = 'scaled';
        clim = [0 1];
        if (isa(cdata,'uint8'))
            if (islogical(cdata))
                cdata = double(cdata);
            else
                cdata = double(cdata)/255.0;
            end
        end
        map = gray(defGrayMapLength);
        
    end
    
case 3
    % warp(R,G,B)  GRANDFATHERED
    % warp(z,X,map)
    % warp(z,I,N)
    % warp(z,I,[a b])
    
    if (prod(size(varargin{3})) == 1)
        % warp(z,I,N)
        z = varargin{1};
        cdata = varargin{2};
        map = gray(varargin{3});
        cdatamapping = 'scaled';
        clim = [0 1];
        if (isa(cdata,'uint8'))
            if (islogical(cdata))
                cdata = double(cdata);
            else
                cdata = double(cdata)/255.0;
            end
        end
        
    elseif (isequal(size(varargin{3}), [1 2]))
        % warp(z,I,[a b])
        z = varargin{1};
        cdata = varargin{2};
        cdatamapping = 'scaled';
        clim = varargin{3};
        map = gray(defGrayMapLength);
        if (isa(cdata,'uint8'))
            if (islogical(cdata))
                cdata = double(cdata);
            else
                cdata = double(cdata)/255.0;
                clim = clim/255.0;
            end
        end
        
    elseif (size(varargin{3},2) == 3)
        % warp(z,X,map)
        z = varargin{1};
        cdata = varargin{2};
        map = varargin{3};
        cdatamapping = 'direct';
        if (isa(cdata,'uint8'))
            cdata = double(cdata) + 1;
        end
        
    elseif (isequal(size(varargin{1}), size(varargin{2})) & ...
                isequal(size(varargin{2}), size(varargin{3})))
        % warp(R,G,B)
        cdata = cat(3,varargin{1:3});
        if (isa(cdata,'uint8'))
            cdata = double(cdata)/255.0;
        end
        likeimage = 1;
        
    else
        msg = 'Invalid input arguments';
    end
    
case 4
    % warp(z,R,G,B)    GRANDFATHERED
    % warp(x,y,z,I)
    % warp(x,y,z,RGB)
    
    if ((ndims(varargin{4}) == 3) & (size(varargin{4},3) == 3))
        % warp(x,y,z,RGB)
        x = varargin{1};
        y = varargin{2};
        z = varargin{3};
        cdata = varargin{4};
        if (isa(cdata,'uint8'))
            cdata = double(cdata)/255.0;
        end
        
    elseif (isequal(size(varargin{2}), size(varargin{3})) & ...
                isequal(size(varargin{3}), size(varargin{4})))
        % warp(z,R,G,B)
        z = varargin{1};
        cdata = cat(3,varargin{2:4});
        if (isa(cdata,'uint8'))
            cdata = double(cdata)/255.0;
        end
        
    else
        % warp(x,y,z,I)
        x = varargin{1};
        y = varargin{2};
        z = varargin{3};
        cdata = varargin{4};
        cdatamapping = 'scaled';
        clim = [0 1];
        map = gray(defGrayMapLength);
        if (isa(cdata,'uint8'))
            if (islogical(cdata))
                cdata = double(cdata);
            else
                cdata = double(cdata)/255.0;
            end
        end
        
    end
    
case 5
    % warp(x,y,z,X,map)
    % warp(x,y,z,I,N)
    % warp(x,y,z,I,[a b])
    
    if (prod(size(varargin{5})) == 1)
        % warp(x,y,z,I,N)
        x = varargin{1};
        y = varargin{2};
        z = varargin{3};
        cdata = varargin{4};
        map = gray(varargin{5});
        cdatamapping = 'scaled';
        clim = [0 1];
        if (isa(cdata,'uint8'))
            if (islogical(cdata))
                cdata = double(cdata);
            else
                cdata = double(cdata)/255.0;
            end
        end
        
    elseif (isequal(size(varargin{5}), [1 2]))
        % warp(x,y,z,I,[a b])
        x = varargin{1};
        y = varargin{2};
        z = varargin{3};
        cdata = varargin{4};
        cdatamapping = 'scaled';
        clim = varargin{5};
        map = gray(defGrayMapLength);
        if (isa(cdata,'uint8'))
            if (islogical(cdata))
                cdata = double(cdata);
            else
                cdata = double(cdata)/255.0;
                clim = clim/255.0;
            end
        end
        
    elseif (size(varargin{5},2) == 3)
        % warp(x,y,z,X,map)
        x = varargin{1};
        y = varargin{2};
        z = varargin{3};
        cdata = varargin{4};
        map = varargin{5};
        cdatamapping = 'direct';
        if (isa(cdata,'uint8'))
            cdata = double(cdata) + 1;
        end
        
    else
        msg = 'Invalid input arguments';
    end
    
case 6
    % warp(x,y,z,R,G,B)    GRANDFATHERED
    
    if (~isequal(size(varargin{4}), size(varargin{5})) | ...
                ~isequal(size(varargin{5}), size(varargin{6})))
        msg = 'R, G, and B must have the same size';
    
    else
        x = varargin{1};
        y = varargin{2};
        z = varargin{3};
        cdata = cat(3,varargin{4:6});
        if (isa(cdata,'uint8'))
            cdata = double(cdata)/255.0;
        end
    end
    
otherwise
    msg = 'Too many input arguments';
    
end

if (isempty(msg))
    siz = size(cdata);
    M = siz(1);
    N = siz(2);
    
    if (isempty(z))
        % The surface displays most quickly when we use
        % a simple 2-by-2 z matrix, but that uses up a
        % large amount of printer memory when printed.
        % In IPT v1, the z matrix was the same size as
        % the image; this solution took a long time to
        % display.  The factor of 4 below is a compromise.
        % -sle, September 1996
        p = max(floor(min(size(cdata))/4),2);
        z = zeros(p);
        x = linspace(1,N,p);
        y = linspace(1,M,p);
    end
    
    if (isempty(x))
        [x,y] = meshgrid(1:size(z,2), 1:size(z,1));
    end
    
    if ((length(x) == 2) & (length(y) == 2))
        [x,y] = meshgrid(linspace(x(1),x(2),size(z,2)), ...
                linspace(y(1),y(2),size(z,1)));
    end
end