[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Issue
Volume 2, May
Previous Issue
Volume 2, March
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 2, Issue 4 (April 2010) – 14 articles , Pages 908-1196

  • 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:
398 KiB  
Article
Potential of Using Remote Sensing Techniques for Global Assessment of Water Footprint of Crops
by Mireia Romaguera, Arjen Y. Hoekstra, Zhongbo Su, Maarten S. Krol and Mhd. Suhyb Salama
Remote Sens. 2010, 2(4), 1177-1196; https://doi.org/10.3390/rs2041177 - 26 Apr 2010
Cited by 68 | Viewed by 14120
Abstract
Remote sensing has long been a useful tool in global applications, since it provides physically-based, worldwide, and consistent spatial information. This paper discusses the potential of using these techniques in the research field of water management, particularly for ‘Water Footprint’ (WF) studies. The [...] Read more.
Remote sensing has long been a useful tool in global applications, since it provides physically-based, worldwide, and consistent spatial information. This paper discusses the potential of using these techniques in the research field of water management, particularly for ‘Water Footprint’ (WF) studies. The WF of a crop is defined as the volume of water consumed for its production, where green and blue WF stand for rain and irrigation water usage, respectively. In this paper evapotranspiration, precipitation, water storage, runoff and land use are identified as key variables to potentially be estimated by remote sensing and used for WF assessment. A mass water balance is proposed to calculate the volume of irrigation applied, and green and blue WF are obtained from the green and blue evapotranspiration components. The source of remote sensing data is described and a simplified example is included, which uses evapotranspiration estimates from the geostationary satellite Meteosat 9 and precipitation estimates obtained with the Climatic Prediction Center Morphing Technique (CMORPH). The combination of data in this approach brings several limitations with respect to discrepancies in spatial and temporal resolution and data availability, which are discussed in detail. This work provides new tools for global WF assessment and represents an innovative approach to global irrigation mapping, enabling the estimation of green and blue water use. Full article
(This article belongs to the Special Issue Global Croplands)
Show Figures

Figure 1

Figure 1
<p>Inputs for water footprint estimation.</p>
Full article ">Figure 2
<p>Flowchart proposed for obtaining WF of crops from remote sensing data.</p>
Full article ">Figure 3
<p>(a) Monthly green ET and (b) Monthly blue ET obtained in October 2009, using remote sensing data in Egypt.</p>
Full article ">
6930 KiB  
Article
Remote Sensing of Vegetation Structure Using Computer Vision
by Jonathan P. Dandois and Erle C. Ellis
Remote Sens. 2010, 2(4), 1157-1176; https://doi.org/10.3390/rs2041157 - 21 Apr 2010
Cited by 226 | Viewed by 25813
Abstract
High spatial resolution measurements of vegetation structure in three-dimensions (3D) are essential for accurate estimation of vegetation biomass, carbon accounting, forestry, fire hazard evaluation and other land management and scientific applications. Light Detection and Ranging (LiDAR) is the current standard for these measurements [...] Read more.
High spatial resolution measurements of vegetation structure in three-dimensions (3D) are essential for accurate estimation of vegetation biomass, carbon accounting, forestry, fire hazard evaluation and other land management and scientific applications. Light Detection and Ranging (LiDAR) is the current standard for these measurements but requires bulky instruments mounted on commercial aircraft. Here we demonstrate that high spatial resolution 3D measurements of vegetation structure and spectral characteristics can be produced by applying open-source computer vision algorithms to ordinary digital photographs acquired using inexpensive hobbyist aerial platforms. Digital photographs were acquired using a kite aerial platform across two 2.25 ha test sites in Baltimore, MD, USA. An open-source computer vision algorithm generated 3D point cloud datasets with RGB spectral attributes from the photographs and these were geocorrected to a horizontal precision of <1.5 m (root mean square error; RMSE) using ground control points (GCPs) obtained from local orthophotographs and public domain digital terrain models (DTM). Point cloud vertical precisions ranged from 0.6 to 4.3 m RMSE depending on the precision of GCP elevations used for geocorrection. Tree canopy height models (CHMs) generated from both computer vision and LiDAR point clouds across sites adequately predicted field-measured tree heights, though LiDAR showed greater precision (R2 > 0.82) than computer vision (R2 > 0.64), primarily because of difficulties observing terrain under closed canopy forest. Results confirm that computer vision can support ultra-low-cost, user-deployed high spatial resolution 3D remote sensing of vegetation structure. 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>The Knoll (<b>a</b>) and Herbert Run (<b>b</b>) test sites on the campus of the University of Maryland Baltimore County. Sites and 25 m × 25 m subplots are outlined in red over 2008 leaf-off orthophotograph. Green lines delimit the approximate extent of kite aerial photograph acquisition at each site, blue crosses are GCPs used for geocorrection, and yellow circles are GCPs used in geocorrection accuracy assessment.</p>
Full article ">Figure 2
<p>Ecosynth procedure for vegetation measurements using computer vision.</p>
Full article ">Figure 3
<p>Point clouds produced by LiDAR and Ecosynth across the Knoll (<a href="#remotesensing-02-01157-f001" class="html-fig">Figure 1</a>a) and Herbert Run (<a href="#remotesensing-02-01157-f001" class="html-fig">Figure 1</a>b) test sites, compared with 2008 leaf-off orthophotograph, with 25 m × 25 m subplot grid in red (a and d). Knoll image (<b>a</b>) LiDAR first return (<b>b</b>) and Ecosynth points (<b>c</b>). Herbert Run image (<b>d</b>) LiDAR first return (<b>e</b>) and Ecosynth points (<b>f</b>). Note relief displacement of tree canopy in (d). Height colors have the same scale within each site but not across sites. Black lines delimit tree canopy determined from LiDAR.</p>
Full article ">Figure 4
<p>Oblique views of Ecosynth and LiDAR point clouds at the Knoll (<a href="#remotesensing-02-01157-f001" class="html-fig">Figure 1</a>a) and Herbert Run (<a href="#remotesensing-02-01157-f001" class="html-fig">Figure 1</a>b) test sites. Knoll aerial photograph draped on LiDAR first return (<b>a</b>), LiDAR first return plus bare earth (<b>b</b>), and Ecosynth point cloud (<b>c</b>; RGB colors). Herbert Run aerial photograph draped on LiDAR first return (<b>d</b>), LiDAR first return plus bare earth (<b>e</b>), and Ecosynth point cloud (<b>f</b>; RGB colors). 25 m subplots are outlined in purple at constant 50 m elevation. Heights in (<b>b</b>) and (<b>e</b>) use same colors as <a href="#remotesensing-02-01157-f003" class="html-fig">Figure 3</a>.</p>
Full article ">Figure 5
<p>LiDAR and Ecosynth ground points (<b>a</b>–<b>d</b>), DTMs (<b>e</b>–<b>h</b>) and DTM differences (<b>i</b> and <b>j</b>). Ground points for Knoll site from LiDAR (<b>a</b>) and Ecosynth (<b>b</b>) and Herbert Run LiDAR (<b>c</b>) and Ecosynth (<b>d</b>). DTMs from Knoll LiDAR (<b>e</b>) and Ecosynth (<b>f</b>) and from Herbert Run LiDAR (<b>g</b>) and Ecosynth (<b>h</b>). DTM differences, Ecosynth—LiDAR, for Knoll (<b>i</b>) and Herbert Run (<b>j</b>). Site orientation and height colors in (<b>a</b>) to (<b>h</b>) are same as <a href="#remotesensing-02-01157-f003" class="html-fig">Figure 3</a>. Black lines delimit tree canopy as determined from LiDAR CHM.</p>
Full article ">Figure 6
<p>Results of stepwise multiple linear regressions of subplot canopy height metrics from Ecosynth and LiDAR CHMs on field measured canopy heights. Knoll standard Ecosynth (<b>a</b>), LiDAR (<b>b</b>), and precision Ecosynth with LiDAR DTM (<b>c</b>)<span class="html-italic">.</span> Herbert Run standard Ecosynth (<b>d</b>), LiDAR (<b>e</b>), and precision Ecosynth with LiDAR DTM (<b>f</b>). Dashed lines are regression models; solid line is observed = expected. Model parameters are described in <a href="#remotesensing-02-01157-t003" class="html-table">Table 3</a>.</p>
Full article ">Figure 7
<p>Maps and density plots of differences in Ecosynth CHMs after subtracting LiDAR CHMs, in m. Knoll Ecosynth CHM (<b>a</b>) and Ecosynth CHM with LiDAR DTM (<b>b</b>). Herbert Run Ecosynth CHM (<b>c</b>) and Ecosynth CHM with LiDAR DTM (<b>d</b>). Black lines in difference maps delimit tree canopy determined from LiDAR CHM. Colors are same as <a href="#remotesensing-02-01157-f005" class="html-fig">Figure 5</a>i and <a href="#remotesensing-02-01157-f005" class="html-fig">Figure 5</a>j. Dashed vertical lines in density plots are mean difference and 1 SD from mean, solid vertical lines at 0.</p>
Full article ">
771 KiB  
Review
Ten Years of SeaWinds on QuikSCAT for Snow Applications
by Annett Bartsch
Remote Sens. 2010, 2(4), 1142-1156; https://doi.org/10.3390/rs2041142 - 16 Apr 2010
Cited by 40 | Viewed by 10338
Abstract
The scatterometer SeaWinds on QuikSCAT provided regular measurements at Ku-band from 1999 to 2009. Although it was designed for ocean applications, it has been frequently used for the assessment of seasonal snowmelt patterns aside from other terrestrial applications such as ice cap monitoring, [...] Read more.
The scatterometer SeaWinds on QuikSCAT provided regular measurements at Ku-band from 1999 to 2009. Although it was designed for ocean applications, it has been frequently used for the assessment of seasonal snowmelt patterns aside from other terrestrial applications such as ice cap monitoring, phenology and urban mapping. This paper discusses general data characteristics of SeaWinds and reviews relevant change detection algorithms. Depending on the complexity of the method, parameters such as long-term noise and multiple event analyses were incorporated. Temporal averaging is a commonly accepted preprocessing step with consideration of diurnal, multi-day or seasonal averages. Full article
(This article belongs to the Special Issue Microwave Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>SeaWinds on QuikSCAT timeseries example from autumn 2003 to spring 2004 (Salehard, 66.53<math display="inline"> <msup> <mrow/> <mo>∘</mo> </msup> </math>E, 66.53<math display="inline"> <msup> <mrow/> <mo>∘</mo> </msup> </math>N). Backscatter of all available measurements in dB (green +) compared to daily temperature range (red vertical bars show minimum to maximum in degree Celsius; source: WMO D512 dataset). Blue diamonds represent the difference between average morning and evening backscatter in dB.</p>
Full article ">Figure 2
<p>Estimated standard deviation of noise in dB (<math display="inline"> <msub> <mi>s</mi> <mi>σ</mi> </msub> </math>) of SeaWinds on QuikSCAT above 60<math display="inline"> <msup> <mrow/> <mo>∘</mo> </msup> </math>N (excluding Greenland).</p>
Full article ">Figure 3
<p>Frequency of midwinter daily average backscatter increases of more than 1.5 dB as detected with Ku-band QuikSCAT for the months November to February of winter 2000/1–2008/9 (no masking applied for lakes).</p>
Full article ">Figure 4
<p>Ten years (2000–2009) of spring snow melt dynamics based on diurnal thaw and refreeze detection [<a href="#B12-remotesensing-02-01142" class="html-bibr">12</a>]: (a) mean day of snowmelt start, (b) mean day of snowmelt end, (c) standard deviation of start of snowmelt in days, (d) standard deviation of end of snowmelt in days, (e) mean duration of spring snowmelt period, (f) standard deviation of spring snowmelt duration.</p>
Full article ">
2088 KiB  
Article
Forest Roads Mapped Using LiDAR in Steep Forested Terrain
by Russell A. White, Brian C. Dietterick, Thomas Mastin and Rollin Strohman
Remote Sens. 2010, 2(4), 1120-1141; https://doi.org/10.3390/rs2041120 - 15 Apr 2010
Cited by 86 | Viewed by 17886
Abstract
LiDAR-derived digital elevation models can reveal road networks located beneath dense forest canopy. This study tests the accuracy of forest road characteristics mapped using LiDAR in the Santa Cruz Mountains, CA. The position, gradient, and total length of a forest haul road were [...] Read more.
LiDAR-derived digital elevation models can reveal road networks located beneath dense forest canopy. This study tests the accuracy of forest road characteristics mapped using LiDAR in the Santa Cruz Mountains, CA. The position, gradient, and total length of a forest haul road were accurately extracted using a 1 m DEM. In comparison to a field-surveyed centerline, the LiDAR-derived road exhibited a positional accuracy of 1.5 m, road grade measurements within 0.53% mean absolute difference, and total road length within 0.2% of the field-surveyed length. Airborne LiDAR can provide thorough and accurate road inventory data to support forest management and watershed assessment activities. Full article
(This article belongs to the Special Issue LiDAR)
Show Figures

Figure 1

Figure 1
<p>Little Creek watershed and Swanton Pacific Ranch boundaries shown within the context of the Scotts Creek watershed.</p>
Full article ">Figure 2
<p>Little Creek road with dimensions, note the steep cut- and fill-slopes.</p>
Full article ">Figure 3
<p><b>(a)</b> LiDAR ground returns, <b>(b)</b> Distance to nearest ground return displayed with Euclidian distance grid, <b>(c)</b> 1 m color orthophoto, <b>(d)</b> bare-earth terrain represented with 1 m grayscale slope grid.</p>
Full article ">Figure 4
<p>Cumulative frequency distribution for the distance to the nearest ground return.</p>
Full article ">Figure 5
<p><b>(a)</b> Road topography degraded by over-filtering and removal of true ground returns, and <b>(b)</b> improvements in road topography after re-filtering.</p>
Full article ">Figure 6
<p><b>(a)</b> road centerline digitized from LiDAR shaded-relief grid in ESRI ArcMap, <b>(b)</b> field surveyed centerline (black) and digitized centerline (red), note change in scale.</p>
Full article ">Figure 7
<p>Distribution of LiDAR DEM error computed as DEM elevation—survey elevation at 126 road centerline locations.</p>
Full article ">Figure 8
<p>Left: DEM elevation error <span class="html-italic">versus</span> canopy cover. Grey points denote errors within high slope group. Right: DEM elevation error <span class="html-italic">versus</span> road gradient. Grey points denote errors within the high canopy cover group. Dashed vertical line denotes division between groups.</p>
Full article ">Figure 9
<p>LiDAR-derived roads in the Little Creek watershed.</p>
Full article ">Figure 10
<p>Left: Slope measurement error <span class="html-italic">versus</span> canopy cover. Gray points denote errors within the high road slope group. Right: Slope measurement error <span class="html-italic">versus</span> road gradient. Gray points denote errors within the high canopy cover group. Dashed vertical line denotes division between groups.</p>
Full article ">Figure 11
<p>Longitudinal profiles derived from ground survey and LiDAR data.</p>
Full article ">
914 KiB  
Review
Monitoring Automotive Particulate Matter Emissions with LiDAR: A Review
by Claudio Mazzoleni, Hampden D. Kuhns and Hans Moosmüller
Remote Sens. 2010, 2(4), 1077-1119; https://doi.org/10.3390/rs2041077 - 9 Apr 2010
Cited by 16 | Viewed by 14477
Abstract
Automotive particulate matter (PM) causes deleterious effects on health and visibility. Physical and chemical properties of PM also influence climate change. Roadside remote sensing of automotive emissions is a valuable option for assessing the contribution of individual vehicles to the total PM burden. [...] Read more.
Automotive particulate matter (PM) causes deleterious effects on health and visibility. Physical and chemical properties of PM also influence climate change. Roadside remote sensing of automotive emissions is a valuable option for assessing the contribution of individual vehicles to the total PM burden. LiDAR represents a unique approach that allows measuring PM emissions from in-use vehicles with high sensitivity. This publication reviews vehicle emission remote sensing measurements using ultraviolet LiDAR and transmissometer systems. The paper discusses the measurement theory and documents examples of how these techniques provide a unique perspective for exhaust emissions of individual and groups of vehicles. Full article
(This article belongs to the Special Issue LiDAR)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Diagram (<b>a</b>) and picture (<b>b</b>) of a typical field set-up system of the PM LiDAR and gaseous measurement systems. The main unit contains all the light sources and detectors, while the retro unit incorporates a set of mirrors and retro-reflectors. A trailer hosting the data acquisition and computer system is parked on the side of a single-lane road where the exhaust plume of vehicles passing by is intercepted by the light beams.</p>
Full article ">Figure 2
<p>Diagram and picture of the main unit of the UV LiDAR/transmissometer system.</p>
Full article ">Figure 3
<p>Picture of the calibration procedure and examples of data of the resulting LiDAR returns for filtered air and CO<sub>2</sub>. The insert reports the calculated overlap function.</p>
Full article ">Figure 4
<p>Gasoline vehicle passing through the LiDAR system (the license plate has been obscured for privacy reasons) deployed near Boise, Idaho, USA in winter 2004. On the right panel the excess PM mass concentrations emitted by the gasoline vehicle are represented in an intensity chart. The concentrations were calculated from the LiDAR returns after correcting for the non-linearity of the detector. Range- and extinction-corrections have been applied after the background signal has been subtracted. No correction has been applied to account for the finite temporal pulse width of the laser. The abscissa axis represents the time (in milliseconds) elapsed after the vehicle passed through the system. The ordinate axis represents the distance across the road <span class="html-italic">r</span> from the main unit (in meters). The color scale represents the PM mass concentration (in µg/m<sup>3</sup>) with white being the highest concentration; evident is the plume dispersion.</p>
Full article ">Figure 5
<p>LiDAR returns from a hard target and using different neutral density filters <span class="html-italic">versus</span> the optical density of the filters. The plot shows different curves for different values of the voltage control (gain) of the photomultiplier. The shaded box represents typical working conditions in the filed from the minimum to the maximum signals.</p>
Full article ">Figure 6
<p>Example data from a medium duty Isuzu W-series truck (model year 2000) powered by a 4.75-L diesel high-power engine. Data were taken on May 8, 2002, in Las Vegas, NV, USA (at the northbound on-ramp from N. Eastern Avenue onto I-515; 36.1684°N, 115.2487°W). (<b>a</b>) Excess PM mass column content from UV LiDAR and transmissometer and consumed fuel mass column content from CO<sub>2</sub> as a function of time. (<b>b</b>) Regression between PM and consumed fuel mass column contents yielding the PM emission factors from UV LiDAR and transmissometer as slopes in units of g of PM per kg of fuel [<a href="#B11-remotesensing-02-01077" class="html-bibr">11</a>].</p>
Full article ">Figure 7
<p>Time series for in-plume and VERSS systems for half-hour averaged PM EFs for the 4<sup>th</sup> of March 2004. Overall average PM EFs are 0.78 ± 0.09 gPM/kgFuel and 0.80 ± 0.06 gPM/kgFuel for VERSS and in-plume system, respectively [<a href="#B13-remotesensing-02-01077" class="html-bibr">13</a>].</p>
Full article ">Figure 8
<p>Panel (<b>a</b>) CO emission factors <span class="html-italic">vs.</span> the Vehicle Specific Power, panel (<b>b</b>) PM emission factors <span class="html-italic">vs.</span> the Vehicle Specific Power [<a href="#B17-remotesensing-02-01077" class="html-bibr">17</a>].</p>
Full article ">Figure 9
<p>PM emission factors <span class="html-italic">vs.</span> the vehicle age and PART5 numerical simulation [<a href="#B17-remotesensing-02-01077" class="html-bibr">17</a>].</p>
Full article ">Figure 10
<p>Average PM emission factors <span class="html-italic">vs.</span> CO (<b>a</b>) and <span class="html-italic">vs.</span> HC (<b>b</b>) emission factors [<a href="#B15-remotesensing-02-01077" class="html-bibr">15</a>].</p>
Full article ">Figure 11
<p>Contributions of lowest to highest CO emitters to the CO and PM emission factors for spark-ignition vehicles measured in Las Vegas in 2002. The error bars represent the 95% confidence interval for the respective decile mean. The percent above each column indicates the contribution to the fleet emission factor of CO or PM [<a href="#B14-remotesensing-02-01077" class="html-bibr">14</a>].</p>
Full article ">Figure 12
<p>Quintile distributions of particulate matter (PM) emission factors (EFs) for school buses with hot-stabilized and cold-start engine conditions. Panel (a) represents results obtained in hot engine conditions; panel (b) represents cold-start conditions only [<a href="#B13-remotesensing-02-01077" class="html-bibr">13</a>].</p>
Full article ">Figure 13
<p>Particulate matter (PM) emission factors (EFs) for school buses with cold-start engine conditions stratified by model year (panel (<b>a</b>)) and engine type (panel (<b>b</b>)) [<a href="#B13-remotesensing-02-01077" class="html-bibr">13</a>].</p>
Full article ">
1686 KiB  
Article
Studies on the Rapid Expansion of Sugarcane for Ethanol Production in São Paulo State (Brazil) Using Landsat Data
by Bernardo Friedrich Theodor Rudorff, Daniel Alves Aguiar, Wagner Fernando Silva, Luciana Miura Sugawara, Marcos Adami and Mauricio Alves Moreira
Remote Sens. 2010, 2(4), 1057-1076; https://doi.org/10.3390/rs2041057 - 9 Apr 2010
Cited by 336 | Viewed by 35667
Abstract
This study’s overarching aim is to establish the areal extent and characteristics of the rapid sugarcane expansion and land use change in São Paulo state (Brazil) as a result of an increase in the demand for ethanol, using Landsat type remotely sensed data. [...] Read more.
This study’s overarching aim is to establish the areal extent and characteristics of the rapid sugarcane expansion and land use change in São Paulo state (Brazil) as a result of an increase in the demand for ethanol, using Landsat type remotely sensed data. In 2003 flex fuel automobiles started to enter the Brazilian consumer market causing a dramatic expansion of sugarcane areas from 2.57 million ha in 2003 to 4.45 million ha in 2008. Almost all the land use change, for the sugarcane expansion of crop year 2008/09, occurred on pasture and annual crop land, being equally distributed on each. It was also observed that during the 2008 harvest season, the burned sugarcane area was reduced to 50% of the total harvested area in response to a protocol that aims to cease sugarcane straw burning practice by 2014 for mechanized areas. This study indicates that remote sensing images have efficiently evaluated important characteristics of the sugarcane cultivation dynamic providing quantitative results that are relevant to the debate of sustainable ethanol production from sugarcane in Brazil. Full article
(This article belongs to the Special Issue Global Croplands)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) Location of the South-Central region (highlighting São Paulo state) and the Northeast region of Brazil, (<b>b</b>) the six largest sugarcane producers in the world for the 2007/08 crop year, (<b>c</b>) percentage of sugarcane, ethanol and sugar production from Northeast and South-Central regions of Brazil, with São Paulo state presented separately.</p>
Full article ">Figure 2
<p>Multispectral [4(R)5(G)3(B)] and multi-temporal Landsat-5 images illustrating the basic procedure to identify new sugarcane fields (<b>1a</b>–<b>1d</b>) and renovated sugarcane fields with crop rotation (<b>2a</b>–<b>2d</b>).</p>
Full article ">Figure 3
<p>Examples of evident sugarcane fields harvested without straw burning that appear quite bright (<b>3.1</b> and <b>3.3a</b>) and examples of not so evident sugarcane fields harvested without straw burning that appear in greenish color (<b>3.2</b> and <b>3.3b</b>). Dark blue fields that appear in almost all figures are typical of sugarcane harvested with straw burning.</p>
Full article ">Figure 4
<p>Examples of four different land use conversion: from pasture to sugarcane (<b>4.1a</b>, <b>4.1b</b>); from agricultural crop to sugarcane (<b>4.2a</b>, <b>4.2b</b>); from citrus to sugarcane (<b>4.3a</b>, <b>4.3b</b>); and from reforestation to sugarcane (<b>4.4a</b>, <b>4.4b</b>).</p>
Full article ">Figure 5
<p>Cultivated sugarcane area available in thousands of ha from crop years 2003/04 to 2008/08 for each Administrative Regions (ADRs) and for São Paulo state.</p>
Full article ">Figure 6
<p>Cultivated sugarcane areas by Administrative Regions (ADRs) of São Paulo state for the 2008/09 crop year. The ADR of Ribeirão Preto is detached to exemplify the map of harvest type during the 2008/09 crop year. The ADR of São José do Rio Preto is detached to exemplify the map of prior land use in the expanded sugarcane fields of 2008/09 crop year.</p>
Full article ">Figure 7
<p><b>(a)</b> Percentage classes of cultivated sugarcane area and ratio classes of expansion area to renovation area; (<b>b)</b> Classes of available sugarcane area and proportions of: harvested with straw burning, harvested without straw burning and unharvested; (<b>c)</b> Classes of sugarcane expansion area and proportions of land use change for: pasture, agriculture, citrus, forest (natural vegetation) and reforestation for the Administrative Regions (ADRs) of Araçatuba (AR), Baixada Santista (BS), Barretos (BR), Bauru (BA), Campinas (CA), Central (CE), Franca (FR), Marília (MA), Presidente Prudente (PP), Registro (RE), Ribeirão Preto (RP), São José do Rio Preto (SR), São José dos Campos (SC), São Paulo (SP) and Sorocaba (SO).</p>
Full article ">
1289 KiB  
Article
Per-Field Irrigated Crop Classification in Arid Central Asia Using SPOT and ASTER Data
by Christopher Conrad, Sebastian Fritsch, Julian Zeidler, Gerd Rücker and Stefan Dech
Remote Sens. 2010, 2(4), 1035-1056; https://doi.org/10.3390/rs2041035 - 8 Apr 2010
Cited by 158 | Viewed by 15867
Abstract
The overarching goal of this research was to explore accurate methods of mapping irrigated crops, where digital cadastre information is unavailable: (a) Boundary separation by object-oriented image segmentation using very high spatial resolution (2.5–5 m) data was followed by (b) identification of crops [...] Read more.
The overarching goal of this research was to explore accurate methods of mapping irrigated crops, where digital cadastre information is unavailable: (a) Boundary separation by object-oriented image segmentation using very high spatial resolution (2.5–5 m) data was followed by (b) identification of crops and crop rotations by means of phenology, tasselled cap, and rule-based classification using high resolution (15–30 m) bi-temporal data. The extensive irrigated cotton production system of the Khorezm province in Uzbekistan, Central Asia, was selected as a study region. Image segmentation was carried out on pan-sharpened SPOT data. Varying combinations of segmentation parameters (shape, compactness, and color) were tested for optimized boundary separation. The resulting geometry was validated against polygons digitized from the data and cadastre maps, analysing similarity (size, shape) and congruence. The parameters shape and compactness were decisive for segmentation accuracy. Differences between crop phenologies were analyzed at field level using bi-temporal ASTER data. A rule set based on the tasselled cap indices greenness and brightness allowed for classifying crop rotations of cotton, winter-wheat and rice, resulting in an overall accuracy of 80 %. The proposed field-based crop classification method can be an important tool for use in water demand estimations, crop yield simulations, or economic models in agricultural systems similar to Khorezm. Full article
(This article belongs to the Special Issue Global Croplands)
Show Figures

Figure 1

Figure 1
<p>The study area. (a) Countries and irrigation regions of Central Asia. The Khorezm region is highlighted. (b) Administrative boundaries of the Rayons in the Khorezm region and the extent of the satellite images used in this study. (c) Boundaries of the Water User Associations (WUAs) in Khorezm. The SPOT 5 scenes used in this study are displayed in the background. Highlighted are the WUAs that are investigated in this study (1: Jayhun; 2: Amir Temur; 3: Shomahulum; 4: P. Mahmud; 5: Madir Yop).</p>
Full article ">Figure 2
<p>Idealized cropping calendar of the study region, Khorezm.</p>
Full article ">Figure 3
<p>Schematic workflow of the study. The left part shows the separation of field boundaries for per-field classification, the right part highlights the classification steps. Note: the segmentation of field boundaries can be replaced by any other source of vector field boundaries.</p>
Full article ">Figure 4
<p>Metrics used for segmentation suitability assessment.</p>
Full article ">Figure 5
<p>Tasseled cap brightness and greenness indices for the WUA ‘Jayhun’. Top: June data, Bottom: July data; The scatter plots at the right show brightness (x-axis) and greenness (y-axis) for the respective ASTER data sets.</p>
Full article ">Figure 6
<p>Within-field fractions of basic land cover classes (bar plots between 0 and 1) of all samples used for training the classes ‘fallow’, ‘cotton’, ‘rice’, ‘wheat-rice’ and ‘wheat’.</p>
Full article ">Figure 7
<p>Rule base for final classification at field level. Note: DS = dry soil, WS = wet soil, D&amp;MV = dense and medium vegetation, t1 = first ASTER record (01.06.2007), t2 = second ASTER record (03.07.2007)</p>
Full article ">Figure 8
<p>Assessment of SPOT segmentation settings. Note: Metric 1 and metric 2 are colored red and blue, respectively. The upper plot shows the mean values, the lower plot the standard deviations. The legend explains the labeling of the x-axes for the different settings (scale, color, compactness).</p>
Full article ">Figure 9
<p>Per-field crop map of the WUA Jayhun (close to the Amu Darya River) derived from bi-temporal ASTER data.</p>
Full article ">Figure 10
<p>Map showing the per-field crop distribution over the entire study area derived from bi-temporal ASTER data recorded in 2007.</p>
Full article ">
1384 KiB  
Article
Population Growth and Its Expression in Spatial Built-up Patterns: The Sana’a, Yemen Case Study
by Gunter Zeug and Sandra Eckert
Remote Sens. 2010, 2(4), 1014-1034; https://doi.org/10.3390/rs2041014 - 7 Apr 2010
Cited by 21 | Viewed by 16106
Abstract
In light of rapid global urbanisation, monitoring and mapping of urban and population growth is of great importance. Population growth in Sana’a was investigated for this reason. The capital of the Republic of Yemen is a rapidly growing middle sized city where the [...] Read more.
In light of rapid global urbanisation, monitoring and mapping of urban and population growth is of great importance. Population growth in Sana’a was investigated for this reason. The capital of the Republic of Yemen is a rapidly growing middle sized city where the population doubles almost every ten years. Satellite data from four different sensors were used to explore urban growth in Sana’a between 1989 and 2007, assisted by topographic maps and cadastral vector data. The analysis was conducted by delineating the built-up areas from the various optical satellite data, applying a fuzzy-rule-based composition of anisotropic textural measures and interactive thresholding. The resulting datasets were used to analyse urban growth and changes in built-up density per district, qualitatively as well as quantitatively, using a geographic information system. The built-up area increased by 87 % between 1989 and 2007. Built-up density has increased in all areas, but particularly in the northern and southern suburban districts, also reflecting the natural barrier of surrounding mountain ranges. Based on long-term population figures, geometric population growth was assumed. This hypothesis was used together with census data for 1994 and 2004 to estimate population figures for 1989 and 2007, resulting in overall growth of about 240%. By joining population figures to district boundaries, the spatial patterns of population distribution and growth were examined. Further, urban built-up growth and population changes over time were brought into relation in order to investigate changes in population density per built-up area. Population densities increased in all districts, with the greatest density change in the peripheral areas towards the North. The results reflect the pressure on the city’s infrastructure and natural resources and could contribute to sustainable urban planning in the city of Sana’a. Full article
(This article belongs to the Special Issue Multi-Temporal Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Population growth of Sana’a from 1950 to 2010 [<a href="#B8-remotesensing-02-01014" class="html-bibr">8</a>].</p>
Full article ">Figure 2
<p>Location of the study area.</p>
Full article ">Figure 3
<p>The processing and analysis workflow. (<b>a</b>) Delineation of built-up areas and change analysis. (<b>b</b>) Population distribution and change analysis. (<b>c</b>) Spatial population change analysis related to built-up expansion.</p>
Full article ">Figure 4
<p>Input and output of the PanTex workflow. (<b>a</b>) Quickbird (2004) panchromatic image resampled to 5 m spatial resolution. (<b>b</b>) Built-up area contours (red) delineated using the PanTex workflow. Grey areas represent high contrast statistics (built-up), and dark areas represent low contrast statistics (non built-up).</p>
Full article ">Figure 5
<p>Suburban (a) and urban (b) subsets from Quickbird (2004) showing correctly classified built-up (BU) (dark green) and non-built-up (NBU) (light green) areas. Errors of omission are presented in (red) and errors of commission in (orange).</p>
Full article ">Figure 6
<p>Urban expansion (1978–2007) and district map. Note that district 1306 consists of two sub-districts.</p>
Full article ">Figure 7
<p>Urban density [km<sup>2</sup>/district] in 2007, (<b>a</b>) changes in urban density by district between 1994 and 2007, (<b>b</b>) displayed as a ratio.</p>
Full article ">Figure 8
<p>Population growth between 1989 and 2007 per district together with the compound annual growth rates.</p>
Full article ">Figure 9
<p>Population density per district [Population / km<sup>2</sup>] for 2007 (<b>a</b>) and changes of population density per district between 2007 and 1994, displayed as ratio (<b>b</b>).</p>
Full article ">Figure 10
<p>Population density per built-up area of each district [Population / km<sup>2</sup>] for 2007 (<b>a</b>), and changes of population density per built-up area of each district between 2007 and 1994 displayed as ratio (<b>b</b>).</p>
Full article ">Figure 11
<p>Comparison of population densities and density change during the observed period per district. Each curve represents a district. The ordinate shows the observed period from 1989 to 2007. The abscissa indicates the number of people per district area [km<sup>2</sup>].</p>
Full article ">Figure 12
<p>Comparison of built-up densities and density change during the observed period per district. Each curve represents a district. The ordinate shows the observed period from 1989 to 2007. The abscissa indicates the built-up area per district area. A logarithmic scale was chosen to allow for a proportional comparison.</p>
Full article ">
1198 KiB  
Article
Monitoring Vegetation Phenological Cycles in Two Different Semi-Arid Environmental Settings Using a Ground-Based NDVI System: A Potential Approach to Improve Satellite Data Interpretation
by Malika Baghzouz, Dale A. Devitt, Lynn F. Fenstermaker and Michael H. Young
Remote Sens. 2010, 2(4), 990-1013; https://doi.org/10.3390/rs2040990 - 6 Apr 2010
Cited by 31 | Viewed by 11810
Abstract
In semi-arid environmental settings with sparse canopy covers, obtaining remotely sensed information on soil and vegetative growth characteristics at finer spatial and temporal scales than most satellite platforms is crucial for validating and interpreting satellite data sets. In this study, we used a [...] Read more.
In semi-arid environmental settings with sparse canopy covers, obtaining remotely sensed information on soil and vegetative growth characteristics at finer spatial and temporal scales than most satellite platforms is crucial for validating and interpreting satellite data sets. In this study, we used a ground-based NDVI system to provide continuous time series analysis of individual shrub species and soil surface characteristics in two different semi-arid environmental settings located in the Great Basin (NV, USA). The NDVI system was a dual channel SKR-1800 radiometer that simultaneously measured incident solar radiation and upward reflectance in two broadband red and near-infrared channels comparable to Landsat-5 TM band 3 and band 4, respectively. The two study sites identified as Spring Valley 1 site (SV1) and Snake Valley 1 site (SNK1) were chosen for having different species composition, soil texture and percent canopy cover. NDVI time-series of greasewood (Sarcobatus vermiculatus) from the SV1 site allowed for clear distinction between the main phenological stages of the entire growing season during the period from January to November, 2007. NDVI time series values were significantly different between sagebrush (Artemisia tridentata) and rabbitbrush (Chrysothamnus viscidiflorus) at SV1 as well as between the two bare soil types at the two sites. Greasewood NDVI from the SNK1 site produced significant correlations with chlorophyll index (r = 0.97), leaf area index (r = 0.98) and leaf xylem water potential (r = 0.93). Whereas greasewood NDVI from the SV1 site produced lower correlations (r = 0.89, r = 0.73), or non significant correlations (r = 0.32) with the same parameters, respectively. Total percent cover was estimated at 17.5% for SV1 and at 63% for SNK1. Results from this study indicated the potential capabilities of using this ground-based NDVI system to extract spatial and temporal details of soil and vegetation optical properties not possible with satellite derived NDVI. Full article
Show Figures

Figure 1

Figure 1
<p>False color Landsat 5 TM image (2007) showing the location of the two study sites (Spring Valley 1 site and Snake Valley 1 site) in the Great Basin.</p>
Full article ">Figure 2
<p>Time series of average midday (11:30–13:30 h) NDVI values of a single greasewood plant at the Spring Valley 1 site. Data were acquired for the period from January 19 to November 30, 2007. I: dormancy phase (January 19–April 2); II: active growth and canopy development phase (April 3–June 6); III: full canopy development and stable physiological status phase (June 7–August 2); IV: water limitation and stress response phase (August 3–September 26); V: leaf senescence phase (September 27–October 22); I: dormancy phase (October 23–November 30).</p>
Full article ">Figure 3
<p>Time series of average midday (11:30–13:30h) NDVI values of a single greasewood plant at Spring Valley 1 site and at Snake Valley 1 site. Data were acquired during the experimental period from May 5 to September 30, 2007. The area under the curve represents an integrated greasewood growth total of 59.31 for the Spring Valley 1 site and 65.47 for the Snake Valley 1 site.</p>
Full article ">Figure 4
<p>Time series of average midday (11:30-13:30h) NDVI values of a single sagebrush plant and a single rabbitbrush plant at Spring Valley 1 site. Data were acquired from May 5 to November 30, 2007. I: active growth and canopy development (May 5–June 4 for sagebrush) and (May 5–June 14 for rabbitbrush); II: downward adjustment entering summer period (June 5–July 10 for sagebrush) and (June 15–July 2 for rabbitbrush); III: stable physiological status during summer period (July 11–October 7 for sagebrush) and (July 3–October 7 for rabbitbrush); IV: downward adjustment entering winter period (October 8–November 30 for sagebrush) and (October 8-November 4 for rabbitbrush); V: leaf senescence phase (November 5–November 30 for rabbitbrush).</p>
Full article ">Figure 5
<p>Time series of average midday (11:30-13:30h) NDVI values of a bare soil surface at Spring Valley 1 site and at Snake Valley 1 site. Data were acquired during the experimental period from May 5 to September 30, 2007. Average daily rainfall data acquired from a weather station at the two sites are illustrated for the same experimental period.</p>
Full article ">Figure 6
<p>Comparison of plant measurements relationships with SKR-1800 NDVI values between the Snake Valley 1 site (left panel: A, B, C and D) and the Spring Valley 1 site (right panel: E, F, G and H).</p>
Full article ">Figure 6 Cont.
<p>Comparison of plant measurements relationships with SKR-1800 NDVI values between the Snake Valley 1 site (left panel: A, B, C and D) and the Spring Valley 1 site (right panel: E, F, G and H).</p>
Full article ">Figure 7
<p>Leaf moisture content for greasewood, sagebrush and rabbitbrush.</p>
Full article ">Figure 8
<p>Average midday (11:30–13:30) NDVI values and average daily ET<sub>a</sub> values between May 5 and September 30, 2007. Average midday NDVI values are both for a single greasewood plant from the Snake Valley 1 site (closed squares) and from the Spring Valley 1 site (open squares). Average daily ET<sub>a</sub> values are from Snake Valley 1 site (closed circles) and from the Spring Valley 1 site (open circles).</p>
Full article ">Figure 9
<p>Comparison of ground-based NDVI values and Landsat-NDVI acquired during satellite overpasses. Average midday ground-based NDVI values are for (A): greasewood and bare soil (closed symbols) from the Snake Valley 1 site and for (B): greasewood, sagebrush, rabbitbrush and bare soil (open symbols) from the Spring Valley 1 site. Landsat-NDVI values (closed inverted triangle) are from (A): Snake Valley 1 site and (open inverted triangle) (B): Spring Valley 1 site.</p>
Full article ">Figure 9 Cont.
<p>Comparison of ground-based NDVI values and Landsat-NDVI acquired during satellite overpasses. Average midday ground-based NDVI values are for (A): greasewood and bare soil (closed symbols) from the Snake Valley 1 site and for (B): greasewood, sagebrush, rabbitbrush and bare soil (open symbols) from the Spring Valley 1 site. Landsat-NDVI values (closed inverted triangle) are from (A): Snake Valley 1 site and (open inverted triangle) (B): Spring Valley 1 site.</p>
Full article ">
996 KiB  
Article
Comparative Analysis of Clustering-Based Approaches for 3-D Single Tree Detection Using Airborne Fullwave Lidar Data
by Sandeep Gupta, Holger Weinacker and Barbara Koch
Remote Sens. 2010, 2(4), 968-989; https://doi.org/10.3390/rs2040968 - 1 Apr 2010
Cited by 91 | Viewed by 13044
Abstract
In the past, many algorithms have been applied for three-dimensional (3-D) single tree extraction using Airborne Laser Scanner (ALS) data. Clustering based algorithms are widely used in different applications but rarely being they used in the field of forestry using ALS data as [...] Read more.
In the past, many algorithms have been applied for three-dimensional (3-D) single tree extraction using Airborne Laser Scanner (ALS) data. Clustering based algorithms are widely used in different applications but rarely being they used in the field of forestry using ALS data as an input. In this paper, a comparative qualitative study was conducted using the iterative partitioning and hierarchical clustering based mechanisms and full waveform ALS data as an input to extract the individual trees/tree crowns in their most appropriate shape. The full waveform LIght Detection And Ranging (LIDAR) data was collected from the Waldkirch black forest area in the south-western part of Germany in August 2005 with density of 4–5 points/m2. Both the clustering algorithms were used in their original and modified form for a comparative qualitative analysis of the results obtained in the form of individual clusters containing 3-D points for each tree/tree crown. A total of 378 trees were found in all the 1.2 ha area with height ranging from 15 m to 50.9 m. The forest contains mainly older trees with deciduous, coniferous and mixed stands. The findings showed that among the three kind of clustering methods applied (normal k-means, modified k-means and hierarchical clustering), the modified k-means algorithm using external seed points and scaling down the height for initialization of the clustering process was the most promising method for the extraction of clusters of individual trees/tree crowns. A 3-D reconstruction of extracted individual clusters was carried out using QHull algorithm. In this study, the result was not possible to validate quantitatively due to lack of the field inventory data. Full article
Show Figures

Figure 1

Figure 1
<p>LIDAR raw data overlaid on DSM.</p>
Full article ">Figure 2
<p>Methodology flow chart.</p>
Full article ">Figure 3
<p>LIDAR raw and normalized points. (a) LIDAR raw points of the whole test area projected above DTM. (b) LIDAR raw points above DTM in a closed view. (c) Normalized points above zero height in a closed view.</p>
Full article ">Figure 4
<p>Subdivision of study area (in progress) into a 20 m × 20 m grid.</p>
Full article ">Figure 5
<p>Divided cells (30) in yellow color overlaid on the DSM.</p>
Full article ">Figure 6
<p>Extracted local maxima overlaid on DSM (red colored points).</p>
Full article ">Figure 7
<p>Cell 6 - distribution of normalized 3-D LIDAR points, projected into a freely chosen vertical plane, at 3 different height levels (0–2 m, 2–16 m, and above 16 m) shown in 3 different colors (red, green and blue, respectively). The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 8
<p>Result after running N <span class="html-italic">k</span>-means without scaling down the height value on two different datasets of Cell 6 at different height levels. (a) Cell 6—clusters from N <span class="html-italic">k</span>-means in 2 height classes (between 2 and 16 m and above 16 m). (b) Cell 6—clusters from N <span class="html-italic">k</span>-means above 16 m height. The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 8 Cont.
<p>Result after running N <span class="html-italic">k</span>-means without scaling down the height value on two different datasets of Cell 6 at different height levels. (a) Cell 6—clusters from N <span class="html-italic">k</span>-means in 2 height classes (between 2 and 16 m and above 16 m). (b) Cell 6—clusters from N <span class="html-italic">k</span>-means above 16 m height. The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 9
<p>Result after running N <span class="html-italic">k</span>-means by scaling down the height value on two datasets of Cell 6 at different height levels. (<b>a</b>) Cell 6 – tree clusters from N <span class="html-italic">k</span>-means in two height classes (between 2 and 16 m and above 16 m). (<b>b</b>) Cell 6 – tree clusters from N <span class="html-italic">k</span>-means above 16 m height. The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 10
<p>Cluster of an individual tree from Cell 6 after running N <span class="html-italic">k</span>-means on the dataset above 16 m height and respective convex polytope, projected into a freely chosen vertical plane. (<b>a</b>) Cell 6—an individual tree cluster by applying N <span class="html-italic">k</span>-means without scaling down the height value. (<b>b</b>) 3-D convex polytope reconstructed from an individual tree cluster as shown in (<b>a</b>). (<b>c</b>) Cell 6—an individual tree cluster by applying N <span class="html-italic">k</span>-means after scaling down the height value. (<b>d</b>) 3-D convex polytope reconstructed from an individual tree cluster as shown in (c). The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 11
<p>Result after running M <span class="html-italic">k</span>-means without scaling down the height value on Cell 6 datasets of height above 16 m. The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 12
<p>Result after running M <span class="html-italic">k</span>-means after scaling down the height value on Cell 6 datasets of height above 16 m. The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 13
<p>Cluster of an individual tree from Cell 6 after running M <span class="html-italic">k</span>-means without scaling down the height value on the dataset above 16 m height and respective convex polytope. (<b>a</b>) Cell 6—an individual tree cluster above 16 m height. (<b>b</b>) 3-D Convex polytope reconstructed from tree cluster as shown in (<b>a</b>). The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 14
<p>Cluster of an individual tree from Cell 6 by applying M <span class="html-italic">k</span>-means after scaling down the height value on the dataset above 16 m height and respective convex polytope. (<b>a</b>) Cell 6—an individual tree cluster above 16 m height. (<b>b</b>) 3-D Convex polytope reconstructed from an individual tree cluster as shown in (<b>a</b>). The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 15
<p>Result after running hierarchical tree clustering without scaling down the height value on two datasets of Cell 6 at different height levels. (<b>a</b>) Cell 6 clusters after hierarchical clustering in 2 height classes (between 2 and 16 m height and above 16 m height). (<b>b</b>) Cell 6 clusters after hierarchical clustering performed on dataset above 16 m height. The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 16
<p>Result after running hierarchical tree clustering and scaling down the height value on two datasets of Cell 6 at different height levels. (<b>a</b>) Cell 6 clusters after hierarchical clustering in 2 height classes (between 2 and 16 m height and above 16 m height). (<b>b</b>) Cell 6 clusters after hierarchical clustering performed on dataset above 16 m height. The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">Figure 16 Cont.
<p>Result after running hierarchical tree clustering and scaling down the height value on two datasets of Cell 6 at different height levels. (<b>a</b>) Cell 6 clusters after hierarchical clustering in 2 height classes (between 2 and 16 m height and above 16 m height). (<b>b</b>) Cell 6 clusters after hierarchical clustering performed on dataset above 16 m height. The x and y coordinate values (in meters) are displayed horizontally and the z value (in meters) is displayed vertically.</p>
Full article ">
466 KiB  
Article
Towards Multidecadal Consistent Meteosat Surface Albedo Time Series
by Alexander Loew and Yves Govaerts
Remote Sens. 2010, 2(4), 957-967; https://doi.org/10.3390/rs2040957 - 31 Mar 2010
Cited by 38 | Viewed by 10249
Abstract
Monitoring of land surface albedo dynamics is important for the understanding of observed climate trends. Recently developed multidecadal surface albedo data products, derived from a series of geostationary satellite data, provide the opportunity to study long term surface albedo dynamics at the regional [...] Read more.
Monitoring of land surface albedo dynamics is important for the understanding of observed climate trends. Recently developed multidecadal surface albedo data products, derived from a series of geostationary satellite data, provide the opportunity to study long term surface albedo dynamics at the regional to global scale. Reliable estimates of temporal trends in surface albedo require carefully calibrated and homogenized long term satellite data records and derived products. The present paper investigates the long term consistency of a new surface albedo product derived from Meteosat First Generation (MFG) geostationary satellites for the time period 1982–2006. The temporal consistency of the data set is characterized. The analysis of the long term homogeneity reveals some discrepancies in the time series related to uncertainties in the characterization of the sensor spectral response of some of the MFG satellites. A method to compensate for uncertainties in the current data product is proposed and evaluated. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Spectral response functions of the visible channel of Meteosat First Generation radiometers.</p>
Full article ">Figure 2
<p>Surface albedo (<math display="inline"> <mrow> <mi>B</mi> <mi>H</mi> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>s</mi> <mi>o</mi> <mo>,</mo> <mi>b</mi> <mi>b</mi> </mrow> </msub> </mrow> </math>) time series for dark and bright desert targets for different spectral conversion approaches: (a) broadband albedo with standard conversion coefficients, (b) broadband with new conversion coefficients; different lines correspond to different targets.</p>
Full article ">Figure 3
<p>Time-latitude (Hovmoeller) diagrams of monthly Meteosat surface albedo anomalies based on different spectral conversion approaches: (a) broadband visible, based on original spectral conversion coefficients [<a href="#B18-remotesensing-02-00957" class="html-bibr">18</a>], (b) broadband visible based on new coefficients. Reference period for anomaly calculation: 1982–2006.</p>
Full article ">Figure 4
<p>Density plot of the relationship between <math display="inline"> <mrow> <mi>B</mi> <mi>H</mi> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>s</mi> <mi>o</mi> <mo>,</mo> <mi>λ</mi> </mrow> </msub> </mrow> </math> and <math display="inline"> <mrow> <mi>B</mi> <mi>H</mi> <msub> <mi>R</mi> <mrow> <mi>i</mi> <mi>s</mi> <mi>o</mi> <mo>,</mo> <mi>b</mi> <mi>b</mi> </mrow> </msub> </mrow> </math> in the Meteosat VIS band for all Meteosat satellites. The red dashed line corresponds to the theoretical relationship as given in the product user manual [<a href="#B18-remotesensing-02-00957" class="html-bibr">18</a>] and the green solid line corresponds to the new empirically estimated relationship. While strong corrections are applied for the Met-2/3-4 sensors, the curves for Met-5/6/7 are nearly identical for the albedo range <math display="inline"> <mrow> <mn>0</mn> <mi>…</mi> <mn>0</mn> <mo>.</mo> <mn>6</mn> </mrow> </math>.</p>
Full article ">
1516 KiB  
Article
Eucalyptus Biomass and Volume Estimation Using Interferometric and Polarimetric SAR Data
by Fábio Furlan Gama, João Roberto Dos Santos and José Claudio Mura
Remote Sens. 2010, 2(4), 939-956; https://doi.org/10.3390/rs2040939 - 31 Mar 2010
Cited by 75 | Viewed by 13133
Abstract
This work aims to establish a relationship between volume and biomass with interferometric and radiometric SAR (Synthetic Aperture Radar) response from planted Eucalyptus saligna forest stands, using multi-variable regression techniques. X and P band SAR images from the airborne OrbiSAR-1 sensor, were acquired [...] Read more.
This work aims to establish a relationship between volume and biomass with interferometric and radiometric SAR (Synthetic Aperture Radar) response from planted Eucalyptus saligna forest stands, using multi-variable regression techniques. X and P band SAR images from the airborne OrbiSAR-1 sensor, were acquired at the study area in the southeast region of Brazil. The interferometric height (Hint = difference between interferometric digital elevation model in X and P bands), contributed to the models developed due to fact that Eucalyptus forest is composed of individuals whose structure is predominantly cylindrical and vertically oriented, and whose tree heights have great correlation with volume and biomass. The volume model showed that the stand volume was highly correlated with the interferometric height logarithm (Log10Hint), since Eucalyptus tree volume has a linear relationship with the vegetation height. The biomass model showed that the combination of both Hint2 and Canopy Scattering Index—CSI (relation of s°VV by the sum of s°VV and s°HH, which represents to the canopy interaction) were used in this model, due to the fact that the Eucalyptus biomass and the trees height relationship is not linear. Both models showed a prediction error of around 10% to estimate the Eucalyptus biomass and volume, which represents a great potential to use this kind of technology to help establish Eucalyptus forest inventory for large areas. Full article
Show Figures

Figure 1

Figure 1
<p>Study area and the inventory stands.</p>
Full article ">Figure 2
<p>Methodological procedure diagram.</p>
Full article ">Figure 3
<p>Biomass.</p>
Full article ">Figure 4
<p>Tree heights.</p>
Full article ">Figure 5
<p>Total, Commercial and Canopy mean volume.</p>
Full article ">Figure 6
<p>(a) Volume and tree heights; (b) Biomass and tree heights.</p>
Full article ">Figure 7
<p>(a) Gaps and ground truth (Htopo + Hinv); (b) DEM and gaps.</p>
Full article ">Figure 8
<p>(a) Flow diagram of variables selection; (b) Graph of volume regression model.</p>
Full article ">Figure 9
<p>(a) Comparative graph between inventory volume, volume regression model results and difference between the ground truth and the model (Difference); (b) graphic of regression residues and percentage of stand gaps.</p>
Full article ">Figure 10
<p>(a) IHS image of vegetation volume model (I = X band image, H = vegetation volume, S = vegetation mask); (b) Standard deviation image.</p>
Full article ">Figure 11
<p>(a) Flow diagram of variables selection; (b) biomass regression model graph.</p>
Full article ">Figure 12
<p>(a) Biomass and CSI; (b) Biomass and Hint<sup>2</sup>.</p>
Full article ">Figure 13
<p>(a) Comparative graph of inventoried biomass, modeled biomass, and the correspondent difference; (b) Graph of stand gaps percentage and model regression residues.</p>
Full article ">Figure 14
<p>(a) IHS image from biomass model (I = X-band image, H = vegetation biomass, S = vegetation mask); (b) Standard deviation image related to biomass model.</p>
Full article ">
677 KiB  
Article
Non-Lambertian Corrected Albedo and Vegetation Index for Estimating Land Evapotranspiration in a Heterogeneous Semi-Arid Landscape
by Isabella Mariotto and Vincent P. Gutschick
Remote Sens. 2010, 2(4), 926-938; https://doi.org/10.3390/rs2040926 - 30 Mar 2010
Cited by 21 | Viewed by 9094
Abstract
The application of energy balance algorithms to remotely sensed imagery often fails to account for surface roughness variation with diverse land cover, resulting in poor resolution of evapotranspiration (ET) variations. Furthermore, the assumption of a horizontally homogeneous Lambertian surface reflecting energy equally in [...] Read more.
The application of energy balance algorithms to remotely sensed imagery often fails to account for surface roughness variation with diverse land cover, resulting in poor resolution of evapotranspiration (ET) variations. Furthermore, the assumption of a horizontally homogeneous Lambertian surface reflecting energy equally in all directions affects the calculations of albedo and vegetation index. The primary objective of this study is to improve the accuracy of the estimation and discrimination of ET among different land cover types in Southern New Mexico from ASTER datasets, by formulating the spatial variation of non-Lambertian reflectance using a wavelength-dependent Minnaert function. Full article
Show Figures

Figure 1

Figure 1
<p>ASTER surface reflectance VNIR band 3N over Las Cruces area on 08 June 2005.</p>
Full article ">Figure 2
<p>Eddy-Covariance ET flux data at the JER in June 2005.</p>
Full article ">Figure 3
<p>Minnaert regressions for each land cover class and each ASTER reflectance band at the JER.</p>
Full article ">Figure 4
<p>Comparison of reflectance (a) before and (b) after correction.</p>
Full article ">Figure 5
<p>Comparison of MSAVI (<b>a</b>) before and (<b>b</b>) after reflectance correction.</p>
Full article ">Figure 6
<p>Final (a) non-reflectance corrected and (b) reflectance corrected ET (mm/day) maps, which show a pattern with (c) surface temperature and (d) land cover type, respectively.</p>
Full article ">Figure 7
<p>TIN from contours. Playa grassland is located at the lowest elevation of the area.</p>
Full article ">Figure 8
<p>Scatterplot of ET single pixels values (mm/day) of corrected SEBAL (y) <span class="html-italic">versus</span> non-corrected SEBAL (x).</p>
Full article ">Figure 9
<p>Scatterplots of ET single pixels values (mm/day) for (a) playa grassland, (b) creosotebush scrub, (c) tarbush shrubland, and (d,e) mesquite duneland for corrected model (y) and non-corrected (x).</p>
Full article ">Figure 10
<p>Mean ET (mm/day) of each land cover type.</p>
Full article ">
590 KiB  
Article
Using Spatial Structure Analysis of Hyperspectral Imaging Data and Fourier Transformed Infrared Analysis to Determine Bioactivity of Surface Pesticide Treatment
by Christian Nansen, Noureddine Abidi, Amelia Jorge Sidumo and Ali Hosseini Gharalari
Remote Sens. 2010, 2(4), 908-925; https://doi.org/10.3390/rs2040908 - 26 Mar 2010
Cited by 16 | Viewed by 11565
Abstract
Many food products are subjected to quality control analyses for detection of surface residue/contaminants, and there is a trend of requiring more and more documentation and reporting by farmers regarding their use of pesticides. Recent outbreaks of food borne illnesses have been a [...] Read more.
Many food products are subjected to quality control analyses for detection of surface residue/contaminants, and there is a trend of requiring more and more documentation and reporting by farmers regarding their use of pesticides. Recent outbreaks of food borne illnesses have been a major contributor to this trend. With a growing need for food safety measures and “smart applications” of insecticides, it is important to develop methods for rapid and accurate assessments of surface residues on food and feed items. As a model system, we investigated detection of a miticide applied to maize leaves and its miticidal bioactivity over time, and we compared two types of reflectance data: fourier transformed infrared (FTIR) data and hyperspectral imaging (HI) data. The miticide (bifenazate) was applied at a commercial field rate to maize leaves in the field, with or without application of a surfactant, and with or without application of a simulated “rain event”. In addition, we collected FTIR and HI from untreated control leaves (total of five treatments). Maize leaf data were collected at seven time intervals from 0 to 48 hours after application. FTIR data were analyzed using conventional analysis of variance of miticide-specific vibration peaks. Two unique FTIR vibration peaks were associated with miticide application (1,700 cm−1 and 763 cm−1). The integrated intensities of these two peaks, miticide application, surfactant, rain event, time between miticide application, and rain event were used as explanatory variables in a linear multi-regression fit to spider mite mortality. The same linear multi-regression approach was applied to variogram parameters derived from HI data in five selected spectral bands (664, 683, 706, 740, and 747 nm). For each spectral band, we conducted a spatial structure analysis, and the three standard variogram parameters (“sill”, “range”, and “nugget”) were examined as possible “indicators” of miticide bioactivity. We demonstrated that both FTIR peaks and standard variogram parameters could be used to accurately predict spider mite mortality, but linear multi-regression fits based on standard variogram parameters had the highest accuracy and were successfully validated with independent data. Based on experimental manipulation of HI data, the use of spatial structure analysis in classification of HI data was discussed. Full article
(This article belongs to the Special Issue Ecological Status and Change by Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Variogram illustrates relationship of distance between paired observations (lag distance) and variance, and variogram analysis is used to determine three parameters (“Nugget”, “Range”, and “Sill”).</p>
Full article ">Figure 2
<p>Average miticide bioactivity (spider mite mortality) from different treatments. Statistical analyses of spider mite bioassay results are presented in <a href="#remotesensing-02-00908-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 3
<p>Fourier transformed infrared analysis (FTIR) spectra from maize leaves. Untreated control maize leaf (in green), maize leaf treated with miticide before rain event (blue), and maize leaf treated with miticide after simulated rain (red).</p>
Full article ">Figure 4
<p>Average hyperspectral profiles acquired from untreated and miticide-treated maize leaves over time.</p>
Full article ">Figure 5
<p>Variogram analysis of reflectance data from single spectral band (683 nm) of an untreated maize leaf were experimentally manipulated in four ways and compared with actual: multiplying all reflectance values with either 1.025 or 1.050 to simulate a 2.5% and 5.0% increase in all reflectance values or multiplying half the reflectance values (random selection) with 1.025 or one-third of the reflectance values (random selection) with 1.050.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop