matlab分析地形,基于Matlab的地形数据处理

基于Matlab的地形数据处理

第 3 0卷第 1 期 2 0 1 0年 3月 测绘科学 与工 程 Ge o ma t i c S c i e nc e a n d En g i n e e r i ng Vo1 ． 3 0． No． 1 Ma r ． 2 01 0 基于 Ma t l a b的地形数据处理 杨韫澜① 胡海彦 陈 虹 江振 治 【 摘 要】 地形特征点散点、 规则网格、 三角网格及等高线等多种地形数据的存储形式及数据结构特殊 性很适合 m a t l a b的数据 处理 特点 。本 文研究如何 在 m a t l a b中通过 简单 的编程 ， 从 而对 地形数据 进行 必要 的 处理 ， 以满足地形测量要求 ， 并给出了数据处理实例 。 【 关键词】 角网格 ； 规则网格； 等高线 ； 内插； 曲线拟合； Ma t l a b 【 分类号】 1 2 3 1 El e v a t i o n Da t a Pr o c e s s i n g ba s e d o n M a t l a b Ya n g Yu n l a n Hu Ha i y a n Ch e n Ho ng J i a n g Z he n z h i Abs t r a c t Th e s t o r a g e for ms o f t e r r a i n d a t a s u c h a s c ha r a c t e r i s t i c p o i n t da t a，g r i ds， t r i a n g u l a t i o n n e t ，e t c ．，a nd t h e s p e c i a l d a t a s t r uc t u r e a r e q ui t e s u i t e d t o b e pr o c es s e d i n ma t l a b．Thi s pa pe r s t u di e s t h e me t h o d me e t i n g t he r e — q u i r eme n t s o f t e r r a i n me a s u r e me n t f o r t e r r a i n d a t a p r o c e s s i ng i n Ma t l a b b y s i mpl e p r o g r a mmi n g ． Fi na l l y，t h e p a p e r g i v e s t h e e x a mp l e o f d a t a p r o c e s s i n g． Ke y wo r ds t r i a ng l a t i on n e t wo r k s； g r l ds ； c o n t o ur l i n e s； i n t e r p o l a t i o n； C H I V e fit t i n g； ma t l a b 1 引言 M a t l a b是适合多学科 、 多种工作平 台的功能 强大 、 界面友好且开放性很强 的大型优秀应用软 件， 是数值分析 、 数字信号处理 、 自动控制及工程 应用的有效试验仿真工具 j 。Ma t l a b中的所有数 据都是用数组表示 和存储 的， 并提供 数组运算功 能用于对各种数据集合进行处理 ， 这种特点很符 合对地形数据的处理。本文讨论如何运用 Ma t l a b 对地形特征散点数据 、 规则网格 、 三角网及等高线 等地形数据的具体处理方法， 其 主要技术环节包 括 ： 散点数据构建 D e l a u n a y三角网格 ； 规则 网格 的内插 ； 等高线绘制等。整个试验演示所用数据 为陕北某 区域 ( 大型冲沟 ) 实测地形数据 ， 散点为 1 5 2 9 3 0个 ， 在全数字工作站上通过立体观测和影 像匹配而得。 ①测绘信息技术总站工程师 图 1 地 形 特 征 散 列点 数 据 2 散点数据构建 D e l a u n a y 三角网 地形特征散点数据是最 “ 原始 ” 的地形表现 形式 ， 它一般通过立体像对 的解析和影像 匹配而 获得散点的( X ， Y ， z ) 坐标。由于这些散点坐标数 据的排列是随机的， 规律性差 ， 所 以直接应用较为 3 0 杨韫 澜 、胡海彦、 陈虹、 江振治： 基 于 M a t 1 a b 的地形数据处理 第 1 期 复杂 ’ 。通 常 的 做 法 是 对 该 数 据 集 构 建 D e l a u n a y 三角网格 ， 而 D e l a u n a y三角网格是表示 地形信息公认 的数据冗余度最小、 表现力最强 的 地形数据形式之一。 2 ． 1 构建 D e l a u n a y三角网 D e l a u n a y方法是采用一组符合 D e l a u n a y原则 的三角形将所有 的数据点 连接起来 ， 并且没有任 何一个散点落在任何一个三角形之内。在 M a t l a b 中， 函数 d e l a u n a y ( X ， Y ) 返回一列数据索引用于标 明各个三角形 的顶点。例如 ， 导入散点数据后 ， 试 验中返 回如下结果 ： 》 t r i =d e l a u n a y ( x ， y ) t r i = 图 2 Ma t l a b中的 De l a u n a y函数 其中， 每一行 由三角形 的 3个 顶点 的 x和 Y 的索 引构 成。例 如， 第一 个 三角 形 是 以位 于 X [ 9 6 3 6 ， 3 8 4 2， 3 8 4 1 ] 和 Y [ 9 6 3 6 ， 3 8 4 2 ， 3 8 4 1 ] 的三个 数据点为顶点构成 的。继而可 以用 函数 t r i m e s h 将这些三角形绘制出来 ： 》 h o l d o n ， t r i me s h ( t r i ， x ， Y ， z ) ， h o l d o ff 》h i d d e n o ff 》 t i t l e ( I ) e l a u n a y T ri a n g u l a t i o n ) 2 ． 2 D e l a u n a y三角网的应 用 获得 D e l a u n a y三角网格 t r i 之后 ， 即可以用函 数 t s e a r c h和 d s e a r c h对散列数据进行插值 。例如 可 以用 t s e a r c h函数 查 看平 面 坐标 ( 4 ． 7 6 6 3 e+ 0 0 5 ， 3 ． 4 3 0 8 e+0 0 6 ) 位于那个三角形 内， 代码如 下 ： 》 t s e a r e h ( x ， Y ， t r i ， 4 ． 7 6 6 3 e+0 0 5 ， 3 ． 4 3 0 8 e+ 0 0 6) 图 3 d e l a u n a y三角 网绘 制 ( 顶 视 图 ) anS = 1 5 0 9 3 》 t r i ( a n s ， ： ) anS