多时相MODIS影像水田信息提取研究

水稻种植及其分布信息是土地覆被变化、作物估产、甲烷排放、粮食安全和水资源管理分析的重要数据源。基于遥感的水田利用监测中,通常采用时序NDVI植被指数法和影像分类法分别进行AVHRR和TM影像的水田信息获取。针对8天合成MODIS陆地表面反射比数据的特点和水稻生长特征,选取水稻种植前的休耕期、秧苗移植期、秧苗生长期和成熟期等多时相MODIS地表反射率影像数据,通过归一化植被指数、增强植被指数及利用对土壤湿度和植被水分含量较敏感的短波红外波段计算得到的陆表水指数进行水田信息获取。将提取结果与基于ETM＋影像的国土资源调查水田数据,通过网格化计算处理并进行对比分析,结果表明,利用MODIS影像的8天合成地表反射率数据,进行区域甚至全国的水田利用监测是可行的。     

融合多源时序遥感数据大尺度不透水面覆盖率估算

不透水面信息是监测城市扩张及区域生态环境变化研究的重要指标,基于遥感技术对地表不透水面信息进行快速提取具有重要意义。传统大范围不透水面覆盖率估算模型主要基于单一遥感信息与不透水面比例之间的相关性,通过单因子回归模型实现不透水面覆盖率的估算。受限于单一遥感信息的信息量及普适性等影响,这类方法在大尺度不透水面提取中具有较大局限性,估算结果的区域适应性存在较大差异。针对该问题,本文提出基于多特征遥感信息进行不透水面估算的方法,以弥补单一特征在大范围不透水面提取中的不确定性。该方法首先以多时相MOD13Q1、MOD09A1产品、夜间灯光数据（NPP-VIIRS）和Landsat 8 OLI为遥感数据源,从不同角度构建突出不透水面信息的多个指数特征;在此基础上利用多元回归模型建立多因子不透水面覆盖率估算模型,进而实现大尺度不透水面覆盖率的遥感估算。本研究选择分布于全国范围内13个典型城市作为主要研究区对提出的模型进行了验证,结果表明：该方法能够适应不同区域不透水面覆盖率的估算,在复杂城市区域表现出较传统方法更好的效果,明显改善了城市内部不透水面覆盖率的估算精度。     

像素工厂中真正射影像的生产及其精度分析

像素工厂是当今世界一流的遥感影像自动化处理系统,集自动化、并行处理、多种影像兼容、远程管理等特点于一身,代表了当前遥感影像数据处理技术的发展方向,主要用于地形图测绘、城市规划等。这里主要分析像素工厂生产数字表面模型(DSM)和真正射影像的工作流程,并结合实际生产进行精度分析。     

基于多角度遥感物理模型的叶面积指数反演

与传统遥感观测相比,多角度对地观测通过对林木多个方向的观察,可得到丰富的森林三维空结构信息,为定量遥感提供新的途径。本文根据不同遥感数据选择合适的模型并建立相应的查找表。TM/ETM+数据采用混合像元分解模型,CHRIS数据考虑不同森林场景选择不同的遥感物理模型,在选择合适的模型基础上,根据模型的不同敏感参数和试验区选择的特点设计查找表参数,并由模型正演建立查找表,再根据插值的方法由遥感图像的反射率值反演LAI。两种数据反演得到的LAI与MODIS 15A2 LAI产品比较分析。结果表明多角度遥感反演准确率有一定的提高,具有一定的理论研究和实际应用价值。     

Assessing the optimal preprocessing steps of MODIS time series to map cropping systems in Mato Grosso,Brazil

The adoption of new cropping practices such as integrated Crop-Livestock systems (iCL) aims at improving the land use sustainability of the agricultural sector in the Brazilian Amazon. The emergence of such integrated systems, based on crop and pasture rotations over and within years, challenges the remote sensing community who needs to implement accurate and efficient methods to process satellite image time series (SITS) in order to come up with a monitoring protocol. These methods generally include a SITS preprocessing step which can be time consuming. The aim of this study is to assess the importance of preprocessing operations such as temporal smoothing and computation of phenological metrics on the mapping of main cropping systems (i.e. pasture, single cropping, double cropping and iCL), with a special emphasis on the iCL class. The study area is located in the state of Mato Grosso, an important producer of agriculture commodities located in the Southern Brazilian Amazon. SITS were composed of a set of 16-day composites of MODIS Vegetation Indices (MOD13Q1 product) covering a one year period between 2014 and 2015. Two widely used classifiers, i.e. Random Forest (RF) and Support Vector Machine (SVM), were tested using five data sets issued from a same SITS but with different preprocessing levels: (i) raw NDVI; (ii) raw NDVI + raw EVI; (iii) smoothed NDVI; (iv) NDVI-derived phenometrics; (v) raw NDVI + phenometrics. Both RF and SVM classification results showed that the “raw NDVI + raw EVI” data set achieved the highest performance (RF OA = 0.96, RF Kappa = 0.94, SVM OA = 0.95, SVM Kappa = 0.93), followed closely by the “raw NDVI” and the “raw NDVI + phenometrics” datasets. The “NDVI-derived phenometrics” alone achieved the lowest accuracies (RF OA = 0.58 and SVM OA = 0.66). Considering that the implementation of preprocessing steps is computationally expensive and does not provide significant gains in terms of classification accuracy, we recommend to use raw vegetation indices for mapping cropping practices in Mato Grosso, including the integrated Crop-Livestock systems.     

利用多光谱卫星遥感和深度学习方法进行青藏高原积雪判识

青藏高原积雪对全球气候变化十分重要,针对已有积雪遥感判识方法中普遍采用的可见光与红外光谱数据易受复杂地形与高海拔影响,导致青藏高原地区积雪判识精度较低的问题,提出了一种基于多光谱遥感与地理信息数据特征级融合的积雪遥感判识方法:以风云三号卫星可见光与红外多光谱遥感资料与多要素地理信息作为数据源,由地面实测雪深数据与现有积雪产品交叉筛选出样本标签,构建并训练基于层叠去噪自编码器(SDAE)的特征融合与分类网络,从而有效辨识青藏高原遥感图像中的云、积雪以及无雪地表。经地面实测雪深数据验证,该方法分类精度显著高于使用相同数据源的FY-3A/MULSS积雪产品,略高于国际主流积雪产品MOD10A1与MYD10A1,并且年均云覆盖率最低。试验结果表明该方法可有效地减少云层对积雪判识的干扰,提升分类精度。     

无人机平台的沙漠场地BRDF特征观测方法及建模

目标双向反射分布函数BRDF（Bidirectional Reflectance Distribution Function）不仅是陆面遥感的关键地球物理参数,也是星载光学遥感仪器基于地面目标的场地辐射校正重要参量,是影响定标精度的关键要素。传统野外地物多角度测量使用的观测设备,一般其结构较为复杂,重量体积较大,而且运输和组装过程繁琐,观测目标时容易受地形和交通限制,难以进行高效快速精确的野外测量。近年来,无人机由于其设备操作简便、运输和观测方式灵活等方面的优点,可作为新的观测平台应用于当前遥感试验中。本文设计了一种基于无人机平台的地表BRDF测量装置、观测方案和数据处理流程。利用多旋翼低空无人机和云台的组合,搭载野外地物光谱仪和跟拍相机,通过对地面目标多角度观测和高精度定位及角度控制,实现针对固定目标的多方位角和天顶角观测。本文采用上述设计方案和观测流程,在敦煌辐射校正场开展多次稳定均匀沙漠目标的多角度光谱观测试验,并利用实验观测数据,基于Ross-Li核驱动模型推算了场地BRDF模型参数,并与MODIS的陆表BRDF产品（MCD43C1）及反射率产品（MOD/MYD09）进行对比验证。通过开展野外实验,核验了这种新的BRDF观测手段的可靠性,获取的敦煌地表BRDF参数与MODIS遥感产品有良好的一致性,各波段的相对偏差在5%以内。本研究表明,基于多旋翼无人机的BRDF观测系统,提供了一种全新的地物目标方向反射特性观测方法,可用于自动化高频次场地特性观测以及卫星同步定标等野外实验活动。在保证观测精度的同时,极大地减轻人力物力的投入,值得广泛推广应用。     

基于无人机遥感的信息提取研究

无人机遥感系统可以获取矿区高分辨率遥感影像,具有较强的灵活性和现势性,成图速度快,能广泛应用于各种调查任务。本文以江西省定南县稀土矿区为研究区,介绍了无人机遥感系统的组成、航拍流程和数据处理的方法,通过面向对象的多尺度分割技术,建立不同地物类别对应的分类层次,然后根据地物特征进行信息提取,最后总结出一套基于无人机遥感的矿产信息提取方法。     

利用MODIS数据识别水稻关键生长发育期

利用遥感方法提取中国范围内的水稻关键生长发育期。首先, 对时间序列Terra MODIS-EVI(Enhanced Vegetation Index)进行傅里叶和小波低通滤波平滑处理, 然后, 根据水稻在移栽期、分蘖初期、抽穗期和成熟期的EVI变化特征, 实现对各个生长发育期的识别。通过将利用2005年MODIS数据识别的结果与当年气象台站的地面观测资料进行比较, 采用本研究中的识别方法得出的水稻各个生长发育期的绝对误差大部分小于16d, 经过F检验表明提取的结果与地面观测资料在0.05水平下具有显著一致性。研究中的信息提取方法可被用于其他年份的水稻生长发育期识别, 根据其他作物的生长发育特点, 也可能适合于提取其他作物的生长发育期。     

Comparison of four UAV georeferencing methods for environmental monitoring purposes focusing on the combined use with airborne and satellite remote sensing platforms

This work is aimed at the environmental remote sensing community that uses UAV optical frame imagery in combination with airborne and satellite data. Taking into account the economic costs involved and the time investment, we evaluated the fit-for-purpose accuracy of four positioning methods of UAV-acquired imagery: 1) direct georeferencing using the onboard raw GNSS (GNSSNAV) data, 2) direct georeferencing using Post-Processed Kinematic single-frequency carrier-phase without in situ ground support (PPK1), 3) direct georeferencing using Post-Processed Kinematic double-frequency carrier-phase GNSS data with in situ ground support (PPK2), and 4) indirect georeferencing using Ground Control Points (GCP). We tested a multispectral sensor and an RGB sensor, onboard multicopter platforms. Orthophotomosaics at <0.05 m spatial resolution were generated with photogrammetric software. The UAV image absolute accuracy was evaluated according to the ASPRS standards, wherein we used a set of GCPs as reference coordinates, which we surveyed with a differential GNSS static receiver. The raw onboard GNSSNAV solution yielded horizontal (radial) accuracies of RMSEr≤1.062 m and vertical accuracies of RMSEz≤4.209 m; PPK1 solution gave decimetric accuracies of RMSEr≤0.256 m and RMSEz≤0.238 m; PPK2 solution, gave centimetric accuracies of RMSEr≤0.036 m and RMSEz≤0.036 m. These results were further improved by using the GCP solution, which yielded accuracies of RMSEr≤0.023 m and RMSEz≤0.030 m. GNSSNAV solution is a fast and low-cost option that is useful for UAV imagery in combination with remote sensing products, such as Sentinel-2 satellite data. PPK1, which can register UAV imagery with remote sensing products up to 0.25 m pixel size, as WorldView-like satellite imagery, airborne lidar or orthoimagery, has a higher economic cost than the GNSSNAV solution. PPK2 is an acceptable option for registering remote sensing products of up to 0.05 m pixel size, as with other UAV images. Moreover, PPK2 can obtain accuracies that are approximate to the usual UAV pixel size (e.g. co-register in multitemporal studies), but it is more expensive than PPK1. Although indirect georeferencing can obtain the highest accuracy, it is nevertheless a time-consuming task, particularly if many GCPs have to be placed. The paper also provides the approximate cost of each solution.     

黑河流域遥感物候产品验证与分析

植被物候遥感产品对全球变化响应、农业生产管理、生态学的应用等多领域研究具有重要意义。但现有植被物候遥感产品还有较多问题,主要包括一方面使用不同参数的时间序列数据以及不同提取算法导致的产品结果差异较大,另一方面在地面验证中地面观测数据与遥感反演数据的物理含义不一致导致的验证方法的系统性误差。本文以黑河流域为研究区,对比验证基于EVI(Enhanced Vegetation Index)时间序列数据提取的MLCD(MODIS global land cover dynamics product)植被遥感物候产品和基于LAI(Leaf Area Index)时间序列数据提取的UMPM(product by universal multi-life-cycle phenology monitoring method)植被遥感物候产品的有效性及精度等。同时,通过验证分析进一步评估基于EVI和LAI时间序列提取的物候特征的差异及特点,探讨由于地面观测植被物候与遥感提取植被物候的物理意义的不一致问题导致的直接验证结果偏差。结果表明:UMPM产品有效性整体高于MLCD产品,但在以草地和灌木为主的稀疏植被区,由于LAI取值精度的原因,UMPM产品存在较多缺失数据,且时空稳定性较低;基于玉米地面观测数据表明,EVI对植被开始生长的信号比LAI更加敏感,更适合提取生长起点,但植被指数易饱和,峰值起点普遍提前,基于LAI提取的峰值起点更加合理。由于地面观测的物候期在后期更加关注果实生长,遥感观测仅关注叶片的生长,遥感定义的峰值终点和生长终点与玉米的乳熟期和成熟期差异较大。     

Object-based urban land cover classification using rule inheritance over very high-resolution multisensor and multitemporal data

Very high spatial and temporal resolution remote sensing data facilitate mapping highly complex and diverse urban environments. This study analyzed and demonstrated the usefulness of combined high-resolution aerial digital images and elevation data, and its processing using object-based image analysis for mapping urban land covers and quantifying buildings. It is observed that mapping heterogeneous features across large urban areas is time consuming and challenging. This study presents and demonstrates an approach for formulating an optimal land cover classification rule set over small representative training urban area image, and its subsequent transfer to the multisensor, multitemporal images. The classification results over the training area showed an overall accuracy of 96%, and the application of rule set to different sensor images of other test areas resulted in reduced accuracies of 91% for the same sensor, 90% and 86% for the different sensors temporal data. The comparison of reference and classified buildings showed ±4% detection errors. Classification through a transferred rule set reduced the classification accuracy by about 5%–10%. However, the trade-off for this accuracy drop was about a 75% reduction in processing time for performing classification in the training area. The factors influencing the classification accuracies were mainly the shadow and temporal changes in the class characteristics.     

基于QR码的瓦片遥感影像数字水印算法

针对瓦片遥感影像数据,提出一种基于快速响应矩阵码（quick response code,QR code）的强鲁棒性盲水印算法。算法的设计充分顾及了瓦片遥感影像的数据特征,将包含版权认证信息的QR码水印图像进行数据压缩和预处理,有效降低了嵌入数据量;通过在水印嵌入单元和水印信息位之间建立稳定的映射关系确保水印信息的可靠嵌入;采用改变像素值位平面统计特征的方法实现对含水印瓦片遥感影像的快速定位。最后,对所提算法进行了试验验证。结果表明,QR码的引入有效提升了水印信息容量,所提算法对于常规图像处理操作和拼接攻击具备良好的保真度和鲁棒性,可实际应用于瓦片遥感影像数据的版权保护。     

结合像元级和目标级的高分辨率遥感影像建筑物变化检测

快速、精准的建筑物变化检测对城市规划建设等业务管理具有重要意义。随着卫星遥感技术的快速发展,基于高分辨率遥感影像的建筑物变化检测得到了广泛关注。针对像元级建筑物变化检测方法往往精度不足而目标级建筑物变化检测方法过程烦琐等问题,本文提出结合像元级和目标级的高分辨率遥感影像建筑物变化检测方法。首先综合高分辨率遥感影像的多维特征,利用随机森林分类器进行影像集分类,以获取像元级建筑物变化检测结果;然后对后时相遥感影像进行图像分割,获得影像对象;最后融合像元级建筑物变化检测结果和影像对象,识别变化的建筑物目标。利用双时相QuickBird高分辨率遥感影像进行建筑物变化检测试验,结果表明：本文提出的方法能够削弱光照、观测角度等环境差异对建筑物变化检测的影响,显著改善建筑物变化的检测精度。     

Comparison of eMODIS and MOD/MYD13A2 NDVI products during 2012–2014 spring green-up periods in Alaska and northwest Canada

Accurate monitoring of vegetation dynamics is required to understand the inter-annual variability and long term trends in terrestrial carbon exchange in tundra and boreal ecoregions. In western North America, two Normalized Vegetation Index (NDVI) products based on spectral reflectance data from the Moderate Resolution Imaging Spectroradiometer (MODIS) are available. The MOD/MYD13A2 NDVI product is available as a 16-day composite product in a sinusoidal projection as global hdf tiles. The eMODIS Alaska NDVI product is available as a 7-day composite geotif product in a regional equal area conic projection covering Alaska and the entire Yukon River Basin. These two NDVI products were compared for the 2012–2014 late May–late June spring green-up periods in Alaska and the Yukon Territory. Relative to the MOD/MYD13A2 NDVI product, it is likely that the eMODIS NDVI product contained more cloud-contaminated NDVI values. For example, the MOD/MYD13A2 product flagged substantially fewer pixels as “good quality” in each 16-day composite period compared to the corresponding MODIS Alaska NDVI product from a 7-day composite period. During the spring green-up period, when field-based NDVI increases, the eMODIS NDVI product averaged 43 percent of pixels that declined by at least 0.05 NDVI between 2 composite periods, consistent with cloud-contamination problems, while the MOD/MYD13A2 NDVI averaged only 6 percent of pixels. Based on a cloudy Landsat-8 scene, the eMODIS compositing process selected 23 percent pixels, while the MOD/MYD13A2 compositing process selected less than 0.003 percent pixels. Based on the results, it appears that the MOD/MYD13A2 NDVI product is superior for scientific applications based on NDVI phenology in the tundra and boreal regions of northwestern North America.     

Hyperspectral Satellite Data in Mapping Salt-Affected Soils Using Linear Spectral Unmixing Analysis

Development of salt-affected soils in the irrigated lands of arid and semi-arid region is major cause of land degradation. Hyperion hyperspectral remote sensing data (EO-1) was used in the present study for characterization and mapping of salt-affected soils in a part of irrigation command area of Indo-Gangetic alluvial plains. Linear spectral mixture analysis approach was used to map various categories of salt affected soils represented by spectral endmembers of slightly, moderately and highly salt-affected soils. These endmembers were related to surface expression of various categories of salt-affected soils in the area. The endmembers were selected by performing minimum noise fraction (MNF) transformation and pixel purity index (PPI) on Hyperion (EO-1) data with reference to high resolution LISS IV data and field data. The results showed that various severity classes of salt-affected soils could be reliably mapped using linear spectral unmixing analysis. A low RMSE value (0.0193) over the image was obtained that revealed a good fit of the model in identification and classification of endmembers of various severities of salt affected soils. The overall classification accuracies for slight, moderate and highly salt-affected soils were estimated of 78.57, 79.81 and 84.43% respectively.     

Geographical model for precise agriculture monitoring with real-time remote sensing

This paper presents a remote sensing model for crop monitoring that was developed by the authors in a multi-year study. It also presents two experiments conducted for testing a newly developed application. The model combines remote sensing models using mapping of the spatial distribution of vegetation in an agricultural field, with precision agricultural models that maximize the output (yield) while minimizing the input (cost). This combination enables one to operate a monitoring and management process that includes every sub-unit of the field using remote sensing mapping.The model consists of five steps: (1) Preparing information layers that map the crop-affecting elements, e.g. irrigation and topography; (2) Collecting spectral and plant data simultaneously; (3) Processing and analyzing the data in order to prepare vegetation maps; (4) Decision-making in accordance with the above-mentioned maps or with predicted-yield maps; and (5) Quality control.The experiments showed that although the results were not statistically significant, the application of the proposed model enables one to draw recommendations within 45 h, and that remote sensing monitoring results in more benefits than do traditional control methods. The quality control was not ideal, due to the narrow range of the spectrum used in the remote sensing monitoring.     

光学遥感影像智能化处理研究进展

随着对地观测技术的发展,海量遥感影像不断传输到地面。传统的遥感信息处理方法在处理效率、精度上的不足,限制了遥感信息的挖掘及利用,亟需发展智能化方法满足遥感影像处理的需求。受自然界中生物进化机制的启发,基于进化计算的遥感影像智能化处理方法具有以下特点:(1)拥有全局优化能力,对目标函数的优化能力更强;(2)具有自组织、自学习的特点,能够从遥感数据本身学习,不依赖数据分布等先验信息;(3)拥有处理多目标问题的能力,同时考虑多个目标函数而不需要人工确定它们之间的权重。因此,智能化遥感信息处理方法能够在海量遥感影像中有效地提取适用于不同应用目的的信息。本文主要介绍智能化遥感信息处理方法的典型应用包括遥感影像智能化分类(监督分类、聚类)、遥感影像亚像素信息提取(高光谱影像混合像元分解、亚像元制图),并讨论了遥感信息智能化处理方法的发展方向。     

陆地总初级生产力遥感估算精度分析

准确估算陆地总初级生产力GPP(Gross Primary Productivity)数值对碳循环过程模拟有重要影响。本文介绍了多种基于植被指数以及基于光能利用率的遥感GPP算法,综述了不同算法在其研究区域的估算精度;并分析了MODIS/GPP以及BESS/GPP两种遥感GPP产品在不同植被类型的估算精度。通过对比全球碳通量站网络GPP数据表明,MODIS/GPP产品在全球估算结果具显著相关性(R2=0.59)及中等标准误差(RMSE=2.86 g C/m2/day),估算精度较高的植被类型有落叶阔叶林,草地等;估算精度较低类型包括常绿阔叶林,稀树草原等。本文对GPP产品中存在的不确定性进行分析,通过综述前人研究中发现的遥感估算GPP方法中存在的问题,指出可能的提高卫星遥感GPP产品估算精度的方法及发展趋势。     

Evaluation of ICA and CEM algorithms with Landsat-8/ASTER data for geological mapping in inaccessible regions

Many regions remain poorly studied in terms of geological mapping and mineral exploration in inaccessible regions especially in the Arctic and Antarctica due to harsh conditions and logistic difficulties. Application of specialized image processing techniques is capable of revealing the hidden linear mixing spectra in multispectral and hyperspectral satellite images. In this study, the applications of Independent component analysis (ICA) and Constrained Energy Minimization (CEM) algorithms were evaluated for Landsat-8 and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) remote sensing data for geological mapping in Morozumi Range and Helliwell Hills areas, Northern Victoria Land (NVL), Antarctica. The results of this investigation demonstrate the capability of the two algorithms in distinguishing pixel and subpixel targets in the multispectral satellite data. The application of the methods for identifying poorly exposed geologic materials and subpixel exposures of alteration minerals has invaluable implications for geological mapping and mineral exploration in inaccessible regions.