comp.lang.idl-pvwave archive: archive » Re: fft and least sqaues problem

Home » Public Forums » archive » Re: fft and least sqaues problem

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: fft and least sqaues problem [message #81141 is a reply to message #81135]

Tue, 14 August 2012 11:22

rogass
Messages: 200
Registered: April 2008

Senior Member

Finally, I found the mistake - sitting in front of this pc :)

Here is the code:

s=size(im,/dim)*1.
im1 =im
dx=.25
dy=.3
im2=image_shift(im1,dx,dy);more precise then interpolate(im,findgen(s[0])+dx,findgen(s[1])+dy,/grid,/cubi c)
fim1=fft(im1,-1)
fim2=fft(im2,-1)
corr=fim1*conj(fim2)/abs(fim1*fim2)
corr=shift(corr,s/2)
phase=atan(imaginary(corr)/real_part(corr))
wx=(findgen(s) mod s[0])
wy=(rebin(findgen(1,s[1]),s))
wx-=s[0]/2
wy-=s[1]/2
wx*=2.*!pi/s[0]
wy*=2.*!pi/s[1]
r=5;fitting radius
phase2=phase[s[0]/2 - r : s[0]/2 + r,s[1]/2 - r : s[1]/2 + r]
wx2=wx[s[0]/2 - r : s[0]/2 + r,s[1]/2 - r : s[1]/2 + r]
wy2=wy[s[0]/2 - r : s[0]/2 + r,s[1]/2 - r : s[1]/2 + r]
print,la_least_squares(transpose([[wx2[*]],[wy2[*]]]),phase2 [*])

Thanks Helder!

Cheers

CR

Report message to a moderator

[Message index]

		Re: fft and least sqaues problem By: rogass on Fri, 17 August 2012 07:40
		Re: fft and least sqaues problem By: Craig Markwardt on Thu, 16 August 2012 21:58
		Re: fft and least sqaues problem By: rogass on Thu, 16 August 2012 02:13
		Re: fft and least sqaues problem By: Craig Markwardt on Tue, 14 August 2012 16:53
		Re: fft and least sqaues problem By: rogass on Tue, 14 August 2012 11:22
		Re: fft and least sqaues problem By: Helder Marchetto on Tue, 14 August 2012 07:50
		Re: fft and least sqaues problem By: rogass on Tue, 14 August 2012 07:13
		Re: fft and least sqaues problem By: Craig Markwardt on Fri, 17 August 2012 08:45

Previous Topic:	Coyote Library Updates
Next Topic:	RESOLVE_ROUTINE not finding MEDIAN function?

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sun Oct 12 12:09:35 PDT 2025

Total time taken to generate the page: 1.52017 seconds