Mater`s Thesis

SLOVAK UNIVERSITY OF TECHNOLOGY IN BRATISLAVA
FACULTY OF CIVIL ENGINEERING
Reg. No.: SvF-5330- 67547
MONITORING LANDSCAPE CHANGES
USING SATELLITE RADAR IMAGERY
MASTER'S THESIS
2015
Bc. Jakub Vanko
SLOVAK UNIVERSITY OF TECHNOLOGY IN BRATISLAVA
FACULTY OF CIVIL ENGINEERING
Reg. No.: SvF-5330- 67547
MONITORING LANDSCAPE CHANGES
USING SATELLITE RADAR IMAGERY
MASTER'S THESIS
Study programme:
Study field number:
Study field:
Training workplace:
Thesis supervisor:
Consultant :
Bratislava 2015
Geodesy and Cartography
3636
5.1.3. Geodesy and Cartography
Department of Theoretical Geodesy
Ing. Juraj Papčo, PhD.
Doc. RNDr. Zuzana Krivá, PhD.
Bc. Jakub Vanko
Affidavit
I hereby declare that this master's thesis was written by myself on the basis of
theoretical and practical knowledge obtained during my studies. No further sources were
used than those explicitly stated. All citations quoted from publications in this thesis are
labelled as such.
Bratislava 21. May 2015
.....................
Bc. Jakub Vanko
Acknowledgements
TerraSAR-X radar imagery for this master's thesis was kindly provided by the
German Aerospace Centre (DLR) within the project ID LAN1583: Object Recognition
Based on High-Resolution Radar Imagery and the project LAN2833: Utilization of the
High Spatio-temporal Resolution of TerraSAR-X Observations for Recent Ground
Deformation Monitoring and SAR Image Processing. Sentinel-1 dataset was provided by
the Sentinel-1 Scientific Data Hub of the European Space Agency.
I would like to thank my supervisor and now a fellow colleague from the insar.sk
team Ing. Juraj Papčo, PhD. for introducing me to the beautiful science of SAR, for
providing me hints and new insights and steering me on the right track to completing this
master's thesis.
I also appreciate the help of Doc. RNDr. Zuzana Krivá, PhD. from the Department
of Mathematics and Descriptive Geometry SvF STU for answering my beginner questions
about image processing, hours of consultations and also for providing the Perona-Malik
equation based filter for denoising SAR images which were used in my thesis.
I am also very thankful to Ing. Matúš Bakoň for accepting me to the insar.sk team,
for helping me with the final stages of this thesis and mostly for inspiring me to work
harder and to seek more and more knowledge.
Last but not least I would like to thank my family for supporting me thru all the
hard times of my studies. I am forever in your debt and I would like to dedicate this thesis
to you.
Abstract
The general aim of this master's thesis is to test the application of satellite radar
imagery on detection of landscape elements and their changes in time. The experiment was
conducted in the vicinity of Bratislava, the capital city of Slovak Republic. Radar imagery
was provided by the TerraSAR-X mission and by the new European Space Agency
initiative, the Sentinel-1 mission. All radar images are inherently degraded by granular
noise known as speckle and before application the very first step is to reduce the amount of
this effect. For this purpose several speckle reduction filters were quantitatively compared
and evaluated on real radar imageries. After speckle reduction, it was possible to perform
object detection. The obtained results were compared with in situ measurements using
global navigation satellite systems. The results have been presented by filtered images,
statistical tables and diagrams.
Key words: Synthetic Aperture Radar, speckle noise reduction, spatial filtering, edge
detection, TerraSAR-X, Sentinel-1
Abstrakt
Hlavným cieľom tejto diplomovej práce je testovanie využitia družicových
radarových snímok na detekciu prvkov krajiny a ich zmien v čase. Experiment bol
zrealizovaný v okolí Bratislavy, hlavného mesta Slovenskej republiky. Radarové dáta boli
poskytnuté z misie TerraSAR-X a novej misie Európskej vesmírnej agentúry, Sentinel-1.
Každá radarová snímka je do určitej miery znehodnotená speckle šumom a pred jej
použitím musí byť tento efekt zredukovaný. K tomuto účelu práca obsahuje kvantitatívne
porovnanie a zhodnotenie viacero speckle filtrov otestovaných na reálnych radarových
snímkach. Po redukcii šumu bolo možné vykonať detekciu prvkov krajiny. Pre
zhodnotenie presnosti detegovaných prvkov boli výsledky porovnané s in situ meraním
pomocou globálnych navigačných satelitných systémov. Výsledky sú prezentované
odšumenými snímkami, štatistickými tabuľkami a grafmi.
Kľúčové slová: Radar so syntetickou apertúrou, redukcia speckle šumu, filtrácia
digitálnych obrazov, detekcia hrán, TerraSAR-X, Sentinel-1
Contents
List of Appendices ............................................................................................................... 11
List of Acronyms and Abbreviations ................................................................................... 12
List of Symbols ................................................................................................................... 14
Introduction .......................................................................................................................... 15
1.
2.
3.
Theoretical Aspects of SAR Technology .................................................................... 16
1.1.
SAR Introduction .................................................................................................. 16
1.2.
Complex SAR Image ............................................................................................ 17
1.3.
Detected SAR Image ............................................................................................. 17
1.4.
Commonly Used Radar Frequencies in RS ........................................................... 18
1.5.
SAR Polarization................................................................................................... 19
1.6.
SAR Viewing Geometry ....................................................................................... 20
Speckle Noise Reduction ............................................................................................. 21
2.1.
Speckle Noise Definition ...................................................................................... 21
2.2.
Speckle Suppression Techniques .......................................................................... 22
2.3.
Tested Filters ......................................................................................................... 22
2.4.
Speckle Filters Evaluation..................................................................................... 23
2.4.1.
Equivalent Number of Looks ......................................................................... 24
2.4.2.
Speckle Suppressions Index........................................................................... 24
2.4.3.
Speckle Suppression and Mean Preservation Index ...................................... 24
2.4.4.
Edge-Enhancing Index ................................................................................... 25
2.4.5.
Feature-Preserving Index ............................................................................... 25
2.5.
Testing SAR Dataset ............................................................................................. 26
2.6.
Speckle Filter Evaluation Results ......................................................................... 27
Object Detection on SAR Images ................................................................................ 31
3.1.
Area of Interest ...................................................................................................... 31
3.2.
TSX Imagery ......................................................................................................... 32
3.2.1.
TSX Mission Description .............................................................................. 32
9
3.2.2.
LAN 1583 Dataset ......................................................................................... 33
3.2.3.
LAN 2833 Dataset ......................................................................................... 34
3.3.
Sentinel-1 Imagery ................................................................................................ 35
3.3.1.
Sentinel-1 Mission Description ..................................................................... 36
3.3.2.
Sentinel 1 Dataset .......................................................................................... 37
3.4.
Edge Detection ...................................................................................................... 38
3.5.
In Situ Measurements ............................................................................................ 43
3.6.
Comparison with In Situ Measurements ............................................................... 44
Conclusion ........................................................................................................................... 52
Resume................................................................................................................................. 53
Bibliography ........................................................................................................................ 56
Appendices........................................................................................................................... 60
10
List of Appendices
Appendix A
Speckle Evaluation Results Part 1
Appendix B
Speckle Evaluation Results Part 2
Appendix C
Speckle Evaluation Results Sorted by ENL
Appendix D
Speckle Evaluation Results Sorted by SSI
Appendix E
Speckle Evaluation Results Sorted by SMPI
Appendix F
Speckle Evaluation Results Sorted by EEI
Appendix G
Speckle Evaluation Results Sorted by FPI
11
List of Acronyms and Abbreviations
AOI
Area of Interest
ASCII
American Standard Code for Information Interchange
C-band
Microwave band, frequency 4-8 GHz, wavelength 3.75-7.5 cm
CEOS
The Committee on Earth Observation Satellites
DLR
German Aerospace Centre
EEI
Edge-Enhancing Index
EM
Electromagnetic
ENL
Equivalent Number of Looks
EPSG
European Petroleum Survey Group
ESA
European Space Agency
EU
European Union
FPI
Feature-Preserving Index
GMES
Global Monitoring for Environment and Security (Copernicus)
GMT
Greenwich Mean Time
GNSS
Global Navigation Satellite Systems
GNU
GNU's not Unix
GPL
General Public License
GUI
Graphic User Interface
HH
SAR polarization, horizontal transmission, horizontal reception
HV
SAR polarization, horizontal transmission, vertical reception
IW
Interferometric Wide Swath Mode
K-band
Microwave band, frequency 18-26.5 GHz, wavelength 11.3-16.7 mm
Ka -band
Microwave band, frequency 26.5-40 GHz, wavelength 5-11.3 mm
Ku-band
Microwave band, frequency 12-18 GHz, wavelength 16.7-25 mm
L-band
Microwave band, frequency 1-2 GHz, wavelength 15-30 cm
MAP
Maximum a Posteriori
NE
Northeast
NMS
Non-maximum Suppression
NW
Northwest
P-band
Microwave band, frequency 0.3-1 GHz, wavelength 30-100 cm
PGM
Portable Graymap Format
PGS
Payload Ground Segment
12
RADAR
Radio Detection and Ranging
RAR
Real Aperture Radar
RS
Remote Sensing
SAR
Synthetic Aperture Radar
S-band
Microwave band, frequency 2-4 GHz, wavelength 7.5-15 cm
SE
Southeast
STBX
Sentinel-1 Toolbox
SM
StripMap imaging mode of the TerraSAR-X satellite
SMPI
Speckle Suppression and Mean Preservation Index
SSI
Speckle Suppression Index
ST
SpotLight imaging mode of the TerraSAR-X satellite
STD
Standard Deviation
STU
Slovak University of Technology in Bratislava
SVF
Faculty of Civil Engineering
SW
Southwest
TIFF
Tagged Image File Format
TSX
TerraSAR-X
VAR
Variance
VH
SAR polarization, vertical transmission, horizontal reception
VV
SAR polarization, vertical transmission, vertical reception
WGS 84
World Geodetic System 1984
X-band
Microwave band, frequency 8-12 GHz, wavelength 25-37.5 mm
XLS
Excel Binary File Format
13
List of Symbols
C
Number of columns in SAR image
f
SAR image of true radiometric values
g
Original image value
g1
Original image pixel values on one side of the edge
g2
Original image pixel values on the other side of the edge
gf1
Filtered image pixel values on one side of the edge
gf2
Filtered image pixel values on the other side of the edge
gf
Filtered amplitude image
n
Kernel dimension
Q
SMPI coefficient
R
Number of rows in SAR image
u
Speckle noise
y
Image column coordinate
x
Image line coordinate
β
Ration of the standard deviation to the mean of a filtered image
σu
Speckle noise standard deviation
14
Introduction
Currently our planet faces an increased threat of natural disasters due to climate
changes caused by either natural processes or human activities. This is the reason why
many countries and intergovernmental organizations such as the European Space Agency
(ESA) are interested in monitoring the state of the environment, oceans and atmosphere.
The goal is to gather accurate and timely information to improve the management of the
environment, understand and mitigate the effects of climate changes and ensure civil
security. The single location where we can learn the most about our planet is found
nowhere on Earth but high up above it (unknown ESA's spokesperson, 2013). A powerful
tool for Earth observation is remote sensing (RS). RS is the science and art of obtaining
information about an object, area, or phenomenon through the analysis of data acquired by
a device that is not in contact with the object, area, or phenomenon under investigation
(Lillesand, Kiefer, Chipman, 2004). Based on the wavelength in which the system
operates, RS is divided into optical and microwave RS.
Microwave RS using Synthetic Aperture Radar (SAR) offers many advantages and
new applications over optical RS but it also has its weaknesses. Microwave radiation can
penetrate through cloud cover and most weather conditions. Because radar is an active
sensor, it can also be used to image the surface at any time, day or night.
All SAR images are inherently degraded by granular noise known as radar speckle.
Speckle is caused by random constructive and destructive interference from multiple
backscatters and causes difficulties for analysis, classification and image interpretation for
various RS applications, including object detection. Before using SAR images, the very
first step is to reduce the effect of speckle. One of the common methods which are used to
reduce speckle noise is spatial filtering. For the purpose of this research, five speckle filters
were quantitatively evaluated. The filter showing the best results was latter used to reduce
the speckle effect on SAR images used for detection of landscape elements.
This thesis is divided into three main chapters. The first chapter contains theoretical
aspects of SAR. The second chapter deals with speckle noise and its reduction techniques.
The goal of this chapter is to choose the most suitable speckle filter, denoise the SAR
images and prepare them for edge detection. Finally, chapter three describes edge detection
and edge extraction from despeckled SAR images. The detected edges were then compared
with more conventional data collection method, the global navigation satellite systems
(GNSS).
15
1. Theoretical Aspects of SAR Technology
SAR is an active microwave imaging system. The microwave portion of the
electromagnetic spectrum covers the range from approximately 1 cm to 1m in wavelength.
These relatively long wavelengths are independent of atmospheric scattering and can
penetrate through clouds, dust, fog or heavy rainfall. Active sensor supplies its own source
of energy to illuminate the target. The advantage gained from the use of active sensors, is
the possibility to acquire data also during the night.
1.1. SAR Introduction
The word radar is an acronym for RAdio Detection And Ranging. As its name
implies, radar was developed as a means of using radio waves to detect the presence of
objects and to determine their distance and sometimes their angular position. The process
entails transmitting short bursts, or pulses, of microwave energy in the direction of interest
and recording the strength and origin of "echoes" or ''reflections'' received from objects
within the system's field of view. Radar systems may or may not produce images, and they
may be ground based, or mounted in aircraft or spacecraft (Lillesand, Kiefer, Chipman,
2004).
SAR transmits pulses of microwave energy towards the Earth's surface and detects
the reflected backscatter from each target. The received backscattered signals are separated
into two components which are carrying information about the amplitude and the phase of
the returning signal. Amplitude depends on target properties (structure and dielectric
properties). Phase is a function of the distance between the sensor and the target as well as
target properties. Together they form a complex number. To create an image, the returning
signal of a single pulse is sampled and these samples are stored in an image line. As the
radar platform moves forward, recording and processing of the backscattered signals build
up a two-dimensional image of the surface.
SAR systems were developed to overcome the limitations of Real Aperture Radar
(RAR) systems. Compared to RAR, SAR synthetically increases the antenna's size to
achieve finer resolution in the platform motion direction. RAR systems have their physical
limit to the length of the antenna that can be carried on an aircraft or satellite. SAR
increases the size of the antenna or aperture electronically.
16
1.2. Complex SAR Image
A SAR image is a two-dimensional array of pixels formed by columns and rows
where a pixel is associated with a small area of Earth's surface whose size depends only on
the SAR system characteristics. Each pixel provides a complex number (amplitude and
phase information) associated to the reflectivity of all the scatters contained in SAR
resolution cell (Lee, Pottier, 2009).
1.3. Detected SAR Image
The detected SAR Image contains a measurement of the amplitude of the radiation
backscattered toward the radar by the objects (scatters) contained in each SAR resolution
cell (Ferreti et al., 2007).
The detected amplitude depends on how the radar energy interacts with the surface,
which is a function of several variables or parameters. These parameters include the
particular characteristics of the radar system (viewing geometry, frequency, polarization,
etc.) as well as characteristics of the surface (land cover type, topography, relief, etc.).
Generally the dominant factor in determining the tones in radar images is surface
roughness. Exposed rocks and urban areas show strong amplitudes, whereas smooth
surfaces like water show low amplitudes, since the radiation is reflected away from the
radar. Amplitude images are generally visualised by means of grey scale levels.
Assume that a detected SAR image g(x,y) is sampled so that the resulting digital
image has R rows and C columns. Figure 1.1 shows the coordinate convention used
throughout this thesis. The values of the coordinates (x,y) now become discrete numbers
quantities. The origin coordinates are (x,y) = (0,0). The next coordinate value along the
first row of the image are represented as (x,y) = (1,0). This knowledge is essential for
proper transformation of the detected landscape elements from image coordinates to
projected map coordinates.
Figure 1.1: Coordinate convention used to represent digital images in this thesis.
17
1.4. Commonly Used Radar Frequencies in RS
Similarly to optical RS, even SAR sensors operate with one or more different
bands. Radar waves fall under the microwave region of the electromagnetic (EM)
spectrum. The microwave region of the EM spectrum is quite large, relatively to the visible
and infrared. There are several wavelength ranges or bands commonly used, which were
given code letters during World War II and remained to this day (CCRS, 2014):
•
Ka, K and Ku bands: very short wavelengths used in early airborne radar systems
but uncommon today,
•
X-band: used extensively on airborne systems for military reconnaissance and
terrain mapping, X-band is also used by the TerraSAR-X satellite, whose data is
utilized in this thesis,
•
C-band: common on many airborne (CCRS Convair-580 and NASA AirSAR) and
spaceborne systems including ERS-1 and 2, RADARSAT , Sentinel-1,
•
S-band: used on board of the Russian ALMAZ satellite,
•
L-band: used on board of the American SEASAT, Japanese JERS-1 satellites,
NASA airborne system, ALOS, ALOS-2, L-band will also be used in the planned
Tandem-L mission,
•
P-band: longest radar wavelengths, used on NASA experimental airborne research
system and also planned for the ESA Biomass mission.
Figure 1.2: Microwave region of the EM spectrum.
(Adapted from (CCRS, 2014))
18
1.5. SAR Polarization
The microwave polarization refers to the orientation of the electric field vector of
the transmitted beam with respect to the horizontal direction. If the electric field vector
oscillates along a direction parallel to the horizontal direction, the beam is said to be
horizontally polarized. On the other hand, if the electric field vector oscillates along a
direction perpendicular to the horizontal direction, the beam is vertically polarized
(CRISP, 2001).
Using different polarizations and wavelengths, it is possible to collect information
that is useful for particular applications, e.g. to classify agricultural fields. In radar system
descriptions the reader can come across the following abbreviations (Tempfli et al., 2009):
•
HH: horizontal transmission and horizontal reception,
•
VV: vertical transmission and vertical reception,
•
HV: horizontal transmission and vertical reception,
•
VH: vertical transmission and horizontal reception.
Older SAR sensors were only capable to provide single polarization channel
images. Modern SAR sensors are often dual-polarized or even capable of providing full
polarization datasets (HH, HV, VV, VH).
Figure 1.3: Example of few SAR satellites and their polarization capabilities.
(Adapted from (MDA, 2014))
19
1.6. SAR Viewing Geometry
SAR imaging geometry is different from systems employed for optical RS. SAR
sensors are side-looking instruments while most optical instruments are nadir-looking.
This difference exists because optical sensors are able to distinguish among targets using
angular distance from the sensor. SAR can only distinguish the returns from various targets
based upon the arrival time of the received signals. A nadir-looking radar could not
distinguish between two scatters with equal distances from the sensor because a single
incident wave front illuminates both points at the same instant, so the scatter returns from
both points arrive at the same receiving antenna simultaneously (Austin Center for Space
Research, 2015).
SAR viewing geometry is explained in Figure 1.4. The SAR platform (A) travels
forward along the flight path or orbit (B) with the nadir directly beneath the platform (C).
The projection of the flight path on the surface is called the ground track (D). The
microwave beam illuminates an area or swath (E). The along-track dimension parallel to
the flight direction refers to azimuth (F), while the cross-track dimension refers to range
(G). The portion of the image swath closest to the ground track of the radar platform is
called the near range while the portion of the swath farthest from the ground track is called
far range. The look angle (H) is the angle at which the radar looks at the surface. The
incidence angle (I) of the system is defined as the angle between the radar beam and the
local surface vertical. Moving from near range to far range, the incidence angle increases.
At all ranges the radar antenna measures the radial line of sight distance between the radar
and each target on the surface. This line is called slant range (J). The true horizontal
distance along the ground corresponding to each measured point in slant range is called
ground range (K).
Figure 1.4: SAR viewing geometry.
20
2. Speckle Noise Reduction
Unlike optical RS images, characterized by very neat and uniform features, SAR
images are affected by speckle noise. Speckle forms a main obstacle to analyse, interpret
and classify SAR images for various RS applications. Figure 2.1 displays a subset of a
detected SAR image degraded by speckle noise.
Figure 2.1: Detected SAR image: the speckle effect on the homogenous
fields is clearly visible.
2.1. Speckle Noise Definition
Speckle is a direct result of the fact that the incident energy is coherent - that is, it
can be assumed to have a single frequency and the wavefront arrives at a pixel with a
single phase. If there were a single large dominant scatter in the pixel, such as a corner
reflector or a building, then the returned signal would be largely determined by the
response of that dominant element, and any scattering from the background would be
negligible. More often, though the pixel will be a sample of very large number of
incremental scatters; their returns combine to give the resultant received signal for that
pixel (Richards, 2009). Such situation is illustrated in Figure 2.2.
Figure 2.2: Generation of speckle through the interference of a very large number of rays
scattered from within a pixel. (Adapted from (Richards,2009))
21
Speckle can be statistically characterized by multiplicative noise model. Degraded
image with speckle noise is given by the equation (Sarode, Deshmukh, 2011):
g ( x , y ) = f ( x, y ) ⋅ u ( x, y )
(2.1)
Where, g (x,y) is the observed amplitude image with speckle, f (x,y) is the image of
the true radiometric values and u (x,y) is the speckle noise. (x,y) denotes the pixel location.
The noise is characterized by a distribution with a unit mean (E[u]=1) and a standard
deviation σu. The multiplicative nature of speckle complicates the noise reduction process.
2.2. Speckle Suppression Techniques
Generally speaking, speckle noise can be reduced by multi-look processing or
spatial filtering (Raney, 1998). Multi-looking or multiple-look processing is usually done
during data acquisition and speckle reduction by spatial filtering is performed on the output
image in a digital image analysis environment. Spatial filtering is a neighbourhood
operation which works with the values of the image pixels in the neighbourhood and the
corresponding values of a subimage that has the same dimensions as the neighbourhood.
The subimage is called a filter, mask, kernel, template, or window, with the first three
terms being the most prevalent terminology. The values in a filter subimage are referred to
as coefficients, rather than pixels (Gonzales, 2008). The mask of a few pixels in dimension
(e.g. 3x3, 5x5, etc.) is moving over each pixel in the image, applying a mathematical
calculation using the pixel values under that mask, and replacing the central pixel with the
new value. The kernel is moved along the image one pixel at a time until the entire image
has been covered. By applying the filter a smoothing effect can be achieved and the visual
appearance of the speckle is reduced. The ideal filter should preserve object edges, textural
and radiometric information. This thesis evaluates the effect of Frost, Lee, Refined Lee,
Gamma MAP and Perona-Malik equation based filter.
2.3. Tested Filters
This subchapter contains basic information about the evaluated speckle filters. For
more detailed description and mathematical principles, readers are referred to relevant
publications quoted in the text.
The Frost filter replaces the pixel of interest with a weighted sum of the values
within the n x n moving kernel. It is based on the local statistics and the multiplicative
noise model. The weighting factors decrease with distance from the pixel of interest. They
22
increase for the central pixels as variance within the kernel increases with the edge areas to
preserve edge structure (Mansourpour, Rajabi, Blais, 2006).
Gamma or maximum a posteriori (MAP) filter minimizes the loss of texture
information in gamma-distributed scenes which are the images of forested areas,
agricultural land and oceans. In this filtering technique, the smoothing process is
determined by the coefficient of variation and contrast ratios of probability density
functions (Shanthi, Valarmathi, 2013).
The Lee filter utilizes the statistical distribution of digital number within the
moving kernel to estimate the value of pixel of interest. This filter assumes Gaussian
distribution for the noise in image data (Lee, 1981). The Lee filter removes the noise by
minimizing either the mean square error or the weighted least square estimation.
Refined Lee Filter is based on the Lee filter. It considers a Local Linear Minimum
Mean Square Error estimation within edge aligned windows, achieving a better
preservation of the image details.
Perona-Malik filter was tested with the cooperation of Department of Mathematics,
Faculty of Civil Engineering (SVF), Slovak University of Technology in Bratislava (STU).
Perona-Malik nonlinear image selective smoothing equation (called anisotropic diffusion
in image processing) was proposed in 1990 by Perona and Malik. It usually works well for
images disturbed by additive noise, provided that the gradients of noise are smaller than
the edge gradients. The multiplicative character of the speckle is changed into additive one
by logarithming. The explored algorithm works on an adaptive grid, i.e., for homogenous
areas it uses larger elements of the computational grid. This filter does not preserve
radiometric properties, its aim is to support edge detection.
2.4. Speckle Filters Evaluation
To assess the capability of the tested filters to remove speckle noise, and their
effectiveness in successfully preserving the real structure of the scene backscatter, several
performance measures were used. The following measures were also used to assess the
performance of several iterations of filtering and also different kernel sizes 3x3, 5x5, and
7x7. This doesn't apply to Perona-Malik filter, which works on an adaptive grid. Also
Refined Lee filter was unable to set a different kernel size than 7x7 pixels. In addition to
these measures visual inspection also provides a good assessment of the speckle filters
23
performance. It was an efficient and easy way to assess the speckle reduction capability
and also effectiveness in edge preserving.
2.4.1. Equivalent Number of Looks
Equivalent Number of Looks (ENL) for amplitude image is calculated using the
following equation (Lee et al. 1994), (Lee, Pottier, 2009):
ENL = ( 0.5527 / β )
2
(2.2)
Where 0.5527 is the value of the σu of a 1-look amplitude SAR image. This definition is
consistent with the multi-look processing using amplitude averaging of individual looks.
According to (Lee et al., 1994) and (Lee, Pottier, 2009) β is the ration of the standard
deviation to the mean of the image and it is defined as the coefficient of variance (speckle
index):
var( g f )
β=
E  g f 
(2.3)
Where gf is the pixel value of the filtered SAR image.
2.4.2. Speckle Suppressions Index
The Speckle Suppression Index (SSI) is the coefficient of variance of the filtered
image normalized by that of the original image, which is defined as (Sheng, Xia, 1996):
SSI =
var( g f )
E  g f 
⋅
E[g]
var( g )
(2.4)
Where gf is the filtered value, g is the original value. For most cases SSI < 1.0, which
means speckle is suppressed. The lower the SSI is, the stronger suppression ability the
filter has.
2.4.3. Speckle Suppression and Mean Preservation Index
According to (Shamsoddini, Trinder, 2010) ENL and SSI are not reliable when the
filter overestimates the mean value. They have developed an index called Speckle
Suppression and Mean Preservation Index (SMPI). The equation of this index is as follow:
SMPI = Q ⋅
24
var( g f )
var( g )
(2.5)
And Q is calculated as:
Q = 1 + E [ g ] − E [ gˆ 0 ]
(2.6)
The lower values indicate better performance of the filter in terms of mean preservation
and noise reduction.
2.4.4. Edge-Enhancing Index
The Edge-Enhancing Index (EEI) is an important parameter in assessing the
enhancement of step functions in an image. It is often used to evaluate a filters's ability to
preserve edges, such as boundaries between water bodies and land. The EEI is defined as
(Qiu et al., 2004):
EEI =
∑g
∑g
f1
1
− gf2
− g2
(2.7)
Where g1 and g2 are the original values of the pixels on either side of the edge, whereas gf1
and gf2 are the corresponding filtered values. The numerator is the absolute difference in
intensity of the pixels on the two sides of the edge in the filtered image and the
denominator is that difference in the original image. Therefore, the EEI is usually smaller
than 1.0 and higher EEI values correspond to a better edge preserving capability.
2.4.5. Feature-Preserving Index
Feature-Preserving Index is a measure for assessing a filter's ability to preserve
linear features and subtle structures. For a one-pixel wide linear feature of n-pixel length,
the FPI is given by (Qiu et al., 2004):
∑ (2⋅ g
n
FPI =
1
f
− gf1 − gf 2 )
n
∑ (2 ⋅ g − g
1
− g2 )
(2.8)
1
Where g is the original value of a pixel on the linear feature. g1 and g2 are the
original values of the neighbouring pixels on both sides of the feature, and gf , gf1 and gf2
are the filtered values of corresponding pixels. In most cases, the values of FPI measure are
also lower than 1. The higher the FPI measure, the better the linear feature.
25
2.5. Testing SAR Dataset
The test image is a 4000 x 4000 pixel subset from a detected SAR image. It was
selected from the dataset intended for project ID LAN1583: Object Recognition
Recogni
Based on
High-Resolution
Resolution Radar Imagery. The dataset was kindly provided by the German
Aerospace Center (DLR). The image was acquired by the TerraSAR-X
TerraSAR (TSX) mission
(X-band,
band, 31 mm wavelength) on December 3, 2008 with Spotlight acquisition mode, VV
polarization
ation and sample spacing 0.75 m. In order to test the particular filters with ENL, SSI
SS
and SMPI, a homogenous segment was selected from the SAR image. The area is
represented with the purple rectangle in Figure 2.3.
2.3 For edge preservation testing with EEI
and FPI indices, a cross section of the agricultural field's edge was selected. It is
represented in Figure 2.3 with the small yellow line.
Figure 2.3: Subset of the TerraSAR-X
TerraSAR image from December 3, 2008.
2008 The purple rectangle
shows the selected homogeneous area, the yellow line represents the edge cross section.
26
2.6. Speckle Filter Evaluation Results
Initial processing and speckle reduction was performed in Sentinel-1 Toolbox
(S1TBX). S1TBX was developed under contract to ESA by Array Systems Computing Inc.
and contributors. The toolbox and full source code is distributed freely under the GNU
GPL license. The S1TBX is consisting of a collection of processing tools, data product
readers, writers and a display and analysis applications to support data from ESA SAR
missions including Sentinel-1 and third party SAR data from TerraSAR-X (S1TBX Help,
2015). S1TBX is programmed in Java, it provides user friendly graphic user interface
(GUI) and supports Windows, Linux and MAC OS X operating systems. S1TBX was used
mostly for reading Sentinel-1 and TSX data, creating subsets, speckle filtering, map
projection and GeoTiff export.
The Perona-Malik filter was a standalone program written in C programming
language. The program enables to read ASCII images in Portable Graymap Format (PGM)
and set different parameters of the smoothing effect.
Manually calculating all parameters for each smoothed image would be time
consuming with a greater possibility of calculation errors. A simple MATLAB function
was created to overcome this problem. The function is called Speckleeval. This function
quickly and easily automates the calculation. Speckleeval reads both the original image and
the smoothed image, calculates all evaluation parameters and writes the results into Excel
Binary File Format (XLS).
Mean values, variations (VAR), standard deviations (STD) and all performance
measures mentioned above were computed for 55 smoothed images. Generally increasing
the filters kernel size and the number of iterations improved speckle reduction but
decreased edge preservation. Numerical results for each filter and variable kernel sizes are
presented in Appendix A and Appendix B. For example label 3Frost1 means that the SAR
image was smoothed by one iteration of Frost filter using a 3x3 pixel size kernel. Same
labelling was used in the following Appendices.
After speckle reduction all smoothed SAR images were submitted to edge detection
with the Canny Edge Detection Algorithm. The most suitable results from edge detection
were achieved on images smoothed by applying a 7x7 pixel size kernel, so only these
results were compared including Perona-Malik filter with adaptive grid. All results sorted
by each measure are listed in Appendices C to G. Best speckle reduction effect was
achieved by Gamma MAP filter according to ENL, SSI and SMPI indices. Lee and Frost
27
filter's ability to reduce speckle after 5 iterations was nearly identical to Gamma MAP
filter. Poorest results were achieved by Refined Lee filter and Perona-Malik filter after 5
iterations. Overall, all filters except Perona-Malik are capable of preserving the mean of a
homogenous segment well. Emphasis of this thesis is on edge detection, therefore EEI and
FPI indices have greater weight. Appendix F shows numerical values of the mean, VAR,
STD and all performance measures of all 5 adaptive speckle reduction filters sorted by EEI
in descending order. Appendix G shows results sorted by FPI also in descending order. The
higher EEI and FPI indices, the better edge preservation. From Appendix F and Appendix
G it is clear that best results in preserving edges were obtained by Refined Lee and PeronaMalik filter. Although according to EEI the Perona-Malik filter after 5 iterations performed
worse than all other filters, on average it still performs exceptionally well.
To demonstrate the smoothing effect of the speckle filters, close-up of the original
image and the same image smoothed by 5 iterations of Lee, Frost, Gamma MAP, Refined
Lee (with kernel size 7x7 pixels) and by 40 iterations of Perona-Malik filter are displayed
in figure 2.4. First of all, it can be clearly seen that speckle noise is significantly reduced in
the filtered images. Secondly, although the filtered images appear to be blurry within a
land cover segment, differences between neighbouring land covers are visibly enhanced.
28
(b)
(a)
(c)
(d)
(e)
(f)
Figure 2.4: Original image (a) and smoothed images by Lee (b), Gamma MAP (c),
Refined Lee (d), Perona-Malik after 40 iterations (e) and Frost filter (f).
This thesis focuses on landscape objects detection, so the main criterium is to
choose the filter with the best edge preserving capability. It would seem that there are two
obvious candidates Refined Lee and Perona-Malik filter. However applying the edge
detector and comparing the detected edges revealed some interesting new facts. The
situation is shown in Figure 2.5. For example, the detected edge on Gamma MAP filtered
image was much smoother than the edge detected on the Perona-Malik smoothed image.
29
(b)
(a)
(c)
(d)
Figure 2.5: Smoothed images by Gamma MAP filter (a), Perona-Malik filter (b), and
detected edges on Gamma MAP smoothed image (c) and Perona-Malik smoothed image
(d).
Despite the numerical evaluation of all filters, visual comparison proofed to be the
most reliable measure. For this reason, the final decision was made in favour of the
Gamma MAP filter. Except the tests presented in this thesis, numerous other test were
conducted on different SAR images bringing different results. The conclusion is that there
is no single best filter for every possible scenario. The most suitable filter depends on
many parameters and variables e. g. SAR viewing geometry, polarization, land cover,
prerequisites for certain applications etc.. Also the most reliable and quickest measure for
selecting the most suitable filter for images being processed in practical use is visual
assessment.
30
3. Object Detection on SAR Images
This experiment is focused on edge detection and edge extraction of landscape
elements in SAR images. The detection was performed on multiple SAR images including
TSX and Sentinel-1 radar satellite missions within the vicinity of Bratislava, the capital
city of Slovakia.
3.1. Area of Interest
The area of interest (AOI) is situated northeast of the capital city Bratislava next to
class II road 503 between the D1 highway and the city of Pezinok. This location was
selected due to multiple available SAR images of this area and for its good accessibility via
the D1 highway. The study area is located in the Danubian Lowland at the intersection of
Danubian Hills and Danubian Flat. The AOI's center geographic coordinates are 48° 14´
52´´ N and 17° 19´ 20´´ E in World geodetic system 1984 (WGS 84). According to the
territorial and administrative system of Slovakia it is situated in Bratislava Region, Pezinok
District and it falls under the cadastral areas of Slovenský Grob and Viničné.
The AOI (Figure 3.1) consists of an agricultural field with the size of 219 hectares.
The field is bounded on the northeast (NE) by an asphalt road, on the northwest (NW) by
an alley of trees, on the southwest (SW) by a field road and on the southeast (SE) side it is
bounded by various border types such as trees, field roads and more agricultural fields.
These different borders are essential for evaluating the SAR's and edge detector's ability to
recognize edges of landscape elements.
Figure 3.1: Area of interest. (Adapted from Google Earth)
31
3.2. TSX Imagery
The TerraSAR-X images analyzed in this experiment were obtained through
German Aerospace Center (DLR) project ID LAN1583: Object Recognition Based on
High-Resolution Radar Imagery and project ID LAN2833: Utilization of the High SpatioTemporal Resolution of TerraSAR-X Observations for Recent Ground Deformation
Monitoring and SAR Image Processing. More information about these projects is available
at the TerraSAR-X Science Service System. For the purposes of the experiment one image
from project LAN1583 and two images from LAN2833 were selected.
3.2.1. TSX Mission Description
TerraSAR-X is a joint project between the DLR and the German industry
(ASTRIUM). Operational data acquisition for the TerraSAR-X mission is executed by two
satellites TSX-1 (in orbit since June 2007) and TDX-1 (launch in June 2010). Both
instruments will also fulfill the TanDEM-X mission in a joint operation. DLR owns and
operates the satellites and the payload ground segment (PGS) and holds the rights for the
scientific exploitation of the TerraSAR-X mission data (Fritz, Eineder, 2013).
The TSX-1 is a side-looking X-band SAR based on active phased array antenna
technology. Figure 3.2 illustrates an artist's view of the satellite. Characteristic values of
the platform and SAR instrument are summarized in Table 3.1. In June 2010, TSX-1 was
supplemented in orbit by its twin, the TanDEM-X instrument.
Figure: 3.2: Artist's view of TerraSAR-X.
Adapted from (DLR, 2015)
32
Table 3.1: TerraSAR-X technical facts.
TerrSAR-X technical facts
5 years, for both satellites (TerraSAR-X and
TanDEM-X), an extended lifetime of that least
Operational life
another 5 years (beyond 2018) is expected by
the operator DLR (Status: April 2014).
Orbit
Sun-synchronous repeat orbit
Repeat period
11 days
Inclination
97.44°
Altitude at the equator
514 km (319.8 miles)
Centre Frequency
9.65 GHz (X band)
Nominal acquisition direction
Right side
Single, dual - depending on imaging mode
Polarization
quadruple is available as advanced polarization
mode for dedicated acquisition campaigns
(Adapted from TerraSAR-X Image Product Guide, 2014)
3.2.2. LAN 1583 Dataset
The TSX data was delivered as basic image products. They correspond to the
CEOS Level 1b quality. TSX imagery from project LAN 1583 was acquired with the
SpotLight (ST) imaging mode. Spotlight imaging modes uses phased array beam steering
in azimuth direction to increase the illumination time, i.e. the size of the synthetic aperture.
(TerraSAR-X Image Product Guide, 2014). Figure 3.3 shows the image subset of the AOI.
and Table 3.2 summarizes the basic metadata of the image.
Figure 3.3: Preview of the TSX image of the AOI from December 3, 2008.
33
Table 3.2: Basic metadata of the TSX image acquired on December 3, 2008.
Image metadata
Product type
Acquisition mode
Beams
Acquisition date
Cycle
Track
Orbit
Pass
Sample type
Polarization
Range sample spacing
Azimuth sample spacing
Pulse repetition frequency
Radar frequency
Total size
Raster height
Raster width
Average scene height ellipsoid
Map projection applied
SE_SL_S
Spotlight
Spot_066
03-Dec-2008
49
146
8162
Ascending
Detected
VV
0.75 m
0.75 m
3757 Mhz
9650 MHz
922 Mb
23333 lines
20666 samples
268.394 m
UTM zone 33N
3.2.3. LAN 2833 Dataset
LAN 2833 was acquired with the StripMap (SM) mode. Stripmap is the basic SAR
imaging mode known, e.g. ERS-1 and other radar satellites. The ground swath is
illuminated with continuous sequence of pulses while the antenna beam is fixed in
elevation and azimuth. This results in an image strip with a continuous image quality (in
flight direction) (TerraSAR-X Image Product Guide, 2014). Previews of the two images
used in the practical experiment are shown in Figure 3.4 and Figure 3.5. Table 3.3
summarizes the basic metadata of the images. The fist image from this dataset was
acquired five years after the initial image from LAN1583 and the second image after
another three months, both with HH polarization.
Figure 3.4: Preview of the TSX image of the AOI from January 21, 2013.
34
Figure 3.5: Preview of the TSX image of the AOI from June 2, 2013.
Table 3.3: Basic metadata of the TSX images from the LAN2833 dataset.
Image metadata
Product type
EEC_SE_SM_S
Acquisition mode
StripMap
Beams
Strip_007
Acquisition date
21-JAN-2013
Cycle
187
Track
17
Orbit
31079
Pass
Descending
Sample type
Detected
Polarization
HH
Range sample spacing
1.25 m
Azimuth sample spacing
1.25 m
Pulse repetition frequency
3437 Mhz
Radar frequency
9650 MHz
Total size
3439 Mb
Raster height
52400 lines
Raster width
34800 samples
Average Scene height ellipsoid
228.715 m
Map projection applied
UTM zone 33N
EEC_SE_SM_S
StripMap
Strip_007
02-JUN-2013
199
17
33083
Descending
Detected
HH
1.25 m
1.25 m
3437 Mhz
9650 MHz
3479 Mb
52400 lines
34800 samples
228.728 m
UTM zone 33N
3.3. Sentinel-1 Imagery
Sentinel 1 data was acquired from the Sentinel-1 Scientific Data Hub. The Sentinel
1 Scientific data Hub provides free and open access to a Rolling Archive of Sentinel-1
Level-0 and Level-1 use products. Access to the data hub is granted after free one time
registration.
35
3.3.1. Sentinel-1 Mission Description
Sentinel-1 is a polar-orbiting satellite for operational SAR applications. The
constellation of two C-band radar satellites will provide continuous all weather day/night
imagery for user services, especially those identified in ESA's GMES (Copernicus) service
elements programme and on projects funded by the European Union (EU) (ESA Bulletin,
2007). Sentinel-1 satellites are being built by an industrial consortium led by Thales Alenia
Space (Italy) as the prime Contractor, while ASTRIUM (Germany) is responsible for the
C-band SAR payload. Figure 3.6 illustrates an artist's view of the Sentinel-1 satellite.
Technical parameters of the satellite are summarized in Table 3.4.
Figure: 3.6: Artist's view of Sentinel-1.
Adapted from (www.esa.int, 2015)
Table 3.4: Sentinel-1 technical facts.
Sentinel-1 technical facts
Operational life
Orbit
7 years with consumables for 12 years
Near-polar, Sun-synchronous, circular
12 days repeat cycle (1 satellite), 6 days
Repeat period
(constellation)
Inclination
98.18°
Altitude at the equator
698 km (433.7 miles)
Centre Frequency
5.405 GHz (C band)
Nominal acquisition direction
Right side
Polarization
Single and dual polarization for all modes
(Adapted from Sentinel-1 User Handbook, 2014)
36
3.3.2. Sentinel 1 Dataset
The downloaded product specification is Level-1 Ground Range Detected (GRD)
acquired by the Interferometric Wide Swath Mode (IW). The product consists of focused
SAR data that has been detected, multi-looked and projected to ground range. Ground
range coordinates are the slant range coordinates projected onto the ellipsoid of the Earth.
Pixel values represent detected amplitude. Phase information is lost. The resulting product
has approximately square resolution pixels and square pixel spacing with reduced speckle
at a cost of reduced geometric resolution. (Sentinel Online, 2015). Tested data vas
available in both VV (Figure 3.7) and VH (Figure 3.8) polarization with 10 m spacing. The
data was acquired on March 7, 2015.
Figure 3.7: Sentinel-1 image of the AOI from March 7, 2015, VV polarization.
Figure 3.8: Sentinel-1 image of the AOI from March 7, 2015, VH polarization.
37
Table 3.5: Basic Metadata of the Sentinel-1 images acquired on March 7, 2015.
Image metadata
Product type
Acquisition mode
Beams
Processed time
Cycle
Track
Orbit
Pass
Sample type
Polarization
Range sample spacing
Azimuth sample spacing
Pulse repetition frequency
Radar frequency
Total size
Raster height
Raster width
Average scene height ellipsoid
Map projection applied
GRD
IW
07-MAR-2015
42
124
4921
Descending
Detected
VH, VV
10 m
9.996 m
5.40500045433435E9 MHz
3285 Mb
16697 lines
25764 samples
268.394 m
-
3.4. Edge Detection
Edge detection refers to operation or multiple operations performed on an image to
detect edges. Edges are significant local changes of intensity in an image. Usually they
characterize boundaries between two different regions and are therefore a problem of
fundamental importance in image processing.
Since edge detection has been studied, several techniques for edge enhancement
were proposed. Robust edge detection techniques are essentially based on the following
two steps: edge enhancement and decision. Unlike optical images, in SAR data, which is
highly heterogeneous, a robust edge enhancement phase is critical in providing acceptable
detection rates. This phase is usually performed through techniques that are related to
derivation, namely simple differences, Sobel filter (Pratt, 1978), Roberts filter, Prewitt
filter (Prewitt, 1970), morphological gradients etc., possibly combined with smoothing.
However there are three problems by the use of derivative as an edge detector (Prasad,
Rao, Chandrasekhar, 2013):
•
False edges produced by noise are enhanced,
•
ramp edges produce lower derivative magnitude,
•
the location of the ramp edges are not accurate.
38
Canny edge detection algorithm (Canny, 1986) provides tools to overcome these
edge detection problems by applying the derivative to edge detection.
Edge detection and edge extraction was achieved by MATLAB Image Processing
Toolbox. For the purpose of edge detection and extraction from georeferenced SAR
images in GeoTiff format a simple MATLAB function was proposed. The function is
called Saredge. This function enables reading images in GeoTiff format (both optical and
SAR), perform edge detection on greyscale images and export map coordinates of every
single pixel belonging to the detected edge in text file format.
The function stores both image data and a spatial referencing object. Edge
detection is executed by the built-in MATLAB Canny Edge Detection Algorithm. The
Canny edge detector first applies Gaussian smoothing on the image to eliminate remaining
noise. Secondly it computes the squared gradient magnitude. Then it applies nonmaximum suppression (NMS) to get rid of spurious response to edge detection. After NMS
it applies double thresholding to determine potential edges. Local maxima above the first
threshold are identified as weak edges and local maxima above the second threshold are
identified as strong edges. It finalizes the detection by hysteresis. Edges that are weak and
not connected to strong edges are suppressed.
Edge detection returns a binary image containing ''1'' where the function finds edges
in the input image and ''0'' elsewhere. Image coordinates of every pixel containing ''1'' are
then transformed to map coordinates using the stored reference object.
Saredge was gradually applied to all test SAR images previously smoothed by
Gamma MAP filter. The criterion was to detect maximum number of true edges while
detecting minimum number of false edges. Every image required a different threshold
setting to achieve the best possible results. The best results for each SAR image are shown
in figures 3.9 - 3.13.
39
Figure 3.9: Detected edges on the TSX image from December 3, 2008.
Figure 3.10: Detected edges on the TSX image from January 21, 2013.
40
Figure 3.11: Detected edges on the TSX image from June 2, 2013.
Figure 3.12: Detected edges on the Sentinel-1 image with VV polarization from
March 7, 2015.
41
Figure 3.13: Detected edges on the Sentinel-1 image with VH polarization from
March 7, 2015.
Edge detection on the TSX SAR image from December 2008 brought the best
results. All desired edges were identified correctly with a small number of false edges
identified in the homogenous areas due to remaining speckle noise.
Detection on the TSX image from January 2013 did not bring such effective
results. During the sensing period in January 2013 the AOI was apparently covered with
snow. Snow cover is a smooth surface and therefore most of the radiation was reflected
away from the SAR sensor. Because of this the image appears with dark tones and dull
edges, which were not recognisable by Saredge. On the other hand edges detected on the
TSX image from June 2013 are comparable with results obtained from December 2008. By
comparing figures 3.9 and 3.11 it is clearly visible that the original single agricultural field
is now divided into two new agricultural fields.
The attempt to detect edges in Sentinel-1 imagery met with the poorest results and
the edges were incomparable with the TSX dataset. The main drawback of the Sentinel-1
imagery was the lower resolution then the TSX imagery resolution. Also it is clearly
visible that the VH polarised image contains lower amplitudes. The best and only usable
results were achieved on the TSX SAR images from December 2008 and June 2013.
42
3.5. In Situ Measurements
For evaluation of edge detection accuracy it was necessary to perform in situ
measurements in the AOI. The borders of the AOI were mapped with Real Time
Kinematic (RTK) GNSS method on March 20, 2015 with three TOPCON GNSS receivers.
Two receivers TOPCON Hiper GD were used in base/rover setup and the receiver
TOPCON Hiper II was connected to the Slovakian Spatial Observation Service (SKPOS)
which represents active geodetic controls in Slovakia. The border was mapped with
approximately 50 meter step. Single types of borders are shown in Figure 3.14 and Figure
3.15 from (a) to (d). The point coordinates were observed in the geographic coordinate
system WGS 84 (EPSG: 4326). The provided SAR products were projected in the
coordinate system UTM zone 33N (EPSG: 32633). Projection of the observed points was
carried out in Arcmap. The observed points are displayed on the background of the TSX
image in Figure 3.15.
(a)
(b)
(c)
(d)
Figure 3.14: Representative images of the NW border (a), SW border (b),
SE border (c), NE border (d).
43
Figure 3.15:: Points observed with Real-Time
Real Time Kinematic method displayed on the background
of the TSX SAR image.
3.6. Comparison with In Situ Measurements
The points observed by GNSS served as a flawless reference for the detected edges.
The detected edges consist of many single points. These points represent pixels that were
identified as edges. The nearest detected points to the in situ observed points were
we selected
and the differences between them were computed. Note that the SAR images
image were acquired in
2008 and 2013 while the in situ measurements took place in 2015. Since 2008
20 the boundaries
of the AOI may have slightly changed due to agricultural activity. The differences between
GNSS observations and edges detected on the TSX imageries are displayed in
i figures 3.16 3.23. Overlay of the differences from both images is displayed in figures 3.24
3.2 - 3.27.
During the results analysis of the northeast border,, a slight systematic shift was
recorded in both SAR images. Five
Fiv meter shift in image from 2008 and 6 meter shift in image
from 2013. The shift may have been caused by the SAR imaging geometry and by the
microwave's physical nature. The main road is surrounded by trees and ditches on both sides.
side
One possible explanation is that the detected amplitude may have been detected from the
treetops. These systematic shifts were eliminated and the results are shown in figures 3.16
3.1 and
3.17.
44
A slight shift can be also observed on the northwest border shown in figures 3.18 and
3.19. This border consists from alley of trees and this shift may have been caused by the
subjective determination of the border during in situ measurements.
The southeast and southwest borders consisted of different land cover types varying
from field roads, trees and agricultural fields. The larger differences on the southeast and
southwest may have been caused by information loss durinng speckle reduction. We can not
also exclude that the agricultural fileds were differently plowed, however it is highly unlikely
that the differences were so significant.
The results imply that edge detection on SAR images can't compete with automatic
detection on optical remote sensing images due to the distortive effect of the speckle noise.
However it can still find its use in aplications that don't require fine accuracy such as
agriculture or forestry. Such an example can be provided by figures 3.9 and 3.11. The first
image from December 2008 shows an uncultivated single large field. The second image from
June 2013 shows that the field was divided into two sepperated fields. By utilizing multi
polarised datasets in radar polarimetry it is also possible to identify the crops.
Another possible application that doesn't require fine accuracy is monitoring river
banks, floads, icerberg shifts or even pre-earthquake and post-earthquake evaulation. Also
very commonly are detected SAR images utilized in oil spill detection or ship detection.
NE Border, December 2008
1.5
Difference (m)
1
0.5
0
-0.5
-1
-1.5
-2
-2.5
1
3
5
7
9
11
13
15
17
19
21
23
25
27
29
Point ID
Figure 3.16: Differences between observed points and detected edges of the NE border from
December 2008, Mean = 0.00 m, STD = 0.75 m.
45
NE Border, June 2013
5
4
Difference (m)
3
2
1
0
-1
-2
-3
-4
-5
1
3
5
7
9
11
13
15
17
19
21
23
25
27
29
Point ID
Figure 3.17: Differences between observed points and detected edges of the NE border from
June 2013, Mean = 0.00 m, STD = 2.46 m.
NW Border, December 2008
10
Difference (m)
8
6
4
2
0
-2
-4
1
3
5
7
9
11
13
15
17
19
21
23
25
27
29
31
33
Point ID
Figure 3.18: Differences between observed points and detected edges of the NW border from
December 2008, Mean = 3.62 m, STD = 1.98 m.
46
NW Border, June 2013
20
Difference (m)
15
10
5
0
-5
1
3
5
7
9
11
13
15
17
19
21
23
25
27
29
31
33
Point ID
Figure 3.19: Differences between observed points and detected edges of the NW border from
June 2013, Mean = 5.65 m, STD = 4.36 m.
SE Border, December 2008
8
Difference (m)
6
4
2
0
-2
-4
-6
-8
-10
1
3
5
7
9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47
Point ID
Figure 3.20: Differences between observed points and detected edges of the SE border from
December 2008, Mean = -1.61 m, STD = 3.48 m.
47
SE Border, June 2013
20
15
Difference (m)
10
5
0
-5
-10
-15
-20
-25
1
3
5
7
9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47
Point ID
Figure 3.21: Differences between observed points and detected edges of the SE border from
June 2013, Mean = -4.30 m, STD = 9.39 m.
SW Border, December 2008
6
5
Difference (m)
4
3
2
1
0
-1
-2
-3
1
3
5
7
9
11
13
15
17
19
21
23
25
27
29
31
33
Point ID
Figure 3.22: Differences between observed points and detected edges of the SW border from
December 2008, Mean = 1.45 m, STD = 1.77 m.
48
SW Border, June 2013
15
Difference (m)
10
5
0
-5
-10
-15
-20
1
3
5
7
9
11
13
15
17
19
21
23
25
27
29
31
33
Point ID
Figure 3.23: Differences between observed points and detected edges of the SW border from
June 2013, Mean = -1.60 m, STD = 4.94 m.
NE Border
TSX December 3, 2008
TSX June 2, 2013
5
4
3
Difference (m)
2
1
0
-1
-2
-3
-4
-5
1
3
5
7
9
11
13
15
17
19
21
23
25
Point ID
Figure 3.24: Overlay of the differences between observed points
and detected edges of the NE border.
49
27
29
NW Border
TSX December 3, 2008
TSX June 2, 2013
20
Difference (m)
15
10
5
0
-5
1
3
5
7
9
11
13
15
17
19
21
23
25
27
29
31
33
Point ID
Figure 3.25: Overlay of the differences between observed points
and detected edges of the NW border.
SE Border
TSX December 3, 2008
TSX June 2, 2013
20
15
Difference (m)
10
5
0
-5
-10
-15
-20
-25
1
3
5
7
9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47
Point ID
Figure 3.26: Overlay of the differences between observed points
and detected edges of the SE border.
50
SW Border
TSX December 3, 2008
TSX June 2, 2013
15
10
Difference (m)
5
0
-5
-10
-15
-20
1
3
5
7
9
11
13
15
17
19
21
23
25
27
29
Point ID
Figure 3.27: Overlay of the differences between observed points
and detected edges of the SW border.
51
31
33
Conclusion
The goal of this thesis was to test object detection on SAR images. Due to the physical
nature of microwave radiation and SAR imaging properties all radar images are inherently
degraded by speckle noise. Five adaptive speckle filters were explored to find the most
acceptable one for denoising speckle and preparing SAR images for edge detection.
The main criterium was to reduce speckle noise and to preserve sharpness of the edges
as much as possible. Quantitative assessment showed interesting results. According to ENL,
SSI and SMPI measures the best speckle reduction was achieved by Gamma MAP filter, note
that the differences between Gamma MAP, Frost and Lee filter where negligible. The filters
with the best edge preservation ability according to EEI and FPI were Refined Lee and
Perona-Malik. Practical edge detection questioned these results. Best edge detection results
were achieved on images that have been previously smoothed by Gamma MAP filter however
results on different SAR datasets may vary. The best filter depends on many parameters and
variables of the SAR system and the observed surface.
A simple MATLAB function for extracting edges from images in GeoTiff format was
proposed. The function was tested on TSX and Sentinel-1 datasets. The results have shown
that the main prerequisite for edge detection is an effective speckle reduction filter and high
resolution of the SAR images. While the TSX datasets were available in 0.75 m and 1.25 m
sample spacing, the freely available Sentinel-1 dataset had only 10 m sample spacing. The
results of the edge detection between TSX and Sentinel-1 were incomparable. The detected
edges from TSX images were compared with edges mapped by GNSS. The obtained results
can't compete with automatic detection on optical images, however due to the SAR image
degradation by speckle noise, but the results are still acceptable for certain applications like
agriculture or forestry, monitoring large scale events such as earthquakes, floods or oil spill.
In addition SAR imagery carries more information that can be utilized besides object
detection in different applications such as radar interferometry or radar polarimetry.
52
Resume
V súčasnosti náš svet čelí zvýšenej hrozbe prírodných katastrof zapríčinených
klimatickými zmenami. Tieto zmeny sú podmienené prirodzeným alebo antropogénnym
pôvodom. Z toho dôvodu je v záujme mnohých svetových organizácií, ako je napr. Európska
vesmírna agentúra (ESA),
monitorovať
stav
životného
prostredia
na pevnine,
v oceánoch a v atmosfére. Cieľom je zaistiť dáta pre zistenie aktuálneho stavu životného
prostredia, predikciu vývoja klimatických zmien a tiež poskytnúť dáta pre krízový manažment
za účelom zvýšenia bezpečnosti obyvateľstva.
Účinným nástrojom na monitorovanie životného prostredia Zeme sa stal diaľkový
prieskum Zeme s využitím radaru so syntetickou apertúrou (SAR). Radar na rozdiel
od optických senzorov je aktívny senzor vybavený vlastným zdrojom mikrovlnného žiarenia.
Radar ožaruje oblasť záujmu a po interakcii žiarenia s cieľom zaznamenáva odrazené
žiarenie. Produktom spracovania tohto odrazeného signálu je 2D obraz snímanej oblasti.
Vlastný zdroj žiarenia a použitie vĺn z mikrovlnnej časti spektra umožňuje radaru
snímkovať zemský povrch cez deň, v noci a bez závislosti od počasia. Jediná radarová snímka
má tak potenciál zachytiť aktuálny stav ľadovcov alebo škody napáchané zemetrasením,
či lesnými požiarmi na území väčšom ako je 1000 km2. Dlhodobé kontinuálne snímkovanie
umožňuje
sledovanie
environmentálnych
zmien
a porovnaním
aktuálnych
snímok
s archivovanými snímkami sa dá zistiť napr. rozsah výrubu pralesov alebo nárast hladiny
oceánov.
Hlavným cieľom tejto práce je testovanie využitia družicových radarových snímok
na detekciu prvkov krajiny a ich zmien v čase. Hlavnou prekážkou je, že všetky radarové
snímky sú znehodnotené tzv. speckle šumom, ktorý vzniká interferenciou odrazeného
radarového signálu. Tento šum sťažuje akúkoľvek interpretáciu radarových snímok vrátane
detekcie objektov. Z tohto dôvodu musia byť pred extrakciou nových informácií radarové
snímky odšumené. K tomuto účelu práca obsahuje kvantitatívne porovnanie a zhodnotenie
účinnosti Frostovho filtra, Leeho filtra, vylepšeného Leeho filtra, Gamma MAP filtra
a vyhladzovacieho filtra založeného na Perona-Malikovej rovnici. Okrem jednotlivých filtrov
boli testované aj rôzne veľkosti másk filtrov a tiež aj viacero iterácií filtrovania.
Na posúdenie účinnosti filtrov bolo vybratých päť testovacích kritérií: ENL, SSI,
SMPI, EEI a FPI. ENL, SSI a SMPI indexy slúžili na zhodnotenie schopnosti redukovať
speckle šum na veľkých homogénnych oblastiach a pritom zachovávať rádiometrické
vlastnosti radarových snímok. EEI a FPI indexy slúžili na porovnanie schopnosti zachovávať
53
hrany. Manuálny výpočet pre každú snímku by bol zdĺhavý a s väčšou pravdepodobnosťou
omylu, preto bola navrhnutá jednoduchú funkciu Speckleeval, ktorá dokázala načítať spolu
s originálnou snímkou aj vyhladenú snímku, vypočítať všetky testovacie parametre a zapísať
ich do súboru XLS.
V zásade platilo, že s pribúdajúcim množstvom iterácií sa zlepšovala redukcia speckle
šumu, ale zhoršovala sa schopnosť zachovávania hrán. Najlepšie výsledky v redukcii speckle
šumu boli dosiahnuté Gamma MAP filtrom. Zároveň treba spomenúť, že po piatich iteráciách
filtrovania, boli výsledky Leeho a Frostovho filtra takmer totožné. Najlepšie zachovávanie
hrán podľa EEI a FPI indexov dosiahli Leeho vylepšený filter a Perona-Malikov filter.
Aplikovanie hranového detektora na vyhladené snímky prinieslo ďalšie výsledky. Detegované
hrany na snímkach vyhladených Gamma MAP filtrom boli oveľa hladšie a neprerušované ako
napr. pri Perona-Malikovom filtri. Z tohto dôvodu bol na odšumenie a pripravenie snímok na
detekciu hrán v záujmovom území zvolený Gamma MAP filter.
Záujmové územie je situované severovýchodne od hlavného mesta Slovenskej
republiky, Bratislavy. Záujmovú oblasť tvorí poľnohospodárske pole s rozlohou 219 hektárov.
Toto pole je na severovýchode ohraničené asfaltovou cestou, na severozápade stromoradím,
na juhozápade poľnou cestou a juhovýchodnú hranicu poľa tvorí stromoradie, poľná cesta
a ďalšie poľnohospodárske polia. Táto lokalita bola vybratá z dôvodu dobrého prístupu,
vhodným rôznorodým typom ohraničenia a zároveň veľkého množstva radarových snímok
z tohto územia.
Radarové dáta poskytnuté k experimentu pochádzali z misie TerraSAR-X (TSX)
a novej misie ESA, Sentinel-1. Snímky TerraSAR-X boli poskytnuté v rámci projektov
LAN1583 a LAN2833. Z projektu LAN1583 bola vybratá radarová snímka záujmového
územia zachytená 3. decembra 2008 s VV polarizáciou a so štvorcovým rozlíšením pixla
s rozmerom 0.75 m. Zo snímok určených pre projekt LAN 2833 boli vybraté dve snímky
z 21. januára 2013 a 2. júna 2013 s HH polarizáciou a rozlíšením pixla 1.25m
Dáta z družice Sentinel-1 sú voľne k dispozícií na portáli Sentinel-1 Scientific Data
Hub. Použitý produkt vo formáte GRD bol dostupný vo VV aj VH polarizácii so štvorcovým
rozlíšením pixla 10 m. Snímky boli zachytené 7. marca 2015.
Pre účel detekcie a extrakcie hrán bola vytvorená funkcia s názvom Saredge. Funkcia
dokáže načítať radarovú (alebo aj optickú) snímku vo formáte GeoTiff, vykonať detekciu
hrán s nastavením rôznych prahov, pretransformovať obrazové súradnice na mapové
súradnice za pomoci referencie uloženej v pôvodnej snímke a výsledok zapísať do textového
súboru. Detekcia hrán prebehla na všetkých uvedených radarových snímkach. Kritériom
54
bolo zaznamenať čo najviac skutočných hrán s súčasne minimalizovať detegované falošné
hrany.
Najlepšie výsledky boli dosiahnuté na snímke TSX z decembra 2008. Všetky dôležité
hrany boli zaznamenané s minimálnym počtom falošných hrán. Snímka z januára 2013
nepriniesla tak dobré výsledky. V dobe snímkovania bola záujmová oblasť zrejme
pod snehovou pokrývkou. Sneh je hladký povrch, preto sa väčšina žiarenia odrazila smerom
od radaru. Z toho dôvodu snímka obsahuje nízke amplitúdy (tmavé pixle) a nevýrazné hrany.
Detekcia hrán na snímke TSX z júna 2013 priniesla porovnateľné výsledky so snímkou
z decembra 2012. Porovnaním snímok bol zistené, že pôvodne jedno veľké poľnohospodárske
pole bolo rozdelené na dve menšie. Pokus detekcie na snímkach Sentinel-1 priniesol
najslabšie výsledky, ktoré nie sú porovnateľné so snímkami z družice TSX. Hlavným
dôvodom tohto neúspechu bolo nižšie rozlíšenie snímok.
Kvôli porovnaniu presnosti detegovaných hrán bolo v záujmovom území vykonané
meranie prístrojmi GNSS kinematickou metódou v reálnom čase. Jednotlivé hrany
pozostávali s veľkého množstva bodov. Body hrán, ktoré boli najbližšie k odmeraným bodom
prístrojmi GNSS boli vybraté a ich rozdiely medzi odmeranými bodmi boli vypočítané.
Rozdiely medzi GNSS meraním a detegovanými hranami boli prezentované vo forme grafov.
Z výsledkov je zrejmé, že táto metóda sa nemôže rovnať presnosti detekcie prvkov na
optických snímkach s vysokým rozlíšením, ale stále je použiteľná pre aplikáciu
v poľnohospodárstve alebo v lesníctve. Okrem amplitúdovej zložky, radarové dáta obsahujú
aj fázovú zložku, ktorá sa dá použiť v ďalších aplikáciách ako je napr. radarová
interferometreia, čo robí z technológie SAR užitočný nástroj pri sledovaní našej planéty.
55
Bibliography
Austin Center for Space Research. http://www.csr.utexas.edu/ (2015-04-29).
ALI, M., CLAUSI, D.: Using The Canny Edge Detector for Feature Extraction and
Enhancement of Remote Sensing Images. Canada: University of Waterloo 2001.
ArcGIS 10.3 Help. http://resources.arcgis.com/en/help/ (2015-04-28).
CANNY, J.: A Computational Approach to Edge Detection in Pattern Analysis and Machine
Inteligence. United States of America: IEEE 1986.
European Space Agency. http://www.esa.int/ (2015-04-25).
ESA Bulletin 131.
http://www.esa.int/esapub/bulletin/bulletin131/bul131a_attema.pdf (2015-04-29).
FERRETI, A., MONTI-GUARNIERI, A., PRATI, C., ROCCA, F., MASSONET, D.: InSAR
Principles: Guidelines for SAR Interferometry Processing and Interpretation. The
Netherlands: ESA Publications 2007, ISBN 92-9092-233-8.
FRITZ, T., EINEDER, M.: TerraSAR-X Ground Segment Basic Product Specification
Document. Germany: DLR 2013.
Centre
for
Remote
Imaging,
Sensing
and
Processing
(CRISP).
http://www.crisp.nus.edu.sg/~research/tutorial/freqpol.htm (2015-04-15).
German Aerospace Center. http://www.dlr.de/ (2015-04-29).
Google Earth. https://www.google.com/earth/ (2015-03-15).
GONZALES, R., C.: Digital Image Processing. New Jersey: Pearson Education International
2008, ISBN 978-0-13-505267-9.
Slovakian InSAR team. http://insar.sk/ (2015-03-23).
LEE, J.-S.: Speckle Analysis and Smoothing of Synthetic Aperture Radar Images in
Computer Graphics and Image Processing. Washington, D.C.: Elsevier 1981.
56
LEE, J.-S., JURKEVICH, I., DEWAELE, P., WARMBACQ, P., OOSTERLINC,A.: Speckle
Filtering of Synthetic Aperture Radar Images: A Review. In: Remote Sensing Reviews.
United States of America: Harwood Academic Publishers 1994.
LEE, J.-S., POTTIER, L.: Polarimetric Radar Imaging: from Basics to Applications.
Boca Raton: CRC Press, Taylor & Francis Group 2009, ISBN 978-1-4200-5497-2.
LILLESAND, T. M., KIEFER, R. W., CHIPMAN, J. W.: Remote Sensing and Image
Interpretation. Hoboken: John Wiley & Sons 2008, ISBN: 978-471-15227-7.
MacDonald, Dettwiller and Associates Ltd. (MDA).
http://gs.mdacorporation.com/ (2015-04-17).
MANSOURPOUR, M., RAJABI, M. A., BLAIS, J. A. R.: Effects and Performance of
Speckle Noise Reduction Filters on Active Radar and SAR Images. Ankara: ISPRS 2006.
Matlab Help. http://www.mathworks.com/help/matlab/index.html (2015-03-13).
Matlab Image Processing Toolbox Help.
http://www.mathworks.com/help/images/ (2015-03-13).
Matlab Mapping Toolbox Help.
http://www.mathworks.com/help/map/index.html (2015-03-13).
Microsoft Office Excel Help.
https://support2.microsoft.com/find-solutions/office/excel (215-03-13).
PAPČO, J., BAKOŇ, M.: Object Recognition Based on High Resolution Radar Imagery.
German Aerospace Center (DLR), TSX-Archive-2012 project ID LAN1583 2012.
PAPČO, J.: Utilization of the High Spatio-Temporal Resolution of TerraSAR-X Observations
for Recent Ground Deformation Monitoring and SAR Image Processing. German
Aerospace Center (DLR), TSX-Archive-2015 project ID LAN2833 2015.
PERONA, P., MALIK, J.: Scale-Space and Edge Detection Using Anisotropic Diffusion in
Pattern Analysis and Machine Intelligence. USA: IEEE Transactions 1990.
57
PRASAD, P.S, RAO, B. P, CHANDRASEKHAR, P.N.R.L.: Enhanced Edge Detection
Technique for SAR Images. In: International Journal of Computer Applications. New
York: IJCA 2013.
PRATT, W.K.: Digital image Processing. New York: Wiley, 1978.
PREWITT, J.M.S: Object Enhancement and Extraction. In: Picture Processing nad
Psychopictorics. New York: Academic, 1970.
QIU, F., BERGLUND, J., JENSEN, J. R., THAKKAR, P., DIANWEI, R.: Speckle Noise
Reduction in SAR Imagery Using A Local Adaptive Median Filter. In: GIScience and
Remote Sensing. Maryland: V. H. Winston & Son 2004.
RANEY, R. K.: Radar Fundamentals: Technical Perspective. Chapter 2. In: Principles and
Applications of Imaging Radar, Manual of Remote Sensing. Toronto: John Wiley & Sons,
1998. ISBN 978-0-471-29406-1.
RICHARDS, J.A.: Remote Sensing with Imaging Radar. Berlin: Springer 2009,
ISBN 978-3-642-02019-3.
SARODE, M. V., DESHMUKH, P. R.: Reduction of Speckle Noise and Image Enhancement
of Images Using Filtering Technique. India: IJACT 2011.
Sentinel-1 Toolbox Help.
https://sentinel.esa.int/web/sentinel/toolboxes/sentinel-1 (2015-01-25).
Sentinel online: https://sentinel.esa.int (2015-01-25).
Sentinel-1 Scientific Data Hub: https://scihub.esa.int/ (2015-03-15).
Sentinel-1 User Handbook.
https://sentinel.esa.int/documents/247904/685163/Sentinel-1_User_Handbook
(2015-05-10)
SHAMSODDINI, A., TRINDER, J. C.: Image Texture Preservation in Speckle Noise
Suppression. In: ISPRS TC VII Symposium - 100 Years ISPRS. Vienna: IAPRS 2010.
58
SHANTHI, L., VALARMATHI, M., L.: SAR Image Despeckling using Possibilistic Fuzzy
C-means Clustering and Edge Detection in Bandelet Domain. London: Springer 2013.
SHENG, Y., XIA, Z.-G.: A comprehensive Evaluation of Filters for Radar Speckle
Suppression. Boston: University of Massachusetts at Boston 1996.
Slovakian Spatial Observation Service. http://www.skpos.gku.sk/ (2015-03-20).
Spatial Reference. http://spatialreference.org/ (2015-03-20).
TEMPFLI, K., HUURNEMAN, G. C., BAKKER, W. H., JANSSEN, L. L. F.: Principles of
Remote Sensing: An Introductory textbook. Enschede: ITC 2009.
ISBN 978-90-6164-270-1.
TerraSAR-X Image Product Guide, Basic and Enhanced Radar Satellite Imagery.
http://www2.geo-airbusds.com/files/pmedia/public/r459_9_201408_tsxx-itd-ma-0009_tsxproductguide_i2.00.pdf (2015-04-29).
TerraSAR-X Science Service System. http://sss.terrasar-x.dlr.de (2015-02-01).
Tutorial: Fundamentals of Remote Sensing. Canada Centre for Remote Sensing: Natural
Resources Canada (CCRS). http://www.nrcan.gc.ca/earth-sciences/geomatics/satelliteimagery-air-photos/satellite-imagery-products/educational-resources/9309 (2015-04-23).
59
Appendices
60
Appendix A
Speckle Evaluation Results Part 1
Frost filter evaluation results
NAME
MEAN
STD
VAR
ENL
SSI
SMPI
EEI
FPI
3Frost1
155.6345
47.8899
2293.4438
3.2263
0.6837
1.0052
0.8224
0.8988
3Frost2
155.5941
38.8406
1508.5917
4.9022
0.5546
0.8376
0.7014
0.8232
3Frost3
155.5833
33.5823
1127.7728
6.5567
0.4796
0.7294
0.6175
0.7713
3Frost4
155.5798
30.0370
902.2218
8.1954
0.4289
0.6538
0.5548
0.7130
3Frost5
155.5772
27.4432
753.1301
9.8175
0.3919
0.5984
0.5113
0.6675
5Frost1
155.4685
34.3807
1182.0294
6.2465
0.4913
0.8029
0.6214
0.7826
5Frost2
155.4273
25.6567
658.2652
11.2107
0.3668
0.6142
0.4756
0.6417
5Frost3
155.4191
21.4931
461.9534
15.9731
0.3073
0.5170
0.3956
0.5499
5Frost4
155.4156
18.9233
358.0923
20.6050
0.2705
0.4562
0.3476
0.4824
5Frost5
155.4123
17.1365
293.6583
25.1251
0.2450
0.4139
0.3242
0.4373
7Frost1
155.6699
27.0595
732.2166
10.1099
0.3862
0.5544
0.5040
0.7074
7Frost2
155.6372
19.2665
371.1975
19.9343
0.2750
0.4037
0.3557
0.5053
7Frost3
155.6294
15.9727
255.1257
29.0007
0.2280
0.3364
0.3195
0.4136
7Frost4
155.6248
14.0116
196.3238
37.6846
0.2000
0.2960
0.3095
0.3670
7Frost5
155.6221
12.6711
160.5576
46.0777
0.1809
0.2682
0.3064
0.3433
Lee filter evaluation results
NAME
MEAN
STD
VAR
ENL
SSI
SMPI
EEI
FPI
3Lee1
156.0630
47.0906
2217.5234
3.3551
0.6704
0.7012
0.8180
0.9010
3Lee2
156.0611
38.3620
1471.6450
5.0555
0.5461
0.5723
0.6964
0.8267
3Lee3
156.0599
33.2633
1106.4455
6.7241
0.4736
0.4968
0.6139
0.7661
3Lee4
156.0575
29.8075
888.4841
8.3733
0.4244
0.4462
0.5520
0.7077
3Lee5
156.0562
27.2668
743.4802
10.0063
0.3882
0.4087
0.5077
0.6614
5Lee1
156.1017
32.4273
1051.5328
7.0790
0.4615
0.4650
0.5949
0.7610
5Lee2
156.0967
24.7935
614.7196
12.1085
0.3529
0.3573
0.4569
0.6167
5Lee3
156.0956
23.3398
544.7467
13.6636
0.3322
0.3367
0.4330
0.5821
5Lee4
156.0922
20.1087
404.3604
18.4066
0.2862
0.2911
0.3705
0.5127
5Lee5
156.0886
17.9803
323.2927
23.0211
0.2559
0.2612
0.3348
0.4556
7Lee1
156.0987
24.5578
603.0831
12.3424
0.3495
0.3532
0.4371
0.6141
7Lee2
156.0931
18.3072
335.1531
22.2077
0.2606
0.2648
0.3348
0.4572
7Lee3
156.0882
15.4482
238.6475
31.1862
0.2199
0.2245
0.3122
0.3954
7Lee4
156.0851
13.6735
186.9637
39.8056
0.1946
0.1993
0.3064
0.3615
7Lee5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
61
Appendix B
Speckle Evaluation Results Part 1
Refined Lee filter evaluation results
NAME
MEAN
STD
VAR
ENL
SSI
SMPI
EEI
FPI
RefLee1
157.3106
43.6266
1903.2788
3.9718
0.6162
1.3668
0.6844
0.8410
RefLee2
157.5995
34.6859
1203.1088
6.3064
0.4890
1.2293
0.5682
0.7688
RefLee3
157.5961
27.1274
735.8940
10.3099
0.3824
0.9601
0.4536
0.6847
RefLee4
157.5578
25.0301
626.5037
12.1042
0.3530
0.8722
0.4212
0.6486
RefLee5
157.5211
23.4333
549.1184
13.8035
0.3305
0.8043
0.3978
0.6228
ENL
SSI
EEI
FPI
Gamma MAP filter evaluation results
NAME
MEAN
STD
VAR
SMPI
3Gamma1
156.0265
47.1540
2223.4971
3.3446
0.6715
0.7267
0.8180
0.9010
3Gamma2
156.0201
38.4182
1475.9576
5.0381
0.5471
0.5956
0.6964
0.8267
3Gamma3
156.0171
33.3139
1109.8189
6.6999
0.4744
0.5179
0.6139
0.7661
3Gamma4
156.0132
29.8555
891.3493
8.3417
0.4252
0.4657
0.5520
0.7077
3Gamma5
156.0103
27.3142
746.0678
9.9657
0.3890
0.4272
0.5077
0.6614
5Gamma1
156.0998
32.4315
1051.8019
7.0770
0.4616
0.4660
0.5949
0.7610
5Gamma2
156.0946
24.7967
614.8767
12.1051
0.3529
0.3581
0.4569
0.6167
5Gamma3
156.0913
20.9949
440.7865
16.8853
0.2988
0.3042
0.3881
0.5345
5Gamma4
156.0892
18.5933
345.7098
21.5284
0.2647
0.2700
0.3426
0.4679
5Gamma5
156.0874
16.8969
285.5053
26.0676
0.2405
0.2457
0.3226
0.4325
7Gamma1
156.0987
24.5578
603.0842
12.3424
0.3495
0.3532
0.4371
0.6141
7Gamma2
156.0931
18.3072
335.1548
22.2075
0.2606
0.2648
0.3348
0.4572
7Gamma3
156.0882
15.4482
238.6478
31.1862
0.2199
0.2245
0.3122
0.3954
7Gamma4
156.0851
13.6735
186.9640
39.8055
0.1946
0.1993
0.3064
0.3615
7Gamma5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
Perona-Malik filter evaluation results
NAME
MEAN
STD
VAR
ENL
SSI
SMPI
EEI
FPI
PM_1
153.3511
52.3490
2740.4207
2.6214
0.7584
2.8000
0.7140
0.8638
PM_2
151.7369
37.9394
1439.3985
4.8863
0.5555
2.9009
0.5489
0.6711
PM_3
150.8968
27.5230
757.5182
9.1822
0.4052
2.4335
0.4123
0.5443
PM_4
150.5074
20.8780
435.8890
15.8752
0.3082
1.9617
0.3357
0.4635
PM_5
150.3294
16.9216
286.3391
24.1094
0.2501
1.6328
0.3005
0.4045
62
Appendix C
Speckle Evaluation Results Sorted by ENL
NAME
MEAN
STD
VAR
ENL
SSI
SMPI
EEI
FPI
7Gamma5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
7Lee5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
7Frost5
155.6221
12.6711
160.5576
46.0777
0.1809
0.2682
0.3064
0.3433
7Lee4
156.0851
13.6735
186.9637
39.8056
0.1946
0.1993
0.3064
0.3615
7Gamma4
156.0851
13.6735
186.9640
39.8055
0.1946
0.1993
0.3064
0.3615
7Frost4
155.6248
14.0116
196.3238
37.6846
0.2000
0.2960
0.3095
0.3670
7Lee3
156.0882
15.4482
238.6475
31.1862
0.2199
0.2245
0.3122
0.3954
7Gamma3
156.0882
15.4482
238.6478
31.1862
0.2199
0.2245
0.3122
0.3954
7Frost3
155.6294
15.9727
255.1257
29.0007
0.2280
0.3364
0.3195
0.4136
PM_5
150.3294
16.9216
286.3391
24.1094
0.2501
1.6328
0.3005
0.4045
7Lee2
156.0931
18.3072
335.1531
22.2077
0.2606
0.2648
0.3348
0.4572
7Gamma2
156.0931
18.3072
335.1548
22.2075
0.2606
0.2648
0.3348
0.4572
7Frost2
155.6372
19.2665
371.1975
19.9343
0.2750
0.4037
0.3557
0.5053
PM_4
150.5074
20.8780
435.8890
15.8752
0.3082
1.9617
0.3357
0.4635
RefLee5
157.5211
23.4333
549.1184
13.8035
0.3305
0.8043
0.3978
0.6228
7Lee1
156.0987
24.5578
603.0831
12.3424
0.3495
0.3532
0.4371
0.6141
7Gamma1
156.0987
24.5578
603.0842
12.3424
0.3495
0.3532
0.4371
0.6141
RefLee4
157.5578
25.0301
626.5037
12.1042
0.3530
0.8722
0.4212
0.6486
RefLee3
157.5961
27.1274
735.8940
10.3099
0.3824
0.9601
0.4536
0.6847
7Frost1
155.6699
27.0595
732.2166
10.1099
0.3862
0.5544
0.5040
0.7074
PM_3
150.8968
27.5230
757.5182
9.1822
0.4052
2.4335
0.4123
0.5443
RefLee2
157.5995
34.6859
1203.1088
6.3064
0.4890
1.2293
0.5682
0.7688
PM_2
151.7369
37.9394
1439.3985
4.8863
0.5555
2.9009
0.5489
0.6711
RefLee1
157.3106
43.6266
1903.2788
3.9718
0.6162
1.3668
0.6844
0.8410
PM_1
153.3511
52.3490
2740.4207
2.6214
0.7584
2.8000
0.7140
0.8638
original
156.1093
70.2631
4936.9075
1.5079
-
-
-
-
The higher ENL, the better speckle reduction.
63
Appendix D
Speckle Evaluation Results Sorted by SSI
NAME
MEAN
STD
VAR
ENL
SSI
SMPI
EEI
FPI
original
156.1093
70.2631
4936.9075
1.5079
-
-
-
-
7Gamma5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
7Lee5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
7Frost5
155.6221
12.6711
160.5576
46.0777
0.1809
0.2682
0.3064
0.3433
7Lee4
156.0851
13.6735
186.9637
39.8056
0.1946
0.1993
0.3064
0.3615
7Gamma4
156.0851
13.6735
186.9640
39.8055
0.1946
0.1993
0.3064
0.3615
7Frost4
155.6248
14.0116
196.3238
37.6846
0.2000
0.2960
0.3095
0.3670
7Lee3
156.0882
15.4482
238.6475
31.1862
0.2199
0.2245
0.3122
0.3954
7Gamma3
156.0882
15.4482
238.6478
31.1862
0.2199
0.2245
0.3122
0.3954
7Frost3
155.6294
15.9727
255.1257
29.0007
0.2280
0.3364
0.3195
0.4136
PM_5
150.3294
16.9216
286.3391
24.1094
0.2501
1.6328
0.3005
0.4045
7Lee2
156.0931
18.3072
335.1531
22.2077
0.2606
0.2648
0.3348
0.4572
7Gamma2
156.0931
18.3072
335.1548
22.2075
0.2606
0.2648
0.3348
0.4572
7Frost2
155.6372
19.2665
371.1975
19.9343
0.2750
0.4037
0.3557
0.5053
PM_4
150.5074
20.8780
435.8890
15.8752
0.3082
1.9617
0.3357
0.4635
RefLee5
157.5211
23.4333
549.1184
13.8035
0.3305
0.8043
0.3978
0.6228
7Lee1
156.0987
24.5578
603.0831
12.3424
0.3495
0.3532
0.4371
0.6141
7Gamma1
156.0987
24.5578
603.0842
12.3424
0.3495
0.3532
0.4371
0.6141
RefLee4
157.5578
25.0301
626.5037
12.1042
0.3530
0.8722
0.4212
0.6486
RefLee3
157.5961
27.1274
735.8940
10.3099
0.3824
0.9601
0.4536
0.6847
7Frost1
155.6699
27.0595
732.2166
10.1099
0.3862
0.5544
0.5040
0.7074
PM_3
150.8968
27.5230
757.5182
9.1822
0.4052
2.4335
0.4123
0.5443
RefLee2
157.5995
34.6859
1203.1088
6.3064
0.4890
1.2293
0.5682
0.7688
PM_2
151.7369
37.9394
1439.3985
4.8863
0.5555
2.9009
0.5489
0.6711
RefLee1
157.3106
43.6266
1903.2788
3.9718
0.6162
1.3668
0.6844
0.8410
PM_1
153.3511
52.3490
2740.4207
2.6214
0.7584
2.8000
0.7140
0.8638
The lower SSI, the better speckle reduction.
64
Appendix E
Speckle Evaluation Results Sorted by SMPI
NAME
MEAN
STD
VAR
ENL
SSI
SMPI
EEI
FPI
original
156.1093
70.2631
4936.9075
1.5079
-
-
-
-
7Gamma5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
7Lee5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
7Lee4
156.0851
13.6735
186.9637
39.8056
0.1946
0.1993
0.3064
0.3615
7Gamma4
156.0851
13.6735
186.9640
39.8055
0.1946
0.1993
0.3064
0.3615
7Lee3
156.0882
15.4482
238.6475
31.1862
0.2199
0.2245
0.3122
0.3954
7Gamma3
156.0882
15.4482
238.6478
31.1862
0.2199
0.2245
0.3122
0.3954
7Lee2
156.0931
18.3072
335.1531
22.2077
0.2606
0.2648
0.3348
0.4572
7Gamma2
156.0931
18.3072
335.1548
22.2075
0.2606
0.2648
0.3348
0.4572
7Frost5
155.6221
12.6711
160.5576
46.0777
0.1809
0.2682
0.3064
0.3433
7Frost4
155.6248
14.0116
196.3238
37.6846
0.2000
0.2960
0.3095
0.3670
7Frost3
155.6294
15.9727
255.1257
29.0007
0.2280
0.3364
0.3195
0.4136
7Lee1
156.0987
24.5578
603.0831
12.3424
0.3495
0.3532
0.4371
0.6141
7Gamma1
156.0987
24.5578
603.0842
12.3424
0.3495
0.3532
0.4371
0.6141
7Frost2
155.6372
19.2665
371.1975
19.9343
0.2750
0.4037
0.3557
0.5053
7Frost1
155.6699
27.0595
732.2166
10.1099
0.3862
0.5544
0.5040
0.7074
RefLee5
157.5211
23.4333
549.1184
13.8035
0.3305
0.8043
0.3978
0.6228
RefLee4
157.5578
25.0301
626.5037
12.1042
0.3530
0.8722
0.4212
0.6486
RefLee3
157.5961
27.1274
735.8940
10.3099
0.3824
0.9601
0.4536
0.6847
RefLee2
157.5995
34.6859
1203.1088
6.3064
0.4890
1.2293
0.5682
0.7688
RefLee1
157.3106
43.6266
1903.2788
3.9718
0.6162
1.3668
0.6844
0.8410
PM_5
150.3294
16.9216
286.3391
24.1094
0.2501
1.6328
0.3005
0.4045
PM_4
150.5074
20.8780
435.8890
15.8752
0.3082
1.9617
0.3357
0.4635
PM_3
150.8968
27.5230
757.5182
9.1822
0.4052
2.4335
0.4123
0.5443
PM_1
153.3511
52.3490
2740.4207
2.6214
0.7584
2.8000
0.7140
0.8638
PM_2
151.7369
37.9394
1439.3985
4.8863
0.5555
2.9009
0.5489
0.6711
The lower SMPI, the better speckle reduction.
65
Appendix F
Speckle Evaluation Results Sorted by EEI
NAME
MEAN
STD
VAR
ENL
SSI
SMPI
EEI
FPI
PM_1
153.3511
52.3490
2740.4207
2.6214
0.7584
2.8000
0.7140
0.8638
RefLee1
157.3106
43.6266
1903.2788
3.9718
0.6162
1.3668
0.6844
0.8410
RefLee2
157.5995
34.6859
1203.1088
6.3064
0.4890
1.2293
0.5682
0.7688
PM_2
151.7369
37.9394
1439.3985
4.8863
0.5555
2.9009
0.5489
0.6711
7Frost1
155.6699
27.0595
732.2166
10.1099
0.3862
0.5544
0.5040
0.7074
RefLee3
157.5961
27.1274
735.8940
10.3099
0.3824
0.9601
0.4536
0.6847
7Lee1
156.0987
24.5578
603.0831
12.3424
0.3495
0.3532
0.4371
0.6141
7Gamma1
156.0987
24.5578
603.0842
12.3424
0.3495
0.3532
0.4371
0.6141
RefLee4
157.5578
25.0301
626.5037
12.1042
0.3530
0.8722
0.4212
0.6486
PM_3
150.8968
27.5230
757.5182
9.1822
0.4052
2.4335
0.4123
0.5443
RefLee5
157.5211
23.4333
549.1184
13.8035
0.3305
0.8043
0.3978
0.6228
7Frost2
155.6372
19.2665
371.1975
19.9343
0.2750
0.4037
0.3557
0.5053
PM_4
150.5074
20.8780
435.8890
15.8752
0.3082
1.9617
0.3357
0.4635
7Lee2
156.0931
18.3072
335.1531
22.2077
0.2606
0.2648
0.3348
0.4572
7Gamma2
156.0931
18.3072
335.1548
22.2075
0.2606
0.2648
0.3348
0.4572
7Frost3
155.6294
15.9727
255.1257
29.0007
0.2280
0.3364
0.3195
0.4136
7Lee3
156.0882
15.4482
238.6475
31.1862
0.2199
0.2245
0.3122
0.3954
7Gamma3
156.0882
15.4482
238.6478
31.1862
0.2199
0.2245
0.3122
0.3954
7Frost4
155.6248
14.0116
196.3238
37.6846
0.2000
0.2960
0.3095
0.3670
7Frost5
155.6221
12.6711
160.5576
46.0777
0.1809
0.2682
0.3064
0.3433
7Lee4
156.0851
13.6735
186.9637
39.8056
0.1946
0.1993
0.3064
0.3615
7Gamma4
156.0851
13.6735
186.9640
39.8055
0.1946
0.1993
0.3064
0.3615
7Gamma5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
7Lee5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
PM_5
150.3294
16.9216
286.3391
24.1094
0.2501
1.6328
0.3005
0.4045
original
156.1093
70.2631
4936.9075
1.5079
-
-
-
-
The higher EEI, the better edge preservation.
66
Appendix G
Speckle Evaluation Results Sorted by FPI
NAME
MEAN
STD
VAR
ENL
SSI
SMPI
EEI
FPI
PM_1
153.3511
52.3490
2740.4207
2.6214
0.7584
2.8000
0.7140
0.8638
RefLee1
157.3106
43.6266
1903.2788
3.9718
0.6162
1.3668
0.6844
0.8410
RefLee2
157.5995
34.6859
1203.1088
6.3064
0.4890
1.2293
0.5682
0.7688
7Frost1
155.6699
27.0595
732.2166
10.1099
0.3862
0.5544
0.5040
0.7074
RefLee3
157.5961
27.1274
735.8940
10.3099
0.3824
0.9601
0.4536
0.6847
PM_2
151.7369
37.9394
1439.3985
4.8863
0.5555
2.9009
0.5489
0.6711
RefLee4
157.5578
25.0301
626.5037
12.1042
0.3530
0.8722
0.4212
0.6486
RefLee5
157.5211
23.4333
549.1184
13.8035
0.3305
0.8043
0.3978
0.6228
7Lee1
156.0987
24.5578
603.0831
12.3424
0.3495
0.3532
0.4371
0.6141
7Gamma1
156.0987
24.5578
603.0842
12.3424
0.3495
0.3532
0.4371
0.6141
PM_3
150.8968
27.5230
757.5182
9.1822
0.4052
2.4335
0.4123
0.5443
7Frost2
155.6372
19.2665
371.1975
19.9343
0.2750
0.4037
0.3557
0.5053
PM_4
150.5074
20.8780
435.8890
15.8752
0.3082
1.9617
0.3357
0.4635
7Lee2
156.0931
18.3072
335.1531
22.2077
0.2606
0.2648
0.3348
0.4572
7Gamma2
156.0931
18.3072
335.1548
22.2075
0.2606
0.2648
0.3348
0.4572
7Frost3
155.6294
15.9727
255.1257
29.0007
0.2280
0.3364
0.3195
0.4136
PM_5
150.3294
16.9216
286.3391
24.1094
0.2501
1.6328
0.3005
0.4045
7Lee3
156.0882
15.4482
238.6475
31.1862
0.2199
0.2245
0.3122
0.3954
7Gamma3
156.0882
15.4482
238.6478
31.1862
0.2199
0.2245
0.3122
0.3954
7Frost4
155.6248
14.0116
196.3238
37.6846
0.2000
0.2960
0.3095
0.3670
7Lee4
156.0851
13.6735
186.9637
39.8056
0.1946
0.1993
0.3064
0.3615
7Gamma4
156.0851
13.6735
186.9640
39.8055
0.1946
0.1993
0.3064
0.3615
7Frost5
155.6221
12.6711
160.5576
46.0777
0.1809
0.2682
0.3064
0.3433
7Gamma5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
7Lee5
156.0827
12.4316
154.5449
48.1541
0.1770
0.1816
0.3036
0.3282
The higher FPI, the better edge preservation.
67