Nothing Special   »   [go: up one dir, main page]

Next Issue
Volume 2, January
Previous Issue
Volume 1, September
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
remotesensing-logo

Journal Browser

Journal Browser

Remote Sens., Volume 1, Issue 4 (December 2009) – 39 articles , Pages 620-1394

  • Issues are regarded as officially published after their release is announced to the table of contents alert mailing list.
  • You may sign up for e-mail alerts to receive table of contents of newly released issues.
  • PDF is the official format for papers published in both, html and pdf forms. To view the papers in pdf format, click on the "PDF Full-text" link, and use the free Adobe Reader to open them.
Order results
Result details
Section
Select all
Export citation of selected articles as:
1170 KiB  
Article
Leaf Area Index (LAI) Estimation of Boreal Forest Using Wide Optics Airborne Winter Photos
by Terhikki Manninen, Lauri Korhonen, Pekka Voipio, Panu Lahtinen and Pauline Stenberg
Remote Sens. 2009, 1(4), 1380-1394; https://doi.org/10.3390/rs1041380 - 22 Dec 2009
Cited by 26 | Viewed by 14092
Abstract
A new simple airborne method based on wide optics camera is developed for leaf area index (LAI) estimation in coniferous forests. The measurements are carried out in winter, when the forest floor is completely snow covered and thus acts as a light background [...] Read more.
A new simple airborne method based on wide optics camera is developed for leaf area index (LAI) estimation in coniferous forests. The measurements are carried out in winter, when the forest floor is completely snow covered and thus acts as a light background for the hemispherical analysis of the images. The photos are taken automatically and stored on a laptop during the flights. The R2 value of the linear regression of the airborne and ground based LAI measurements was 0.89. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The camera system attached to the helicopter support landing gear. The arrow points to the hole, through which the camera was mechanically attached.</p>
Full article ">Figure 2
<p>Use of principal components in removal of shadows from the images (just a small part of entire image is shown).</p>
Full article ">Figure 3
<p>Image thresholding methods validated with LAI-2000 data.</p>
Full article ">Figure 4
<p>Example of photos of one vertical profile taken at Kommattivaara (67.440° N, 26.745° E) in March 13, 2009 using the wide optics airborne camera. (a) Altitude above ground about 25 m, (b) altitude above ground about 210 m. The corresponding thresholded binary images are shown below.</p>
Full article ">Figure 4 Cont.
<p>Example of photos of one vertical profile taken at Kommattivaara (67.440° N, 26.745° E) in March 13, 2009 using the wide optics airborne camera. (a) Altitude above ground about 25 m, (b) altitude above ground about 210 m. The corresponding thresholded binary images are shown below.</p>
Full article ">Figure 5
<p>A detail of image (a) in <a href="#remotesensing-01-01380-f004" class="html-fig">Figure 4</a> demonstrating the resolution of the camera.</p>
Full article ">Figure 6
<p>The airborne LAI profile corresponding to <a href="#remotesensing-01-01380-f004" class="html-fig">Figure 4</a>.</p>
Full article ">Figure 7
<p>The relationship of the canopy height and the average altitude of the two GPS systems corresponding to the images to be used for LAI retrieval comparison with ground measurements. The missing GPS height values of the camera system were interpolated before taking the average of the two GPS systems. The terrain height was subtracted from the GPS height values. The GPS heights of the two systems are shown as error bars of the average value. The outlier (blue marker) not included in the regression corresponds to a plot, which was so close to the road, that the LAI values derived from the higher altitudes were deteriorated by that.</p>
Full article ">Figure 8
<p>Relationship of the ground based and airborne LAI values. The airborne values were measured in March 13, 2009 and corresponding ground measurements were carried out during August 26–September 5, 2008 The points for airborne measurements 2008 and ground measurements 2007, for which the airborne location is closer than 30 m to that of the ground measurements and which don’t have deciduous LAI contribution, are shown as well.</p>
Full article ">Figure 9
<p>Standwise values of the mean projection of unit foliage area (G) as a function of sun zenith angle, calculated from ground based hemispherical photos of 2008.</p>
Full article ">Figure 10
<p>LAI is underestimated in forests with large percentage of deciduous trees.</p>
Full article ">Figure 11
<p>The observed LAI value versus the flight altitude above ground in March 13, 2009. The values are scaled using the regression line of <a href="#remotesensing-01-01380-f008" class="html-fig">Figure 8</a> for all data points of 2009.</p>
Full article ">
1249 KiB  
Article
Using Urban Landscape Trajectories to Develop a Multi-Temporal Land Cover Database to Support Ecological Modeling
by Jeffrey Hepinstall-Cymerman, Stefan Coe and Marina Alberti
Remote Sens. 2009, 1(4), 1353-1379; https://doi.org/10.3390/rs1041353 - 22 Dec 2009
Cited by 12 | Viewed by 14576
Abstract
Urbanization and the resulting changes in land cover have myriad impacts on ecological systems. Monitoring these changes across large spatial extents and long time spans requires synoptic remotely sensed data with an appropriate temporal sequence. We developed a multi-temporal land cover dataset for [...] Read more.
Urbanization and the resulting changes in land cover have myriad impacts on ecological systems. Monitoring these changes across large spatial extents and long time spans requires synoptic remotely sensed data with an appropriate temporal sequence. We developed a multi-temporal land cover dataset for a six-county area surrounding the Seattle, Washington State, USA, metropolitan region. Land cover maps for 1986, 1991, 1995, 1999, and 2002 were developed from Landsat TM images through a combination of spectral unmixing, image segmentation, multi-season imagery, and supervised classification approaches to differentiate an initial nine land cover classes. We then used ancillary GIS layers and temporal information to define trajectories of land cover change through multiple updating and backdating rules and refined our land cover classification for each date into 14 classes. We compared the accuracy of the initial approach with the landscape trajectory modifications and determined that the use of landscape trajectory rules increased our ability to differentiate several classes including bare soil (separated into cleared for development, agriculture, and clearcut forest) and three intensities of urban. Using the temporal dataset, we found that between 1986 and 2002, urban land cover increased from 8 to 18% of our study area, while lowland deciduous and mixed forests decreased from 21 to 14%, and grass and agriculture decreased from 11 to 8%. The intensity of urban land cover increased with 252 km2 in Heavy Urban in 1986 increasing to 629 km2 by 2002. The ecological systems that are present in this region were likely significantly altered by these changes in land cover. Our results suggest that multi-temporal (i.e., multiple years and multiple seasons within years) Landsat data are an economical means to quantify land cover and land cover change across large and highly heterogeneous urbanizing landscapes. Our data, and similar temporal land cover change products, have been used in ecological modeling of past, present, and likely future changes in ecological systems (e.g., avian biodiversity, water quality). Such data are important inputs for ecological modelers, policy makers, and urban planners to manage and plan for future landscape change. Full article
(This article belongs to the Special Issue Ecological Status and Change by Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Six county study area in western Washington, USA showing the 2002 Urban Growth Areas, elevation, water, and county boundaries.</p>
Full article ">Figure 2
<p>Classification steps, landscape trajectory analysis, and ancillary GIS-derived classes used to develop the 14 land cover classes for each year of imagery.</p>
Full article ">Figure 3
<p>Urban land cover from 1986 to 2002 showing when land transitioned to urban classes with respect to the 2002 Urban Growth Boundaries (UGB).</p>
Full article ">
386 KiB  
Article
“Group Inversion Approach” for Detection of Soil Moisture Temporal-Invariant Locations
by Claudia Notarnicola
Remote Sens. 2009, 1(4), 1338-1352; https://doi.org/10.3390/rs1041338 - 21 Dec 2009
Viewed by 10875
Abstract
This paper presents an approach denominated Group Inversion Approach (GIA) which aims at detecting soil moisture temporal invariants, i.e., the stable temporal soil moisture locations, by using mainly remotely sensed data. The soil moisture temporal invariants are those locations where [...] Read more.
This paper presents an approach denominated Group Inversion Approach (GIA) which aims at detecting soil moisture temporal invariants, i.e., the stable temporal soil moisture locations, by using mainly remotely sensed data. The soil moisture temporal invariants are those locations where independently of the absolute value changes, the relative spatial distribution of soil moisture remains almost constant. In this procedure, the soil moisture values estimated from different inversion approaches and sensor configurations are compared among themselves and with the ground data. The procedure has been tested in a watershed of around 7,000 km2 with data collected during the SMEX’02 experiment in Iowa (USA). The results indicate that fields with inversion errors lower than five times the soil moisture variability detected with ground measurements represent well the mean watershed soil moisture values. The GIA technique has been also found in good agreement with the classical technique used to detect the stable soil moisture features, based exclusively on ground measurements. Full article
(This article belongs to the Special Issue Microwave Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>The overall procedure used to detect soil stable characteristics.</p>
Full article ">Figure 2
<p>The error <math display="inline"> <semantics> <mrow> <msubsup> <mi mathvariant="sans-serif">ξ</mi> <mi mathvariant="normal">t</mi> <mn>2</mn> </msubsup> </mrow> </semantics> </math>variability due to the different algorithms and acquisition dates (5 July, 7 July, 8 July and 9 July).</p>
Full article ">Figure 3
<p>The error variability <math display="inline"> <semantics> <mrow> <msubsup> <mi mathvariant="sans-serif">ξ</mi> <mi mathvariant="normal">t</mi> <mn>2</mn> </msubsup> </mrow> </semantics> </math>for each field and for each analyzed date. The errors have been averaged on all inversion algorithms. The error bars are ±1.0 standard deviation with respect to the mean value. The graph reports also temporal trend indicating the increased variability after the rain event of 5 July.</p>
Full article ">Figure 4
<p>The error ξ<sup>2</sup> variability for each field is reported. The errors have been averaged over all inversion algorithms and all four dates. The error bars are ±1.0 standard deviation with respect to the mean value.</p>
Full article ">Figure 5
<p>Comparison between the watershed soil moisture mean values and the field soil moisture mean values for fields that were considered as stable (WC 01, 09, 13) and as non stable (WC05, 12).</p>
Full article ">Figure 6
<p>Comparison between the errors calculated with the GIA method present in this paper and the classical methodology introduced by Vachaud <span class="html-italic">et al.</span>, [<a href="#B4-remotesensing-01-01338" class="html-bibr">4</a>]. The case “No_ground” indicates when the ground measurements are skipped in the GIA approach.</p>
Full article ">
739 KiB  
Article
Direct Georeferencing of Stationary LiDAR
by Ahmed Mohamed and Benjamin Wilkinson
Remote Sens. 2009, 1(4), 1321-1337; https://doi.org/10.3390/rs1041321 - 17 Dec 2009
Cited by 17 | Viewed by 16441
Abstract
Unlike mobile survey systems, stationary survey systems are given very little direct georeferencing attention. Direct Georeferencing is currently being used in several mobile applications, especially in terrestrial and airborne LiDAR systems. Georeferencing of stationary terrestrial LiDAR scanning data, however, is currently performed indirectly [...] Read more.
Unlike mobile survey systems, stationary survey systems are given very little direct georeferencing attention. Direct Georeferencing is currently being used in several mobile applications, especially in terrestrial and airborne LiDAR systems. Georeferencing of stationary terrestrial LiDAR scanning data, however, is currently performed indirectly through using control points in the scanning site. The indirect georeferencing procedure is often troublesome; the availability of control stations within the scanning range is not always possible. Also, field procedure can be laborious and involve extra equipment and target setups. In addition, the conventional method allows for possible human error due to target information bookkeeping. Additionally, the accuracy of this procedure varies according to the quality of the control used. By adding a dual GPS antenna apparatus to the scanner setup, thereby supplanting the use of multiple ground control points scattered throughout the scanning site, we mitigate not only the problems associated with indirect georeferencing but also induce a more efficient set up procedure while maintaining sufficient precision. In this paper, we describe a new method for determining the 3D absolute orientation of LiDAR point cloud using GPS measurements from two antennae firmly mounted on the optical head of a stationary LiDAR system. In this paper, the general case is derived where the orientation angles are not small; this case completes the theory of stationary LiDAR direct georeferencing. Simulation and real world field experimentation of the prototype implementation suggest a precision of about 0.05 degrees (~1 milli-radian) for the three orientation angles. Full article
(This article belongs to the Special Issue LiDAR)
Show Figures

Figure 1

Figure 1
<p>LiDAR survey crew (left) and manual survey crew installing the reflex targets on the slip surface of the February 17, 2006 landslide in Guinsaugon, Southern Leyte, Philippines, photos courtesy Professor Guitierrez and Riegl USA.</p>
Full article ">Figure 2
<p>Experimental Setup of the Dual-Antenna System (DAS): Two GPS antennae fixed on a metal bar mounted on top of the LiDAR rotating sensor head. The Inertial Measuring Unit (IMU) in the middle of the setup is for check purpose only and is not part of the DAS system.</p>
Full article ">Figure 3
<p>Simulation Results of the Single-Antenna Case showing precision of calculated orientation angles <span class="html-italic">vs.</span> number of device stops; the curve practically starts to dampen at about 20 stops, however the achieved accuracy is not close to what is required for stationary LiDAR.</p>
Full article ">Figure 4
<p>Simulation Results of the Dual-Antenna Case showing precision of calculated orientation angles <span class="html-italic">vs.</span> number of device stops; the curve's actual dampening happens at close to 100 stops; the achieved accuracy is twice as high as what is required for stationary LiDAR.</p>
Full article ">Figure 5
<p>Simulation Results of the Dual-Antenna Case showing precision of calculated orientation angles <span class="html-italic">vs.</span> their magnitude; as the orientation angle gets larger than 5 degrees, the assumption of small angle rotation gets violated and small-angle model becomes invalid.</p>
Full article ">Figure 6
<p>Simulation Results of the Dual-Antenna Case showing precision of calculated orientation angles <span class="html-italic">vs.</span> their magnitude over the first 5 degrees; the 1 mrad precision is possible with orientation angles smaller than 3 degrees.</p>
Full article ">Figure 7
<p>Simulation Results of the Dual-Antenna Case showing precision of calculated orientation angles <span class="html-italic">vs.</span> mounting bar length for systematic two device stops; the curve's actual dampening happens at about 1 m bar length which show the practicality of the method; it can also be noticed that the achieved accuracy at the 1 m bar length value is close to what is required for stationary LiDAR but not enough with two device stops.</p>
Full article ">Figure 8
<p>Simulation Results of the Dual-Antenna Case showing precision of calculated orientation angles <span class="html-italic">vs.</span> mounting bar length for systematic ten device stops; the curve's actual dampening happens again at about 1 m bar length; it can be noticed that the achieved accuracy at the 1 m bar length value with ten device stops is better than what is required for stationary LiDAR.</p>
Full article ">Figure 9
<p>Field Test Results of the DAS showing values of calculated orientation angles <span class="html-italic">vs.</span> number of device stops for different symmetrical configurations of the device stops; the chart shows that configuration plays a role in determining the value of the calculated angle.</p>
Full article ">Figure 10
<p>Field Test Results of the DAS showing precision of calculated orientation angles <span class="html-italic">vs.</span> number of device stops for different symmetrical configurations of the device stops; the chart shows the asymptotic nature of the precision curve as the number of device stops approaches 10 stops.</p>
Full article ">Figure 11
<p>Field Test Results of the DAS showing precision of calculated orientation angle Omega <span class="html-italic">vs.</span> number of device stops for different symmetrical configurations of the device stops; the chart shows the asymptotic nature of the precision curve as the number of device stops approaches 10 stops.</p>
Full article ">Figure 12
<p>GPS Tilt and Azimuth Errors Compared to calculated instantaneous orientation angles using the mounted Inertial Navigation System (INS); one milli-radian accuracy is not achievable from a single GPS DAS observation justifying the need for multiple setups and the need for the algorithm proposed.</p>
Full article ">
827 KiB  
Article
Measurement of Crown Cover and Leaf Area Index Using Digital Cover Photography and Its Application to Remote Sensing
by Burak Pekin and Craig Macfarlane
Remote Sens. 2009, 1(4), 1298-1320; https://doi.org/10.3390/rs1041298 - 15 Dec 2009
Cited by 74 | Viewed by 18263
Abstract
Digital cover photography (DCP) is a high resolution, vertical field-of-view method for ground-based estimation of forest metrics, and has advantages over fisheye sensors owing to its ease of application and high accuracy. We conducted the first thorough technical appraisal of DCP using both [...] Read more.
Digital cover photography (DCP) is a high resolution, vertical field-of-view method for ground-based estimation of forest metrics, and has advantages over fisheye sensors owing to its ease of application and high accuracy. We conducted the first thorough technical appraisal of DCP using both single-lens-reflex (DSLR) and point-and-shoot cameras and concluded that differences result primarily from the better quality optics available for the DSLR camera. File compression, image size and ISO equivalence had little or no effect on estimates of forest metrics. We discuss the application of DCP for ground truthing of remotely sensed canopy metrics, and highlight its strengths over fisheye photography for testing and calibration of vertical field-of-view remote sensing. Full article
Show Figures

Figure 1

Figure 1
<p>Comparison of wide-angle fisheye and 57° photographic images with narrow FOV cover images.</p>
Full article ">Figure 2
<p>Example of a cover image that has been classified into canopy (black), small within crown gaps pixels (white) and large between-crown gaps (grey). Crown cover (0.51) is the fractional cover of black+white pixels. Foliage cover (0.29) is the fractional cover of black pixels. Crown porosity (Equation 1) is the ratio of white to black+white pixels (0.44).</p>
Full article ">Figure 3
<p>Indicative sampling positions for fisheye photography (crosses) and cover images (circles) in the sample plots. Cover images are spaced 10 m apart but sample outside the 40 m × 40 m, effectively sampling an area of 50 m × 50 m, represented by the square box.</p>
Full article ">Figure 4
<p>Comparison of openness, and leaf area index calculated after Lang and Xiang [<a href="#B40-remotesensing-01-01298" class="html-bibr">40</a>], obtained from the Nikon D80 DSLR camera and from the Coolpix 4500 camera using full-frame fisheye photography. Each data point is the mean of 5 images from the 4500 or 10 (medium and small combined) images from the D80. Standard errors (not shown) were between 0.01 and 0.04 for openness and 0.06 and 0.18 for leaf area index. The dashed line is the 1:1 line.</p>
Full article ">Figure 5
<p>Comparison of crown porosity and leaf area index obtained from the Nikon D80 DSLR camera and from the Coolpix 4500 camera using DCP. Each data point is the mean of 25 images from the 4500 or 50 (medium and small) images from the D80. Standard errors (not shown) were between 0.01 and 0.02 for crown porosity and 0.08 and 0.15 for leaf area index. The dashed line is the 1:1 line.</p>
Full article ">Figure 6
<p>The percentage change of crown cover, foliage cover, crown porosity and leaf area index in response to a changed threshold used to segment images into sky and foliage. The Threshold change is the absolute difference between the automatic threshold selected by WinSCANOPY and the manual threshold used for comparison.</p>
Full article ">Figure 7
<p>Comparison of leaf area index calculated after Lang and Xiang [<a href="#B40-remotesensing-01-01298" class="html-bibr">40</a>], and gap fraction, from full-frame fisheye photography with leaf area index and either crown cover or foliage cover obtained from DCP. Open circles represent data plotted against foliage cover and closed circles represent data plotted against crown cover. The dashed line is the 1:1 line. All data are from large JPEGs collected using the Nikon D80 DSLR. Each data point is the mean of 5 full-frame fisheye images or 25 cover images. Standard errors (not shown) were between 0.07 and 0.16 for leaf area index and 0.01 and 0.05 for crown and foliage cover.</p>
Full article ">
452 KiB  
Article
An Empirical Algorithm for Estimating Agricultural and Riparian Evapotranspiration Using MODIS Enhanced Vegetation Index and Ground Measurements of ET. I. Description of Method
by Pamela L. Nagler, Kiyomi Morino, R. Scott Murray, John Osterberg and Edward P. Glenn
Remote Sens. 2009, 1(4), 1273-1297; https://doi.org/10.3390/rs1041273 - 10 Dec 2009
Cited by 67 | Viewed by 13379
Abstract
We used the Enhanced Vegetation Index (EVI) from MODIS to scale evapotranspiration (ETactual) over agricultural and riparian areas along the Lower Colorado River in the southwestern US. Ground measurements of ETactual by alfalfa, saltcedar, cottonwood and arrowweed were expressed as [...] Read more.
We used the Enhanced Vegetation Index (EVI) from MODIS to scale evapotranspiration (ETactual) over agricultural and riparian areas along the Lower Colorado River in the southwestern US. Ground measurements of ETactual by alfalfa, saltcedar, cottonwood and arrowweed were expressed as fraction of potential (reference crop) ETo (EToF) then regressed against EVI scaled between bare soil (0) and full vegetation cover (1.0) (EVI*). EVI* values were calculated based on maximum and minimum EVI values from a large set of riparian values in a previous study. A satisfactory relationship was found between crop and riparian plant EToF and EVI*, with an error or uncertainty of about 20% in the mean estimate (mean ETactual = 6.2 mm d−1, RMSE = 1.2 mm d−1). The equation for ETactual was: ETactual = 1.22 × ETo-BC × EVI*, where ETo-BC is the Blaney Criddle formula for ETo. This single algorithm applies to all the vegetation types in the study, and offers an alternative to ETactual estimates that use crop coefficients set by expert opinion, by using an algorithm based on the actual state of the canopy as determined by time-series satellite images. Full article
(This article belongs to the Special Issue Global Croplands)
Show Figures

Figure 1

Figure 1
<p>Sap flow measurement sites for saltcedar at Cibola National Wildlife Refuge. SL = Slitherin; HS = Hot Springs; SW = Swamp; DT = Diablo Tower; DSW = Diablo Southwest; DE = Diablo East.</p>
Full article ">Figure 2
<p>Climate conditions measured at the Parker, Arizona, AZMET station during periods in which sap flow measurements were made at Cibola National Wildlife Refuge in 2007 (closed circles) and 2008 (open circles). Values are means over the June-August measurement intervals of air temperature (A), vapor pressure deficit (B), and solar radiation (R<sub>S</sub>) (C).</p>
Full article ">Figure 2 Cont.
<p>Climate conditions measured at the Parker, Arizona, AZMET station during periods in which sap flow measurements were made at Cibola National Wildlife Refuge in 2007 (closed circles) and 2008 (open circles). Values are means over the June-August measurement intervals of air temperature (A), vapor pressure deficit (B), and solar radiation (R<sub>S</sub>) (C).</p>
Full article ">Figure 3
<p>Saltcedar plant-specific leaf area index (LAPS) (A), fractional cover (B), leaf-area transpiration (C) and ground-area transpiration (D) at six sites at Cilobla National Wildlife Refuge on the Lower Colorado River.</p>
Full article ">Figure 4
<p>Leaf-area transpiration (E<sub>L</sub>) (closed circles) and stomatal conductance (G<sub>S</sub>) (open circles) of saltcedar at six sites at Cibola National Wildlife Refuge on the Lower Colorado River. Results are hourly mean values over each measurement period. Error bars are standard errors. Horizontal white rectangles above the x-axis denotes daylight hours.</p>
Full article ">Figure 5
<p>Evaporative fraction (EF), defined as ET<sub>actual</sub>/R<sub>n</sub>, for saltcedar at six sites at Cibola National Wildlife Refuge on the Lower Colorado River. Sites are Slitherin (closed circles), Swamp (open circles), Diablo Tower (closed triangles), Diablo Southwest (open triangles), Hot Springs (closed squares) and Diablo East (open squares).</p>
Full article ">Figure 6
<p>Alfalfa ET<sub>actual</sub> at a field in the Palo Verde Irrigation District calculated from the difference in soil moisture levels measured 48 hours after an irrigation (closed circles) and 6–7 days latter (open circles). Each data point is the mean of five probe measurements. SE is the standard error of ET<sub>actual</sub> over all five ports. P is the probability that mean moisture contents across soil depths are equal at the two measurement intervals by paired t-test.</p>
Full article ">Figure 6 Cont.
<p>Alfalfa ET<sub>actual</sub> at a field in the Palo Verde Irrigation District calculated from the difference in soil moisture levels measured 48 hours after an irrigation (closed circles) and 6–7 days latter (open circles). Each data point is the mean of five probe measurements. SE is the standard error of ET<sub>actual</sub> over all five ports. P is the probability that mean moisture contents across soil depths are equal at the two measurement intervals by paired t-test.</p>
Full article ">Figure 7
<p>Ratio of ET<sub>actual</sub> or E<sub>G</sub> to reference crop ET (ET<sub>o</sub>F) for plants on the Lower Colorado River, using the Blaney Criddle method (A) and the Penman Monteith method (B) for ET<sub>o</sub>. Plants and locations are: saltcedar at Slitherin (SL), Swamp (SW), Diablo East (DE), Diablo Southwest (DSW) and Diablo Tower (DT) at Cibola National Wildlife Refuge (closed circles); saltcedar at Hot Springs at Cibola National Wildlife Refuge (cross); saltcedar at Havasu National Wildlife Refuge in 2002 and 2003 (open triangles); arrowweed at Havasu National Wildlife Refuge in 2003 (closed square); and alfalfa at Palo Verde Irrigation District on three dates (open circles). Hot Springs was not included in the regression analyses. Regression equations were passed through the origin and dashed lines denote 95% confidence intervals. Error bars are standard errors of means.</p>
Full article ">Figure 8
<p>ET<sub>o</sub> calculated by the Penman Monteith method (ET<sub>o-PM</sub>), the Blaney Criddle method (ET<sub>o-BC</sub>) and saltcedar ET<sub>actual</sub> measured by Bowen ratio moisture flux towers at Havasu National Wildlife Refuge on the Lower Colorado River in 2002 and 2003. Data are from [<a href="#B7-remotesensing-01-01273" class="html-bibr">7</a>]. Correlation coefficients (r) between each variable are shown under the curves.</p>
Full article ">
593 KiB  
Article
Improving Landsat and IRS Image Classification: Evaluation of Unsupervised and Supervised Classification through Band Ratios and DEM in a Mountainous Landscape in Nepal
by Krishna Bahadur K.C.
Remote Sens. 2009, 1(4), 1257-1272; https://doi.org/10.3390/rs1041257 - 8 Dec 2009
Cited by 61 | Viewed by 18363
Abstract
Modification of the original bands and integration of ancillary data in digital image classification has been shown to improve land use land cover classification accuracy. There are not many studies demonstrating such techniques in the context of the mountains of Nepal. The objective [...] Read more.
Modification of the original bands and integration of ancillary data in digital image classification has been shown to improve land use land cover classification accuracy. There are not many studies demonstrating such techniques in the context of the mountains of Nepal. The objective of this study was to explore and evaluate the use of modified band and ancillary data in Landsat and IRS image classification, and to produce a land use land cover map of the Galaudu watershed of Nepal. Classification of land uses were explored using supervised and unsupervised classification for 12 feature sets containing the LandsatMSS, TM and IRS original bands, ratios, normalized difference vegetation index, principal components and a digital elevation model. Overall, the supervised classification method produced higher accuracy than the unsupervised approach. The result from the combination of bands ration 4/3, 5/4 and 5/7 ranked the highest in terms of accuracy (82.86%), while the combination of bands 2, 3 and 4 ranked the lowest (45.29%). Inclusion of DEM as a component band shows promising results. Full article
(This article belongs to the Special Issue Ecological Status and Change by Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of the study area.</p>
Full article ">Figure 2
<p>Satellite images over Galaudu watershed (a) Landsat MSS 1976 (b) Landsat TM 1990 (c) Landsat TM 2000 (d) IRS LISS III 2002.</p>
Full article ">Figure 2 Cont.
<p>Satellite images over Galaudu watershed (a) Landsat MSS 1976 (b) Landsat TM 1990 (c) Landsat TM 2000 (d) IRS LISS III 2002.</p>
Full article ">Figure 3
<p>Steps of digital image processing to produce land use land cover map.</p>
Full article ">Figure 4
<p>Land use in Galaudu watershed in 1976, 1990, 2000 and 2002.</p>
Full article ">
833 KiB  
Article
Upliftment Estimation of the Zagros Transverse Fault in Iran Using Geoinformatics Technology
by Hojjat Ollah Safari, Saeid Pirasteh and Biswajeet Pradhan
Remote Sens. 2009, 1(4), 1240-1256; https://doi.org/10.3390/rs1041240 - 8 Dec 2009
Cited by 16 | Viewed by 13302
Abstract
The Izeh fault zone is a transverse fault zone with dextral strike slip (and some reverse component) in the Zagros Mountains (Iran). It causes some structural deformations. This fault zone is acting as eastern boundary of Dezful Embayment and forms subsidence of the [...] Read more.
The Izeh fault zone is a transverse fault zone with dextral strike slip (and some reverse component) in the Zagros Mountains (Iran). It causes some structural deformations. This fault zone is acting as eastern boundary of Dezful Embayment and forms subsidence of the embayment. The fault has been recognized using remote sensing techniques in conjunction with surface and subsurface analyses. The stratigraphic columns have been prepared in 3D form using Geographical Information System (GIS) tools on the basis of structural styles and thickness of lithologic units. Height differences for erosion levels have been calculated in stratigraphic columns with respect to the subsidence in the Dezful Embayment, which is related to Izeh zone. These height differences have been estimated to be 5,430 m in the central part (and 5,844 m in the northern part) from the Eocene to recent times. This study shows that comparison of the same erosion levels in Asmari-Pabdeh formation boundaries for interior and eastern block of the Izeh fault zone with the absolute uplifting due to the fault activity which is about 533 m per million years in the Izeh zone. The present study reveals that subtracting the absolute uplifting from total subsidence; the real subsidence of Dezful embayment from Eocene to Recent is 0.13 mm/year. The mean rate of uplifting along the Izeh fault zone is 0.015 mm/year. Full article
(This article belongs to the Special Issue Remote Sensing in Seismology)
Show Figures

Figure 1

Figure 1
<p>Structural zones of the Zagros fold-thrust belt and location of the study area together with a structural map of part of the Izeh fault zone.</p>
Full article ">Figure 2
<p>The study area with structural features and AB cross section on ETM+ image with FCC 7-4-1.</p>
Full article ">Figure 3
<p>Geology map of the study area with vertical displacement along the Izeh Fault Zone and thrust.</p>
Full article ">Figure 4
<p>DEM of the study area.</p>
Full article ">Figure 5
<p>The different stages for 3D sketches generation: (a) The geometric corrected Landsat ETM+ satellite image overlaid on DEM; (b) surveyed cross sections in opposite directions.</p>
Full article ">Figure 6
<p>Cross sections together with subsidence-uplifting.</p>
Full article ">Figure 7
<p>Field photos indicate active faulting and uplifting during the sedimentation and erosion processes. (a) thrust faulting in mountain front; (b) uplifting along fault zone; (c) active faulting in quaternary sediments; (d) bending in mountain front.</p>
Full article ">
2645 KiB  
Review
Antarctic Ice Sheet and Radar Altimetry: A Review
by Frédérique Rémy and Soazig Parouty
Remote Sens. 2009, 1(4), 1212-1239; https://doi.org/10.3390/rs1041212 - 7 Dec 2009
Cited by 80 | Viewed by 16934
Abstract
Altimetry is probably one of the most powerful tools for ice sheet observation. Our vision of the Antarctic ice sheet has been deeply transformed since the launch of the ERS1 satellite in 1991. With the launch of ERS2 and Envisat, the series of [...] Read more.
Altimetry is probably one of the most powerful tools for ice sheet observation. Our vision of the Antarctic ice sheet has been deeply transformed since the launch of the ERS1 satellite in 1991. With the launch of ERS2 and Envisat, the series of altimetric observations now provides 19 years of continuous and homogeneous observations that allow monitoring of the shape and volume of ice sheets. The topography deduced from altimetry is one of the relevant parameters revealing the processes acting on ice sheet. Moreover, altimeter also provides other parameters such as backscatter and waveform shape that give information on the surface roughness or snow pack characteristics. Full article
(This article belongs to the Special Issue Microwave Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Map of the Antarctic ice sheet with the main places cited in the text. Ice shelf names are in red.</p>
Full article ">Figure 2
<p>Ice sheet mechanisms and time-scale outline adapted from [<a href="#B7-remotesensing-01-01212" class="html-bibr">7</a>]. Surface melting, snow precipitation, sublimation, wind-driven sublimation and basal melting or refreezing are assumed to instantaneously react to climate change. Ice streams, outlet glaciers are assumed to react quickly (between 1 yr to 10 yr) or slowly as ice-shelf dynamics to climate change (meaning between 100 and 1,000 yr). On the contrary, ice flow, basal temperature and fusion isostasy take a very long time to react (time lag is between 10,000 and 100,000 years).</p>
Full article ">Figure 3
<p>Altimetric waveform shape. Altimetric observations also provide the return waveform that can be seen as the histogram of the backscattered energy with respect to the return time. The signal is the sum of a surface echo (in light grey) and of a volume echo (in dark grey). The altimeter provides then the surface altitude, the waveform shape (with the parameters such as leading edge width and trailing edge slope) and the total backscattered energy from the surface.</p>
Full article ">Figure 4
<p>Mean altimetric waveform parameters mapped over Antarctica in Ku-band from 2003 to 2007. (a) The backscattering coefficient expressed in dB, (b) the leading edge width expressed in meters, (c) the trailing edge slope expressed in 10<sup>6</sup> s<sup>−1</sup>.</p>
Full article ">Figure 5
<p>Difference in the altimetric waveform parameters between the Ku-band and the S-band, from the Envisat altimeter. (a) The backscattering coefficient expressed in dB, (b) the leading edge width expressed in meters, (c) the trailing edge slope expressed in 10<sup>6</sup> s<sup>−1</sup>.</p>
Full article ">Figure 6
<p>Difference in altimetric height between the Ku-band and the S-band of the altimeter from [<a href="#B22-remotesensing-01-01212" class="html-bibr">22</a>]. This map illustrates very well the effect of the penetration of the radar wave within the snow pack because it depends on the wave frequency. This map can also be used to derive the principal characteristics of the snow pack and of the surface roughness.</p>
Full article ">Figure 7
<p>High-resolution map of the Antarctic Ice Sheet topography from the ERS-1 geodetic mission [<a href="#B7-remotesensing-01-01212" class="html-bibr">7</a>]. The elevation reaches 4000 m. Note the numerous details: ice shelves surrounding most of the continent, surface undulations, flat areas reflecting subglacial lakes, elongated scars due to hydrological networks.</p>
Full article ">Figure 8
<p>Map of the balance velocity expressed in m/yr from [<a href="#B15-remotesensing-01-01212" class="html-bibr">15</a>]. Note the presence of numerous ice streams that can be followed far away in the upslope direction. The name and location of the main outlet glaciers are shown.</p>
Full article ">Figure 9
<p>Enlargement of the changes in height in the west Antarctic ice sheet during the period 2002–2008. The Pine Island, Thwaites and Smith Glaciers sectors (noted A, B, C) decrease at a rate greater than −0.15 m/yr while other glaciers (D or E) increase during this period. D is found to thinning only during the Envisat period, this shows that the upslope signal is propagated in few years over several hundred kilometres. E corresponds to the upslope part of the Kamb ice stream that flows out of the Envisat coverage. Its thinning up to the coast is confirmed with the IceSat data [<a href="#B114-remotesensing-01-01212" class="html-bibr">114</a>].</p>
Full article ">Figure 10
<p>Mass balance of the Antarctic ice sheet (in m/yr) updated from [<a href="#B18-remotesensing-01-01212" class="html-bibr">18</a>], for the ERS period (1992-2003) and the Envisat period (2002-2006). Note the important thinning of the West Antarctic ice sheet over a large sector. The Eastern part of the continent exhibits less impressive signals, but is also shown local fluctuations that depend on the period.</p>
Full article ">Figure 11
Full article ">
872 KiB  
Article
HF Radar Bistatic Measurement of Surface Current Velocities: Drifter Comparisons and Radar Consistency Checks
by Belinda Lipa, Chad Whelan, Bill Rector and Bruce Nyden
Remote Sens. 2009, 1(4), 1190-1211; https://doi.org/10.3390/rs1041190 - 1 Dec 2009
Cited by 28 | Viewed by 12872
Abstract
We describe the operation of a bistatic HF radar network and outline analysis methods for the derivation of the elliptical velocity components from the radar echo spectra. Bistatic operation is illustrated by application to a bistatic pair: Both remote systems receive backscattered echo, [...] Read more.
We describe the operation of a bistatic HF radar network and outline analysis methods for the derivation of the elliptical velocity components from the radar echo spectra. Bistatic operation is illustrated by application to a bistatic pair: Both remote systems receive backscattered echo, with one remote system in addition receiving bistatic echoes transmitted by the other. The pair produces elliptical velocity components in addition to two sets of radials. Results are compared with drifter measurements and checks performed on internal consistency in the radar results. We show that differences in drifter/radar current velocities are consistent with calculated radar data uncertainties. Elliptical and radial velocity components are demonstrated to be consistent within the data uncertainties. Inclusion of bistatic operation in radar networks can be expected to increase accuracy in derived current velocities and extend the coverage area. Full article
(This article belongs to the Special Issue Ocean Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Schematic diagram of a bistatic radar configuration. Radar at transmit site is a transmitter only. Radar at the receive site receives bistatic echo (green) from an elliptical range contour. It also transmits and receives backscatter echo (purple) from a circular range ring.</p>
Full article ">Figure 2
<p>Received sea echo power vs. range at Doppler frequency 0.375 Hz measured by the an HF radar receiver located at Commonweal, California on April 04, 2008, 22:40, showing the offset between received backscatter echo (range &lt; 90 km) and bistatic echo (160 km &lt; range &lt; 230 km) Red: Loop 1 Yellow: Loop 2 Blue: monopole. The channel powers have been artificially separated by 20 dB.</p>
Full article ">Figure 3
<p>Schematic diagrams of radar configuration (a) Backscatter from a circle of radius R around the transmitter/receiver location, showing azimuth angle ⎕⎕⎕of the received signal; (b) Bistatic scatter from an ellipse with receiver and transmitter at the focal points separated by distance <span class="html-italic">F</span>. Ellipse major and minor semi-axes <span class="html-italic">a</span>, <span class="html-italic">b</span>, are shown, as well as the azimuth angle <span class="html-italic">φ</span>⎕⎕⎕ of the received signal, the angles of incidence and reflection⎕⎕, the distances <span class="html-italic">P</span>, <span class="html-italic">Q</span> from the receiver, transmitter to a scatter point on the ellipse, and the slope⎕⎕ of the normal to the ellipse at the scatter point.</p>
Full article ">Figure 4
<p>Radar echo spectra measured by an HF radar receiver located at Commonweal, California on April 02, 2008 23:10 (a) backscatter spectra, range cell 2; (b) bistatic echo spectra from the transmitter at Fort Funston, range cell 59.</p>
Full article ">Figure 5
<p>Radar current velocity maps produced by an HF radar located at Commonweal, California, on April 02, 2008, 23:10 from (a) backscatter; (b) bistatic echo from the transmitter at Fort Funston.</p>
Full article ">Figure 6
<p>Drifter tracks recorded from April 1 to 4, 2008 off the Northern California coast, with positions plotted as triangles, with a different color for each drifter.</p>
Full article ">Figure 7
<p>Drifter tracks recorded from April 1, 2008 14:20 to 18:40 showing the currents veering from an easterly to a westerly direction.</p>
Full article ">Figure 8
<p>Radar total current velocity map produced from the radial velocities measured by the HF radars at COMM, FORT on April 02, 2008 23:10.</p>
Full article ">Figure 9
<p>A comparison map of radial velocities measured by COMM (Black) and drifters (Red) on April 02, 2008 23:10. Time series of radial velocities at the point shown marked <b>X</b> are plotted in <a href="#remotesensing-01-01190-f010" class="html-fig">Figure 10</a>.</p>
Full article ">Figure 10
<p>Time series of radial velocities measured at the point shown in <a href="#remotesensing-01-01190-f009" class="html-fig">Figure 9</a> by COMM (Black) and Drifter (Red).</p>
Full article ">Figure 11
<p>(a) Scatter plot of radial velocity data points: Drifter <span class="html-italic">vs</span>. radar; (b) Comparison statistics averaged over the map, plotted vs. time: Blue—standard deviation between radar and drifter radial velocities. Green—bias. Red—spatial standard deviation in COMM radials Black—temporal standard deviation in COMM radials</p>
Full article ">Figure 12
<p>A comparison map of elliptical velocities, measured by COMM (Black) and drifters (Red) on April 02, 2008, 23:10. Time series of elliptical velocities at the point shown marked <b>X</b> are plotted in <a href="#remotesensing-01-01190-f013" class="html-fig">Figure 13</a>. Normalized radar antenna patterns are shown: Loop 1 (black) ; Loop 2 (red).</p>
Full article ">Figure 13
<p>Time series of elliptical velocities measured at the point shown in <a href="#remotesensing-01-01190-f012" class="html-fig">Figure 12</a> by COMM (Black) and Drifter (Red).</p>
Full article ">Figure 14
<p>(a) Scatter plot of elliptical velocity data points: Drifter vs radar; (b) Comparison statistics averaged over the map, plotted vs. time: Blue—standard deviation between radar and drifter elliptical velocities. Green—bias. Red—spatial standard deviation in COMM ellipticals. Black—temporal standard deviation in COMM ellipticals.</p>
Full article ">Figure 15
<p>A comparison map of radial velocities. Black—calculated from total velocity vectors produced by COMM, FORT; Red—measured by DRAK. Time series of radial velocities at the point shown marked <b>X</b> are plotted in <a href="#remotesensing-01-01190-f016" class="html-fig">Figure 16</a>.</p>
Full article ">Figure 16
<p>Time series of radial velocities measured at the point shown in <a href="#remotesensing-01-01190-f015" class="html-fig">Figure 15</a>. Black—calculated from total velocity vectors produced by COMM, FORT. Red—measured by DRAK.</p>
Full article ">Figure 17
<p>(a) Scatter plot of radial velocity data points: Radial velocity from COMM/FORT totals vs DRAK radial velocity; (b) Comparison statistics averaged over the map, plotted vs. time: Blue –standard deviation between DRAK radials and radials from COMM/FORT total velocities. Green—bias. Red—spatial standard deviation in DRAK radials Black—temporal standard deviation in DRAK radials.</p>
Full article ">Figure 18
<p>Elliptical velocity comparison map. Black—elliptical component of total velocity vectors from COMM, FORT. Red—measured by COMM. Time series of elliptical velocities at the point shown marked <b>X</b> are plotted in <a href="#remotesensing-01-01190-f019" class="html-fig">Figure 19</a>.</p>
Full article ">Figure 19
<p>Time series of elliptical velocities measured at the point shown in <a href="#remotesensing-01-01190-f018" class="html-fig">Figure 18</a>. Black—calculated from total velocity vectors produced by COMM, FORT. Red—measured by COMM.</p>
Full article ">Figure 20
<p>(a) Scatter plot of elliptical velocity data points: Elliptical velocity from COMM/FORT totals <span class="html-italic">vs.</span> COMM elliptical velocity; (b) Comparison statistics averaged over the map, plotted <span class="html-italic">vs.</span> time: Blue—standard deviation between COMM ellipticals and ellipticals from COMM/FORT total velocities. Green—bias. Red—spatial standard deviation in COMM elliptical.</p>
Full article ">
857 KiB  
Article
A Class-Oriented Strategy for Features Extraction from Multidate ASTER Imagery
by Nicola Crocetto and Eufemia Tarantino
Remote Sens. 2009, 1(4), 1171-1189; https://doi.org/10.3390/rs1041171 - 30 Nov 2009
Cited by 28 | Viewed by 11726
Abstract
In this paper we propose a hybrid classification method, adopting the best features extraction strategy for each land cover class on multidate ASTER data. To enable an effective comparison among images, Multivariate Alteration Detection (MAD) transformation was applied in the pre-processing phase, because [...] Read more.
In this paper we propose a hybrid classification method, adopting the best features extraction strategy for each land cover class on multidate ASTER data. To enable an effective comparison among images, Multivariate Alteration Detection (MAD) transformation was applied in the pre-processing phase, because of its high level of automation and reliability in the enhancement of change information among different images. Consequently, different features identification procedures, both spectral and object-based, were implemented to overcome problems of misclassification among classes with similar spectral response. Lastly, a post-classification comparison was performed on multidate ASTER-derived land cover (LC) maps to evaluate the effects of change in the study area. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area and original subsets of ASTER data: acquisition dates of 2003-06-24 (left) and 2007-07-14 (right).</p>
Full article ">Figure 2
<p>NEM processed ASTER data (left), the on-demand Level-2 standard product AST-08 (centre) of (a) 24 June 2003 and (b) 14 July 2007 and the scatter plots showing their correlation (right).</p>
Full article ">Figure 3
<p>Spatial distribution of Land Surface Temperature (TS) in the subsets ASTER 2003 and 2007 data.</p>
Full article ">Figure 4
<p>Workflow of the procedures in the class-oriented approach to features extraction.</p>
Full article ">Figure 5
<p>LC classes after Class-oriented classification on the ASTER data subsets of 24 June, 2003 (left) and 14 July, 2007 (right).</p>
Full article ">Figure 6
<p>(a) Image difference resulting from post-classification comparison of two independently classified images (Black-Change, White-No Change); (b) Subset area of land cover transformations identified on ASTER 2003 and 2007 data.</p>
Full article ">
1422 KiB  
Article
Enhanced Automated Canopy Characterization from Hyperspectral Data by a Novel Two Step Radiative Transfer Model Inversion Approach
by Wouter Dorigo, Rudolf Richter, Frédéric Baret, Richard Bamler and Wolfgang Wagner
Remote Sens. 2009, 1(4), 1139-1170; https://doi.org/10.3390/rs1041139 - 27 Nov 2009
Cited by 57 | Viewed by 17842
Abstract
Automated, image based methods for the retrieval of vegetation biophysical and biochemical variables are often hampered by the lack of a priori knowledge about land cover and phenology, which makes the retrieval a highly underdetermined problem. This study addresses this problem by presenting [...] Read more.
Automated, image based methods for the retrieval of vegetation biophysical and biochemical variables are often hampered by the lack of a priori knowledge about land cover and phenology, which makes the retrieval a highly underdetermined problem. This study addresses this problem by presenting a novel approach, called CRASh, for the concurrent retrieval of leaf area index, leaf chlorophyll content, leaf water content and leaf dry matter content from high resolution solar reflective earth observation data. CRASh, which is based on the inversion of the combined PROSPECT+SAILh radiative transfer model (RTM), explores the benefits of combining semi-empirical and physically based approaches. The approach exploits novel ways to address the underdetermined problem in the context of an automated retrieval from mono-temporal high resolution data. To regularize the inverse problem in the variable domain, RTM inversion is coupled with an automated land cover classification. Model inversion is based on a two step lookup table (LUT) approach: First, a range of possible solutions is selected from a previously calculated LUT based on the analogy between measured and simulated reflectance. The final solution is determined from this subset by minimizing the difference between the variables used to simulate the spectra contained in the reduced LUT and a first guess of the solution. This first guess of the variables is derived from predictive semi-empirical relationships between classical vegetation indices and the single variables. Additional spectral regularization is obtained by the use of hyperspectral data. Results show that estimates obtained with CRASh are significantly more accurate than those obtained with a tested conventional RTM inversion and semi-empirical approach. Accuracies obtained in this study are comparable to the results obtained by various authors for better constrained inversions that assume more a priori information. The completely automated and image-based nature of the approach facilitates its use in operational chains for upcoming high resolution airborne and spaceborne imaging spectrometers. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Overview of CRASh inversion scheme. The upper part shows the general sequence, the grey box shows the inversion approach itself.</p>
Full article ">Figure 2
<p>Template vegetation spectra used in SPECL and typical vegetation types falling into the respective class (in brackets).</p>
Full article ">Figure 3
<p>Location of study site within the catchment of Lake Waging-Taching in Southeast Germany. Field boundaries are superimposed on a false colour HyMap hyperspectral image recorded at June 30, 2003 (Red = band 29 (849nm); green = band 15 (646nm); blue = band 9 (555nm).</p>
Full article ">Figure 4
<p>Inversion results of CRASh approach versus conventional spectral RMSE based approach.</p>
Full article ">Figure 5
<p><span class="html-italic">In situ</span> measurements of <span class="html-italic">Cw</span>, <span class="html-italic">Cdm</span> and <span class="html-italic">LAI</span> versus estimates obtained from spectral reflectance using the CRASh approach (above) and the spectral RMSE (below). Circles refer to the measurements taken at MEA1, bullets to the measurements taken at MEA2. The dotted line indicates the 1:1 line.</p>
Full article ">
470 KiB  
Article
An Empirical Algorithm for Estimating Agricultural and Riparian Evapotranspiration Using MODIS Enhanced Vegetation Index and Ground Measurements of ET. II. Application to the Lower Colorado River, U.S.
by R. Scott Murray, Pamela L. Nagler, Kiyomi Morino and Edward P. Glenn
Remote Sens. 2009, 1(4), 1125-1138; https://doi.org/10.3390/rs1041125 - 20 Nov 2009
Cited by 43 | Viewed by 13171
Abstract
Large quantities of water are consumed by irrigated crops and riparian vegetation in western U.S. irrigation districts. Remote sensing methods for estimating evaporative water losses by soil and vegetation (evapotranspiration, ET) over wide river stretches are needed to allocate water for agricultural and [...] Read more.
Large quantities of water are consumed by irrigated crops and riparian vegetation in western U.S. irrigation districts. Remote sensing methods for estimating evaporative water losses by soil and vegetation (evapotranspiration, ET) over wide river stretches are needed to allocate water for agricultural and environmental needs. We used the Enhanced Vegetation Index (EVI) from MODIS sensors on the Terra satellite to scale ET over agricultural and riparian areas along the Lower Colorado River in the southwestern U.S., using a linear regression equation between ET of riparian plants and alfalfa measured on the ground, and meteorological and remote sensing data, with an error or uncertainty of about 20%. The algorithm was applied to irrigation districts and riparian areas from Lake Mead to the U.S./Mexico border. The results for agricultural crops were similar to results produced by crop coefficients developed for the irrigation districts along the river. However, riparian ET was only half as great as crop coefficient estimates set by expert opinion, equal to about 40% of reference crop evapotranspiration. Based on reported acreages in 2007, agricultural crops (146,473 ha) consumed 2.2 × 109 m3 yr−1 of water. All riparian shrubs and trees (47,014 ha) consumed 3.8 × 108 m3 yr−1, of which saltcedar, the dominant riparian shrub (25,044 ha), consumed 1.8 × 108 m3 yr−1, about 1% of the annual flow of the river. This method could supplement existing protocols for estimating ET by providing an estimate based on the actual state of the canopy as determined by frequent-return satellite data. Full article
(This article belongs to the Special Issue Global Croplands)
Show Figures

Figure 1

Figure 1
<p>Wide area riparian and agricultural areas for which ET was estimated on the Lower Colorado River. MID-Riparian = riparian zone north of Mohave Irrigation District (MID); HNWR = Havasu National Wildlife Refuge; BWNWR = Bill Williams National Wildlife Refuge; PVID = Palo Verde Irrigation District; CNWR = Cibola National Wildlife Refuge; INWR = Imperial National Wildlife Refuge; MLWA = Mittry Lake Wildlife Area; YID = Yuma Valley Irrigation District agricultural fields (images from Quickbird).</p>
Full article ">Figure 2
<p>MODIS estimates (closed circles) of alfalfa (A); cottonwood (B); and arrowweed (C) ET plotted against ground estimates (open circles). Alfalfa was measured in 2006 and 2007 [<a href="#B15-remotesensing-01-01125" class="html-bibr">15</a>]; arrowweed was measured in 2003 [<a href="#B18-remotesensing-01-01125" class="html-bibr">18</a>]; and cottonwood was measured in 2005 [<a href="#B33-remotesensing-01-01125" class="html-bibr">33</a>].</p>
Full article ">Figure 3
<p>MODIS estimates (closed circles) of saltcedar ET at six sites at Cibola National Wildlife Refuge on the Lower Colorado River plotted against ground estimates from sap flow sensors (open circles). Havasu data were collected in 2003–2003 [<a href="#B18-remotesensing-01-01125" class="html-bibr">18</a>]; other sites were measured in 2007 and 2008 [<a href="#B14-remotesensing-01-01125" class="html-bibr">14</a>,<a href="#B15-remotesensing-01-01125" class="html-bibr">15</a>].</p>
Full article ">
324 KiB  
Article
Evaluating the Effects of Environmental Changes on the Gross Primary Production of Italian Forests
by Fabio Maselli, Marco Moriondo, Marta Chiesi, Gherardo Chirici, Nicola Puletti, Anna Barbati and Piermaria Corona
Remote Sens. 2009, 1(4), 1108-1124; https://doi.org/10.3390/rs1041108 - 19 Nov 2009
Cited by 10 | Viewed by 14282
Abstract
A ten-year data-set descriptive of Italian forest gross primary production (GPP) has been recently constructed by the application of Modified C-Fix, a parametric model driven by remote sensing and ancillary data. That data-set is currently being used to develop multivariate regression models which [...] Read more.
A ten-year data-set descriptive of Italian forest gross primary production (GPP) has been recently constructed by the application of Modified C-Fix, a parametric model driven by remote sensing and ancillary data. That data-set is currently being used to develop multivariate regression models which link the inter-year GPP variations of five forest types (white fir, beech, chestnut, deciduous and evergreen oaks) to seasonal values of temperature and precipitation. The five models obtained, which explain from 52% to 88% of the inter-year GPP variability, are then applied to predict the effects of expected environmental changes (+2 °C and increased CO2 concentration). The results show a variable response of forest GPP to the simulated climate change, depending on the main ecosystem features. In contrast, the effects of increasing CO2 concentration are always positive and similar to those given by a combination of the two environmental factors. These findings are analyzed with reference to previous studies on the subject, particularly concerning Mediterranean environments. The analysis confirms the plausibility of the scenarios obtained, which can cast light on the important issue of forest carbon pool variations under expected global changes. Full article
(This article belongs to the Special Issue Ecological Status and Change by Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Spatial distribution of the five forest types considered in Italy.</p>
Full article ">Figure 2
<p>Monthly GPP image of August 2003 obtained by the application of Modified C-Fix (see text for details).</p>
Full article ">Figure 3
<p>Thermo-pluviometric diagram descriptive of the present and future climate scenarios for the areas covered by deciduous oaks forests (FT 8), which are the most widespread forest type over the Italian territory.</p>
Full article ">Figure 4
<p>Correlation coefficients found for the five forest types between annual GPP estimated by Modified C-Fix and seasonal temperatures (A) and rainfall (B).</p>
Full article ">Figure 5
<p>Annual GPP predicted by Modified C-Fix for the five Italian forest types in the environmental scenarios considered (present scenario, climate change, increased atmospheric CO<sub>2</sub> and combination of the two factors).</p>
Full article ">
213 KiB  
Article
A Simple Method to Determine the Timing of Snow Melt by Remote Sensing with Application to the CO2 Balances of Northern Mire and Heath Ecosystems
by Janne Rinne, Mika Aurela and Terhikki Manninen
Remote Sens. 2009, 1(4), 1097-1107; https://doi.org/10.3390/rs1041097 - 19 Nov 2009
Cited by 8 | Viewed by 11413
Abstract
The timing of the disappearance of the snow cover in spring, or snow melt day (SMD), is a key parameter controlling the carbon dioxide balance between the northern mire and heath ecosystems and the atmosphere. We present a simple method for the determination [...] Read more.
The timing of the disappearance of the snow cover in spring, or snow melt day (SMD), is a key parameter controlling the carbon dioxide balance between the northern mire and heath ecosystems and the atmosphere. We present a simple method for the determination of the SMD using a satellite-based surface albedo product (SAL). The method is based on the local change of albedo from higher wintertime values towards the lower summertime values. The satellite SMD timing correlates well with the SMD determined from snow depth measurements at Finnish weather stations (r = 0.86, slope 1.05). In 50% of the cases the error was 3.4 days or less and bias less than half a day. This would lead to a moderate uncertainty in the annual CO2 balance of mire and heath ecosystems, if the published SMD—CO2 balance relations are valid. However, due to the limited data sets available a systematic validation is left for the future. Full article
(This article belongs to the Special Issue Ecological Status and Change by Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Time series of albedo values in Southern Finland in 2006. The coordinates of the selected SAL grid square are 61°48’N, 24°11’E. Siikaneva (61°50’N, 24°12’E) and Hyytiälä (61°51’N, 24°17’E) are surface observation sites within the selected SAL grid square that are located in an open wetland and a Scots pine forest, respectively.</p>
Full article ">Figure 2
<p>Snow melt day in Northern Europe as derived using the SAF surface albedo product in 2005 (Panel A) and 2006 (Panel B).</p>
Full article ">Figure 3
<p><b>(</b>A) Snow melt day (SMD) as derived using the SAL surface albedo product plotted against that derived from surface snow depth measurements. The grey dots are observations from the year 2005 and the black ones from the year 2006. The solid line indicates the 1:1 relation; (B) Cumulative distribution of the absolute differences between SAL and snow-depth-derived SMD.</p>
Full article ">Figure 4
<p>The difference between snow-depth-derived SMD and SAL-derived SMD at Finnish weather stations in the year 2005 against that in the year 2006.</p>
Full article ">Figure 5
<p>Time series of albedo at the Kaamanen mire measured locally (black solid line) and remotely (circles) in 2006. The error bars mark the periods from which the remotely- measured albedo are determined.</p>
Full article ">
12663 KiB  
Article
Remote Sensing of Channels and Riparian Zones with a Narrow-Beam Aquatic-Terrestrial LIDAR
by Jim McKean, Dave Nagel, Daniele Tonina, Philip Bailey, Charles Wayne Wright, Carolyn Bohn and Amar Nayegandhi
Remote Sens. 2009, 1(4), 1065-1096; https://doi.org/10.3390/rs1041065 - 19 Nov 2009
Cited by 169 | Viewed by 17213
Abstract
The high-resolution Experimental Advanced Airborne Research LIDAR (EAARL) is a new technology for cross-environment surveys of channels and floodplains. EAARL measurements of basic channel geometry, such as wetted cross-sectional area, are within a few percent of those from control field surveys. The largest [...] Read more.
The high-resolution Experimental Advanced Airborne Research LIDAR (EAARL) is a new technology for cross-environment surveys of channels and floodplains. EAARL measurements of basic channel geometry, such as wetted cross-sectional area, are within a few percent of those from control field surveys. The largest channel mapping errors are along stream banks. The LIDAR data adequately support 1D and 2D computational fluid dynamics models and frequency domain analyses by wavelet transforms. Further work is needed to establish the stream monitoring capability of the EAARL and the range of water quality conditions in which this sensor will accurately map river bathymetry. Full article
(This article belongs to the Special Issue LiDAR)
Show Figures

Figure 1

Figure 1
<p>Schematic of EAARL data acquisition [<a href="#B16-remotesensing-01-01065" class="html-bibr">16</a>]. (A) Example transmitted pulses; (B) Received waveform of energy (I) as a function of time (t) reflected from vegetation; (C) Received waveform of energy reflected through water. The upper and lower peaks are reflections from the water surface and the bed, respectively.</p>
Full article ">Figure 2
<p>Typical bathymetric waveforms for three photo-detectors with a strong water surface reflection and weaker bottom reflection. Each detector independently records the time history of energy from each laser pulse. Here the three records are also plotted independently, and thus have small offsets on the vertical axis. Normally the offsets are removed during calibration. Digital counts are inverted in the photodetector and lower counts correspond to higher energy return.</p>
Full article ">Figure 3
<p>CIR digital imagery from Elk Creek, tributary to the upper Middle Fork Salmon River, ID, USA. Pixel resolution is about 15 cm.</p>
Full article ">Figure 4
<p>Airborne LIDAR Processing Software (ALPS) user interface. All inset windows display information from the position of the red dot in the lower middle of the flight track map.</p>
Full article ">Figure 5
<p>Typical bathymetric contours mapped by the EAARL in Bear Valley Creek, ID, USA. Red dots are measurement points made by the EAARL during multiple passes over this reach.</p>
Full article ">Figure 6
<p>Bathymetric LIDAR and ground survey comparison at a meandering channel reach. (A) Example cross sections at locations shown in map view: black = ground survey, green = LIDAR survey; (B) Map view comparison. The contours are from ground survey data. The colored intervals are elevation differences calculated as ground survey grid node elevation minus LIDAR grid node elevation. These colored intervals are georeferenced to the contour map and show that the largest errors tend to occur along the stream banks. The dashed line is the location of channel elevation profile; (C) Channel elevation profiles: black = ground survey points (most accurate), red = profile through ground survey raster, green = profile through LIDAR raster.</p>
Full article ">Figure 7
<p>Comparison of cross-sectional wet areas below an arbitrary water surface elevation in the meandering and straight channel study reaches, calculated from ground survey and EAARL-derived rasters. Dashed line is a 1:1 correspondence.</p>
Full article ">Figure 8
<p>Bathymetric LIDAR and ground survey comparison at a straight channel reach. (A) Example cross sections at locations shown in map view: black = ground survey, green = LIDAR survey; (B) Map view comparison. The contours are from ground survey data. The colored intervals are elevation differences calculated as ground survey grid node elevation minus LIDAR grid node elevation. These colored intervals are georeferenced to the contour map and show that the largest errors tend to occur along the stream banks. The dashed line is the location of channel elevation profile; (C) Channel elevation profiles: black = ground survey points (most accurate), red = profile through ground survey raster, green = profile through LIDAR raster.</p>
Full article ">Figure 8 Cont.
<p>Bathymetric LIDAR and ground survey comparison at a straight channel reach. (A) Example cross sections at locations shown in map view: black = ground survey, green = LIDAR survey; (B) Map view comparison. The contours are from ground survey data. The colored intervals are elevation differences calculated as ground survey grid node elevation minus LIDAR grid node elevation. These colored intervals are georeferenced to the contour map and show that the largest errors tend to occur along the stream banks. The dashed line is the location of channel elevation profile; (C) Channel elevation profiles: black = ground survey points (most accurate), red = profile through ground survey raster, green = profile through LIDAR raster.</p>
Full article ">Figure 9
<p>Effects of bathymetric errors on model predictions of bed mobility in a straight, plane bed channel reach. (A) Percent of flow mesh nodes that switched mobility state when flow was the maximum possible inside the channel (6 m<sup>3</sup>/sec); (B) Percent of flow mesh nodes that switched mobility state at low flow conditions (1 m<sup>3</sup>/sec); (C) Spatial distribution of changes in mobility state at high flow and using the critical shear stress for bed mobility of 41 Pa that is applicable for this field site. Brown nodes had no change in predicted mobility, red nodes switched from immobile to mobile and blue nodes switched from mobile to immobile.</p>
Full article ">Figure 10
<p>(A) Detrended DEM of lower Bear Valley Creek, Idaho, USA. The solid field of high elevations in the downstream half of the image is mountains adjacent to the channel which is flowing in a deep canyon through this area; (B) Detailed view of a portion of the detrended DEM. Elevations are relative to an arbitrary datum set during detrending.</p>
Full article ">Figure 11
<p>Bathymetric LIDAR Toolkit (BLT) graphical user interface. A) Detrended elevation raster with 1 m grid spacing. Flow is left to right. B) The Cross Section Layout tool provides a variety of methods to define cross section locations. C) The Cross Section Explorer displays the channel hydraulic geometry at single or multiple cross sections.</p>
Full article ">Figure 12
<p>Histogram of pool distribution in upper Bear Valley Creek, Middle Fork Salmon River drainage, ID, USA. The volume of pools was calculated in each consecutive 300 m reach of this 10 km-long stream segment. EAARL data which cover the channel and floodplain are inset in a standard 10 m DEM. Red dots are locations of individual Chinook salmon spawning nests in 2003.</p>
Full article ">Figure 13
<p>Histogram of distribution of off-channel habitat that is hydraulically connected to the main channel at bankfull flow conditions. Histogram bars show the area of off-channel habitat connected to each 300 m reach of the 10 km stream segment.</p>
Full article ">Figure 14
<p>Gaussian order 6 wavelet.</p>
Full article ">Figure 15
<p>Gaussian (order 6) transform of a channel profile using a wavelet dilation of 100 m. (A) Detrended elevation profile along the channel thalweg; (B) Spatial distribution of spectral power in the thalweg profile. The transition from a meandering pool-riffle channel with undulating bed topography to a straight channel with plane bed topography is noted at the dashed arrow.</p>
Full article ">Figure 16
<p>Change in a reach of Elk Creek, tributary to Middle Fork Salmon River, Idaho, between EAARL surveys in 2004 and 2007. Red and green ellipses and rectangles refer to areas discussed in text. (A) Bathymetry in 2004; (B) 2004 CIR photograph; (C) Bathymetry in 2007; (D) 2007 CIR photograph; (E) Colored intervals map the change in channel topography over the three year period, calculated as 2004 minus 2007 grid node elevations. Contours are from the 2004 bathymetric data and are included to give topographic context to the elevation difference values.</p>
Full article ">
1755 KiB  
Communication
Investigating the Impacts of Landuse-landcover (LULC) Change in the Pearl River Delta Region on Water Quality in the Pearl River Estuary and Hong Kong’s Coast
by Yuanzhi Zhang, Yufei Wang, Yunpeng Wang and Hongyan Xi
Remote Sens. 2009, 1(4), 1055-1064; https://doi.org/10.3390/rs1041055 - 17 Nov 2009
Cited by 24 | Viewed by 14392
Abstract
Water quality information in the coastal region of Hong Kong and the Pearl River Estuary (PRE) is of great concern to the local community. Due to great landuse-landcover (LULC) changes with rapid industrialization and urbanization in the Pearl River Delta (PRD) region, water [...] Read more.
Water quality information in the coastal region of Hong Kong and the Pearl River Estuary (PRE) is of great concern to the local community. Due to great landuse-landcover (LULC) changes with rapid industrialization and urbanization in the Pearl River Delta (PRD) region, water quality in the PRE has worsened during the last 20 years. Frequent red tide and harmful algal blooms have occurred in the estuary and its adjacent coastal waters since the 1980s and have caused important economic losses, also possibly threatening to the coastal environment, fishery, and public health in Hong Kong. In addition, recent literature shows that water nutrients in Victoria Harbor of Hong Kong have been proven to be strongly influenced by both the Pearl River and sewage effluent in the wet season (May to September), but it is still unclear how the PRE diluted water intrudes into Victoria Harbor. Due to the cloudy and rainy conditions in the wet season in Hong Kong, ASAR images will be used to monitor the PRE river plumes and track the intruding routes of PRE water nutrients. In this paper, we first review LULC change in the PRD and then show our preliminary results to analyze water quality spatial and temporal information from remote observations with different sensors in the coastal region and estuary. The study will also emphasizes on time series of analysis of LULC trends related to annual sediment yields and critical source areas of erosion for the PRD region since the 1980s. Full article
(This article belongs to the Special Issue Ecological Status and Change by Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>The study area in the PRE and coastal region of Hong Kong (adopted from MapPoint).</p>
Full article ">Figure 2
<p>Flowchart of the proposed research.</p>
Full article ">Figure 3
<p>Spatial distribution of SPM and Chl-a in the study area using MERIS data.</p>
Full article ">
871 KiB  
Article
MODIS Hotspot Validation over Thailand
by Veerachai Tanpipat, Kiyoshi Honda and Prayoonyong Nuchaiya
Remote Sens. 2009, 1(4), 1043-1054; https://doi.org/10.3390/rs1041043 - 17 Nov 2009
Cited by 29 | Viewed by 18698
Abstract
To ensure remote sensing MODIS hotspot (also known as active fire products or hotspots) quality and precision in forest fire control and management in Thailand, an increased level of confidence is needed. Accuracy assessment of MODIS hotspots utilizing field survey data validation is [...] Read more.
To ensure remote sensing MODIS hotspot (also known as active fire products or hotspots) quality and precision in forest fire control and management in Thailand, an increased level of confidence is needed. Accuracy assessment of MODIS hotspots utilizing field survey data validation is described. A quantitative evaluation of MODIS hotspot products has been carried out since the 2007 forest fire season. The carefully chosen hotspots were scattered throughout the country and within the protected areas of the National Parks and Wildlife Sanctuaries. Three areas were selected as test sites for validation guidelines. Both ground and aerial field surveys were also conducted in this study by the Forest Fire Control Division, National Park, Wildlife and Plant Conversation Department, Ministry of Natural Resources and Environment, Thailand. High accuracy of 91.84 %, 95.60% and 97.53% for the 2007, 2008 and 2009 fire seasons were observed, resulting in increased confidence in the use of MODIS hotspots for forest fire control and management in Thailand. Full article
Show Figures

Figure 1

Figure 1
<p>Thailand–the three validation test sites of MODIS hotspots in 2007 are circled.</p>
Full article ">Figure 2
<p>187 protected areas (in 130 national parks and 57 wildlife sanctuaries) in Thailand.</p>
Full article ">Figure 3
<p>Field validation pictures.</p>
Full article ">
412 KiB  
Article
Detection of Cypress Canopies in the Florida Panhandle Using Subpixel Analysis and GIS
by Jialing Wang and Paul A. Lang
Remote Sens. 2009, 1(4), 1028-1042; https://doi.org/10.3390/rs1041028 - 17 Nov 2009
Cited by 6 | Viewed by 12770
Abstract
In this study, multitemporal subpixel analysis was used to identify cypress canopies from Landsat 7 ETM+ imagery. One spring and one fall image were selected for each of two sites, an eastern one centered on Tallahassee, FL and a western one centered on [...] Read more.
In this study, multitemporal subpixel analysis was used to identify cypress canopies from Landsat 7 ETM+ imagery. One spring and one fall image were selected for each of two sites, an eastern one centered on Tallahassee, FL and a western one centered on Panama City, FL. Signatures derived from the two eastern images were applied on the two western images that served as the control images for accuracy assessment. Results indicated that multitemporal subpixel analysis greatly improved the classification accuracy and signatures developed from one scene could be used to the subpixel classification of another scene with caution. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The location of the study area. The two boxes with thick black lines indicate the coverage of the two images.</p>
Full article ">Figure 2
<p>Flowchart of the research methodology.</p>
Full article ">Figure 3
<p>Cypress distribution in the Florida Panhandle.</p>
Full article ">Figure 4
<p>Examples of cypress occurrence with three probabilities on public lands.</p>
Full article ">
814 KiB  
Article
Assessing the Impact of Canopy Structure Simplification in Common Multilayer Models on Irradiance Absorption Estimates of Measured and Virtually Created Fagus sylvatica (L.) Stands
by Dimitrios Biliouris, Dimitry Van der Zande, Willem W. Verstraeten, Bart Muys and Pol Coppin
Remote Sens. 2009, 1(4), 1009-1027; https://doi.org/10.3390/rs1041009 - 16 Nov 2009
Cited by 7 | Viewed by 12860
Abstract
Multilayer canopy representations are the most common structural stand representations due to their simplicity. Implementation of recent advances in technology has allowed scientists to simulate geometrically explicit forest canopies. The effect of simplified representations of tree architecture (i.e., multilayer representations) of [...] Read more.
Multilayer canopy representations are the most common structural stand representations due to their simplicity. Implementation of recent advances in technology has allowed scientists to simulate geometrically explicit forest canopies. The effect of simplified representations of tree architecture (i.e., multilayer representations) of four Fagus sylvatica (L.) stands, each with different LAI, on the light absorption estimates was assessed in comparison with explicit 3D geometrical stands. The absorbed photosynthetic radiation at stand level was calculated. Subsequently, each geometrically explicit 3D stand was compared with three multilayer models representing horizontal, uniform, and planophile leaf angle distributions. The 3D stands were created either by in situ measured trees or by modelled trees generated with the AMAP plant growth software. The Physically Based Ray Tracer (PBRT) algorithm was used to simulate the irradiance absorbance of the detailed 3D architecture stands, while for the three multilayer representations, the probability of light interception was simulated by applying the Beer-Lambert’s law. The irradiance inside the canopies was characterized as direct, diffuse and scattered irradiance. The irradiance absorbance of the stands was computed during eight angular sun configurations ranging from 10° (near nadir) up to 80° sun zenith angles. Furthermore, a leaf stratification (the number and angular distribution of leaves per LAI layer inside a canopy) analysis between the 3D stands and the multilayer representations was performed, indicating the amount of irradiance each leaf is absorbing along with the percentage of sunny and shadow leaves inside the canopy. The results reveal that a multilayer representation of a stand, using a multilayer modelling approach, greatly overestimated the absorbed irradiance in an open canopy, while it provided a better approximation in the case of a closed canopy. Moreover, the actual stratification of leaves differed significantly between a multilayer representation and a 3D architecture canopy of the same LAI. The deviations in irradiance absorbance were caused by canopy structure, clumping and positioning of leaves. Although it was found that the use of canopy simplifications for modelling purposes in closed canopies is demonstrated as a valid option, special care should be taken when considering forest stands irradiance simulation for sparse canopies and particularly on higher sun zenith angles where the surrounding trees strongly affect the absorbed irradiance and results can highly deviate from the multilayer assumptions. Full article
Show Figures

Figure 1

Figure 1
<p>Trees used for explicit simulation of the <span class="html-italic">Fagus sylvatica</span> (L.) canopy. Left (A), a 9 year old <span class="html-italic">in situ</span> measured tree; right (B) a 15 year old AMAP simulated tree.</p>
Full article ">Figure 2
<p>Simulated <span class="html-italic">Fagus sylvatica</span> (L.) stand comprised of 9 trees. Four LAI scenarios are presented. Upper stands (9 years old Fagus) are derived from in situ measured trees; (A) with the observed plantation distance of 1.8 m × 2.0 m (LAI = 0.58), and (B) with a changed plantation distance of 1.0 m × 1.0 m (LAI = 1.29). Bottom stands (15 years old Fagus) are derived from AMAP library and have an LAI of 2.76 (C) and 3.65 (D) respectively. The axes of the figures indicate the stand dimensions in meters.</p>
Full article ">Figure 3
<p>PBRT generated APAR per leaf for all 3D simulated stands. The four figures present the nadir view of the stands with LAI 0.58, 1.29, 2.76, and 3.65, respectively. The leaves are separated in percentage classes of absorbed light ranging from yellow (80%–100% absorbed light) to black (0%–5% absorbed light). Sun is at 10° from nadir. The axes of the figures indicate the stand dimensions in meters.</p>
Full article ">Figure 4
<p>The two figures present the side view for the stands of LAI of 1.29, and 2.76 left and right respectively. The leaves are separated in percentage classes of absorbed light ranging from yellow (80%–100% absorbed light) to black (0%–5% absorbed light). Sun is at 10° from nadir. The axes of the figures indicate the stand dimensions in meters.</p>
Full article ">Figure 5
<p>Diurnal sun angular variation of canopy APAR (w m<sup>−2</sup>). The four different graphs show absorbed irradiance from all canopy representations and LAIs.</p>
Full article ">Figure 6
<p>Diurnal variation of Absorbed Photosynthetically Active Radiation irradiance (APAR) (w m<sup>−2</sup>). Data were collected for a single day during June 2004 depicting direct, diffuse, and total irradiance for a partly cloudy sky.</p>
Full article ">Figure 7
<p>The PACL per cumulative LAI (increments in steps of one sixth of total LAI) is presented for the 3D architecture stands. At the same time the absorbed irradiance of the multilayer representations are depicted for the equivalent LAI increments.</p>
Full article ">
1231 KiB  
Article
Estimating Flow Resistance of Wetlands Using SAR Images and Interaction Models
by Mercedes Salvia, Mariano Franco, Francisco Grings, Pablo Perna, Roman Martino, Haydee Karszenbaum and Paolo Ferrazzoli
Remote Sens. 2009, 1(4), 992-1008; https://doi.org/10.3390/rs1040992 - 13 Nov 2009
Cited by 13 | Viewed by 12431
Abstract
The inability to monitor wetland drag coefficients at a regional scale is rooted in the difficulty to determine vegetation structure from remote sensing data. Based on the fact that the backscattering coefficient is sensitive to marsh vegetation structure, this paper presents a methodology [...] Read more.
The inability to monitor wetland drag coefficients at a regional scale is rooted in the difficulty to determine vegetation structure from remote sensing data. Based on the fact that the backscattering coefficient is sensitive to marsh vegetation structure, this paper presents a methodology to estimate the drag coefficient from a combination of SAR images, interaction models and ancillary data. We use as test case a severe fire event occurred in the Paraná River Delta (Argentina) at the beginning of 2008, when 10% of the herbaceous vegetation was burned up. A map of the reduction of the wetland drag coefficient is presented. Full article
(This article belongs to the Special Issue Microwave Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Landcover map of the Paraná River Delta at the beginning of 2008 [<a href="#B9-remotesensing-01-00992" class="html-bibr">9</a>]. Field sites areas are marked in red.</p>
Full article ">Figure 2
<p>Envisat ASAR Wide Swath image of the Paraná River Delta prior to the fire events (2008.02.04, left) and after them (2008.09.01, right). Burned junco marsh areas are shown in red, while selected not burned areas are shown in green.</p>
Full article ">Figure 3
<p>Photographs of one of the test sites, prior to the burning events (2008.02.04, left) and after them (2008.09.01, right). A strong decrease in the emerged biomass was observed. Also, it can be seen that the soil condition is similar in both dates.</p>
Full article ">Figure 4
<p>Box plots of C HH σ<sup>0</sup> values of junco marsh sites in normal condition (green), and for the burned condition (red). The associated mean junco plant density (JPD, plants/m<sup>2</sup>) measured in the field sites is also informed.</p>
Full article ">Figure 5
<p>Water level at Zárate city, used to infer water level inside field sites (see [<a href="#B6-remotesensing-01-00992" class="html-bibr">6</a>]). The dates corresponding to ASAR acquisitions are indicated.</p>
Full article ">Figure 6
<p>Box plots of C HH σ<sup>0</sup> values of junco marsh sites in normal condition (green) and burned condition (red) superimposed with a model simulation of the HH σ<sup>0</sup> values of the junco marsh as a function of JPD. The input model parameters used in this simulation are summarized in <a href="#remotesensing-01-00992-t002" class="html-table">Table 2</a>.</p>
Full article ">Figure 7
<p>Plot of the mean values of <span class="html-italic">a<sub>x</sub></span> (black) and <span class="html-italic">a<sub>y</sub></span> (red) as a function of JPD. A linear fit of mean <span class="html-italic">a<sub>x</sub></span>, <span class="html-italic">a<sub>y</sub></span> as a function of JPD is also presented (R<sup>2</sup>(a<sub>x</sub>) = 0.94, R<sup>2</sup>(a<sub>y</sub>) = 0.92).</p>
Full article ">Figure 8
<p>Map of junco marsh drag coefficient derived from the ASAR image of 2008.09.01 using the methodology described in (4.3). Only burned junco marsh area is shown. The segmentation criterion is summarized in <a href="#remotesensing-01-00992-t003" class="html-table">Table 3</a>. Classes’ colours and names are informed in the figure legend. <span class="html-italic">C<sub>d</sub></span> range values are informed in the figure legend (Coastal, &lt;25, 25–34, 34–43, 43–54, 54–65, 65–76, and &gt;76).</p>
Full article ">
1349 KiB  
Article
An Improved ASTER Index for Remote Sensing of Crop Residue
by Guy Serbin, E. Raymond Hunt, Jr., Craig S. T. Daughtry, Gregory W. McCarty and Paul C. Doraiswamy
Remote Sens. 2009, 1(4), 971-991; https://doi.org/10.3390/rs1040971 - 11 Nov 2009
Cited by 103 | Viewed by 16476
Abstract
Unlike traditional ground-based methodology, remote sensing allows for the rapid estimation of crop residue cover (fR). While the Cellulose Absorption Index (CAI) is ideal for fR estimation, a new index, the Shortwave Infrared Normalized Difference Residue Index (SINDRI), utilizing [...] Read more.
Unlike traditional ground-based methodology, remote sensing allows for the rapid estimation of crop residue cover (fR). While the Cellulose Absorption Index (CAI) is ideal for fR estimation, a new index, the Shortwave Infrared Normalized Difference Residue Index (SINDRI), utilizing ASTER bands 6 and 7, is proposed for future multispectral sensors and would be less costly to implement. SINDRI performed almost as well as CAI and better than other indices at five locations in the USA on multiple dates. A minimal upgrade from one broad band to two narrow bands would provide fR data for carbon cycle modeling and tillage verification. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Photographs of <b>(a)</b> conventional (intensive) and <b>(b)</b> conservation (no-till) tillage from the area surrounding Ames, Iowa, USA.</p>
Full article ">Figure 2
<p>Visible, NIR, and SWIR spectra acquired from a ground-based spectrophotoradiometer of soil, crop residue, and live corn canopy along with ASTER, Landsat TM, and CAI bands. 2031, 2101, and 2211 denote CAI bands. Soil and residue spectra were acquired in the lab, corn canopy outdoors. Relative spectral response functions for ASTER and Landsat TM are courtesy USGS [<a href="#B45-remotesensing-01-00971" class="html-bibr">45</a>]. Abbreviations in legend: Cla is Clarion loam (0.8% SOC) from Ames, IA; CrA, Hm, and Res are Crosier loam (0.5% SOC), Houghton muck (44.9% SOC), and 7 month old corn residue from a field with standing stubble from Fulton, IN, respectively; LC is live corn canopy at silking (R1) [<a href="#B46-remotesensing-01-00971" class="html-bibr">46</a>] stage in Beltsville, MD.</p>
Full article ">Figure 3
<p>Composite coefficient of determination (<span class="html-italic">r<sup>2</sup></span>) values for all data sets mentioned in <a href="#remotesensing-01-00971-t001" class="html-table">Table 1</a> for general normalized difference indices (gNDI<span class="html-italic"><sub>i,j</sub></span>) in comparison with line-point transect crop residue cover (<span class="html-italic">f<sub>R</sub></span>) estimates.</p>
Full article ">Figure 4
<p>Comparison of crop residue cover (<span class="html-italic">f<sub>R</sub></span>) and <span class="html-italic">f<sub>R</sub></span> index values with regression lines for the Beltsville, MD data sets [<a href="#B30-remotesensing-01-00971" class="html-bibr">30</a>]. Index acronyms: CAI–Cellulose Absorption Index; SINDRI–Shortwave Infrared Normalized Difference Residue Index; LCA–Lignin-Cellulose Absorption Index; and NDTI–Normalized Difference Tillage Index. Straight lines denote regression lines. Dotted lines denote tillage class boundaries. Data were acquired with an ASD Inc. Fieldspec Pro FR spectroradiometer (Boulder, CO).</p>
Full article ">Figure 5
<p>Comparison of crop residue cover (<span class="html-italic">f<sub>R</sub></span>) and <span class="html-italic">f<sub>R</sub></span> index values with regression lines for Fulton, IN, 29 May 2006. Index acronyms: CAI – Cellulose Absorption Index; SINDRI–Shortwave Infrared Normalized Difference Residue Index; LCA–Lignin-Cellulose Absorption Index; and NDTI–Normalized Difference Tillage Index. Straight lines denote regression lines. Dotted lines denote tillage class boundaries. Data were acquired from aircraft by SpecTIR LLC (Sparks, NV).</p>
Full article ">Figure 6
<p>Comparison of crop residue cover (<span class="html-italic">f<sub>R</sub></span>) and <span class="html-italic">f<sub>R</sub></span> index values with regression lines for the Ames, IA area. Index acronyms: CAI–Cellulose Absorption Index; SINDRI–Shortwave Infrared Normalized Difference Residue Index; LCA–Lignin-Cellulose Absorption Index; and NDTI–Normalized Difference Tillage Index. Straight lines denote regression lines. Dotted lines denote tillage class boundaries. 22 May 2005 and 19 May 2007 multispectral data were acquired by the ASTER sensor. 27 May 2007 hyperspectral data were acquired from aircraft by SpecTIR LLC (Sparks, NV). CAI is only available for 27 May 2007 imagery.</p>
Full article ">Figure 7
<p>Comparison of crop residue cover (<span class="html-italic">f<sub>R</sub></span>) and <span class="html-italic">f<sub>R</sub></span> index values with regression lines for the Pesotum, IL (PIL) and Centreville, MD (CMD) areas. Index acronyms: CAI–Cellulose Absorption Index; SINDRI–Shortwave Infrared Normalized Difference Residue Index; LCA–Lignin-Cellulose Absorption Index; and NDTI–Normalized Difference Tillage Index. Straight lines denote regression lines. Dotted lines denote tillage class boundaries. Data were acquired from aircraft by SpecTIR LLC (Sparks, NV).</p>
Full article ">Figure 8
<p>Range of ASTER derived NDVI values for each data set. BMD, CMD, AIA, FIN, and PIL denote data from the Beltsville, MD, Centreville, MD, Ames, IA, Fulton, IN, and Pesotum, IL sites, respectively. Black dots denote outlier values above the 10th and below and the 90th percentile values (whiskers), the box denotes values between the 25th and 75th percentiles, with the median values being the bars in the middle of the boxes.</p>
Full article ">Figure 9
<p>Comparison of ASTER and Landsat TM NDTI values for the Ames, IA area on 27 May 2007.</p>
Full article ">Figure 10
<p>Comparison of crop residue indices and classification results for Fulton, IN, 29 May 2006. The false-color NIR image utilized narrow spectral bands (red–860.7 nm, green–648.2 nm, and blue–548.2 nm) and is shown to help the reader differentiate between green vegetation pixels and non-vegetated fields that were undergoing tillage and planting operations. CAI is calculated from hyperspectral data, and the other indices via hyperspectral data that were convolved to equivalent ASTER (SINDRI and LCA) and Landsat TM (NDTI) bands. Classification results were post processed using a 5 × 5 majority analysis to minimize noise. Circles denote ground-truth locations. Color-coded circles in the false-color NIR image denote ground truth tillage classes. Data were acquired from aircraft by SpecTIR LLC (Sparks, NV).</p>
Full article ">Figure 11
<p>Comparison of crop residue indices and classification results for the area southwest of Ames, IA. The top three images were calculated from 19 May 2007 multispectral ASTER imagery; the bottom five are derived from airborne hyperspectral imagery acquired by SpecTIR LLC (Sparks, NV) on 27 May 2007. 27 May SINDRI, LCA, and NDTI were calculated from convolved equivalent ASTER and Landsat TM bands. Circles denote ground-truth locations. Color-coded circles in the 27 May 2007 false-color NIR image denote ground truth tillage classes.</p>
Full article ">
1218 KiB  
Article
Analysis of Land Use/Cover Changes and Animal Population Dynamics in a Wildlife Sanctuary in East Africa
by Charles Ndegwa Mundia and Yuji Murayama
Remote Sens. 2009, 1(4), 952-970; https://doi.org/10.3390/rs1040952 - 11 Nov 2009
Cited by 33 | Viewed by 19408
Abstract
Changes in wildlife conservation areas have serious implications for ecological systems and the distribution of wildlife species. Using the Masai Mara ecosystem as an example, we analyzed long-term land use/cover changes and wildlife population dynamics. Multitemporal satellite images, together with physical and social [...] Read more.
Changes in wildlife conservation areas have serious implications for ecological systems and the distribution of wildlife species. Using the Masai Mara ecosystem as an example, we analyzed long-term land use/cover changes and wildlife population dynamics. Multitemporal satellite images, together with physical and social economic data were employed in a post classification analysis with GIS to analyze outcomes of different land use practices and policies. The results show rapid land use/cover conversions and a drastic decline for a wide range of wildlife species. Integration of land use/cover monitoring data and wildlife resources data can allow for the analysis of changes, and can be used to project trends to provide knowledge about potential land use/cover change scenarios and ecological impacts. Full article
(This article belongs to the Special Issue Ecological Status and Change by Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Map of the study area showing the Masai Mara National Reserve and the surrounding privately owned group ranches.</p>
Full article ">Figure 2
<p>Wildlife movements in the study area (modified from Maddock, [<a href="#B13-remotesensing-01-00952" class="html-bibr">13</a>]).</p>
Full article ">Figure 3
<p>Study approach adopted for the analysis of land use/cover changes in Masai Mara Ecosystem.</p>
Full article ">Figure 4
<p>Land use/cover maps of Masai Mara Ecosystem derived from satellite data for 1975, 1986 and 2007.</p>
Full article ">Figure 5
<p>High concentrations of different wildlife species are common in the Masai Mara Ecosystem.</p>
Full article ">Figure 6
<p>Wildlife and livestock population trends in Masai Mara Ecosystem, 1975–2007. Source: Aerial survey by Department of Resource Surveys and Remote Sensing (DRSRS).</p>
Full article ">Figure 7
<p>Livestock and wildlife grazing together. There is increased competition for pastures due to increasing livestock production in Masai Mara.</p>
Full article ">Figure 8
<p>Maps showing the extent of agricultural expansion and mushrooming tourism facilities in Masai Mara Ecosystem. Tourist facilities have increased from five in 1975 to 140 in 2007: (a) 1975 scenario; (b) 2007 scenario.</p>
Full article ">Figure 8 Cont.
<p>Maps showing the extent of agricultural expansion and mushrooming tourism facilities in Masai Mara Ecosystem. Tourist facilities have increased from five in 1975 to 140 in 2007: (a) 1975 scenario; (b) 2007 scenario.</p>
Full article ">Figure 9
<p>Off-road driving in the Masai Mara National Reserve as tourist vehicles track wild animals. The resulting road tracks, which eventually lead to habitat degradation, were digitized from year 2000 aerial photographs.</p>
Full article ">Figure 10
<p>Effects of off road driving in the Masai Mara National Reserve and surrounding areas.</p>
Full article ">Figure 11
<p>Wet and dry rainfall variation in Masai Mara between 1975 and 2007. Source: Kenya Meteorological Department.</p>
Full article ">Figure 12
<p>Conceptual model depicting the factors contributing to habitat loss and wildlife decline in the Masai Mara Ecosystem.</p>
Full article ">
131 KiB  
Review
LiDAR Utility for Natural Resource Managers
by Andrew Thomas Hudak, Jeffrey Scott Evans and Alistair Matthew Stuart Smith
Remote Sens. 2009, 1(4), 934-951; https://doi.org/10.3390/rs1040934 - 11 Nov 2009
Cited by 141 | Viewed by 18819
Abstract
Applications of LiDAR remote sensing are exploding, while moving from the research to the operational realm. Increasingly, natural resource managers are recognizing the tremendous utility of LiDAR-derived information to make improved decisions. This review provides a cross-section of studies, many recent, that demonstrate [...] Read more.
Applications of LiDAR remote sensing are exploding, while moving from the research to the operational realm. Increasingly, natural resource managers are recognizing the tremendous utility of LiDAR-derived information to make improved decisions. This review provides a cross-section of studies, many recent, that demonstrate the relevance of LiDAR across a suite of terrestrial natural resource disciplines including forestry, fire and fuels, ecology, wildlife, geology, geomorphology, and surface hydrology. We anticipate that interest in and reliance upon LiDAR for natural resource management, both alone and in concert with other remote sensing data, will continue to rapidly expand for the foreseeable future. Full article
(This article belongs to the Special Issue LiDAR)
802 KiB  
Article
Remote Sensing and Spectral Characteristics of Desert Sand from Qatar Peninsula, Arabian/Persian Gulf
by Abdulali Sadiq and Fares Howari
Remote Sens. 2009, 1(4), 915-933; https://doi.org/10.3390/rs1040915 - 11 Nov 2009
Cited by 18 | Viewed by 16721
Abstract
Remote sensing data can provide valuable information about the surface expression of regional geomorphologic and geological features of arid regions. In the present study, several processing techniques were applied to reveal such in the Qatar Peninsula. Those included preprocessing for radiometric and geometric [...] Read more.
Remote sensing data can provide valuable information about the surface expression of regional geomorphologic and geological features of arid regions. In the present study, several processing techniques were applied to reveal such in the Qatar Peninsula. Those included preprocessing for radiometric and geometric correction, various enhancement methods, classification, accuracy assessment, contrast stretching, color composition, and principal component analyses. Those were coupled with field groundtruthing and lab analyses. Field groundtruthing included one hundred and forty measurements of spectral reflectance for various sediment exposures representing main sand types in the four studied parts in Qatar. Lab investigations included grain size analysis, X-ray diffraction and laboratory measurements of spectral reflectance. During the course of this study three sand types have been identified: (i) sabkha-derived salt-rich, quartz sand, and (ii) beach-derived calcareous sand and (iii) aeolian dune quartz. Those areas are spectrally distinct in the VNIR, suggesting that VNIR spectral data can be used to discriminate them. The study found that the main limitation of the ground spectral reflectance study is the difficulty of covering large areas. The study also found that ground and laboratory spectral radiance are generally higher in reflectance than those of Landsat TM. This is due to several factors such as atmospheric conditions, the low altitude or different scales. Whereas for areas with huge size of dune sand, the Landsat TM spectral has higher reflectance than those from field and laboratory. The study observed that there is a good correspondence or correlation of the wavelengths maximum sensitivity between the three spectral measurements i.e lab, field and space-borne measurements. Full article
Show Figures

Figure 1

Figure 1
<p>Locator map of Qatar and surrounding countries.</p>
Full article ">Figure 2
<p>Natural color image of Qatar as well as RGB and PC images of the western part of the study area.</p>
Full article ">Figure 3
<p>RGB and PC images of the north’ tip, northeast coast and southeast coast of Qatar.</p>
Full article ">Figure 4
<p>Unsupervised classification of: (a) west Qatar (1. Inland sabkha, 2. Upper Dammam Fm., 3. Rus Fm., 4. Salty sand, 5. Lower Dammam Fm.,6. Cherty sand and 7. Coastal sabkha ) (b) north east coast (1. Coastal sabkha, 2. Beach sand mixed with sabkha, 3. Blown beach sand, 4. Upper Dammam Fm, 5. Depression sand and Rodah soil and 6. Rus Fm.), (c) north tip (1. Beach-derived sand mixed with sabkha, 2. Uppper Dammam Fm., 3. Blown sand, and 4. Coastal Sabkha) and (d) southeast of Qatar Peninsula (Coastal sabkha cover with water, 2. Upper Dammam Formation, 3. Dry sabkha, 4. Beach-derived sand, 5. Sabkha deposit mixed with dune sand, 6. Sand dune, and 7. Exposed sedimentary outcrop).</p>
Full article ">Figure 5
<p>Spectral profiles of: (a) Landsat TM raw digital number (DN) values for the study area and (b) Landsat TM mean spectral reflectance of various sand types.</p>
Full article ">Figure 6
<p>Example of Landsat TM spectral reflectance from the studied areas.</p>
Full article ">Figure 7
<p>Field and laboratory spectral profiles of the main sand fields in Qatar.</p>
Full article ">Figure 8
<p>Mean spectral reflectances of: (a) Sabkha derived sand of west of Qatar and (b) beach-derived sand, northeast coast of Qatar Peninsula, (c) beach derived sand, north tip, and (d) dune sand, southeast coast of Qatar Peninsula.</p>
Full article ">Figure 9
<p>Average of grain size analyses of the three sand types.</p>
Full article ">Figure 10
<p>XRD analysis for: (a) west sabkha-derived sand, (b) northeast beach-derived sand and (c) southeast sand dunes of Qatar Peninsula (Gyp: Gypsum; Bass: Bassanite; Q: Quartz , P-PF: Plagoclase-Feldespar; Cal: Calcite; Dol: Dolomite; Arag: Aragonite).</p>
Full article ">
970 KiB  
Article
Regional Assessment of Aspen Change and Spatial Variability on Decadal Time Scales
by Temuulen Tsagaan Sankey
Remote Sens. 2009, 1(4), 896-914; https://doi.org/10.3390/rs1040896 - 10 Nov 2009
Cited by 11 | Viewed by 11820
Abstract
Quaking aspen (Populus tremuloides) is commonly believed to be declining throughout western North America. Using a historical vegetation map and Landsat TM5 imagery, this study detects changes in regional aspen cover over two different time periods of 85 and 18 years [...] Read more.
Quaking aspen (Populus tremuloides) is commonly believed to be declining throughout western North America. Using a historical vegetation map and Landsat TM5 imagery, this study detects changes in regional aspen cover over two different time periods of 85 and 18 years and examines aspen change patterns with biophysical variables in the Targhee National Forest of eastern Idaho, USA. A subpixel classification approach was successfully used to classify aspen. The results indicate greater spatial variability in regional aspen change patterns than indicated by local-scale studies. The observed spatial variability appears to be an inherent pattern in regional aspen dynamics, which interacts with biophysical variables, but persists over time. Full article
(This article belongs to the Special Issue Ecological Status and Change by Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Study area and 300 randomly-generated sample polygon locations in the Targhee National Forest in Idaho, USA.</p>
Full article ">Figure 2
<p>Mean spectral reflectance of aspen and other dominant vegetation cover types within the study area in green (G = 0.52–0.60 μm), red (R = 0.63–0.69 μm), near infrared (NIR = 0.76–0.90 μm), and middle infrared (Mid IR = 1.55–1.75 μm) portions of the electromagnetic spectrum in the fall (F) and summer (S) seasons. Error bars are standard errors.</p>
Full article ">Figure 3
<p>Aspen presence and absence classification of 2005 Landsat TM5 multitemporal composite using Mixture Tuned Matched Filtering (MTMF) technique with a regression approach. The exponential regression model was fitted to the MTMF-produced matched filtering scores and infeasibility values (<span class="html-italic">R<sup>2</sup></span> = 0.57, <span class="html-italic">p</span> &lt; 0.0001). Pixels that fell under the regression curve (solid grey line) that had matched filtering scores of 0.5–1 (dashed grey lines) and infeasibility values of &lt;5 (dotted grey line) were classified as aspen presence. All other pixels were classified as aspen absence.</p>
Full article ">Figure 4
<p>The Mixture Tuned Matched Filtering (MTMF) classification images and final aspen map for 2005. The image of matched filtering scores (a) estimates the abundance of target cover within each pixel, while the image of infeasibility values (b) indicates the relative accuracy of the matched filtering score in each pixel. Aspen presence and absence map (c) is produced after the regression integration of the two images.</p>
Full article ">Figure 5
<p>Examples of local-scale aspen changes between 1920 and 2005. Simple image differencing was performed using 1920 (a) and 2005 (b) aspen presence and absence maps, which resulted in three different classes: aspen decrease, no-change, and aspen increase.</p>
Full article ">Figure 6
<p>Aspen change patterns in the 1920–2005 time period. Aspen increase (positive grey bars) and aspen decrease (negative black bars) were simultaneously analyzed as two response variables in a MANOVA model. Grazing, forest harvest, and vegetation cover types were significant predictor variables in aspen increase (p &lt; 0.05), while all predictor variables were significant in aspen decrease (p &lt; 0.05). (a) Aspen changes patterns and grazing; (b) Aspen change patterns and forest harvest; (c) Aspen change patterns and vegetation cover type (LP pine = Lodgepole pine); (d) Aspen change patterns and forest stand age.</p>
Full article ">Figure 7
<p>Aspen change patterns in the 1987–2005 time period. Aspen increase (positive grey bars) and aspen decrease (negative black bars) were simultaneously analyzed as two response variables in a MANOVA model. Grazing and vegetation cover types were statistically significant predictor variables in aspen increase (p &lt; 0.05), while forest harvest and stand age were significant in aspen decrease (p &lt; 0.05). (a) Aspen changes patterns and grazing; (b) Aspen change patterns and forest harvest; (c) Aspen change patterns and vegetation cover type (LP pine = Lodgepole pine); (d) Aspen change patterns and forest stand age.</p>
Full article ">
1570 KiB  
Article
Supervised Classification of Agricultural Land Cover Using a Modified k-NN Technique (MNN) and Landsat Remote Sensing Imagery
by Luis Samaniego and Karsten Schulz
Remote Sens. 2009, 1(4), 875-895; https://doi.org/10.3390/rs1040875 - 9 Nov 2009
Cited by 48 | Viewed by 17028
Abstract
Nearest neighbor techniques are commonly used in remote sensing, pattern recognition and statistics to classify objects into a predefined number of categories based on a given set of predictors. These techniques are especially useful for highly nonlinear relationship between the variables. In most [...] Read more.
Nearest neighbor techniques are commonly used in remote sensing, pattern recognition and statistics to classify objects into a predefined number of categories based on a given set of predictors. These techniques are especially useful for highly nonlinear relationship between the variables. In most studies the distance measure is adopted a priori. In contrast we propose a general procedure to find an adaptive metric that combines a local variance reducing technique and a linear embedding of the observation space into an appropriate Euclidean space. To illustrate the application of this technique, two agricultural land cover classifications using mono-temporal and multi-temporal Landsat scenes are presented. The results of the study, compared with standard approaches used in remote sensing such as maximum likelihood (ML) or k-Nearest Neighbor (k-NN) indicate substantial improvement with regard to the overall accuracy and the cardinality of the calibration data set. Also, using MNN in a soft/fuzzy classification framework demonstrated to be a very useful tool in order to derive critical areas that need some further attention and investment concerning additional calibration data. Full article
Show Figures

Figure 1

Figure 1
<p>Distribution and location of the test sites.</p>
Full article ">Figure 2
<p>Chronogram of the vegetation periods and management activities. The dates of the landsat scenes are also depicted.</p>
Full article ">Figure 3
<p>Scatterplots depicting the location of the land use classes in the predictor space <math display="inline"> <mi mathvariant="bold">x</mi> </math> at three different points in time as indicated at each panel. In this case, the reflectance measurements represent bands 3-4-5 of the training set whose sample size is <math display="inline"> <mrow> <mi>n</mi> <mo>=</mo> <mn>270</mn> </mrow> </math>.</p>
Full article ">Figure 4
<p>Effect of the optimization on the variance function <math display="inline"> <msub> <mi>G</mi> <msub> <mi>B</mi> <mn>1</mn> </msub> </msub> </math>. Here <math display="inline"> <mrow> <mi>h</mi> <mo>=</mo> <mn>1</mn> </mrow> </math>, i.e., the variance function of class 1 for the training set <math display="inline"> <mi mathvariant="script">T</mi> </math>.</p>
Full article ">Figure 5
<p>Scatterplots depicting the location of the land use classes in the embedded space <math display="inline"> <msub> <mi mathvariant="bold">u</mi> <mn>1</mn> </msub> </math>. Class 1 (encircled by a continuous line) is completely isolated from the others. Sample size <math display="inline"> <mrow> <mi>n</mi> <mo>=</mo> <mn>270</mn> </mrow> </math>.</p>
Full article ">Figure 6
<p>Graph showing the variation of the 90% confidence interval depending on the number of nearest neighbors. As reference the confidence intervals for ML and k-NN (<math display="inline"> <mi mathvariant="bold">B</mi> </math> not optimized) are also shown.</p>
Full article ">Figure 7
<p>Graph showing the influence of the sample size per class <math display="inline"> <msub> <mi>n</mi> <mi>h</mi> </msub> </math> on the overall accuracy of the classification. The dimension of the transformation matrix <math display="inline"> <mi mathvariant="bold">B</mi> </math> in the mono- and multi-temporal classifications is <math display="inline"> <mrow> <mo>[</mo> <mn>3</mn> <mo>×</mo> <mn>6</mn> <mo>]</mo> </mrow> </math> and <math display="inline"> <mrow> <mo>[</mo> <mn>3</mn> <mo>×</mo> <mn>9</mn> <mo>]</mo> </mrow> </math> respectively.</p>
Full article ">Figure 8
<p>Graph showing the influence of the type of transformation and the number of nearest neighbors on the overall accuracy of the classification.</p>
Full article ">Figure 9
<p>Ambiguity index [panel (a)] and land cover map obtained with the MNN classifier using <math display="inline"> <mrow> <mi>N</mi> <mo>=</mo> <mn>5</mn> </mrow> </math> neighbors [panel (b)].</p>
Full article ">Figure 10
<p>Histograms depicting the ambiguity index <math display="inline"> <mrow> <mi>b</mi> <mo>(</mo> <msub> <mi mathvariant="bold">y</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </math> as a function of the sample size of the training set <math display="inline"> <msub> <mi>n</mi> <mi>h</mi> </msub> </math>. (All histograms were obtained with the MNN classifier using <math display="inline"> <mrow> <mi>N</mi> <mo>=</mo> <mn>5</mn> </mrow> </math> neighbors.)</p>
Full article ">
555 KiB  
Article
Hyperspectral Reflectance and Fluorescence Imaging to Detect Scab Induced Stress in Apple Leaves
by Stephanie Delalieux, Annemarie Auwerkerken, Willem W. Verstraeten, Ben Somers, Roland Valcke, Stefaan Lhermitte, Johan Keulemans and Pol Coppin
Remote Sens. 2009, 1(4), 858-874; https://doi.org/10.3390/rs1040858 - 6 Nov 2009
Cited by 64 | Viewed by 15420
Abstract
Apple scab causes significant losses in the production of this fruit. A timely and more site-specific monitoring and spraying of the disease could reduce the number of applications of fungicides in the fruit industry. The aim of this leaf-scale study therefore lies in [...] Read more.
Apple scab causes significant losses in the production of this fruit. A timely and more site-specific monitoring and spraying of the disease could reduce the number of applications of fungicides in the fruit industry. The aim of this leaf-scale study therefore lies in the early detection of apple scab infections in a non-invasive and non-destructive way. In order to attain this objective, fluorescence- and hyperspectral imaging techniques were used. An experiment was conducted under controlled environmental conditions, linking hyperspectral reflectance and fluorescence imaging measurements to scab infection symptoms in a susceptible apple cultivar (Malus x domestica Borkh. cv. Braeburn). Plant stress was induced by inoculation of the apple plants with scab spores. The quantum efficiency of Photosystem II (PSII) photochemistry was derived from fluorescence images of leaves under light adapted conditions. Leaves inoculated with scab spores were expected to have lower PSII quantum efficiency than control (mock) leaves. However, besides scab-induced, also immature leaves exhibited low PSII quantum efficiency. Therefore, this study recommends the simultaneous use of fluorescence imaging and hyperspectral techniques. A shortwave infrared narrow-waveband ratio index (R1480/R2135) is presented in this paper as a promising tool to identify scab stress before symptoms become visible to the naked eye. Low PSII quantum efficiency attended by low narrow waveband R1480/R2135 index values points out scab stress in an early stage. Apparent high PSII quantum efficiency together with high overall reflectance in VIS and SWIR spectral domains indicate a severe, well-developed scab infection. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Histograms depicting the distribution of the Φ<sub>PSII</sub> pixel values throughout time for all measured scab inoculated (grey) and control leaves (black). Histograms were constructed for images taken at 3, 4, 5, 6, 7, 10, 14 and 17 days after inoculation.</p>
Full article ">Figure 2
<p>The average of the continuum removed SWIR areas for each measurement day for scab-inoculated (×) and control (mock) (▪) leaves. Ratio of the average areas of both treatments were plotted on the right Y-axis and represented by connected dots.</p>
Full article ">Figure 3
<p>(top) The average spectral behavior of wavelengths 1,480 nm and 2,135 nm for scab inoculated (––) and mock leaves (-) as obtained from the first until the 17<sup>th</sup> day after inoculation. (bottom) The average spectral behavior of ratio 1480/2135 for the scab inoculated (––) and mock leaves (-) as obtained from the first until the 17<sup>th</sup> day after inoculation.</p>
Full article ">Figure 4
<p>Boxplots representing the behavior of the index R<sub>1480</sub>/R<sub>2135</sub> for each measurement day. The measurement days are shown on the X-axis with * representing measurements on scab inoculated leaves and without * the control leaves.</p>
Full article ">Figure 5
<p>(a) c-Values representing the discriminatory performance of all possible two spectral ratio indices throughout the full spectral region (extracted from [<a href="#B29-remotesensing-01-00858" class="html-bibr">29</a>]). Logistic regression was used for binary classification of scab-inoculated and control leaves. (b) Determination coefficients of linear models fitted through protein content and all possible two spectral ratio indices. Data taken from LOPEX.</p>
Full article ">Figure 6
<p>Boxplots depicting the outcomes of the destructive pigment extraction on different measurement days. Measurement days indicated with * represent results of scab inoculated leaves and without * the control leaves.</p>
Full article ">
207 KiB  
Article
Derivation of Soil Line Influence on Two-Band Vegetation Indices and Vegetation Isolines
by Hiroki Yoshioka, Tomoaki Miura, José A. M. Demattê, Karim Batchily and Alfredo R. Huete
Remote Sens. 2009, 1(4), 842-857; https://doi.org/10.3390/rs1040842 - 3 Nov 2009
Cited by 14 | Viewed by 12842
Abstract
This paper introduces derivations of soil line influences on two-band vegetation indices (VIs) and vegetation isolines in the red and near infra-red reflectance space. Soil line variations are described as changes in the soil line parameters (slope and offset) and the red reflectance [...] Read more.
This paper introduces derivations of soil line influences on two-band vegetation indices (VIs) and vegetation isolines in the red and near infra-red reflectance space. Soil line variations are described as changes in the soil line parameters (slope and offset) and the red reflectance of the soil surface. A general form of a VI model equation written as a ratio of two linear functions (e.g., NDVI and SAVI) was assumed. It was found that relative VI variations can be approximated by a linear combination of the three soil parameters. The derived expressions imply the possibility of estimating and correcting for soil-induced bias errors in VIs and their derived biophysical parameters, caused by the assumption of a general soil line, through the use of external data sources such as regional soil maps. Full article
Show Figures

Figure 1

Figure 1
<p>Illustrations of (a) soil line variation and (b) vegetation isoline variation.</p>
Full article ">
563 KiB  
Article
Photogrammetric Methodology for the Production of Geomorphologic Maps: Application to the Veleta Rock Glacier (Sierra Nevada, Granada, Spain)
by Javier De Matías, José Juan De Sanjosé, Gonzalo López-Nicolás, Carlos Sagüés and José Jesús Guerrero
Remote Sens. 2009, 1(4), 829-841; https://doi.org/10.3390/rs1040829 - 28 Oct 2009
Cited by 19 | Viewed by 13242
Abstract
In this paper we present a stereo feature-based method using SIFT (Scale-invariant feature transform) descriptors. We use automatic feature extractors, matching algorithms between images and techniques of robust estimation to produce a DTM (Digital Terrain Model) using convergent shots of a rock glacier.The [...] Read more.
In this paper we present a stereo feature-based method using SIFT (Scale-invariant feature transform) descriptors. We use automatic feature extractors, matching algorithms between images and techniques of robust estimation to produce a DTM (Digital Terrain Model) using convergent shots of a rock glacier.The geomorphologic structure observed in this study is the Veleta rock glacier (Sierra Nevada, Granada, Spain). This rock glacier is of high scientific interest because it is the southernmost active rock glacier in Europe and it has been analyzed every year since 2001. The research on the Veleta rock glacier is devoted to the study of its displacement and cartography through geodetic and photogrammetric techniques. Full article
(This article belongs to the Special Issue Geomorphological Processes and Natural Hazards)
Show Figures

Figure 1

Figure 1
<p>(a). Veleta rock glacier situation; (b). Close view of the glacier from Mount Veleta.</p>
Full article ">Figure 2
<p>(a) Photogrammetric control point; (b) Photogrametric shots from mount Veleta towards the rock glacier.</p>
Full article ">Figure 3
<p>(a) Geometry of normal photogrammetry; (b) Geometry of convergent images. O<sub>1</sub> and O<sub>2</sub> are the centres of projection of the photographs.</p>
Full article ">Figure 4
<p>Map of the Veleta rock glacier (scale: 1/1,000, contour lines: 1 m).</p>
Full article ">Figure 5
<p>Photograph of the distribution of control points (E1–E9) denoted with marks, around the Veleta rock glacier (2008).</p>
Full article ">Figure 6
<p>SIFT matched points denoted with crosses between two images of the Veleta rock glacier area.</p>
Full article ">Figure 7
<p>Camera pose estimation (front view and top view).</p>
Full article ">Figure 8
<p>Example of Delaunay triangulation of the rock glacier (the view is oriented as in <a href="#remotesensing-01-00829-f009" class="html-fig">Figure 9</a>).</p>
Full article ">Figure 9
<p>Comparison of our reconstructed surface cross sections (black polylines) and total station cross sections (blue polylines) and their situation in a Veleta rock glacier sketch.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop