imageryintro.1 9.39 KB
Newer Older
xuebingbing's avatar
xuebingbing committed
1
.TH imageryintro 1 "" "GRASS 7.8.3" "GRASS GIS User's Manual"
xuebingbing's avatar
xuebingbing committed
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252
.SH Image processing in GRASS GIS
.SS Image processing in general
\fBDigital numbers and physical values (reflection/radiance\-at\-sensor):\fR
.PP
Satellite imagery is commonly stored in Digital Numbers (DN) for
minimizing the storage volume, i.e. the originally sampled analog
physical value (color, temperature, etc) is stored a discrete
representation in 8\-16 bits. For example, Landsat data are stored in
8bit values (i.e., ranging from 0 to 255); other satellite data may
be stored in 10 or 16 bits. Having data stored in DN, it implies
that these data are not yet the observed ground reality. Such data
are called \(dqat\-satellite\(dq, for example the amount of energy sensed
by the sensor of the satellite platform is encoded in 8 or more
bits. This energy is called radiance\-at\-sensor. To obtain physical
values from DNs, satellite image providers use a linear transform
equation (y = a * x + b) to encode the radiance\-at\-sensor
in 8 to 16 bits. DNs can be turned back into physical values by
applying the reverse formula (x = (y \- b) / a).
.PP
The GRASS GIS module i.landsat.toar
easily transforms Landsat DN to radiance\-at\-sensor (top of atmosphere,
TOA). The equivalent module for ASTER data is
i.aster.toar.
For other satellites, r.mapcalc can
be employed.
.PP
\fBReflection/radiance\-at\-sensor and surface reflectance\fR
.PP
When radiance\-at\-sensor has been obtained, still the atmosphere
influences the signal as recorded at the sensor. This atmospheric
interaction with the sun energy reflected back into space by
ground/vegetation/soil needs to be corrected. The need of
removing atmospheric artifacts stems from the fact that the
atmosphericic conditions are changing over time. Hence, to gain
comparability between Earth surface images taken at different
times, atmospheric need to be removed converting at\-sensor values
which are top of atmosphere to surface reflectance values.
.PP
In GRASS GIS, there are two ways to apply atmospheric correction for
satellite imagery. A simple, less accurate way for Landsat is with
i.landsat.toar,
using the DOS correction method. The more accurate way is using
i.atcorr (which supports many satellite
sensors). The atmospherically corrected sensor data represent
surface reflectance,
which ranges theoretically from 0% to 100%. Note that this level of
data correction is the proper level of correction to calculate
vegetation indices.
.PP
In GRASS GIS, image data are identical to raster data.
However, a couple of commands are explicitly dedicated to image
processing. The geographic boundaries of the raster/imagery file are
described by the north, south, east, and west fields. These values
describe the lines which bound the map at its edges. These lines do
NOT pass through the center of the grid cells at the edge of the
map, but along the edge of the map itself.
.PP
As a general rule in GRASS:
.IP
.IP \fB1\fR
Raster/imagery output maps have their bounds and resolution equal
to those of the current region.
.IP \fB2\fR
Raster/imagery input maps are automatically cropped/padded and
rescaled (using nearest\-neighbor resampling) to match the current
region.
.PP
.SS Imagery import
The module r.in.gdal offers a common
interface for many different raster and satellite image
formats. Additionally, it also offers options such as on\-the\-fly
location creation or extension of the default region to match the
extent of the imported raster map.  For special cases, other import
modules are available. Always the full map is imported. Imagery data
can be group (e.g. channel\-wise) with i.group.
.PP
For importing scanned maps, the user will need to create a
x,y\-location, scan the map in the desired resolution and save it into
an appropriate raster format (e.g. tiff, jpeg, png, pbm) and then use
r.in.gdal to import it. Based on
reference points the scanned map can be rectified to obtain geocoded
data.
.SS Image processing operations
GRASS raster/imagery map processing is always performed in the current
region settings (see g.region), i.e. the
current region extent and current raster resolution is used. If the
resolution differs from that of the input raster map(s), on\-the\-fly
resampling is performed (nearest neighbor resampling). If this is not
desired, the input map(s) has/have to be resampled beforehand with one
of the dedicated modules.
.SS Geocoding of imagery data
GRASS is able to geocode raster and image data of various types:
.RS 4n
.IP \(bu 4n
unreferenced scanned maps by defining four corner points
(i.group, i.target,
g.gui.gcp, i.rectify)
.IP \(bu 4n
unreferenced satellite data from optical and Radar sensors by
defining a certain number of ground control points
(i.group, i.target,
g.gui.gcp, i.rectify)
.IP \(bu 4n
interactive graphical Ground Control Point (GCP) manager
.IP \(bu 4n
orthophoto generation based on DEM: i.ortho.photo
.IP \(bu 4n
digital handheld camera geocoding: modified procedure for
i.ortho.photo
.RE
.SS Visualizing (true) color composites
To quickly combine the first three channels to a near natural color
image, the GRASS command d.rgb can be used or
the graphical GIS manager (wxGUI). It assigns
each channel to a color which is then mixed while displayed. With a
bit more work of tuning the grey scales of the channels, nearly
perfect colors can be achieved. Channel histograms can be shown with
d.histogram.
.SS Calculation of vegetation indices
An example for indices derived from multispectral data is the NDVI
(normalized difference vegetation index). To study the vegetation
status with NDVI, the Red and the Near Infrared channels (NIR) are
taken as used as input for simple map algebra in the GRASS command
r.mapcalc
(ndvi = 1.0 * (nir \- red)/(nir + red)). With
r.colors an optimized \(dqndvi\(dq color table
can be assigned afterward. Also other vegetation indices can be
generated likewise.
.SS Calibration of thermal channel
The encoded digital numbers of a thermal infrared channel can be
transformed to degree Celsius (or other temperature units) which
represent the temperature of the observed land surface. This requires
a few algebraic steps with r.mapcalc
which are outlined in the literature to apply gain and bias values
from the image metadata.
.SS Image classification
Single and multispectral data can be classified to user defined land
use/land cover classes. In case of a single channel, segmentation will
be used.
GRASS supports the following methods:
.RS 4n
.IP \(bu 4n
Radiometric classification:
.RS 4n
.IP \(bu 4n
Unsupervised classification (i.cluster,
i.maxlik) using the Maximum Likelihood
classification method
.IP \(bu 4n
Supervised classification (i.gensig
or g.gui.iclass, i.maxlik)
using the Maximum Likelihood classification method
.RE
.IP \(bu 4n
Combined radiometric/geometric (segmentation based) classification:
.RS 4n
.IP \(bu 4n
Supervised classification (i.gensigset,
i.smap)
.RE
.IP \(bu 4n
Object\-oriented classification:
.RS 4n
.IP \(bu 4n
Unsupervised classification (segmentation based:
i.segment)
.RE
.RE
Kappa statistic can be calculated to validate the results
(r.kappa). Covariance/correlation matrices can be
calculated with r.covar.
.SS Image fusion
In case of using multispectral data, improvements of the resolution
can be gained by merging the panchromatic channel with color
channels. GRASS provides the HIS (i.rgb.his,
i.his.rgb) and the Brovey and PCA transform
(i.pansharpen) methods.
.SS Radiometric corrections
Atmospheric effects can be removed with i.atcorr.
Correction for topographic/terrain effects is offered in i.topo.corr.
Clouds in LANDSAT data can be identified and removed with
i.landsat.acca.
Calibrated digital numbers of LANDSAT and ASTER imagery may be converted to
top\-of\-atmosphere radiance or reflectance and temperature (i.aster.toar,
i.landsat.toar).
.SS Time series processing
GRASS also offers support for time series processing (r.series). Statistics can be derived from a
set of coregistered input maps such as multitemporal satellite
data. The common univariate statistics and also linear regression can
be calculated.
.SS Evapotranspiration modeling
In GRASS, several types of evapotranspiration (ET) modeling methods
are available:
.RS 4n
.IP \(bu 4n
Reference ET: Hargreaves (i.evapo.mh),
Penman\-Monteith (i.evapo.pm);
.IP \(bu 4n
Potential ET: Priestley\-Taylor (i.evapo.pt);
.IP \(bu 4n
Actual ET: i.evapo.time.
.RE
Evaporative fraction: i.eb.evapfr,
i.eb.hsebal01.
.SS Energy balance
Emissivity can be calculated with i.emissivity.
Several modules support the calculation of the energy balance:
.RS 4n
.IP \(bu 4n
Actual evapotranspiration for diurnal period  (i.eb.eta);
.IP \(bu 4n
Evaporative fraction and root zone soil moisture (i.eb.evapfr);
.IP \(bu 4n
Sensible heat flux iteration (i.eb.hsebal01);
.IP \(bu 4n
Net radiation approximation (i.eb.netrad);
.IP \(bu 4n
Soil heat flux approximation (i.eb.soilheatflux).
.RE
.SS See also
.RS 4n
.IP \(bu 4n
GRASS GIS Wiki page: Image processing
.IP \(bu 4n
The GRASS 4
\fIImage
Processing manual\fR
.IP \(bu 4n
Introduction into raster data processing
.IP \(bu 4n
Introduction into 3D raster data (voxel) processing
.IP \(bu 4n
Introduction into vector data processing
.IP \(bu 4n
Introduction into temporal data processing
.IP \(bu 4n
Database management
.IP \(bu 4n
Projections and spatial transformations
.RE
.SH SOURCE CODE
.PP
Available at: Image processing in GRASS GIS source code (history)
.PP
Main index |
Imagery index |
Topics index |
Keywords index |
Graphical index |
Full index
.PP
xuebingbing's avatar
xuebingbing committed
253
© 2003\-2020
xuebingbing's avatar
xuebingbing committed
254
GRASS Development Team,
xuebingbing's avatar
xuebingbing committed
255
GRASS GIS 7.8.3 Reference Manual