Fourier-Mellin based Image Registration (with GUI)

GUI-implementation of the Fourier-Mellin transform for simple image rotation, scale and translation.

Description

This GUI implementation of the Fourier-Mellin transform allows registration of an image with respect to a base image. Translation, rotation and scale shift values are retrieved and used to form a registered reconstruction.

• Required Products: Image Processing Toolbox
• MATLAB release: MATLAB 6.0 (R12)
• HTML Files: Fourier-Mellin Transform GUI Help

Other Files

imlogpolar.m,
lena2.bmp,
lena_org.bmp,
parse_inputs.m,
readpgm.m,
transformImage.m,
window1d.m,
window2d.m,
zeropad.m,
hipass_filter.m,
crosspowerspectrum.m,
fft_resize_test.m,
fm_gui.m,
fm_guifun.m,
fm_parse_inputs.m,
fourier_mellin.m,
image_scale.m

[Gang Pan]

Doesn't work with my images

[Adam Wilmer]

How about you tell me what happens and I try and sort it out for you. Use the CONTACT AUTHOR bit above. Cheers. Adam.

[neel bhattacharya]

excellent help from adam

[Sic Volo]

Good stuff. True, it does not work for non square images, but still it was very helpful

[Kishore Mogatadakala]

Eccellet work!

[a c]

Excellent ,Thanks

[pavan kumar]

Its a very useful article

[Romeil Sandhu]

Hey, I couldn't run this program because it gave the following error when trying to register the images:

Change directory to the location of fm_gui.m...
...or change the value of FM_PATH in fm_gui.m.
Cannot rotate/scale a non-existent image
PROBLEM ALERT: phasecorrelation contains a lot of NaNs (check FFTs exist as they cannot cope with certain inputs for some reason)
??? Index exceeds matrix dimensions.

Error in ==> C:\MATLAB6p5\work\fm_gui_v2\fourier_mellin.m
On line 56 ==> sorted_scale = scale_neighbourhood(2);

Error in ==> C:\MATLAB6p5\work\fm_gui_v2\fm_guifun.m
On line 337 ==> [data.registered,data.input1registered,data.input2registered,reg_output,data.pc_rs,data.pc_trans] = fourier_mellin(data);

??? Error while evaluating uicontrol Callback.

Can someone help me please!

[pau soler]

very good

[Tomer Ish shalom]

Excelent work

However, The resampling algorithms can be improved and instead of using matlab standard resampling routines.
The routines suggested by yeroslavsky in "digital holography" can be used.

[vivek reddy]

 

[Robin Tucker]

Any hints on how to convert this to work with Matlab 7?

[Ji-Young Lim]

The transformImage.m has a few small bugs that do not let the following FFT be successful. GUI is wonderful!

[klesk chonkin]

Excellent job!

[Moti Freiman]

Is there any way to use this code for generating fourier mellin rotation, translation, scale invariant descriptors?

[Andrei Luk]

 

[shirzad Kalhori]

Hello !
Very nice page.
Best regards

[li haichao]

 

[bereket haile]

helo I need matlab codes for log-polar transform

[Filmon Abraham]

 

[rajavel pandeeswaran]

Really good

[Riber Kaban]

 

[Lin Zhang]

I have found out that the Fourier-Mellin transform (FMT) can only do image registration for two images with the same size and the same background. The FMT cannot do the template matching where the pattern template only contains an object in the input image and it is much smaller than the input image, but the input image contains many other things in addition to the object contained in the template.
In this case, I have found out that the FMT can not do the template matching to locate the pattern in the input image.

I would like to know if anyone else has the similar problem and how to solve this problem.

Thanks.

[wan mariana]

 

[Muharram Mansoorizadeh]

Dear Zhang, this is not surprising , because the FM method is built upon a set of transforms that apply on image structures(including background and objects).
To do template matching, in some cases, you can apply preprocessing to eliminate or soften background. Example of such cases are removing constant backgrounds, filtering out certain colors(e.g. skin color), morphological filtering and segmentation.

[Antai Guo]

I have tested the FMT Gui, while it seems that the matlab application doesn't work well, and when processing 'lena.bmp', it dumpped to dummy.

[Rudy Slingerland]

This program doesn't work, even with the test image provided in the zip file. Upon attempting to lad the first image it returns an error:
??? Undefined function or method 'isind' for input arguments of type 'uint8'.

Error in ==> fm_guifun at 77
            if isind(M1) & ~isempty(ma)

??? Error using ==> fm_guifun loadA; fm_guifun RotateScaleCropInput2; fm_guifun input1Select; fm_guifun input2Select;
Undefined function or method 'isind' for input arguments of type 'uint8'.

??? Error while evaluating uicontrol Callback

For debugging purposes, please note that I am using Matlab v. 7.4 (R2007a)

Thanks,

R Slingerland

[Christian Recknagel]

I have found an error in the subfunction "transformImage.m". In line 46,47,48,61 you want to set specific values to black. The value for black is zero and not NAN!!!! If you set the values to NAN the fft2 function used in "crosspowerspectrum.m" does not work.

[Sohan Ranjan]

The zip file is missing. Can anyone share?

[Greg Baxter]

While the GUI appears to function as designed, there were numerous errors in the underlying code.

NOTE: I'm using R2008A

??? Error using ==> get
There is no 'zdata' property in the 'image' class.

Error in ==> objbounds at 65
        zd = get(h(i), 'zdata');

Error in ==> axis>LocSetTight at 266
limits = objbounds(findall(ax));

Error in ==> axis>LocSetImage at 393
LocSetTight(ax);

Error in ==> axis at 139
                LocSetImage(ax(j),pbarlimit);

Error in ==> fm_guifun at 448
    axis image; % fit axis box tightly around image

??? Error while evaluating uicontrol Callback

Warning: isind is obsolete and may be removed in the future.
See product release notes for more information.
> In isind at 32
  In fm_guifun at 126
PROBLEM ALERT: phasecorrelation contains a lot of NaNs (check FFTs exist as they cannot cope with certain inputs for some reason)
??? Attempted to access scale_neighbourhood(2); index out of bounds because
numel(scale_neighbourhood)=0.

Error in ==> fourier_mellin at 56
    sorted_scale = scale_neighbourhood(2);

Error in ==> fm_guifun at 337
        [data.registered,data.input1registered,data.input2registered,reg_output,data.pc_rs,data.pc_trans]
        = fourier_mellin(data);

??? Error while evaluating uicontrol Callback

[JIAO …]

I found that it can work for the common images, but it can not work for the astronomic image.And anybody else can make a file without GUI?