R
Ray Schumacher
Hi all,
I'm trying to port Hongjie Xie's IDL/ENVI Implementation of the FFT Based
Algorithm for Automatic Image Registration
ftp://ftp.iamg.org/VOL29/v29-08-10.zip
http://www.nmt.edu/~hjxie/xie-paper.pdf
http://www.iamg.org/CGEditor/cg2003.htm
to Python/numarray/Numeric and running into some subtle problems.
I don't have IDL, so I could only read the PDF docs.
For instance, the IDL ROT() requires some additional steps to emulate, and
the MAX(array, position) as well.
And, I'm not at all familiar with the correspondence of the FFT
libraries... even the simple shift-only routine
http://rjs.org/astro/1004x/Python/register/shift_idl.py
produces different (incorrect) x shift values than the IDL.
My current code
http://rjs.org/astro/1004x/Python/register/
conforms to Xie's flow as much as possible for now, but should and will be
re-factored when functional. His LogPolar function, to convert rectangular
coordinate arrays to polar, is particularly slow.
If anyone has some insight as to possible, critical differences, please let
me know!
Thank you,
Ray Schumacher
user:ray at the domain rjs daaat oh-arr-gee
I'm trying to port Hongjie Xie's IDL/ENVI Implementation of the FFT Based
Algorithm for Automatic Image Registration
ftp://ftp.iamg.org/VOL29/v29-08-10.zip
http://www.nmt.edu/~hjxie/xie-paper.pdf
http://www.iamg.org/CGEditor/cg2003.htm
to Python/numarray/Numeric and running into some subtle problems.
I don't have IDL, so I could only read the PDF docs.
For instance, the IDL ROT() requires some additional steps to emulate, and
the MAX(array, position) as well.
And, I'm not at all familiar with the correspondence of the FFT
libraries... even the simple shift-only routine
http://rjs.org/astro/1004x/Python/register/shift_idl.py
produces different (incorrect) x shift values than the IDL.
My current code
http://rjs.org/astro/1004x/Python/register/
conforms to Xie's flow as much as possible for now, but should and will be
re-factored when functional. His LogPolar function, to convert rectangular
coordinate arrays to polar, is particularly slow.
If anyone has some insight as to possible, critical differences, please let
me know!
Thank you,
Ray Schumacher
user:ray at the domain rjs daaat oh-arr-gee