云南省降水量時空分布圖

1、 本科學(xué)生綜合性、設(shè)計性 實驗報告姓名 趙藝淞 學(xué)號 104130138 專業(yè) 地理信息系統(tǒng) 班級 2010級地信班 實驗課程名稱 專題地圖編制 指導(dǎo)教師及職稱 何云燕 開課學(xué)期 2012 至 2013 學(xué)年第二 學(xué)期 云南師范大學(xué)旅游與地理科學(xué)學(xué)院編印一、實驗設(shè)計實驗名稱：云南省降水量動態(tài)分析圖實驗時間：2013年5月20日6月28日實驗類型： 驗證實驗 綜合實驗 設(shè)計實驗小組合作： 是 否小組成員：趙藝淞 尹文姣1、實驗內(nèi)容（含實驗原理介紹）：（1）研究背景：云南省全省129個縣降水量的時空分布分析（以2000年，2005年，2010年三年數(shù)據(jù)為例），分析比較同一年限降水量在整個云南省的空

2、間分布格局，不同年限降水量的差異和變化趨勢。（2）一幅完整的專題地圖由三個大的方面構(gòu)成，即專題地圖的數(shù)學(xué)要素、專題要素和地理底圖要素。專題地圖的數(shù)學(xué)要素包括坐標(biāo)網(wǎng)、比例尺和地圖定向等內(nèi)容。次幅專題地圖是反映人文經(jīng)濟(jì)現(xiàn)象的專題圖，在云南省行政區(qū)域范圍內(nèi)比較各個地州市的農(nóng)業(yè)發(fā)展和組成情況，采用較小比例尺1：oooo，在這里坐標(biāo)網(wǎng)為地理坐標(biāo)網(wǎng)也就是經(jīng)緯網(wǎng)。地圖的定向以中央經(jīng)線為正北方向。專題要素是專題地圖的主體。此次的專題圖采用了利用分區(qū)和分級統(tǒng)計方法配合進(jìn)行制圖，兩種方法配合使用，用分級統(tǒng)計圖作為背景，在圖上每一分區(qū)內(nèi)描繪統(tǒng)計圖表，從而使得地圖的主題內(nèi)容更加的全面。地理底圖要素：采用云南省行政區(qū)劃

3、圖作為地理底圖，人文經(jīng)濟(jì)專題圖主要是考慮研究區(qū)域的等面積性，和各研究區(qū)劃分正確性。適當(dāng)添加河流道路山峰等自然地理要素，可以幫讀圖者定位，原則是簡潔適當(dāng)，底圖顏色要淺淡，不能喧賓奪主。數(shù)學(xué)基礎(chǔ)：地圖投影是“高斯克呂格投影”，地理坐標(biāo)系統(tǒng)是“西安1980”。具體參數(shù)如下：地圖投影地理坐標(biāo)系統(tǒng)Projection: Gauss_KrugerFalse_Easting: 500000.000000False_Northing: 0.000000Central_Meridian: 102.000000Scale_Factor: 1.000000Latitude_Of_Origin: 0.000000Li

4、near Unit: Meter (1.000000)Geographic Coordinate System: GCS_Xian_1980Angular Unit: Degree (0.017453292519943299)Prime Meridian: Greenwich (0.000000000000000000)Datum: D_Xian_1980Spheroid: Xian_1980Semimajor Axis: 6378140.000000000000000000Semiminor Axis: 6356755.288157528300000000Inverse Flattening

5、: 298.257000000000010000地圖的地理基礎(chǔ)要素：基礎(chǔ)要素有各州市行政區(qū)劃邊界，河流，道理，山峰等地圖的專題內(nèi)容：專題是云南省2011年各州市農(nóng)林牧漁業(yè)總產(chǎn)值分布，顯示方式是以分區(qū)三維柱狀圖加分級設(shè)色顯示。圖例設(shè)計：圖例設(shè)計以出圖時排版為依據(jù)來設(shè)計圖面配置與整飾原則：圖面配置要考慮整幅圖的平衡問題，主圖要放人眼的視覺中心，即幾何中心偏左上位置，圖名，圖例，比例尺，附圖和說明圖表等要根據(jù)主圖的位置做出合理的配置。而地圖整飾要考慮各種顏色的搭配問題，使讓人看起來協(xié)調(diào)舒暢統(tǒng)計數(shù)據(jù)：數(shù)據(jù)是通過查閱云南省統(tǒng)計局的云南省統(tǒng)計數(shù)據(jù)庫，然后下載的到的， 地圖配準(zhǔn)數(shù)據(jù)：地圖配準(zhǔn)的坐標(biāo)是在網(wǎng)絡(luò)上

6、查找的一部分州市在西安1980坐標(biāo)系下的坐標(biāo)值。地圖的打印要求：本專題地圖采用A3紙打印2、實驗?zāi)康暮鸵螅簩ｎ}地圖所反映的主題內(nèi)容十分廣泛，從自然地理現(xiàn)象到社會經(jīng)濟(jì)等人文地理現(xiàn)象無所不容，既能表示各種具體有形的事物和現(xiàn)象又能表現(xiàn)抽象無形的事物現(xiàn)象。專題地圖還具有地理基礎(chǔ)和專題要素雙重內(nèi)容的層面的特性，它表現(xiàn)方法多樣，圖面配置靈活。學(xué)會專題地圖的制作對于我們GIS專業(yè)的學(xué)生是十分重要的。本次實驗的目的是進(jìn)一步加強(qiáng)對專題地圖的理解和認(rèn)識，特別是對專題地圖的表現(xiàn)方法和圖例設(shè)計的理解和操作。本次實驗是基于ArcGIS軟件的專題地圖制作，將所學(xué)知識付諸實踐，對鞏固專題地圖編制和地圖學(xué)的理論知識，加深對

7、ArcGIS軟件的理解和認(rèn)識，增強(qiáng)動手能力，提高制圖水平都有很大意義。3、實驗材料及相關(guān)設(shè)備：實驗數(shù)據(jù): （1）查閱云南統(tǒng)計年鑒獲取2000年，2005年，2010年云南省各縣的年均降水量數(shù)據(jù)。（2）云南省行政區(qū)劃圖一張，用于云南省省界的數(shù)字化。制圖軟件：ArcGIS軟件（ArcCatolog、ArcMap 、ArcView）相關(guān)輔助：Excel、Word等。4、理論依據(jù)（1）專題圖的優(yōu)勢：專題地圖是指突出而盡可能完善、詳盡地表示制圖區(qū)內(nèi)的一種或幾種自然或經(jīng)濟(jì)（人文）要素的地圖。專題地圖的制圖領(lǐng)域?qū)拸V，凡具有空間屬性的信息數(shù)據(jù)都可用其來表示。其內(nèi)容、形式多種多樣，能夠廣泛應(yīng)用與國民經(jīng)濟(jì)建設(shè)、教

8、學(xué)和科學(xué)研究、國防建設(shè)等行業(yè)部門。專題地圖和普通地圖相比，具有獨特的特征。（1）地統(tǒng)計（Geostatistics）又稱地質(zhì)統(tǒng)計，也可以稱為空間統(tǒng)計分析，其實統(tǒng)計學(xué)的一個分支，地統(tǒng)計學(xué)于20世紀(jì)50年代初開始形成，60年代在法國著名統(tǒng)計學(xué)家G.Matheron的大量理論研究工作基礎(chǔ)上，形成的一門新的統(tǒng)計學(xué)分支。ArcGIS地統(tǒng)計分析模塊主要由三個功能模塊組成：探索性數(shù)據(jù)分析、地統(tǒng)計分析向?qū)б约吧蓴?shù)據(jù)子集利用這些基本功能模塊，可以方便的完成多種地統(tǒng)計分析創(chuàng)建完善的專題地圖。（3）利用分區(qū)和分級統(tǒng)計方法配合進(jìn)行制圖分區(qū)統(tǒng)計圖法是在制圖區(qū)域各個區(qū)劃內(nèi)，按其相應(yīng)的統(tǒng)計數(shù)據(jù)，描繪不同形式的統(tǒng)計圖表，以

9、表示并比較各個區(qū)劃單位內(nèi)現(xiàn)象的總和、構(gòu)成及動態(tài)的方法。分級統(tǒng)計圖法是在整個制圖區(qū)域的若干個小的區(qū)劃單位內(nèi)（行政區(qū)劃或其他區(qū)劃單位），根據(jù)各分區(qū)資料的數(shù)量（相對）指標(biāo)進(jìn)行分級，并用相應(yīng)色級或不同疏密的暈線，反映各區(qū)現(xiàn)象的集中程度或發(fā)展水平的分布差別。分區(qū)統(tǒng)計圖表法與分級統(tǒng)計圖法均是以統(tǒng)計資料為基礎(chǔ)的表示方法。他們都能反映個區(qū)劃單位之間的數(shù)量差別，但不能反映每個區(qū)劃單元內(nèi)部的具體差異。兩種方法配合使用，用分級統(tǒng)計圖作為背景，在圖上每一分區(qū)內(nèi)描繪統(tǒng)計圖表，從而使得地圖的主題內(nèi)容更加的全面。4、實驗方法步驟及注意事項：數(shù)據(jù)準(zhǔn)備通過查閱2001年云南省統(tǒng)計年鑒、2006年云南省統(tǒng)計年鑒、2011年云南省

10、統(tǒng)計年鑒資料分別獲取2000年、2005年、2010年云南省129個縣區(qū)的年平均降水量數(shù)據(jù)，并制作成Excel表格的形式保存。下載JPG格式的云南省縣區(qū)行政區(qū)劃圖一份，并進(jìn)行數(shù)字化，并在屬性表添加縣名、所屬州市、周長、面積等屬性值，得到一幅云南省縣區(qū)行政區(qū)劃電子地圖。將Excel表格中的降水量數(shù)據(jù)鏈接導(dǎo)入到數(shù)字化好的地圖里，完成數(shù)據(jù)的前期準(zhǔn)備工作。步驟:地圖數(shù)字化，編制與成果輸出。地圖數(shù)字化以及屬性的添加（1）在ArcMap中添加底圖，建立地理坐標(biāo)系統(tǒng)，選控制點進(jìn)行配準(zhǔn)，選擇投影方式（2）在ArcCatolog創(chuàng)建所需圖層，面狀元素包括州市級的行政區(qū)劃、縣級行政區(qū)劃、云南省邊界，點狀元素包括各

11、測站點的位置。坐標(biāo)投影信息一定要和原底圖匹配。（3）創(chuàng)建好的圖層加載到ArcMap，調(diào)用Editor工具進(jìn)行編輯，根據(jù)底圖信息分別進(jìn)行數(shù)字化跟蹤，得到需要的各圖層數(shù)據(jù)并進(jìn)行拓?fù)錂z查（5）完善屬性信息。聯(lián)接屬性：將云南省降水量數(shù)據(jù)表在Excel中進(jìn)行編輯，將編輯好的數(shù)據(jù)另存為.DBF的數(shù)據(jù)格式。在屬性表中進(jìn)行屬性關(guān)聯(lián)，連接文件和被連接文件的關(guān)鍵字一定要一致。專題地圖編制（1）為數(shù)字化好的個圖層（各州市點，國界線等）設(shè)置顏色、大小。（2）利用分區(qū)和分級統(tǒng)計方法配合進(jìn)行制圖分級統(tǒng)計圖法是在整個制圖區(qū)域的若干個小的區(qū)劃單位內(nèi)（行政區(qū)劃或其他區(qū)劃單位），根據(jù)各分區(qū)資料的數(shù)量（相對）指標(biāo)進(jìn)行分級，并用相應(yīng)

12、色級或不同疏密的暈線，反映各區(qū)現(xiàn)象的集中程度或發(fā)展水平的分布差別。分區(qū)統(tǒng)計圖發(fā)是在制圖區(qū)域各個區(qū)劃內(nèi)，按其相應(yīng)的統(tǒng)計數(shù)據(jù)，描繪不同形式的統(tǒng)計圖表，以表示并比較各個區(qū)劃單位內(nèi)現(xiàn)象的總和、構(gòu)成及動態(tài)的方法。分區(qū)統(tǒng)計圖可以顯示現(xiàn)象的內(nèi)部結(jié)構(gòu)。兩種方法配合使用，用分級統(tǒng)計圖作為背景，在圖上每一分區(qū)內(nèi)描繪統(tǒng)計圖表，從而使得地圖的主題內(nèi)容更加的全面。本次實驗的數(shù)據(jù)是云南省2000年降水量的統(tǒng)計數(shù)據(jù)。分區(qū)統(tǒng)計圖表法與分級統(tǒng)計圖法均是以統(tǒng)計資料為基礎(chǔ)的表示方法。他們都能反映個區(qū)劃單位之間的數(shù)量差別，但不能反映每個區(qū)劃單元內(nèi)部的具體差異。先用分級統(tǒng)計法根據(jù)各州市的降水量進(jìn)行分級，對底色進(jìn)行分層設(shè)色。分級統(tǒng)計圖則

13、按照各個地州市的農(nóng)業(yè)總產(chǎn)值劃分為6個等級，并用不同的顏色來區(qū)分。然后再用分區(qū)統(tǒng)計法根據(jù)2000、2005、2010三年的降水量數(shù)據(jù)分別的產(chǎn)值生成柱狀圖。在這里分區(qū)統(tǒng)計圖以地州市級行政區(qū)域為劃分依據(jù)根據(jù)每個地州市在全省各地州年平均降水量的生成折線圖，作為附圖插入在地圖上，進(jìn)行補(bǔ)充說明。圖面設(shè)計點擊view-layout 添加各種地圖要素（1）添加地圖名稱，圖名圖廓作者日期（2）添加圖例（顏色） （3）添加指北針 （4）添加比例尺（5）對整圖進(jìn)行整飾 選擇landscape 進(jìn)行版面設(shè)計,對指北針、圖例、比例尺等各部分進(jìn)行修改調(diào)整。輸出地圖（1）先數(shù)值化云南省區(qū)劃圖（2）對數(shù)值化后的底圖進(jìn)行屬性賦

14、值（3）分區(qū)和分級統(tǒng)計方法配合進(jìn)行專題地圖編制 （4）利用Excel輔助成圖（5）圖面設(shè)計（6）輸出地圖二、實驗內(nèi)容、步驟和結(jié)果（要求：詳細(xì)寫清楚本次實驗的完成的主要內(nèi)容、具體實施步驟和實驗結(jié)果。紙張不夠可以自行添加。）1、地圖數(shù)字化以及屬性的添加1.1、數(shù)字化云南省縣區(qū)行政區(qū)劃地圖，得到一個命名為“云南省縣區(qū)行政區(qū)邊界”的面圖層，并在屬性表中添加字段名為“縣名-縮”、“縣名-全”“所屬州市”、“周長”、“面積”等屬性值。1.2、將2000年、2005年、2010年云南省129個縣區(qū)的三年的年平均降水量數(shù)據(jù)制作成Excel表格的形式保存，建立縣名的字段名為“縣名-縮”使之與屬性表的字段名相匹配

15、。1.3、打開面圖層“云南省縣區(qū)行政區(qū)劃圖”的屬性表，點擊Options按鈕，使用Join and Relates-Join工具，打開對話框，在打開的對話框里面選擇要Join的表，并選擇相對應(yīng)的字段。這里地圖數(shù)據(jù)使用的是“縣名-縮字段,Excel數(shù)據(jù)選擇的也是縣名-縮字段,因為它們含有相同的內(nèi)容. 點擊OK,然后再打開地圖數(shù)據(jù)的屬性表就可以看到Excel數(shù)據(jù)已經(jīng)導(dǎo)入到ArcGIS地圖數(shù)據(jù)里面了 . 因為Join后的屬性表只保存在內(nèi)存中，要是關(guān)掉程序再打開地圖數(shù)據(jù)，Excel里的數(shù)據(jù)又沒有了，所以需要將地圖數(shù)據(jù)導(dǎo)出保存。方法也是在地圖數(shù)據(jù)上右鍵，選擇Data-Export Data，選擇保存位置

16、即可。加載剛保存的數(shù)據(jù)，這時候Excel里的數(shù)據(jù)就已經(jīng)導(dǎo)入到地圖數(shù)據(jù)里面了。1.4、面圖層向點圖層轉(zhuǎn)化。地統(tǒng)計分析需要點圖層進(jìn)行內(nèi)插，從而得到連續(xù)的表面，我們現(xiàn)有數(shù)據(jù)為行政區(qū)界的面圖層，采用多邊形質(zhì)心的方法（面轉(zhuǎn)為幾何中心點）來把面圖層轉(zhuǎn)為點圖層。在ArcToolboxData Management Tools-Features-Feature to Point 得到如下結(jié)果1.5、將圖層“云南省各縣區(qū)行政邊界”合并成一個多邊形，只保留云南省省界。在ArcToolBox下選擇Data Management Tools- -Dissolve 。1.6、設(shè)置邊界可見，背景透明。單擊目錄表中的“云南

17、省各縣”圖層的圖例，打開Symbol Selector對話框。單擊Fill Color下拉箭頭，然后點擊No Color，在Symbol Selector對話框中單擊OK按鈕。使圖層酒杯透明顯示，只是有輪廓可見，以便看見在后面所要創(chuàng)建的其他圖層。2、基于空間插值方法，對云南省全省的平均降水量進(jìn)行分級2.1、激活地統(tǒng)計分析模塊。在ArcMap中，單擊Tools,再單擊Extensions,選中Geostatistical Analyst復(fù)選框，單擊Close按鈕。2.2、添加Geostatistical Analyst工具條到ArcMap中。單擊View菜單，光標(biāo)指向Toolbars，然后單擊G

18、eostatistical Analyst。2.3、在創(chuàng)建表面之前對數(shù)據(jù)進(jìn)行檢查，為了找出數(shù)據(jù)中那些離群值并且發(fā)現(xiàn)數(shù)據(jù)中存在的趨勢。用ESDA(空間數(shù)據(jù)探索分析)工具檢查數(shù)據(jù)我們對數(shù)據(jù)已經(jīng)有如下了解：2.3.1、降水量數(shù)據(jù)不符合正態(tài)分布，不過在經(jīng)歷了LOG數(shù)據(jù)變換之后，數(shù)據(jù)近似的符合了正態(tài)分布，平均值和中值大致相等，在直方圖中，變換后的數(shù)據(jù)呈單峰且在平均值/中值兩側(cè)有較好的對稱性。2.3.2、正態(tài)QQ圖中再次證實數(shù)據(jù)不符合正態(tài)分布，LOG變換后，QQ圖中的點構(gòu)成了一條近似的直線。3、利用Trend Analysis工具可以看到數(shù)據(jù)顯示出某種趨勢，還趨勢提取后可以是一個二階多項式對其進(jìn)行最佳擬合

19、。4、半變異函數(shù)表面表明數(shù)據(jù)中存在空間自相關(guān)。4.1、創(chuàng)建連續(xù)表面（克立格插值）結(jié)合上述數(shù)據(jù)分析發(fā)現(xiàn)的空間關(guān)系，根據(jù)降水量的樣本點，進(jìn)行空間內(nèi)插，創(chuàng)建剔除趨勢并模擬空間自相關(guān)的連續(xù)表面。Arcgis 的地統(tǒng)計分析有一個地統(tǒng)計向?qū)?，按照這個向?qū)б徊揭徊骄涂梢詫崿F(xiàn)空間插值值。單擊Geostatistical Analyst工具條，然后點擊地統(tǒng)計向?qū)eostatistical Wizard。 選擇輸入的數(shù)據(jù)和屬性，選擇的插值模型，這里選擇 kriging方法，如圖，然后點擊Next即可。選擇插值方法和數(shù)據(jù)預(yù)處理方法選擇克里金插值法。云南省各個縣區(qū)的平均降水量并不符合正態(tài)分布，而且在趨勢分析中存在明

20、顯的拋物線特征，“u”型曲線，所以這里選擇通過二階多項式擬合對其預(yù)處理，在Order of Trend Removal下拉箭頭選擇Second，比較插值結(jié)果的精度，然后 next。 這一步對插值的精度影響很大。根據(jù)左上角的半變異函數(shù)云圖，在右邊模型欄內(nèi)，選擇對其擬合較好的模型，通過可以通過調(diào)整搜尋角度和搜尋半徑（左下），最終確定一個合適的模型，該模型的參數(shù)在左下角藍(lán)色條框內(nèi)給出。單擊Finish按鈕，會首先給出這個過程的一個總結(jié) Summary。Output layer information 對話框顯示了用于創(chuàng)建表面的模型的信息摘要。例如，數(shù)據(jù)預(yù)處理過程，采用的模型等。4.2、得到插值的結(jié)果

21、，即圖層“Ordinary Kriging”將預(yù)測表面限制在云南省邊界內(nèi)。目錄表中右鍵點擊圖層“Ordinary Kriging”再單擊Properties,在Properties對話框中點擊General標(biāo)簽，將圖層名改為“云南省年平均降水量分布”由于預(yù)測到的降水分布圖并不能覆蓋云南省全部地理范圍（圖層“市區(qū)面積”），要解決這個問題，需要對表面進(jìn)行外推值（即預(yù)測邊界之外的點的值），點擊Extent標(biāo)簽，在Set the extent to框中選擇a custom extent entered below,在Visible Extent 區(qū)域中輸入下面的值，Left（左）：97,Right（右

22、）：107，Top（上）：30，Bottom（下）：21，然后點擊確定按鈕。在目錄中右擊Layer,然后單擊Properties，在對話框里選擇Data Frame標(biāo)簽，選中Enable復(fù)選框，而后點擊Specify Shape按鈕，在彈出對話框中單擊Outline of Features按鈕，在Layer下拉框中單擊“云南省邊界”。點擊OK,關(guān)閉窗口，完成整個裁剪過程。 最終得到裁剪后的分級圖如下圖所示：5、在每個地州市上生成柱狀圖，進(jìn)行分區(qū)顯示。5.1、將圖層“云南省各縣區(qū)行政邊界”按照“州市”字段合并成多邊形，只保留云南省各地州市邊界。在ArcToolBox下選擇Data Managem

23、ent Tools- -Dissolve 。改名為“云南省地州行政邊界”5.2、在Excel中，將2000年、2005年、2010年云南省129個縣的降水量數(shù)據(jù)以地州為單位進(jìn)行匯總，并建立新的各個地州市的三年的年平均降水量數(shù)據(jù)制作成Excel表格的形式保存，建立縣名的字段名為“州市”使之與屬性表的字段名相匹配。5.3、打開面圖層“云南省地州行政邊界”的屬性表，點擊Options按鈕，使用Join and Relates-Join工具，打開對話框，在打開的對話框里面選擇要Join的表，并選擇相對應(yīng)的字段。這里地圖數(shù)據(jù)使用的是“縣名-縮字段,Excel數(shù)據(jù)選擇的也是縣名-縮字段,因為它們含有相同的

24、內(nèi)容. 點擊OK,然后再打開地圖數(shù)據(jù)的屬性表就可以看到Excel數(shù)據(jù)已經(jīng)導(dǎo)入到ArcGIS地圖數(shù)據(jù)里面了 .5.4、將地圖數(shù)據(jù)導(dǎo)出保存。方法也是在地圖數(shù)據(jù)上右鍵，選擇Data-Export Data，選擇保存位置即可。加載剛保存的數(shù)據(jù)，這時候Excel里的數(shù)據(jù)就已經(jīng)導(dǎo)入到地圖數(shù)據(jù)里面了。5.5、右鍵點擊圖層“云南省地州行政邊界”，在下拉菜單中選擇Properties,在Symbology進(jìn)行如下選擇，制作一個反映每個地州三年的降水量變化的柱狀圖。5.6、將背景顏色設(shè)為無色，使得上一步得到的分級圖能夠顯示在圖上。6、國界線制作：將圖層“云南省各縣區(qū)行政邊界”合并成一個多邊形，只保留云南省省界。在

25、ArcToolBox下選擇Data Management Tools- -Dissolve。然后將得到的面狀圖層轉(zhuǎn)變?yōu)榫€狀圖層。在ArcToolboxData Management Tools-Features-Feature to line 制作邊界線符號：右鍵單擊該線狀圖層的符號，在打開的窗口Properyies中進(jìn)行設(shè)置，制作邊界線。7、添加昆明的位置，把云南省行政中心添加上8、添加地圖要素，完成專題地圖的輸出。在主菜單上單擊View下的Layout View，然后單擊Insert下的相應(yīng)模塊進(jìn)行添加，地圖顯示：最后將制作好的專題地圖轉(zhuǎn)化為JPG格式并輸出成圖實驗結(jié)果：參考文獻(xiàn)：1劉湘南，黃方，王平等.2008.GIS空間分析原理與方法(第二版)M.北京.科學(xué)出版社2宋小冬，鈕心毅.2007.地理信息系統(tǒng)實習(xí)教程（Arc GIS 9.x）（第一版）M.北京.科學(xué)出版社3黃杏元.馬勁松.2008.地理信息系統(tǒng)概論（第三版）M.北京.高等教育出版社4云南省統(tǒng)計局，國家統(tǒng)計局云南省調(diào)查總隊.2012.云南省統(tǒng)計年鑒 M中國統(tǒng)計出版社5黃仁濤，龐小平等.2003.專題地圖編制.(第一版)