2012年4月地理空間信息Apr, 2012第10卷第2期GEOSPATIAL INFORMATIONVol. 10,No.2基于Matlab的空間分析胡圣武(河南理工大學(xué)測(cè)繪學(xué)院，河南焦作454000)摘要:空間分析被認(rèn)為是地理信息系統(tǒng)區(qū)別于一般信息系統(tǒng)、CAD或電子地圖系統(tǒng)的主要功能特征，也是評(píng)價(jià)一個(gè)地理信息系統(tǒng)的主要指標(biāo)。從介紹Matlab強(qiáng)大的數(shù)值計(jì)算功能和數(shù)據(jù)可視化功能著手，分析了Matlab在空間分析方面的技術(shù)應(yīng)用，重點(diǎn)研究了Matlab 在空間數(shù)據(jù)分析中插值法的應(yīng)用以及空間數(shù)據(jù)的三維表達(dá)。關(guān)鍵詞: Matlab; 空間分析;插值分析;三維表達(dá)中圖分類(lèi)號(hào): P208文獻(xiàn)標(biāo)志碼: B文章編號(hào): 1672-4623 (2012) 02-0022-03空間分析是GIS的核心和重要功能之一, GIs具1.1 Matlab 的特點(diǎn)有強(qiáng)大的空間分析功能，是區(qū)別計(jì)算機(jī)制圖系統(tǒng)和數(shù)Matlab的特點(diǎn)比較多，其主要特點(diǎn)為D!:①友好據(jù)庫(kù)管理系統(tǒng)的顯著特征之一n。空間分析使GIS不的工作平臺(tái)和編程環(huán)境;②簡(jiǎn)單易用的程序語(yǔ)言;③僅體現(xiàn)在地圖制圖上，用戶還可以通過(guò)與系統(tǒng)交互而強(qiáng)大的科學(xué)計(jì)算機(jī)數(shù)據(jù)處理能力;④出色的圖形處理將地理數(shù)據(jù)經(jīng)過(guò)分析轉(zhuǎn)換為對(duì)自己有用的信息。同時(shí)，功能;⑤應(yīng)用廣泛的模塊集合工具箱;⑥實(shí)用的程序利用空間分析技術(shù)，對(duì)原始數(shù)據(jù)模型進(jìn)行觀察和實(shí)驗(yàn)，接口和發(fā)布平臺(tái);⑦應(yīng)用軟件開(kāi)發(fā)(包括用戶界面)。用戶可以獲得新的經(jīng)驗(yàn)和知識(shí)，并以此為空間行為的1.2 Matlab 用于空間分析的優(yōu)點(diǎn)決策依據(jù)?？臻g分析對(duì)空間信息的提取和傳輸，已使1)強(qiáng)大的空間模型成圖能力。Matlab能夠?qū)⒂?jì)算GIS成為區(qū)別于一般信息系統(tǒng)的主要功能特征，也是的結(jié)果以二維、三維乃至四維的圖形表現(xiàn)出來(lái)，并且評(píng)價(jià)GIS功能強(qiáng)弱的重要指標(biāo)之一?？蔀橛脩籼峁╈`可以通過(guò)對(duì)圖形線性、立面、色彩、渲染、光線、視活解決各類(lèi)專(zhuān)門(mén)問(wèn)題的有效工具。隨著空間分析技術(shù)角等特性的處理，把計(jì)算結(jié)果表現(xiàn)得淋漓盡致。的不斷發(fā)展,GIS也將從一-般的空間事務(wù)處理向空間決2)強(qiáng)大的圖像處理能力。圖像處理工具包是由一策支持方向邁進(jìn)”。系列支持圖像處理操作的函數(shù)組成的，所支持的圖像處理操作有:幾何操作、區(qū)域操作和塊操作，線性濾1 Matlab 理論波和濾波器設(shè)計(jì)，變換(DCT變換),圖像分析和增Matlab是由美國(guó)Mathworks公司發(fā)布的主要面對(duì)強(qiáng)，二值圖像操作等?？茖W(xué)計(jì)算、可視化以及交互式程序設(shè)計(jì)的高科技計(jì)算Matlab不但提供非常方便的繪圖功能和強(qiáng)大的圖環(huán)境。它將數(shù)值分析、矩陣計(jì)算、科學(xué)數(shù)據(jù)可視化以形圖像處理能力，而且提供了增強(qiáng)的交互處理能力。顏及非線性動(dòng)態(tài)系統(tǒng)的建模和仿真等諸多強(qiáng)大功能集成色映射、線形、標(biāo)記、顏色、圖例、坐標(biāo)軸、多窗口.在一個(gè)易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)制圖根據(jù)用戶要求都可以定制;三維圖形可以縮放、旋計(jì)以及必須進(jìn)行有效數(shù)值計(jì)算的眾多科學(xué)領(lǐng)域提供了轉(zhuǎn)和改變視角4.一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非3)豐富的庫(kù)函數(shù)。工具箱"tolbox” 是Matlab 擴(kuò)交互式程序設(shè)計(jì)語(yǔ)言的編輯模式，代表了當(dāng)今國(guó)際科展自身功能的開(kāi)放式的開(kāi)發(fā)模式。map和image是Matlab學(xué)計(jì)算軟件的先進(jìn)水平。提供的2個(gè)與GIS密切相關(guān)的工具箱。Matlab的基本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式2基于Matlab的空間插值分析與數(shù)學(xué)、工程中常用的形式十分相似，故用Matlab來(lái)解算問(wèn)題要比用C、FORTRAN等語(yǔ)言完成相同的事情本文重點(diǎn)研究了二維插值空間分析方法，在插值簡(jiǎn)捷得多，并且Matlab也吸收了像Maple等軟件的優(yōu)廣]中國(guó)煤花士以某種方法描述數(shù)據(jù)點(diǎn)，使Matlab成為-個(gè)強(qiáng)大的數(shù)學(xué)軟件.點(diǎn)之間TMHCNMH了4種插值方法，在分析成西雙不后，叮以丈對(duì)地用這幾種插值方法。收稿日期: 2011-08-22項(xiàng)目來(lái)源:國(guó)家白然科學(xué)基金資助項(xiàng)目(40474003)。第10卷第2期胡圣武:基于Matlab的空間分析2.1數(shù)據(jù)準(zhǔn)備要在某山區(qū)方圓大約27 000m2的范圍內(nèi)修一條公路,從山下經(jīng)過(guò)一個(gè)居民區(qū)，然后通向一個(gè)礦區(qū)。橫向縱向區(qū)城分別每隔400 m測(cè)量- -次，得到一些點(diǎn)的高程，需要用Matlab做出地貌圖和等高線，進(jìn)行空間分析，表1為空間數(shù)據(jù)表51.表1空間數(shù)據(jù)表圖3等高線圖1200 1600 2000 2400 2800 3200 3600 4 0002.34種插值方法的演示比較1200 1130 1250 1280 1230 1040 90500700以上數(shù)據(jù)演示了空間數(shù)據(jù)在Matlab的實(shí)現(xiàn)，但對(duì)于1600 1320 1490 1420 1400 1300 709085空間插值的演示，可以通過(guò)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可以對(duì)粗糙的表面插值光滑，Method有 4種方法:1) 'nearest': 最近點(diǎn)插值法，它根據(jù)巳知兩點(diǎn)間2800 1500 1200 1100 1550 1600 1550 1380 1 070320015001550160015501600160016001550的插值點(diǎn)和這兩點(diǎn)間位置的遠(yuǎn)近進(jìn)行插值。2) 'linear': 線性插值法，它把相鄰的數(shù)據(jù)點(diǎn)用直3600 1480 1500 1550 1540 1450 1300 1200 980線相連，按所生成的直線進(jìn)行插值。2.2 Matlab 的實(shí)現(xiàn)) 'spline' :樣條插值法，它用已知數(shù)據(jù)求出樣Matlab實(shí)現(xiàn)首先需要規(guī)定數(shù)據(jù)的輸入條函數(shù)后，按照樣條函數(shù)進(jìn)行插值。>>x=1200: 400: 4000; y=1200: 400: 3600;4) 'cubic': 立方插值法，它根據(jù)已知數(shù)據(jù)求出立>>z= [ 1130 1250 1280 1230 1040 900 500 700;方數(shù)據(jù)立方函數(shù)后，按照立方函數(shù)進(jìn)行插值。1320 1490 1420 1400 1300 700 900 850;對(duì)圖2進(jìn)行插值處理，我們分別按照4種方式進(jìn)1390 1500 1500 1400 900 1100 1060 950;行(結(jié)果如圖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);得到網(wǎng)格圖(見(jiàn)圖1)>>zi4=interp2 (x, y, z, xi, yi' cubic' );>>surf (x, y, z);得到曲面圖(見(jiàn)圖2)a)最近點(diǎn)插值法效果圖b)線性插值法效果圖圖1網(wǎng)格圖圖2曲面圖等高線在GIS領(lǐng)域中有著很重要的地位，Matlab中國(guó)煤化工強(qiáng)大的數(shù)學(xué)分析功能使等高線的形成變得簡(jiǎn)單，contourHCNMHG語(yǔ)句可以很快地形成等高線(如圖3所示)。在Matlabc)樣條插值法效果圖d)立方條插值法效果圖中鍵入:圖44種效果圖的比較>>contour (x, y, z); %形 成等高線由圖4不難看出，插值后的圖像要比插值前圖像24.地理空間信息第10卷第2期表面光滑得多，但比較上面4個(gè)插值效果圖可以發(fā)現(xiàn)，U, V),它輸出的參量X, Y是數(shù)據(jù)點(diǎn)的坐標(biāo)，U和V是立方插值法和樣條插值法生成的圖形確實(shí)要比雙線性相應(yīng)數(shù)據(jù)點(diǎn)處的想x, y偏導(dǎo)數(shù)值，quiver 最后繪制出插值法和最近點(diǎn)插值法圖形效果好得多。在應(yīng)用中，可的梯度向量箭頭就是根據(jù)函數(shù)gradient求得的偏導(dǎo)數(shù)求以根據(jù)不同的插值方法的優(yōu)劣以及對(duì)插值結(jié)果的要求，得的，生成梯度圖如圖8所示。選擇不同的插值方法。上大13基于Matlab的三維數(shù)據(jù)表達(dá)Matlab不僅具有強(qiáng)大的數(shù)值計(jì)算功能，還可以利用二維和三維繪圖功能，將運(yùn)算結(jié)果可視化，將雜亂的數(shù)據(jù)通過(guò)圖形表示，從中觀察出其內(nèi)在關(guān)系，圖形命令格式簡(jiǎn)單，可以使用不同的線性、色彩、數(shù)據(jù)點(diǎn)標(biāo)記和標(biāo)注等來(lái)修飾圖形。許多GIS都提供專(zhuān)業(yè)的三維分析功能，但使用和支持開(kāi)發(fā)功能很弱1an]。在Matlab圖7帶基準(zhǔn)平面的網(wǎng)格圖圖8梯度圖中分析三維數(shù)據(jù)，不但提供很多現(xiàn)成的函數(shù)，而且其鍵入命令為:靈活運(yùn)用矩陣運(yùn)算的方法，可以實(shí)現(xiàn)許多其他軟件無(wú)》[x, y] =meshgrid (-3: 0.5: 3);法實(shí)現(xiàn)的功能l7.1。>> z-peaks (x, y); surf (x, y, z);》[px, y] =gradient (z, 0.5,0.5);3.1三維圖的生成Matlab提供函數(shù)peaks來(lái)生成高斯分布的數(shù)據(jù)。首》contour (x, y, z);> hold on;先生成原始數(shù)據(jù)，并繪制數(shù)據(jù)圖形，其中用到了函數(shù)surf,繪制的圖形如圖5所示。》quiver (x, y, px, py);利用Matlab來(lái)實(shí)現(xiàn)空間模型，不僅使空間模型的>> [x, y] =meshgrid (-3: 0.5: 3);成圖、轉(zhuǎn)換、空間分析變得可行又簡(jiǎn)單，而且能更加>>z=peaks Gx, ); surf (x, y, z);“meshc (x, y, z)”生成具有基本等高線的網(wǎng)格圖，清晰地表達(dá)成果。“meshz (x, y, z)” 生成帶有基準(zhǔn)面的網(wǎng)格圖。輸入以參考文獻(xiàn)下命令，可以得到如圖6所示的帶基本等高線和如圖[1]鄔倫,劉瑜,張晶,等.地理信息系統(tǒng):原理、方法和應(yīng)用[M].北7所示的帶基準(zhǔn)面的平面網(wǎng)格圖。京:科學(xué)出版社,2005>>meshc (x, y, z); %繪制帶基 本等高線的網(wǎng)格圖2] 王勁峰,廖一蘭,劉鑫.空間數(shù)據(jù)分析教程[M].北京:科學(xué)出版社,2010>>meshz(x,y,z);%繪制帶基準(zhǔn)平面的網(wǎng)格圖3] Duane Hanselman,Bruce Litefield著, 朱仁峰譯.精通Matlab7[M].北京:清華大學(xué)出版社,20064] 曾文華,鈄春紅.Matlab在地理信息系統(tǒng)中的應(yīng)用[]計(jì)算機(jī).應(yīng)用,2022);91-9395] 崔利,武文波.測(cè)繪領(lǐng)城中Matlab的應(yīng)用[].遼寧工程技術(shù)大學(xué)學(xué)報(bào)2005,24(SI):47-496] 楊馳.GIS空間分析建模構(gòu)想[].測(cè)繪通報(bào),06011:22227] 黃光東,管建和,李響,等.基于Matlab平臺(tái)的插值法技術(shù)實(shí)現(xiàn)與應(yīng)用[J].電腦編程技巧與維護(hù),2006,2:30-328]田雨,盧秀山,獨(dú)知行,等.GIS空間分析實(shí)驗(yàn)中的Matlab應(yīng)用圖5三維高斯分布的原始空間圖6帶基本等高線 的.[]地理空間信息,2009, 7()1113數(shù)據(jù)圖網(wǎng)格圖3.2梯度圖的生成作者簡(jiǎn)介:胡圣武，博士，副教授,現(xiàn)主要從事GIS基礎(chǔ)理論在Matlab中產(chǎn)生梯度箭頭的函數(shù)是quiver (X, Y,和圖像處理技術(shù)研究。中國(guó)煤化工MHCNMHG