地理信息系統(tǒng)實(shí)習(xí)教程第九章

本文格式為Word版，下載可任意編輯——地理信息系統(tǒng)實(shí)習(xí)教程第九章第九章再分類與疊合計(jì)算

1中學(xué)選址的依據(jù)

本練習(xí)要求對現(xiàn)有中學(xué)、人口分布、土地使用進(jìn)行評定，為新建中學(xué)的選址提供依據(jù)，可利用的數(shù)據(jù)有三種：

（1）現(xiàn)有中學(xué)。點(diǎn)狀空間數(shù)據(jù)，新建中學(xué)不應(yīng)離現(xiàn)有中學(xué)太近，為此建立離開現(xiàn)有中學(xué)的距離圖。（2）人口分布。點(diǎn)狀空間數(shù)據(jù)，新建中學(xué)應(yīng)在人口比較密集的地方，為此利用點(diǎn)狀人口調(diào)查數(shù)據(jù)產(chǎn)生密度圖。

（3）規(guī)劃土地使用。矢量多邊形空間數(shù)據(jù)，某些用地不應(yīng)當(dāng)建中學(xué)，如工業(yè)用地，某些用地不太適合建中學(xué)，如商業(yè)用地，而居住用地適合建中學(xué)。

根據(jù)上述數(shù)據(jù)和基本依據(jù)，將離開現(xiàn)有中學(xué)的距離分為4類，分別賦予評定指標(biāo)：0～500米：0，500～1000米：1，1000～1500米：2，>1500米：3圖層①人口密度分為4類（人/公頃），分別賦予評定指標(biāo)：

0～50：0，50～100：1，100～200：2，>200：3圖層②土地使用再分為3類，分別賦予評定指標(biāo)：

工業(yè)或綠地：0，商業(yè)：1，居?。?圖層③

新建中學(xué)選址的綜合評定指標(biāo)：圖層①×圖層②×圖層③

2產(chǎn)生離開現(xiàn)有中學(xué)的距離圖、再分類

啟動地圖文檔\\gis_ex09\\ex09\\ex09.mxd，激活dataframe1，有4個圖層（見圖9-1）：點(diǎn)狀圖層“現(xiàn)有中學(xué)〞、點(diǎn)狀圖層“人口調(diào)查〞、面狀圖層“規(guī)劃土地使用〞、線狀圖層“道路〞（只起限定柵格邊界和背景圖的作用）。鼠標(biāo)雙擊dataframe1，在彈出的DataFrameProperties對話框的General標(biāo)簽中將MapUnites和DisplayUnits從Unknown（未定義）改為Meters（米），完成后按“確定〞鍵關(guān)閉。選用菜單Tools/Extension…，加載SpatialAnalyst，在菜單View/Toolbars下勾選SpatialAnalyst，調(diào)出柵格分析工具條。

圖9-1dataframe1的顯示

設(shè)置SpatialAnalyst的初始參數(shù)，選用菜單SpatialAnalyst/Options…：

（1）General標(biāo)簽

Working：D:\\gis_ex09\\ex09\\temp\\鼠標(biāo)展開選擇SpatialAnalyst的工作路徑Analysismask：不選，本練習(xí)暫不考慮AnalysisCoordinateSystem：

●Analysisoutputwillbesavedinthesamecoordinatesystemastheinput(orfirstrasterinputiftherearemultiple…點(diǎn)選上側(cè)，產(chǎn)生柵格的坐標(biāo)系和輸入數(shù)據(jù)一致

（2）Extents標(biāo)簽

Analysisextent：SameasLayer：“道路〞下拉選擇圖層，限定分析空間的范圍（3）Cellsize標(biāo)簽

Analysiscell：AsSpecifiedBelow下拉選擇

Cellsize：50鍵盤輸入柵格單元的大小NumberofRows：82邊界和柵格單元大小確定后，自動確定柵格行數(shù)NumberofColumns：136邊界和柵格單元大小確定后，自動確定柵格列數(shù)

按“確定〞鍵，完成初始化設(shè)置。選用菜單SpatialAnalyst/Distance/StraightLine…，出現(xiàn)StraightLine參數(shù)設(shè)置對話框：

Distanceto：現(xiàn)有中學(xué)下拉選擇圖層名，消防站為距離的參照點(diǎn)Maximumdistance：不限定最大的計(jì)算范圍，保持空白Outputcellsize：50默認(rèn)，使用初始化設(shè)置

Outputraster：school鍵盤輸入柵格數(shù)據(jù)名稱，存放路徑為初始化設(shè)置

按OK鍵，產(chǎn)生離開現(xiàn)有中學(xué)的距離分布圖，點(diǎn)擊圖層名school，選用菜單SpatialAnalysis/Reclassify…，系統(tǒng)出現(xiàn)Reclassify對話框，點(diǎn)擊Classify…按鈕，出現(xiàn)Classification對話框：

Method:EqualInterval按等間距方式分類Classes：4共分4類

按OK鍵，返回Reclassify對話框，柵格數(shù)據(jù)被分為4類，還要按分析之前確定的要求修改分類的間距：

OldValues（原來值）NewValue（修改值）0–5000500–100011000–150021500–22003

在Outputraster中規(guī)定文件名及輸出路徑，這里為R_school，按OK鍵返回，系統(tǒng)產(chǎn)生新的分類圖層R_school（見圖9-2）。

圖9-2離開現(xiàn)有中學(xué)的距離再分類圖

3產(chǎn)生人口分布密度圖、再分類

選用菜單SpatialAnalysis/Density…，出現(xiàn)Density參數(shù)設(shè)置對話框：

Inputdata：人口調(diào)查下拉選擇圖層名，得到人口調(diào)查樣本點(diǎn)PopulationField：POPU下拉選擇，該字段在“人口調(diào)查〞要素屬性表中

DensityType：Kernel采用核心式計(jì)算方法，離樣本點(diǎn)近的單元密度大，遠(yuǎn)的密度小，大于探尋半徑時，單元的密度值為零

SearchRadius：800計(jì)算密度時的搜尋半徑AreaUnits:Squaremapunits面積單位用地圖單位的平方Outputcellsize：50已經(jīng)在初始化時設(shè)定，默認(rèn)Outputraster：popu_den鍵盤輸入柵格數(shù)據(jù)名稱，路徑按初始化設(shè)置

按OK鍵，產(chǎn)生估計(jì)的人口密度分布圖popu_den（產(chǎn)生密度圖的原理已在前章介紹）。選用菜單SpatialAnalysis/Reclassify…，出現(xiàn)Reclassify對話框，在Inputraster下拉菜單中選擇popu_den，點(diǎn)擊Classify…按鈕，出現(xiàn)Classification對話框：

Method:EqualIntervalClasses：4

按OK鍵后，返回Reclassify對話框，柵格數(shù)據(jù)被自動分為4類，再按要求修改：

OldValues（原來值）NewValue（修改值）

0–0.00500.005–0.0110.01–0.0220.02–0.033

在Outputraster中指定產(chǎn)生柵格數(shù)據(jù)的名稱及路徑，這里名稱為R_popu，路徑按初始設(shè)定，按OK鍵返回，

系統(tǒng)產(chǎn)生新的分類圖層“R_popu〞（見圖9-3）。

本章開始，人口密度的面積按公頃（萬平方米）考慮，練習(xí)中地圖單位是米，計(jì)算密度時按地圖單位的平方即平方米計(jì)算，因此柵格單元的取值和當(dāng)時確定的指標(biāo)之間是1：10000的關(guān)系，再分類后和原來的定義一致。

圖9-3重新分類后的人口密度圖

4規(guī)劃土地使用矢量圖層轉(zhuǎn)換為柵格圖、再分類

鼠標(biāo)右鍵開啟“Attributeof土地使用〞，可以看到該要素屬性表有字段LANDUSE，為每個多邊形地塊的規(guī)劃土地使用性質(zhì)，按城市規(guī)劃專業(yè)習(xí)慣，土地使用分類編碼的意義為：C：商業(yè)，G：綠地，M：工業(yè)，R1：一類居住，R2：二類居住。選用菜單SpatialAnalysis/Convert/FeaturestoRaster…（矢量要素轉(zhuǎn)換為柵格），系統(tǒng)出現(xiàn)FeaturestoRaster參數(shù)設(shè)置對話框：

Inputfeatures：土地使用下拉選擇圖層名

Field：landuse下拉選擇字段名，該字段在“土地使用〞要素屬性表中，決定柵格單元取值Outputcellsize：50已經(jīng)在初始化時設(shè)定

Outputraster：ld_use鍵盤輸入柵格數(shù)據(jù)名稱，路徑按初始設(shè)置

按OK鍵后，產(chǎn)生柵格狀土地使用圖層ld_use（見圖9-4），開啟圖層屬性表“Attributeofld_use〞，該表的內(nèi)容如下：

ValueCountlanduse

柵格單元取值取該值的單元數(shù)累計(jì)所對應(yīng)的轉(zhuǎn)換前的屬性

11568M（工業(yè)用地）2854C（商業(yè)用地）34120R2（一類居住用地）41079R1（二類居住用地）5658G（綠地）

圖9-4土地使用矢量圖轉(zhuǎn)化成柵格圖

按評價的需要對landuse字段再分類。選用菜單SpatialAnalysis/Reclassify…，系統(tǒng)出現(xiàn)Reclassify對話框，在Inputraster下拉菜單中選擇ld_use，在Reclassfield中下拉選擇landuse，表示對landuse字段進(jìn)行重新賦值。注意，這里應(yīng)直接修改NewValues，具體操作如下：

OldValues（原來值）NewVa