利用BRDF原型和单方向反射率数据估算地表反照率

地表反照率是影响地表能量收支平衡的决定性参数之一,精确反演地表反照率需要考虑地表各向异性反射特征。本文尝试以双向反射分布函数BRDF原型为地表各向异性反射的先验知识,通过单方向反射率反演地表反照率。首先根据地面实测及MODIS多角度反射率数据对反演方法进行分析和精度评价,然后借助MODIS BRDF产品统计出研究区的主导BRDF原型,并联合环境一号卫星(HJ-1)单方向反射率数据反演30 m地表反照率,最终将结果与地表实测数据进行比较。结果表明:BRDF原型对BRDF的变化进行了约束,且能够适用于几十米尺度的遥感数据反照率的反演;不同级别的各向异性反射特征的分布是不均一的,借助于主导BRDF原型能够使大部分样本的地表反照率满足精度要求;利用研究区MODIS BRDF产品统计得到的主导BRDF原型为先验知识,通过HJ-1数据反演得到的地表反照率与地表实测反照率有较高的一致性,而朗伯假定条件下的反照率高于实测结果。本文算法简单高效,可为产生全国范围的中高分辨卫星反照率产品提供有价值的算法参考。     

Evaluation of directional normalization methods for Landsat TM/ETM+ over primary Amazonian lowland forests

Biodiversity mapping in extensive tropical forest areas poses a major challenge for the interpretation of Landsat images, because floristically clearly distinct forest types may show little difference in reflectance. In such cases, the effects of the bidirectional reflection distribution function (BRDF) can be sufficiently strong to cause erroneous image interpretation and classification. Since the opening of the Landsat archive in 2008, several BRDF normalization methods for Landsat have been developed. The simplest of these consist of an empirical view angle normalization, whereas more complex approaches apply the semi-empirical Ross–Li BRDF model and the MODIS MCD43-series of products to normalize directional Landsat reflectance to standard view and solar angles. Here we quantify the effect of surface anisotropy on Landsat TM/ETM+ images over old-growth Amazonian forests, and evaluate five angular normalization approaches. Even for the narrow swath of the Landsat sensors, we observed directional effects in all spectral bands. Those normalization methods that are based on removing the surface reflectance gradient as observed in each image were adequate to normalize TM/ETM+ imagery to nadir viewing, but were less suitable for multitemporal analysis when the solar vector varied strongly among images. Approaches based on the MODIS BRDF model parameters successfully reduced directional effects in the visible bands, but removed only half of the systematic errors in the infrared bands. The best results were obtained when the semi-empirical BRDF model was calibrated using pairs of Landsat observation. This method produces a single set of BRDF parameters, which can then be used to operationally normalize Landsat TM/ETM+ imagery over Amazonian forests to nadir viewing and a standard solar configuration.     

Consistency of vegetation index seasonality across the Amazon rainforest

Vegetation indices (VIs) calculated from remotely sensed reflectance are widely used tools for characterizing the extent and status of vegetated areas. Recently, however, their capability to monitor the Amazon forest phenology has been intensely scrutinized. In this study, we analyze the consistency of VIs seasonal patterns obtained from two MODIS products: the Collection 5 BRDF product (MCD43) and the Multi-Angle Implementation of Atmospheric Correction algorithm (MAIAC). The spatio-temporal patterns of the VIs were also compared with field measured leaf litterfall, gross ecosystem productivity and active microwave data. Our results show that significant seasonal patterns are observed in all VIs after the removal of view-illumination effects and cloud contamination. However, we demonstrate inconsistencies in the characteristics of seasonal patterns between different VIs and MODIS products. We demonstrate that differences in the original reflectance band values form a major source of discrepancy between MODIS VI products. The MAIAC atmospheric correction algorithm significantly reduces noise signals in the red and blue bands. Another important source of discrepancy is caused by differences in the availability of clear-sky data, as the MAIAC product allows increased availability of valid pixels in the equatorial Amazon. Finally, differences in VIs seasonal patterns were also caused by MODIS collection 5 calibration degradation. The correlation of remote sensing and field data also varied spatially, leading to different temporal offsets between VIs, active microwave and field measured data. We conclude that recent improvements in the MAIAC product have led to changes in the characteristics of spatio-temporal patterns of VIs seasonality across the Amazon forest, when compared to the MCD43 product. Nevertheless, despite improved quality and reduced uncertainties in the MAIAC product, a robust biophysical interpretation of VIs seasonality is still missing.     

MISR和MODIS二向性反射数据产品的对比分析

利用4种地表类型共16个地面站点的MISR 9个角度反射率数据产品和MODIS反射率数据产品(MOD09A1), 对比两种传感器获取的地表方向反射观测数据和用BRDF模型参数产品计算的方向反射数据, 分析了这两种BRDF模型参数的外推能力。研究结果表明: (1) MISR和MODIS BRDF模型参数数据产品都具有一定的对观测数据以外方向反射的外推能力, 相比用MISR BRDF模型参数数据外推到相应MODIS观测角度方向反射的结果一致性较强。(2) 采用两种BRDF模型参数推算的方向反射率与观测数据的接近程度有随观测天顶角的增大而减弱的趋势, 在观测天顶角较小时一致性较好。(3) 所用数据处理结果显示, MODIS BRDF模型参数数据产品在近垂直主平面方向反射率推算结果相比近主平面方向的推算结果更接近观测数据。     

利用POLDER数据验证MODIS BRDF模型参数产品及Ross-Li模型

将MODIS BRDF模型参数产品(MCD43A1)模拟的近红外波段反射率与从POLDER-3/PARASOL BRDF全球数据集中筛选的9961个像元的BRDF观测数据进行对比,验证了MCD43A1所采用的RossThick-LiSparseR BRDF 模型(Ross-Li模型)拟合二向反射的能力。结果表明,Ross-Li模型总体上可以有效地模拟地物的二向反射,所有像元的近红外波段反射率模拟值与POLDER-3观测数据之间的R²达到0.943,RMSE为0.016,模拟反射率比POLDER-3数据总体偏低5.2%。但Ross-Li模型明显低估了热点反射率,热点模拟结果比POLDER-3数据平均偏低14%,模拟值的R2为0.824,RMSE为0.07。热点反射率模拟误差与地表覆盖类型有关,针叶林热点反射率模拟值偏低最严重,其次是阔叶林、草地与农田,灌木与裸地热点反射率模拟值偏低相对较小。通过修正Ross-Li模型中的体散射核,可以明显改善热点反射率的模拟效果(R²＝0.839,RMSE＝0.043)。Ross-Li模型对天底、暗点等特征方向反射率的模拟较为准确。Ross-Li模型的模拟精度随太阳天顶角和观测天顶角的增大而降低。对于农田与草地而言,Ross-Li模型的模拟精度随NDVI的增加而降低;但在森林与灌木覆盖条件下,当NDVI约为0.5时,Ross-Li模型的模拟效果最差。     

MODIS二向反射模型方差函数的研究与应用

地表二向性反射分布函数(BRDF)是表征地物反射随太阳和观测方向变化的物理量。在统计意义上,BRDF表示均值统计量,BRVF(Bidirectional Reflectance Variance Function)表示方差统计量,它们对研究地表各向异性反射特征有着重要意义。本文首先采用误差传播理论,推导出基于MODIS BRDF模型的BRVF表达形式。研究结果表明,BRVF的空间分布模式主要由几何光学核Kgeo和体散射核Kvol的一次项和二次项权重和决定。然后利用EOS地面验证核心站点(EOS Land Validation Core Sites)的MODIS BRDF产品,对BRVF空间分布模式随地表类型、光谱波段和观测角度范围进行验证。验证结果表明,基于MODIS BRDF产品的验证与理论推导有较好的一致性。BRVF空间分布模式和地表类型有关,通常在热点处有一个峰值。在大观测天顶角(60°)下,BRVF随着角度的增大而增大。BRVF在近红外波段整体上大于红波段,表明其波段依赖性。最后,将上述理论成果初步应用于69组地表测量数据的模拟中。模拟结果表明,当大角度缺少观测数据时,模型外延所引起的方向反射方差显著增大,对地表反照率的反演精度和不确定性有较大影响。其中,红波段的白天空反照率的相对误差最大可达38.26%。本研究对利用小角度观测数据进行地表反照率反演的不确定性分析有指导意义;对大角度观测数据缺失情况下,先验知识在地表反照率的反演应用可提供有意义的理论支撑。     

机载WIDAS地表观测的BRDF原型反演算法验证

地表反照率表征地面对太阳辐射的反射能力,在地表能量平衡中起着重要的作用。2008年黑河的WATER实验中,针对机载红外广角双模式成像仪(WIDAS)观测角度小的特点,发展了一种基于MODIS二向性反射分布函数(BRDF)原型估算WIDAS反照率的算法。然而,在2012年黑河HiWATER实验中,WIDAS的观测角度范围由早期30°观测天顶角升级为最大观测天顶角52°,这一改造对数据预处理(辐射定标、大气校正和多角度配准)产生影响,使WIDAS数据出现明显的噪声,为WIDAS地表反照率的反演带来新挑战。针对新问题,采用新的地表观测数据,进一步验证了BRDF原型算法反演机载WIDAS反照率的能力。为了获取准实时BRDF先验知识,基于核驱动模型和各向异性平整指数(AFX)首先准实时提取了实验区5种MODIS的BRDF原型;然后,将其作为先验知识应用到黑河WIDAS机载数据的反照率反演中;最后,利用实验区的反照率实测数据进行算法验证。本研究比较了3个算法处理该问题的能力:(1)基于BRDF原型作为先验知识的算法;(2)完全基于核驱动的全反演算法;(3)基于朗伯假设的方法,通过算法对比分析,可更好地分析验证算法(1)的性能。结果表明:(1)BRDF原型算法反演精度最高且稳定,尤其是当观测数据位于垂直主平面时,全反演算法由于缺少大角度观测数据的约束,更容易出现异常值;BRDF原型算法的最大绝对误差为0.034,相对于全反演算法和朗伯假设方法精度分别提高了约18%和71%;(2)针对2012年黑河机载WIDAS观测噪声较大的新问题,BRDF原型反演算法表现出了良好的抗噪声性能。通过对WIDAS反照率反演和验证进一步表明,在多角度观测数据信息量不足以进行全反演,且观测受噪声影响较大时,BRDF原型作为先验知识,提供了一种有效解决该问题的手段,对提高反照率反演的精度、增强反演的稳定性有重要作用。     

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

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

一种改进的MODIS影像BRDF辐射校正方法

针时二向性造成的辐射失真，在核驱动模型的基础上提出了一种改进的二向性反射分布模型，并通过时MODIS影像的试验，得到了减弱辐射差异的较好效果，并对试验结果进行了分析和评价。     

新几何光学核的验证以及用核驱动模型反演地表反照率(之二)

基于核驱动模型的AMBRALS程序将用于提供MODIS的反照率和二向反射产品;该文提出一个新的几何光学核-LiTransit核,用以取代目前AMBRALS运行程序中的LiSparseR核。当植被及阴蔚总面积在视场中比例接近于;时,新核由LiSparseR核向LiDens。核自然诫。因而新核兼有公Sparse·核面LiDense核的优点,比LiSparseR核更符合几何光学模型的基本原理。对于29组地面观测数据的反演试验表明;LiTrarisit核有与LiSparse核相当的数据拟合能力,,外推到大的天顶角时有.ense核的稳定性;与LiSparseR核比较,RossThick-LiTransit的核组合更能反映直入扇出反照率随太阳天顶角变化的趋势。用NOAA AVHRR多时相数据反演的结果,也表明RossThick-LiTransit组合在反演甲更为稳定。     

改善MODIS BRDF产品热点效应的方法研究

核驱动的Ross Thick-LiSparse Reciprocal(RTLSR)双向反射函数(BRDF)模型已广泛地应用于MODIS等星载传感器的业务化产品处理中。但是,对于多年MODIS二向反射产品历史数据,如何基于RTMLSR模型发展一种简单有效快速的方法,进行热点效应的校正是一个迫切需要解决的问题。本文提出了一种简单有效的方法,不需要对观测数据重新反演,直接在现有MODIS二向反射产品的基础上进行热点校正,方便用户对历史MODIS二向反射产品的使用。该方法应用POLDER-3/BRDF数据库和部分经MODIS业务化算法筛选的反射率数据进行验证,并与RTLSR模型和RTMLSR模型的结果进行比较,结果表明:(1)该方法比现有的MODIS业务化RTLSR算法,对热点反射率有明显改善,拟合相对误差平均降低了10.12%;(2)该方法相对于RTMLSR模型在热点反射率拟合效果上差别不大,相对误差相差2.10%;(3)该方法对热点和冷点归一化的植被指数(NDHD)的估算效果有一定程度的改善,相对于RTLSR模型降低了约4.99%,与RTMLSR模型的相对误差相差1.32%,该方法对直接应用现有MODIS BRDF产品,基于热点方向反射率反演植被结构参数(如植被聚集指数)的精度提高有现实应用价值。     

正则化滤波在反演地表光谱反照率中的应用

利用有限的卫星观测数据,通过反演地物二向性反射分布函数(BRDF)模型提取复杂地表的真实信息,是一个不适定的地球科学反演问题。为了克服离散不适定性带来的困难,在分析引起模型解不稳定的原因和数据误差传递机制的基础上,使用构建滤波函数的方法实现BRDF模型的正则化约束反演。实验结果表明,正则化滤波算法与MODIS AMBRALS算法精度相当,完全适用于具有充足观测和稀疏观测情况下的地物参数反演。     

Evaluating the potential of multi-view data extraction from small Unmanned Aerial Systems (UASs) for object-based classification for Wetland land covers

Unmanned Aerial Systems (UASs) have the potential to provide multi-view data, but the approaches used to extract the multi-view data from UAS and investigation of their use in image classification are currently unavailable in publications to our best knowledge. This study presents a method that combines collinearity equations and a two-phase optimization procedure to automatically project a point from real world coordinate system of an orthoimage to UAS image coordinate system (row and column numbers) to be used in multi-view data extraction. The results show average errors for the computed UAS column and row numbers were 1.6 and 1.8 pixels respectively evaluated with leave-one-out method. Based on this algorithm, it’s also for the first time that object-based multi-view data were extracted and presented, and the potential of using the multi-view data to aid Geographic Object-Based Image Analysis(GEOBIA) through bidirectional reflectance distribution function (BRDF) modelling was evaluated with two representatives of BRDFs, the Rahman-Pinty-Verstraete(RPV) and Ross-Thick-LiSparse (RTLS). Our results indicate the RPV model tends to overestimate the bidirectional reflectance for land cover types with high reflectance, while perform well for those with relatively low reflectance in our study area. To test the impact of using multi-view data on image classification, we extracted parameters from BRDF models and used these parameters as object features for object-based classification. The 10-fold cross validation results show that the 3-parameter RTLS significantly improved overall accuracy compared to the classifications relying only on the orthoimage features, while other BRDF models did not show significant improvements, raising the needs to develop new methods to better utilize the multi-view information in GEOIBA in the future.     

Estimation of forest canopy structural parameters using kernel-driven bi-directional reflectance model based multi-angular vegetation indices

Near-surface bi-directional reflectance and high-spatial resolution true-color imagery of several forested canopies were acquired using an unmanned helicopter. The observed reflectance from multiple view-zenith angles were simulated with a kernel-driven bidirectional reflectance model, and the BRDF parameters were retrieved. Based on the retrieved BRDF parameters, kernel-derived multi-angular vegetation indices (KMVIs) were computed. The potential of KMVI for prediction of canopy structural parameters such as canopy fraction and canopy volume was assessed. The performance of each KMVI was tested by comparison to field measured canopy fraction and canopy volume. For the prediction of canopy fraction, the KMVI that included the nadir-based NDVI performed better than other KMVI emphasizing the importance of nadir observation for remote estimation of the canopy fraction. The Nadir BRDF-adjusted NDVI was found to be superior for the prediction of canopy fraction, which could explain 77% variation of the canopy fraction. However, none of the existing KMVI predicted the canopy volume better than Nadir BRDF-adjusted NDVI and Nadir-view NDVI. The Canopy structural index (CSI) was proposed with the combination of normalized difference between dark-spot near infrared reflectance and hot-spot red reflectance. The CSI could establish an improved relationship with the canopy volume over Nadir BRDF-adjusted NDVI and Nadir-view NDVI, explaining 72% variation in canopy volume. In addition, MODIS based KMVI were evaluated for the prediction of canopy fraction and canopy volume. MODIS based KMVI also showed similar results to the helicopter based KMVI. The promising results shown by the CSI suggest that it could be an appropriate candidate for remote estimation of three-dimensional canopy structure.     

评估MODIS的BRDF角度指数产品

应用地表观测的二向性反射数据集和多种MODIS数据产品,通过统计分析,对MODIS的二向性反射角度指数产品进行综合评估,结果表明:(1)MODIS角度指数包含了地表三维结构信息,有望用来反演地表的物理结构参数;(2)MODIS角度指数是内在的三维关系,各向异性因子(Anisotropic Factor:ANIF)和各向异性指数(Anisotropic Index:ANIX)高相关,建议去掉ANIF以精炼MODIS角度指数产品;(3)各向异性平整指数(Anisotropic FlatIndex:AFX)较好地指示了地表基本散射类型的变化,且具有较小的类内方差,对改善特定地表分类精度可能会更有用.     

Monitoring land surface albedo and vegetation dynamics using high spatial and temporal resolution synthetic time series from Landsat and the MODIS BRDF/NBAR/albedo product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warming/cooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500 m Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF/NBAR/albedo products and 30 m Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDF/Albedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30 m Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30 m albedos for the intervening daily time steps in this study. These enhanced daily 30 m spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of ±0.006. These synthetic time series provide much greater spatial detail than the 500 m gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 km by 14 km) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30 m resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces.     

Inter- and intra-annual variations of clumping index derived from the MODIS BRDF product

Clumping index quantifies the level of foliage aggregation, relative to a random distribution, and is a key structural parameter of plant canopies and is widely used in ecological and meteorological models. In this study, the inter- and intra-annual variations in clumping index values, derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) BRDF product, are investigated at six forest sites, including conifer forests, a mixed deciduous forest and an oak-savanna system. We find that the clumping index displays large seasonal variation, particularly for the deciduous sites, with the magnitude in clumping index values at each site comparable on an intra-annual basis, and the seasonality of clumping index well captured after noise removal. For broadleaved and mixed forest sites, minimum clumping index values are usually found during the season when leaf area index is at its maximum. The magnitude of MODIS clumping index is validated by ground data collected from 17 sites. Validation shows that the MODIS clumping index can explain 75% of variance in measured values (bias = 0.03 and rmse = 0.08), although with a narrower amplitude in variation. This study suggests that the MODIS BRDF product has the potential to produce good seasonal trajectories of clumping index values, but with an improved estimation of background reflectance.     

基于TM/ETM和DEM数据的BRDF特征提取

巴丹吉林沙漠地区地物类型单一,地形起伏,形成了天然的二向反射数据集;因此,本研究利用巴丹吉林沙漠地区的ASTER GDEM产品提供的地面高程数据,计算出每个坡元所对应的太阳-观测几何信息(包括太阳天顶角与方位角和观测天顶角与方位角),假设沙丘上每个坡元的表面结构不随其坡度和坡向变化,加上Landsat-TM/ETM+对地观测的信息,就形成了对同一地物的多角度观测数据集,从而可以提取该地区的BRDF特征。为了检验该方法,利用该方法获取的BRDF特征信息模拟了25景Landsat-TM/ETM+数据,并与实际的Landsat-TM/ETM+图像进行对比分析。结果表明, Landsat-TM/ETM+前4个波段的模拟图像与真实图像地表平均反射率相比,平均误差分别为2.80%、1.92%、2.68%和2.32%,高于一般辐射定标中5%—7%的误差要求,因此本研究方法可为高分辨率数据的交叉辐射定标等应用提供参考。     

针叶林冠层热点及多角度无人机遥感观测及特征分析

多角度遥感观测是研究植被冠层BRDF (Bidirectional Reflectance Distribution Function)特性的重要手段,但目前对森林冠层进行连续间隔采样的多角度遥感观测及数据较少,热点方向的观测尤为缺乏。本研究基于无人机多角度高光谱成像系统,在主平面上对针叶林冠层以等角度连续间隔采样进行多角度观测,获取了主平面上多角度(包括热点和暗点)高光谱影像,并将观测结果与四尺度几何光学模型模拟结果进行对比分析。多角度观测获取的植被冠层反射率呈现出典型的植被方向反射特征,后向大部分角度观测的冠层反射率高于前向,在热点处出现峰值,在暗点附近方向出现最低值,观测天顶角VZA (View Zenith Angle)较大时表现出明显的"碗边效应"。结果表明:(1)观测的针叶林冠层反射率及BRDF特性与四尺度模型模拟基本一致,但红光波段模拟的热点反射率稍低于观测,前向观测VZA较大时模拟与观测结果差异稍大;(2)冠层结构及叶片光学特性的差异会导致观测到的BRDF特征不同;(3)观测的针叶林冠层BRDF呈现明显的光谱效应,不同波段呈现的各向异性特性不同,红光波段各向异性最强,近红外波段最弱;(4) BRDF的光谱效应差异导致观测到的植被指数也表现出各向异性,NDVI (Normalized Difference Vegetation Index)、PRI (Photochemical Reflectance Index)和MTCI(MERIS Terrestrial Chlorophyll Index)在热点方向最低,EVI (Enhanced Vegetation Index)在热点方向最高。本研究中无人机多角度成像观测提供的角度和高光谱信息,尤其是热点和暗点信息,在地物识别及分类、植被冠层结构反演及碳循环关键参数获取等研究方面具有较好的应用前景,在其它地物反射或热辐射等方向性特性研究中也具有较大的潜力。     

A cloud detection method based on a time series of MODIS surface reflectance images

The Moderate Resolution Imaging Spectroradiometer (MODIS)-Terra surface reflectance product (MOD09A1), with bands 1 to 7, is a gridded, eight-day composite product derived from the MODIS-Terra top of atmosphere reflectance swaths. It performs cloud detection and corrects for the effects of atmospheric gases and aerosols. The cloud mask (CM) algorithms for MODIS are based on empirical thresholds on spectral reflectance and brightness temperature. Since the spatial resolution of the thermal band is 1000 m, while that of MOD09A1 is 500 m, many undetected and false clouds are observed in MOD09A1. These errors always result in temporal and spatial inconsistencies in higher-level products. In this paper, a cloud detection algorithm (TSCD) based on a MOD09A1 time series is introduced. Time series cloud detection (TSCD) algorithm is based on the relative stability of ground reflectance and the sudden variations in reflectance that result from cloud cover. The algorithm first searches the clear-sky reference data, and then discriminates clouded and unclouded pixels by detecting a sudden change of reflectance in the blue wavelength and spectral correlation coefficient at the pixel level. Compared with cloud cover assessments obtained from MODIS' original CM, TSCD provides similar or better discrimination in most situations when the land surface changes slowly.