Developing detailed age-specific thematic maps for coffee (Coffea arabica L.) in heterogeneous agricultural landscapes using random forests applied on Landsat 8 multispectral sensor

Coffee is a commodity of international trade significance, and its value chain can benefit from age-specific thematic maps. This study aimed to assess the potential of Landsat 8 OLI to develop these maps. Using field-collected samples with the random forest classifier, splitting coffee into three age classes (Scheme A) was compared with running the classification with one compound coffee class (Scheme B). Higher overall classification accuracy was obtained in Scheme B (90.3% for OLI and 86.8% for ETM+) than in Scheme A (86.2% for OLI and 81.0% for ETM+). The NIR band of OLI was the most important band in intra-class discrimination of coffee. Landsat 8 OLI mapped area closely matched farm records (R^2?=?0.88) compared to that of Landsat 7 ETM+ (R^2?=?0.78). It was concluded that Landsat 8 OLI data can be used to produce age-specific thematic maps in coffee production areas although disaggregating coffee classes reduces overall accuracy.     

Detection and mapping the spatial distribution of bracken fern weeds using the Landsat 8 OLI new generation sensor

Bracken fern is an invasive plant that presents serious environmental, ecological and economic problems around the world. An understanding of the spatial distribution of bracken fern weeds is therefore essential for providing appropriate management strategies at both local and regional scales. The aim of this study was to assess the utility of the freely available medium resolution Landsat 8 OLI sensor in the detection and mapping of bracken fern at the Cathedral Peak, South Africa. To achieve this objective, the results obtained from Landsat 8 OLI were compared with those derived using the costly, high spatial resolution WorldView-2 imagery. Since previous studies have already successfully mapped bracken fern using high spatial resolution WorldView-2 image, the comparison was done to investigate the magnitude of difference in accuracy between the two sensors in relation to their acquisition costs. To evaluate the performance of Landsat 8 OLI in discriminating bracken fern compared to that of Worldview-2, we tested the utility of (i) spectral bands; (ii) derived vegetation indices as well as (iii) the combination of spectral bands and vegetation indices based on discriminant analysis classification algorithm. After resampling the training and testing data and reclassifying several times (n = 100) based on the combined data sets, the overall accuracies for both Landsat 8 and WorldView-2 were tested for significant differences based on Mann-Whitney U test. The results showed that the integration of the spectral bands and derived vegetation indices yielded the best overall classification accuracy (80.08% and 87.80% for Landsat 8 OLI and WorldView-2 respectively). Additionally, the use of derived vegetation indices as a standalone data set produced the weakest overall accuracy results of 62.14% and 82.11% for both the Landsat 8 OLI and WorldView-2 images. There were significant differences {U (100) = 569.5, z = −10.8242, p < 0.01} between the classification accuracies derived based on Landsat OLI 8 and those derived using WorldView-2 sensor. Although there were significant differences between Landsat and WorldView-2 accuracies, the magnitude of variation (9%) between the two sensors was within an acceptable range. Therefore, the findings of this study demonstrated that the recently launched Landsat 8 OLI multispectral sensor provides valuable information that could aid in the long term continuous monitoring and formulation of effective bracken fern management with acceptable accuracies that are comparable to those obtained from the high resolution WorldView-2 commercial sensor.     

Sentinel-2A与Landsat 8O LI逐像元辐射归一化方法研究

考虑不同传感器光谱响应函数差异及不同地物类型反射率光谱的差异,提出了一种逐像元辐射归一化方法,并以2017年7月17日内蒙古达里诺尔湖地区准同步过境的Sentinel-2A及Landsat 8数据为例,对两类数据可见-近红外波段（VNIR）地表反射率结果进行归一化。首先采用Sen2cor方法及NASA官方提供大气校正算法,分别对Sentinel-2A及Landsat 8 OLI影像进行大气校正并重采样到同一空间分辨率;然后基于光谱库计算匹配因子并构建图像与光谱库之间的匹配转换模型,实现像元尺度上从Sentinel-2影像到Landsat 8影像地表反射率相似波段之间的转换。结果表明,经逐像元归一化的影像相比原始影像及经HLS光谱归一化的影像,与Landsat 8 VNIR波段的相关性明显提高,辐射一致性增强。该转换模型为多源中高分辨率遥感图像高精度辐射归一化提供了新思路。     

Comprehensive study of the biophysical parameters of agricultural crops based on assessing Landsat 8 OLI and Landsat 7 ETM+ vegetation indices

For three agricultural crop types, winter wheat (Triticum aestivum L.), barley (Hordeum vulgare L.), and canola (Brassica napus L.), we estimated biophysical parameters including fresh and dry biomass, leaf area index (LAI), and vegetation water content, for which we found the equivalent water thickness (EWT), fuel moisture content per fresh weight (FMC_FW), and fuel moisture content per dry weight (FMC_DW). We performed these estimations using data from the newly launched Landsat 8 Operational Land Imager (OLI) sensor, as well as its predecessor the Landsat 7 Enhanced Thematic Mapper Plus (ETM+). Progress in the design of the new sensor (i.e., Landsat 8), including narrower near-infrared (NIR) wavebands, higher signal-to-noise ratio (SNR), and greater radiometric resolution highlights the necessity to investigate the biophysical parameters of agricultural crops, especially compared to data from its predecessor. This study aims to evaluate vegetation indices (VIs) derived from the Landsat 8 OLI and the Landsat 7 ETM+. Both the Landsat 8 OLI and Landsat 7 ETM+ VIs agreed well with in-situ data measurements. However, the Landsat 8 OLI-derived VIs were generally more consistent with in situ data than the Landsat 7 ETM+ VIs. We also note that the Landsat 8 OLI is better able to capture the small variability of the VIs because of its higher SNR and wider radiometric range; in addition, the saturation phenomenon occurred earlier for the Landsat 7 ETM+ than for the Landsat 8 OLI. This indicates that the new sensor is better able to estimate the biophysical parameters of crops.     

Using multi-date satellite imagery to monitor invasive grass species distribution in post-wildfire landscapes: An iterative,adaptable approach that employs open-source data and software

Among the most pressing concerns of land managers in post-wildfire landscapes are the establishment and spread of invasive species. Land managers need accurate maps of invasive species cover for targeted management post-disturbance that are easily transferable across space and time. In this study, we sought to develop an iterative, replicable methodology based on limited invasive species occurrence data, freely available remotely sensed data, and open source software to predict the distribution of Bromus tectorum (cheatgrass) in a post-wildfire landscape. We developed four species distribution models using eight spectral indices derived from five months of Landsat 8 Operational Land Imager (OLI) data in 2014. These months corresponded to both cheatgrass growing period and time of field data collection in the study area. The four models were improved using an iterative approach in which a threshold for cover was established, and all models had high sensitivity values when tested on an independent dataset. We also quantified the area at highest risk for invasion in future seasons given 2014 distribution, topographic covariates, and seed dispersal limitations. These models demonstrate the effectiveness of using derived multi-date spectral indices as proxies for species occurrence on the landscape, the importance of selecting thresholds for invasive species cover to evaluate ecological risk in species distribution models, and the applicability of Landsat 8 OLI and the Software for Assisted Habitat Modeling for targeted invasive species management.     

Land cover classification using Landsat 8 Operational Land Imager data in Beijing,China

The successful launch of Landsat 8 provides a new data source for monitoring land cover, which has the potential to significantly improve the characterization of the earth’s surface. To assess data performance, Landsat 8 Operational Land Imager (OLI) data were first compared with Landsat 7 ETM + data using texture features as the indicators. Furthermore, the OLI data were investigated for land cover classification using the maximum likelihood and support vector machine classifiers in Beijing. The results indicated that (1) the OLI data quality was slightly better than the ETM + data quality in the visible bands, especially the near-infrared band of OLI the data, which had a clear improvement; clear improvement was not founded in the shortwave-infrared bands. Moreover, (2) OLI data had a satisfactory performance in terms of land cover classification. In summary, OLI data were a reliable data source for monitoring land cover and provided the continuity in the Landsat earth observation.     

Assessment of Landsat atmospheric correction methods for water color applications using global AERONET-OC data

With the longest archive of satellite remote sensing images, the Landsat series of satellites have demonstrated their great potential in aquatic environmental studies. However, although various atmospheric correction (AC) methods have been developed for Landsat observations in water color applications, a comprehensive assessment of their accuracies across different AC methods and instruments has yet to be performed. Using in situ spectral data collected by Aerosol Robotic Network-Ocean Color (AERONET-OC) sites, the performances of five types of AC methods over three different Landsat missions (i.e., Landsat 5/7/8) were evaluated. The Landsat 8 Operational Land Imager (OLI) showed more accurate AC retrievals than the other two instruments, and the results for its green and red bands appeared more reliable than those for the other wavelengths (uncertainty levels of ∼30 %). The iterative NIR algorithm with 2-bands (NIR-SWIR2) model selection embedded in SeaDAS showed the best performances for OLI in two blue bands. Moreover, larger residual errors were found for most Landsat 5/7 bands regardless of the AC methods and spectral bands employed with an uncertainty of >50 %. Interestingly, a simple aerosol subtraction method over the Rayleigh-corrected reflectance (Rrc) outperformed the exponential extrapolation (EXP) algorithms, especially for Landsat 5/7. Neither the image-based AC algorithm nor the surface reflectance (SR) products provided by the United States Geological Survey (USGS) showed acceptable performances over coastal environments. The uncertainties in the various Landsat reflectance products over water surfaces could be associated with a relatively poor signal-to-noise ratio (SNR) in addition to radiometric calibration uncertainties, imperfect aerosol removal methods. Future research is required to collect in situ data across a wider range of water optical properties (particularly more turbid inland waters) to examine the corresponding applicability of Landsat-series observations.     

基于随机森林的土地利用分类与景观格局分析

基于景观格局分析方法,分别选取8个反映景观格局类型和6个反映景观水平格局指数用于探索土地利用变化状况.实验结果表明:①基于选择合适的训练样本下,对3期遥感影像采用随机森林的监督分类方法,总体精度均在94％以上,kappa系数均为90％以上;②综合来看,研究期内建设用地为优势景观类型.研究区斑块数量增加,破碎化程度有所加...     

An experimental comparison between KELM and CART for crop classification using Landsat-8 OLI data

The operational land imager (OLI) is the latest instrument in the Landsat series of satellite imagery, which officially began normal operations on 30 May 2013. The OLI includes two bands that are not on the thematic mapper series of sensors aboard Landsat-5 and 7; a cirrus band and a coastal/aerosol band. This paper compares the classification and regression tree and the kernel-based extreme learning machine (KELM) for mapping crops in Hokkaido, Japan, using OLI data, except the cirrus band and the pan band. The OLI data acquired on 8 July 2013 was used for crop classification of beans, beets, grassland, maize, potatoes and winter wheat. The KELM algorithm performed better in this study and achieved overall accuracies of 90.1%. According to the Jeffries–Matusita (J–M) distances, the short wavelength infrared band provides the greater contribution (the highest value was observed for band 6 in OLI data).     

Assessing the potential of integrated Landsat 8 thermal bands,with the traditional reflective bands and derived vegetation indices in classifying urban landscapes

Reliable and up-to-date urban land cover information is valuable in urban planning and policy development. Due to the increasing demand for reliable land cover information there has been a growing need for robust methods and datasets to improve the classification accuracy from remotely sensed imagery. This study sought to assess the potential of the newly launched Landsat 8 sensor’s thermal bands and derived vegetation indices in improving land cover classification in a complex urban landscape using the support vector machine classifier. This study compared the individual and combined performance of Landsat 8’s reflective, thermal bands and vegetation indices in classifying urban land use-land cover. The integration of Landsat 8 reflective bands, derived vegetation indices and thermal bands overall produced significantly higher accuracy classification results than using traditional bands as standalone (i.e. overall, user and producer accuracies). An overall accuracy above 89.33% and a kappa index of 0.86, significantly higher than the one obtained with the use of the traditional reflective bands as a standalone data-set and other analysis stages. On average, the results also indicate high producer and user accuracies (i.e. above 80%) for most of the classes with a McNemar’s Z score of 9.00 at 95% confidence interval showing significant improvement compared with classification using reflective bands as standalone. Overall, the results of this study indicate that the integration of the Landsat 8’s OLI and TIR data presents an invaluable potential for accurate and robust land cover classification in a complex urban landscape, especially in areas where the availability of high resolution datasets remains a challenge.     

Detection of coffee berry necrosis by digital image processing of landsat 8 oli satellite imagery

Coffee berry necrosis is a fungal disease that, at a high level, significantly affects coffee productivity. With the advent of surface mapping satellites, it was possible to obtain information about the spectral signature of the crop on a time scale pertinent to the monitoring and detection of plant phenological changes. The objective of this paper was to define the best machine learning algorithm that is able to classify the incidence CBN as a function of Landsat 8 OLI images in different atmospheric correction methods. Landsat 8 OLI images were acquired at the dates closest to sampling anthracnose field data at three times corresponding to grain filling period and were submitted to atmospheric corrections by DOS, ATCOR, and 6SV methods. The images classified by the algorithms of machine learning, Random Forest, Multilayer Perceptron and Naive Bayes were tested 30 times in random sampling. Given the overall accuracy of each test, the algorithms were evaluated using the Friedman and Nemenyi tests to identify the statistical difference in the treatments. The obtained results indicated that the overall accuracy and the balanced accuracy index were on an average around 0.55 and 0.45, respectively, for the Naive Bayes and Multilayer Perceptron algorithms in the ATCOR atmospheric correction. According to the Friedman and Nemenyi tests, both algorithms were defined as the best classifiers. These results demonstrate that Landsat 8 OLI images were able to identify an incidence of the coffee berry necrosis by means of machine learning techniques, a fact that cannot be observed by the Pearson correlation.     

一种拓展的半物理时空融合算法及其初步应用

Landsat 5卫星较低的时间分辨率(16天)使得其很难获得大区域的、时相一致的清晰影像数据集。本文发展了一种基于半物理模型的时空融合算法-即乘性调制融合算法,并借助多时序的MODIS反射率数据来生成多时相的Landsat TM/ETM+反射率合成影像,经镶嵌后得到区域尺度的高时空分辨率地表反射率数据集(Landsat TM/ETM+)。本文利用吉林省2006年—2011年的Landsat 5 TM地表反射率数据以及500 m的MOD09A1反射率产品来生成3个时相的Landsat 5 TM反射率合成数据,从而获得研究区在上述时相下地表反射率数据的镶嵌图。初步分析表明,所生成的Landsat 5 TM反射率数据的光谱分布特征与MOD09A1反射率数据较为一致,且图像在整体上光谱特征的连续性较好。     

Object-based crop classification using multi-temporal SPOT-5 imagery and textural features with a Random Forest classifier

In this study, an object-based image analysis (OBIA) approach was developed to classify field crops using multi-temporal SPOT-5 images with a random forest (RF) classifier. A wide range of features, including the spectral reflectance, vegetation indices (VIs), textural features based on the grey-level co-occurrence matrix (GLCM) and textural features based on geostatistical semivariogram (GST) were extracted for classification, and their performance was evaluated with the RF variable importance measures. Results showed that the best segmentation quality was achieved using the SPOT image acquired in September, with a scale parameter of 40. The spectral reflectance and the GST had a stronger contribution to crop classification than the VIs and GLCM textures. A subset of 60 features was selected using the RF-based feature selection (FS) method, and in this subset, the near-infrared reflectance and the image acquired in August (jointing and heading stages) were found to be the best for crop classification.     

应用时间序列EVI的MERSI多光谱混合像元分解

针对风云3数据的特点,本文将EVI生长曲线引入多光谱混合像元的分解。首先,利用Landsat8 OLI影像,采用支持向量机的分类方法,提取研究区域的耕地信息,利用该信息对风云MERSI数据进行掩膜处理,获得研究区域的耕地影像。接着,利用MERSI时序影像,计算像元EVI值,通过SG滤波,构建农作物(端元)和混合像元的EVI生长曲线。通过实地调查,获取研究区的农作物端元,尤其对主要的农作物玉米,在空间上均匀选取了14个端元。然后,采用传统的方法,将14种玉米端元生长曲线分别与其它端元组合,进行混合像元分解。发现分解的效果差异很大,提取的玉米种植面积从191.90 km2到574.83 km2不等。为提高分解精度,借用光谱匹配(光谱夹角最小)的方法(用生长曲线代替光谱曲线)自适应选择与混合像元EVI曲线最相似的玉米端元作为组合端元,进行混合像元分解。结果得到玉米的种植面积为589.95 km2,比传统方法的最好(相对)精度提高了2%。     

北京地区Landsat 8 OLI高空间分辨率气溶胶光学厚度反演

卫星气溶胶光学厚度(AOD)反演中,传统暗目标方法在反射率较低的水体、浓密植被覆盖区域取得了较好效果,在反射率较高且结构复杂的高反射地表上空目前多采用深蓝算法,但存在空间分辨率较低,对细节分布描述性较差等问题。为解决这一问题,本文首先以5年(2008年—2012年)长时间序列MODIS地表反射率产品为基础,采用最小值合成法建立500 m分辨率逐月地表反射率产品数据集,然后利用地物波谱库中典型地物波谱数据,分析建立MODIS与Landsat 8 OLI传感器蓝光波段反射率转换模型,最后北京地区AERONET地基观测数据确定了气溶胶光学物理参数,并反演获取了北京地区上空500 m分辨率的AOD分布。为验证反演算法的精度,分别将反演结果同AERONET及MODIS/Terra气溶胶产品(MOD04)进行交叉对比,同时利用相关系数R,均方根误差RMSE,平均绝对误差MAE以及MODIS AOD产品预期误差EE共4个指标进行衡量。结果表明:算法反演获取的AOD与AERONET观测值具有较高的一致性,各指标分别为R=0.963,RMSE=0.156,MAE=0.097,EE=85.3%,稍优于MOD04产品(R=0.962,RMSE=0.158,MAE=0.101,EE=75.8%),并且有效的对比点数也高于MOD04。通过与地基观测相比,卫星遥感获取的高分辨率城市地区AOD精度可作为定量评估城市空气质量的有效依据。     

基于遥感的区域尺度森林地上生物量估算研究

森林是陆地生态系统最大的碳库,精确估算森林生物量是陆地碳循环研究的关键。首先从机载LiDAR数据中提取高度和密度统计量,采用逐步回归模型进行典型样区生物量估算;然后利用机载LiDAR数据估算的生物量作为样本数据,与多光谱遥感数据Landsat8 OLI的波段反射率及植被指数建立回归模型,实现区域尺度森林地上生物量估算。实验结果显示,机载LiDAR数据估算的鼎湖山样区生物量与地面实测生物量的相关性R2达0.81,生物量RMSE为40.85 t/ha,说明机载LiDAR点云数据的高度和密度统计量与生物量存在较高的相关性。以机载LiDAR数据估算的生物量为样本数据,结合多光谱遥感数据Landsat8 OLI估算粤西北地区的森林地上生物量,精度验证结果为:R2为0.58,RMSE为36.9 t/ha;针叶林、阔叶林和针阔叶混交林等3种不同森林类型生物量的估算结果为:R2分别为0.51(n=251)、0.58(n=235)和0.56(n=241),生物量RMSE分别为24.1 t/ha、31.3 t/ha和29.9 t/ha,估算精度相差不大。总体上看,利用遥感数据可以开展区域尺度的森林地上生物量估算,为森林固碳监测提供有力的参考数据。     

Performance evaluation of object based greenhouse detection from Sentinel-2 MSI and Landsat 8 OLI data: A case study from Almería (Spain)

This paper shows the first comparison between data from Sentinel-2 (S2) Multi Spectral Instrument (MSI) and Landsat 8 (L8) Operational Land Imager (OLI) headed up to greenhouse detection. Two closely related in time scenes, one for each sensor, were classified by using Object Based Image Analysis and Random Forest (RF). The RF input consisted of several object-based features computed from spectral bands and including mean values, spectral indices and textural features. S2 and L8 data comparisons were also extended using a common segmentation dataset extracted form VHR World-View 2 (WV2) imagery to test differences only due to their specific spectral contribution. The best band combinations to perform segmentation were found through a modified version of the Euclidian Distance 2 index. Four different RF classifications schemes were considered achieving 89.1%, 91.3%, 90.9% and 93.4% as the best overall accuracies respectively, evaluated over the whole study area.     

Improved coastal wetland mapping using very-high 2-meter spatial resolution imagery

Accurate wetland maps are a fundamental requirement for land use management and for wetland restoration planning. Several wetland map products are available today; most of them based on remote sensing images, but their different data sources and mapping methods lead to substantially different estimations of wetland location and extent. We used two very high-resolution (2 m) WorldView-2 satellite images and one (30 m) Landsat 8 Operational Land Imager (OLI) image to assess wetland coverage in two coastal areas of Tampa Bay (Florida): Fort De Soto State Park and Weedon Island Preserve. An initial unsupervised classification derived from WorldView-2 was more accurate at identifying wetlands based on ground truth data collected in the field than the classification derived from Landsat 8 OLI (82% vs. 46% accuracy). The WorldView-2 data was then used to define the parameters of a simple and efficient decision tree with four nodes for a more exacting classification. The criteria for the decision tree were derived by extracting radiance spectra at 1500 separate pixels from the WorldView-2 data within field-validated regions. Results for both study areas showed high accuracy in both wetland (82% at Fort De Soto State Park, and 94% at Weedon Island Preserve) and non-wetland vegetation classes (90% and 83%, respectively). Historical, published land-use maps overestimate wetland surface cover by factors of 2–10 in the study areas. The proposed methods improve speed and efficiency of wetland map production, allow semi-annual monitoring through repeat satellite passes, and improve the accuracy and precision with which wetlands are identified.     

Use of landsat ETM+ SLC-off segment-based gap-filled imagery for crop type mapping

Failure of the Scan Line Corrector (SLC) on the Landsat ETM+ sensor has had a major impact on many applications that rely on continuous medium resolution imagery to meet their objectives. The United States Department of Agriculture (USDA) Cropland Data Layer (CDL) program uses Landsat imagery as the primary source of data to produce crop-specific maps for 20 states in the USA. A new method has been developed to fill the image gaps resulting from the SLC failure to support the needs of Landsat users who require coincident spectral data, such as for crop type mapping and monitoring. We tested the new gap-filled method for a CDL crop type mapping project in eastern Nebraska. Scan line gaps were simulated on two Landsat 5 images (spring and late summer 2003) and then gap-filled using landscape boundary models, or segment models, that were derived from 1992 and 2002 Landsat images (used in the gap-fill process). Various date combinations of original and gap-filled images were used to derive crop maps using a supervised classification process. Overall kappa values were slightly higher for crop maps derived from SLC-off gap-filled images compared to crop maps derived from the original imagery (0.3–1.3% higher). Although the age of the segment model used to derive the SLC-off gap-filled product did not negatively impact the overall agreement, differences in individual cover type agreement did increase (?0.8%–1.6% using the 2002 segment model to ?5.0–5.1% using the 1992 segment model). Classification agreement also decreased for most of the classes as the size of the segment used in the gap-fill process increased.     

Crop classification and acreage estimation in North Korea using phenology features

In North Korea, reliable and timely information on crop acreage and spatial distribution is hard to obtain. In this study, we developed a fast and robust method to estimate crop acreage in North Korea using time-series normalized difference vegetation index (NDVI) derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) data. We proposed a method to identify crop type based on NDVI phenology features using data collected in other areas with similar agri-environmental conditions to mitigate the shortage of ground truth data. Eventually the classification map (MODIScrop) was assessed using the Food and Agriculture Organization (FAO) statistical data and high-resolution crop classification maps derived from one Landsat scene (LScrop). The Pareto boundary method was used to assess the accuracy and crop distribution of the MODIScrop maps. Results showed that acreage derived from the MODIScrop maps was generally consistent with that reported in the FAO data (a relative error <4.1% for rice and <6.1% for maize, and <9.0% for soybean except for in 2004, 2008, and 2009) and the maps derived from the LScrop (a relative error about 5% in 2013, and 7% in 2008 and 2014). The classification accuracy reached 74.4%, 69.8%, and 73.1% of the areas covered by the Landsat images in 2008, 2013, and 2014, respectively. This indicates that features derived from NDVI profiles were able to characterize major crops, and the approaches developed in this study are feasible for crop mapping and acreage estimation in regions with limited ground truth data.