MATLAB: Error message in Intensity-Based Automatic Image Registration

Question

Hi

I am trying to compare and analyze blank and dermatome maps (attached) in order to quantify number of shaded pixels and percentage of shaded pixels within the template body outline. Ultimately this will be used with pain and paresthesia maps. My understanding is that I wil want to use Image Registration in the Image processing toolbox, and that 'Intensity-Based Automatic Image Registration' would serve my purposes. While following the online tutorial I keep running into this error when I try to use imRegister:

Error using imregtform>parseInputs (line 252) The value of 'MovingImage' is invalid. All dimensions of the moving image should be greater than 4.

Error in imregtform (line 123) parsedInputs = parseInputs(varargin{:});

Error in imregister (line 119) tform = imregtform(varargin{:});

Any help on how to solve this problem or any general information on image registration would be greatly appreciated. My entire coding attempt is below for your viewing pleasure:

>> fixed = imread ('blankmap.jpg');

>> moving = imread ('dermatomemap.jpg');

>> figure, imshowpair (moving, fixed, 'montage')

Warning: Image is too big to fit on screen; displaying at 13% > In imuitools/private/initSize at 71 In imshow at 282 In imshowpair at 124

>> title ('unregistered')

>> [optimizer, metric] = imregconfig ('multimodal'); >> movingRegisteredDefault = imregister (moving, fixed, 'affine', optimizer, metric);

Error using imregtform>parseInputs (line 252) The value of 'MovingImage' is invalid. All dimensions of the moving image should be greater than 4.

Error in imregtform (line 123) parsedInputs = parseInputs(varargin{:});

Error in imregister (line 119) tform = imregtform(varargin{:});

Matt

Answer 1

It might be that you need to use grayscale images (using rgb2gray). Eg:

movingRegisteredDefault = imregister (rgb2gray(moving), rgb2gray(fixed), 'affine', optimizer, metric);

or you can register using a single channel :

movingRegisteredDefault = imregister (moving(:,:,Channel), fixed(:,:,Channel), 'affine', optimizer, metric);