Error using imwarp - ippgeotrans requestest a huge array

While using imwarp on an 777x3161 image I get the following error:

Error using ippgeotrans
Requested 50995x174010 (33.1GB) array exceeds maximum array size preference. Creation of arrays greater than this limit may take a long time and
cause MATLAB to become unresponsive. See array size limit or preference panel for more information.
Error in imwarp>ippWarpAffine (line 566)
    outputImage = ippgeotrans(single(inputImage),double(T),Rout.ImageSize,interp,single(fillVal));
Error in imwarp>remapPointsAndResample (line 223)
        outputImage = ippWarpAffine(inputImage,R_A,tform,outputRef,method,fillValues);
Error in imwarp (line 210)
outputImage = remapPointsAndResample(parsedInputs.InputImage,R_A,tform,outputRef,method,fillValues);

I use imwarp inside a loop and the error only occurs for one particular image. For the other images imwarp returns an image with the expected size (777x3161)

There is no documentation on ippgeotrans so I have no idea why it is requesting an array of this astronomical size... Would it help if I perhaps chose another interpolation method?

Can anyone help/explain what is going on here?

Best regards

Silja

Answer 1

[Walter Roberson]

Is it possible that the one image is 3 dimensional when the others are all 2 dimensional?

At the command line, issue the command

dbstop if error

and run the program for that one problem image. When it stops in the error, ask it to show size(inputImage) . Then give the command dbup and ask for the size again. Keep using dbup and asking for the size until you get to remapPointsAndResample and in there look at size(parsedInputs.InputImage)

[Silja Heilmann]

Dear Walter Thank you for your answer. The problem image in definitely not 3-dim (I checked).

I tried following your instructions, the output I got is cp'ed below:

>> dbstop if error

ans =
         777        3161

(this was me asking for the size of the problem image before passing it to imwarp!)

Error using ippgeotrans
Requested 50995x174010 (33.1GB) array exceeds maximum array size preference. Creation of arrays greater than this limit may take a long time and cause MATLAB to become unresponsive. See array size limit or preference panel for more information.

Error in imwarp>ippWarpAffine (line 566)
    outputImage = ippgeotrans(single(inputImage),double(T),Rout.ImageSize,interp,single(fillVal));

Error in imwarp>remapPointsAndResample (line 223)
        outputImage = ippWarpAffine(inputImage,R_A,tform,outputRef,method,fillValues);

Error in imwarp (line 210)
outputImage = remapPointsAndResample(parsedInputs.InputImage,R_A,tform,outputRef,method,fillValues);

566     outputImage = ippgeotrans(single(inputImage),double(T),Rout.ImageSize,interp,single(fillVal));

Now I start asking for the size of "inputImage"

K>> size(inputImage) 

ans =

           781        3165

K>> dbup
In workspace belonging to imwarp>remapPointsAndResample (line 223)
K>> size(inputImage)

ans =

           777        3161

K>> dbup
In workspace belonging to imwarp (line 210)
K>> size(inputImage)
Undefined function or variable 'inputImage'.

K>> size(parsedInputs.InputImage)

ans =

           777        3161

K>> dbup
In base workspace.

Hope that means something to you? (I tried with a different interpolation method and got the same error so that is not it). Also I got the error for one other image in my batch. But most of them run with no problem!

[Walter Roberson]

Sorry, I do not see how those input sizes lead to attempting to allocate that size of array. I would need access to the toolbox code to follow it.

[Silja Heilmann]

Hey Walter I just found out that when I generate the transform using 'pure' images instead of pre-processing them using imhistmatch() I do not get the error!

The code was:

fixed = timestack_BF(:,:,1,tt);
moving = imhistmatch(timestack_BF(:,:,1,tt+1),fixed);

[optimizer,metric] = imregconfig('monomodal');
tform = imregtform(moving,fixed,'similarity',optimizer,metric);  

moving_reg = imwarp(moving,tform,'nearest');

(I got the error in the last line).

This code however:

fixed = timestack_BF(:,:,1,tt);
moving = timestack_BF(:,:,1,tt+1);

[optimizer,metric] = imregconfig('monomodal');
tform = imregtform(moving,fixed,'similarity',optimizer,metric);  

moving_reg = imwarp(moving,tform,'nearest');

runs without error for the two problem images!