著录项信息
专利名称 | 基于曲面论基本定律对地形曲面建立数字模型的方法 |
申请号 | CN200710064595.7 | 申请日期 | 2007-03-21 |
法律状态 | 撤回 | 申报国家 | 中国 |
公开/公告日 | 2008-09-24 | 公开/公告号 | CN101271596 |
优先权 | 暂无 | 优先权号 | 暂无 |
主分类号 | G06T17/50 | IPC分类号 | G;0;6;T;1;7;/;5;0查看分类表>
|
申请人 | 中国科学院地理科学与资源研究所 | 申请人地址 | 北京市安外大屯路3号917大楼
变更
专利地址、主体等相关变化,请及时变更,防止失效 |
权利人 | 中国科学院地理科学与资源研究所 | 当前权利人 | 中国科学院地理科学与资源研究所 |
发明人 | 岳天祥;宋敦江;杜正平 |
代理机构 | 中科专利商标代理有限责任公司 | 代理人 | 周国城 |
摘要
本发明属于地球观测与导航技术,涉及地理信息系统处理,公开一种基于曲面论基本定律对地形曲面建立数字模型的方法,将曲面进行空间离散化,然后根据曲面论中的高斯方程建立曲面离散化格网方程,求解超定线性系统,即可获得曲面的数字模型。本发明提高了数字地面模型的精度,可应用于卫星影像校正、双星定位系统以及国防信息化建设。
1. 一种基于曲面论基本定律对地形曲面建立数字模型的方法,步骤 如下:
步骤S1:采集地形表面曲面数据;
步骤S2:对地形曲面的采样数据进行空间离散化处理,得到离散化 的曲面数据Q(0);
步骤S3:根据曲面论,计算地形曲面数据Q(K),Q(K)基本要素为: 地形曲面上每个点的第一类基本量的系数E、F、G和第二类基本量的系 数L、M、N;
步骤S4:判断迭代次数K,第一次迭代时,K等于0;当K等于0 时,执行步骤S6;当K不等于0时,执行步骤S5;
步骤S5:判断迭代过程中相临两次曲面数据的差异Q(K)-Q(K-1), 当差异大于某一个阈值时,执行步骤S6;当Q(K)-Q(K-1)等于小于或等 于某一个阈值时,执行步骤S7;
步骤S6:根据曲面论中的高斯方程,结合地形曲面上点的数据参数, 建立地形曲面的差分微分方程为:S×f=T,求解差分微分方程,得到地 形曲面的模拟值为:f=S-1×T,完成一次地形曲面的迭代过程,并且K=K +1;其中S是曲面方程的系数矩阵,f是待求的曲面未知数,T为计算 出的曲面常数项;
步骤S7:输出地形曲面数据Q(K),得到高精度的数字地形模型。
2. 根据权利要求1所述的方法,其特征在于还包括,采用离散化曲 面的方法建立高斯方程,并通过迭代计算的方法进行求解。
3. 根据权利要求1所述的方法,其特征在于,所述曲面由第一类基 本量和第二类基本量决定。
4. 根据权利要求1所述的方法,其特征在于,所述地形曲面为蒙格 面片形式:
z=f(x,y),其中:x,y分别表示x坐标和y坐标。
5. 根据权利要求1所述的方法,其特征在于,所述地形曲面的第一 类基本量的系数为:
F=fxfy ,第一类基本量的系数E表示当 变化x时,变化将离曲面的多远,G表示当变化y时,变化将离曲面多 远,F表示当变化x,y同时变化时,他们之间的夹角是多大。
6. 根据权利要求1、4所述的方法,其特征在于,所述地形曲面的第 二类基本量的系数为:
第二类基本量的系数L、M、N与曲面的法向量曲率紧密相关。
7. 根据权利要求4所述的方法,其特征在于还包括,选择地形数学 曲面f(x,y)=2sin(πx)sin(πy)+1作为标准曲面,对标准曲面进行数值 模拟分析,当标准曲面顶部没有采样点时,经典模型的模拟结果有峰值 削平和边界振荡,曲面整体扭曲;将曲面整体扭曲经过HASM迭代后, 获得被削平的峰值,消除边界振荡,使模拟曲面与标准曲面相同。
8. 根据权利要求1所述的方法,其特征在于还包括,对地形数学曲 面f进行数值模拟时,先根据采样点的值求解出基本量 E、F、G、L、M、N,再利用基本量E、F、G、L、M、N求出曲面的第二类 克里斯托弗尔变量Γi j,i=11,12,22;j=1,2,在给定的边界条件下,模拟得 到曲面的数值解。
9. 权利要求1所述的方法用于空间图像数据处理。
引用专利(该专利引用了哪些专利)
序号 | 公开(公告)号 | 公开(公告)日 | 申请日 | 专利名称 | 申请人 | 该专利没有引用任何外部专利数据! |
被引用专利(该专利被哪些专利引用)
序号 | 公开(公告)号 | 公开(公告)日 | 申请日 | 专利名称 | 申请人 | 1 | | 2011-12-08 | 2011-12-08 | | |
2 | | 2011-01-19 | 2011-01-19 | | |
3 | | 2011-12-08 | 2011-12-08 | | |