[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (9)

Search Parameters:
Keywords = ground canopy thermal imagery

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
25 pages, 12024 KiB  
Article
Linking High-Resolution UAV-Based Remote Sensing Data to Long-Term Vegetation Sampling—A Novel Workflow to Study Slow Ecotone Dynamics
by Fabian Döweler, Johan E. S. Fransson and Martin K.-F. Bader
Remote Sens. 2024, 16(5), 840; https://doi.org/10.3390/rs16050840 - 28 Feb 2024
Viewed by 1466
Abstract
Unravelling slow ecosystem migration patterns requires a fundamental understanding of the broad-scale climatic drivers, which are further modulated by fine-scale heterogeneities just outside established ecosystem boundaries. While modern Unoccupied Aerial Vehicle (UAV) remote sensing approaches enable us to monitor local scale ecotone dynamics [...] Read more.
Unravelling slow ecosystem migration patterns requires a fundamental understanding of the broad-scale climatic drivers, which are further modulated by fine-scale heterogeneities just outside established ecosystem boundaries. While modern Unoccupied Aerial Vehicle (UAV) remote sensing approaches enable us to monitor local scale ecotone dynamics in unprecedented detail, they are often underutilised as a temporal snapshot of the conditions on site. In this study in the Southern Alps of New Zealand, we demonstrate how the combination of multispectral and thermal data, as well as LiDAR data (2019), supplemented by three decades (1991–2021) of treeline transect data can add great value to field monitoring campaigns by putting seedling regeneration patterns at treeline into a spatially explicit context. Orthorectification and mosaicking of RGB and multispectral imagery produced spatially extensive maps of the subalpine area (~4 ha) with low spatial offset (Craigieburn: 6.14 ± 4.03 cm; Mt Faust: 5.11 ± 2.88 cm, mean ± standard error). The seven multispectral bands enabled a highly detailed delineation of six ground cover classes at treeline. Subalpine shrubs were detected with high accuracy (up to 90%), and a clear identification of the closed forest canopy (Fuscospora cliffortioides, >95%) was achieved. Two thermal imaging flights revealed the effect of existing vegetation classes on ground-level thermal conditions. UAV LiDAR data acquisition at the Craigieburn site allowed us to model vegetation height profiles for ~6000 previously classified objects and calculate annual fine-scale variation in the local solar radiation budget (20 cm resolution). At the heart of the proposed framework, an easy-to-use extrapolation procedure was used for the vegetation monitoring datasets with minimal georeferencing effort. The proposed method can satisfy the rapidly increasing demand for high spatiotemporal resolution mapping and shed further light on current treeline recruitment bottlenecks. This low-budget framework can readily be expanded to other ecotones, allowing us to gain further insights into slow ecotone dynamics in a drastically changing climate. Full article
(This article belongs to the Section Ecological Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Overview of the workflow proposed to monitor fine-scale changes at alpine ecosystem boundaries. The workflow is based on UAV surveys, collecting RGB (2 cm resolution), multispectral, thermal and LiDAR data. These datasets are then used to create ancillary datasets (vegetation classification, vegetation height, solar radiation and seedling position), and all results are then linked to an existing field transect, here the monitoring of <span class="html-italic">Fuscospora cliffortioides</span> recruitment.</p>
Full article ">Figure 2
<p>Comparative display of RGB and multispectral imagery for Mt Faust, New Zealand. This image presents a side-by-side comparison of 5 cm resolution multispectral imagery (displayed in false colour; <b>top</b> section) and 2 cm resolution RGB imagery (<b>bottom</b> section). The purpose is to assess the precision of initial outputs by identifying vegetation boundaries at the study site.</p>
Full article ">Figure 3
<p>Accuracy assessment (in terms of Root Mean Square Error, RMSE) of RGB vs. multispectral imagery for the Craigieburn site, New Zealand, by manually selecting GCP centres in ArcGIS Pro. Images (<b>a</b>,<b>c</b>) show the RGB imagery, while (<b>b</b>,<b>d</b>) depict the multispectral imagery with marked GCPs. Image (<b>e</b>) presents an overlay of both spectral layers, enabling the measurement of spatial displacement. The mean offset across all 14 GCP’s between RGB and multispectral imagery was 5.29 ± 2.34 cm.</p>
Full article ">Figure 4
<p>An RGB image and a point cloud generated from LiDAR data for the Craigieburn site, New Zealand. The image shows (<b>a</b>) an RGB image (2 cm resolution), (<b>b</b>) a vegetation height model (20 cm resolution), and (<b>c</b>) a highlighted vegetation height model for a selected area within the site. The vegetation height data for both topographic features and vegetation are obtained by computing the difference between the DTM and the DSM.</p>
Full article ">Figure 5
<p>A 2D representation of seedling transects in the Southern Alps, New Zealand (−43.064°S, 171.425°E, Craigieburn), originally established by P. Wardle in 1991. The positions of seedlings along (<span class="html-italic">X</span>-axis, purple lines) and perpendicular to (<span class="html-italic">Y</span>-axis, teal lines) the transect are derived from the GPS coordinates of marked transect poles (indicated by yellow points).</p>
Full article ">Figure 6
<p>Results of the vegetation classification for both subalpine treeline ecotones. The figure presents RGB images of (<b>a</b>) Craigieburn and (<b>b</b>) Mt Faust (2 cm resolution), as well as vegetation classification via seven spectral reflectance bands of (<b>c</b>) Craigieburn and (<b>d</b>) Mt Faust (5 cm resolution). Bare ground, unvegetated ground (scree), mountain beech (<span class="html-italic">Fuscospora cliffortioides</span>) and tussock grasses (<span class="html-italic">Chionochloa</span> spp.), and three types of alpine shrubs (<span class="html-italic">Dracophyllum uniflorum</span>, <span class="html-italic">Leucopogon colensoi</span>, <span class="html-italic">Podocarpus nivalis</span>) are clearly distinguished. At Mt Faust, scree areas were lacking.</p>
Full article ">Figure 7
<p>Calculated solar radiation budget (Wh/m<sup>2</sup>) for the Craigieburn site, New Zealand (−43.064°S, 171.425°E), based on a DSM (20 cm resolution) generated from LiDAR data (20 May 2019). (<b>a</b>) RGB visualisation (2 cm resolution) and (<b>b</b>) LiDAR generated point cloud visualisation of the site. (<b>c</b>) Predicted solar radiation in Wh/m<sup>2</sup> for an entire summer and (<b>d</b>) late autumn day using an (<b>e</b>) adjusted colour scale.</p>
Full article ">Figure 8
<p>Effect of DEM resolution in predicting solar radiation regimes. Solar radiation map based on the LiDAR DSM (<b>left</b>, 20 cm resolution) compared with the standard DTM (<b>right</b>, 8 m resolution; LINZ, 2016) source material for Craigieburn site in late austral autumn (20 May 2019).</p>
Full article ">Figure 9
<p>Results of the UAV thermal surveys at Craigieburn site, New Zealand (late austral autumn, 20 May 2019). Temperature recordings (°C) in the morning (<b>a</b>), at noon (<b>b</b>) and the differences between both surveys (<b>c</b>).</p>
Full article ">Figure 10
<p>Box-and-violin plots of the maximum surface temperature range during daylight hours (late autumn) by ground cover type. Maximum differences in surface temperature were derived from UAV-based thermal imaging surveys in the morning and at noon on east-facing slopes of Craigieburn site, New Zealand. Different lower-case letters indicate statistically significant differences at <span class="html-italic">α</span> = 0.05 (post hoc test based on Tukey contrasts).</p>
Full article ">Figure 11
<p>A 3D-visualisation of the Craigieburn site, New Zealand, with <span class="html-italic">Fuscospora cliffortioides</span> seedling position and height at the treeline. The illustration used UAV-based RGB imagery (2 cm resolution) over the DTM (subalpine zone) and the DSM. Seedling positions were aligned along a long-term vegetation monitoring transect. The symbol size is proportional to seedling height recordings in the summer of 2019.</p>
Full article ">Figure 12
<p>High-quality LiDAR DSM layer (20 cm resolution) used for calculation of shading potential of individual trees at the treeline (<b>a</b>,<b>b</b>) of the Craigieburn site, New Zealand. The procedure enabled (<b>c</b>) 3D visualisations of exposed <span class="html-italic">Fuscospora cliffortioides</span> trees (green cylinder), measured as part of a treeline monitoring programme. Solar radiation projections (Wh/m<sup>2</sup>) illustrate the shading potential of individual trees in (<b>d</b>) late austral autumn and (<b>e</b>) summer.</p>
Full article ">Figure 13
<p>A 3D-visualisation of the Craigieburn site, New Zealand, with tree seedling positions and heights (indicated by symbol size). (<b>a</b>,<b>c</b>) Thermal differences (late austral autumn, 20 May 2019, morning vs. noon) above canopy draped over the LiDAR DSM (20 cm resolution). Seedling positions from a treeline census dataset are located along the transect, and their heights (summer 2019) are displayed and compared with (<b>b</b>) calculated solar radiation (Wh/m<sup>2</sup>) for the same late austral autumn day.</p>
Full article ">
19 pages, 5157 KiB  
Article
Measuring Evapotranspiration Suppression from the Wind Drift and Spray Water Losses for LESA and MESA Sprinklers in a Center Pivot Irrigation System
by Behnaz Molaei, R. Troy Peters, Abhilash K. Chandel, Lav R. Khot, Claudio O. Stockle and Colin S. Campbell
Water 2023, 15(13), 2444; https://doi.org/10.3390/w15132444 - 2 Jul 2023
Viewed by 2573
Abstract
Wind drift and evaporation loss (WDEL) of mid-elevation spray application (MESA) and low-elevation spray application (LESA) sprinklers on a center pivot and linear-move irrigation machines are measured and reported to be about 20% and 3%, respectively. It is important to estimate the fraction [...] Read more.
Wind drift and evaporation loss (WDEL) of mid-elevation spray application (MESA) and low-elevation spray application (LESA) sprinklers on a center pivot and linear-move irrigation machines are measured and reported to be about 20% and 3%, respectively. It is important to estimate the fraction of WDEL that cools and humidifies the microclimate causing evapotranspiration (ET) suppression, mitigating the measured irrigation system losses. An experiment was conducted in 2018 and 2019 in a commercial spearmint field near Toppenish, Washington. The field was irrigated with an 8-span center pivot equipped with MESA but had three spans that were converted to LESA. All-in-one weather sensors (ATMOS-41) were installed just above the crop canopy in the middle of each MESA and LESA span and nearby but outside of the pivot field (control) to record meteorological parameters on 1 min intervals. The ASCE Penman–Monteith (ASCE-PM) standardized reference equations were used to calculate grass reference evapotranspiration (ETo) from this data on a one-minute basis. A comparison was made for the three phases of before, during, and after the irrigation system passed the in-field ATMOS-41 sensors. In addition, a small unmanned aerial system (UAS) was used to capture 5-band multispectral (ground sampling distance [GSD]: 7 cm/pixel) and thermal infrared images (GSD: 13 cm/pixel) while the center pivot irrigation system was irrigating the field. This imagery data was used to estimate crop evapotranspiration (ETc) using a UAS-METRIC energy balance model. The UAS-METRIC model showed that the estimated ETc under MESA was suppressed by 0.16 mm/day compared to the LESA. Calculating the ETo by the ASCE-PM method showed that the instantaneous ETo rate under the MESA was suppressed between 8% and 18% compared to the LESA. However, as the time of the ET suppression was short, the total amount of the estimated suppressed ET of the MESA was less than 0.5% of the total applied water. Overall, the total reduction in the ET due to the microclimate modifications from wind drift and evaporation losses were small compared to the reported 17% average differences in the irrigation application efficiency between the MESA and the LESA. Therefore, the irrigation application efficiency differences between these two technologies were very large even if the ET suppression by wind drift and evaporation losses was accounted for. Full article
(This article belongs to the Special Issue Evapotranspiration Measurements and Modeling II)
Show Figures

Figure 1

Figure 1
<p>Scheme of the experimental design in 2018 and 2019.</p>
Full article ">Figure 2
<p>ATMOS-41 weather stations installed at 1.8 m above ground under spans of MESA and LESA configurations of the center pivot as it approached the sensors in August 2018.</p>
Full article ">Figure 3
<p>ATMOS-41 weather stations installed at 1 m above ground under spans of MESA and LESA configurations of the center pivot as it passed over the sensors in August 2019.</p>
Full article ">Figure 4
<p>Relative humidity and precipitation of the 20 days of recorded data in 2019 for the MESA, LESA, and Control treatments to find days with irrigation events (the dates with red borders) that were used in this study.</p>
Full article ">Figure 5
<p>Wind rose of the MESA and LESA treatments for day of year 220 and 226 in year 2019.</p>
Full article ">Figure 6
<p>Visual flowchart for processing UAS-based imagery to obtain mapping evapotranspiration at high resolution with internalized calibration (METRIC) energy balance model. The manually selected hot and cold pixels on this flowchart only applies to this field of view.</p>
Full article ">Figure 7
<p>(<b>A</b>) Thermal map of the spearmint field irrigated with MESA (spans 5 and 7) and LESA (6 and 8) irrigation systems. The green triangles on pivot road and the spearmint canopy indicated the hot and cold pixels, respectively, that were used for the internal calibration of the UAS-METRIC model. (<b>B</b>) Spearmint evapotranspiration from the UAS-METRIC model using aerial images taken 100 m above the ground level. The ET<sub>c</sub> map was imported in MATLAB for extracting ROI samples. ROIs were divided into three phases: before, during, and after irrigation by the center pivot system.</p>
Full article ">Figure 8
<p>Measurements of relative humidity and air temperature (1 min interval) relative to the ATMOS-41 installed under MESA and LESA treatments in three phases of before (BI), during (DI), and after irrigation (AI) in 2018. The measured depth of irrigation (mm) is plotted to identify these phases.</p>
Full article ">Figure 9
<p>(<b>A</b>) One-minute basis calculation of reference evapotranspiration (ET<sub>o</sub>) under MESA and LESA treatments relative to the ATMOS-41 in three phases of before (BI), during (DI), and after irrigation (AI). (<b>B</b>) Spread of calculated ET<sub>o</sub> with statistical analysis between MESA and LESA treatments in each phase of irrigation. Different letters indicate that mean value of the measured parameters was significantly different after ANOVA test with a single factor (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 10
<p>One-minute measurements of relative humidities and air temperatures under the MESA and LESA treatments for day of year (DOY) 220 and 226. The phases are shown relative to the irrigation event for 3 h before irrigation (3 h BI) until 3 h after irrigation (3 h AI) the ATMOS-41. The depth of irrigation (mm) is plotted to identify these phases.</p>
Full article ">Figure 11
<p>(<b>A</b>) Monitoring changes in the reference evapotranspiration (ET<sub>o</sub>) under center pivot spans with MESA and LESA irrigation systems while the irrigation system is passing the ATMOS-41. (<b>B</b>) Spread of calculated ET<sub>o</sub> for DOY 220 and 226 with statistical analysis between MESA and LESA treatments in each phase of irrigation. Different letters indicate that mean value of the measured parameters was significantly different after ANOVA test with a single factor (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">
22 pages, 4997 KiB  
Article
Drone-Based Hyperspectral and Thermal Imagery for Quantifying Upland Rice Productivity and Water Use Efficiency after Biochar Application
by Hongxiao Jin, Christian Josef Köppl, Benjamin M. C. Fischer, Johanna Rojas-Conejo, Mark S. Johnson, Laura Morillas, Steve W. Lyon, Ana M. Durán-Quesada, Andrea Suárez-Serrano, Stefano Manzoni and Monica Garcia
Remote Sens. 2021, 13(10), 1866; https://doi.org/10.3390/rs13101866 - 11 May 2021
Cited by 13 | Viewed by 6242
Abstract
Miniature hyperspectral and thermal cameras onboard lightweight unmanned aerial vehicles (UAV) bring new opportunities for monitoring land surface variables at unprecedented fine spatial resolution with acceptable accuracy. This research applies hyperspectral and thermal imagery from a drone to quantify upland rice productivity and [...] Read more.
Miniature hyperspectral and thermal cameras onboard lightweight unmanned aerial vehicles (UAV) bring new opportunities for monitoring land surface variables at unprecedented fine spatial resolution with acceptable accuracy. This research applies hyperspectral and thermal imagery from a drone to quantify upland rice productivity and water use efficiency (WUE) after biochar application in Costa Rica. The field flights were conducted over two experimental groups with bamboo biochar (BC1) and sugarcane biochar (BC2) amendments and one control (C) group without biochar application. Rice canopy biophysical variables were estimated by inverting a canopy radiative transfer model on hyperspectral reflectance. Variations in gross primary productivity (GPP) and WUE across treatments were estimated using light-use efficiency and WUE models respectively from the normalized difference vegetation index (NDVI), canopy chlorophyll content (CCC), and evapotranspiration rate. We found that GPP was increased by 41.9 ± 3.4% in BC1 and 17.5 ± 3.4% in BC2 versus C, which may be explained by higher soil moisture after biochar application, and consequently significantly higher WUEs by 40.8 ± 3.5% in BC1 and 13.4 ± 3.5% in BC2 compared to C. This study demonstrated the use of hyperspectral and thermal imagery from a drone to quantify biochar effects on dry cropland by integrating ground measurements and physical models. Full article
(This article belongs to the Special Issue Ecohydrological Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location of biochar plots at an upland rice experimental site in Costa Rica. The experiment consisted of three groups (BC1: bamboo biochar, BC2: sugarcane biochar, and C: control without biochar) in triplicate (Plot 1, 2, and 3). Each plot was divided into 3 sub-plots used in statistical analysis (indicated as dashed lines). The true-color image (4.5 mm resolution) is from a Cubert FireflEYE 185 VNIR camera (band 56, 29, and 12) taken on 14 November 2018, orthorectified and mosaicked using Agisoft Metashape software. The star (☆) denotes the location of the weather station, and circles (○) the locations of soil sensors.</p>
Full article ">Figure 2
<p>(<b>a</b>) False-color composite image from Cubert hyperspectral camera band 104, 29, and 12. Red color denotes vegetation and blue color bare soil. (<b>b</b>) Comparison of rice canopy reflectance measured by ASD spectroradiometer and Cubert camera after the empirical line correction using (<b>c</b>) reflectance of dark soil object, and (<b>d</b>) reflectance of bright metal box object (black line) in four days measurement (colored lines).</p>
Full article ">Figure 3
<p>Land surface temperature (LST, 2.25 cm resolution) at 12:30 on 16 November 2018 estimated from the FLIR tau2 thermal camera and the Cubert-derived land surface emissivity. The yellow color indicates hot pixels and the dark color cool pixels.</p>
Full article ">Figure 4
<p>(<b>a</b>) Map of volumetric soil water content (θ) of three experimental groups from UAV on 14 November 2018. The circle (○) denote sensor locations. (<b>b</b>) Violin plots of the θ from UAV in each plot and each treatment. The shaded areas of the violin plot denote the distribution of the day’s measurement, the white dots denote the median values, and the bottom and top edges of the black bars denote the 25th and 75th percentiles respectively. The circles (○) denote sensor measurements during flight. (<b>c</b>) The comparison of plot-averaged θ in three groups BC1, BC2, and C from in situ sensors during the four-day campaign. Whiskers denote the 95% confidence interval of each group on each day.</p>
Full article ">Figure 5
<p>(<b>a</b>) Map of soil matric potential (ψ) of three experimental groups from UAV on 14 November 2018. The circle (○) denote sensor locations. (<b>b</b>) Violin plots of the ψ from UAV in each plot and each treatment. The shaded areas of the violin plot denote the distribution of the day’s measurement, the white dots denote the median values, and the bottom and top edges of the black bars denote the 25th and 75th percentiles respectively. The circles (○) denote sensor measurements during flight. (<b>c</b>) The comparison of plot-averaged ψ in three groups BC1, BC2, and C from in situ sensors during the four-day campaign. Whiskers denote the 95% confidence interval of each group on each day.</p>
Full article ">Figure 6
<p>(<b>a</b>) Map of canopy chlorophyll content (CCC) over three experimental groups. The comparison of three groups BC1, BC2, and C is shown in (<b>b</b>) CCC from Cubert, and (<b>c</b>) normalized difference vegetation index (NDVI). Whiskers denote the 95% confidence interval of each group on each day.</p>
Full article ">Figure 7
<p>(<b>a</b>) Map of latent heat flux (λET) over three experimental groups. The comparison of three groups BC1, BC2, and C is shown in (<b>b</b>) latent heat flux (λET), and (<b>c</b>) sensible heat flux (H). Whiskers denote the 95% confidence interval of each group on each day.</p>
Full article ">Figure 8
<p>Overall relationships between soil matric potential (ψ) and (<b>a</b>) canopy chlorophyll content (CCC), (<b>b</b>) normalized difference vegetation index (NDVI), and (<b>c</b>) leaf area index (LAI) in the three groups BC1, BC2, and C.</p>
Full article ">Figure 9
<p>Relationships of evaporative fraction (EF) to (<b>a</b>) soil moisture content (θ) and (<b>b</b>) soil matric potential (ψ) derived from UAV for the three groups BC1, BC2, and C. Dash lines in (<b>a</b>) indicate the extrapolation of the linear fitting lines (solid lines) to maximum EF (EF<sub>max</sub> ≤ 1). Dash lines in (<b>b</b>) indicate the extrapolation of the exponential fitting lines (solid lines) to EF = EF<sub>max</sub>, using sigmoidal curves.</p>
Full article ">
17 pages, 24772 KiB  
Article
High-Resolution Spatiotemporal Water Use Mapping of Surface and Direct-Root-Zone Drip-Irrigated Grapevines Using UAS-Based Thermal and Multispectral Remote Sensing
by Abhilash K. Chandel, Lav R. Khot, Behnaz Molaei, R. Troy Peters, Claudio O. Stöckle and Pete W. Jacoby
Remote Sens. 2021, 13(5), 954; https://doi.org/10.3390/rs13050954 - 4 Mar 2021
Cited by 21 | Viewed by 3977
Abstract
Site-specific irrigation management for perennial crops such as grape requires water use assessments at high spatiotemporal resolution. In this study, small unmanned-aerial-system (UAS)-based imaging was used with a modified mapping evapotranspiration at high resolution with internalized calibration (METRIC) energy balance model to map [...] Read more.
Site-specific irrigation management for perennial crops such as grape requires water use assessments at high spatiotemporal resolution. In this study, small unmanned-aerial-system (UAS)-based imaging was used with a modified mapping evapotranspiration at high resolution with internalized calibration (METRIC) energy balance model to map water use (UASM-ET approach) of a commercial, surface, and direct-root-zone (DRZ) drip-irrigated vineyard. Four irrigation treatments, 100%, 80%, 60%, and 40%, of commercial rate (CR) were also applied, with the CR estimated using soil moisture data and a non-stressed average crop coefficient of 0.5. Fourteen campaigns were conducted in the 2018 and 2019 seasons to collect multispectral (ground sampling distance (GSD): 7 cm/pixel) and thermal imaging (GSD: 13 cm/pixel) data. Six of those campaigns were near Landsat 7/8 satellite overpass of the field site. Weather inputs were obtained from a nearby WSU-AgWeatherNet station (1 km). First, UASM-ET estimates were compared to those derived from soil water balance (SWB) and conventional Landsat-METRIC (LM) approaches. Overall, UASM-ET (2.70 ± 1.03 mm day−1 [mean ± std. dev.]) was higher than SWB-ET (1.80 ± 0.98 mm day−1). However, both estimates had a significant linear correlation (r = 0.64–0.81, p < 0.01). For the days of satellite overpass, UASM-ET was statistically similar to LM-ET, with mean absolute normalized ET departures (ETd,MAN) of 4.30% and a mean r of 0.83 (p < 0.01). The study also extracted spatial canopy transpiration (UASM-T) maps by segmenting the soil background from the UASM-ET, which had strong correlation with the estimates derived by the standard basal crop coefficient approach (Td,MAN = 14%, r = 0.95, p < 0.01). The UASM-T maps were then used to quantify water use differences in the DRZ-irrigated grapevines. Canopy transpiration (T) was statistically significant among the irrigation treatments and was highest for grapevines irrigated at 100% or 80% of the CR, followed by 60% and 40% of the CR (p < 0.01). Reference T fraction (TrF) curves established from the UASM-T maps showed a notable effect of irrigation treatment rates. The total water use of grapevines estimated using interpolated TrF curves was highest for treatments of 100% (425 and 320 mm for the 2018 and 2019 seasons, respectively), followed by 80% (420 and 317 mm), 60% (391 and 318 mm), and 40% (370 and 304 mm) of the CR. Such estimates were within 5% to 11% of the SWB-based water use calculations. The UASM-T-estimated water use was not the same as the actual amount of water applied in the two seasons, probably because DRZ-irrigated vines might have developed deeper or lateral roots to fulfill water requirements outside the irrigated soil volume. Overall, results highlight the usefulness of high-resolution imagery toward site-specific water use management of grapevines. Full article
Show Figures

Figure 1

Figure 1
<p>Soil moisture for surface-drip-irrigated (100% of the commercial rate (CR)) grapevine treatment blocks pertinent to the (<b>a</b>) 2018 and (<b>b</b>) 2019 growth seasons.</p>
Full article ">Figure 2
<p>Unmanned aerial system (UAS) used for multispectral and thermal infrared imaging of grapevines.</p>
Full article ">Figure 3
<p>Small UAS imagery-derived, high-resolution normalized difference vegetation index (NDVI) orthomosaic layer (season: 2019, day before harvest (DBH): 71) showing the surface and direct-root-zone (DRZ) irrigation treatment layout in Block 2.</p>
Full article ">Figure 4
<p>Sample daily ET maps from (<b>a</b>) UASM and (<b>b</b>) LM approaches for grapevines irrigated at the commercial rate (day of the year: 207, season: 2018).</p>
Full article ">Figure 5
<p>Comparison plots between daily ET estimates from UASM and LM approaches for grapevines irrigated at the commercial rate on selected days for the (<b>a</b>) 2018 and (<b>b</b>) 2019 growing seasons. DOY: day of the year.</p>
Full article ">Figure 6
<p>Comparison between daily ET estimates from UASM and SWB approaches on imaging days for the (<b>a</b>) 2018 and (<b>b</b>) 2019 growing seasons.</p>
Full article ">Figure 7
<p>Daily UASM-T for DRZ-irrigated grapevines at different rate and depth treatments on selected days in the 2018 (<b>a</b>–<b>d</b>) and 2019 (<b>e</b>–<b>h</b>) growing seasons. DBH: days before harvest.</p>
Full article ">Figure 7 Cont.
<p>Daily UASM-T for DRZ-irrigated grapevines at different rate and depth treatments on selected days in the 2018 (<b>a</b>–<b>d</b>) and 2019 (<b>e</b>–<b>h</b>) growing seasons. DBH: days before harvest.</p>
Full article ">Figure 8
<p>Basal crop coefficient curves derived from UASM-T maps (non-dotted lines) for grapevines irrigated at different rates in the (<b>a</b>) 2018, and (<b>b</b>) 2019 growth seasons (Ts: treatment start day; S represents locally adjusted standard basal crop coefficient curve pertinent to non-stressed canopy; 40%, 60%, and 80% of the CR were DRZ, and 100% of the CR was surface drip irrigation).</p>
Full article ">Figure 9
<p>Crop water use (or total T) and actual water applied during treatments in the (<b>a</b>) 2018 and (<b>b</b>) 2019 growing seasons (C. ET<sub>r</sub>: cumulative ET<sub>r</sub> for the treatment period).</p>
Full article ">
13 pages, 5839 KiB  
Article
Smartphone Application-Enabled Apple Fruit Surface Temperature Monitoring Tool for In-Field and Real-Time Sunburn Susceptibility Prediction
by Bin Wang, Rakesh Ranjan, Lav R. Khot and R. Troy Peters
Sensors 2020, 20(3), 608; https://doi.org/10.3390/s20030608 - 22 Jan 2020
Cited by 13 | Viewed by 4192
Abstract
Heat stress and resulting sunburn is a major abiotic stress in perineal specialty crops. For example, such stress to the maturing fruits on apple tree canopies can cause several physiological disorders that result in considerable crop losses and reduced marketability of the produce. [...] Read more.
Heat stress and resulting sunburn is a major abiotic stress in perineal specialty crops. For example, such stress to the maturing fruits on apple tree canopies can cause several physiological disorders that result in considerable crop losses and reduced marketability of the produce. Thus, there is a critical technological need to effectively monitor the abiotic stress under field conditions for timely actuation of remedial measures. Fruit surface temperature (FST) is one of the stress indicators that can reliably be used to predict apple fruit sunburn susceptibility. This study was therefore focused on development and in-field testing of a mobile FST monitoring tool that can be used for real-time crop stress monitoring. The tool integrates a smartphone connected thermal-Red-Green-Blue (RGB) imaging sensor and a custom developed application (‘AppSense 1.0’) for apple fruit sunburn prediction. This tool is configured to acquire and analyze imagery data onboard the smartphone to estimate FST. The tool also utilizes geolocation-specific weather data to estimate weather-based FST using an energy balance modeling approach. The ‘AppSense 1.0’ application, developed to work in the Android operating system, allows visual display, annotation and real-time sharing of the imagery, weather data and pertinent FST estimates. The developed tool was evaluated in orchard conditions during the 2019 crop production season on the Gala, Fuji, Red delicious and Honeycrisp apple cultivars. Overall, results showed no significant difference (t110 = 0.51, p = 0.6) between the mobile FST monitoring tool outputs, and ground truth FST data collected using a thermal probe which had accuracy of ±0.4 °C. Upon further refinements, such tool could aid growers in real-time apple fruit sunburn susceptibility prediction and assist in more effective actuation of apple fruit sunburn preventative measures. This tool also has the potential to be customized for in-field monitoring of the heat stressors in some of the sun-exposed perennial and annual specialty crops at produce maturation. Full article
(This article belongs to the Special Issue Advanced Sensor Technologies for Crop Phenotyping Application)
Show Figures

Figure 1

Figure 1
<p>Mobile fruit surface temperature monitoring tool with ‘AppSense 1.0’ Android application designed for rapid and real-time apple sunburn prediction.</p>
Full article ">Figure 2
<p>The ‘AppSense 1.0’ application data analysis flow (blue and red lines represent the weather data based fruit surface temperature (FST) and imaging based FST models, respectively). Indicated icons on the work-flow represent various controls (CAMERA-image acquisition; WEATHER-weather data collection form nearest weather station; FST-estimation of imagery and weather based FST; SAVE-save the data on-board the smartphone, DATABASE-summarize the result and exports to ‘*.xlse’ file).</p>
Full article ">Figure 3
<p>Weather data based (<b>a</b>) FST estimation process flow and (<b>b</b>) location specific weather information and FST estimation results.</p>
Full article ">Figure 4
<p>(<b>a</b>) Flow chart for thermal-Red-Green-Blue (RGB) imaging based FST estimation and (<b>b</b>) steps involved in color feature based fruit segmentation.</p>
Full article ">Figure 5
<p>A typical thermal-RGB image analysis outputs (<b>a</b>) an RGB image, (<b>b</b>) a segmented and binarized image, (<b>c</b>) a binary image after blur filtering, (<b>d</b>) a binary image with the largest contour, (<b>e</b>) an apple blob extraction, (<b>f</b>) a raw thermal image, (<b>g</b>) a pixel offset between the raw overlapped thermal and RGB images, and (<b>h</b>) an offset-adjusted image overlay.</p>
Full article ">Figure 6
<p>(<b>a</b>) Fruit surface temperature monitoring using developed tool and pertinent (<b>b</b>) ground-truthing using a thermal probe.</p>
Full article ">Figure 7
<p>(<b>a</b>) The ‘AppSense 1.0’ application display interface depicting segmentation output and summary of FST estimates and representative segmentation results obtained for (<b>b</b>) Fuji, (<b>c</b>) Gala, (<b>d</b>) Red delicious, and (<b>e</b>) Honeycrisp fruits through real-time imagery data analysis on-board the smartphone running the mobile FST monitoring tool.</p>
Full article ">Figure 8
<p>Thermal-RGB imaging sensor calibration equation for temperature data correction.</p>
Full article ">Figure 9
<p>(<b>a</b>) Thermal-RGB imagery and weather data derived FST estimates and corresponding ground truth measured FST and ambient air temperature data (FST<sub>i</sub>: Thermal-RGB imagery derived FST, FST<sub>i-max</sub>: imagery based maximum FST, FST<sub>w</sub>: weather data derived FST, T<sub>air</sub>: ambient air temperature, FST<sub>G</sub>: ground truth measured FST, and n: sample size), and (<b>b</b>) typical fruit thermal map for <span class="html-italic">cv.</span> Honeycrisp.</p>
Full article ">Figure 10
<p>Relationships between FST<sub>i</sub> and (<b>a</b>) FST<sub>w</sub> and (<b>b</b>) T<sub>air</sub> (FST<sub>i</sub>: Thermal-RGB imagery derived FST, FST<sub>w</sub>: weather data derived FST, T<sub>air</sub>: ambient air temperature).</p>
Full article ">Figure 11
<p>(<b>a</b>) Thermal RGB imagery and weather data derived FST estimates and corresponding ambient air temperature with temperature map for <span class="html-italic">cv.</span> (<b>b</b>) Fuji, (<b>c</b>) Gala and (<b>d</b>) Red delicious (FST<sub>i</sub>: Thermal-RGB imagery derived FST, FST<sub>i-max</sub>: imagery based maximum FST, FST<sub>w</sub>: weather data derived FST, T<sub>air</sub>: ambient air temperature, FST<sub>G</sub>: ground truth measured FST, and n: sample size).</p>
Full article ">
15 pages, 10781 KiB  
Article
A Non-Reference Temperature Histogram Method for Determining Tc from Ground-Based Thermal Imagery of Orchard Tree Canopies
by Arachchige Surantha Ashan Salgadoe, Andrew James Robson, David William Lamb and Derek Schneider
Remote Sens. 2019, 11(6), 714; https://doi.org/10.3390/rs11060714 - 25 Mar 2019
Cited by 13 | Viewed by 4688
Abstract
Obtaining average canopy temperature (Tc) by thresholding canopy pixels from on-ground thermal imagery has historically been undertaken using ‘wet’ and ‘dry’ reference surfaces in the field (reference temperature thresholding). However, this method is extremely time inefficient and can suffer inaccuracies if [...] Read more.
Obtaining average canopy temperature (Tc) by thresholding canopy pixels from on-ground thermal imagery has historically been undertaken using ‘wet’ and ‘dry’ reference surfaces in the field (reference temperature thresholding). However, this method is extremely time inefficient and can suffer inaccuracies if the surfaces are non-standardised or unable to stabilise with the environment. The research presented in this paper evaluates non-reference techniques to obtain average canopy temperature (Tc) from thermal imagery of avocado trees, both for the shaded side and sunlit side, without the need of reference temperature values. A sample of 510 thermal images (from 130 avocado trees) were acquired with a FLIR B250 handheld thermal imaging camera. Two methods based on temperature histograms were evaluated for removing non-canopy-related pixel information from the analysis, enabling Tc to be determined. These approaches included: 1) Histogram gradient thresholding based on temperature intensity changes (HG); and 2) histogram thresholding at one or more standard deviation (SD) above and below the mean. The HG method was found to be more accurate (R2 > 0.95) than the SD method in defining canopy pixels and calculating Tc from each thermal image (shaded and sunlit) when compared to the standard reference temperature thresholding method. The results from this study present an alternative non-reference method for determining Tc from ground-based thermal imagery without the need of calibration surfaces. As such, it offers a more efficient and computationally autonomous method that will ultimately support the greater adoption of non-invasive thermal technologies within a precision agricultural system. Full article
(This article belongs to the Section Remote Sensing in Agriculture and Vegetation)
Show Figures

Figure 1

Figure 1
<p>Flow diagram of methodologies used in the study; HG—histogram gradient thresholding, SD—standard deviation envelope thresholding, and RT—reference temperature thresholding.</p>
Full article ">Figure 2
<p>‘Dry’ and ‘wet’ reference surfaces. (<b>a</b>) Practical implementation of the panel with reference surfaces; (<b>b</b>) extracting temperature by manually drawn polygons on thermal image.</p>
Full article ">Figure 3
<p>Sequence of canopy thermal image processing; (<b>a</b>) original thermal imagery (blue: Sky area; yellow: Tree canopy, bright yellow: Non-leaf material), (<b>b</b>) extracted temperature matrix, (<b>c</b>) unimodal histogram of the temperature matrix, and (<b>d</b>) pictorial interpretation of a threshold thermal image (black: Excluded sky area; white area: Excluded hot stem parts; and colored: Leaf material).</p>
Full article ">Figure 4
<p>Graphical interpretation of the histogram gradient thresholding (HG) technique; (<b>a</b>) calculating ratio pixel change (RPC) in the lower ‘tail’, (<b>b</b>) RPC points with corresponding temperature bin values.</p>
Full article ">Figure 5
<p>Root mean squared error (RMSE) and R<sup>2</sup> of T<sub>c</sub> reported by the histogram gradient (HG) method at different RPC gradients tested in the histogram for (<b>a</b>) shaded-side and (<b>b</b>) sunlit-side images.</p>
Full article ">Figure 6
<p>Shaded-side image thresholding for T<sub>c</sub>, HG vs reference temperature (RT) method.</p>
Full article ">Figure 7
<p>Sunlit-side image thresholding for T<sub>c</sub>, HG vs reference temperature (RT) method.</p>
Full article ">Figure 8
<p>Pictorial representation of outliers (HG vs RT thresholding); (<b>a</b>) original thermal image (coloured: canopy and sky), (<b>b</b>) incorrect threshold by RT method (white: canopy incorrectly excluded; black: excluded sky), (<b>c</b>) correct threshold by HG method (coloured: canopy; black: excluded sky).</p>
Full article ">Figure 9
<p>Thresholding for T<sub>c</sub>. Standard deviation (SD) method vs reference tempererature (RT) thresholding; (<b>a</b>) shaded-side images, (<b>b</b>) sunlit-side images.</p>
Full article ">Figure 10
<p>Comparison of outliers (SD method vs RT method); smaller canopy situation (a, b, and c); (<b>a</b>) original thermal image (coloured: canopy and sky), (<b>b</b>) correct threshold by RT method (coloured: canopy; black: excluded sky), (<b>c</b>) incorrect threshold by SD method (white: canopy incorrectly excluded; black: excluded sky). Larger canopy situation (d, e, and f); (<b>d</b>) original thermal image (coloured: canopy and sky), (<b>e</b>) correct threshold by RT method (coloured: canopy; black: excluded sky), (<b>f</b>) correct threshold by SD method (coloured: canopy; black: excluded sky).</p>
Full article ">Figure 11
<p>Validation results of T<sub>c</sub> from the HG method; (<b>a</b>) shaded-side images, (<b>b</b>) sunlit-side images.</p>
Full article ">
18 pages, 5260 KiB  
Article
Determining the Mechanisms that Influence the Surface Temperature of Urban Forest Canopies by Combining Remote Sensing Methods, Ground Observations, and Spatial Statistical Models
by Shudi Zuo, Shaoqing Dai, Xiaodong Song, Chengdong Xu, Yilan Liao, Weiyin Chang, Qi Chen, Yaying Li, Jianfeng Tang, Wang Man and Yin Ren
Remote Sens. 2018, 10(11), 1814; https://doi.org/10.3390/rs10111814 - 15 Nov 2018
Cited by 10 | Viewed by 3930
Abstract
The spatiotemporal distribution pattern of the surface temperatures of urban forest canopies (STUFC) is influenced by many environmental factors, and the identification of interactions between these factors can improve simulations and predictions of spatial patterns of urban cool islands. This quantitative research uses [...] Read more.
The spatiotemporal distribution pattern of the surface temperatures of urban forest canopies (STUFC) is influenced by many environmental factors, and the identification of interactions between these factors can improve simulations and predictions of spatial patterns of urban cool islands. This quantitative research uses an integrated method that combines remote sensing, ground surveys, and spatial statistical models to elucidate the mechanisms that influence the STUFC and considers the interaction of multiple environmental factors. This case study uses Jinjiang, China as a representative of a city experiencing rapid urbanization. We build up a multisource database (forest inventory, digital elevation models, population, and remote sensing imagery) on a uniform coordinate system to support research into the interactions that influence the STUFC. Landsat-5/8 Thermal Mapper images and meteorological data were used to retrieve the temporal and spatial distributions of land surface temperature. Ground observations, which included the forest management planning inventory and population density data, provided the factors that determine the STUFC spatial distribution on an urban scale. The use of a spatial statistical model (GeogDetector model) reveals the interaction mechanisms of STUFC. Although different environmental factors exert different influences on STUFC, in two periods with different hot spots and cold spots, the patch area and dominant tree species proved to be the main factors contributing to STUFC. The interaction between multiple environmental factors increased the STUFC, both linearly and nonlinearly. Strong interactions tended to occur between elevation and dominant species and were prevalent in either hot or cold spots in different years. In conclusion, the combining of multidisciplinary methods (e.g., remote sensing images, ground observations, and spatial statistical models) helps reveal the mechanism of STUFC on an urban scale. Full article
(This article belongs to the Section Forest Remote Sensing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Flow chart describing the proposed methodology. The data were processed by combining Landsat-5/8 TM remote sensing images, FMPI, and the GeogDetector model.</p>
Full article ">Figure 2
<p>Map of study area. The inset shows the location in China.</p>
Full article ">Figure 3
<p>Changes in FMPI data in Jinjiang City between 2004 and 2014.</p>
Full article ">Figure 4
<p>Comparison of STUFC between on-site observations and Landsat data.</p>
Full article ">Figure 5
<p>Spatial distribution of temperature in 2004 and 2014 obtained using a 2250 m threshold and local Getis–Ord Gi* statistics. Red (GiZscore ≥ 1.65) indicates a hot spot, blue (GiZscore ≤ −1.65) indicates a cold spot, yellow (1.65 &gt; GiZscore &gt; −1.65) indicates a nonsignificant area, and gray indicates altitude. The numbers indicate the main regions of mutual conversion of cold spots, hot spots, and nonsignificant areas.</p>
Full article ">Figure 6
<p>Population density in different cluster regions classified with local Getis–Ord G<sub>i</sub>* statistics, which indicate population changes between 2004 and 2014.</p>
Full article ">Figure 7
<p>Degree of influence of forest attributes, soil, topography, and population on STUFC in Jinjiang City. PA, patch area (i.e., forest-stand area); DS, dominant species; CD, canopy density; SA, stand age (i.e., average forest-stand age); ShI, shape index; SI, site index; SD, soil depth; HD, humus depth; ELE, elevation; SDe, slope degree; SPo, slope position; SDi, slope direction; PD, population; ISP100 percent of impervious surface within 100 m radius.</p>
Full article ">Figure 8
<p>Influence (<span class="html-italic">q</span> value) of interactions in cold spot areas in 2004 and 2014.</p>
Full article ">
1110 KiB  
Article
A Three-Dimensional Index for Characterizing Crop Water Stress
by Jessica A. Torrion, Stephan J. Maas, Wenxuan Guo, James P. Bordovsky and Andy M. Cranmer
Remote Sens. 2014, 6(5), 4025-4042; https://doi.org/10.3390/rs6054025 - 2 May 2014
Cited by 10 | Viewed by 9290
Abstract
The application of remotely sensed estimates of canopy minus air temperature (Tc-Ta) for detecting crop water stress can be limited in semi-arid regions, because of the lack of full ground cover (GC) at water-critical crop stages. Thus, soil background may restrict [...] Read more.
The application of remotely sensed estimates of canopy minus air temperature (Tc-Ta) for detecting crop water stress can be limited in semi-arid regions, because of the lack of full ground cover (GC) at water-critical crop stages. Thus, soil background may restrict water stress interpretation by thermal remote sensing. For partial GC, the combination of plant canopy temperature and surrounding soil temperature in an image pixel is expressed as surface temperature (Ts). Soil brightness (SB) for an image scene varies with surface soil moisture. This study evaluates SB, GC and Ts-Ta and determines a fusion approach to assess crop water stress. The study was conducted (2007 and 2008) on a commercial scale, center pivot irrigated research site in the Texas High Plains. High-resolution aircraft-based imagery (red, near-infrared and thermal) was acquired on clear days. The GC and SB were derived using the Perpendicular Vegetation Index approach. The Ts-Ta was derived using an array of ground Ts sensors, thermal imagery and weather station air temperature. The Ts-Ta, GC and SB were fused using the hue, saturation, intensity method, respectively. Results showed that this method can be used to assess water stress in reference to the differential irrigation plots and corresponding yield without the use of additional energy balance calculation for water stress in partial GC conditions. Full article
(This article belongs to the Special Issue Analysis of Remote Sensing Image Data)
Show Figures

Graphical abstract

Graphical abstract
Full article ">
<p>(<b>Top</b>) Cumulative crop evapotranspiration (ET), rainfall and applied total water (irrigation plus rainfall) for the high irrigated (HighIrr), low irrigated (LowIrr) and non-irrigated (dryland) treatments plotted <span class="html-italic">vs.</span> the day of year (DOY); (<b>Bottom</b>) the amount and timing of daily precipitation events. The results are presented for the 2007 and 2008 growing seasons from planting (vertical line with date) to crop maturity. Arrows indicate the dates on which airborne remote sensing imagery was acquired.</p>
Full article ">
<p>Array of sensors used to monitor surface temperature (<span class="html-italic">Ts</span>): (<b>left</b>) high irrigation treatment; (<b>right</b>) dryland treatment.</p>
Full article ">
<p>Typical scatter plot of image pixel reflectances in the red and near-infrared bands.</p>
Full article ">
<p>Red <span class="html-italic">vs.</span> NIR scatter plot of pixel reflectance from an image scene taken on 8 August 2007 (<b>A</b>), and the soil line evaluated from multi-temporal image data for georeferenced bare soil surfaces (near the pivot pump and adjacent to the plots) with varying degrees of moisture (<b>B</b>).</p>
Full article ">
<p>Percent of ground cover (GC) of cotton (C) and sorghum, <span class="html-italic">Sorghum bicolor</span> (S), percent soil brightness (SB) and canopy minus air temperature (<span class="html-italic">Tc-Ta</span>) (°C), for two dates in 2007 (<b>A</b>,<b>B</b>) and 2008 (<b>C</b>) along with the respective irrigation plots for the highly irrigated (HighIrr), low irrigated (LowIrr) and dryland treatment.Note: sorghum was part of a 6-year (2003–2008) crop rotation and irrigation project and is out of the scope of this 2-year (2007–2008) study.</p>
Full article ">
<p>Estimated <span class="html-italic">vs.</span> measured ground cove (GC) for the various irrigation plots in the research site (2007–2008).</p>
Full article ">
<p>Color composites produced by hue, saturation and intensity (HIS) fusion of <span class="html-italic">Ts-Ta</span> (−30 °C to 30 °C), GC (0–100%) and SB (0–100%) for 2007 (<b>A</b>,<b>B</b>) and 2008 (<b>C</b>).</p>
Full article ">
<p>Lint yield for the three irrigation treatments in each year [<a href="#b29-remotesensing-06-04025" class="html-bibr">29</a>]. The standard error of estimates was calculated using the Procmixed model in SAS [<a href="#b30-remotesensing-06-04025" class="html-bibr">30</a>] with the irrigation treatment as the fixed effect and the year as a random effect. The yield difference between treatments greater than the least significant difference (LSD<sub>0.05</sub>) is significant.</p>
Full article ">
<p>GC, SB and <span class="html-italic">Ts-Ta</span> plotted on three orthogonal axes. GC is associated with crop growth, SB with apparent soil moistness and <span class="html-italic">Ts-Ta</span> with the surface energy balance. Plotted data are for the 8 August 2007, image acquisition date.</p>
Full article ">
1381 KiB  
Article
Deciphering the Precision of Stereo IKONOS Canopy Height Models for US Forests with G-LiHT Airborne LiDAR
by Christopher S. R. Neigh, Jeffrey G. Masek, Paul Bourget, Bruce Cook, Chengquan Huang, Khaldoun Rishmawi and Feng Zhao
Remote Sens. 2014, 6(3), 1762-1782; https://doi.org/10.3390/rs6031762 - 26 Feb 2014
Cited by 41 | Viewed by 9197
Abstract
Few studies have evaluated the precision of IKONOS stereo data for measuring forest canopy height. The high cost of airborne light detection and ranging (LiDAR) data collection for large area studies and the present lack of a spaceborne instrument lead to the need [...] Read more.
Few studies have evaluated the precision of IKONOS stereo data for measuring forest canopy height. The high cost of airborne light detection and ranging (LiDAR) data collection for large area studies and the present lack of a spaceborne instrument lead to the need to explore other low cost options. The US Government currently has access to a large archive of commercial high-resolution imagery, which could be quite valuable to forest structure studies. At 1 m resolution, we here compared canopy height models (CHMs) and height data derived from Goddard’s airborne LiDAR Hyper-spectral and Thermal Imager (G-LiHT) with three types of IKONOS stereo derived digital surface models (DSMs) that estimate CHMs by subtracting National Elevation Data (NED) digital terrain models (DTMs). We found the following in three different forested regions of the US after excluding heterogeneous and disturbed forest samples: (1) G-LiHT DTMs were highly correlated with NED DTMs with R2 > 0.98 and root mean square errors (RMSEs) < 2.96 m; (2) when using one visually identifiable ground control point (GCP) from NED, G-LiHT DSMs and IKONOS DSMs had R2 > 0.84 and RMSEs of 2.7 to 4.1 m; and (3) one GCP CHMs for two study sites had R2 > 0.7 and RMSEs of 2.6 to 3 m where data were collected less than four years apart. Our results suggest that IKONOS stereo data are a useful LiDAR alternative where high-quality DTMs are available. Full article
Show Figures


<p>Upper left shows the NASA’s Blue Marble 500 m imagery of the A–C study area locations. Panels (<b>A</b>) Harvard Forest, Massachusetts; (<b>B</b>) Hoquiam, Washington; and (<b>C</b>) Jamison, South Carolina show I-cubed 1 m seamless color mosaic of commercial and government imagery with tiled G-LiHT products indicated with white quadrangles and IKONOS stereo footprints indicated with red quadrangles.</p>
Full article ">
<p>Schema of how data were compared to decipherer components that impact IKONOS canopy height model (CHM) error. The left portion describes data used in analysis, center describes derived data products, and the right portion describes the comparison of products.</p>
Full article ">
<p>Scatter plots of three regions using a standard normal distribution of 25,000 randomly sampled points after filtering from <a href="#t2-remotesensing-06-01762" class="html-table">Table 2</a> comparing NED DTM <span class="html-italic">vs.</span> G-LiHT DTM. The black line indicates a one-to-one relationship and the thin gray line indicates a least squares linear fit. The total number of samples is displayed in the plot title. Black points indicate no filtering, red points indicate data filtered with high LiDAR reflectance and blue points indicate data filtered with high LiDAR reflectance and slopes less than 10°.</p>
Full article ">
<p>Scatter plots of three study regions using a standard normal distribution of pixels randomly sampled with totals indicated in plot title. Three types of DSMs from IKONOS were compared to G-LiHT LiDAR first returns. Values not meeting criteria in <a href="#t2-remotesensing-06-01762" class="html-table">Table 2</a> were excluded from analysis. The colors from left to right indicate no ground control points (GCPs) in black, one GCP in red, and 16 GCPs in blue. The black line indicates a one-to-one relationship and the gray line indicates a least squares linear fit.</p>
Full article ">
<p>Height difference of LiDAR DSMs minus IKONOS DSMs. Data are presented as a histogram count percent of total from the standard normal distribution of randomly sampled pixels. The colors indicate no GCP DSMs with dashed black lines, one GCP DSMs with solid red lines, and 16 GCP DSMs are shown with dashed blue lines. Values not meeting criteria in <a href="#t2-remotesensing-06-01762" class="html-table">Table 2</a> were excluded from analysis.</p>
Full article ">
<p>CHM scatter plots of three study regions with total number of random samples indicated in plot title. Values not meeting valid criteria listed in <a href="#t2-remotesensing-06-01762" class="html-table">Table 2</a> were excluded from analysis. The solid black line indicates a one-to-one relationship and the solid gray line indicates a least squares linear fit. Colors from left to right indicate no GCPs in black, one GCP in red, and 16 GCPs in blue.</p>
Full article ">
<p>A comparison between IKONOS CHM and G-LiHT CHM near Hoquiam Washington south of Grays Harbor. (<b>A</b>) Pan sharpened multi-spectral image (MSI) true color IKONOS with patterns of even aged forest and clear cut harvest. (<b>B</b>) Landsat vegetation change tracker annual disturbance history from 1984 to 2010. (<b>C</b>) IKONOS CHM with 16 GCPs. (<b>D</b>) G-LiHT CHM. (<b>E</b>) Oblique 3-D image of IKONOS MSI fused with G-LiHT CHM.</p>
Full article ">
Back to TopTop