# Difference between revisions of "Wavenumber"

(added stub content) |
(→2 Open your image: remark about (c)) |
||

(14 intermediate revisions by 2 users not shown) | |||

Line 1: | Line 1: | ||

− | Wavenumber, k, sometimes called the ''propagation number'', is a quantification of spatial scale. It can be thought of as a spatial analog to the temporal [[frequency]]. It is often defined as the number of | + | {{A2Z}} |

+ | Wavenumber, k, sometimes called the ''propagation number'', is a quantification of spatial scale. It can be thought of as a spatial analog to the temporal [[frequency]], and is often called ''spatial frequency''. It is often defined as the number of wavelengths per unit distance, or in terms of wavelength, ''λ'': | ||

:<math>k = \frac{1}{\lambda}\ </math> | :<math>k = \frac{1}{\lambda}\ </math> | ||

− | In | + | This is analogous to frequency ''f'', which is the reciprocal of period ''T''; that is, ''f'' = 1/''T''. In a sense, period can be thought of as a temporal 'wavelength—the length of an oscillation in time. |

+ | |||

+ | ==Angular wavenumber== | ||

+ | If you've explored the applications of frequency in geophysics, you'll have noticed that we sometimes don't use ordinary frequency f, in Hertz. Because geophysics deals with oscillating waveforms, ones that vary around a central value (think of a wiggle trace of seismic data), we often use the '''circular''' or '''angular frequency'''. In this way, we can also express the close relationship between frequency and phase, which is an angle. So in many geophysical applications, we want the ''angular wavenumber''. It is expressed in radians per metre: | ||

:<math>k = \frac{2\pi}{\lambda}\ </math> | :<math>k = \frac{2\pi}{\lambda}\ </math> | ||

+ | |||

+ | The relationship between angular wavenumber and angular frequency is analogous to that between wavelength and ordinary frequency—they are related by the velocity ''V'': | ||

+ | |||

+ | :<math>k = \frac{\omega}{V}\ </math> | ||

+ | |||

+ | It's unfortunate that there are two definitions of ''wavenumber''. Some people reserve the term ''spatial frequency'' for the ordinary wavenumber, or use ν (that's a Greek ''nu'', not a vee — another potential source of confusion!), or even ''σ'' for it. But just as many call it the wavenumber and use ''k'', so the only sure way through the jargon is to specify what you mean by the terms you use. As usual! | ||

+ | |||

+ | Here are two spectra: the FFTs of an image of some waves, and a binary image of some particles: | ||

+ | |||

+ | [[file:waves-particles_FFT.png|600px]] | ||

+ | |||

+ | ==Computing 2D Fourier transform in FIJI== | ||

+ | You can use FIJI to compute the 2D FFT of an image. FIJI is a popular open-source, scriptable scientific image-processing tool. '''F'''IJI '''I'''s '''J'''ust (a more awesome) [http://rsbweb.nih.gov/ij/ '''I'''mageJ]. | ||

+ | |||

+ | ===1 Download and install FIJI=== | ||

+ | [http://fiji.sc/wiki/index.php/Fiji Get the software from here.] | ||

+ | |||

+ | ===2 Open your image=== | ||

+ | This image is a low-res version of a photomicrograph of some shale from Schieber et al (2010)<ref>Juergen Schieber, John B. Southard, and Arndt Schimmelmann (2010). Lenticular shale fabrics resulting from intermittent erosion of water-rich muds—interpreting the rock record in the light of recent flume experiments. ''Journal of Sedimentary Research'' '''80''' (1), p 119–128. DOI [http://dx.doi.org/10.2110/jsr.2010.005 10.2110/jsr.2010.005].</ref>. | ||

+ | |||

+ | [[File:FIJI_FFT_Image.png|400px]] | ||

+ | |||

+ | '''Please note:''' This image is not covered by the CC-BY license on the rest of this wiki. It is copyright of the authors. I'm calling the use of a low-res version 'fair use'. I will try to build an example with an open image at some point — feel free to change it! | ||

+ | |||

+ | ===3 Set the scale=== | ||

+ | No matter what you're doing in FIJI, it's often a good idea to tell FIJI the real-world units of your image, so you don't have to think in pixels. To do this, use the rectangle selection tool to measure a length in pixels (e.g. using a scale in the image), then select '''Analyze > Set scale...''' | ||

+ | |||

+ | [[File:FIJI_Set_scale_example.png|800px]] | ||

+ | |||

+ | # Select the rectangle select tool | ||

+ | # Draw a rectangle around your scale | ||

+ | # Note the length of the scale in pixels and select '''Analyze > Set Scale...''' | ||

+ | # Give the length, and set the units (mm in this case), and note the dimensions now appear for your image | ||

+ | |||

+ | ===4 Crop the image=== | ||

+ | You need to remove any decoration like scales, annotation, borders, etc. These would show up in the wavenumber domain. Select a good area, then '''Image > Crop''' | ||

+ | |||

+ | [[File:FIJI_crop_image.png]] [[File: FIJI_Crop_menu.png|180px]] | ||

+ | |||

+ | ===5 Run the FFT=== | ||

+ | There are two ways: '''Process > FFT > FFT''' or '''Plugins > Process > Fast FFT'''. The methods are basically the same, just slightly different implementations. | ||

+ | |||

+ | [[file:FIJI_Analyze_FFT.png|300px]] [[File:FIJI_Example_2DFFT.png|300px]] | ||

+ | |||

+ | ===6 Play around with filters=== | ||

+ | To better understand what you're looking at, try running some spatial filters from '''Process > FFT > Bandpass Filter...'''. Be sure to check ''Display filter''. This applies a mask in the 2D wavenumber domain, which results in removing some wavenumbers from the image. For example, to eliminate vertical stripes from an image, we want to eliminate all of the wavenumbers in the ''x'' direction that have a value of 0 in the ''y'' direction — this requires us to cut out a horizontal band in the 2D FFT display: | ||

+ | |||

+ | [[File:FIJI_Bandpass-filter_example.png|300px]] | ||

+ | |||

+ | ==References== | ||

+ | <references /> | ||

+ | |||

+ | ==External links== | ||

+ | * [[:Wikipedia:Wavenumber|Wavenumber]] — Wikipedia entry |

## Latest revision as of 12:53, 2 May 2012

Wavenumber, k, sometimes called the *propagation number*, is a quantification of spatial scale. It can be thought of as a spatial analog to the temporal frequency, and is often called *spatial frequency*. It is often defined as the number of wavelengths per unit distance, or in terms of wavelength, *λ*:

This is analogous to frequency *f*, which is the reciprocal of period *T*; that is, *f* = 1/*T*. In a sense, period can be thought of as a temporal 'wavelength—the length of an oscillation in time.

## Contents

## Angular wavenumber

If you've explored the applications of frequency in geophysics, you'll have noticed that we sometimes don't use ordinary frequency f, in Hertz. Because geophysics deals with oscillating waveforms, ones that vary around a central value (think of a wiggle trace of seismic data), we often use the **circular** or **angular frequency**. In this way, we can also express the close relationship between frequency and phase, which is an angle. So in many geophysical applications, we want the *angular wavenumber*. It is expressed in radians per metre:

The relationship between angular wavenumber and angular frequency is analogous to that between wavelength and ordinary frequency—they are related by the velocity *V*:

It's unfortunate that there are two definitions of *wavenumber*. Some people reserve the term *spatial frequency* for the ordinary wavenumber, or use ν (that's a Greek *nu*, not a vee — another potential source of confusion!), or even *σ* for it. But just as many call it the wavenumber and use *k*, so the only sure way through the jargon is to specify what you mean by the terms you use. As usual!

Here are two spectra: the FFTs of an image of some waves, and a binary image of some particles:

## Computing 2D Fourier transform in FIJI

You can use FIJI to compute the 2D FFT of an image. FIJI is a popular open-source, scriptable scientific image-processing tool. **F**IJI **I**s **J**ust (a more awesome) **I**mageJ.

### 1 Download and install FIJI

### 2 Open your image

This image is a low-res version of a photomicrograph of some shale from Schieber et al (2010)^{[1]}.

**Please note:** This image is not covered by the CC-BY license on the rest of this wiki. It is copyright of the authors. I'm calling the use of a low-res version 'fair use'. I will try to build an example with an open image at some point — feel free to change it!

### 3 Set the scale

No matter what you're doing in FIJI, it's often a good idea to tell FIJI the real-world units of your image, so you don't have to think in pixels. To do this, use the rectangle selection tool to measure a length in pixels (e.g. using a scale in the image), then select **Analyze > Set scale...**

- Select the rectangle select tool
- Draw a rectangle around your scale
- Note the length of the scale in pixels and select
**Analyze > Set Scale...** - Give the length, and set the units (mm in this case), and note the dimensions now appear for your image

### 4 Crop the image

You need to remove any decoration like scales, annotation, borders, etc. These would show up in the wavenumber domain. Select a good area, then **Image > Crop**

### 5 Run the FFT

There are two ways: **Process > FFT > FFT** or **Plugins > Process > Fast FFT**. The methods are basically the same, just slightly different implementations.

### 6 Play around with filters

To better understand what you're looking at, try running some spatial filters from **Process > FFT > Bandpass Filter...**. Be sure to check *Display filter*. This applies a mask in the 2D wavenumber domain, which results in removing some wavenumbers from the image. For example, to eliminate vertical stripes from an image, we want to eliminate all of the wavenumbers in the *x* direction that have a value of 0 in the *y* direction — this requires us to cut out a horizontal band in the 2D FFT display:

## References

- ↑ Juergen Schieber, John B. Southard, and Arndt Schimmelmann (2010). Lenticular shale fabrics resulting from intermittent erosion of water-rich muds—interpreting the rock record in the light of recent flume experiments.
*Journal of Sedimentary Research***80**(1), p 119–128. DOI 10.2110/jsr.2010.005.

## External links

- Wavenumber — Wikipedia entry