用于Matlab的图像处理库.zip

classdef Slicer < handle %GUI for exploration of 3D images, using Object Oriented Programming. % % SLICER is an graphical interface to explore slices of a 3D image. % Index of the current slice is given under the slider, mouse position as % well as cursor value are indicated when mouse is moved over image, and % scrollbars allow to navigate within image. % % SLICER should work with any kind of 3D images: binary, gray scale % (integer or floating-point) or color RGB. % % slicer(IMG); % where IMG is a preloaded M*N*P matrix, opens the slicer GUI, % initialized with image IMG. % User can change current slice with the slider to the left, X and Y % position with the two corresponding sliders, and change the zoom in the % View menu. % % Slicer(IMGNAME, ...); % Load the stack specified by IMGNAME. It can be either a tif bundle, the % first file of a series, or a 3D image stored in one of the medical % image format: % * DICOM (*.dcm) % * Analyze (*.hdr) % * MetaImage (*.mhd, *.mha) % It is also possible to import a raw data file, from the File->Import % menu. % % Slicer % without argument opens with an empty interface, and allows to open an % image through the menu. % % Slicer(..., PARAM, VALUE); % Specifies one or more display options as name-value parameter pairs. % Available parameter names are: % % * 'Slice' the display uses slice given by VALUE as current slice % % * 'Name' gives a name to the image (for display in title bar) % % * 'Spacing' specifies the size of voxel elements. VALUE is a 1-by-3 % row vector containing spacing in x, y and z direction. % % * 'Origin' specifies coordinate of first voxel in user space % % * 'UnitName' the name of the unit used for spatial calibration (can % be any string, default is empty). % % * 'DisplayRange' the values of min and max gray values to display. The % default behaviour is to use [0 255] for uint8 images, or to % compute bounds such as 95% of the voxels are converted to visible % gray levels for other image types. % % * 'ColorMap' The colormap used for displaying grayscale images % (default is gray). Should be a N-by-3 array of double, with N=256 % for grayscale or intensity images, as N=the number of labels for % label images, and N=2 for binary images. % % * 'BackgroundColor' the color used as background for label images. % % * 'ImageType' The type of image, used for adapting display. Can be % one of 'binary', 'grayscale', 'intensity', 'label', 'color', % 'vector', 'none'. Default value is assessed from data type and % size. % % * 'Parent' another instance of Slicer, that is used to initialize % several parameters like spatial resolution, display range, LUT % for display... % % % Example % % Explore human brain MRI % metadata = analyze75info('brainMRI.hdr'); % I = analyze75read(metadata); % Slicer(I); % % % show the 10-th slice, and add some setup % Slicer(I, 'Slice', 10, 'Spacing', [1 1 2.5], 'Name', 'Brain', 'DisplayRange', [0 90]); % % See also % imStacks, imscrollpanel % % Requires % GUI Layout Toolbox version 2.0, or 1.17 (try to switch depending on % Matlab version) % % ------ % Author: David Legland % e-mail: david.legland@inrae.fr % Created: 2011-04-12, using Matlab 7.9.0.529 (R2009b) % https://github.com/mattools/matImage % Copyright 2011 INRAE - BIA - BIBS %% Properties properties % Image data stored as a 3D or 4D array, in YX(C)Z order. ImageData; % type of image. Can be one of: % 'binary' % 'grayscale' % 'intensity' % 'label' % 'color' % 'vector' % 'none' ImageType; % size of the reference image (1-by-3 row vector, in XYZ order) ImageSize; % extra info for image, such as the result of imfinfo ImageInfo; % extra info for image, such as the result of imfinfo ImageName; % displayed data (2D image) Slice; % z-index of slice within 3D image SliceIndex; % used to adjust constrast of the slice DisplayRange; % Look-up table for display of uint8, label and binary images. % can be empty, in obj case gray colormap is assumed ColorMap = []; % background color for label to RGB conversion. Given as RGB triplet of % values bewteen 0 and 1. Default is white = (1,1,1). BgColor = [1 1 1]; % calibration information for image VoxelOrigin; VoxelSize; VoxelSizeUnit; % shortcut for avoiding many tests. Should be set to true when either % voxelOrigin, voxelsize or voxelSizeUnit is different from its default % value. Calibrated = false; % keep last path for opening new images LastPath = pwd; % list of handles to the widgets. % Structure with following fields: % * Figure: the main figure % * SubFigures: a list of figures that can be closed when the main % figure is closed % * Image: the image display % * ImageAxis: the axis containing the image % * ZSlider: the slider used to change the slice index % * ZEdit: the edit uicontrol used to change slice index % * ZProfileFigure used to display the profile along z % * ZProfileAxis used to display the profile along z Handles; % the last position of mouse click, in user coordinates. LastClickedPoint = []; end %% Constructor methods function obj = Slicer(varargin) % call parent constructor obj = obj@handle(); % initialize structure containing handles obj.Handles = struct(); % keep pointer to current path obj.LastPath = pwd; % initialize using image given as argument if ~isempty(varargin) var = varargin{1}; if isa(var, 'Image') setupImageFromClass(obj, var); elseif ischar(var) setupImageFromFile(obj, var); else if ndims(var) < 3 %#ok<ISMAT> error('Requires input image with at least 3 dimensions'); end setupImageData(obj, var, inputname(1)); end varargin(1) = []; if isempty(obj.ImageData) return; end else obj.ImageData = []; obj.ImageType = 'none'; end % parses input arguments, given as list of name-value pairs parsesInputArguments(); % add checkup on visible image slice if ~isempty(obj.ImageData) obj.SliceIndex = min(obj.SliceIndex, obj.ImageSize(3)); end updateCalibrationFlag(obj); % create default figure fig = figure(); set(fig, 'Menubar', 'none'); set(fig, 'NumberTitle', 'off'); set(fig, 'Name', 'Slicer'); obj.Handles.Figure = fig; obj.Handles.SubFigures = []; % create main figure menu and layout setupMenuBar(obj); setupLayout(fig); % setup new image display if ismember(obj.ImageType, {'label', 'binary'}) maxi = max(obj.ImageData(:)); obj.DisplayRange = [0 maxi]; if isempty(obj.ColorMap) obj.ColorMap = jet(double(maxi)); end colormap([obj.BgColor ; obj.ColorMap]); end updateSlice(obj); displayNewImage(obj); updateTitle(obj); % setup listeners