基于Matlab的空间分析

2012年4月地理空间信息Apr, 2012第10卷第2期GEOSPATIAL INFORMATIONVol. 10,No.2基于Matlab的空间分析胡圣武(河南理工大学测绘学院，河南焦作454000)摘要:空间分析被认为是地理信息系统区别于一般信息系统、CAD或电子地图系统的主要功能特征，也是评价一个地理信息系统的主要指标。从介绍Matlab强大的数值计算功能和数据可视化功能着手，分析了Matlab在空间分析方面的技术应用，重点研究了Matlab 在空间数据分析中插值法的应用以及空间数据的三维表达。关键词: Matlab; 空间分析;插值分析;三维表达中图分类号: P208文献标志码: B文章编号: 1672-4623 (2012) 02-0022-03空间分析是GIS的核心和重要功能之一, GIs具1.1 Matlab 的特点有强大的空间分析功能，是区别计算机制图系统和数Matlab的特点比较多，其主要特点为D!:①友好据库管理系统的显著特征之一n。空间分析使GIS不的工作平台和编程环境;②简单易用的程序语言;③仅体现在地图制图上，用户还可以通过与系统交互而强大的科学计算机数据处理能力;④出色的图形处理将地理数据经过分析转换为对自己有用的信息。同时，功能;⑤应用广泛的模块集合工具箱;⑥实用的程序利用空间分析技术，对原始数据模型进行观察和实验，接口和发布平台;⑦应用软件开发(包括用户界面)。用户可以获得新的经验和知识，并以此为空间行为的1.2 Matlab 用于空间分析的优点决策依据。空间分析对空间信息的提取和传输，已使1)强大的空间模型成图能力。Matlab能够将计算GIS成为区别于一般信息系统的主要功能特征，也是的结果以二维、三维乃至四维的图形表现出来，并且评价GIS功能强弱的重要指标之一。可为用户提供灵可以通过对图形线性、立面、色彩、渲染、光线、视活解决各类专门问题的有效工具。随着空间分析技术角等特性的处理，把计算结果表现得淋漓尽致。的不断发展,GIS也将从一-般的空间事务处理向空间决2)强大的图像处理能力。图像处理工具包是由一策支持方向迈进”。系列支持图像处理操作的函数组成的，所支持的图像处理操作有:几何操作、区域操作和块操作，线性滤1 Matlab 理论波和滤波器设计，变换(DCT变换),图像分析和增Matlab是由美国Mathworks公司发布的主要面对强，二值图像操作等。科学计算、可视化以及交互式程序设计的高科技计算Matlab不但提供非常方便的绘图功能和强大的图环境。它将数值分析、矩阵计算、科学数据可视化以形图像处理能力，而且提供了增强的交互处理能力。颜及非线性动态系统的建模和仿真等诸多强大功能集成色映射、线形、标记、颜色、图例、坐标轴、多窗口.在一个易于使用的视窗环境中，为科学研究、工程设制图根据用户要求都可以定制;三维图形可以缩放、旋计以及必须进行有效数值计算的众多科学领域提供了转和改变视角4.一种全面的解决方案，并在很大程度上摆脱了传统非3)丰富的库函数。工具箱"tolbox” 是Matlab 扩交互式程序设计语言的编辑模式，代表了当今国际科展自身功能的开放式的开发模式。map和image是Matlab学计算软件的先进水平。提供的2个与GIS密切相关的工具箱。Matlab的基本数据单位是矩阵，它的指令表达式2基于Matlab的空间插值分析与数学、工程中常用的形式十分相似，故用Matlab来解算问题要比用C、FORTRAN等语言完成相同的事情本文重点研究了二维插值空间分析方法，在插值简捷得多，并且Matlab也吸收了像Maple等软件的优广]中国煤花士以某种方法描述数据点，使Matlab成为-个强大的数学软件.点之间TMHCNMH了4种插值方法，在分析成西双不后，叮以丈对地用这几种插值方法。收稿日期: 2011-08-22项目来源:国家白然科学基金资助项目(40474003)。第10卷第2期胡圣武:基于Matlab的空间分析2.1数据准备要在某山区方圆大约27 000m2的范围内修一条公路,从山下经过一个居民区，然后通向一个矿区。横向纵向区城分别每隔400 m测量- -次，得到一些点的高程，需要用Matlab做出地貌图和等高线，进行空间分析，表1为空间数据表51.表1空间数据表图3等高线图1200 1600 2000 2400 2800 3200 3600 4 0002.34种插值方法的演示比较1200 1130 1250 1280 1230 1040 90500700以上数据演示了空间数据在Matlab的实现，但对于1600 1320 1490 1420 1400 1300 709085空间插值的演示，可以通过interp2命令绘图比较插值的2000 1390 1500 1500 1400 900 1100 1060 9504种方法在效果上的差异。interp2 (x, y, z, xi, yi, method)2400 1500 1200 1 100 1350 1450 1 2001 150 1 010可以对粗糙的表面插值光滑，Method有 4种方法:1) 'nearest': 最近点插值法，它根据巳知两点间2800 1500 1200 1100 1550 1600 1550 1380 1 070320015001550160015501600160016001550的插值点和这两点间位置的远近进行插值。2) 'linear': 线性插值法，它把相邻的数据点用直3600 1480 1500 1550 1540 1450 1300 1200 980线相连，按所生成的直线进行插值。2.2 Matlab 的实现) 'spline' :样条插值法，它用已知数据求出样Matlab实现首先需要规定数据的输入条函数后，按照样条函数进行插值。>>x=1200: 400: 4000; y=1200: 400: 3600;4) 'cubic': 立方插值法，它根据已知数据求出立>>z= [ 1130 1250 1280 1230 1040 900 500 700;方数据立方函数后，按照立方函数进行插值。1320 1490 1420 1400 1300 700 900 850;对图2进行插值处理，我们分别按照4种方式进1390 1500 1500 1400 900 1100 1060 950;行(结果如图4所示)，命令如下:1500 1200 1100 1350 1450 1200 1150 1010;[xi, yi] =meshgrid (1200: 100: 4000,1 200: 100:1500 1200 1100 1550 1600 1550 1380 1070;3 600);1500 1550 1600 1550 1600 1500 1600 1550;>>zi1=interp2 (x, y, z, xi, yi,' neareast' );1480 1500 1550 1540 1450 1300 1200 980];>>zi2=interp2 (x, y, z, xi, yi,' linear' );然后在Matlab中键入:>>zi3=interp2 (x, y, z, xi, yi,’ spline' );>>mesh (x, y, z);得到网格图(见图1)>>zi4=interp2 (x, y, z, xi, yi' cubic' );>>surf (x, y, z);得到曲面图(见图2)a)最近点插值法效果图b)线性插值法效果图图1网格图图2曲面图等高线在GIS领域中有着很重要的地位，Matlab中国煤化工强大的数学分析功能使等高线的形成变得简单，contourHCNMHG语句可以很快地形成等高线(如图3所示)。在Matlabc)样条插值法效果图d)立方条插值法效果图中键入:图44种效果图的比较>>contour (x, y, z); %形 成等高线由图4不难看出，插值后的图像要比插值前图像24.地理空间信息第10卷第2期表面光滑得多，但比较上面4个插值效果图可以发现，U, V),它输出的参量X, Y是数据点的坐标，U和V是立方插值法和样条插值法生成的图形确实要比双线性相应数据点处的想x, y偏导数值，quiver 最后绘制出插值法和最近点插值法图形效果好得多。在应用中，可的梯度向量箭头就是根据函数gradient求得的偏导数求以根据不同的插值方法的优劣以及对插值结果的要求，得的，生成梯度图如图8所示。选择不同的插值方法。上大13基于Matlab的三维数据表达Matlab不仅具有强大的数值计算功能，还可以利用二维和三维绘图功能，将运算结果可视化，将杂乱的数据通过图形表示，从中观察出其内在关系，图形命令格式简单，可以使用不同的线性、色彩、数据点标记和标注等来修饰图形。许多GIS都提供专业的三维分析功能，但使用和支持开发功能很弱1an]。在Matlab图7带基准平面的网格图图8梯度图中分析三维数据，不但提供很多现成的函数，而且其键入命令为:灵活运用矩阵运算的方法，可以实现许多其他软件无》[x, y] =meshgrid (-3: 0.5: 3);法实现的功能l7.1。>> z-peaks (x, y); surf (x, y, z);》[px, y] =gradient (z, 0.5,0.5);3.1三维图的生成Matlab提供函数peaks来生成高斯分布的数据。首》contour (x, y, z);> hold on;先生成原始数据，并绘制数据图形，其中用到了函数surf,绘制的图形如图5所示。》quiver (x, y, px, py);利用Matlab来实现空间模型，不仅使空间模型的>> [x, y] =meshgrid (-3: 0.5: 3);成图、转换、空间分析变得可行又简单，而且能更加>>z=peaks Gx, ); surf (x, y, z);“meshc (x, y, z)”生成具有基本等高线的网格图，清晰地表达成果。“meshz (x, y, z)” 生成带有基准面的网格图。输入以参考文献下命令，可以得到如图6所示的带基本等高线和如图[1]邬伦,刘瑜,张晶,等.地理信息系统:原理、方法和应用[M].北7所示的带基准面的平面网格图。京:科学出版社,2005>>meshc (x, y, z); %绘制带基 本等高线的网格图2] 王劲峰,廖一兰,刘鑫.空间数据分析教程[M].北京:科学出版社,2010>>meshz(x,y,z);%绘制带基准平面的网格图3] Duane Hanselman,Bruce Litefield著, 朱仁峰译.精通Matlab7[M].北京:清华大学出版社,20064] 曾文华,钭春红.Matlab在地理信息系统中的应用[]计算机.应用,2022);91-9395] 崔利,武文波.测绘领城中Matlab的应用[].辽宁工程技术大学学报2005,24(SI):47-496] 杨驰.GIS空间分析建模构想[].测绘通报,06011:22227] 黄光东,管建和,李响,等.基于Matlab平台的插值法技术实现与应用[J].电脑编程技巧与维护,2006,2:30-328]田雨,卢秀山,独知行,等.GIS空间分析实验中的Matlab应用图5三维高斯分布的原始空间图6带基本等高线 的.[]地理空间信息,2009, 7()1113数据图网格图3.2梯度图的生成作者简介:胡圣武，博士，副教授,现主要从事GIS基础理论在Matlab中产生梯度箭头的函数是quiver (X, Y,和图像处理技术研究。中国煤化工MHCNMHG