GIS在区域地质应用的两个数据处理技术问题.pdfVIP

工程地质地理信息系统研讨会资料汇编 2000·北京 GI S在区域地质应用的两个数据处理技术问题 严泰来 朱德海 (中国农业大学100094) 【摘要】本文就区域地质研究两个常见的数据处理问题，即数值分析与数据加密填图问题，提供 开发计算机G1S栅格格式图形处理功能的两个算法。此算法引进了遥感图像处理的数字滤波理论与 图形分形理论，对于GIS中栅格格式数据处理这一较为薄弱的部分做一技术补充，为拓宽GIS在 区域地质应用研究领域提供可资借鉴的思想。 【关键剐数字滤波、翻、．盐毯面插值 区域地质研究在地学中是一个基础性的研究领域，具有图件比例尺跨度大，数据 属性庞杂，空间分析复杂及要求高等诸多特点。地理信息系统是将空间信息数字化， 使用计算机存储量大、数据处理速度快的强大功能，进行空间分析的系统，因而地理 信息系统出现后不久，这一技术很快就应用于区域地质研究领域，并获得很大成功： 为了适应区域地质研究深入的需要，这里介绍在GIS开发应用中具有普遍应用意义的 两个数据处理的实际问题。 地质普查数据如重力场分析、磁场分析、电导率分布等等都是以方里网采样得到 的数据，这些数据在GIS中称之为栅格数据，在小比例尺图件中方里网格就是栅格， 每个方里网格中的重力加速度G、磁场强度、磁偏角、电导率就是其栅格的属性。栅 格格式数据是GIS空间信息的重要表达方式，尽管在目前这种表达方式还不是主流表 达方式，但是不少人认为按照GIS发展的趋势，相信栅格格式将会取代矢量格式成为 地学空间信息的主流表达格式，其原因是栅格格式数据处理简单，空间分析便捷，特 别适宜于数据场的动态分析。栅格格式的缺点诸如数据精度差、存储空间占用量巨大 等缺点，随着计算机内存与磁盘存储空间的扩展将逐渐消失。作为GIS在区域地质上 的应用，我们必须对栅格格式的数据处理给予充分的重视，充分发挥GIS在这种数据 格式下的功能优势，以便从海量数据中挖掘、提取我们感兴趣的信息。 一、数字滤波技术 所谓的数字滤波是基于这样一个思想，即：将栅格格式下的属性数据空间分布A (I)看作是一个随时间变化的非周期的复杂波函数， G(，)可(t)……………………………………………………～(1) 上式中的I为网格编号。 工程地质地理信息系统研讨会资料汇编 2000·北京 按照傅里叶(Fourier)分析的思想，任何一个函数，包括周期函数与非周期函数。 都可以分解为无穷个三角谐波的迭加，即： J．………………………………(2) sin(k09+rpt ，川=A+∑Ah j=J 上式无穷三角级数称为傅里叶级数。u，203,……nol各次三角波频率为1次、2次…… n次正弦谐波圆频率。 选择适当的函数作为作用函数，对原函数进行作用，并在某一范围内积分，这就是 所谓的函数卷积(faltung)处理，其解析数学表达式为： ^rr)o止ff J=J：：^f卜’{J，2f{彬善………………………………(3) 上式中原函数f。(t)与作用函数f2(t)都是时间变量t的函数，串是卷积运算的算符。 而式(3)在离散变量情况下有下面的表达式： ‘o#一m』Ⅷ 这里，矩阵P称之为原函数矩阵，矩阵C称之为卷积矩阵，或卷积模板。 选择不同的卷积模板，作用于原函数会对原函数产生不同的特殊效果。