一种由采样定理生成较高分辨率dem的方法.docxVIP

一种由采样定理生成较高分辨率dem的方法 dem的规模效应对dem应用的适应性和可靠性以及dem是否真的适合实践，是当前研究的热点。为了验证分辨率，有必要连续修改dem，并将空间定位的一致性作为数据源。其中之一是使用不同的尺度地形图生成不同的分辨率dem。然而，在此过程中，存在空间定位基础匹配的精度问题，该误差对后续分析具有一定的不确定性。图1显示了不同分辨率的dem由专业插值软件在洛杉矶高原上的重叠软件中生成的，并且具有相应的堆叠特性。结果表明，这种不确定性存在于一侧的数据集，并且该方法无法生成分辨率变化的dem数据集。 另一途径是首先基于大比例尺地形图建立高分辨率DEM,然后通过对其表面综合概括,生成空间定位基础一致的、分辨率较低的DEM.常规方法有双线性、双三次等内插方法,但这些方案对于同一地区不同分辨率的DEM均需新构建,各种分辨率之间无法进行转换,同时缺乏简捷有效的逆过程.小波变换等方法具备时频多分辨率特性,能够在不同的时频分辨率下表达对象不同详细程度的特征信息,较好地弥补了上述缺陷,受到学者关注.然而,对于如何基于该方法进一步生成分辨率可连续变化或能自由指定的DEM,目前研究较少,且主要是通过直接对小波高频系数作阈值处理,保留不同百分比数量的较大系数来实现. DEM是一类二维信号,其空间域的采样率(或采样间隔,即水平分辨率)与其频域的最高频率之间存在着一定关系,而现有研究尚未对其深入分析.本文首先从采样定理出发分析了这种关系,在此基础上,探讨了如何通过定量降低DEM的有效最高频率来实现降低其采样率至指定值,即如何通过运用两者关系生成指定分辨率的综合DEM;然后借助于分数进制小波变换实现了这一方法;最后,通过目视对比分析、高程统计指标分析、坡度与剖面曲率频率分布分析、等高线与河流总长度统计分析等,对所得的综合DEM进行了实验论证,结果显示了本方法在一定范围内的适用性. 本文探讨的DEM均为格网型.对DEM的分析、综合处理在Matlab平台中参照相关算法编程实现;之后对所得综合DEM的实验论证在Matlab与ArcGIS平台中联合实现. 1 方法分析 1.1 信号采样理论指导下的采样 地面高程模型(DEM)是实际地表高程表面的逼近.具有不同逼近程度或详细程度信息的DEM所适合的栅格尺寸有所不同:栅格尺寸过大将导致某些细节信息无法表达,过小则造成数据冗余.此最适合的栅格尺寸即为DEM的水平分辨率.图2显示了同一地区两幅具有不同详细程度地表高程信息的DEM,两者各自所适合的栅格尺寸也不相同.DEM作为一种二维信号,其分辨率的选取问题即是其重采样时采样间隔的选取问题,因此在这一过程中需要考虑到采样定理的约束. 假定二维信号g(x,y),x,y∈(-∞,+∞)具有有限带宽,G(u,v)为其傅里叶变换,fmax_u和fmax_v分别为其在u、v方向上的带宽或最高频率(此处频率的含义为单位长度内高程值的周期性变化次数),则有:当|u|?[0,fmax_u]或v?[0,fmax_v]时,G(u,v)=0,即原信号在频域(uv域)的频谱信息集中在长宽为2fmax_u、2fmax_v的矩形范围内,如图3所示.数字技术的广泛应用需要连续信号的数字化,采样是数字化的首要过程.为确保采样时信息不丢失,采样时所采用的采样间隔或采样率必须满足采样定理,即 {1Δx≥2fmax_u或fsam_x≥2fmax_u1Δy≥2fmax_v或fsam_y≥2fmax_v{1Δx≥2fmax_u或fsam_x≥2fmax_u1Δy≥2fmax_v或fsam_y≥2fmax_v 当信号的频谱信息并没有集中分布于一个有限的范围时,可以采用一个低通抗混叠滤波器对其作预处理,使其最高频率降至某指定值.采样定理在采样过程中确保信息不丢失方面起着理论上的保证作用,其是否能得到满足往往关系到后续处理的成败与质量. DEM水平、垂直两个方向的采样间隔或采样率通常相等,即有fsam_x=fsam_y.根据离散傅里叶变换(DFT)理论,在一个主值周期内,有fmax_u=fmax_v(鉴于DFT的快捷有效性,本文中采用DFT进行傅里叶分析).由此以下用fsam表示空间域两个方向的采样率,fmax表示为频域两个方向的最高频率(简称最高频率).采样定理要求fsam与fmax的比值n必须大于等于2;但n也不可过大,以避免不必要的数据冗余. 在建立了采样率与最高频率之间的关系后,如能通过某种方法将具有较高采样率的DEM(以下称为原DEM)的最高频率降低至指定值,即可尝试获得所适合的采样率降至指定值的综合DEM. 1.2 原代地线域划分 对于如何由原DEM生成指定分辨率的综合DEM,基本思路为: (1)从频域能量分布的角度分析如何更有效地定义与计算DEM的频域最高频率; (2)先假定通过