1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Fast Fourier Transform (FFT) power spectrum angle

  1. Mar 18, 2014 #1
    Dear Physics Buddies,

    How are well all, okay I hope. I was wondering if I might browse all your infinite intellects and ask you a very simple question.

    I am working with some medical images in matlab and my collaborators would like to know the orientation of the fibres that it contains.

    So far I have:

    Tiled the images (with an increasing tile size of the power of 2)
    Calculated the regions ft

    Then I was going to use PCA to find the largest eigen vector and declare that as the most dominant orientation for this region. However, I am a little stuck.

    Someone suggested using the interia matrix of the ft power spectrum. This involves computing:


    I am using brocploc in matlab and so I was wondering if any of you guys new of quick ways to compute this for image?

    Hoping you can help

  2. jcsd
  3. Mar 25, 2014 #2


    User Avatar
    Science Advisor

    I guess the method used would be highly dependent on the patterns in the image and the colour / grey scale format. Are you trying to measure the principal axis in each subsections of the image or over the whole of the image? The fibres must be significantly more than 2 pixels wide to process the image.

    I would consider a quick and dirty arithmetic technique such as sampling the circumference of a circle described on each subsection. By accumulating the absolute differences of sequential pixel values around the circle and saving all the values you will have a ramp with variable slope. Now subtract a linear fit to put the final value at zero. The resulting plot will have two peaks and two troughs. The positive slope will be where the circle is crossing the fibres, the negative slope will be where the fibres are tangent. Fit a 2nd harmonic to the plot, the phase reveals the principle direction of fibres.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook