comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Butterworth Band-Pass Filter
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: Butterworth Band-Pass Filter [message #73914 is a reply to message #73913] Sun, 12 December 2010 08:56 Go to previous messageGo to previous message
David Fanning is currently offline  David Fanning
Messages: 11724
Registered: August 2001
Senior Member
burton449 writes:

> According to your documentation, a high-pass butterworth filter will
> be defined like that:
> filter = 1 / [1 + C(Ro/R)^2n]
>
> So in IDL, for a cuttoff of 15, and freqImage = FFT(image, -1), the
> filter will be defined like that:
> filter = 1.0 / ( 1.0d + (15.0/freqImage)^2 )
>
> Is it right?

No, what you are calling freqImage is called freqDomainImage
in that example. What you want is something like this:

s = Size(image, /Dimensions)
filter = 1.0D / (1.0D + 15.0/Dist(s[0],s[1]))^2)

In other words, what you are calling "freqImage" in
the filter is actually built with the DIST function,
and is different from the freqDomainImage, which is
defined like this:

freqDomainImage = FFT(image, -1)

Cheers,

David



--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
Sepore ma de ni thui. ("Perhaps thou speakest truth.")
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Over-PLOT() on MAP()?
Next Topic: Run envi routine in IDL Virtual MAchine

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Thu Oct 09 22:54:20 PDT 2025

Total time taken to generate the page: 0.23489 seconds