Assigment 8 in GEG2210 - Data Collection - Land Surveying, Remote Sensing and Digital Photogrammetry

Image enhancement, filtering and sharpening

By Petter Reinholdtsen and Shanette Dallyn, 2005-05-01.

This exercise was performed by logging into jern.uio.no using ssh and running ERDAS Imagine. Started by using 'imagine' on the command line. The images were loaded from /mn/geofag/gggruppe-data/geomatikk/

We tried to use svalbard/tm87.img, but it only have 5 bands. We decided to switch, and next tried jotunheimen/tm.img, which had 7 bands.

Some notes on the digital images

The pixel values in a given band is only a using a given range of values. This is because sensor data in a single image rarely extend over the entire range of possible values.

The peak values of the histograms represent the the spectral sensitivity values that occure the most often with in the image band being analysed.

Evaluation of the different bands

This image show the "true colour" version, with the blue range assigned to the blue colour, green range to green colour and red range to red colour.

band 1, blue (0.45-0.52 micrometer - um)

Visible light, and will display a broad range of values both over land and water. Reflected from ice, as those are visible white and reflect all visible light waves. Histogram show most values between 30 and 136. Mean values of 66.0668. There are one wide peak with center around 50. There are two peaks at 0 and 255.

band 2, green (0.52-0.60 um)

Visible light, and will display a broad range of values both over land and water. Reflected from ice, as those are visible white and reflect all visible light waves. Histogram show most values from 8 to 120. The mean value is 30.9774. There are two main peaks at 20 and 27. There is also a pie at 0.

band 3, red (0.60-0.69 um)

Visible light, and will display a broad range of values both over land and water. Reflected from ice, as those are visible white and reflect all visible light waves. Histogram show most values from 33 t 135, with one wide peak around 52. There are also seem to be two peaks at 0 and 255. The mean value is 34.3403.

band 4, near-infraread (0.76-0.90 um)

Water acts as an absorbing body so in the near infrared spectrum, water features will appear dark or black meaning that all near infrared bands are absorbed. On the other hand, land features including ice, act as reflector bodies in this band. The histogram show most values between 7 and 110. The mean is 40.1144. There are two peaks at 7 and 40.

band 5, mid-infrared (1.55-1.75 um)

The ice, glaciers and water do not reflect any mid-infrared light. The histogram show most values between 1 and 178. The mean is 49.8098 and there are two peaks at 6 and 78, in addition to two peaks at 0 and 255.

band 6, thermal infrared (10.4-12.5 um)

Display the temperature on earth. We can for example see that the ice is colder than the surrounding areas. The histogram show most values between 36 to 122. The mean is 102.734. There are one wide peak around 53, in addition to two peaks at 0 and 255.

band 7, mid-infrared (2.08-2.35 um)

The ice, glaciers and water do not reflect any mid-infrared frequencies. The histogram show most values between 77 and 150. The mean is 24.04, and there are one wide peak at 130 and a smaller peak at 83, in addition to one peak at 0.

Image enhancement

When we look at the linear contrast functions, we can move the slope and shift values increasing or decreasing the contrast of the image. For example, in the linear contrasting we moved the slope value from 1.00 to 3.00 to obtain a brighter appearing image, and then we moved the shift from 0 to 10 to recieve a sharper image.


Next we tried the piecewise linear stretching for contrast. In this image we tried to make all of the histograms in the red, blue and green spectrum as similar as possible so we could detect a change in the image.(insert histogram change)
We tried to break the slope and move the break point to slightly after each histogram peak. This resulted in the image obtaining a slightly blue tint and dullness. (put ugly blueish picture here) As this result was not really increasing the contrast, we tried another variation to try to spread out the histogram peak to use a wider range. This setting gave an improved image, were it is easier to see the red vegetation and the white ice.


We also tried to do histogram equilization on the standard infrared composition. This changed the colours in the image, making the previously green areas red, and the brown areas more light blue. In this new image, we can clearly see the difference between two kind of water, one black and one green. We suspect the green water might be deeper, but do not know for sure.

We can get best contrast stretch by using the histogram equalisation. This gave us the widest range of visible separation between features.

Displaying colour images

Comparing a map we found on the web, and the standard infrared image composition, we can identify some features from the colors used:

Next, we tried to shift the frequencies displayed to use blue for the red band, green for the near ir band and red for the mid ir (1.55-1.75 um). With this composition, we get some changes in the colours of different features:

Filtering and image sharpening

'

We decided to work on the grey scale version of the near infrared (band4). We changed the colour assignment to use this band for all three colours, giving us a gray scale image.

'

We applied the 3x3 low pass filter on this image, and this gave us almost the same image as the original. If you look closely you can see that some white dots in the original disapper, and some of the water edges seem to blur very slightly.

'

We also tried the 3x3 high pass filter on the band4 grey scale image. This gave a very noisy image. Edges of vallies and ice are not well defined. The black waters are still obvious.

'

We also tried the 3x3 edge detection, and this gave us an image that makes it difficult to distinguish elevation features such as the valleys. Rather, edge detection allows us to study main features in an area like the lakes. (insert band4 edge 3 image) '

We tried a gradient filter using this 3x3 matrix. The matrix was chosen to make sure the sum of all the weights were zero, and to make sure the sum of horizontal, vertical and diagonal numbers were zero too.

12-1
20-2
1-2-1

The gradient filter used gave us enhancement on lines in the vertical, horizontal and diagonal directions. This is seen by the white lines that outline certain areas of main features like the rivers within the vallies and some of the lakes.

'

When we rework the matrix to equal negative one, we end up with a lot of noise in the image that also seems to blurr the image. Using a negative one matrix is not optimal if you are trying to obtain sharpness.

-1-1-1
-17-1
-1-1-1

'

We then tried with a 3x3 matrix were the sum of all values equals 1, to enhance the high frequency parts of the image.

-1-1-1
-19-1
-1-1-1

This gave us a sharper looking image compared to the result of the negative 1 filter. This is not really obvious unless one is comparing the two images carefully. In order to see more differences the matrix sums would have to be more then plus/minus one.

References


Petter Reinholdtsen
Last modified: Sun May 1 14:28:48 CEST 2005