function freezeColors(varargin)
% freezeColors Lock colors of plot, enablingmultiple colormaps per figure. (v2.3)
%
% Problem: There is only onecolormap per figure. This function provides
% an easy solution when plots using different colomaps aredesired
% in the same figure.
%
% freezeColors freezes thecolors of graphics objects in the current axis so
% that subsequent changes to the colormap (or caxis) will not changethe
% colors of these objects. freezeColors works on any graphicsobject
% with CData in indexed-color mode: surfaces, images,scattergroups,
% bargroups, patches, etc. It works by converting CData to true-colorrgb
% based on the colormap active at the time freezeColors iscalled.
%
% The original indexed colordata is saved, and can be restored using
% unfreezeColors, making the plot once again subject to the colormapand
% caxis.
%
%
% Usage:
% freezeColors applies to all objects in current axis (gca),
% freezeColors(axh) same, butworks on axis axh.
%
% Example:
% subplot(2,1,1); imagesc(X); colormap hot; freezeColors
% subplot(2,1,2); imagesc(Y); colormap hsv; freezeColors etc...
%
% Note: colorbars must also be frozen. Due to Matlab 'improvements'this can
% nolonger be done with freezeColors. Instead, please
% usethe function CBFREEZE by Carlos Adrian Vargas Aguilera
% thatcan be downloaded from the MATLAB File Exchange
% (http://www.mathworks.com/matlabcentral/fileexchange/24371)
%
% h=colorbar; cbfreeze(h), or simply cbfreeze(colorbar)
%
% For additional examples, see test/test_main.m
%
% Side effect on render mode:freezeColors does not work with the painters
% renderer, because Matlab doesn't support rgb color data in
% painters mode. If the current renderer is painters,freezeColors
% changes it to zbuffer. This may have unexpected effects on otheraspects
% of your plots.
%
% See also unfreezeColors, freezeColors_pub.html, cbfreeze.
%
%
% John Iversen (iversen@nsi.edu) 3/23/05
%
% Changes:
% JRI (iversen@nsi.edu)4/19/06 Correctly handlesscaled integer cdata
% JRI9/1/06 should now handle allobjects with cdata: images, surfaces,
% scatterplots. (v 2.1)
% JRI 11/11/06 Preserves NaNcolors. Hidden option (v 2.2, not uploaded)
% JRI 3/17/07 Preserve caxis after freezing--maintains colorbar scale (v2.3)
% JRI 4/12/07 Check for painters mode as Matlab doesn't support rgb in it.
% JRI4/9/08 Fix preserving caxis forobjects within hggroups (e.g. contourf)
% JRI4/7/10 Change documentation forcolorbars
% Hidden option for NaN colors:
% Missing data are oftenrepresented by NaN in the indexed color
% data, which renderstransparently. This transparency will be preserved
% when freezing colors. Ifinstead you wish such gaps to be filled with
% a real color, add'nancolor',[r g b] to the end of the arguments. E.g.
% freezeColors('nancolor',[r gb]) or freezeColors(axh,'nancolor',[r g b]),
% where [r g b] is a colorvector. This works on images & pcolor, but noton
% surfaces.
% Thanks to Fabiano Busdraghiand Jody Klymak for the suggestions. Bugfixes
% attributed in the code.
% Free for all uses, but please retain the following:
% Original Author:
% John Iversen, 2005-10
% john_iversen@post.harvard.edu
appdatacode = 'JRI__freezeColorsData';
[h, nancolor] = checkArgs(varargin);
%gather all children with scaled or indexed CData
cdatah = getCDataHandles(h);
%current colormap
cmap = colormap;
nColors = size(cmap,1);
cax = caxis;
% convert object color indexes into colormap to true-color datausing
% current colormap
for hh = cdatah',
g =get(hh);
%preserveparent axis clim
parentAx =getParentAxes(hh);
originalClim= get(parentAx,'clim');
% Note: Special handling ofpatches: For some reason, setting
% cdata on patches created bybar() yields an error,
% so instead we'll setfacevertexcdata instead for patches.
if~strcmp(g.Type,'patch'),
cdata = g.CData;
else
cdata = g.FaceVertexCData;
end
%get cdatamapping (most objects (except scattergroup) have it)
ifisfield(g,'CDataMapping'),
scalemode = g.CDataMapping;
else
scalemode = 'scaled';
end
%saveoriginal indexed data for use with unfreezeColors
siz =size(cdata);
setappdata(hh, appdatacode, {cdata scalemode});
%convertcdata to indexes into colormap
ifstrcmp(scalemode,'scaled'),
%4/19/06 JRI, Accommodate scaled display of integer cdata:
% in MATLAB, uint * double = uint, so must coerce cdata todouble
% Thanks to O Yamashita for pointing this need out
idx = ceil( (double(cdata) - cax(1)) / (cax(2)-cax(1)) *nColors);
else %directmapping
idx = cdata;
/8/09 in case direct data is non-int (e.g.image;freezeColors)
% (Floor mimics how matlab converts data into colormapindex.)
% Thanks to D Armyr for the catch
idx = floor(idx);
end
%clamp to[1, nColors]
idx(idx<1) = 1;
idx(idx>nColors) = nColors;
%handlenans in idx
nanmask =isnan(idx);
idx(nanmask)=1; %temporarily replace w/ a valid colormap index
%maketrue-color data--using current colormap
realcolor =zeros(siz);
for i =1:3,
c = cmap(idx,i);
c = reshape(c,siz);
c(nanmask) = nancolor(i); %restore Nan (or nancolor ifspecified)
realcolor(:,:,i) = c;
end
%apply newtrue-color color data
%true-coloris not supported in painters renderer, so switch out of that
ifstrcmp(get(gcf,'renderer'), 'painters'),
set(gcf,'renderer','zbuffer');
end
%replaceoriginal CData with true-color data
if~strcmp(g.Type,'patch'),
set(hh,'CData',realcolor);
else
set(hh,'faceVertexCData',permute(realcolor,[1 3 2]))
end
%restoreclim (so colorbar will show correct limits)
if~isempty(parentAx),
set(parentAx,'clim',originalClim)
end
end %loop on indexed-color objects
%============================================================================%
% Local functions
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% getCDataHandles -- get handles of all descendents with indexedCData
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function hout = getCDataHandles(h)
% getCDataHandles Find all objects with indexedCData
%recursively descend object tree, finding objects with indexedCData
% An exception: don't include children of objects that themselveshave CData:
% for example, scattergroupsare non-standard hggroups, with CData. Changing
% such a group's CDataautomatically changes the CData of its children,
% (as well as the children'shandles), so there's no need to act on them.
error(nargchk(1,1,nargin,'struct'))
hout = [];
if isempty(h),return;end
ch = get(h,'children');
for hh = ch'
g =get(hh);
ifisfield(g,'CData'), %does object have CData?
%is it indexed/scaled?
if ~isempty(g.CData) &&isnumeric(g.CData) &&size(g.CData,3)==1,
hout = [hout; hh]; %#ok<AGROW> %yes,add to list
end
else %noCData, see if object has any interesting children
hout = [hout; getCDataHandles(hh)];%#ok<AGROW>
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% getParentAxes -- return handle of axes object to which a givenobject belongs
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function hAx = getParentAxes(h)
% getParentAxes Return enclosing axes of a givenobject (could be self)
error(nargchk(1,1,nargin,'struct'))
%object itself may be an axis
if strcmp(get(h,'type'),'axes'),
hAx =h;
return
end
parent = get(h,'parent');
if (strcmp(get(parent,'type'), 'axes')),
hAx =parent;
else
hAx =getParentAxes(parent);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% checkArgs -- Validate input arguments
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [h, nancolor] = checkArgs(args)
% checkArgs Validate input arguments tofreezeColors
nargs = length(args);
error(nargchk(0,3,nargs,'struct'))
%grab handle from first argument if we have an odd number ofarguments
if mod(nargs,2),
h =args{1};
if~ishandle(h),
error('JRI:freezeColors:checkArgs:invalidHandle',...
'The first argument must be a valid graphics handle (to anaxis)')
end
% 4/2010check if object to be frozen is a colorbar
ifstrcmp(get(h,'Tag'),'Colorbar'),
if ~exist('cbfreeze.m'),
warning('JRI:freezeColors:checkArgs:cannotFreezeColorbar',...
['You seem to be attempting to freeze a colorbar. This nolonger'...
'works. Please read the help for freezeColors for thesolution.'])
else
cbfreeze(h);
return
end
end
args{1} =[];
nargs =nargs-1;
else
h =gca;
end
%set nancolor if that option was specified
nancolor = [nan nan nan];
if nargs == 2,
ifstrcmpi(args{end-1},'nancolor'),
nancolor = args{end};
if ~all(size(nancolor)==[1 3]),
error('JRI:freezeColors:checkArgs:badColorArgument',...
'nancolor must be [r g b] vector');
end
nancolor(nancolor>1) = 1;nancolor(nancolor<0) = 0;
else
error('JRI:freezeColors:checkArgs:unrecognizedOption',...
'Unrecognized option (%s). Only ''nancolor'' isvalid.',args{end-1})
end
end
function CBH = cbfreeze(varargin)
�FREEZE Freezes the colormap ofa colorbar.
%
% SYNTAX:
% cbfreeze
% cbfreeze('off')
% cbfreeze(H,...)
% CBH = cbfreeze(...);
%
% INPUT:
% H - Handles of colorbars to be freezed, or from figures tosearch
% for them or from peer axes (see COLORBAR).
% DEFAULT: gcf (freezes all colorbars from the current figure)
% 'off' - Unfreezes the colorbars, other options are:
% 'on' Freezes
% 'un' same as'off'
% 'del' Deletes thecolormap(s).
% DEFAULT: 'on' (of course)
%
% OUTPUT (all optional):
% CBH - Color bar handle(s).
%
% DESCRIPTION:
% MATLAB works with a unique COLORMAP by figure which is a big
% limitation. Function FREEZECOLORS by John Iversen allows touse
% different COLORMAPs in a single figure, but it fails freezingthe
% COLORBAR. This program handles this problem.
%
% NOTE:
% * Optional inputs use its DEFAULT value when not given or [].
% * Optional outputs may or not be called.
% * If no colorbar is found, one is created.
% * The new frozen colorbar is an axes object and does notbehaves
% as normally colorbars when resizing the peer axes. Although,some
% time the normal behavior is not that good.
% * Besides, it does not have the 'Location' property anymore.
% * But, it does acts normally: no ZOOM, no PAN, no ROTATE3D andno
% mouse selectable.
% * No need to say that CAXIS and COLORMAP must be defined beforeusing
% this function. Besides, the colorbar location. Anyway, 'off'or
% 'del' may help.
% * The 'del' functionality may be used whether or not thecolorbar(s)
% is(are) froozen. The peer axes are resized back. Try:
% >> colorbar, cbfreeze del
%
% EXAMPLE:
% surf(peaks(30))
% colormap jet
% cbfreeze
% colormap gray
% title('What...?')
%
% SEE ALSO:
% COLORMAP, COLORBAR, CAXIS
% and
% FREEZECOLORS by John Iversen
% at http://www.mathworks.com/matlabcentral/fileexchange
%
%
% ---
% MFILE: cbfreeze.m
% VERSION: 1.1 (Sep 02, 2009)(<ahref="matlab:web('http://www.mathworks.com/matlabcentral/fileexchange/authors/11258')">download</a>)
% MATLAB: 7.7.0.471 (R2008b)
% AUTHOR: Carlos Adrian Vargas Aguilera (MEXICO)
% CONTACT: nubeobscura@hotmail.com
% REVISIONS:
% 1.0 Released. (Jun 08, 2009)
% 1.1 Fixed BUG with image handle on MATLAB R2009a. Thanks toSergio
% Muniz. (Sep 02, 2009)
% DISCLAIMER:
% cbfreeze.m is provided "asis" without warranty of any kind, under the
% revised BSD license.
% Copyright (c) 2009 CarlosAdrian Vargas Aguilera
% INPUTS CHECK-IN
%-------------------------------------------------------------------------
% Parameters:
cbappname ='Frozen'; % Colorbar application data with fields:
% 'Location' from colorbar
% 'Position' from peer axes befor colorbar
%'pax' handle from peer axes.
axappname = 'FrozenColorbar'; % Peer axes application data withfrozen
% colorbar handle.
% Set defaults:
S ='on'; Sopt = {'on','un','off','del'};
H = get(0,'CurrentFig');
% Check inputs:
if nargin==2 &&(~isempty(varargin{1}) &&all(ishandle(varargin{1})) &&...
isempty(varargin{2}))
% Check for CallBacks functionalities:
% ------------------------------------
varargin{1} = double(varargin{1});
ifstrcmp(get(varargin{1},'BeingDelete'),'on')
% Working as DeletFcn:
if (ishandle(get(varargin{1},'Parent'))&& ...
~strcmpi(get(get(varargin{1},'Parent'),'BeingDeleted'),'on'))
% The handleinput is being deleted so do the colorbar:
S ='del';
if~isempty(getappdata(varargin{1},cbappname))
% The frozencolorbar is being deleted:
H =varargin{1};
else
% The peeraxes is being deleted:
H =ancestor(varargin{1},{'figure','uipanel'});
end
else
% The figure is gettingclose:
return
end
elseif (gca==varargin{1}&& ...
gcf==ancestor(varargin{1},{'figure','uipanel'}))
% Working as ButtonDownFcn:
cbfreezedata =getappdata(varargin{1},cbappname);
if ~isempty(cbfreezedata)
ifishandle(cbfreezedata.ax)
% Turns thepeer axes as current (ignores mouse click-over):
set(gcf,'CurrentAxes',cbfreezedata.ax);
return
end
else
% Clears applicationdata:
rmappdata(varargin{1},cbappname)
end
H = varargin{1};
end
else
% Checks for normal calling:
% --------------------------
% Looks for H:
if nargin &&~isempty(varargin{1}) &&all(ishandle(varargin{1}))
H = varargin{1};
varargin(1) = [];
end
% Looks for S:
if ~isempty(varargin)&& (isempty(varargin{1}) ||ischar(varargin{1}))
S = varargin{1};
end
end
% Checks S:
if isempty(S)
S = 'on';
end
S = lower(S);
iS = strmatch(S,Sopt);
if isempty(iS)
error('CVARGAS:cbfreeze:IncorrectStringOption',...
['Unrecognized ''' S ''' argument.' ])
else
S = Sopt{iS};
end
% Looks for CBH:
CBH = cbfreeze(H); �H = cbhandle(H);
if ~strcmp(S,'del') &&isempty(CBH)
% Creates a colorbar and peer axes:
pax = gca;
CBH = colorbar('peer',pax);
else
pax = [];
end
%-------------------------------------------------------------------------
% MAIN
%-------------------------------------------------------------------------
% Note: only CBH and S are necesary, but I use pax to avoid the useof the
% "hidden" 'Axes' COLORBAR's property. Why... �
% Saves current position:
fig = get( 0,'CurrentFigure');
cax = get(fig,'CurrentAxes');
% Works on every colorbar:
for icb = 1:length(CBH)
% Colorbar axes handle:
h = double(CBH(icb));
% This application data:
cbfreezedata = getappdata(h,cbappname);
% Gets peer axes:
if ~isempty(cbfreezedata)
pax = cbfreezedata.pax;
if ~ishandle(pax) % just in case
rmappdata(h,cbappname)
continue
end
elseif isempty(pax) % not generated
try
pax =double(get(h,'Axes')); % NEW feature inCOLORBARs
catch
continue
end
end
% Choose functionality:
switch S
case 'del'
% Deletes:
if~isempty(cbfreezedata)
% Returnsaxes to previous size:
oldunits =get(pax,'Units');
set(pax,'Units','Normalized');
set(pax,'Position',cbfreezedata.Position)
set(pax,'Units',oldunits)
set(pax,'DeleteFcn','')
ifisappdata(pax,axappname)
rmappdata(pax,axappname)
end
end
ifstrcmp(get(h,'BeingDelete'),'off')
delete(h)
end
case {'un','off'}
% Unfrozes:
if~isempty(cbfreezedata)
delete(h);
set(pax,'DeleteFcn','')
ifisappdata(pax,axappname)
rmappdata(pax,axappname)
end
oldunits =get(pax,'Units');
set(pax,'Units','Normalized')
set(pax,'Position',cbfreezedata.Position)
set(pax,'Units',oldunits)
CBH(icb) =colorbar(...
'peer' ,pax,...
'Location',cbfreezedata.Location);
end
otherwise % 'on'
% Freezes:
% Gets colorbar axesproperties:
cb_prop =get(h);
% Gets colorbar image handle.Fixed BUG, Sep 2009
hi =findobj(h,'Type','image');
% Gets image data andtransform it in a RGB:
CData = get(hi,'CData');
if size(CData,3)~=1
% It'salready frozen:
continue
end
% Gets image tag:
Tag = get(hi,'Tag');
% Deletes previous colorbarpreserving peer axes position:
oldunits =get(pax,'Units');
set(pax,'Units','Normalized')
Position =get(pax,'Position');
delete(h)
cbfreezedata.Position =get(pax,'Position');
set(pax,'Position',Position)
set(pax,'Units',oldunits)
% Generates new colorbaraxes:
% NOTE: this is needed becauseeach time COLORMAP or CAXIS is used,
% MATLAB generates a new COLORBAR! This eliminates thatbehaviour
% and is the central point on this function.
h = axes(...
'Parent' ,cb_prop.Parent,...
'Units' ,'Normalized',...
'Position',cb_prop.Position...
);
% Save location for futurecall:
cbfreezedata.Location =cb_prop.Location;
% Move ticks because IMAGEdraws centered pixels:
XLim = cb_prop.XLim;
YLim = cb_prop.YLim;
if isempty(cb_prop.XTick)
%Vertical:
X = XLim(1)+ diff(XLim)/2;
Y =YLim +diff(YLim)/(2*length(CData))*[+1 -1];
else % isempty(YTick)
%Horizontal:
Y = YLim(1)+ diff(YLim)/2;
X =XLim +diff(XLim)/(2*length(CData))*[+1 -1];
end
% Draws a new RGB image:
image(X,Y,ind2rgb(CData,colormap),...
'Parent' ,h,...
'HitTest' ,'off',...
'Interruptible' ,'off',...
'SelectionHighlight','off',...
'Tag' ,Tag...
)
% Removesall '...Mode' properties:
cb_fields =fieldnames(cb_prop);
indmode =strfind(cb_fields,'Mode');
for k=1:length(indmode)
if~isempty(indmode{k})
cb_prop = rmfield(cb_prop,cb_fields{k});
end
end
% Removes special COLORBARsproperties:
cb_prop =rmfield(cb_prop,{...
'CurrentPoint','TightInset','BeingDeleted','Type',... % read-only
'Title','XLabel','YLabel','ZLabel','Parent','Children',... % handles
'UIContextMenu','Location',... % colorbars
'ButtonDownFcn','DeleteFcn',... % callbacks
'CameraPosition','CameraTarget','CameraUpVector','CameraViewAngle',...
'PlotBoxAspectRatio','DataAspectRatio','Position',...
'XLim','YLim','ZLim'});
% And now, set new axesproperties almost equal to the unfrozen
% colorbar:
set(h,cb_prop)
% CallBack features:
set(h,...
'ActivePositionProperty','position',...
'ButtonDownFcn' ,@cbfreeze,... % mhh...
'DeleteFcn' ,@cbfreeze) % again
set(pax,'DeleteFcn' ,@cbfreeze) % and again!
% Do not zoom or pan orrotate:
setAllowAxesZoom (zoom ,h,false)
setAllowAxesPan (pan ,h,false)
setAllowAxesRotate(rotate3d,h,false)
% Updates data:
CBH(icb) =h;
% Saves data for futureundo:
cbfreezedata.pax = pax;
setappdata( h,cbappname,cbfreezedata);
setappdata(pax,axappname,h);
end % switchfunctionality
end % MAIN loop
% OUTPUTS CHECK-OUT
%-------------------------------------------------------------------------
% Output?:
if ~nargout
clear CBH
else
CBH(~ishandle(CBH)) = [];
end
% Returns current axes:
if ishandle(cax)
set(fig,'CurrentAxes',cax)
end
% [EOF] cbfreeze.m
大家只要把这个.m文件matlab工作文件夹中即可,此函数用来在一个figure上作出多个立体图形。
e.g.
subplot(1,2,1)
sphere(10)
colorbar
freezeColors
subplot(1,2,2)
peaks(20)
colorbar
cbfreeze
图形如下:
这样在同一个figure中就显示出两个colorbar。
()