带先验约束的地表参数提取的有效反演方法.pdfVIP

中国科学 D 辑 :地球科学 2009 年 第 39 卷 第 3 期 : 360 ~ 369 SCIENCE IN CHINA PRESS 带先验约束的地表参数提取的有效反演方法 ①* ④ ②③ ②③ ②③ 王彦飞 , Shiqian Ma , 杨华 , 王锦地 , 李小文 ① 中国科学院地质与地球物理研究所油气综合地球物理重点实验室, 北京 100029; ② 北京师范大学遥感与 GIS 中心, 北京 100875; ③ 北京师范大学与中国科学院遥感应用研究所联合遥感科学国家重点实验室, 北京 100875 ； ④ Department of Industrial Engineering and Operations Research, Columbia University, New York, NY 10027-6902, USA * E-mail: yfwang_ucf@ 收稿日期: 2008-06-12; 接受日期: 2008-10-22 国家自然科学基金项目(批准号:和国家重点基础研究发展计划(编号: 2007CB714400, 2005CB422104)资助 摘要 地球表面的各向异性特性可以用地表二向反射函数(BRDF)恰当地描述. BRDF 的核 关键词 心是利用线性核驱动模型, 数学上表述为各向同性核、体散射核和几何光学核的线性组合. 随 不适定问题 着多角度遥感领域的发展, BRDF 模型越来越被看作是可以反演重要的有关地表生物的或气候 地表参数提取 的参数, 比如说叶面积指数和地表反照率. 一个线性逼近的核驱动BRDF 模型通常可以写成下 最优化 正则化 述形式(Roujean 等, 1992):f +k (t , t , φ)f +k (t , t , φ)f =r(t , t , φ), 其中r 表示地表的二向反 iso vol i v vol geo i v geo i v 射; kvol 和kgeo 为通常所说的核函数, 即为已知的入射和观测几何特性的函数, 分别描述了体散 射和几何散射(包括折射和反射); t 是太阳方向天顶角, t 是观测方向天顶角; φ表示太阳-观测方 i v 向的相对方位角;f iso,f vol 和f geo 为未知的待反演参数, 可以用来拟合观测. 计算过程的稳定性是 由核矩阵的代数算子特征谱和观测噪音/误差来刻画的. 因此为了计算地表反照率, 成功反演 模型参数是至关重要的环节. 我们首先考虑了为计算 BRDF 模型反演的光滑解方法. 业已知 道, 这是一个不适定的反问题. 不适定性是由线性核驱动BRDF 模型的欠定性表征的, 比如说 观测严重不足或观测方向范围有限, 或者是观测数据高度线性相关以及噪音的污染等. 例如, 一次单角度观测可以导致一个欠定的系统(核算子的零空间含有非零向量)或者系统无解(系数 矩阵的秩不等于增广矩阵的秩). 因此, 光滑性或正则化技巧应当加以利用来压制不适定性. Li 等(2001)应用先验知识把原始模型转换为一个超定的模型并求得最小二乘解. Pokrovsky 等 (2002)应用QR 分解反演BRDF 模型. Wang 等(2007)考虑到了反演的正则化策略并提出了不适 定地表参数反演的一个完整的正则化理论. 在文中,