应用MATLAB构造数字高程模型工具箱.pdf

第２３卷　第４期　　　　　　　　　　　　　 桂 林 工 学 院 学 报 　　　　　　　　　　　　Ｖｏｌ．２３Ｎｏ．４ ２００３年 １０月　　　　　　　　　ＪＯＵＲＮＡＬＯＦＧＵＩＬＩＮＩＮＳＴＩＴＵＴＥＯＦＴＥＣＨＮＯＬＯＧＹ　　　　　　　　Ｏｃｔ２００３  文章编号：１００６－５４４Ｘ（２００３）０４－０４５３－０４ 应用 ＭＡＴＬＡＢ构造数字高程模型工具箱 陈天伟 （桂林工学院土木工程系，广西 桂林　 ５４１００４） 摘　要：基于方格网的数字地面模型常面临不规则图幅边界处理等问题．设计了应用ＭＡＴ ＬＡＢ建立ＤＥＭ的算法、函数的调用格式及其工具箱界面．讨论了规则与不规则图幅边界处 理和等高线遇地物断开的问题：将不规则图幅边界和地物都看作多边形，求得各边界线段 经过的网格单元，相应编码为１，将边界外的单元编码为０；将地物轮廓线经过的单元编码 为２，其包含的单元编码为０，设编码为０的各单元的高程为非数ＮａＮ，则其等高线不会绘 出．对于图幅拼接，设置各图幅的 ‘Ｐｏｓｉｔｉｏｎ’属性进行定位即可．通过实例运行，验证了 方案的正确性． 关键词：工具箱；函数；插值；网格；编码 中图分类号：Ｐ２０８，Ｐ２０９　　　　　　　　　　文献标识码：Ａ① 　　ＭＡＴＬＡＢ作为一种大型数学软件，适合于多 ２　工具箱运行界面回调程序的设计 种学科，它包含了丰富的函数资源，采用直观的 符合人们思维习惯的代码，不但使人们解决各种 ２１　运行界面 学科的数学问题的编程简易可行，也使得开发工 　　运行界面包括命令窗和工具箱界面．工具箱 程软件的效率大大提高，开发周期大大缩短．随 界面利用ＭＡＴＬＡＢ的ＧＵＩ设计向导来设计，主要 着计算机技术的发展，数字高程模型 （ＤＥＭ）无 命令按钮有：高程数据、等高线、剖面图、三维 论是算法还是软件的实现都已有成果报道，本文 表面图、图例、退出、帮助信息，还有色彩菜单 就如何利用ＭＡＴＬＡＢ平台构造ＤＥＭ作一探讨． 及其他一些功能按钮．界面如图２所示． １　数字高程模型建立过程 　　 （１）读入数据：采用２种方式：文件输入或 屏幕数字化．屏幕数字化即用鼠标在打开的地图 ［１］ 影像上单击离散点 ，获得坐标并输入高程值． 　　 （２）划分格网：确定合适的坐标范围，由离 散点密度确定加密网格点． 　　 （３）调用Ｍａｔｌａｂ插值函数ｇｒｉｄｄａｔａ计算网格 点的高程，从而建立ＤＥＭ模型． 　　 （４）在此基础上，分别调用等高线、剖面 图、三维图等子程序输出相应信息． 图１　建立ＤＥＭ的流程图 Ｆｉｇ１　ＦｌｏｗｃｈａｒｔｏｆｃｏｎｓｔｒｕｃｔｉｎｇＤＥＭ 　　工作流程如图１所示． ① 收稿日期：２００３－０５－１５；修订日期：２００３－０７－０８ 基金项目：广西区科技攻关项目 （桂科攻０３３０００８－７） 作者简介：陈天伟 （１９６５－），男，广西北流人，讲师，工程测量专业． ４５４ 桂　林　工　学　院　学　报　　　　　　　　　　　　　　　２００３年 图２　工具箱界面和输出等高线 Ｆｉｇ２　Ｉｎｔｅｒｆａｃｅｏｆｔｏｏｌｂｏｘａｎｄｅｘｐｏｒｔｃｏｎｔｏｕｒ ２２　回调程序的算法设计 高程，格式为：［Ｃ，ｈ］＝ｃｏｎｔｏｕｒ（ｘｘ，ｙｙ，ｚｚ，［ｊ　 ［２］ ２２１　网格点高程插值 　 （１）获得离散点信 ｊ］）．其中ｘｘ，ｙｙ，ｚｚ为网格点坐标高程数组；ｊ为等 息后，调用Ｌｉｎｓｐａｃｅ函数计算网格的ｘ坐标数组， 高线高程值，ｈ为该等高线句柄．Ｃ表示在鼠标单击 格式为：