[go: up one dir, main page]

Next Issue
Volume 7, May
Previous Issue
Volume 7, March
 
 
remotesensing-logo

Journal Browser

Journal Browser

Remote Sens., Volume 7, Issue 4 (April 2015) – 70 articles , Pages 3426-4972

  • 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:
1442 KiB  
Editorial
Innovative Technologies for Terrestrial Remote Sensing
by Paul Aplin and Doreen S. Boyd
Remote Sens. 2015, 7(4), 4968-4972; https://doi.org/10.3390/rs70404968 - 22 Apr 2015
Cited by 1 | Viewed by 6540
Abstract
Characterizing and monitoring terrestrial, or land, surface features, such as forests, deserts, and cities, are fundamental and continuing goals of Earth Observation (EO). EO imagery and related technologies are essential for increasing our scientific understanding of environmental processes, such as carbon capture and [...] Read more.
Characterizing and monitoring terrestrial, or land, surface features, such as forests, deserts, and cities, are fundamental and continuing goals of Earth Observation (EO). EO imagery and related technologies are essential for increasing our scientific understanding of environmental processes, such as carbon capture and albedo change, and to manage and safeguard environmental resources, such as tropical forests, particularly over large areas or the entire globe. This measurement or observation of some property of the land surface is central to a wide range of scientific investigations and industrial operations, involving individuals and organizations from many different backgrounds and disciplines. However, the process of observing the land provides a unifying theme for these investigations, and in practice there is much consistency in the instruments used for observation and the techniques used to map and model the environmental phenomena of interest. There is therefore great potential benefit in exchanging technological knowledge and experience among the many and diverse members of the terrestrial EO community. [...] Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Earth Observation Technology Cluster word cloud, showing the initiative’s principal EO technology connections and collaborations.</p>
Full article ">
2116 KiB  
Article
Knowledge-Based Detection and Assessment of Damaged Roads Using Post-Disaster High-Resolution Remote Sensing Image
by Jianhua Wang, Qiming Qin, Jianghua Zhao, Xin Ye, Xiao Feng, Xuebin Qin and Xiucheng Yang
Remote Sens. 2015, 7(4), 4948-4967; https://doi.org/10.3390/rs70404948 - 22 Apr 2015
Cited by 24 | Viewed by 8660
Abstract
Road damage detection and assessment from high-resolution remote sensing image is critical for natural disaster investigation and disaster relief. In a disaster context, the pairing of pre-disaster and post-disaster road data for change detection and assessment is difficult to achieve due to the [...] Read more.
Road damage detection and assessment from high-resolution remote sensing image is critical for natural disaster investigation and disaster relief. In a disaster context, the pairing of pre-disaster and post-disaster road data for change detection and assessment is difficult to achieve due to the mismatch of different data sources, especially for rural areas where the pre-disaster data (i.e., remote sensing imagery or vector map) are hard to obtain. In this study, a knowledge-based method for road damage detection and assessment solely from post-disaster high-resolution remote sensing image is proposed. The road centerline is firstly extracted based on the preset road seed points. Then, features such as road brightness, standard deviation, rectangularity, and aspect ratio are selected to form a knowledge model. Finally, under the guidance of the road centerline, the post-disaster roads are extracted and the damaged roads are detected by applying the knowledge model. In order to quantitatively assess the damage degree, damage assessment indicators with their corresponding standard of damage grade are also proposed. The newly developed method is evaluated using a WorldView-1 image over Wenchuan, China acquired three days after the earthquake on 15 May 2008. The results show that the producer’s accuracy (PA) and user’s accuracy (UA) reached about 90% and 85%, respectively, indicating that the proposed method is effective for road damage detection and assessment. This approach also significantly reduces the need for pre-disaster remote sensing data. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Flow chart of road damage detection and assessment.</p>
Full article ">Figure 2
<p>The diagram of road. (<b>a</b>) The road before disaster; (<b>b</b>) The road after disaster.</p>
Full article ">Figure 3
<p>Schematic diagram of post-disaster road extraction. <span class="html-italic">P</span> is a group of pixels located in the road centerline (<math display="inline"> <semantics> <mrow> <mi>P</mi> <mo>=</mo> <mo>{</mo> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>p</mi> <mn>2</mn> </msub> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msub> <mi>p</mi> <mi>n</mi> </msub> <mo>}</mo> </mrow> </semantics> </math>, <span class="html-italic">n</span> is the total number of pixels). <span class="html-italic">l</span> represents the searching line which moves along the road centerline and the length is D (<span class="html-italic">D</span> ≥ <span class="html-italic">w<sub>road</sub></span>).</p>
Full article ">Figure 4
<p>Road extraction. (<b>a</b>) The test image, which is an urban image without damage; (<b>b</b>) The hypothetic roads; (<b>c</b>) The roads after verification.</p>
Full article ">Figure 5
<p>Schematic diagram of road damage detection.</p>
Full article ">Figure 6
<p>The spatial location of the study area.</p>
Full article ">Figure 7
<p>Image of WorldView-1 in the study area.</p>
Full article ">Figure 8
<p>Result of road centerline extraction. The road centerline is shown as red line and the road seed points are shown as yellow crosses.</p>
Full article ">Figure 9
<p>Results of post-disaster road extraction. Green regions are the post-disaster roads.</p>
Full article ">Figure 10
<p>Results of damage detection. Red regions are the damaged road segments.</p>
Full article ">Figure 11
<p>Tree and vehicle that are mistaken as damaged roads. (<b>a</b>) Tree shadow; (<b>b</b>) Tree shadow is mistaken as damaged road; (<b>c</b>) Vehicle; (<b>d</b>) Vehicle is mistaken as damaged road.</p>
Full article ">Figure 12
<p>Sensitivity test of free parameters. (<b>a</b>) the brightness threshold <math display="inline"> <semantics> <mrow> <msub> <mi>b</mi> <mn>1</mn> </msub> </mrow> </semantics> </math>; (<b>b</b>) the brightness threshold <math display="inline"> <semantics> <mrow> <msub> <mi>b</mi> <mn>2</mn> </msub> </mrow> </semantics> </math>; (<b>c</b>) the standard deviation threshold <math display="inline"> <semantics> <mrow> <msub> <mi>s</mi> <mn>1</mn> </msub> </mrow> </semantics> </math>; (<b>d</b>) the standard deviation threshold <math display="inline"> <semantics> <mrow> <msub> <mi>s</mi> <mn>2</mn> </msub> </mrow> </semantics> </math>; (<b>e</b>) the rectangularity threshold <math display="inline"> <semantics> <mrow> <msub> <mi>r</mi> <mn>1</mn> </msub> </mrow> </semantics> </math>; (<b>f</b>) the length-to-width ratio threshold <math display="inline"> <semantics> <mrow> <msub> <mi>w</mi> <mn>1</mn> </msub> </mrow> </semantics> </math>.</p>
Full article ">Figure 13
<p>Result of damage grade identification. The damage grades of the green, blue, yellow and red road segments are basic, minor, moderate and major, respectively.</p>
Full article ">
716 KiB  
Article
A Test of the New VIIRS Lights Data Set: Population and Economic Output in Africa
by Xi Chen and William Nordhaus
Remote Sens. 2015, 7(4), 4937-4947; https://doi.org/10.3390/rs70404937 - 22 Apr 2015
Cited by 98 | Viewed by 10783
Abstract
The present study analyses the new Visible Infrared Imaging Radiometer Suite (VIIRS) lights data to determine whether it can provide more accurate proxies for socioeconomic data in areas with poor quality data than proxies based on stable lights. Our analysis indicates that VIIRS [...] Read more.
The present study analyses the new Visible Infrared Imaging Radiometer Suite (VIIRS) lights data to determine whether it can provide more accurate proxies for socioeconomic data in areas with poor quality data than proxies based on stable lights. Our analysis indicates that VIIRS lights are a promising supplementary source for standard measures on population and economic output at a small scale, especially for low population and economic density areas in Africa. The current analysis also suggests that in comparison to stable lights generated by the DMSP-OLS system, data generated by the VIIRS system provide more information to estimate population than output index. However, further analysis and formal statistical models are needed to evaluate the usefulness of VIIRS lights versus other lights products. With more advanced methods, there is also a potential to generate a synthetic index by combining different lights products to produce a better proxy measure for other indexes. Full article
(This article belongs to the Special Issue Remote Sensing with Nighttime Lights)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Scatter plot of cell population, stable lights, and Visible Infrared Imaging Radiometer Suite (VIIRS) lights for cells with a population &lt; 1000.</p>
Full article ">
4028 KiB  
Article
Users’ Assessment of Orthoimage Photometric Quality for Visual Interpretation of Agricultural Fields
by Agnieszka Tarko, Sytze De Bruin, Dominique Fasbender, Wim Devos and Arnold K. Bregt
Remote Sens. 2015, 7(4), 4919-4936; https://doi.org/10.3390/rs70404919 - 21 Apr 2015
Cited by 7 | Viewed by 5813
Abstract
Land cover identification and area quantification are key aspects of implementing the European Common Agriculture Policy. Legitimacy of support provided to farmers is monitored using the Land Parcel Identification System (LPIS), with land cover identification performed by visual image interpretation. While the geometric [...] Read more.
Land cover identification and area quantification are key aspects of implementing the European Common Agriculture Policy. Legitimacy of support provided to farmers is monitored using the Land Parcel Identification System (LPIS), with land cover identification performed by visual image interpretation. While the geometric orthoimage quality required for correct interpretation is well understood, little is known about the photometric quality needed for LPIS applications. This paper analyzes the orthoimage quality characteristics chosen by authors as being most suitable for visual identification of agricultural fields. We designed a survey to assess users’ preferred brightness and contrast ranges for orthoimages used for LPIS purposes. Survey questions also tested the influence of a background color on the preferred orthoimage brightness and contrast, the preferred orthoimage format and color composite, assessments of orthoimages with shadowed areas, appreciation of image enhancements and, finally, consistency of individuals’ preferred brightness and contrast settings across multiple sample images. We find that image appreciation is stable at the individual level, but preferences vary across respondents. We therefore recommend that LPIS operators be enabled to personalize photometric settings, such as brightness and contrast values, and to choose the displayed band combination from at least four spectral bands. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The orthoimages selected for the online survey are from the zones outlined in yellow (indicated with a black pin).</p>
Full article ">Figure 2
<p>Two questions concerning brightness and contrast using the same images but presenting them on different background colors—white (question 3, <b>Left</b>) and black (question 11, <b>Right</b>). The orthoimages displayed have the following properties (brightness and contrast): A (-30, 0), B (0, 0) and C (15, 0).</p>
Full article ">Figure 3
<p>Preferred modification levels for brightness (<b>Left</b>) and contrast (<b>Right</b>).</p>
Full article ">Figure 4
<p>Respondents’ choice of the best and the worst images against the white and the black background (the brightness and contrast combinations used in the sample images are specified). Figures A, B and C each represent a pair of questions using the same sample images against a different background color (best on the left and worst on the right; white column = white background, black column = black background).</p>
Full article ">Figure 5
<p>Estimated Shannon entropy (points) and its standard deviation times the 97.5% quantile of the Student distribution with 196 degrees of freedom (whiskers) for choices of the most and the least appreciated image triplets, paired with the white background (gray line) and black background (black line). The pairs of questions referred to as A, B and C are the same as those in <a href="#remotesensing-07-04919-f004" class="html-fig">Figure 4</a>.</p>
Full article ">Figure 6
<p>Influence of best image format choice on the subsequent worst one (FCC = false color composite, NCC = natural color composite, TIFF = tagged image file format, ECW = Enhanced Compression Wavelet).</p>
Full article ">Figure 7
<p>Number of choices (indicated by the marker size) for most and least preferred brightness and contrast. The gray dashed line indicates the theoretically most consistent choice possibilities.</p>
Full article ">
5772 KiB  
Article
Comparative Assessment of Satellite-Retrieved Surface Net Radiation: An Examination on CERES and SRB Datasets in China
by Xin Pan, Yuanbo Liu and Xingwang Fan
Remote Sens. 2015, 7(4), 4899-4918; https://doi.org/10.3390/rs70404899 - 21 Apr 2015
Cited by 33 | Viewed by 7019
Abstract
Surface net radiation plays an important role in land–atmosphere interactions. The net radiation can be retrieved from satellite radiative products, yet its accuracy needs comprehensive assessment. This study evaluates monthly surface net radiation generated from the Clouds and the Earth’s Radiant Energy System [...] Read more.
Surface net radiation plays an important role in land–atmosphere interactions. The net radiation can be retrieved from satellite radiative products, yet its accuracy needs comprehensive assessment. This study evaluates monthly surface net radiation generated from the Clouds and the Earth’s Radiant Energy System (CERES) and the Surface Radiation Budget project (SRB) products, respectively, with quality-controlled radiation data from 50 meteorological stations in China for the period from March 2000 to December 2007. Our results show that surface net radiation is generally overestimated for CERES (SRB), with a bias of 26.52 W/m2 (18.57 W/m2) and a root mean square error of 34.58 W/m2 (29.49 W/m2). Spatially, the satellite-retrieved monthly mean of surface net radiation has relatively small errors for both CERES and SRB at inland sites in south China. Substantial errors are found at northeastern sites for two datasets, in addition to coastal sites for CERES. Temporally, multi-year averaged monthly mean errors are large at sites in western China in spring and summer, and in northeastern China in spring and winter. The annual mean error fluctuates for SRB, but decreases for CERES between 2000 and 2007. For CERES, 56% of net radiation errors come from net shortwave (NSW) radiation and 44% from net longwave (NLW) radiation. The errors are attributable to environmental parameters including surface albedo, surface water vapor pressure, land surface temperature, normalized difference vegetation index (NDVI) of land surface proxy, and visibility for CERES. For SRB, 65% of the errors come from NSW and 35% from NLW radiation. The major influencing factors in a descending order are surface water vapor pressure, surface albedo, land surface temperature, NDVI, and visibility. Our findings offer an insight into error patterns in satellite-retrieved surface net radiation and should be valuable to improving retrieval accuracy of surface net radiation. Moreover, our study on radiation data of China provides a case example for worldwide validation. Full article
(This article belongs to the Special Issue Remote Sensing of Solar Surface Radiation)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Distribution of Chinese meteorological sites used in the present study.</p>
Full article ">Figure 2
<p>Comparison of monthly net radiation derived from the CERES and the SRB products with surface observations at 50 meteorological sites in China. (<b>a</b>) Surface observations <span class="html-italic">versus</span> CERES data; (<b>b</b>) surface observations <span class="html-italic">versus</span> SRB data; and (<b>c</b>) SRB <span class="html-italic">versus</span> CERES.</p>
Full article ">Figure 3
<p>Error distribution of monthly surface net radiation generated from the CERES and SRB data for the period from March 2000 to December 2007.</p>
Full article ">Figure 4
<p>Statistical comparison of CERES, SRB, and observed (OBS) surface-measured monthly net radiation.</p>
Full article ">Figure 5
<p>Spatial distribution of multi-year averaged errors in monthly surface net radiation generated from the CERES and SRB products.</p>
Full article ">Figure 6
<p>Monthly variation of errors in surface net radiation generated from the CERES and the SRB products.</p>
Full article ">Figure 7
<p>Spatial distribution of annual mean errors in surface net radiation generated from the CERES and the SRB products from 2000 to 2007.</p>
Full article ">Figure 8
<p>Annual mean errors in surface net radiation generated from the CERES and the SRB products from 2000 to 2007.</p>
Full article ">Figure 9
<p>Correlation coefficients (R<sup>2</sup>) between environmental parameters and errors in NSW and NLW for CERES and SRB at 11 level-1 sites.</p>
Full article ">Figure 10
<p>Spatial distributions of correlation coefficients (R<sup>2</sup>) between environmental parameters and errors in NSW and NLW for CERES and SRB at 50 sites.</p>
Full article ">
11535 KiB  
Article
Assessing the Impacts of Urbanization-Associated Land Use/Cover Change on Land Surface Temperature and Surface Moisture: A Case Study in the Midwestern United States
by Yitong Jiang, Peng Fu and Qihao Weng
Remote Sens. 2015, 7(4), 4880-4898; https://doi.org/10.3390/rs70404880 - 20 Apr 2015
Cited by 98 | Viewed by 11378
Abstract
Urbanization-associated land use and land cover (LULC) changes lead to modifications of surface microclimatic and hydrological conditions, including the formation of urban heat islands and changes in surface runoff pattern. The goal of the paper is to investigate the changes of biophysical variables [...] Read more.
Urbanization-associated land use and land cover (LULC) changes lead to modifications of surface microclimatic and hydrological conditions, including the formation of urban heat islands and changes in surface runoff pattern. The goal of the paper is to investigate the changes of biophysical variables due to urbanization induced LULC changes in Indianapolis, USA, from 2001 to 2006. The biophysical parameters analyzed included Land Surface Temperature (LST), fractional vegetation cover, Normalized Difference Water Index (NDWI), impervious fractions evaporative fraction, and soil moisture. Land cover classification and changes and impervious fractions were obtained from the National Land Cover Database of 2001 and 2006. The Temperature-Vegetation Index (TVX) space was created to analyze how these satellite-derived biophysical parameters change during urbanization. The results showed that the general trend of pixel migration in response to the LULC changes was from the areas of low temperature, dense vegetation cover, and high surface moisture conditions to the areas of high temperature, sparse vegetation cover, and low surface moisture condition in the TVX space. Analyses of the T-soil moisture and T-NDWI spaces revealed similar changed patterns. The rate of change in LST, vegetation cover, and moisture varied with LULC type and percent imperviousness. Compared to conversion from cultivated to residential land, the change from forest to commercial land altered LST and moisture more intensively. Compared to the area changed from cultivated to residential, the area changed from forest to commercial altered 48% more in fractional vegetation cover, 71% more in LST, and 15% more in soil moisture Soil moisture and NDWI were both tested as measures of surface moisture in the urban areas. NDWI was proven to be a useful measure of vegetation liquid water and was more sensitive to the land cover changes comparing to soil moisture. From a change forest to commercial land, the mean soil moisture changed 17%, while the mean NDWI changed 90%. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Three sample areas for assessing the impact of urbanization in the City of Indianapolis, USA. The vegetation fraction was produced by the Landsat TM image that acquired on 16 June 2001. (<b>a</b>) The whole study area, (<b>b</b>) Area 1 is characterized by conversions from cultivated to residential lands, (<b>c</b>) Area 2 and Area 3 represent changes from forest to commercial, and from open area to commercial respectively.</p>
Full article ">Figure 2
<p>The flowchart of the study. The blue polygons with white text are the resultant biophysical parameters that were compared between 2001 and 2006.</p>
Full article ">Figure 3
<p>Method to interpolate <math display="inline"><semantics> <mrow><msub><mtext>φ</mtext><mtext>i</mtext></msub></mrow> </semantics></math> for each pixel. <math display="inline"> <semantics> <mrow> <msub> <mtext>φ</mtext> <mtext>i</mtext> </msub> </mrow> </semantics> </math> is the <math display="inline"> <semantics> <mi>φ</mi> </semantics> </math> value for each random pixel in the TVX space. In the similar triangles ABC and ADE, AC/AE = BC/DE. <math display="inline"> <semantics> <mrow> <msub> <mtext>φ</mtext> <mtext>i</mtext> </msub> </mrow> </semantics> </math> can be calculated using the values of <math display="inline"> <semantics> <mrow> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> <mo>,</mo> <mo> </mo> <msub> <mi>T</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>,</mo> <mo> </mo> <msub> <mi>T</mi> <mi>i</mi> </msub> <mo>,</mo> <mo> </mo> <msub> <mi>φ</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo> </mo> </mrow> </semantics> </math> and <math display="inline"> <semantics> <mrow> <msub> <mi>φ</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> </semantics> </math>.</p>
Full article ">Figure 4
<p>Land cover, the percentage of impervious surface, scaled vegetation fraction, scaled Land Surface Temperature (LST), soil moisture, and Normalized Difference Water Index (NDWI) [<a href="#B15-remotesensing-07-04880" class="html-bibr">15</a>] in 2001 and 2006, and the percentage of changes of each parameter from 2001 to 2006. Land cover change and percent of imperviousness were acquired from National Land Cover Database (NLCD), while other parameters were generated from this study.</p>
Full article ">Figure 5
<p>The characteristics measured by scaled vegetation fraction (Scaled Fr), scaled LST, soil moisture, NDWI and percentage of imperviousness for different land types in three sample areas on 17 June 2001 and 1 July 2006.</p>
Full article ">Figure 6
<p>Scatterplot of scaled LST (x axis) <span class="html-italic">versus</span> scaled fractional vegetation cover (y axis) for Landsat TM images that was acquired on 17 June 2001 and 1 July 2006. Compared to the shape of the scatter plot in 2001(<b>upper</b>), 2006 one (<b>lower</b>) became “shorter” and “wider”, which indicated the general trend of the surface condition changed to lower vegetation cover, lower moisture availability, and higher temperature. The scaled LST was transformed from LST by Equation (8) using maximum, minimum, and average LST. The scaled fractional vegetation cover was transformed by the same method.</p>
Full article ">Figure 7
<p>(<b>a</b>) Pixel trajectories in the TVX space, (<b>b</b>) Temperature-soil moisture space, and (<b>c</b>) Temperature-NDWI space from 17 June 2001 to 1 July 2006. Cultivated to residential was represented by Area 1, Forest to commercial was represented by Area 1, and open area to commercial was represented by Area 3.</p>
Full article ">
46047 KiB  
Article
Hydrodynamic and Inundation Modeling of China’s Largest Freshwater Lake Aided by Remote Sensing Data
by Peng Zhang, Jianzhong Lu, Lian Feng, Xiaoling Chen, Li Zhang, Xiongwu Xiao and Honggao Liu
Remote Sens. 2015, 7(4), 4858-4879; https://doi.org/10.3390/rs70404858 - 20 Apr 2015
Cited by 34 | Viewed by 8761
Abstract
China’s largest freshwater lake, Poyang Lake, is characterized by rapid changes in its inundation area and hydrodynamics, so in this study, a hydrodynamic model of Poyang Lake was established to simulate these long-term changes. Inundation information was extracted from Moderate Resolution Imaging Spectroradiometer [...] Read more.
China’s largest freshwater lake, Poyang Lake, is characterized by rapid changes in its inundation area and hydrodynamics, so in this study, a hydrodynamic model of Poyang Lake was established to simulate these long-term changes. Inundation information was extracted from Moderate Resolution Imaging Spectroradiometer (MODIS) remote sensing data and used to calibrate the wetting and drying parameter by assessing the accuracy of the simulated inundation area and its boundary. The bottom friction parameter was calibrated using current velocity measurements from Acoustic Doppler Current Profilers (ADCP). The results show the model is capable of predicting the inundation area dynamic through cross-validation with remotely sensed inundation data, and can reproduce the seasonal dynamics of the water level, and water discharge through a comparison with hydrological data. Based on the model results, the characteristics of the current velocities of the lake in the wet season and the dry season of the lake were explored, and the potential effect of the current dynamic on water quality patterns was discussed. The model is a promising basic tool for prediction and management of the water resource and water quality of Poyang Lake. Full article
(This article belongs to the Special Issue Remote Sensing in Flood Monitoring and Management)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Study area and measurement sites.</p>
Full article ">Figure 2
<p>Flow chart of the methodologies.</p>
Full article ">Figure 3
<p>Model grids (<b>A</b>) and bathmetry (<b>B</b>).</p>
Full article ">Figure 4
<p>RMSE values of the simulated inundation areas, the mean Kappa coefficient for the maximum remotely sensed inundation area (MKMAX) and mean Kappa coefficient for the minimum remotely sensed inundation area (MKMIN) during the model calibration period (2001–2005) for calibration experiments with different critical water depth.</p>
Full article ">Figure 5
<p>Comparison of simulated and remotely sensed inundation areas in the calibration period (2001–2005) and the validation period (2006–2010).</p>
Full article ">Figure 6
<p>Model simulated water boundary (Model_Boundary) and remotely sensed inundation area (RS_Area) at the time of maximum (Max) and minimum (Min) inundation area for each year from 2006 to 2010.</p>
Full article ">Figure 7
<p>Comparison of simulated and measured velocities along two separated ADCP tracks.</p>
Full article ">Figure 8
<p>Comparison of water level and water discharge from the model and the hydrological stations during the calibration periods (2001–2005) and the validation periods (2006–2010). Note that the water level measurements at Tangyin and Longkou are not available in the year 2010.</p>
Full article ">Figure 9
<p>Monthly mean depth-averaged current vectors and magnitudes in the dry season (June, 2005) and the wet season (November, 2005). (<b>A</b>) and (<b>B</b>) are current vectors in the dry season and wet season, respectively, and (<b>C</b>) and (<b>D</b>) are current magnitudes in the dry season and wet season, respectively.</p>
Full article ">
1624 KiB  
Article
Comparing the Dry Season In-Situ Leaf Area Index (LAI) Derived from High-Resolution RapidEye Imagery with MODIS LAI in a Namibian Savanna
by Manuel J. Mayr and Cyrus Samimi
Remote Sens. 2015, 7(4), 4834-4857; https://doi.org/10.3390/rs70404834 - 20 Apr 2015
Cited by 24 | Viewed by 7618
Abstract
The Leaf Area Index (LAI) is one of the most frequently applied measures to characterize vegetation and its dynamics and functions with remote sensing. Satellite missions, such as NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS) operationally produce global datasets of LAI. Due to their [...] Read more.
The Leaf Area Index (LAI) is one of the most frequently applied measures to characterize vegetation and its dynamics and functions with remote sensing. Satellite missions, such as NASA’s Moderate Resolution Imaging Spectroradiometer (MODIS) operationally produce global datasets of LAI. Due to their role as an input to large-scale modeling activities, evaluation and verification of such datasets are of high importance. In this context, savannas appear to be underrepresented with regards to their heterogeneous appearance (e.g., tree/grass-ratio, seasonality). Here, we aim to examine the LAI in a heterogeneous savanna ecosystem located in Namibia’s Owamboland during the dry season. Ground measurements of LAI are used to derive a high-resolution LAI model with RapidEye satellite data. This model is related to the corresponding MODIS LAI/FPAR (Fraction of Absorbed Photosynthetically Active Radiation) scene (MOD15A2) in order to evaluate its performance at the intended annual minimum during the dry season. Based on a field survey we first assessed vegetation patterns from species composition and elevation for 109 sites. Secondly, we measured in situ LAI to quantitatively estimate the available vegetation (mean = 0.28). Green LAI samples were then empirically modeled (LAImodel) with high resolution RapidEye imagery derived Difference Vegetation Index (DVI) using a linear regression (R2 = 0.71). As indicated by several measures of model performance, the comparison with MOD15A2 revealed moderate consistency mostly due to overestimation by the aggregated LAImodel. Model constraints aside, this study may point to important issues for MOD15A2 in savannas concerning the underlying MODIS Land Cover product (MCD12Q1) and a potential adjustment by means of the MODIS Burned Area product (MCD45A1). Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Map of Northern-Central Namibia illustrating the study area and the Cuvelai catchment.</p>
Full article ">Figure 2
<p>Schematic illustration of spatial sampling for ESU I66: measurements were made at 5 m intervals along two perpendicularly intersecting transects. Points (orange) indicate the 15 measurements below the canopy in ESU I66 (green). Directions of the sampling process are indicated by the white arrows.</p>
Full article ">Figure 3
<p>Exemplary sites of green vegetation in the study area: (<b>left</b>) <span class="html-italic">Colophospermum mopane</span> shrub lands. (<b>right</b>) Open woodlands, mainly containing Makalani palms (<span class="html-italic">Hyphaene petersiana</span>).</p>
Full article ">Figure 4
<p>Conditional plots for LAI<sub>2200</sub> and estimated total plant cover in %, per terrain position class. The black circles denote the samples from an ESU, the red lines show the respective LOESS smoothing lines (NOTE: the “bottom-middle”-plot only shows 20 from 21 samples due to presentation purposes).</p>
Full article ">Figure 5
<p>Bivariate plots of <span class="html-italic">in situ</span> LAI<sub>2200</sub> and DVI derived from RapidEye imagery. The linear regression model (R<sup>2</sup> = 0.71) is indicated by the solid line, whereas green points represent samples of the classes “C. mopane shrub land”and “open woodland” (n = 17). Note the different origins of axes in the figure.</p>
Full article ">Figure 6
<p>LAI maps of the study region: (<b>a</b>) High-resolution map of LAI<sub>model</sub> (5 × 5 m) based on the transfer function given in Equation (6). For cartographic reasons, no differentiation for pixels with a LAI<sub>model</sub> &gt; 1.3 was made. (<b>b</b>) Aggregated map of the LAI<sub>model</sub>. (<b>c</b>) 8-day mean MODIS LAI (MOD15A2) map (spatial resolution: 1 × 1 km). (<b>d</b>) Absolute difference between (<b>b</b>) and <b>(c)</b>, where positive values indicate the aggregated LAI<sub>model</sub> to exceed MOD15A2, and vice versa. (NOTE: For (<b>b</b>) and (<b>d</b>), urban areas, as classified by (<b>c</b>), were a priori excluded from processing).</p>
Full article ">Figure 7
<p>(<b>a</b>) Standard deviation of MODIS LAI (MOD15A2) and non-vegetated pixels (grey in <a href="#remotesensing-07-04834-f006" class="html-fig">Figure 6</a>c) separated into desert (black) and water (blue). <b>(b)</b> Monthly MODIS Burned Area (MCD45A1) product from September 2010.</p>
Full article ">
4348 KiB  
Article
Analysis of the Urban Heat Island Effect in Shijiazhuang, China Using Satellite and Airborne Data
by Kai Liu, Hongbo Su, Lifu Zhang, Hang Yang, Renhua Zhang and Xueke Li
Remote Sens. 2015, 7(4), 4804-4833; https://doi.org/10.3390/rs70404804 - 20 Apr 2015
Cited by 75 | Viewed by 11562
Abstract
The urban heat island (UHI) effect resulting from rapid urbanization generally has a negative impact on urban residents. Shijiazhuang, the capital of Hebei Province in China, was selected to assess surface thermal patterns and its correlation with Land Cover Types (LCTs). This study [...] Read more.
The urban heat island (UHI) effect resulting from rapid urbanization generally has a negative impact on urban residents. Shijiazhuang, the capital of Hebei Province in China, was selected to assess surface thermal patterns and its correlation with Land Cover Types (LCTs). This study was conducted using Landsat TM images on the mesoscale level and airborne hyperspectral thermal images on the microscale level. Land surface temperature (LST) was retrieved from four scenes of Landsat TM data in the summer days to analyze the thermal spatial patterns and intensity of surface UHI (SUHI). Surface thermal characteristics were further examined by relating LST to percentage of imperious surface area (ISA%) and four remote sensing indices (RSIs), the Normalized Difference Vegetation Index (NDVI), Universal Pattern Decomposition method (VIUPD), Normalized Difference Built-up Index (NDBI) and Biophysical Composition Index (BCI). On the other hand, fives scenes of airborne TASI (Thermal Airborne Spectrographic Imager sensor) images were utilized to describe more detailed urban thermal characteristics of the downtown of Shijiazhuang city. Our results show that an obvious surface heat island effect existed in the study area during summer days, with a SUHI intensity of 2–4 °C. The analyses reveal that ISA% can provide an additional metric for the study of SUHI, yet its association with LST is not straightforward and this should a focus in future work. It was also found that two physically based indices, VIUPD and BCI, have the potential to account for the variation in urban LST. The results concerning on TASI indicate that diversity of impervious surfaces (rooftops, concrete, and mixed asphalt) contribute most to the SUHI, among all of the land cover features. Moreover, the effect of impervious surfaces on LST is complicated, and the composition and arrangement of land cover features may play an important role in determining the magnitude and intensity of SUHI. Overall, the analysis of urban thermal signatures at two spatial scales complement each other and the use of airborne imagery data with higher spatial resolution is helpful in revealing more details for understanding urban thermal environments. Full article
Show Figures

Figure 1

Figure 1
<p>Study area in Shijiazhuang, China. (<b>a</b>) WordView-2 RGB color composite of the study area. Line 1 (Left) and 2 (Right) are two north–south direction profile transects. (<b>b</b>) Location of the study area and Landsat TM imagery.</p>
Full article ">Figure 2
<p>Overall process flow.</p>
Full article ">Figure 3
<p>LCTs patterns and percent impervious surface area distribution derived from Landsat TM image in the study area on 5 September 2006 and 15 August 2010, (<b>a</b>) LCTs patterns on 5 September 2006; (<b>b</b>) LCT patterns on 15 August 2010; (<b>c</b>) percent impervious surface area on 5 September 2006; (<b>d</b>) percent impervious surface area on 15 August 2010.</p>
Full article ">Figure 4
<p>LST distribution in the study area derived from Landsat TM image on (<b>a</b>) 5 September 2006; (<b>b</b>) 23 August 2007; (<b>c</b>) 12 August 2009; and (<b>d</b>) 15 August 2010.</p>
Full article ">Figure 5
<p>(<b>a</b>) Accuracy assessment of estimated ISA% on August 15, 2010; (<b>b</b>) Scatterplot of LST <span class="html-italic">vs.</span> ISA% on August 15, 2010; (<b>c</b>) Scatterplot of Mean LST <span class="html-italic">vs.</span> ISA% on 15 August 2010.</p>
Full article ">Figure 6
<p>Spectral indices in the study area derived from Landsat TM image on 15 August 2010, (<b>a</b>) Normalized Difference Vegetation Index (NDVI); (<b>b</b>) Universal Pattern Decomposition (VIUPD); (<b>c</b>) Normalized Difference Built-up Index (NDBI) and (<b>d</b>) Biophysical Composition Index (BCI).</p>
Full article ">Figure 7
<p>Scatterplot of LST <span class="html-italic">versus</span> the remote sensing indexes on 15 August 2010 (<b>a</b>) Scatterplot of LST <span class="html-italic">vs.</span> NDVI; (<b>b</b>) Scatterplot of mean LST <span class="html-italic">vs.</span> NDVI; (<b>c</b>) Scatterplot of LST <span class="html-italic">vs.</span> VIUPD; and (<b>d</b>) Scatterplot of mean LST <span class="html-italic">vs.</span> VIUPD.</p>
Full article ">Figure 8
<p>Scatterplot of LST <span class="html-italic">vs.</span> remote sensing indexes on 15 August 2010 (a) Scatterplot of LST <span class="html-italic">vs.</span> NDBI; (<b>b</b>) Scatterplot of mean LST <span class="html-italic">vs.</span> NDBI; (<b>c</b>) Scatterplot of LST <span class="html-italic">vs.</span> BCI; and (<b>d</b>) Scatterplot of mean LST <span class="html-italic">vs.</span> BCI.</p>
Full article ">Figure 9
<p>LCTs map derived from WV2 image and SASI image using SVM and watershed segmentation method.</p>
Full article ">Figure 10
<p>LST derived from TASI imagery on (<b>a</b>) the morning of 7 August 2010; (<b>b</b>) the noon of 7 August 2010; (<b>c</b>) the noon of 15 August 2010; (<b>d</b>) the evening of 25 July 2010; and (<b>e</b>) the evening of 27 July 2010. (Note that the white regions are those could not pass the data quality control.)</p>
Full article ">Figure 11
<p>Histogram of LST over all land covers on (<b>a</b>) the morning of 7 August 2010; (<b>b</b>) the noon of 7 August 2010; (<b>c</b>) the noon of 15 August 2010; (<b>d</b>) the evening of 25 July 2010; (<b>e</b>) the evening of 27 July 2010.</p>
Full article ">Figure 12
<p>North–south land surface temperature profiles derived from TASI LST images, (<b>a</b>) Line 1; (<b>b</b>) Line 2 (locations marked in <a href="#remotesensing-07-04804-f001" class="html-fig">Figure 1</a>).</p>
Full article ">Figure 13
<p>Comparison of the LST from TASI and Landsat TM on 15 August 2010 (<b>a</b>) Scatter plot of the LST; (<b>b</b>) the LST profile on two line transects (Lines marked in <a href="#remotesensing-07-04804-f001" class="html-fig">Figure 1</a>).</p>
Full article ">
37001 KiB  
Article
Polarimetric Calibration of CASMSAR P-Band Data Affected by Terrain Slopes Using a Dual-Band Data Fusion Technique
by Lu Liao, Jie Yang, Pingxiang Li and Fenfen Hua
Remote Sens. 2015, 7(4), 4784-4803; https://doi.org/10.3390/rs70404784 - 20 Apr 2015
Cited by 2 | Viewed by 6259
Abstract
For airborne synthetic aperture radar (SAR) polarimetric calibration (PolCAL) based on distributed targets, it is important to ensure the removal of both the polarimetric distortion and terrain slope effect. This paper proposes a new technique for PolCAL in mountainous areas, without the use [...] Read more.
For airborne synthetic aperture radar (SAR) polarimetric calibration (PolCAL) based on distributed targets, it is important to ensure the removal of both the polarimetric distortion and terrain slope effect. This paper proposes a new technique for PolCAL in mountainous areas, without the use of corner reflectors (CRs). The technique based on dual-band data fusion consists of two steps. First, the polarization orientation angle shift (POAS), as a priori asymmetry information, is derived from X-band interferometry and applied to P-band fully-polarimetric data. Second, the crosstalk and cross-polarization (cross-pol) channel imbalance are iteratively determined using the POAS after dual-band data fusion. The performance and feasibility of the technique was evaluated by CRs. It was demonstrated that the proposed technique is capable of deriving the distortion parameters and performs better than the methods presented in Quegan and Ainsworth et al. The signal-to-noise ratio (SNR) and pedestal height have been investigated in polarimetric signatures. The proposed technique is useful for PolCAL in mountainous areas and for monitoring systems without CRs in long-term operation. Full article
(This article belongs to the Special Issue Remote Sensing Dedicated to Geographical Conditions Monitoring)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) Antennas mounted under the aircraft body; (<b>b</b>) Corner reflector deployed in the calibration site.</p>
Full article ">Figure 2
<p>Schematic diagram of inducing <span class="html-italic">θ</span> in the geometry.</p>
Full article ">Figure 3
<p>Process flow of the proposed polarimetric calibration (PolCAL) technique. POAS, polarization orientation angle shift (POAS); CR, corner reflector.</p>
Full article ">Figure 4
<p>(<b>a</b>) Manual matching, with the left amplitude image acquired in the X-band HH-polarization, and the right amplitude acquired in the P-band HH-polarization; (<b>b</b>) Overlapped image of the X-band single-pass interferometric DEM image and the P-band HH amplitude image; (<b>c</b>) The POAS derived from the X-band interferometric DEM. The near range is at the top of image.</p>
Full article ">Figure 5
<p>Comparison of the four methods. (<b>a</b>–<b>b</b>) The amplitude and phase of α; (<b>c</b>–<b>f</b>) The amplitude of <span class="html-italic">u</span>, <span class="html-italic">v</span>, <span class="html-italic">w</span>, and <span class="html-italic">z</span>.</p>
Full article ">Figure 5 Cont.
<p>Comparison of the four methods. (<b>a</b>–<b>b</b>) The amplitude and phase of α; (<b>c</b>–<b>f</b>) The amplitude of <span class="html-italic">u</span>, <span class="html-italic">v</span>, <span class="html-italic">w</span>, and <span class="html-italic">z</span>.</p>
Full article ">Figure 6
<p>Polarimetric signatures of the four CRs in several cases. (<b>a</b>–<b>c</b>) The Co-pol signature of CR 2; (<b>d</b>–<b>f</b>) The Cross-pol signature of CR 2; (<b>g</b>–<b>i</b>) The Co-pol signature of CR 3; (<b>j</b>–<b>l</b>) The Cross-pol signature of CR 3; (<b>m</b>–<b>o</b>) The Co-pol signature of CR 5; (<b>p</b>–<b>r</b>) The Cross-pol signature of CR 5; (<b>s</b>–<b>u</b>) The Co-pol signature of CR 6; (<b>v</b>–<b>x</b>) The Cross-pol signature of CR 6. From left to right, the signatures represent the results of un-calibrated, CR method and the proposed technique, respectively.</p>
Full article ">Figure 6 Cont.
<p>Polarimetric signatures of the four CRs in several cases. (<b>a</b>–<b>c</b>) The Co-pol signature of CR 2; (<b>d</b>–<b>f</b>) The Cross-pol signature of CR 2; (<b>g</b>–<b>i</b>) The Co-pol signature of CR 3; (<b>j</b>–<b>l</b>) The Cross-pol signature of CR 3; (<b>m</b>–<b>o</b>) The Co-pol signature of CR 5; (<b>p</b>–<b>r</b>) The Cross-pol signature of CR 5; (<b>s</b>–<b>u</b>) The Co-pol signature of CR 6; (<b>v</b>–<b>x</b>) The Cross-pol signature of CR 6. From left to right, the signatures represent the results of un-calibrated, CR method and the proposed technique, respectively.</p>
Full article ">Figure 7
<p>(<b>a</b>–<b>f</b>) RGB images by Pauli decomposition. (<b>g</b>–<b>k</b>) Col-pol signatures of the five typical distributed targets.</p>
Full article ">Figure 8
<p>(<b>a</b>) Optical image from Google Earth. (<b>b</b>) X-band HH-polarization amplitude image. (<b>c</b>) P-band HH-polarization amplitude image. CR 5 is displayed inside the red circle; bright points indicate the CR.</p>
Full article ">
2317 KiB  
Article
Object-Based Approach for Multi-Scale Mangrove Composition Mapping Using Multi-Resolution Image Datasets
by Muhammad Kamal, Stuart Phinn and Kasper Johansen
Remote Sens. 2015, 7(4), 4753-4783; https://doi.org/10.3390/rs70404753 - 17 Apr 2015
Cited by 155 | Viewed by 15613
Abstract
Providing accurate maps of mangroves, where the spatial scales of the mapped features correspond to the ecological structures and processes, as opposed to pixel sizes and mapping approaches, is a major challenge for remote sensing. This study developed and evaluated an object-based approach [...] Read more.
Providing accurate maps of mangroves, where the spatial scales of the mapped features correspond to the ecological structures and processes, as opposed to pixel sizes and mapping approaches, is a major challenge for remote sensing. This study developed and evaluated an object-based approach to understand what types of mangrove information can be mapped using different image datasets (Landsat TM, ALOS AVNIR-2, WorldView-2, and LiDAR). We compared and contrasted the ability of these images to map five levels of mangrove features, including vegetation boundary, mangrove stands, mangrove zonations, individual tree crowns, and species communities. We used the Moreton Bay site in Australia as the primary site to develop the classification rule sets and Karimunjawa Island in Indonesia to test the applicability of the rule sets. The results demonstrated the effectiveness of a conceptual hierarchical model for mapping specific mangrove features at discrete spatial scales. However, the rule sets developed in this study require modification to map similar mangrove features at different locations or when using image data acquired by different sensors. Across the hierarchical levels, smaller object sizes (i.e., tree crowns) required more complex classification rule sets. Incorporation of contextual information (e.g., distance and elevation) increased the overall mapping accuracy at the mangrove stand level (from 85% to 94%) and mangrove zonation level (from 53% to 59%). We found that higher image spatial resolution, larger object size, and fewer land-cover classes result in higher mapping accuracies. This study highlights the potential of selected images and mapping techniques to map mangrove features, and provides guidance for how to do this effectively through multi-scale mangrove composition mapping. Full article
(This article belongs to the Special Issue Remote Sensing of Mangroves: Observation and Monitoring)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Conceptual temporal and spatial hierarchical organization of mangroves features identifiable from remotely-sensed images, and the required image pixel resolution for mapping the features. (Symbols are courtesy of the Integration and Application Network, University of Maryland Center for Environmental Science—ian.umces.edu/symbols/).</p>
Full article ">Figure 2
<p>Study sites showing the major land-cover types and the field transects across the mangrove zonations.</p>
Full article ">Figure 3
<p>Flowchart of the mapping process applied in this study.</p>
Full article ">Figure 4
<p>Image objects hierarchy for multi-scale mangrove mapping, objects relationships, and the levels of information at each hierarchy level.</p>
Full article ">Figure 5
<p>Spectral reflectance profiles of the major land cover types extracted from WV-2 image within the Moreton Bay site.</p>
Full article ">Figure 6
<p>Subset of Whyte Island maps, showing (<b>a</b>) WV-2 standard false color composite, (<b>b</b>) vegetation (V) and non-vegetation (NV) discrimination using FDI, (<b>c</b>) spectral-based mangroves (M) and non-mangroves (NM) separation, and (<b>d</b>) band combination image to enhance the mangrove zonations. Tree crown delineation process showing (<b>e</b>) color composite of PC1, PC2, PC1 (RGB), (<b>f</b>) masked canopy gaps (white), (<b>g</b>) tree canopy seeds (red) on top of NIR band, and (<b>h</b>) tree crown polygons derived from region growing.</p>
Full article ">Figure 7
<p>Example of mangroves and non-mangroves area-based accuracy assessment; (<b>a</b>) reference map, (<b>b</b>) classified map from WV-2 image, and (<b>c</b>) classes produced from area intersection process.</p>
Full article ">Figure 8
<p>Subsets of green band semi-variograms, with up to 50 m lag distance, showing the variation of vegetation structure at different sites: (<b>a</b>) Fisherman Island, (<b>b</b>) Whyte Island, (<b>c</b>) Boondall wetlands and (<b>d</b>) Karimunjawa Island. Coordinates represent the approximate center of each image.</p>
Full article ">Figure 9
<p>Example subsets of mangrove composition maps at Whyte Island, Moreton Bay, showing all hierarchy levels produced from different image sources: (<b>a</b>) level 1 TM, (<b>b</b>) level 1 AVNIR-2, (<b>c</b>) level 1 WV-2, (<b>d</b>) WV-2 image (RGB 753), (<b>e</b>) level 2 TM, (<b>f</b>) level 2 AVNIR-2, (<b>g</b>) level 2 WV-2, (<b>h</b>) level 2 WV-2+LiDAR, (<b>i</b>) level 3 AVNIR-2, (<b>j</b>) level 3 WV-2, (<b>k</b>) level 3 WV-2+LiDAR, (<b>l</b>) level 4 pan-sharpened WV-2, (<b>m</b>) level 4 pan-sharpened WV-2+LiDAR, (<b>n</b>) WV-2 PC1,2,1, (<b>o</b>) level 5 pan-sharpened WV-2, (<b>p</b>) level 5 pan-sharpened WV-2+LiDAR and (<b>q</b>) true color aerial photograph.</p>
Full article ">Figure 10
<p>Comparison of the overall accuracy result of (<b>a</b>) different images and (<b>b</b>) different levels, and (<b>c</b>) area of the produced maps for level 1, 2, and 3 in Whyte Island, Moreton Bay.</p>
Full article ">
30530 KiB  
Article
A Quantitative Inspection on Spatio-Temporal Variation of Remote Sensing-Based Estimates of Land Surface Evapotranspiration in South Asia
by Ainong Li, Wei Zhao and Wei Deng
Remote Sens. 2015, 7(4), 4726-4752; https://doi.org/10.3390/rs70404726 - 17 Apr 2015
Cited by 13 | Viewed by 8494
Abstract
Evapotranspiration (ET) plays a key role in water resource management. It is important to understand the ET spatio-temporal pattern of South Asia for understanding and anticipating serious water resource shortages. In this study, daily ET in 2008 was estimated over South Asia by [...] Read more.
Evapotranspiration (ET) plays a key role in water resource management. It is important to understand the ET spatio-temporal pattern of South Asia for understanding and anticipating serious water resource shortages. In this study, daily ET in 2008 was estimated over South Asia by using MODerate Resolution Imaging Spectroradiometer (MODIS) products combined with field observations and Global Land Data Assimilation System (GLDAS) product through Surface Energy Balance System (SEBS) model. Monthly ET data were calculated based on daily ET and evaluated by the GLDAS ET data. Good agreements were found between two datasets for winter months (October to February) with R2 from 0.5 to 0.7. Spatio-temporal analysis of ET was conducted. Ten specific sites with different land cover types at typical climate regions were selected to analyze the ET temporal change pattern, and the result indicated that the semi-arid or arid areas in the northwest had the lowest average daily ET (around 0.3 mm) with a big fluctuation in the monsoon season, while the sites in the Indo-Gangetic Plain and in southern India has bigger daily ET (more than 3 mm) due to a large water supplement. It is suggested that the monsoon climate has a large impact on ET spatio-temporal variation in the whole region. Full article
(This article belongs to the Special Issue Recent Advances in Thermal Infrared Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Core countries and the locations of available meteorological stations around South Asia.</p>
Full article ">Figure 2
<p>Comparison between the observed air temperature and the fitted air temperature before and after the modification on the calibration method (Equation (7)).</p>
Full article ">Figure 3
<p>(<b>a</b>) Evaporative Fraction (EF); (<b>b</b>) Daily net radiation and (<b>c</b>) Daily ET estimation results of South Asia on 2 March 2008.</p>
Full article ">Figure 3 Cont.
<p>(<b>a</b>) Evaporative Fraction (EF); (<b>b</b>) Daily net radiation and (<b>c</b>) Daily ET estimation results of South Asia on 2 March 2008.</p>
Full article ">Figure 4
<p>Monthly ET estimation results in South Asia in 2008.</p>
Full article ">Figure 5
<p>Scatter plot between estimated monthly ET and the NOAH GLDAS monthly product in 2008.</p>
Full article ">Figure 6
<p>Monthly ET difference between the estimation result and the GLDAS product in March and June 2008. (<b>a</b>) presents the estimated monthly ET, (<b>b</b>) presents the GLDAS monthly ET product, (<b>c</b>) shows their difference in March, and (<b>d</b>) shows the NDVI spatial distribution of day of year (DOY) 81, 2008. (<b>e</b>) presents the estimated monthly ET, (<b>f</b>) presents the GLDAS monthly ET product, (<b>g</b>) shows their difference in June, and (<b>h</b>) shows the NDVI spatial distribution of DOY 161, 2008.</p>
Full article ">Figure 7
<p>Total precipitation of (<b>a</b>) dry season (March to May) and (<b>b</b>) monsoon season (June to August) acquired from TRMM precipitation product.</p>
Full article ">Figure 8
<p>Land cover map of South Asia acquired from the MODIS land cover product and the locations of the selected ten sites.</p>
Full article ">Figure 9
<p>Temporal variation of daily ET (<b>cycle</b>) and NDVI (<b>green symbol line</b>) of the sampling sites A-J with different land cover types listed in <a href="#remotesensing-07-04726-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 9 Cont.
<p>Temporal variation of daily ET (<b>cycle</b>) and NDVI (<b>green symbol line</b>) of the sampling sites A-J with different land cover types listed in <a href="#remotesensing-07-04726-t001" class="html-table">Table 1</a>.</p>
Full article ">Figure 10
<p>Annual ET spatial distribution in 2008 in South Asia.</p>
Full article ">Figure 11
<p>Scatter plot between annual precipitation and annual ET in South Asia in 2008.</p>
Full article ">
11355 KiB  
Article
Vertical Height Errors in Digital Terrain Models Derived from Airborne Laser Scanner Data in a Boreal-Alpine Ecotone in Norway
by Erik Næsset
Remote Sens. 2015, 7(4), 4702-4725; https://doi.org/10.3390/rs70404702 - 17 Apr 2015
Cited by 17 | Viewed by 6241
Abstract
It has been suggested that airborne laser scanning (ALS) could be used for operational monitoring of vegetation changes in the alpine tree line caused by climate change. Because the vegetation is low in such tree-less areas close to the alpine zone, the accuracy [...] Read more.
It has been suggested that airborne laser scanning (ALS) could be used for operational monitoring of vegetation changes in the alpine tree line caused by climate change. Because the vegetation is low in such tree-less areas close to the alpine zone, the accuracy of the digital terrain model (DTM) becomes crucial for early detection of, e.g., pioneer trees representing an ongoing tree migration given that the height of the vegetation may be on the same order of magnitude as the DTM uncertainty. The goal of this study was to assess and exemplify the vertical height errors of DTMs derived from ALS data under varying flying altitudes and pulse repetition frequencies (PRF). Important effects in the analysis were local terrain form, terrain surface, ground vegetation height, and terrain slope, because they may be correlated with recruitment patterns of pioneer trees. Based on 426 ground control points collected in a boreal-alpine ecotone, a standard deviation of 0.07–0.08 m was found for the lowest flying altitudes and lowest PRFs. For the highest PRF the standard deviation was 0.13 m. There were statistically significant mean errors for the different terrain forms and ground vegetation heights (?0.11 to 0.13 m). Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) Location of the study area in southern Norway (black square); (<b>b</b>) Design of the trial. The predefined sample points are marked as circles (<span class="html-italic">n</span> = 40) and those that were identified in field (<span class="html-italic">n</span> = 39) and terrain control points (<span class="html-italic">n</span> = 387) are shown as black dots. The green area is defined as forest according to the official N50 topographic map series. Accordingly, the light area is above the tree line. The black triangle shows the location of the national reference point of the Norwegian Mapping Authority (60°00′21.10929″N 9°01′43.97760″E, 951.904 m above sea level). Contour interval is 5 m.</p>
Full article ">Figure 2
<p>Examples of three ground control points with different characteristics. (<b>a</b>) A flat location with green vegetation surface and vegetation height 10–20 cm; (<b>b</b>) A flat location with rock/bare surface and vegetation height &lt; 10 cm; (<b>c</b>) A concave location with green vegetation surface and vegetation height &gt; 20 cm.</p>
Full article ">Figure 3
<p>Photographs taken from a central part of the study area towards south (<b>a</b>) and north (<b>b</b>) showing the stage of vegetation development as observed in the field at 11:00 UTC on 10 June 2007. The ALS data were acquired at 15:20–15:55 UTC on 11 June 2007. Apart from a few wet spots with some dead grass from the previous year still visible, the ground vegetation appeared fully developed.</p>
Full article ">Figure 4
<p>Illustration by backscatter intensity of the unexpected automatic shift in pulse repetition frequency (PRF) from 166 to 142 kHz for ACQ<sub>3</sub>. (<b>a</b>) Backscatter intensity for the two strips of ACQ<sub>3</sub> over the study area. Low backscatter intensity (dark grey) is shown for the western strip flown at 166 kHz, as planned. High backscatter intensity (light grey) is seen for the eastern strip (142 kHz); (<b>b</b>) Illustration of where the initial PRF (166 kHz) was restored in the eastern strip, approximately 1.2 km outside the study area, by a significant drop in backscatter intensity.</p>
Full article ">Figure 5
<p>Normal Q-Q plot for the distribution of the individual errors (ε) between TIN surface elevation and ground control elevation for the four different TIN models.</p>
Full article ">Figure 5 Cont.
<p>Normal Q-Q plot for the distribution of the individual errors (ε) between TIN surface elevation and ground control elevation for the four different TIN models.</p>
Full article ">
3067 KiB  
Article
Building Deformation Assessment by Means of Persistent Scatterer Interferometry Analysis on a Landslide-Affected Area: The Volterra (Italy) Case Study
by Silvia Bianchini, Fabio Pratesi, Teresa Nolesini and Nicola Casagli
Remote Sens. 2015, 7(4), 4678-4701; https://doi.org/10.3390/rs70404678 - 17 Apr 2015
Cited by 97 | Viewed by 9229
Abstract
In recent years, space-borne InSAR (interferometric synthetic aperture radar) techniques have shown their capabilities to provide precise measurements of Earth surface displacements for monitoring natural processes. Landslides threaten human lives and structures, especially in urbanized areas, where the density of elements at risk [...] Read more.
In recent years, space-borne InSAR (interferometric synthetic aperture radar) techniques have shown their capabilities to provide precise measurements of Earth surface displacements for monitoring natural processes. Landslides threaten human lives and structures, especially in urbanized areas, where the density of elements at risk sensitive to ground movements is high. The methodology described in this paper aims at detecting terrain motions and building deformations at the local scale, by means of satellite radar data combined with in situ validation campaigns. The proposed approach consists of deriving maximum settlement directions of the investigated buildings from displacement data revealed by radar measurements and then in the cross-comparison of these values with background geological data, constructive features and on-field evidence. This validation permits better understanding whether or not the detected movements correspond to visible and effective damages to buildings. The method has been applied to the southwestern sector of Volterra (Tuscany region, Italy), which is a landslide-affected and partially urbanized area, through the use of COSMO-SkyMed satellite images as input data. Moreover, we discuss issues and possible misinterpretations when dealing with PSI (Persistent Scatterer Interferometry) data referring to single manufactures and the consequent difficulty of attributing the motion rate to ground displacements, rather than to structural failures. Full article
Show Figures

Figure 1

Figure 1
<p>Methodology flowchart.</p>
Full article ">Figure 2
<p>Step by step procedure for building damage characterization and estimation: (<b>A</b>) PSI data and derived IDW surface displayed as bilinear interpolation layered on building boundaries; (<b>B</b>) IDW surface displayed as a nearest-neighbor discrete raster, with pixel centroids obtained for each raster cell; (<b>C</b>) schematic representation of differential settlement parameters used within the analysis; (<b>D</b>) computation of differential settlement direction within an appropriate buffer around the building.</p>
Full article ">Figure 3
<p>Volterra study area: geological map and section. In order to emphasize the morphology, a vertical exaggeration of 2× is applied to the section. Modified from [<a href="#B34-remotesensing-07-04678" class="html-bibr">34</a>].</p>
Full article ">Figure 4
<p>Velocity rates and spatial distribution of available PSI data overlapped on the landslide inventory map of Volterra area provided by the Tuscany region. The close-up study area is the black-contoured sector: (<b>A</b>) PSI COSMO-SkyMed in descending orbit; (<b>B</b>) PSI COSMO-SkyMed in ascending orbit; (<b>C</b>) sketch representing the two COSMO-SkyMed acquisition geometries combined with the local topography of the study area.</p>
Full article ">Figure 5
<p>(<b>A</b>) Location of the five case studies; (<b>B</b>) geological map of the study area; and (<b>C</b>) two longitudinal sections.</p>
Full article ">Figure 6
<p>Analysis results of building B1: building information, background data, building deformations, settlement parameters, main crack pattern, some photos of the field survey and a PS time series (the gap during year 2011 is due to missing acquisitions).</p>
Full article ">Figure 7
<p>Analysis results of building B3: building information, background data, building deformations, settlement parameters, main crack pattern and photos of field survey.</p>
Full article ">Figure 8
<p>Analysis results of building B3: building information, background data, building deformations, settlement parameters, main crack pattern and photos of field survey.</p>
Full article ">Figure 9
<p>Analysis results of building B4: building information, background data, building deformations, settlement parameters, main crack pattern and photos of the field survey.</p>
Full article ">Figure 10
<p>Analysis results of building B5: building information, background data, building deformations, settlement parameters, main crack pattern, some photos of the field survey and a PS time series (the gap during year 2011 is due to missing acquisitions).</p>
Full article ">
10183 KiB  
Article
Optimal Combination of Classification Algorithms and Feature Ranking Methods for Object-Based Classification of Submeter Resolution Z/I-Imaging DMC Imagery
by Fulgencio Cánovas-García and Francisco Alonso-Sarría
Remote Sens. 2015, 7(4), 4651-4677; https://doi.org/10.3390/rs70404651 - 17 Apr 2015
Cited by 38 | Viewed by 6944
Abstract
Object-based image analysis allows several different features to be calculated for the resulting objects. However, a large number of features means longer computing times and might even result in a loss of classification accuracy. In this study, we use four feature ranking methods [...] Read more.
Object-based image analysis allows several different features to be calculated for the resulting objects. However, a large number of features means longer computing times and might even result in a loss of classification accuracy. In this study, we use four feature ranking methods (maximum correlation, average correlation, Jeffries–Matusita distance and mean decrease in the Gini index) and five classification algorithms (linear discriminant analysis, naive Bayes, weighted k-nearest neighbors, support vector machines and random forest). The objective is to discover the optimal algorithm and feature subset to maximize accuracy when classifying a set of 1,076,937 objects, produced by the prior segmentation of a 0.45-m resolution multispectral image, with 356 features calculated on each object. The study area is both large (9070 ha) and diverse, which increases the possibility to generalize the results. The mean decrease in the Gini index was found to be the feature ranking method that provided highest accuracy for all of the classification algorithms. In addition, support vector machines and random forest obtained the highest accuracy in the classification, both using their default parameters. This is a useful result that could be taken into account in the processing of high-resolution images in large and diverse areas to obtain a land cover classification. Full article
Show Figures


<p>Work flow chart of the methodology.</p>
Full article ">
<p>Irrigation Unit 28. National Topographic Map 1:200,000 from the Spanish National Geographic Institute (IGN).</p>
Full article ">
<p>Kappa indices obtained with three classification algorithms: random forest, weighted k (wk)-NN and SVM and the five ranking methods.</p>
Full article ">
<p>Kappa indices obtained with two classification algorithms: naive Bayes and linear discriminant analysis and the five ranking methods.</p>
Full article ">
<p>Example of the land cover maps obtained with the five classification methods using the optimal feature set obtained with the method based on the Gini index. (<b>a</b>) RF-Gini i., 78 features; (<b>b</b>) wk-NN-Gini i., 46 features; (<b>c</b>) SVM-Gini i., 33 features; (<b>d</b>) nB-Gini i., 11 features; (<b>e</b>) LDA-Gini i., 65 features; and (<b>f</b>) Z/I-Imaging DMC image with 45-cm spatial resolution.</p>
Full article ">
1312 KiB  
Article
Remote Estimation of Leaf and Canopy Water Content in Winter Wheat with Different Vertical Distribution of Water-Related Properties
by Shishi Liu, Yi Peng, Wei Du, Yuan Le and Lu Li
Remote Sens. 2015, 7(4), 4626-4650; https://doi.org/10.3390/rs70404626 - 17 Apr 2015
Cited by 40 | Viewed by 7359
Abstract
This study analyzed the vertical distribution of gravimetric water content (GWC), relative water content (RWC), and equivalent water thickness (EWT) in winter wheat during heading and early ripening stages, and evaluated the position of leaf number at which Vegetation Indexes (VIs) can best [...] Read more.
This study analyzed the vertical distribution of gravimetric water content (GWC), relative water content (RWC), and equivalent water thickness (EWT) in winter wheat during heading and early ripening stages, and evaluated the position of leaf number at which Vegetation Indexes (VIs) can best retrieve canopy water-related properties of winter wheat. Results demonstrated that the vertical distribution of these properties followed a near-bell-shaped curve with the highest values at the intermediate leaf position. GWC of the top three or four leaves during the heading stage and the top two or three leaves during the early ripening stage can represent the GWC of the whole canopy, but the RWC and EWT of the whole canopy should be calculated based on the top four leaves. At leaf level, the analysis demonstrated strong relationships between EWT and VIs for the top leaf layer, but for GWCD, GWCF, and RWC, the strongest relationships with VIs were found in the intermediate leaf layers. At canopy level, VIs provided the most accurate estimation of GWCfor the top three or four leaves. Water absorption-based VIs could estimate canopy EWT of winter wheat for the top four leaves, but the suitable bands sensitive to water absorptions should be carefully selected for the studied species. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The illustration of the sampled wheat with leaves numerically labeled from the bottom (leaf 1) to the top (leaf 5) of the canopy.</p>
Full article ">Figure 2
<p>Leaf GWC<sub>D</sub> (<b>a</b>), GWC<sub>F</sub> (<b>b</b>), RWC (<b>c</b>), and EWT (<b>d</b>) averaged on each measurement date <span class="html-italic">versus</span> leaf number from the top (No.5 refers to the top leaf) to the bottom (No.1 refers to the bottom leaf) of the winter wheat canopy.</p>
Full article ">Figure 2 Cont.
<p>Leaf GWC<sub>D</sub> (<b>a</b>), GWC<sub>F</sub> (<b>b</b>), RWC (<b>c</b>), and EWT (<b>d</b>) averaged on each measurement date <span class="html-italic">versus</span> leaf number from the top (No.5 refers to the top leaf) to the bottom (No.1 refers to the bottom leaf) of the winter wheat canopy.</p>
Full article ">Figure 3
<p>Correlations between reflectance at all wavelengths and GWC<sub>D</sub>, GWC<sub>F</sub>, RWC, EWT at leaf level of winter wheat. The significant correlation (<span class="html-italic">p</span> &lt; 0.001) was indicated by <span class="html-italic">r</span> greater than 0.45 or less than −0.45 (140 samples).</p>
Full article ">Figure 4
<p>Relationships between NDII<sub>1510</sub> and GWC<sub>D</sub>, GWC<sub>F</sub>, RWC, EWT at leaf level of winter wheat.</p>
Full article ">Figure 5
<p>GWC<sub>D</sub> (<b>a</b>), GWC<sub>F</sub> (<b>b</b>), RWC (<b>c</b>), and EWT (<b>d</b>) for the cumulative leaf number within a winter wheat canopy. Cumulative leaf number 1 corresponds to the top leaf (leaf No. 5), and cumulative leaf number 5 corresponds to the top five leaves, from the top leaf (leaf No. 5) to the bottom leaf (leaf No. 1) of the canopy.</p>
Full article ">Figure 6
<p>Correlation between canopy reflectance at all wavelengths and GWC<sub>D</sub>, GWC<sub>F</sub>, RWC, and EWT for the cumulative leaf number. Top one leaf corresponds to leaf No. 5, and top five leaves correspond to the leaf No.5 to leaf No.1. The significant correlation (<span class="html-italic">p</span> &lt; 0.001) was indicated by <span class="html-italic">r</span> greater than 0.52 or less than −0.52 (33 samples).</p>
Full article ">Figure 7
<p>R<sup>2</sup> for the relationships between NDII<sub>2020</sub> and water-related properties for the cumulative leaf number.</p>
Full article ">Figure 8
<p>Mean RMSE between the measured canopy spectra and the canopy spectra simulated with PROSPECT+SAILH model for the top one, two, three, four, and five leaves, respectively.</p>
Full article ">
13766 KiB  
Article
Regional Leaf Area Index Retrieval Based on Remote Sensing: The Role of Radiative Transfer Model Selection
by Gaofei Yin, Jing Li, Qinhuo Liu, Weiliang Fan, Baodong Xu, Yelu Zeng and Jing Zhao
Remote Sens. 2015, 7(4), 4604-4625; https://doi.org/10.3390/rs70404604 - 16 Apr 2015
Cited by 47 | Viewed by 7764
Abstract
Physically-based approaches for estimating Leaf Area Index (LAI) using remote sensing data rely on radiative transfer (RT) models. Currently, many RT models are freely available, but determining the appropriate RT model for LAI retrieval is still problematic. This study aims to evaluate the [...] Read more.
Physically-based approaches for estimating Leaf Area Index (LAI) using remote sensing data rely on radiative transfer (RT) models. Currently, many RT models are freely available, but determining the appropriate RT model for LAI retrieval is still problematic. This study aims to evaluate the necessity of RT model selection for LAI retrieval and to propose a retrieval methodology using different RT models for different vegetation types. Both actual experimental observations and RT model simulations were used to conduct the evaluation. Each of them includes needleleaf forests and croplands, which have contrasting structural attributes. The scattering from arbitrarily inclined leaves (SAIL) model and the four-scale model, which are 1D and 3D RT models, respectively, were used to simulate the synthetic test datasets. The experimental test dataset was established through two field campaigns conducted in the Heihe River Basin. The results show that the realistic representation of canopy structure in RT models is very important for LAI retrieval. If an unsuitable RT model is used, then the root mean squared error (RMSE) will increase from 0.43 to 0.60 in croplands and from 0.52 to 0.63 in forests. In addition, an RT model’s potential to retrieve LAI is limited by the availability of a priori information on RT model parameters. 3D RT models require more a priori information, which makes them have poorer generalization capability than 1D models. Therefore, physically-based retrieval algorithms should embed more than one RT model to account for the availability of a priori information and variations in structural attributes among different vegetation types. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Land cover map of the upper and middle reaches of the Heihe River Basin in China (<b>a</b>) and the locations of the arid region experimental area (AREA) (<b>b</b>) and the forest experimental area (FEA) (<b>c</b>) sub-regions. The land cover map was downloaded from [<a href="#B55-remotesensing-07-04604" class="html-bibr">55</a>]. (b,c) False color images of ASTER/ETM+ acquired on 10 July 2012 and 28 May 2008, respectively. The black boxes and yellow crosses in (b) and (c) display the MODIS 500-m grids and field sampling plots.</p>
Full article ">Figure 2
<p>The regression lines between VI and LAI in AREA (<b>a</b>) and FEA (<b>b</b>). The dots indicate the measurements described in <a href="#sec3dot2dot1-remotesensing-07-04604" class="html-sec">Section 3.2.1</a> and were used to fit the regression lines. The triangles indicate the independent validation samples collected in different field campaigns in the same study area.</p>
Full article ">Figure 3
<p>RMSE and R<sup>2</sup> between the estimated and actual LAI values over the SAIL (<b>a</b>,<b>b</b>) and four-scale (<b>c</b>,<b>d</b>) synthetic test datasets as a function of measurement uncertainty. The SAIL and four-scale synthetic test datasets were simulated by the SAIL and four-scale models, respectively. The four bars for a specific uncertainty level represent the original SAIL and four-scaleNN without noise (SAILNN0 and four-scaleNN0), and the SAIL and four-scaleNN trained using 10% Gaussian white noisy training database (SAILNN10 and four-scaleNN10). Measurement uncertainty represents the corresponding level of Gaussian noise added to the synthetic test dataset.</p>
Full article ">Figure 4
<p>Scatter plots of the estimated and actual values of LAI. Using four-scaleNN to invert the SAIL test dataset (<b>a</b>) and using SAILNN to invert the four-scale test dataset (<b>b</b>). The dashed lines are the regression lines between the estimated and actual LAI.</p>
Full article ">Figure 5
<p>Comparison between SAIL LAI (<b>a</b>) and four-scale LAI (<b>b</b>) in croplands and SAIL LAI (<b>c</b>) and four-scale LAI (<b>d</b>) in forests with the corresponding reference LAI. The solid lines represent the 1:1 lines, and the dotted lines represent the accuracy boundaries (max (0.5, 20%)) specified by the Global Climate Observation System (GCOS).</p>
Full article ">Figure 6
<p>Distribution of the LAI values for each LAI map in AREA (<b>a</b>) and FEA (<b>b</b>). The solid, dashed and dash-doted vertical bars identify the locations of the mean values for the reference, SAIL and four-scale LAI, respectively.</p>
Full article ">
3424 KiB  
Article
Analysis of Geometric Primitives in Quantitative Structure Models of Tree Stems
by Åkerblom Markku, Pasi Raumonen, Mikko Kaasalainen and Eric Casella
Remote Sens. 2015, 7(4), 4581-4603; https://doi.org/10.3390/rs70404581 - 16 Apr 2015
Cited by 68 | Viewed by 10988
Abstract
One way to model a tree is to use a collection of geometric primitives to represent the surface and topology of the stem and branches of a tree. The circular cylinder is often used as the geometric primitive, but it is not the [...] Read more.
One way to model a tree is to use a collection of geometric primitives to represent the surface and topology of the stem and branches of a tree. The circular cylinder is often used as the geometric primitive, but it is not the only possible choice. We investigate various geometric primitives and modelling schemes, discuss their properties and give practical estimates for expected modelling errors associated with the primitives. We find that the circular cylinder is the most robust primitive in the sense of a well-bounded volumetric modelling error, even with noise and gaps in the data. Its use does not cause errors significantly larger than those with more complex primitives, while the latter are much more sensitive to data quality. However, in some cases, a hybrid approach with more complex primitives for the stem is useful. Full article
Show Figures


<p>Geometric primitives. From left to right: circular cylinder <b>(circyl)</b>, elliptic cylinder <b>(ellcyl)</b>, polygon cylinder <b>(polcyl)</b>, truncated circular cone <b>(cone)</b> and polyhedron <b>(trian)</b>. <b>(Top)</b> Perspective side view; <b>(bottom)</b> orthographic top view.</p>
Full article ">
<p>Stages of triangulation. <b>(Left)</b> Each vertex (green) is computed as the mean of points in each cell; <b>(centre)</b> neighbouring cells are connected by edges; <b>(right)</b> The quads are further triangulated, forming the polyhedron.</p>
Full article ">
<p>Curved pipe model used in the sensitivity analysis.</p>
Full article ">
<p>Relative volume error as a function of axis point estimate perturbation. The magnitude of the transition was received by multiplying the radius estimate by the scaling factor shown on the horizontal axis.</p>
Full article ">
<p>Relative volume error as a function of axis direction estimate perturbation.</p>
Full article ">
<p>Flattened portraits of generated stem models (vertical dimension scaled down to 20%).</p>
Full article ">
<p>Taper error curves. Taper curve for the generated Stem 7 (right axis) and the diameter error (left axis) in the taper curves of the reconstructed models.</p>
Full article ">
<p>Cross-sections of the models of Stem 8 along the vertical axis, reconstructed from the noisy measurements.</p>
Full article ">
<p>Effect of vertex count on <b>polcyl</b> and <b>trian</b>.</p>
Full article ">
32217 KiB  
Article
Potential of Space-Borne Hyperspectral Data for Biomass Quantification in an Arid Environment: Advantages and Limitations
by Harald Zandler, Alexander Brenning and Cyrus Samimi
Remote Sens. 2015, 7(4), 4565-4580; https://doi.org/10.3390/rs70404565 - 15 Apr 2015
Cited by 19 | Viewed by 8080
Abstract
In spite of considerable efforts to monitor global vegetation, biomass quantification in drylands is still a major challenge due to low spectral resolution and considerable background effects. Hence, this study examines the potential of the space-borne hyperspectral Hyperion sensor compared to the multispectral [...] Read more.
In spite of considerable efforts to monitor global vegetation, biomass quantification in drylands is still a major challenge due to low spectral resolution and considerable background effects. Hence, this study examines the potential of the space-borne hyperspectral Hyperion sensor compared to the multispectral Landsat OLI sensor in predicting dwarf shrub biomass in an arid region characterized by challenging conditions for satellite-based analysis: The Eastern Pamirs of Tajikistan. We calculated vegetation indices for all available wavelengths of both sensors, correlated these indices with field-mapped biomass while considering the multiple comparison problem, and assessed the predictive performance of single-variable linear models constructed with data from each of the sensors. Results showed an increased performance of the hyperspectral sensor and the particular suitability of indices capturing the short-wave infrared spectral region in dwarf shrub biomass prediction. Performance was considerably poorer in the area with less vegetation cover. Furthermore, spatial transferability of vegetation indices was not feasible in this region, underlining the importance of repeated model building. This study indicates that upcoming space-borne hyperspectral sensors increase the performance of biomass prediction in the world’s arid environments. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Overview of the research area, analyzed satellite images and field sites. The two Hyperion scenes were acquired on 3 August 2012 (western scene) and on 29 July 2013 (eastern scene), the Landsat OLI scene on 28 July 2013 respectively. DEM source: METI &amp; NASA [<a href="#B20-remotesensing-07-04565" class="html-bibr">20</a>].</p>
Full article ">Figure 2
<p>Photographs of (<b>a</b>) dwarf shrub stand located within the eastern Hyperion scene taken in fall 2014, and (<b>b</b>) dwarf shrub stand located within the western Hyperion scene with azonal grass vegetation in the background taken in summer 2013.</p>
Full article ">Figure 3
<p>Boxplots showing dwarf shrub biomass amounts of sites located in Hyperion scenes of August 2012 and July 2013, respectively. Each scene contains 30 field sites.</p>
Full article ">Figure 4
<p>Absolute values of Pearson’s correlation coefficients R of biomass with indices from field sites of feature sets (<b>a</b>) H2012, (<b>b</b>) LS2013a, (<b>c</b>) H2013, and (<b>d</b>) LS2013b. Black lines mark significant values controlled at a FDR &lt; 5%.</p>
Full article ">Figure 5
<p>Occurrence of the 10 most frequently chosen NDIs according to stepwise variable selection from all folds and repetitions (1000 selections) with feature sets (<b>a</b>) H2012, (<b>b</b>) H2012, (<b>c</b>) LS2013a and (<b>d</b>) LS2013b.</p>
Full article ">Figure 5 Cont.
<p>Occurrence of the 10 most frequently chosen NDIs according to stepwise variable selection from all folds and repetitions (1000 selections) with feature sets (<b>a</b>) H2012, (<b>b</b>) H2012, (<b>c</b>) LS2013a and (<b>d</b>) LS2013b.</p>
Full article ">
30947 KiB  
Article
Geo-Positioning Accuracy Using Multiple-Satellite Images: IKONOS, QuickBird, and KOMPSAT-2 Stereo Images
by Jaehoon Jeong, Chansu Yang and Taejung Kim
Remote Sens. 2015, 7(4), 4549-4564; https://doi.org/10.3390/rs70404549 - 15 Apr 2015
Cited by 36 | Viewed by 12062
Abstract
This paper investigates the positioning accuracy of image pairs achieved by integrating images from multiple satellites. High-resolution satellite images from IKONOS, QuickBird, and KOMPSAT-2 for Daejeon, Korea were combined to produce pairs of stereo images. From single-satellite stereo pairs to multiple-satellite image pairs, [...] Read more.
This paper investigates the positioning accuracy of image pairs achieved by integrating images from multiple satellites. High-resolution satellite images from IKONOS, QuickBird, and KOMPSAT-2 for Daejeon, Korea were combined to produce pairs of stereo images. From single-satellite stereo pairs to multiple-satellite image pairs, all available combinations were analyzed via a rational function model (RFM). The positioning accuracy of multiple-satellite pairs was compared to a typical single-satellite stereo pair. The results show that dual-satellite integration can be an effective alternative to single-satellite stereo imagery for horizontal position mapping, but is less accurate for vertical mapping. The integration of additional higher-resolution images can improve the overall accuracy of the existing two images, but, conversely, may result in lower accuracy when very weak convergence or bisector elevation (BIE) angles occur. This highlights that the use of higher resolution images may not ensure improved accuracy, as it can result in very weak geometry. The findings confirm that multiple-satellite images can replace or enhance typical stereo pairs, but also suggest the need for careful verification, including consideration of various geometric elements and image resolution. This paper reveals the potential, limitations, and important considerations for mapping applications using images from multiple satellites. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Representation of convergence and BIE angles on stereo geometry.</p>
Full article ">Figure 2
<p>Image coverage and GCPs’ distribution.</p>
Full article ">Figure 3
<p>Test cases for comparison between single-satellite stereo and dual-satellite integration. (<b>a</b>) single-satellite stereo, (<b>b</b>) dual-satellite stereo.</p>
Full article ">Figure 4
<p>Test cases for comparison before and after the integration of higher resolution single image. (<b>a</b>) single-satellite stereo, (<b>b</b>) single-satellite stereo + higher resolution single image.</p>
Full article ">Figure 5
<p>Test cases for comparison before and after the integration of higher resolution stereo image. (<b>a</b>) single-satellite stereo, (<b>b</b>) single-satellite stereo + higher resolution stereo image.</p>
Full article ">Figure 6
<p>Relationship between pointing error vectors of individual image and 3D position errors.</p>
Full article ">Figure 7
<p>Pointing error vectors and their subtracting vectors for single-satellite stereo and dual-satellite stereo (five ICPs). (<b>a</b>) KOMPSAT-2/IKONOS, (<b>b</b>) KOMPSAT-2/QuickBird, (<b>c</b>) IKONOS/QuickBird.</p>
Full article ">
6832 KiB  
Article
A Comparison of Two Approaches for Estimating the Wheat Nitrogen Nutrition Index Using Remote Sensing
by Pengfei Chen
Remote Sens. 2015, 7(4), 4527-4548; https://doi.org/10.3390/rs70404527 - 15 Apr 2015
Cited by 51 | Viewed by 7656
Abstract
Remote predictions of the nitrogen nutrition index (NNI) are useful for precise nitrogen (N) management in the field. Several studies have recommended two methods for estimating the NNI, which are classified as mechanistic and semi-empirical methods in this study. However, [...] Read more.
Remote predictions of the nitrogen nutrition index (NNI) are useful for precise nitrogen (N) management in the field. Several studies have recommended two methods for estimating the NNI, which are classified as mechanistic and semi-empirical methods in this study. However, no studies have been conducted to thoroughly analyze and compare these two methods. Using winter wheat as an example, this study compared the performances of these two methods for estimating the NNI to determine which method is more suitable for practical use. Field measurements were conducted to determine the above ground biomass, N concentration and canopy spectra during different wheat growth stages in 2012. Nearly 120 samples of data were collected and divided into different calibration and validation datasets (containing data from single or multi-growth stages). Based on the above datasets, the performances of the two NNI estimation methods were compared, and the influences of phenology on the methods were analyzed. All models that used the mechanistic method with different calibration datasets performed well when validated by validation datasets containing single growth or multi-growth stage data. The validation results had R2 values between 0.82 and 0.94, root mean square error (RMSE) values between 0.05 and 0.17, and RMSE% values between 5.10% and 14.41%. Phenology had no effect on this type of NNI estimation method. However, the semi-empirical method was influenced by phenology. The performances of the models established using this method were determined by the type of data used for calibration. Thus, the mechanistic method is recommended as a better method for estimating the NNI. By combining proper N management strategies, it can be used for precise N management. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Correlation matrices between (<b>a</b>) plant nitrogen concentration (%) and SPAD values; (<b>b</b>) biomass and LAI; (<b>c</b>) <span class="html-italic">NNI</span> and SPAD values; and (<b>d</b>) <span class="html-italic">NNI</span> and LAI.</p>
Full article ">Figure 2
<p>Scatterplots for (<b>a</b>) MCARI/MTVI2 <span class="html-italic">vs.</span> N concentration and (<b>b</b>) RTVI <span class="html-italic">vs.</span> biomass.</p>
Full article ">Figure 3
<p>Calibration results of the mechanistic NNI estimation method when using three calibration datasets.</p>
Full article ">Figure 4
<p>Scatterplots for (<b>a</b>) <span class="html-italic">REIP-LI vs. NNI</span> and (<b>b</b>) <span class="html-italic">MTCI vs. NNI</span>.</p>
Full article ">
14386 KiB  
Article
Assessing Field Spectroscopy Metadata Quality
by Barbara A. Rasaiah, Simon. D. Jones, Chris Bellman, Tim J. Malthus and Andreas Hueni
Remote Sens. 2015, 7(4), 4499-4526; https://doi.org/10.3390/rs70404499 - 15 Apr 2015
Cited by 11 | Viewed by 6359
Abstract
This paper presents the proposed criteria for measuring the quality and completeness of field spectroscopy metadata in a spectral archive. Definitions for metadata quality and completeness for field spectroscopy datasets are introduced. Unique methods for measuring quality and completeness of metadata to meet [...] Read more.
This paper presents the proposed criteria for measuring the quality and completeness of field spectroscopy metadata in a spectral archive. Definitions for metadata quality and completeness for field spectroscopy datasets are introduced. Unique methods for measuring quality and completeness of metadata to meet the requirements of field spectroscopy datasets are presented. Field spectroscopy metadata quality can be defined in terms of (but is not limited to) logical consistency, lineage, semantic and syntactic error rates, compliance with a quality standard, quality assurance by a recognized authority, and reputational authority of the data owners/data creators. Two spectral libraries are examined as case studies of operationalized metadata policies, and the degree to which they are aligned with the needs of field spectroscopy scientists. The case studies reveal that the metadata in publicly available spectral datasets are underperforming on the quality and completeness measures. This paper is part two in a series examining the issues central to a metadata standard for field spectroscopy datasets. Full article
(This article belongs to the Special Issue Field Spectroscopy and Radiometry)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Mean Z-scores for completeness by database user.</p>
Full article ">Figure 2
<p>Mean Z-scores for completeness by institute.</p>
Full article ">Figure 3
<p>Cumulative entropy for non-vegetation and mixed groups.</p>
Full article ">
6960 KiB  
Article
Forest Fire Smoke Detection Using Back-Propagation Neural Network Based on MODIS Data
by Xiaolian Li, Weiguo Song, Liping Lian and Xiaoge Wei
Remote Sens. 2015, 7(4), 4473-4498; https://doi.org/10.3390/rs70404473 - 15 Apr 2015
Cited by 71 | Viewed by 13816
Abstract
Satellite remote sensing provides global observations of the Earth’s surface and provides useful information for monitoring smoke plumes emitted from forest fires. The aim of this study is to automatically separate smoke plumes from the background by analyzing the MODIS data. An identification [...] Read more.
Satellite remote sensing provides global observations of the Earth’s surface and provides useful information for monitoring smoke plumes emitted from forest fires. The aim of this study is to automatically separate smoke plumes from the background by analyzing the MODIS data. An identification algorithm was improved based on the spectral analysis among the smoke, cloud and underlying surface. In order to get satisfactory results, a multi-threshold method is used for extracting training sample sets to train back-propagation neural network (BPNN) classification for merging the smoke detection algorithm. The MODIS data from three forest fires were used to develop the algorithm and get parameter values. These fires occurred in (i) China on 16 October 2004, (ii) Northeast Asia on 29 April 2009 and (iii) Russia on 29 July 2010 in different seasons. Then, the data from four other fires were used to validate the algorithm. Results indicated that the algorithm captured both thick smoke and thin dispersed smoke over land, as well as the mixed pixels of smoke over the ocean. These results could provide valuable information concerning forest fire location, fire spreading and so on. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Flowchart of the proposed algorithm for smoke detection.</p>
Full article ">Figure 2
<p>True-color composition RGB images of three cases used to extract seasonal samples: (<b>a</b>) Smoke plumes emitted from fire happened in Daxing’anling area, China in autumn, (<b>b</b>) Smoke plumes emitted from major fires in northeastern Asia in spring, (<b>c</b>) Smoke plumes emitted from fires in Russia in summer.</p>
Full article ">Figure 3
<p>Samples used for spectral analysis are extracted in this area during the forest fire happened on 28 June 2010, (<b>a</b>) True-color composition RGB image acquired over Daxing’anling area, China on 28 June 2010, (<b>b</b>) The extracted smoke samples are shown in the bright white area, (<b>c</b>) The extracted cloud samples are marked with bright white color.</p>
Full article ">Figure 4
<p>Response curves of the four cover types.</p>
Full article ">Figure 5
<p>Normalized distances between smoke and cloud, smoke and water, smoke and vegetation in reflectance of MODIS bands 1–8.</p>
Full article ">Figure 6
<p>3D scatter plot between smoke and cloud in <span class="html-italic">BT11</span>, <span class="html-italic">BTD (3.7–12)</span> and <span class="html-italic">R26</span>.</p>
Full article ">Figure 7
<p>The structure of BP neural networks that designed for identifying smoke plumes. S means smoke, C represents cloud whereas W/V is underlying surface.</p>
Full article ">Figure 8
<p>The evolution of BPNN performance with epoch.</p>
Full article ">Figure 9
<p>The decrease of accuracy and the increase of error by comparing the tests in different ways (one is that the training samples and testing samples extracted in the same season; the other is that the training samples extracted in summer while the testing samples extracted in spring).</p>
Full article ">Figure 10
<p>(<b>a</b>) True-color composition RGB image acquired over Daxing’anling area, China on 16 October 2004, (<b>b</b>) True-color composition RGB image acquired over Daxing’anling area, China and Amur region, Russia on 29 April 2009, (<b>c</b>) True-color composition RGB image acquired over Ryazan region, Russia on 29 July 2010, Panels (<b>d</b>) and (<b>e</b>) are the results of panel (a) and the rectangle area in (a) by using the algorithm, Panels (<b>f</b>) and (<b>g</b>) are the results of panel (b) and the rectangle area in (b), Panels (<b>h</b>) and (<b>i</b>) are the results of panel (c) and the rectangle area in (c).</p>
Full article ">Figure 11
<p>Smoke plumes detection by using the improved algorithm in Russia on 3 August 2012 (summer): (<b>a</b>) True-color composition RGB image of MODIS bands 1, 4 and 3 covering the detected area, (<b>b</b>) The rectangle area shown in panel (<b>a</b>), (<b>c</b>) The detected result of panel (<b>a</b>) and panel (<b>d</b>) is the result of rectangle area. The smoke plumes are depicted in red color.</p>
Full article ">Figure 12
<p>Smoke plumes detection in western Quebec, Canada on 19 June 2013 (spring): (<b>a</b>) True-color composition RGB image of MODIS bands 1, 4 and 3 covering the detected area, (<b>b</b>) The rectangle area shown in panel (a), (<b>c</b>) The detected result of panel (a) and panel (<b>d</b>) is the result of rectangle area, the smoke plumes are depicted in red color.</p>
Full article ">Figure 13
<p>Smoke plumes detection by using the algorithm in Greece on 24 August 2007 (summer): (<b>a</b>) True-color composition RGB image of MODIS bands 1, 4 and 3 covering the detected area, (<b>b</b>) The rectangle area shown in panel (a), (<b>c</b>) The detected result of panel (a) and panel (<b>d</b>) is the result of rectangle area, the smoke plumes are depicted in red color.</p>
Full article ">Figure 14
<p>Smoke plumes detection by using the algorithm around the Alice Spring (Australia) on 30 September 2011 (spring): (<b>a</b>) True-color composition RGB image of MODIS bands 1, 4 and 3 covering the detected area, (<b>b</b>) The rectangle area shown in panel (<b>a</b>), (<b>c</b>) The detected result of panel (<b>a</b>) and panel (<b>d</b>) is the result of rectangle area, the smoke plumes are depicted in red color.</p>
Full article ">Figure 15
<p>Smoke plumes detection by using the multi-threshold method in different locations. Panel (<b>a</b>) and Panel (<b>b</b>) are the results of two groups of thresholds in Russia on 3 August 2012. Panel (<b>c</b>) and Panel (<b>d</b>) are the results of two groups of thresholds in Quebec, Canada on 19 June 2013. Panel (<b>e</b>) and Panel (<b>f</b>) are the results of two groups of thresholds in Greece on 24 August 2007.</p>
Full article ">
8049 KiB  
Article
L-Band SAR Backscatter Related to Forest Cover, Height and Aboveground Biomass at Multiple Spatial Scales across Denmark
by Neha P. Joshi, Edward T. A. Mitchard, Johannes Schumacher, Vivian K. Johannsen, Sassan Saatchi and Rasmus Fensholt
Remote Sens. 2015, 7(4), 4442-4472; https://doi.org/10.3390/rs70404442 - 14 Apr 2015
Cited by 51 | Viewed by 9808
Abstract
Mapping forest aboveground biomass (AGB) using satellite data is an important task, particularly for reporting of carbon stocks and changes under climate change legislation. It is known that AGB can be mapped using synthetic aperture radar (SAR), but relationships between AGB and radar [...] Read more.
Mapping forest aboveground biomass (AGB) using satellite data is an important task, particularly for reporting of carbon stocks and changes under climate change legislation. It is known that AGB can be mapped using synthetic aperture radar (SAR), but relationships between AGB and radar backscatter may be confounded by variations in biophysical forest structure (density, height or cover fraction) and differences in the resolution of satellite and ground data. Here, we attempt to quantify the effect of these factors by relating L-band ALOS PALSAR HV backscatter and unique country-wide LiDAR-derived maps of vegetation penetrability, height and AGB over Denmark at different spatial scales (50 m to 500 m). Trends in the relations indicate that, first, AGB retrieval accuracy from SAR improves most in mapping at 100-m scale instead of 50 m, and improvements are negligible beyond 250 m. Relative errors (bias and root mean squared error) decrease particularly for high AGB values (\(>\)110 Mg ha\(^{-1}\)) at coarse scales, and hence, coarse-scale mapping (\(\ge\)150 m) may be most suited for areas with high AGB. Second, SAR backscatter and a LiDAR-derived measure of fractional forest cover were found to have a strong linear relation (R\(^2\) = 0.79 at 250-m scale). In areas of high fractional forest cover, there is a slight decline in backscatter as AGB increases, indicating signal attenuation. The two results demonstrate that accounting for spatial scale and variations in forest structure, such as cover fraction, will greatly benefit establishing adequate plot-sizes for SAR calibration and the accuracy of derived AGB maps. Full article
Show Figures


<p>(<b>a</b>) Forested areas in Denmark (map provided by Geodatastyrelsen), locations of the species trial plots (STPs) in clusters and the extent of the ALOS PALSAR scenes and LiDAR strips that were removed before analysis. (<b>b</b>,<b>c</b>,<b>d</b>) An example of a forest plantation in Denmark with extracted LiDAR variables, vegetation interception ratio (VIR), mean height (MH) and AGB<sub>L</sub>, at 50 m<span class="html-italic">×</span> 50 m pixel resolution. Built areas and water bodies were removed from all maps prior to analysis.</p>
Full article ">
<p>Example of a cluster of STPs, overlaid on a high-resolution (0.16 m <span class="html-italic">×</span> 0.16 m) airborne optical image obtained during leaf-on season in 2010 by COWI. Each STP contains a single tree species and is placed within the center of a forest stand.</p>
Full article ">
<p>A hypothetical demonstration of the variation of vegetation interception ratio (VIR) with vegetation distribution and aboveground biomass (AGB). (<b>a</b>) Pixels with different AGB density (due to varying tree height), but similar VIR. (<b>b</b>) Pixels with similar AGB density, but varying VIR due to the different distribution of stems.</p>
Full article ">
<p>(<b>a</b>) Results of the LiDAR-to-AGB<sub>STP</sub> model (<a href="#FD2" class="html-disp-formula">Equation (2)</a>) predicting AGB<sub>STP</sub> with the 95th height percentile multiplied by VIR. (<b>b</b>) Normal probability plot of residual distribution. Normality was assessed using the Shapiro–Wilk test (<span class="html-italic">p</span> = 0.17). (<b>c</b>) Predicted <span class="html-italic">vs</span>. measured biomass for STPs of different sizes demonstrate the lack of dependence of bias on plot size.</p>
Full article ">
<p>Distribution of AGB<sub>L</sub> for forested areas at different spatial scales (50 m to 500 m, indicated on each sub–figure). Frequency (<span class="html-italic">y</span>-axis) refers to the number of observations in each range of AGB (<span class="html-italic">x</span>-axis).</p>
Full article ">
<p>SAR-to-VIR<sub>L</sub> model in the power domain (<b>Left panel</b>) and dB (<b>Right panel</b>) for the training dataset at a few of the tested spatial scales in this study (50 m to 500 m, indicated on each sub–figure). Data are represented on a 2D histogram density plot, with values byte-scaled (0 to 255) to the color bars shown.</p>
Full article ">
<p>SAR-to-AGB<sub>L</sub> model (<a href="#FD4" class="html-disp-formula">Equation (4)</a>) (<b>Left panel</b>) and the performance of the inverted model (<a href="#FD5" class="html-disp-formula">Equation (5)</a>) (<b>Right panel</b>) used to retrieve AGB at multiple spatial scales (50 m to 500 m, indicated on each sub–figure). Observations (N) are represented on a 2D histogram density plot, with values byte-scaled (0 to 255). High density observations that run across the residual plots at ~130 Mg ha<sup>−1</sup> show where a maximum retrieval value (MRV) was assigned to predicted AGB when <math display="inline"> <mrow> <msubsup> <mi>σ</mi> <mrow> <mi>H</mi> <mi>V</mi></mrow> <mn>0</mn></msubsup></mrow></math> fell above the predictable range.</p>
Full article ">
<p>SAR-to-MH<sub>L</sub> model at multiple spatial scales (50 m to 500 m, indicated on each sub–figure). Observations (N) are represented on a 2D histogram density plot, with values byte-scaled (0 to 255).</p>
Full article ">
<p>AGB retrieval accuracy ((<b>Left panel</b>), RMSE; (<b>Centre panel</b>), relative RMSE; (<b>Right panel</b>), bias) from <math display="inline"> <mrow> <msubsup> <mi>σ</mi> <mrow> <mi>H</mi> <mi>V</mi></mrow> <mn>0</mn></msubsup></mrow></math> as a function of spatial scale before and after applying a maximum retrieval value (MRV) to predicted AGB.</p>
Full article ">
3007 KiB  
Article
Advancing of Land Surface Temperature Retrieval Using Extreme Learning Machine and Spatio-Temporal Adaptive Data Fusion Algorithm
by Yang Bai, Man Sing Wong, Wen-Zhong Shi, Li-Xin Wu and Kai Qin
Remote Sens. 2015, 7(4), 4424-4441; https://doi.org/10.3390/rs70404424 - 14 Apr 2015
Cited by 47 | Viewed by 9832
Abstract
As a critical variable to characterize the biophysical processes in ecological environment, and as a key indicator in the surface energy balance, evapotranspiration and urban heat islands, Land Surface Temperature (LST) retrieved from Thermal Infra-Red (TIR) images at both high temporal and spatial [...] Read more.
As a critical variable to characterize the biophysical processes in ecological environment, and as a key indicator in the surface energy balance, evapotranspiration and urban heat islands, Land Surface Temperature (LST) retrieved from Thermal Infra-Red (TIR) images at both high temporal and spatial resolution is in urgent need. However, due to the limitations of the existing satellite sensors, there is no earth observation which can obtain TIR at detailed spatial- and temporal-resolution simultaneously. Thus, several attempts of image fusion by blending the TIR data from high temporal resolution sensor with data from high spatial resolution sensor have been studied. This paper presents a novel data fusion method by integrating image fusion and spatio-temporal fusion techniques, for deriving LST datasets at 30 m spatial resolution from daily MODIS image and Landsat ETM+ images. The Landsat ETM+ TIR data were firstly enhanced based on extreme learning machine (ELM) algorithm using neural network regression model, from 60 m to 30 m resolution. Then, the MODIS LST and enhanced Landsat ETM+ TIR data were fused by Spatio-temporal Adaptive Data Fusion Algorithm for Temperature mapping (SADFAT) in order to derive high resolution synthetic data. The synthetic images were evaluated for both testing and simulated satellite images. The average difference (AD) and absolute average difference (AAD) are smaller than 1.7 K, where the correlation coefficient (CC) and root-mean-square error (RMSE) are 0.755 and 1.824, respectively, showing that the proposed method enhances the spatial resolution of the predicted LST images and preserves the spectral information at the same time. Full article
(This article belongs to the Special Issue Recent Advances in Thermal Infrared Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Structure of the single hidden layer feed forward neural network using Extreme Learning Machine (ELM).</p>
Full article ">Figure 2
<p>Flowchart of the proposed fusion model for predicting synthetic LST image at 30 m resolution.</p>
Full article ">Figure 3
<p>Location of the study area.</p>
Full article ">Figure 4
<p>The flowchart of the testing and simulated experiment.</p>
Full article ">Figure 5
<p>Thermal spatial sharpening results of testing experiments via ELM algorithm. False colour images of downscaled Landsat ETM+ multispectral data at 60 m (<b>upper row</b>), downscaled TIR images at 120 m (<b>middle row</b>) and sharpened TIR images at 60 m (<b>lower row</b>). The original TIR images at 60 m refer to the middle row in <a href="#remotesensing-07-04424-f006" class="html-fig">Figure 6</a>. From left to right, they were acquired on 20 October, 7 December and 23 December 2013.</p>
Full article ">Figure 6
<p>Thermal spatial sharpening results of simulated experiments using ELM algorithm. False colour images of observed Landsat ETM+ multispectral data at 30 m (<b>upper row</b>), TIR images at 60 m (<b>middle row</b>) and Enhanced TIR images at 30 m (<b>lower row</b>). From left to right, they were acquired on 20 October, 7 December and 23 December 2013.</p>
Full article ">Figure 7
<p>SADFAT result in the testing experiment. <b>(a)</b> MODIS LST images resampled from 1 km to 60 m; <b>(b)</b> SADFAT-derived image at 60 m; <b>(c)</b> original Landsat ETM+ LST image at 60 m spatial on 7 December 2013.</p>
Full article ">Figure 8
<p>SADFAT result in simulated experiment. <b>(a)</b> MODIS LST images resampled from 1 km to 30 m; <b>(b)</b> SADFAT-derived image at 30 m; <b>(c)</b> up-scaled Landsat ETM+ LST image at 30 m spatial on 7 December 2013.</p>
Full article ">Figure 9
<p>Scatter plots between the predicted and original LSTs at 60 m of the testing experiment on 7 December 2013. <span class="html-italic">x</span>-axis denotes the original LSTs, and <span class="html-italic">y</span>-axis denotes the prediction LSTs.</p>
Full article ">
21160 KiB  
Article
Drought Variability and Land Degradation in Semiarid Regions: Assessment Using Remote Sensing Data and Drought Indices (1982–2011)
by Sergio M. Vicente-Serrano, Daniel Cabello, Miquel Tomás-Burguera, Natalia Martín-Hernández, Santiago Beguería, Cesar Azorin-Molina and Ahmed El Kenawy
Remote Sens. 2015, 7(4), 4391-4423; https://doi.org/10.3390/rs70404391 - 14 Apr 2015
Cited by 122 | Viewed by 14146
Abstract
We analyzed potential land degradation processes in semiarid regions worldwide using long time series of remote sensing images and the Normalized Difference Vegetation Index (NDVI) for the period 1981 to 2011. The objectives of the study were to identify semiarid regions showing a [...] Read more.
We analyzed potential land degradation processes in semiarid regions worldwide using long time series of remote sensing images and the Normalized Difference Vegetation Index (NDVI) for the period 1981 to 2011. The objectives of the study were to identify semiarid regions showing a marked decrease in potential vegetation activity, indicative of the occurrence of land degradation processes, and to assess the possible influence of the observed drought trends quantified using the Standardized Precipitation Evapotranspiration Index (SPEI). We found that the NDVI values recorded during the period of maximum vegetation activity (NDVImax) predominantly showed a positive evolution in the majority of the semiarid regions assessed, but NDVImax was highly correlated with drought variability, and the trends of drought events influenced trends in NDVImax at the global scale. The semiarid regions that showed most increase in NDVImax (the Sahel, northern Australia, South Africa) were characterized by a clear positive trend in the SPEI values, indicative of conditions of greater humidity and lesser drought conditions. While changes in drought severity may be an important driver of NDVI trends and land degradation processes in semiarid regions worldwide, drought did not apparently explain some of the observed changes in NDVImax. This reflects the complexity of vegetation activity processes in the world’s semiarid regions, and the difficulty of defining a universal response to drought in these regions, where a number of factors (natural and anthropogenic) may also affect on land degradation. Full article
(This article belongs to the Special Issue Remote Sensing of Land Degradation in Drylands)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Spatial distribution of the semiarid regions analyzed in the study.</p>
Full article ">Figure 2
<p>Month of the year in which the monthly average maximum Normalized Difference Vegetation Index (NDVI<sub>max</sub>) was recorded.</p>
Full article ">Figure 3
<p>Example of the analysis used to determine the influence of various drought time-scales on the maximum NDVI at 23.1°E, 28.7°S. The left column shows the evolution of the 3-, 6- and 12-month Standardized Precipitation Evapotranspiration Index (SPEI) and the monthly NDVI. The right column shows the evolution of the NDVI (solid line) and the 1-, 3- and 12-month SPEI (dotted line) during the month in which the average maximum NDVI was recorded. The bottom plot shows the correlation between the NDVI and the 1- to 12-month SPEI.</p>
Full article ">Figure 4
<p>Trends in the NDVI in semiarid regions during the month of average maximum vegetation activity, under the assumption of monotonic change (1982–2011). Top: sign and signification (<span class="html-italic">p</span> &lt; 0.1) of the trends. Bottom: magnitude of the trend (in NDVI units decade<sup>−1</sup>).</p>
Full article ">Figure 5
<p>Percentage of semiarid lands according to the observed NDVI magnitude of change (NDVI units decade<sup>−1</sup>) between 1982 and 2011. The inset graph shows the percentages of semiarid lands having negative and positive trends (including signification).</p>
Full article ">Figure 6
<p>Box plots showing the average NDVI, aridity and Water Use Efficiency (WUE) for the four categories of NDVI trend.</p>
Full article ">Figure 7
<p>Pearson’s r correlation coefficients between SPEI and NDVI time series in the month of maximum vegetation activity (1982–2011). Top: magnitude of Pearson’s <span class="html-italic">r</span> correlations. Bottom: areas (in black) having positive and significant correlations between the SPEI and NDVI (<span class="html-italic">p</span> &lt; 0.1).</p>
Full article ">Figure 8
<p>Percentage of semiarid regions according to the correlation coefficient between the NDVI and SPEI (1982–2011). The inset plot shows the percentages of semiarid lands having negative and positive trends (including signification).</p>
Full article ">Figure 9
<p>Box plots showing the average NDVI, aridity and WUE for the four categories of correlation between the SPEI and NDVI.</p>
Full article ">Figure 10
<p>Trends of the SPEI (1982–2011) corresponding to the month of maximum vegetation activity and the SPEI time scale at which the maximum SPEI/NDVI correlation was identified in each pixel. (<b>Top</b>) magnitude of the trends (in SPEI units decade<sup>−1</sup>). (<b>Bottom</b>) sign and signification (<span class="html-italic">p</span> &lt; 0.1) of the trends.</p>
Full article ">Figure 11
<p>Relationship between the magnitude of the NDVI and SPEI trends (units decade<sup>−1</sup>) between 1982 and 2011 for the world’s semiarid regions.</p>
Full article ">Figure 12
<p>Box plot showing the magnitude of change in the SPEI corresponding to the categories of NDVI trend between1982 and 2011.</p>
Full article ">Figure 13
<p>Examples of the evolution of residuals between NDVI observations and predictions, using the SPEI as predictor. Left: USA at −91°E, 31.2°N, Right: Algeria at 2.5°E, 35.5°N.</p>
Full article ">Figure 14
<p>Trend in residual series (difference between observed and predicted NDVI, according to variations of SPEI) in areas having positive and significant correlations between the SPEI and NDVI. Top: sign and signification (<span class="html-italic">p</span> &lt; 0.1) in the trend of the residual series. Bottom: magnitude of the residual series (in NDVI units decade<sup>−1</sup>).</p>
Full article ">Figure 15
<p>Percentage of semiarid regions in which there was a significant correlation between the SPEI and NDVI, according to the magnitude of change in the residual series. The inset plot shows the percentages of the analyzed lands having negative and positive trends of the residuals (including signification).</p>
Full article ">Figure 16
<p>Box plots showing the average NDVI, aridity, WUE and magnitude of change in the SPEI and NDVI in the four categories of trend in the residual series.</p>
Full article ">
27583 KiB  
Article
A Practical Split-Window Algorithm for Retrieving Land Surface Temperature from Landsat-8 Data and a Case Study of an Urban Area in China
by Meijun Jin, Junming Li, Caili Wang and Ruilan Shang
Remote Sens. 2015, 7(4), 4371-4390; https://doi.org/10.3390/rs70404371 - 14 Apr 2015
Cited by 69 | Viewed by 11505
Abstract
This paper proposes a practical split-window algorithm (SWA) for retrieving land surface temperature (LST) from Landsat-8 Thermal Infrared Sensor (TIRS) data. This SWA has a universal applicability and a set of parameters that can be applied when retrieving LSTs year-round. The atmospheric transmittance [...] Read more.
This paper proposes a practical split-window algorithm (SWA) for retrieving land surface temperature (LST) from Landsat-8 Thermal Infrared Sensor (TIRS) data. This SWA has a universal applicability and a set of parameters that can be applied when retrieving LSTs year-round. The atmospheric transmittance and the land surface emissivity (LSE), the essential SWA input parameters, of the Landsat-8 TIRS data are determined in this paper. We also analysed the error sensitivity of these SWA input parameters. The accuracy evaluation of the proposed SWA in this paper was conducted using the software MODTRAN 4.0. The root mean square error (RMSE) of the simulated LST using the mid-latitude summer atmospheric profile is 0.51 K, improving on the result of 0.93 K from Rozenstein (2014). Among the 90 simulated data points, the maximum absolute error is 0.99 °C, and the minimum absolute error is 0.02 °C. Under the Tropical model and 1976 US standard atmospheric conditions, the RMSE of the LST errors are 0.70 K and 0.63 K, respectively. The accuracy results indicate that the SWA provides an LST retrieval method that features not only high accuracy but also a certain universality. Additionally, the SWA was applied to retrieve the LST of an urban area using two Landsat-8 images. The SWA presented in this paper should promote the application of Landsat-8 data in the study of environmental evolution. Full article
(This article belongs to the Special Issue Recent Advances in Thermal Infrared Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The Planck radiance-temperature (180 K–363 K) curve of the Landsat-8 TIRS.</p>
Full article ">Figure 2
<p>(<b>A</b>) The LST evaluated error (K) in the case of undervaluing the atmospheric water vapour content by 0.1 g·cm<sup>−2</sup> for a water vapour content in the range of 1.0–4.0 g·cm<sup>−2</sup> and an interval of 1.0 g·cm<sup>−2</sup> over six gradations of T<sub>10</sub> = 283.0–333.0 K and a T<sub>10</sub>–T<sub>11</sub> value range of −3.0–3.0K. e<sub>10</sub> = 0.967, and e<sub>11</sub> = 0.971. (<b>B</b>) The LST evaluation error (K) for undervaluing the atmospheric water vapour content by 0.2 g·cm<sup>−2</sup> under the same conditions as in (A).</p>
Full article ">Figure 3
<p>The LST evaluation error (K) in the case of undervaluing the atmospheric water vapour content by 0.1 g·cm<sup>−2</sup> (<b>A</b>) and 0.2 g·cm<sup>−2</sup> (<b>B</b>) for water vapour contents in the range of 1.0–4.0 g·cm<sup>−2</sup> at an interval of 1.0 g·cm<sup>−2</sup> over six gradations of T<sub>10</sub> = 283.0–333.0 K and a T<sub>10</sub>–T<sub>11</sub> value of 1 K. The LSE variable is independent.</p>
Full article ">Figure 4
<p>The illustration of the LST calculation error (K) in the case of undervaluing LSE by 0.001 (<b>A</b>) and 0.005 (<b>B</b>) over four gradations of T<sub>10</sub> = 270.0–330.0 K with an interval of 20 K. The atmospheric water vapour content was constant at 1.5 g·cm<sup>−2</sup>, and the T10–T11 values were 1.0, 2.0, and 3.0 K. The LSE variable is independent.</p>
Full article ">Figure 5
<p>The illustration of the LST calculation error (K) in the case of undervaluing LSE by 0.001 (<b>A</b>); and 0.005 (<b>B</b>) over six gradations of T<sub>10</sub> = 283.0–333.0 K with an interval of 10.0 K. The atmospheric water vapour content varies in the range of 1.0–4.0 g·cm<sup>−2</sup> with an interval of 1.0 g·cm<sup>−2</sup>, and the T10–T11 value is 1.0 K. The LSE variable is independent.</p>
Full article ">Figure 6
<p>Workflow for retrieving LST from Landsat-8 data.</p>
Full article ">Figure 7
<p>The LST Spatial distribution of the study urban area retrieved from Landsat-8: (<b>A</b>) on June 27, 2013, and (<b>B</b>) on June 30, 2014.</p>
Full article ">
49581 KiB  
Article
Countrywide Stereo-Image Matching for Updating Digital Surface Models in the Framework of the Swiss National Forest Inventory
by Christian Ginzler and Martina L. Hobi
Remote Sens. 2015, 7(4), 4343-4370; https://doi.org/10.3390/rs70404343 - 13 Apr 2015
Cited by 124 | Viewed by 9959
Abstract
Surface models provide key knowledge of the 3-d structure of forests. Aerial stereo imagery acquired during routine mapping campaigns covering the whole of Switzerland (41,285 km2), offers a potential data source to calculate digital surface models (DSMs). We present an automated [...] Read more.
Surface models provide key knowledge of the 3-d structure of forests. Aerial stereo imagery acquired during routine mapping campaigns covering the whole of Switzerland (41,285 km2), offers a potential data source to calculate digital surface models (DSMs). We present an automated workflow to generate a nationwide DSM with a resolution of 1 × 1 m based on photogrammetric image matching. A canopy height model (CHM) is derived in combination with an existing digital terrain model (DTM). ADS40/ADS80 summer images from 2007 to 2012 were used for stereo matching, with ground sample distances (GSD) of 0.25 m in lowlands and 0.5 m in high mountain areas. Two different image matching strategies for DSM calculation were applied: one optimized for single features such as trees and for abrupt changes in elevation such as steep rocks, and another optimized for homogeneous areas such as meadows or glaciers. The country was divided into 165,500 blocks, which were matched independently using an automated workflow. The completeness of successfully matched points was high, 97.9%. To test the accuracy of the derived DSM, two reference data sets were used: (1) topographic survey points (n = 198) and (2) stereo measurements (n = 195,784) within the framework of the Swiss National Forest Inventory (NFI), in order to distinguish various land cover types. An overall median accuracy of 0.04 m with a normalized median absolute deviation (NMAD) of 0.32 m was found using the topographic survey points. The agreement between the stereo measurements and the values of the DSM revealed acceptable NMAD values between 1.76 and 3.94 m for forested areas. A good correlation (Pearson’s r = 0.83) was found between terrestrially measured tree height (n = 3109) and the height derived from the CHM. Optimized image matching strategies, an automatic workflow and acceptable computation time mean that the presented approach is suitable for operational usage at the nationwide extent. The CHM will be used to reduce estimation errors of different forest characteristics in the Swiss NFI and has high potential for change detection assessments, since an aerial stereo imagery update is available every six years. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Image data used: (<b>a</b>) Coverage of ADS40/ADS80 images acquired in 2007–2012, (<b>b</b>) flight lines for the image acquisition for the entire of Switzerland. The imagery was acquired with a GSD of 0.25 m at lower altitudes (<b>denser flight lines</b>) and a GSD of 0.50 m at higher altitudes in the mountainous regions (<b>less dense flight lines</b>).</p>
Full article ">Figure 2
<p>Image matching workflow. (<b>A</b>) Rasterizing into 0.5 × 0.5 km blocks; (<b>B</b>) selecting the stripe with the smallest distance from nadir to the centre of the block; and (<b>C</b>) selecting the stripe with the second smallest distance from nadir to the center of the block (1) The result after the matching process for the block using the 1st most nadir stereo stripe and the “ADS steep” image matching strategy. Red areas in the image could not be matched. The completeness for this block is 84%. (2) The result after applying the “ADS flat” matching strategy. At 91%, the completeness threshold is still not fulfilled. (3) The result after matching using the 2nd most nadir stripe and the “ADS steep” strategy (completeness 96%) and (4) the final DSM after “ADS flat” strategy with a matching completeness of 96%.</p>
Full article ">Figure 3
<p>Examples of the performance of the two different correlation strategies: “ADS steep”, optimized for steep areas and very heterogeneous surfaces such as forest canopies (<b>a</b>) input CIR image and (<b>b</b>) shaded relief of the resulting DSM and “ADS flat”, optimized for homogeneous areas with little contrast such as glaciers (<b>c</b>) input CIR image and (<b>d</b>) shaded relief of the resulting DSM.</p>
Full article ">Figure 4
<p>Reference data for accuracy assessment. The 198 topographic survey points of the national triangulation network were used for an independent error analysis.</p>
Full article ">Figure 5
<p>Shaded relief of the generated DSM with a colour gradient for canopy height (0–50 m) on the left and the used ADS80 input image data on the right: (<b>a</b>) buildings in Zurich, (<b>b</b>) forests in the Swiss lowlands and (<b>c</b>) glaciers in the Valais.</p>
Full article ">Figure 6
<p>Correlation between tree heights of the CHM and terrestrial tree height measurements. Total number of measurements is 3109. Tree height in the canopy height model refers to the maximum value within a buffer of 2.5 m radius around each individual tree. Tree height in the field was measured using a Vertex ultrasonic hypsometer.</p>
Full article ">Figure 7
<p>Comparison of canopy height from the CHM with terrestrially measured canopy height at the NFI sample plot level. One of the four stages of stand development was assigned to each sample plot. (<b>a</b>) Young growth refers to regeneration areas with a medium tree height of 8 m; (<b>b</b>) pole wood refers to stands with a high number of stems and a mean tree height of 8–20 m; (<b>c</b>) the timber stage is characterized by closed stands with a mean height &gt; 20 m; and (<b>d</b>) mixed refers to stands where there are different stages of stand development.</p>
Full article ">Figure 8
<p>Correlation of tree height extracted from the CHM with terrestrial tree height measurements distinguishing between (<b>a</b>) deciduous and (<b>b</b>) coniferous trees. The total number of measurements was 3109. Tree height in the canopy height model refers to the maximum value within a buffer of 5 m diameter around each individual tree. Tree height in the field was measured using a Vertex ultrasonic hypsometer.</p>
Full article ">Figure 9
<p>Correlation between tree height extracted from the CHM and terrestrial tree height measurements distinguishing between (<b>a</b>) lower and (<b>b</b>) higher elevations. Total number of measurements was 3109. Tree height in the canopy height model refers to the maximum value within a buffer of 5 m diameter around each individual tree. Tree height in the field was measured using a Vertex ultrasonic hypsometer.</p>
Full article ">Figure 10
<p>Comparison of canopy height of the CHM with terrestrially measured canopy height at the sample plot level of the NFI. One of the four species mixture classes (<b>a–d</b>) depending on the percentage of intermixed coniferous trees was assigned to each sample plot.</p>
Full article ">Figure 11
<p>Comparison of canopy height of the CHM with terrestrially measured canopy height at the sample plot level of the NFI. One of the four stand structure classes (<b>a–d</b>), defined by the layering and structure of the canopy was assigned to each sample plot.</p>
Full article ">Figure 12
<p>Sample plots of the Swiss National Forest Inventory spaced in a regular 1.4 km grid, consisting of 25 lattice points each where height and surface cover information was measured.</p>
Full article ">
3013 KiB  
Article
Automatic Case-Based Reasoning Approach for Landslide Detection: Integration of Object-Oriented Image Analysis and a Genetic Algorithm
by Jie Dou, Kuan-Tsung Chang, Shuisen Chen, Ali P. Yunus, Jin-King Liu, Huan Xia and Zhongfan Zhu
Remote Sens. 2015, 7(4), 4318-4342; https://doi.org/10.3390/rs70404318 - 13 Apr 2015
Cited by 134 | Viewed by 12668
Abstract
This paper proposes an automatic method for detecting landslides by using an integrated approach comprising object-oriented image analysis (OOIA), a genetic algorithm (GA), and a case-based reasoning (CBR) technique. It consists of three main phases: (1) image processing and multi-image segmentation; (2) feature [...] Read more.
This paper proposes an automatic method for detecting landslides by using an integrated approach comprising object-oriented image analysis (OOIA), a genetic algorithm (GA), and a case-based reasoning (CBR) technique. It consists of three main phases: (1) image processing and multi-image segmentation; (2) feature optimization; and (3) detecting landslides. The proposed approach was employed in a fast-growing urban region, the Pearl River Delta in South China. The results of detection were validated with the help of field surveys. The experimental results indicated that the proposed OOIA-GA-CBR (0.87) demonstrates higher classification performance than the stand-alone OOIA (0.75) method for detecting landslides. The area under curve (AUC) value was also higher than that of the simple OOIA, indicating the high efficiency of the proposed landslide detection approach. The case library created using the integrated model can be reused for time-independent analysis, thus rendering our approach superior in comparison to other traditional methods, such as the maximum likelihood classifier. The results of this study thus facilitate fast generation of accurate landslide inventory maps, which will eventually extend our understanding of the evolution of landscapes shaped by landslide processes. Full article
(This article belongs to the Special Issue Remote Sensing in Geology)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The spatial location of the study area overlaying the elevation.</p>
Full article ">Figure 2
<p>Examples of landslides in the study area. Principal landslide types observed in Conghua: (<b>a</b>) debris slide threatening a house; (<b>b</b>) soil creep; and (<b>c</b>) rock fall. Evidence of previous landslide scarps in the form of (<b>d</b>) terraced off, but covered with little, sharp boundary in the vegetation, and (<b>e</b>,<b>f</b>) concrete structures built for slope stabilization to arrest future failures.</p>
Full article ">Figure 3
<p>The integrated structure for automatic landslide detection, comprising three processes: (1) multi-segmentation by OOIA; (2) feature selection by GA; and (3) detection by CBR and validation using field work data.</p>
Full article ">Figure 4
<p>Details of QuickBird image segmentation including four scales: (<b>a</b>) scale = 150; (<b>b</b>) scale = 100; (<b>c</b>) scale = 50; (<b>d</b>) scale = 30.</p>
Full article ">Figure 5
<p>Object attributes exported in the OOIA analysis.</p>
Full article ">Figure 6
<p>Illustration of GA: (1) an example of the single point crossover for a binary GA, randomly setting the parents and obtaining the offspring; (2) the mutation for a binary GA, in which the bits are randomly chosen and the allele’s values are altered.</p>
Full article ">Figure 7
<p>Illustration of landslide detection using the CBR method in this study.</p>
Full article ">Figure 8
<p>Examples showing the preparation of reference landslides in the case library, in the form of Quickbird images: (<b>a</b>–<b>d</b>) are young landslides (bright area with no vegetation cover); (<b>e</b>,<b>f</b>) are old landslides (less bright areas covered with concrete structures and vegetation).</p>
Full article ">Figure 9
<p>A schmatic diagram showing the computation of true positives, false positives, false negative, and true negatives for verifying the accuracy of the model. (<b>a</b>) true positive: actual landslides that were correctly classified as landslides; (<b>b</b>) false positive: nonlandslides that were incorreclty classified as landslides; (<b>c</b>) false negative: landslides that were incorrectly classified as nonlandslides; and (<b>d</b>) true neagtive: nonlandslides that were correctly classified as nonlandslides.</p>
Full article ">Figure 10
<p>Relationship between feature number and fitness value in the GA process.</p>
Full article ">Figure 11
<p>Plots of old (<b>a</b>) and young (<b>b</b>) landslide curves. The curves have similar periodic trend changes.</p>
Full article ">Figure 12
<p>Results of landslide detection at different scales (150, 100, and 50).</p>
Full article ">Figure 13
<p>Data for validation: (<b>a</b>) track of GPS route for field survey; (<b>b</b>) examples of landslide detection points overlying the QuickBird image, using CBR.</p>
Full article ">Figure 14
<p>Prediction rate of ROC curve for each landslide class using CBR: (<b>a</b>) young landslides; (<b>b</b>) old landslides.</p>
Full article ">Figure 15
<p>Prediction rate of the ROC curve for each landslide class: (<b>a</b>) young landslides, using OOIA; (<b>b</b>) old landslides, using OOIA; (<b>c</b>) young landslides, using MLC; and (<b>d</b>) old landslides, using MLC.</p>
Full article ">Figure 16
<p>Examples of misclassification: (<b>a</b>) a bare quarry was misclassified as a young landslide, and (<b>b</b>) the related field photo; (<b>c</b>) a rock outcrop with sparse forest was misclassified as an old landslide, and (<b>d</b>) the related field photo. The tones of misclassification were similar, with a bright appearance of the landslide area in the satellite image.</p>
Full article ">
1946 KiB  
Article
Reviewing ALOS PALSAR Backscatter Observations for Stem Volume Retrieval in Swedish Forest
by Maurizio Santoro, Leif E.B. Eriksson and Johan E.S. Fransson
Remote Sens. 2015, 7(4), 4290-4317; https://doi.org/10.3390/rs70404290 - 13 Apr 2015
Cited by 35 | Viewed by 8141
Abstract
Between 2006 and 2011, the Advanced Land Observing Satellite (ALOS) Phased Array L-type Synthetic Aperture Radar (PALSAR) instrument acquired multi-temporal datasets under several environmental conditions and multiple configurations of look angle and polarization. The extensive archive of SAR backscatter observations over the forest [...] Read more.
Between 2006 and 2011, the Advanced Land Observing Satellite (ALOS) Phased Array L-type Synthetic Aperture Radar (PALSAR) instrument acquired multi-temporal datasets under several environmental conditions and multiple configurations of look angle and polarization. The extensive archive of SAR backscatter observations over the forest test sites of Krycklan (boreal) and Remningstorp (hemi-boreal), Sweden, was used to assess the retrieval of stem volume at stand level. The retrieval was based on the inversion of a simple Water Cloud Model with gaps; single estimates of stem volume are then combined to obtain the final multi-temporal estimate. The model matched the relationship between the SAR backscatter and the stem volume under all configurations. The retrieval relative Root Mean Square Error (RMSE) differed depending upon environmental conditions, polarization and look angle. Stem volume was best retrieved in Krycklan using only HV-polarized data acquired under unfrozen conditions with a look angle of 34.3° (relative RMSE: 44.0%). In Remningstorp, the smallest error was obtained using only HH-polarized data acquired under predominantly frozen conditions with a look angle of 34.3° (relative RMSE: 35.1%). The relative RMSE was below 30% for stands >20 ha, suggesting high accuracy of ALOS PALSAR estimates of stem volumes aggregated at moderate resolution. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Map of Sweden showing the location of the two test sites of Remningstorp and Krycklan.</p>
Full article ">Figure 2
<p>Bar chart of stem volume distribution in Remningstorp and Krycklan. Bars were grouped into intervals of 20 m<sup>3</sup>/ha.</p>
Full article ">Figure 3
<p>Panels (<b>a</b>) and (<b>b</b>) illustrate data from Krycklan. Panels (<b>c)</b> and <b>(d)</b> illustrate data from Remningstorp. SAR backscatter with respect to stem volume for unfrozen (Unfr.) conditions (panels (a) and (c)) and frozen (Fr.) conditions (panels (b) and (d)) with among the highest correlation coefficients (see <a href="#remotesensing-07-04290-t005" class="html-table">Table 5</a>). Look angle: 34.3°.</p>
Full article ">Figure 4
<p>Measured and modeled PALSAR backscatter as a function of stem volume for Krycklan. The model curves are based on Equation (1). The crosses and the vertical bars represent the median backscatter and the interquartile range in 25 m<sup>3</sup>/ha large intervals of stem volume. All data acquired under unfrozen conditions unless specified in the legend (fr = frozen).</p>
Full article ">Figure 5
<p>Estimates of the model parameter β at Krycklan with respect to environmental conditions for PALSAR data acquired with 34.3° look angle and HH-polarization. “<span class="html-italic">Unfr.</span>” refers to unfrozen conditions. Frozen conditions refer to images acquired under dry conditions as well as cases with snow fall. If precipitation was recorded at &lt;2 mm, the unfrozen conditions were moist; otherwise, the conditions were wet.</p>
Full article ">Figure 6
<p>Histograms of the estimates of the model parameter β at Krycklan as a function of look angle.</p>
Full article ">Figure 7
<p>Three examples of modeled backscatter as a function of stem volume assuming β unknown (solid curves) and set <span class="html-italic">a priori</span> (dashed curves). The measurements of backscatter are represented by crosses and vertical bars (median and interquartile range) for groups of stem volume, each being 25 m<sup>3</sup>/ha wide. Test site: Krycklan.</p>
Full article ">Figure 8
<p>Scatter plot of single-image RMSEs for a model with three unknowns (horizontal axis) and a model with two unknowns where the parameter β was set <span class="html-italic">a priori</span> equal to 0.006 (vertical axis). Test site: Krycklan.</p>
Full article ">Figure 9
<p>Scatter plot of single-image RMSEs for a model with three unknowns (horizontal axis) and a model with two unknowns where the parameter β was set <span class="html-italic">a priori</span> equal to 0.006 (vertical axis). Test site: Remningstorp.</p>
Full article ">Figure 10
<p>Distribution of single-image retrieval RMSE at Krycklan for combinations of look angle, polarization and environmental conditions for which multi-temporal SAR backscatter observations (at least three) were available.</p>
Full article ">Figure 11
<p>Scatter plot of retrieved stem volume with respect to <span class="html-italic">in situ</span> stem volume in the case of all HV-polarized images acquired during 2008 over Krycklan with a look angle of 34.3°.</p>
Full article ">Figure 12
<p>Scatter plot of retrieved stem volume with respect to <span class="html-italic">in situ</span> stem volume in the case of all HH-polarized images acquired during 2009 over Remningstorp with a look angle of 34.3°.</p>
Full article ">Figure 13
<p>Relative RMSE with respect to minimum stand size for the multi-temporal combination of stem volumes estimated from HH- and HV-polarized images acquired during 2008 over Krycklan with a look angle of 34.3°.</p>
Full article ">
Previous Issue
Next Issue
Back to TopTop