论文部分内容阅读
本课题来源于国家自然科学基金项目“三维表面形状被动重构方法研究”(项目号:60141002)利“十五”国防预研课题“军用图像并行处理算法技术”(项目号:413160103),并且得到了与英国Surrey大学共同申请的英国皇家学会基金项目“三维表面重构(3D SurfaceReconstruction)”(批准号:Q775)的经费支持。在国防和民用等应用领域中,地形模型一直扮演着一个非常重要的角色。尽管大量遥感图像(航片、卫片等)间接提供了丰富的地形信息,但现有三维表面重建技术未能从中有效提取地形表面。另一方面,分形技术已被广泛应用于地形模拟,但却没有与三维重构技术有效结合。为此本文以一种重要的三维重构技术——从明暗恢复形状(Shape-From-Shading,SFS)技术为主,研究和探讨了基于分形的地形图像表面三维重建问题,不仅在方法上克服了传统分形描述工具难以直接与观测图像数据相结合的困难,而且在一定误差范围内可以从观测地形图像中有效提取相应地表形状。本文做出了以下贡献:通过回顾和分析不同约束条件和先验知识在SFS问题中的应用,本文得出以下结论:(1)无论给出确定性约束条件还是随机先验分布,相应的SFS问题均可表示为一个优化问题;(2)SFS问题的约束条件需同时满足先验知识要求和算法稳定性;(3)SFS算法中参数确定和表面形状获取交互迭代或同时迭代;(4)由于反射函数的高度非线性,在SFS算法中采用泰勒线性化方法容易使求解过程陷入局部极小状态。本文在分析分形布朗曲面各分形指标计算的共性基础上,提出了一种利用分形特征变量相关性获取近似分形表面模型的一般性方法。并分别以基于方差图、功率谱、逐层增量以及小波等分析方法为基础,实现了四种基于不同特征变量的具体分形表面模型。构造结果表明,这四种具体近似模型均能有效地反映分形表面的统计特征,较好地逼近了非线性分形表面模型。与传统分形工具相比,本文方法构造的分形表面模型更为直观,易于控制和推广,便于与观测数据相结合并可简化计算。求解高度非线性的SFS问题时,需要构造一系列线性最小二乘子问题。通常采用泰勒线性化并假定其规测噪声为Gauss分布,容易导致求解过程陷入局部极小。为此本文提出了一种改进方法,即应用Pentland线性化结构将SFS问题的非线性部分转化为系统参数,使其观测噪声具有非Gauss分布,并可通过本文分形模型对其进行估计。实验表明,该办法可较好,地避免问题的解陷入局部极小状态。本文利用正则化理论将基于分形的SFS问题表述为一个分形正则化问题,并提出了一种分形正则化SFS算法。其特点是可由迭代次数控制表面分形维数,且具有可控性和灵活性。实验表明,该算法可有效地从地形图像中恢复地表形状。与传统SFS方法相比,该方法的地表恢复结果保持了地形原有的分形特征,而又不易受到噪声影响。对于表面呈现非均匀分形特征的情况,本文利用局部分形分析,分形正则化算法和块Jacobi迭代框架,提出了一种具有内外层迭代结构的局部分形正则化SFS算法。其内层各块迭代次数控制了相应的局部分形维数,外层循环次数则控制了表面的整体一致性。实验表明,该算法可有效地从复杂图像中提取地形表面,可在一定程度上消除某些具体应用中出现的表面拓扑错误。针对复杂的多重分形地形图像的SFS问题,本文提出了一种基于多级优化方法(Multi-Level Optimization,MLO)的多尺度分形正则化SFS算法。该方法利用较大尺度取得的表面结果,对较小尺度中的局部分形迭代常量作修正,达到多尺度分形约束的效果。其各尺度段子块的迭代次数控制了多重分形特征。实验表明,该算法可有效提取复杂地表形状。与经典SFS算法相比,本文算法易于作并行计算和处理大型图像问题,并可灵活控制分形特征,具有可扩展性。