드론을 활용한 위성 지표반사도 산출물 공간 패턴 분석

학위논문(석사) -- 서울대학교대학원 : 농업생명과학대학 생태조경·지역시스템공학부(생태조경학), 2021.8. 조대솔.High-resolution satellites are assigned to monitor land surface in detail. The reliable surface reflectance (SR) is the fundamental in terrestrial ecosystem modeling so the temporal and spatial validation is essential. Usually based on multiple ground control points (GCPs), field spectroscopy guarantees the temporal continuity. Due to limited sampling, however, it hardly illustrates the spatial pattern. As a map, the pixelwise spatial variability of SR products is not well-documented. In this study, we introduced drone-based hyperspectral image (HSI) as a reference and compared the map with Sentinel 2 and Landsat 8 SR products on a heterogeneous rice paddy landscape. First, HSI was validated by field spectroscopy and swath overlapping, which assured qualitative radiometric accuracy within the viewing geometry. Second, HSI was matched to the satellite SRs. It involves spectral and spatial aggregation, co-registration and nadir bidirectional reflectance distribution function (BRDF)-adjusted reflectance (NBAR) conversion. Then, we 1) quantified the spatial variability of the satellite SRs and the vegetation indices (VIs) including NDVI and NIRv by APU matrix, 2) qualified them pixelwise by theoretical error budget and 3) examined the improvement by BRDF normalization. Sentinel 2 SR exhibits overall good agreement with drone HSI while the two NIRs are biased up to 10%. Despite the bias in NIR, the NDVI shows a good match on vegetated areas and the NIRv only displays the discrepancy on built-in areas. Landsat 8 SR was biased over the VIS bands (-9 ~ -7.6%). BRDF normalization just contributed to a minor improvement. Our results demonstrate the potential of drone HSI to replace in-situ observation and evaluate SR or atmospheric correction algorithms over the flat terrain. Future researches should replicate the results over the complex terrain and canopy structure (i.e. forest).원격탐사에서 지표 반사도(SR)는 지표정보를 비파괴적이고 즉각적인 방법으로 전달해주는 매개체 역할을 한다. 신뢰할 수 있는 SR은 육상 생태계 모델링의 기본이고, 이에 따라 SR의 시공간적 검증이 요구된다. 일반적으로 SR은 여러 지상 기준점(GCP)을 기반으로 하는 현장 분광법을 통해서 시간적 연속성이 보장된다. 그러나 현장 분광법은 제한적인 샘플링으로 공간 패턴을 거의 보여주지 않아, 위성 SR의 픽셀 별 공간 변동성은 잘 분석되지 않았다. 본 연구에서는 드론 기반의 초분광 영상(HSI)을 참고자료로 도입하여, 이를 이질적인 논 경관에서 Sentinel 2 및 Landsat 8 SR과 비교하였다. 우선, 드론 HSI는 현장 분광법 및 경로 중첩을 통해서 관측각도 범위 내에서 정성적인 방사 측정을 보장한다고 검증되었다. 이후, 드론 HSI는 위성 SR의 분광반응특성, 공간해상도 및 좌표계를 기준으로 맞춰졌고, 관측 기하를 통일하기 위해서 드론 HIS와 위성 SR은 각각 양방향반사율분포함수 (BRDF) 정규화 반사도 (NBAR)로 변환되었다. 마지막으로, 1) APU 행렬으로 위성 SR과 NDVI, NIRv를 포함하는 식생지수(VI)의 공간변동성을 정량화 했고, 2) 대기보정의 이론적 오차를 기준으로 SR과 VI를 픽셀별로 평가했고, 3) BRDF 정규화를 통한 개선 사항을 검토했다. Sentinel 2 SR은 드론 HSI와 전반적으로 좋은 일치를 보이나, 두 NIR 채널은 최대 10% 편향되었다. NIR의 편향은 식생지수에서 토지 피복에 따라 다른 영향을 미쳤다. NDVI는 식생에서는 낮은 편향을 보여줬고, NIRv는 도시시설물 영역에서만 높은 편향을 보였다. Landsat 8 SR은 VIS 채널에 대해 편향되었다 (-9 ~ -7.6%). BRDF 정규화는 위성 SR의 품질을 개선했지만, 그 영향은 부수적이었다. 본 연구에서는 평탄한 지형에서 드론 HSI가 현장 관측을 대체할 수 있고, 따라서 위성 SR이나 대기보정 알고리즘을 평가하는데 활용될 수 있다는 것을 보였다. 향후 연구에서는 산림으로 대상지를 확대하여, 지형과 캐노피 구조가 드론 HSI 및 위성 SR에 미치는 영향을 분석할 필요가 있다.Chapter 1. Introduction 1 1.1 Background 1 Chapter 2. Method 3 2.1 Study Site 3 2.2 Drone campaign 4 2.3 Data processing 4 2.3.1 Sensor calibration 5 2.3.2 Bidirectional reflectance factor (BRF) calculation 7 2.3.3 BRDF correction 7 2.3.4 Orthorectification 8 2.3.5 Spatial Aggregation 9 2.3.6 Co-registration 10 2.4 Satellite dataset 10 2.4.2 Landsat 8 12 Chapter 3. Result and Discussion 12 3.1 Drone BRF map quality assessment 12 3.1.1 Radiometric accuracy 12 3.1.2 BRDF effect 15 3.2 Spatial variability in satellite surface reflectance product 16 3.2.1 Sentinel 2B (10m) 17 3.2.2 Sentinel 2B (20m) 22 3.2.3 Landsat 8 26 Chapter 4. Conclusion 28 Supplemental Materials 30 Bibliography 34 Abstract in Korean 43석

Mapping the fractional coverage of the invasive shrub Ulex europaeus with multi-temporal Sentinel-2 imagery utilizing UAV orthoimages and a new spatial optimization approach

Mapping the occurrence patterns of invasive plant species and understanding their invasion dynamics is a crucial requirement for preventing further spread to so far unaffected regions. An established approach to map invasive species across large areas is based on the combination of satellite or aerial remote sensing data with ground truth data from fieldwork. Unmanned aerial vehicles (UAV, also referred to as unmanned aerial systems (UAS)) may represent an interesting and low-cost alternative to labor-intensive fieldwork. Despite the increasing use of UAVs in the field of remote sensing in the last years, operational methods to combine UAV and satellite data are still sparse. Here, we present a new methodological framework to estimate the fractional coverage (FC%) of the invasive shrub species Ulex europaeus (common gorse) on Chilo´e Island (south-central Chile), based on ultrahigh- resolution UAV images and a medium resolution intra-annual time-series of Sentinel-2. Our framework is based on three steps: 1) Land cover classification of the UAV orthoimages, 2) reduce the spatial shift between UAV-based land cover classification maps and Sentinel-2 imagery and 3) identify optimal satellite acquisition dates for estimating the actual distribution of Ulex europaeus. In Step 2 we translate the challenging co-registration task between two datasets with very different spatial resolutions into an (machine learning) optimization problem where the UAV-based land cover classification maps obtained in Step 1 are systematically shifted against the satellite images. Based on several Random Forest (RF) models, an optimal fit between varying land cover fractions and the spectral information of Sentinel-2 is identified to correct the spatial offset between both datasets. Considering the spatial shifts of the UAV orthoimages and using optimally timed Sentinel-2 acquisitions led to a significant improvement for the estimation of the current distribution of Ulex europaeus. Furthermore, we found that the Sentinel-2 acquisition from November (flowering time of Ulex europaeus) was particularly important in distinguishing Ulex europaeus from other plant species. Our mapping results could support local efforts in controlling Ulex europaeus. Furthermore, the proposed workflow should be transferable to other use cases where individual target species that are visually detectable in UAV imagery are considered. These findings confirm and underline the great potential of UAV-based groundtruth data for detecting invasive species

Atmospheric Correction Inter-comparison eXercise, ACIX-II Land: An assessment of atmospheric correction processors for Landsat 8 and Sentinel-2 over land

The correction of the atmospheric effects on optical satellite images is essential for quantitative and multitemporal remote sensing applications. In order to study the performance of the state-of-the-art methods in an integrated way, a voluntary and open-access benchmark Atmospheric Correction Inter-comparison eXercise (ACIX) was initiated in 2016 in the frame of Committee on Earth Observation Satellites (CEOS) Working Group on Calibration & Validation (WGCV). The first exercise was extended in a second edition wherein twelve atmospheric correction (AC) processors, a substantially larger testing dataset and additional validation metrics were involved. The sites for the inter-comparison analysis were defined by investigating the full catalogue of the Aerosol Robotic Network (AERONET) sites for coincident measurements with satellites' overpass. Although there were more than one hundred sites for Copernicus Sentinel-2 and Landsat 8 acquisitions, the analysis presented in this paper concerns only the common matchups amongst all processors, reducing the number to 79 and 62 sites respectively. Aerosol Optical Depth (AOD) and Water Vapour (WV) retrievals were consequently validated based on the available AERONET observations. The processors mostly succeeded in retrieving AOD for relatively light to medium aerosol loading (AOD 90% of the results falling within the suggested empirical specifications and with the Root Mean Square Error (RMSE) being mostly <0.25 g/cm2. Regarding Surface Reflectance (SR) validation two main approaches were followed. For the first one, a simulated SR reference dataset was computed over all of the test sites by using the 6SV (Second Simulation of the Satellite Signal in the Solar Spectrum vector code) full radiative transfer modelling (RTM) and AERONET measurements for the required aerosol variables and water vapour content. The performance assessment demonstrated that the retrievals were not biased for most of the bands. The uncertainties ranged from approximately 0.003 to 0.01 (excluding B01) for the best performing processors in both sensors' analyses. For the second one, measurements from the radiometric calibration network RadCalNet over La Crau (France) and Gobabeb (Namibia) were involved in the validation. The performance of the processors was in general consistent across all bands for both sensors and with low standard deviations (<0.04) between on-site and estimated surface reflectance. Overall, our study provides a good insight of AC algorithms' performance to developers and users, pointing out similarities and differences for AOD, WV and SR retrievals. Such validation though still lacks of ground-based measurements of known uncertainty to better assess and characterize the uncertainties in SR retrievals

Application of Sentinel-2 MSI in Arctic Research: Evaluating the Performance of Atmospheric Correction Approaches Over Arctic Sea Ice

Multispectral remote sensing may be a powerful tool for areal retrieval of biogeophysical parameters in the Arctic sea ice. The MultiSpectral Instrument on board the Sentinel-2 (S-2) satellites of the European Space Agency offers new possibilities for Arctic research; S-2A and S-2B provide 13 spectral bands between 443 and 2,202 nm and spatial resolutions between 10 and 60 m, which may enable the monitoring of large areas of Arctic sea ice. For an accurate retrieval of parameters such as surface albedo, the elimination of atmospheric influences in the data is essential. We therefore provide an evaluation of five currently available atmospheric correction processors for S-2 (ACOLITE, ATCOR, iCOR, Polymer, and Sen2Cor). We evaluate the results of the different processors using in situ spectral measurements of ice and snow and open water gathered north of Svalbard during RV Polarstern cruise PS106.1 in summer 2017. We used spectral shapes to assess performance for ice and snow surfaces. For open water, we additionally evaluated intensities. ACOLITE, ATCOR, and iCOR performed well over sea ice and Polymer generated the best results over open water. ATCOR, iCOR and Sen2Cor failed in the image-based retrieval of atmospheric parameters (aerosol optical thickness, water vapor). ACOLITE estimated AOT within the uncertainty range of AERONET measurements. Parameterization based on external data, therefore, was necessary to obtain reliable results. To illustrate consequences of processor selection on secondary products we computed average surface reflectance of six bands and normalized difference melt index (NDMI) on an image subset. Medians of average reflectance and NDMI range from 0.80–0.97 to 0.12–0.18 while medians for TOA are 0.75 and 0.06, respectively

Détection automatisée du réseau routier en forêt boréale par télédétection

Les routes forestières sont essentielles pour l’aménagement forestier durable, il est important pour les gestionnaires des forêts de détenir l’information adéquate du réseau routier dans leur prise de décision. Ce projet a permis d’évaluer l’apport de trois approches orientées objets d’extraction du réseau routier en forêt boréale à partir de; 1) LiDAR aéroporté à 1 mètre de résolution spatiale; 2) l’image satellitaire Sentinel-2 à 10 mètres de résolution et 3) la fusion de deux sources de données précédentes. Le but étant non seulement d’estimer la contribution individuelle de chaque donnée, mais aussi de compléter les informations morphologiques sur les routes forestières afin de mettre à jour les bases de données géographiques disponibles, à partir d’un processus de détection automatique par télédétection. En effet, les bases de données disponibles sont sujettes des incohérences dues aux problèmes liés à la structuration des données reflétant la réalité de terrain ainsi que des limitations sur le géoréférencement qui affectent la prise de décision sur l’ensemble du territoire forestier. Avec l’essor de la technologie de produits de télédétection et de systèmes d’information géographique, nous proposons dans cette étude, une approche de classification automatique basée sur les objets pour l’identification et la caractérisation automatique des chemins en forêt boréale. La segmentation multirésolution a été appliquée aux trois approches sur trois zones d’étude situés au Québec. Les objets linéaires détectés ont été construits de manière itérative en objets linéaires routes par utilisation supplémentaire de la segmentation basée sur la différence spectrale. Les objets linéaires routes ont ainsi été classifiés en réseau routier à partir d’un jeu de règles, définissant ainsi le processus de la classification orientée objet. Les données des composantes morphologiques des routes (largeur de l’emprise et pente longitudinale) ont été extraites du réseau routier pour servir de caractérisation automatique des chemins forestiers. Cette approche méthodologique minimise les incohérences d’informations retrouvées dans les outils cartographiques actuellement disponibles (Routard) et contribue à la documentation sur la morphologie (qualité de l’information sur la géométrie) des données routières disponibles. Les résultats montrent pour les trois zones d’étude que l’approche utilisant la fusion de données Sentinel-2 et LiDAR améliore considérablement les performances de la précision globale (88%) quant à la détection de réseau routier par rapport à l’approche basée sur Sentinel-2 (70%) et celle basée sur LiDAR (63%). Les résultats obtenus sont présentés sous la forme d’une couche vectorielle dans une base de données d’information géographique pour un territoire d’étude

Combination of optical and SAR remote sensing data for wetland mapping and monitoring

Wetlands provide many services to the environment and humans. They play a pivotal role in water quality, climate change, as well as carbon and hydrological cycles. Wetlands are environmental health indicators because of their contributions to plant and animal habitats. While a large portion of Newfoundland and Labrador (NL) is covered by wetlands, no significant efforts had been conducted to identify and monitor these valuable environments when I initiated this project. At that time, there were only two small areas in NL that had been classified using basic Remote Sensing (RS) methods with low accuracies. There was an immediate need to develop new methods for conserving and managing these vital resources using up-to-date maps of wetland distributions. In this thesis, object- and pixel-based classification methods were compared to show the high potential of the former method when medium or high spatial resolution imagery were used to classify wetlands. The maps produced using several classification algorithms were also compared to select the optimum classifier for future experiments. Moreover, a novel Multiple Classifier System (MCS), which combined several algorithms, was proposed to increase the classification accuracy of complex and similar land covers, such as wetlands. Landsat-8 images captured in different months were also investigated to select the time, for which wetlands had the highest separability using the Random Forest (RF) algorithm. Additionally, various spectral, polarimetric, texture, and ratio features extracted from multi-source optical and Synthetic Aperture Radar (SAR) data were assessed to select the most effective features for discriminating wetland classes. The methods developed during this dissertation were validated in five study areas to show their effectiveness. Finally, in collaboration with a team, a website (http://nlwetlands.ca/) and a software package were developed (named the Advanced Remote Sensing Lab (ARSeL)) to automatically preprocess optical/SAR data and classify wetlands using advanced algorithms. In summary, the outputs of this work are promising and can be incorporated into future studies related to wetlands. The province can also benefit from the results in many ways

Derivation of forest inventory parameters from high-resolution satellite imagery for the Thunkel area, Northern Mongolia. A comparative study on various satellite sensors and data analysis techniques.

With the demise of the Soviet Union and the transition to a market economy starting in the 1990s, Mongolia has been experiencing dramatic changes resulting in social and economic disparities and an increasing strain on its natural resources. The situation is exacerbated by a changing climate, the erosion of forestry related administrative structures, and a lack of law enforcement activities. Mongolia’s forests have been afflicted with a dramatic increase in degradation due to human and natural impacts such as overexploitation and wildfire occurrences. In addition, forest management practices are far from being sustainable. In order to provide useful information on how to viably and effectively utilise the forest resources in the future, the gathering and analysis of forest related data is pivotal. Although a National Forest Inventory was conducted in 2016, very little reliable and scientifically substantiated information exists related to a regional or even local level. This lack of detailed information warranted a study performed in the Thunkel taiga area in 2017 in cooperation with the GIZ. In this context, we hypothesise that (i) tree species and composition can be identified utilising the aerial imagery, (ii) tree height can be extracted from the resulting canopy height model with accuracies commensurate with field survey measurements, and (iii) high-resolution satellite imagery is suitable for the extraction of tree species, the number of trees, and the upscaling of timber volume and basal area based on the spectral properties. The outcomes of this study illustrate quite clearly the potential of employing UAV imagery for tree height extraction (R2 of 0.9) as well as for species and crown diameter determination. However, in a few instances, the visual interpretation of the aerial photographs were determined to be superior to the computer-aided automatic extraction of forest attributes. In addition, imagery from various satellite sensors (e.g. Sentinel-2, RapidEye, WorldView-2) proved to be excellently suited for the delineation of burned areas and the assessment of tree vigour. Furthermore, recently developed sophisticated classifying approaches such as Support Vector Machines and Random Forest appear to be tailored for tree species discrimination (Overall Accuracy of 89%). Object-based classification approaches convey the impression to be highly suitable for very high-resolution imagery, however, at medium scale, pixel-based classifiers outperformed the former. It is also suggested that high radiometric resolution bears the potential to easily compensate for the lack of spatial detectability in the imagery. Quite surprising was the occurrence of dark taiga species in the riparian areas being beyond their natural habitat range. The presented results matrix and the interpretation key have been devised as a decision tool and/or a vademecum for practitioners. In consideration of future projects and to facilitate the improvement of the forest inventory database, the establishment of permanent sampling plots in the Mongolian taigas is strongly advised.2021-06-0