MODIS卫星数据地表反照率反演的简化模式

以内蒙西部地区的MODIS遥感图像数据和地表野外同步观测的光谱数据为例,在野外数据量较少且有定标数据的条件下反演地表反照率。使用6S大气1辐射传输模型进行大气校正,并通过MODTRAN4.0模型获取各波段地表入射光通量和窄波段的地表反照率;在窄波段反照率与宽波段反照率之间存在线性关系的前提下,以各波段的入射光通量占总入射通量的比例作为反演参数,实现窄波段到宽波段的反演。反演结果证明此方法简便可行。     

HJ-1A/B星地表反照率生成算法与检验

提出了适合环境与灾害监测预报小卫星-A、B星(简称HJ-1A/B星)CCD相机的大气订正算法,并基于不同地表特性和大气条件下的辐射传输模拟数据,建立HJ-1A/B星的窄波段向宽波段反照率转换的模型.利用多级灰阶靶标实测数据、敦煌检验场实测数据验证了大气订正算法以及转换模型的可靠性和精度,并将HJ-1A/B星影像数据计算的反照率产品与同时相的MODIS反照率产品进行对比分析.结果表明:文章提出的HJ-1A/B星CCD相机大气订正算法可有效校正大气影响;窄波段向宽波段反照率转换模型反演的反照率精度可靠;基于研究成果生成的HJ-1A/B星地表反照率与MODIS反照率产品一致性较好,满足后续遥感数据定量化模型研究的精度需要.     

MODIS、MISR与POLDER 3种全球地表反照率卫星反演产品的比较与分析

对MODIS、MISR和POLDER 3种由多角度卫星观测反演得到的全球地表反照率数据(无冰雪覆盖区域)短波波段(SW,0.3~5 μm)与可见光波段(VIS,0.3~0.7 μm)的黑空地表反照率(DHR)进行了相互比较。3种产品6年平均的全球均值存在显著差异,其值从大到小依次为POLDER\,MISR和MODIS。3种产品的纬向平均在35°N以北区域表现出较大的差异。3种产品彼此之间相关性比较高,其中MODIS与MISR产品的相关性最强,MISR与POLDER产品的相关性最低,短波波段的相关系数(r) 分别为0.939与0.911。3种产品在可见光波段的相关性大于短波波段。在不同地表类型上,3种产品表现出了大致相似的差异,表明其对地表类型并不敏感。对气溶胶的分析表明:MODIS与MISR的550 nm气溶胶光学厚度(AOD)较为相似,其差异不足以解释DHR的差别;但是POLDER的865 nm AOD明显小于MISR,因此可以认为是由于POLDER的AOD估算偏低,导致了POLDER的DHR值大于MODIS与MISR。     

利用MODIS数据计算中国地表短波净辐射通量的研究

地表短波净辐射是辐射能量的重要收入部分,对研究地表辐射平衡、地气能量交换以及各种天气气候的形成都具有决定性的意义.目前国际上利用卫星数据计算地表短波净辐射大多采用空间分辨率较低的宽通道反照率数据(如ERBE卫星数据,空间分辨率为35km),这难以满足局部尺度的能量平衡和蒸散等研究的需要.本文利用空间分辨率达1km的窄通道多光谱MODIS卫星数据,通过相关处理计算,无需气象数据和地面测量数据的参与,就获得了中国的地表短波净辐射通量分布图.经与禹城地区地面实测数据对比可知,晴空时均方根误差小于20 W/m2,有云时均方根误差小于35 W/m2,表明本文使用的地表短波净辐射反演方法适合中国范围地表短波净辐射的计算.     

基于TM影像天山北坡地表反照率反演方法的研究

地表反照率是气候模型和地表能量平衡方程中的重要参数。基于6S模型估算地表反照率,忽略了地形起伏的影响,不同波段组合的地表反照率也有待进一步研究。本文选取干旱区典型流域——天山北坡三工河流域为研究区域,以TM影像为数据源逐步进行地形校正、大气校正等,从而提取窄波段地表反照率。在此基础上,根据亮度、绿度、湿度3个特征变量的物理意义,以各波段能量权重为转换参数对窄波段地表反照率进行组合,实现研究区宽波段地表反照率的反演,得出基于不同波段的物理意义的地表反照率。     

应用MODIS数据监测巢湖蓝藻水华的研究

以巢湖为研究区域,以MODIS 卫星影像为数据源,结合准同步的地面水质监测数据,将MODIS 250 m分辨率的波段反射率与叶绿素a浓度实测值进行相关分析。在此基础上通过回归拟合,构建基于中分辨率成像光谱仪(MODIS) 的叶绿素遥感提取模型。应用模型成功提取出蓝藻爆发水域chl-a的分布。从MODIS遥感图像上可以清晰地反映出巢湖这次蓝藻爆发的强度、地点和分布范围 。研究结果表明:用MODIS影像监测巢湖蓝藻水华是可行的,其中250m分辨率波段1 、2的比值组合r2/r1与叶绿素a浓度实测值高度相关(R=0.909 3),适于反演叶绿素a浓度。     

融合遥感先验信息的叶面积指数反演

借助植被辐射传输模型,利用遥感观测数据估算LAI是一种较为可靠和稳健的反演方法。然而,地表的复杂性、遥感观测的有限性以及自相关性导致遥感数据包含的信息量不足,不能完全支持LAI等地表参数的估算,易造成“病态”反演。在遥感反演过程中引入先验知识能够有效地解决该问题。研究基于遥感数据提取LAI先验信息,并将其用于代价函数的构建,利用PROSAIL辐射传输模型和遗传算法,分别在500 m和250 m尺度反演LAI。将高空间分辨率LAI分别升尺度到500 m和250 m,验证对应尺度LAI结果,评价引入先验信息对于提高LAI反演精度的作用。研究表明,引入先验信息有助于提高不同分辨率下LAI反演精度,且先验信息的质量一定程度上也影响着LAI反演结果。与未加入先验信息的LAI反演结果相比,以MODIS LAI产品作为先验信息反演的500 m尺度LAI结果精度R²由0.55提高至0.65,RMSE由1.29下降至0.38。在250 m尺度,以500 m LAI反演结果作为先验信息反演的叶面积指数,其精度优于以MODIS LAI产品为先验知识的估算结果,验证精度R²增加了0.08,RMSE减少了0.18。研究使用的先验信息主要来自遥感数据本身,没有地面实测数据的参与,在此基础上发展的多分辨率LAI反演方法具有估算大区域尺度LAI的应用潜力。     

融合遥感先验信息的叶面积指数反演

借助植被辐射传输模型,利用遥感观测数据估算LAI是一种较为可靠和稳健的反演方法。然而,地表的复杂性、遥感观测的有限性以及自相关性导致遥感数据包含的信息量不足,不能完全支持LAI等地表参数的估算,易造成"病态"反演。在遥感反演过程中引入先验知识能够有效地解决该问题。研究基于遥感数据提取LAI先验信息,并将其用于代价函数的构建,利用PROSAIL辐射传输模型和遗传算法,分别在500 m和250 m尺度反演LAI。将高空间分辨率LAI分别升尺度到500 m和250 m,验证对应尺度LAI结果,评价引入先验信息对于提高LAI反演精度的作用。研究表明,引入先验信息有助于提高不同分辨率下LAI反演精度,且先验信息的质量一定程度上也影响着LAI反演结果。与未加入先验信息的LAI反演结果相比,以MODIS LAI产品作为先验信息反演的500 m尺度LAI结果精度R2由0.55提高至0.65,RMSE由1.29下降至0.38。在250 m尺度,以500 m LAI反演结果作为先验信息反演的叶面积指数,其精度优于以MODIS LAI产品为先验知识的估算结果,验证精度R2增加了0.08,RMSE减少了0.18。研究使用的先验信息主要来自遥感数据本身,没有地面实测数据的参与,在此基础上发展的多分辨率LAI反演方法具有估算大区域尺度LAI的应用潜力。     

融合遥感先验信息的叶面积指数反演

借助植被辐射传输模型，利用遥感观测数据估算LAI是一种较为可靠和稳健的反演方法。然而，地表的复杂性、遥感观测的有限性以及自相关性导致遥感数据包含的信息量不足，不能完全支持LAI等地表参数的估算，易造成“病态”反演。在遥感反演过程中引入先验知识能够有效地解决该问题。研究基于遥感数据提取LAI先验信息，并将其用于代价函数的构建，利用PROSAIL辐射传输模型和遗传算法，分别在500 m和250 m尺度反演LAI。将高空间分辨率LAI分别升尺度到500 m和250 m，验证对应尺度LAI结果，评价引入先验信息对于提高LAI反演精度的作用。研究表明，引入先验信息有助于提高不同分辨率下LAI反演精度，且先验信息的质量一定程度上也影响着LAI反演结果。与未加入先验信息的LAI反演结果相比，以MODIS LAI产品作为先验信息反演的500 m尺度LAI结果精度R²由0.55提高至0.65，RMSE由1.29下降至0.38。在250 m尺度，以500 m LAI反演结果作为先验信息反演的叶面积指数，其精度优于以MODIS LAI产品为先验知识的估算结果，验证精度R²增加了0.08，RMSE减少了0.18。研究使用的先验信息主要来自遥感数据本身，没有地面实测数据的参与，在此基础上发展的多分辨率LAI反演方法具有估算大区域尺度LAI的应用潜力。     

基于Sentinel-1与FY-3C数据反演植被覆盖地表土壤水分

基于新一代的Sentinel-1SAR数据与FY-3C的MWRI数据,研究植被覆盖地表土壤湿度反演方法。为消除植被对土壤湿度反演影响,首先利用FY-3C/MWRI的微波极化差异指数MPDI,建立植被含水量反演模型;然后,结合植被含水量反演模型和水—云模型,发展一种主被动微波联合反演植被覆盖地表土壤含水量模型;最后,在江淮地区开展反演试验,利用观测的土壤湿度数据进行反演结果的精度验证。结果表明:(1)对于植被覆盖地表土壤湿度反演,由FY3C/MWRI提取的MPDI对于去除植被影响效果较好;(2)相比于VH极化哨兵1号卫星数据,VV极化数据更适用于土壤含水量的反演,能够得到更高的土壤湿度反演精度;(3)哨兵1号卫星数据能够获得较高精度的土壤含水量反演结果,试验反演的土壤湿度值与实测值相关系数为0.561 2,均方根误差为0.044cm~3/cm~3。     

Validation of Narrow-band Surface Albedo Retrieved from FY-3C MERSI Satellite Data

Surface albedo is one of the driving factors in surface radiant energy balance and surface-atmosphere interaction.It is widely used in surface energy balance, medium and long-term weather forecasting and global change research.This study aims to validate the surface albedo retrieved from FY-3C MERSI. This paper selected four regions in Africa and North America as study areas to validate the retrieved albedo from the reflectance data and angle data of FY-3C MERSI at 250 m resolution in 2014. The semi-empirical kernel-driven BRDF(bidirectional reflectance distribution function) model RossThick-LiSparseR and least squares fitting method were used to calculate the parameter of BRDF. Then four narrow-band black-sky albedos and four narrow-band white-sky albedos can be obtained by angle integration. Finally, the cross-validation of FY-3C surface narrow-band albedo products with MODIS albedo products (MCD43A3) was carried out. The results show that theRoot Mean Square Error(RMSE) between the FY-3C narrow-band albedo and the corresponding MODIS narrow-band albedo is in the range of 0.01 ~ 0.04, and the Mean Bias (MBIAS) is 0.09. FY-3C narrow-band albedo has good consistency with the corresponding MODIS narrow-band albedo in the visible and near-infrared bands. So, the methodologyof using the BRDF model to invert the surface albedo of FY-3C medium resolution imaging spectrometer data is feasible and reliable. The further improvement of the inversion accuracy of FY3C-MERSI surface albedo also depends on the improvement of basic data processing quality, including image geometric correction, calibration, and strict data quality control.     

Validation of MISR land surface broadband albedo

Land surface broadband albedo is a critical variable for many scientific applications. Due to the scarcity of spectral albedo measurements of the Earth's surface environments, it is useful to construct broadband albedo from spectral albedo data obtained by multi‐angle satellite observations. The Multi‐angle Imaging SpectroRadiometer (MISR) onboard NASA's Earth Observing System (EOS) Terra satellite provides land surface albedo products from multi‐angular observations; however, the products have not been comprehensively validated. We convert MISR spectral albedos to total shortwave albedos and validate them using ground measurements at different validation sites. For most surface types, a published narrowband to broadband conversion formula was used, but a new conversion formula for snow and ice covered sites is developed in this study where the spectral range of the instrument is different. Several comparisons are made: (1) between MISR directional‐hemispherical reflectance (DHR) or albedo and MODIS (Moderate Resolution Imaging Spectroradiometer) DHR; and (2) between MISR spectral DHR and bi‐hemispherical reflectance (BHR). The results show that: (1) both the value and the temporal trends of the MISR shortwave albedo and the ground measured shortwave albedo are in good agreement, with the exception of the snow and ice sites; (2) the MISR DHR conforms well to MODIS DHR; and (3) the values of MISR DHR and BHR are nearly identical.     

Validation of the MODIS reflectance product under UK conditions

Surface reflectance obtained from remote-sensing data is the main input to almost all remote-sensing applications. The availability and special characteristics of Moderate Resolution Imaging Spectroradiometer (MODIS) products have led to their use worldwide. Validation of the MODIS reflectance product is then crucial to provid information on uncertainty in the reflectance data, and in other MODIS products and in the applied surface–atmosphere models. Compact Airborne Spectrographic Imager (CASI) and Système Pour l'Observation de la Terre (SPOT) data, collected during the Network for Calibration and Validation in Earth Observation (NCAVEO) 2006 Field Campaign, were applied to validate daily MODIS reflectance data over a site in the southern UK. The difference in the view geometry of at-nadir CASI and SPOT data and off-nadir MODIS data was dealt with using a semi-empirical bidirectional reflectance distribution function (BRDF) model. The validation results showed that for our particular study site, the absolute errors in the MODIS reflectance product were too large to allow the albedo data to be used directly in climate models. The errors were mainly related to the uncertainties in the MODIS atmospheric variables, the BRDF model, and undetected clouds and cloud shadows. More generally, the study highlights the extreme difficulty of achieving pixel-level validation of coarse spatial resolution satellite sensor data in an environment in which the atmosphere is constantly changing, and in which the landscape is characterized by high space–time heterogeneity.     

基于ART模型的MODIS积雪反照率反演研究

积雪反照率是研究局地或全球的能量收支平衡和气候变化中的重要参数,遥感反演为积雪反照率的获取提供了便利的手段。积雪反照率大小主要取决于积雪的自身物理属性(雪粒径、形状和污染物等因子)以及天气状况,遥感反演反照率大多基于双向反射模型(BRDF),积雪BRDF模型常使用积雪辐射传输模型获得。采用考虑了雪粒径、粒子形状以及污染物影响的渐进辐射传输理论(ART)模型,建立了MODIS积雪反照率反演算法,得到了MODIS 8d合成积雪反照率产品。将此算法应用于具有均一积雪地表的格陵兰岛地区,并使用GC-Net实测数据进行了验证,反演的总均方根误差(RMSE)为0.018,相关系数(r)为0.83,结果表明考虑了积雪特性的ART模型能够较好地反演积雪反照率,而且反演需要的参数较少。     

Assessment of biases in MODIS surface reflectance due to Lambertian approximation

Using MODIS data and the AERONET-based Surface Reflectance Validation Network (ASRVN), this work studies errors of MODIS atmospheric correction caused by the Lambertian approximation. On one hand, this approximation greatly simplifies the radiative transfer model, reduces the size of the look-up tables, and makes operational algorithm faster. On the other hand, uncompensated atmospheric scattering caused by Lambertian model systematically biases the results. For example, for a typical bowl-shaped bidirectional reflectance distribution function (BRDF), the derived reflectance is underestimated at high solar or view zenith angles, where BRDF is high, and is overestimated at low zenith angles where BRDF is low. The magnitude of biases grows with the amount of scattering in the atmosphere, i.e., at shorter wavelengths and at higher aerosol concentration. The slope of regression of Lambertian surface reflectance vs. ASRVN bidirectional reflectance factor (BRF) is about 0.85 in the red and 0.6 in the green bands. This error propagates into the MODIS BRDF/albedo algorithm, slightly reducing the magnitude of overall reflectance and anisotropy of BRDF. This results in a small negative bias of spectral surface albedo. An assessment for the GSFC (Greenbelt, USA) validation site shows the albedo reduction by 0.004 in the near infrared, 0.005 in the red, and 0.008 in the green MODIS bands.     

Variability in surface BRDF at different spatial scales (30 m-500 m) over a mixed agricultural landscape as retrieved from airborne and satellite spectral measurements

Over the past decade, the role of multiangle remote sensing has been central to the development of algorithms for the retrieval of global land surface properties including models of the bidirectional reflectance distribution function (BRDF), albedo, land cover/dynamics, burned area extent, as well as other key surface biophysical quantities impacted by the anisotropic reflectance characteristics of vegetation. In this study, a new retrieval strategy for fine-to-moderate resolution multiangle observations was developed, based on the operational sequence used to retrieve the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5 reflectance and BRDF/albedo products. The algorithm makes use of a semiempirical kernel-driven bidirectional reflectance model to provide estimates of intrinsic albedo (i.e., directional-hemispherical reflectance and bihemispherical reflectance), model parameters describing the BRDF, and extensive quality assurance information. The new retrieval strategy was applied to NASA's Cloud Absorption Radiometer (CAR) data acquired during the 2007 Cloud and Land Surface Interaction Campaign (CLASIC) over the well-instrumented Atmospheric Radiation Measurement Program (ARM) Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site in Oklahoma, USA. For the case analyzed, we obtained ~ 1.6 million individual surface bidirectional reflectance factor (BRF) retrievals, from nadir to 75° off-nadir, and at spatial resolutions ranging from 3 m to 500 m. This unique dataset was used to examine the interaction of the spatial and angular characteristics of a mixed agricultural landscape; and provided the basis for detailed assessments of: (1) the use of a land cover type-specific a priori knowledge in kernel-driven BRDF model inversions; (2) the interaction between surface reflectance anisotropy and instrument spatial resolution; and (3) the uncertainties that arise when sub-pixel differences in the BRDF are aggregated to a moderate resolution satellite pixel. Results offer empirical evidence concerning the influence of scale and spatial heterogeneity in kernel-driven BRDF models; providing potential new insights into the behavior and characteristics of different surface radiative properties related to land/use cover change and vegetation structure.     

An algorithm for the retrieval of 30-m snow-free albedo from Landsat surface reflectance and MODIS BRDF

We present a new methodology to generate 30-m resolution land surface albedo using Landsat surface reflectance and anisotropy information from concurrent MODIS 500-m observations. Albedo information at fine spatial resolution is particularly useful for quantifying climate impacts associated with land use change and ecosystem disturbance. The derived white-sky and black-sky spectral albedos may be used to estimate actual spectral albedos by taking into account the proportion of direct and diffuse solar radiation arriving at the ground. A further spectral-to-broadband conversion based on extensive radiative transfer simulations is applied to produce the broadband albedos at visible, near infrared, and shortwave regimes. The accuracy of this approach has been evaluated using 270 Landsat scenes covering six field stations supported by the SURFace RADiation Budget Network (SURFRAD) and Atmospheric Radiation Measurement Southern Great Plains (ARM/SGP) network. Comparison with field measurements shows that Landsat 30-m snow-free shortwave albedos from all seasons generally achieve an absolute accuracy of ±0.02-0.05 for these validation sites during available clear days in 2003-2005, with a root mean square error less than 0.03 and a bias less than 0.02. This level of accuracy has been regarded as sufficient for driving global and regional climate models. The Landsat-based retrievals have also been compared to the operational 16-day MODIS albedo produced every 8-days from MODIS on Terra and Aqua (MCD43A). The Landsat albedo provides more detailed landscape texture, and achieves better agreement (correlation and dynamic range) with in-situ data at the validation stations, particularly when the stations include a heterogeneous mix of surface covers.     

FY3A/MERSI地表温度反演

MERSI是我国第二代极轨气象卫星上的重要传感器,可获取高空间分辨率和高时间分辨率的对地观测影像。为使Jimènez-Mu珘nozSobrino算法更适用于FY3A/MERSI传感器通道特性,更新了大气函数的估算系数,并引入观测角度因子,以获取更为精确像元间更为平滑的地表温度。用MODTRAN4模拟验证该算法精度,得引入角度因子后反演精度显著提升,所有角度下平均误差为-0.6±2.2K。用实测的敦煌戈壁地表温度和MODIS地表温度产品评价MERSI反演结果,显示MERSI地表温度的空间分布准确,结果精度也较高。与实测温度对比,平均误差为1.74K,均方根误差小于1.9K。研究区域与MODIS地表温度间差异平均为2.6307K。虽然会受云检测精度和观测亮温偏高的影响,由MERSI反演的高精度地表温度在相关科研和业务方面仍然具有极好的应用前景。     

Theoretical noise sensitivity of BRDF and albedo retrieval from the EOS-MODIS and MISR sensors with respect to angular sampling

The sensitivity of the semiempirical RossThick-LiSparse Ambrals BRDF model to random noise in observed multiangular reflectances was investigated through a study of the impact of angular sampling. The mathematical properties of (linear, additive) kernel-driven BRDF models allow the analytical derivation of so-called weights of determination or noise amplification factors which quantify the uncertainty in retrieved parameters such as nadir-view reflectance or albedo at various solar zenith angles, or in the BRDF model parameters themselves. The study was carried out using simulated angular sampling for the MODIS and MISR instruments to be flown on NASA's Earth Observing System AM-1 platform, as a function of latitude, day of year and sampling period. A similar study was carried out for comparison using the modified RPV BRDF model, a multiplicative model. Results show that the retrieved parameters, reflectance and albedo can be expected to have noise amplification factors that are less than unity, indicating that the retrievals are stable with respect to random noise under the angular sampling schemes occurring. The BRDF model parameters themselves were found to be more susceptible to noise than many of the derived products, especially for the modified RPV model. The effect of different angular sampling regimes on the uncertainty of derived information was further explored. This study provides an indication of the reliability to be expected from the operational BRDF/albedo products from the MODIS and MISR instruments. The findings may qualitatively also apply to AVHRR, SPOT VEGETATION and similar satellite angular sampling regimes.