GIS和空間分析的基本方法[1]

1、#1第一部分 GIS 和空間分析的基本方法第一章 ArcGIS 入門(mén)：數(shù)據(jù)管理及基本的空間分析工具地理信息系統(tǒng)（ GIS ）是一種計(jì)算機(jī)信息系統(tǒng)，用于獲取、存儲(chǔ)、處理、查詢、分析及顯 示地理數(shù)據(jù)。在 GIS 眾多功能中，制圖一直是它的一項(xiàng)主要功能。本章的首要目的是演示如何 借助 GIS 進(jìn)行電腦制圖。主要技巧包括空間與非空間（屬性）數(shù)據(jù)的管理以及二者之間的聯(lián) 系。但是， GIS 遠(yuǎn)不只是一種制圖工具，隨著 GIS 軟件功能越來(lái)越強(qiáng)大，界面越來(lái)越友好，它 在空間分析中的應(yīng)用越來(lái)越廣。本章第二個(gè)目的是介紹 GIS 的一些基本空間分析工具。鑒于 ArcGIS 在教育、商業(yè)及政府機(jī)構(gòu)中的廣泛應(yīng)用，本書(shū)

2、采用它作為主要的軟件平臺(tái)來(lái)完成 GIS 任務(wù)。除非特別說(shuō)明，本書(shū)所有 GIS 操作都是基于 ArcGIS9.0 。各章結(jié)構(gòu)的安排類(lèi) 似：先從基本概念著手，介紹 GIS 方法的基本內(nèi)容；然后用案例來(lái)使者熟悉剛剛介紹的 GIS 方法。本章第 1.1 節(jié)簡(jiǎn)要介紹 ArcGIS 中的空間及屬性數(shù)據(jù)管理，第 1.2 節(jié)為案例 1A ，通過(guò)繪 制庫(kù)婭霍加縣的人口密度分布圖來(lái)演示基于 GIS 作圖的基本過(guò)程。第 1.3 節(jié)介紹 ArcGIS 中基 本的空間分析工具，包括空間查詢、空間連接、地圖疊加等。第 1.4 節(jié)為案例 1B ，演示了一 些空間分析工具：提取克里夫蘭市的普查小區(qū)，生成多邊形鄰接矩陣。在高級(jí)

3、空間統(tǒng)計(jì)研究如 空間聚類(lèi)和空間回歸分析中，我們常常用多邊形鄰接矩陣來(lái)定義空間權(quán)重（參見(jiàn)第九章）。本書(shū)假設(shè)讀者具有初極 GIS 知識(shí)。本章不會(huì)涵蓋所有的 ArcGIS 功能。相反，它只是回顧ArcGIS 的主要功能，給讀者一個(gè) ”熱身賽 ”, 以便引導(dǎo)讀者繼續(xù)學(xué)習(xí)章節(jié)中的一些高級(jí)空間分析 方法。1.1ArcGIS 中的空間和屬性數(shù)據(jù)管理因?yàn)楸緯?shū)選擇 ArcGIS 作為主要的軟件平臺(tái)，所以有必要簡(jiǎn)單介紹一下它的主要模塊及功能。 ArcGIS 是美國(guó)環(huán)境系統(tǒng)研究所（ ESRI）2001 年發(fā)布的一種基于圖形界面（ GUI ）的軟件2平臺(tái)，用以代替以前基于命令行的 ArcInfo 。 ArcGIS 包

4、括 3 個(gè)主要模塊： ArcCatalog 、ArcMap 和 ArcToolbox 。 ArcCatalog 用于查看、管理空間數(shù)據(jù)文件。 ArcMap 用于顯示、分析、編輯空5間及屬性數(shù)據(jù)。 ArcToolbox 整合了各種數(shù)據(jù)管理和數(shù)據(jù)分析工具箱，包括地圖投影管理、數(shù)據(jù) 格式轉(zhuǎn)換以及實(shí)現(xiàn)早期的 ArcInfo 命令。在 ArcGIS9.0 中， ArcToolbox 可以從 ArcMap 或 ArcCatalog 界面中激活。大部分（但不是全部）早期 ArcInfo 命令都可以在 ArcGIS 中實(shí)現(xiàn)。對(duì) 于有些命令或功能，我們作空間分析還是用的較勤,仍然需要用到 ArcInfo 的命令

5、行界面。例如，在附錄 1 中討論了如何用 ArcInfo Workstation 來(lái)讀取、輸出 ASCII 文件。第二章中的案例 2 介紹了如何用 ArcInfo Workstation 來(lái)計(jì)算網(wǎng)絡(luò)距離。1.1.1 地圖投影及空間數(shù)據(jù)模型GIS 不同于其他信息系統(tǒng)的一個(gè)顯著特點(diǎn)是它能夠管理地理或空間（區(qū)位）數(shù)據(jù)。為了理 解這一點(diǎn)，我們需要具備地理坐標(biāo)系統(tǒng)的基礎(chǔ)知識(shí)，例如，我們需要了解經(jīng)緯度及用x、 y 坐標(biāo)在各種平面坐標(biāo)系上的表示方法。將地球的橢球表面轉(zhuǎn)換為平面，或在不同平面坐標(biāo)系之間 進(jìn)行轉(zhuǎn)換的過(guò)程稱(chēng)為地圖投影。在 ArcGIS 中， ArcMap 會(huì)自動(dòng)地將不同坐標(biāo)系的數(shù)據(jù)轉(zhuǎn)化為 最先添加

6、的圖層坐標(biāo)系，這個(gè)過(guò)程一般稱(chēng)為即時(shí)投影變換（ on-the-fly reprojections ）。但是， 如果數(shù)據(jù)量較大的話，這個(gè)過(guò)程可能會(huì)花很多時(shí)間。所以，同一個(gè)項(xiàng)目里面的所有圖層最好用 同一種投影。美國(guó)常用的投影有兩種：通用橫軸墨卡托投影（ UTM ）和州平面坐標(biāo)系統(tǒng) （State Plane Coordinate System， SPCS）。 嚴(yán)格來(lái)說(shuō)， SPCS并不是單一的一種投影，它可能 使用三種投影中的一種：蘭伯特等角圓錐投影、橫軸墨卡托投影、斜軸墨卡托投影。為了盡量 減少變形，南北向長(zhǎng)條形的州或地區(qū)使用橫軸墨卡托投影，東西向長(zhǎng)條形的州或地區(qū)使用蘭伯 特等角圓錐投影。有些州（如阿

7、拉斯加、紐約州）可能用不只一種投影。更多內(nèi)容，讀者可以 參閱 ESRI 公司 ArcGIS 光盤(pán)上的 “理解地圖投影 ”P(pán)DF文件。在 ArcGIS 中查看空間數(shù)據(jù)圖層的投影，可以在 ArcCatalog 中點(diǎn)擊該圖層，然后選擇 Metadata > Spatial ，或者在 ArcMap 右鍵單擊圖層，依次選擇 layer > Properties > Source 。在 ArcToolbox 中進(jìn)行投影相關(guān)操作的辦法為：依次點(diǎn)擊 Data Management Tools > Projections and Transformations ，在投影與變換工具箱（ P

8、rojections and Transformations ）下，投影定義工具（ Define Projection ）將新建一個(gè)包含投影參數(shù)的投影文件（PRJ），或者矯正原來(lái) 的錯(cuò)誤投影。投影定義工具只是標(biāo)記地理數(shù)據(jù)的正確坐標(biāo)系統(tǒng)，坐標(biāo)系統(tǒng)本身不會(huì)被改變。對(duì) 于矢量空間數(shù)據(jù)，選擇 Feature > Project ，將坐標(biāo)系統(tǒng)從一種投影轉(zhuǎn)換到另一種投影，并生成 一個(gè)新的圖層。本工具中提供了如下一些選項(xiàng)：創(chuàng)建新坐標(biāo)系統(tǒng)，使用一個(gè)現(xiàn)存的坐標(biāo)系統(tǒng)， 從一個(gè)地理數(shù)據(jù)中輸入坐標(biāo)系統(tǒng)。對(duì)于柵格數(shù)據(jù)，則選擇 Raster > Project Raster.GIS 可以處理矢量和柵格兩種空間

9、數(shù)據(jù)。在處理矢量數(shù)據(jù)時(shí)， GIS 用地理坐標(biāo)點(diǎn)來(lái)構(gòu)建 點(diǎn)、線、面等要素；在處理柵格數(shù)據(jù)時(shí)， GIS 借助按行列排列柵格單元來(lái)表征空間要素。柵格 數(shù)據(jù)的結(jié)構(gòu)比較簡(jiǎn)單，模型構(gòu)造相對(duì)容易。社會(huì)經(jīng)濟(jì)應(yīng)用中大多使用矢量數(shù)據(jù)，本書(shū)的大部分 GIS 應(yīng)用也使用矢量數(shù)據(jù)。大部分商業(yè) GIS 軟件都提供了矢量、柵格數(shù)據(jù)的互相轉(zhuǎn)換。在 ArcGIS 中，可以通過(guò)調(diào)用 ArcToolbox 中的轉(zhuǎn)換工具（ Conversion Tools ）來(lái)實(shí)現(xiàn)。ESRI 公司的早期 GIS 軟件使用 coverage 數(shù)據(jù)模型。后來(lái)在開(kāi)發(fā) ArcView 軟件包時(shí)，采用 了 shapefiles 數(shù)據(jù)模型。在 ArcGIS 8

10、 及以后的版本中，開(kāi)始使用 geodatabase 模型，代表了面 向?qū)ο髷?shù)據(jù)模型的新趨勢(shì)。面向?qū)ο蟮臄?shù)據(jù)模型把物體的幾何特征（空間數(shù)據(jù)）也當(dāng)作一種屬 性數(shù)據(jù)來(lái)存儲(chǔ)，而傳統(tǒng)的 coverage 和 shapefiles 模型是將空間和屬性兩種數(shù)據(jù)單獨(dú)存儲(chǔ)的。一 般而言，社會(huì)經(jīng)濟(jì)分析中的空間數(shù)據(jù)和屬性數(shù)據(jù)常常來(lái)自不同的渠道，在用 GIS 進(jìn)行制圖或分 析之前的一個(gè)基本任務(wù)是把他們連接在一起。這就涉及到下面將要介紹的屬性數(shù)據(jù)管理了。1.1.2 屬性數(shù)據(jù)管理及屬性連接GIS 數(shù)據(jù)分為兩類(lèi)，即空間數(shù)據(jù)和屬性數(shù)據(jù)。空間數(shù)據(jù)表征地圖要素的幾何特征，屬性數(shù) 據(jù)則描述要素的一些特征。屬性數(shù)據(jù)往往以表格或表格文

11、件的方式存儲(chǔ)。 Shapefile 屬性表使用 dBase 數(shù)據(jù)格式， ArcInfo Workstation 使用 INFO 格式， geodatabase表格使用 Microsoft Access 格式。 ArcGIS 也能夠讀取幾種 ASCII 數(shù)據(jù)格式，包括逗號(hào)分隔和 tab 分隔的文本文件。附錄 1 討介紹了用 ArcGIS 輸入 /輸出 ASCII 數(shù)據(jù)的方法，如果要用 GIS 和其它軟件（如 SAS）進(jìn)行 高級(jí)分析，或者自己編寫(xiě)一些程序來(lái)完成復(fù)雜的任務(wù)，這種數(shù)據(jù)轉(zhuǎn)換是非常重要的操作?；镜臄?shù)據(jù)管理任務(wù)，有些用 ArcCatalog 或 ArcMap 都可以實(shí)現(xiàn)，而其他一些任務(wù)則只

12、有其中一種工具可以完成。例如，新建數(shù)據(jù)表或刪除/拷貝一個(gè)數(shù)據(jù)表都只能在 ArcCatalog 中實(shí)9現(xiàn)（前面我們?cè)岬剑?ArcCatalog 可以用于查看和管理 GIS 數(shù)據(jù)文件）。數(shù)據(jù)表的創(chuàng)建過(guò)程 為：右鍵單擊將要在里面創(chuàng)建數(shù)據(jù)表的文件夾，選擇 new 即可。而要?jiǎng)h除或拷貝數(shù)據(jù)表，只 要在 ArcCatalog 中右鍵單擊該數(shù)據(jù)表，然后單擊 Delete（或 Copy ）即可。如果要在一個(gè)數(shù)據(jù)表中新增一個(gè)變量（可能是在 shape文件屬性表或 dBase 文件中新增一 列或是在 ArcInfo workstation 的 INFO 文件中新增一項(xiàng)）， ArcCatalog 和 ArcMap

13、 都可以勝 任。刪除 INFO 文件中的某一項(xiàng)也可以任選 ArcCatalog 或 ArcMap 之一來(lái)完成；但是，刪除 dBase 文件中的一列則只能用 ArcMap 來(lái)實(shí)現(xiàn)。例如，在 shape 文件屬性表中新增一列，可以 用 ArcCatalog 來(lái)實(shí)現(xiàn)，即右鍵點(diǎn)擊 shape文件 > Properties > Fields ，在空白行里面輸入新建列4 的名字，并定義數(shù)據(jù)類(lèi)型即可。在 ArcMap 中，則要先打開(kāi)數(shù)據(jù)表 > Options > Add Field 。在 ArcMap 中刪除一列，可以打開(kāi)數(shù)據(jù)表，然后右鍵單擊該列，然后選 Delete Field 即

14、可。如果要 進(jìn)行列之間數(shù)據(jù)的運(yùn)算，可以用 ArcMap ：打開(kāi)數(shù)據(jù)表 > 右鍵單擊列 > Calculate Values 。此 外，在 ArcMAP 中，可以通過(guò)右鍵單擊列，然后選擇 Statistics 來(lái)得到一些基本的統(tǒng)計(jì)參數(shù)。在 GIS 中，我們常常使用的屬性連接，是基于某一個(gè)相同列將兩張表的信息連接在一起。 被連接的表可能是一個(gè)與特定地理數(shù)據(jù)庫(kù)有關(guān)的屬性表，也可能是一個(gè)獨(dú)立的數(shù)據(jù)表。在進(jìn)行 屬性連接時(shí)，公共標(biāo)簽的名字可以不同，但它們的數(shù)據(jù)類(lèi)型必須匹配。數(shù)據(jù)表之間的連接關(guān)系 有多種：一對(duì)一、多對(duì)一、一對(duì)多、多對(duì)多。一對(duì)一或多對(duì)一連接是通過(guò) ArcGIS 中的 join 來(lái)

15、完成的。但是，一對(duì)多或多對(duì)多連接則不能用 join 來(lái)實(shí)現(xiàn)，這需要用 ArcGIS 中的 relate 來(lái)關(guān) 聯(lián)兩張表，在關(guān)聯(lián)的同時(shí)會(huì)保持兩張表各自獨(dú)立。在進(jìn)行關(guān)聯(lián)（relate ）時(shí)，是從一張表中選取一條或多條記錄，然后從另一張表中選取若干條關(guān)聯(lián)記錄。表1.1 列出了上述連接的關(guān)系及所用的 ArcGIS 工具。連接（ join ）或關(guān)聯(lián)（ relate）是通過(guò) ArcMap 來(lái)實(shí)現(xiàn)的。在地圖目錄中，右鍵單擊空間數(shù) 據(jù)或目標(biāo)表格，然后依次選擇 Joins and Relates > Join （或 Relate），然后，在連接數(shù)據(jù)對(duì)話框 中選擇 “ Joina ttributes fro

16、m a table ?！边B接只是暫時(shí)的，并不會(huì)新建數(shù)據(jù)表，如果退出項(xiàng)目時(shí)不 保存的話，下次再打開(kāi)時(shí)連接就沒(méi)有了?？梢园堰B接后的結(jié)果輸出為新的數(shù)據(jù)表從而永久地保 存。一旦屬性信息連接到空間圖層，我們就可以用 ArcGIS 方便地制圖了。在 ArcMAP 中，右 鍵單擊圖層，選擇 Properties，在彈出的對(duì)話框中選擇Symbology 。在這里，我們可以選取某一列指標(biāo)來(lái)繪制地圖，可以選擇要素顯示的顏色和圖例，設(shè)置顯示的模板。地圖要素（比例 尺、指北針、圖例）可以從主菜單中插入（Insert） i。1.2 案例 1A：繪制俄亥俄州庫(kù)婭霍加縣人口密度模式圖對(duì)于那些不太熟悉 GIS 的讀者，如果能

17、夠通過(guò)簡(jiǎn)單地點(diǎn)擊幾下鼠標(biāo)就可以親自繪制一張地5 圖，則可以很快嘗 GIS 的甜頭 , 克服對(duì) GIS 復(fù)雜性的神秘感。本節(jié)通過(guò)一個(gè)例子來(lái)演示如何在 GIS 連接空間和非空間信息，并用于繪制地圖。在接下來(lái)的操作中，我們將演示第 1.1 節(jié)里面 介紹的大部分功能。完成一個(gè) GIS 項(xiàng)目開(kāi)始之前先要收集相關(guān)數(shù)據(jù)。一般而言，我們可以使用現(xiàn)存的數(shù)據(jù)。在 美國(guó)社會(huì)經(jīng)濟(jì)應(yīng)用研究中，美國(guó)人口普查局發(fā)布的拓?fù)浼傻乩砭幋a參照文件（ TIGER ）及十 年一度的人口普查數(shù)據(jù)是空間數(shù)據(jù)和屬性數(shù)據(jù)的主要來(lái)源。這兩種數(shù)據(jù)都可以從普查局的網(wǎng)站 下載（ ）。熟練的 ArcGIS 用戶可以直接下載

18、 TIGER 數(shù)據(jù)，然后用 TIGER 轉(zhuǎn)換工具提取所需空間數(shù)據(jù)。依次點(diǎn)擊 ArcToolbox 里面的 Coverage Tools > Conversion > To Coverage > Advanced Tiger Conversion （或 Basic Tiger Conversion ）即可激活 TIGER 轉(zhuǎn)換工 具。轉(zhuǎn)換過(guò)程可能會(huì)花些時(shí)間，而轉(zhuǎn)換之后的數(shù)據(jù)可能還需要進(jìn)一步的加工。所幸的是，一些 網(wǎng)站提供了業(yè)已處理好的 ArcGIS 格式（ shapefiles 或 coverage）的 TIGER 空間數(shù)據(jù)下載服 務(wù)。 ESRI 發(fā)布 ArcGIS 軟件時(shí)，附

19、送有包含這些數(shù)據(jù)的光碟。如果空間數(shù)據(jù)是交換格式 （ e00 ） 的 ， 則 可 以 很 方 便 地 進(jìn) 行 轉(zhuǎn) 換 ： 依 次 點(diǎn) 擊 ： ArcToolbox > Coverage Tools > Conversion > To Coverage > Import from Interchange File 。在下面的案例中，我們將直接使用從 ESRI 網(wǎng)站下載的 shapefile 空間數(shù)據(jù)。盡管讀者可以自己下載數(shù)據(jù)，但為方便起見(jiàn)，本書(shū)光盤(pán)中提供了所需數(shù)據(jù)：1. shape 文件： tgr39035trt002. dBase 文件： tgr39000sf1trt.d

20、bf 。 在本書(shū)中，所有計(jì)算機(jī)文件、變量名以及某些工程中將要用到的命令行都用 Courier New 字體。面是分步介紹的操作過(guò)程。1. 下載空間數(shù)據(jù)登 錄 ESRI 公 司 的 網(wǎng) 站 ， 進(jìn) 入 到 2000 年 人 口 普 查 TIGER/Line Data 頁(yè) 面 ： 選擇 Ohio 州 Cuyahoga 縣，下 載 2000 的普查小區(qū)數(shù)據(jù)，為壓縮后的 shape 文件。解壓后得到名為 tgr39035trt00 的 shapefile 圖層。在 tgr39035trt00 這個(gè)文件名中， tgr 表示它來(lái)源于 TIGER 文件， 39 是州 的 FIPS 編碼， 035 是縣的編碼

21、， trt00 表示 2000 的普查數(shù)據(jù)。 Shapefile 圖層至少包含 3 個(gè) 文件： .dbf 、 .shp、 .shx。有些還包括其他一些文件如 .prj 、 .sbx、 .avl 、 .xml 。為方便起見(jiàn)，本書(shū) 余下部分一律用單數(shù)形式的 “ shapefile指”代一個(gè) shapefile 圖層所有文件。2.轉(zhuǎn)換到 UTM 投影在 ArcCatalog 中，查看 shapefile 文件 tgr39035trt00 的投影，發(fā)現(xiàn)它使用的是地理坐 標(biāo) 系 統(tǒng) 。 在 ArcToolbox 中 ， 依 次 選 擇 Data Management Tools > Project

22、ions and Transformations > Feature > Project ， 激 活 投 影 變 換 對(duì) 話 框 。 在 對(duì) 話 框 中 ， 選 擇 tgr39035trt00.shp 作為輸入數(shù)據(jù)，將輸出數(shù)據(jù)命名為 cuyautm.shp ，定義輸出坐標(biāo)系 統(tǒng)為 UTM(zone 17 ， units meters)。在這里，我們從一個(gè)現(xiàn)存的數(shù)據(jù)集中提取投影文件來(lái)定義輸 出的坐標(biāo)系統(tǒng)：點(diǎn)擊“ Output Coordinate System”旁的圖標(biāo)，激活空間參照屬性對(duì)話框>Import > clevbnd 。圖 1.1 為本任務(wù)的對(duì)話框。點(diǎn)擊 ok

23、執(zhí)行任務(wù)。3. 計(jì)算 shapefile文件中面要素的面積在ArcMap 中(本書(shū)中除非特別說(shuō)明， ArcMap 是我們 ArcGIS 的默認(rèn)工作環(huán)境)，打開(kāi) cuyautm 的屬性表，右鍵單擊 Options 按鈕，選擇 Add Field ，新增一列 area ，設(shè)置數(shù)據(jù)類(lèi)型為 雙精度 Double。右鍵單擊列 area ，選擇 Calculate Values 以計(jì)算面積。在對(duì)話框中，點(diǎn)擊11Advanced ，在第一個(gè)文本框中輸入下述 VBA 命令，Dim dblArea as doubleDim pArea as IAreaSet pArea = shapedblArea = pAr

24、ea.area在第二個(gè)文本框，即 “area= ” 下面的文本框中輸入 dblArea 。點(diǎn)擊 OK 計(jì)算面積。圖 1.2 為計(jì)算面積的對(duì)話框示例。需要說(shuō)明的是，在更新的 ArcGIS 版本（如 9.2）中 , 計(jì)算面積很簡(jiǎn)單：右鍵單擊列 area ， 選擇Calculate Geometry ，然后在對(duì)話框中的 Property項(xiàng)選Area 、Units項(xiàng)選 Square Meters sq m 以計(jì)算面積。4. 下載屬性數(shù)據(jù)7屬性數(shù)據(jù)可以從上面同一個(gè)網(wǎng)站下載，選擇“2000 census tract demographics （SF1） ”即可。屬性數(shù)據(jù) tgr39000sf1trt.db

25、f 是整個(gè)俄亥俄州的數(shù)據(jù)，為 dBase 格式， SF1 表示 Summary File 1 （即基于普查短表的 100%普查的匯總數(shù)據(jù)）。這里的 dBase 文件并不包含所有 的 普 查 變 量 。 如 要 獲 取 完 整 的 普 查 數(shù) 據(jù) ， 可 以 訪 問(wèn) 2000 年 人 口 普 查 網(wǎng) 站 ： /main/www/cen2000.html 。處理 SF1、 SF3 等文件時(shí)，需要知道 2000 年人 口普查數(shù)據(jù)的結(jié)構(gòu)，并要用到一些數(shù)據(jù)分析軟件如SAS、 Access。 SF3 就是基于普查長(zhǎng)表的1%抽樣數(shù)據(jù)。5. 提取凱霍加縣的屬性數(shù)據(jù)在 A

26、rcMap 中，添加并打開(kāi)表 tgr39000sf1trt.dbf 。點(diǎn)擊表右下角的選項(xiàng)按鈕 （Options ） > 選擇 Select By Attributes > 輸入 SQL （結(jié)構(gòu)查詢語(yǔ)言）語(yǔ)句： county= '035' 并按 Apply ，執(zhí)行查詢（第 1.3 節(jié)有更多關(guān)于 ArcGIS 的查詢操作）。所有凱霍加縣的數(shù)據(jù)將被 選中。點(diǎn)擊 Options 鍵，將結(jié)果輸出（ Export ）名為 cuya2k_popu.dbf 的文件。如果愿意，可以把該文件中多余的數(shù)據(jù)列刪去，只保留我們需要的STFID 和 POP2000 兩列。15這一步可以略過(guò)，即

27、直接將表 tgr39000sf1trt.dbf 連接到圖層 cuyautm ，所得結(jié)果 將自動(dòng)去掉其他縣的數(shù)據(jù)。這里設(shè)計(jì)這一步是為了讓讀者練習(xí)屬性查詢這一功能。6. 連接空間數(shù)據(jù)和屬性數(shù)據(jù)右 鍵 單 擊 圖 層 cuyautm > Joins and Relates > Join， 基 于 公 共 標(biāo) 簽 STFID 將 表 cuya2k_popu.dbf 連接到 cuyautm 。 STFID 是每個(gè)普查小區(qū)的唯一標(biāo)志碼，包含州（ 2 位）、縣（ 3 位）及普查小區(qū)（ 6 為）的編碼。圖 1.3 顯示了空間數(shù)據(jù)和屬性數(shù)據(jù)的連接方式 及與地圖要素的關(guān)系。7. 添加并計(jì)算人口密度右鍵

28、單擊圖層 cuyautm ，選擇 Open Attribute Table，打開(kāi)屬性表以檢查連接結(jié)果。在連 接后的表中，列的名稱(chēng)由兩部分組成，第一部分為原表名，第二部分為列名。 例如， cuyautm.area 表示該列為 cuyautm 屬性表中的 area 變量， tgr39000sf1trt.STFID 表示表 tgr39000sf1trt.dbf 中的 STFID （變量名過(guò)長(zhǎng)時(shí)會(huì)自動(dòng)截短）。在本書(shū)后續(xù)講解 （例如計(jì)算公式、表連接等操作）中，除非為了強(qiáng)調(diào)，當(dāng)我們提到列名稱(chēng)時(shí)，將省略源表名。單擊 Options 按鈕 > Add Field ，新增一列 popuden ，此列將顯示

29、在 cuyautm 屬性表的最 后，但位于表 cuya2k_popu.dbf 第一列之前。右鍵單擊列 popuden ，選擇 Calculate Values ，輸入公式 1000000*POP2000/area 。在公式中， POP2000 和 area 都是通 過(guò)直接雙擊最上面對(duì)話框中的列名稱(chēng)來(lái)實(shí)現(xiàn)的，這樣既可以節(jié)省時(shí)間，也可以減少直接輸入可 能出現(xiàn)的錯(cuò)誤。在本書(shū)后續(xù)講解中，為簡(jiǎn)單起見(jiàn)，我們只列出計(jì)算公式如 popuden=1000000*POP2000/area 。注意到地圖投影的單位是米，而人口密度的單位為每 平方公里的人口數(shù) ,所以我們公式中乘上了 1000,000。8. 繪制人口密

30、度模式圖右鍵單擊圖層 cuyautm > Properties > Symbology > Quantities > Graduated Colors ，基于 popuden 繪制人口密度圖。我們可以嘗試不同的分類(lèi)方法、分類(lèi)數(shù)及色彩方案。在主菜單 中，選擇 View > Layout View 來(lái)預(yù)覽地圖。在主菜單中，選擇 Insert > Legend（ Scale Bar, North Arrow 等），可以插入比例尺、指北針、圖例等地圖要素。圖 1.4 為研究區(qū)內(nèi)的人口密度圖，北部的那一大片空白區(qū)域?yàn)橐晾?。本圖中的人口密度 分級(jí)是作圖者自己定義 （C

31、ustomized） 。1.3ArcGIS 中的空間分析工具：查詢、空間連接、地圖疊加 許多空間分析任務(wù)需要利用空間要素之間的位置關(guān)系，進(jìn)行查詢（queries）、空間連接（ spatial joins ）及地圖疊加（ map overlays）。這三種都是空間分析的基本工具。 查詢包括屬性（非空間）查詢和空間查詢。屬性查詢是基于屬性表來(lái)提取在該屬性表中的 屬性信息，或?qū)?yīng)的空間圖層中的空間信息。在 ArcMap 中進(jìn)行屬性查詢，有兩種辦法（ 1） 從主菜單中選擇 Selection by Attributes ；（ 2）在一個(gè)打開(kāi)的數(shù)據(jù)表中，選擇Options 按鈕 >Selectio

32、n by Attributes 。兩種方法都是基于屬性變量用 SQL 查詢語(yǔ)句來(lái)選擇空間要素（或只是 簡(jiǎn)單地從一個(gè)獨(dú)立的屬性表中選取若干條記錄）。案例 1A 第 5 步中已經(jīng)用到了這個(gè)功能。在 主菜單的 Selection 菜單下，另有一個(gè)選項(xiàng)為交互式選擇方法（ Interactive Selection Method ）， 就是用鼠標(biāo)在屏幕（地圖或?qū)傩员恚┥线x擇要素。與其他信息系統(tǒng)相比， GIS 的一個(gè)獨(dú)特之處在于它的空間查詢能力，即能夠基于不同圖層 要素之間的位置關(guān)系進(jìn)行信息查詢。主菜單 Selection 菜單有一個(gè) Selection by Location 選項(xiàng)， 它可以基于一個(gè)圖層

33、中的要素與另一個(gè)圖層中的要素之間的位置關(guān)系進(jìn)行查詢?？晒┎樵兊目?間位置關(guān)系包括相交（ intersect）、在一定范圍之內(nèi)（ are within a distance of ）、完全包含 （ completely contain ）、完全在范圍之內(nèi)（ are completely within ），等等。查詢（屬性查詢或空間查詢）所得結(jié)果可以輸出為新的數(shù)據(jù)文件：（1）用右鍵單擊源圖層，然后選擇 Data > Export Data，可以將查詢所得空間要素輸出為新的圖層文件；（2）單擊數(shù)據(jù)表中的 Options 按鈕 > Export ，可以將結(jié)果保存為數(shù)據(jù)表。屬性連接是基于兩表

34、的共同列進(jìn)行的，而空間連接是基于兩個(gè)圖層空間要素的位置關(guān)系如 重疊、近鄰等來(lái)完成的。我們用源圖層和目標(biāo)圖層來(lái)區(qū)分不同圖層在空間連接時(shí)的作用相異：源圖層的屬性經(jīng)過(guò)空間連接后就轉(zhuǎn)換到目標(biāo)圖層中了。如果源圖層中的一個(gè)對(duì)象對(duì)應(yīng)于目標(biāo)圖 層中的一個(gè)或多個(gè)對(duì)象，則是一個(gè)簡(jiǎn)單的空間連接（ Simple join ）。例如，將縣域多邊形圖層 （源圖層）空間連接到學(xué)校位置（目標(biāo)圖層）的點(diǎn)圖層，縣圖層的屬性（如 FIPS 編碼、縣 名、縣長(zhǎng)名）就轉(zhuǎn)到那些落入縣域邊界內(nèi)的學(xué)校。如果源圖層的多個(gè)對(duì)象對(duì)應(yīng)于目標(biāo)圖層的一 個(gè)對(duì)象，就可能進(jìn)行兩種操作：匯總連接（summarized join ）和距離連接（ distanc

35、e join ）。匯總連接是指將源圖層的數(shù)值屬性進(jìn)行匯總（例如取平均值、求和、最小值、最大值、標(biāo)準(zhǔn)差、 方差）后再將結(jié)果添加到目標(biāo)圖層中。距離連接是從源圖層所有要素中尋找一個(gè)距離目標(biāo)圖層 中的匹配對(duì)象最近的要素，然后將它的屬性及距離值（為二要素之間的距離）添加到目標(biāo)圖層 中。例如，我們可以將地理位置編碼的犯罪數(shù)據(jù)（源圖層）與普查小區(qū)（目標(biāo)圖層）進(jìn)行空間 連接，從而得到按普查小區(qū)匯總的犯罪數(shù)，這就匯總連接；我們也可以將公交站點(diǎn)圖層（源圖 層）與普查街區(qū)重心的點(diǎn)圖層（目標(biāo)圖層）進(jìn)行空間連接，從而得到距離每個(gè)普查街區(qū)最近的 公交站點(diǎn)，這就是距離連接。不同空間要素之間的連接方式多種多樣（Price,

36、2004: 287-288 ）。表 1.2 是 ArcGIS 中各種空間連接的小結(jié)。與屬性連接類(lèi)似的是，空間連接是在 ArcMap 中實(shí)現(xiàn)：右鍵單擊源圖層 > 選擇 Joins and Relates > Join 。在連接數(shù)據(jù)對(duì)話窗口中，選擇 “ Join data from another layer based on spatial location 而不是”“ Join attributes from a table 。后者就是前”面已經(jīng)講述的屬性連接。地圖疊加可以廣義地定義為任何不同圖層改變要素的空間分析?？梢杂?ArcGIS 來(lái)實(shí)現(xiàn) 常用的地圖疊加工具主要有：剪切（ C

37、lip ）、相交（ Intersect ）、合并（ Union ）、緩沖區(qū) （ Buffer ）、多重緩沖區(qū)（ Multiple Ring Buffer ）。剪切是用一個(gè)圖層的邊界來(lái)截取另一個(gè)圖 層。相交是取疊加兩圖層的公共部分即交集。合并是取兩個(gè)圖層的并集。緩沖區(qū)是基于點(diǎn)、線 或面狀要素向外擴(kuò)展一定的緩沖距離形成面狀要素。多重緩沖區(qū)是基于多個(gè)距離同時(shí)生成一系 列的緩沖區(qū)。在 ArcGIS 9.0 中，上述地圖疊加工具分散在 ArcToolbox > Analysis Tools 中的不 同位置：剪切在 Extract 工具庫(kù)中，相交、合并在 Overlay 工具庫(kù)中，緩沖區(qū)、多重緩沖在

38、Erase，參見(jiàn)第 1.4.2 節(jié)第 3Proximity 工具庫(kù)中。本書(shū)用到的其他地圖疊加工具包括：刪除（ 23步）、近鄰（ Near ，參見(jiàn)第 2.3.2 節(jié)第 2 步）、點(diǎn)距離（ Point Distance，參見(jiàn)第 2.3.1 節(jié)第 2 步）、邊界合并（ Dissolve ，參見(jiàn)第 4.3.1 節(jié)第 2步）、附加（ Append ，參加第 4.3節(jié)）等。有讀者可能已經(jīng)注意到空間查詢、空間連接及地圖疊加之間的相似之處。事實(shí)上，許多空 間分析任務(wù)可以用這三種方法中的任何一種來(lái)實(shí)現(xiàn)。表 1.3 列出了他們之間的區(qū)別?？臻g查詢 只是尋找并顯示所需信息，它本身并不創(chuàng)建新的圖層（除非我們將選中的要

39、素輸出成新的文 件）。空間連接總是將連接結(jié)果保存為一個(gè)新的圖層。空間連接與地圖疊加有重要區(qū)別?？臻g 連接只是識(shí)別輸入圖層空間要素之間的位置關(guān)系，它并不改變?cè)瓉?lái)的空間要素，也不創(chuàng)建新的 要素。在地圖疊加過(guò)程中，一些輸入要素被分割、融合或刪去除以創(chuàng)建一個(gè)新的圖層。地圖疊 加比空間連接的運(yùn)算時(shí)間長(zhǎng)，而空間連接的運(yùn)算時(shí)間又比空間查詢長(zhǎng)。1.4 案例 1B：提取克里夫蘭市的普查小區(qū)并進(jìn)行多邊形鄰接分析本例將用到隨書(shū)光盤(pán)中的如下數(shù)據(jù)：1. shapefile 文件 cuyautm ：俄亥俄州庫(kù)亞霍加縣的所有普查小區(qū)；2. coverage 圖層文件 clevbnd ：俄亥俄州克里夫蘭市的邊界。光盤(pán)中的所有

40、 coverage 文件都是 ArcInfo 交換文件格式的（ .e00），需要先進(jìn)行格式轉(zhuǎn) 換：選擇 ArcToolbox 的 Coverage Tools > To Coverage > Import From Interchange File 。這里的 shapefile 文件 cuyautm 是案例 1A 所得結(jié)果，但為了方便在光盤(pán)中也提供了這個(gè)數(shù)據(jù)，以便 讀者可以獨(dú)立于案例 1A 進(jìn)行下面的操作。 Coverage 文件 clevbnd 是從網(wǎng)上下載的。 1.4.1 提取克里夫蘭市的普查小區(qū)1213很多情況下， GIS 使用者需要從一個(gè)較大區(qū)域中提取一個(gè)較小的研究區(qū)。下面

41、我們要做的 就是從庫(kù)亞霍加縣中提取克里夫蘭市的普查小區(qū)。將兩個(gè)圖層疊加后發(fā)現(xiàn)， clevbnd 的邊界 跟 cuyautm 邊緣處的普查小區(qū)邊界不完全重合。 Cuyautm 包含更多的地理細(xì)節(jié)。雖然二者 的邊界不完全一致，但 Cuyautm 中各普查小區(qū)的重心都落在 clevbnd 邊界之內(nèi)。我們的目 的就是要找到 cuyautm 圖層中那些落入 clevbnd 邊界之內(nèi)的重心，從而提取克里夫蘭普查些 cuyautm 中的地理信息。1. 生成庫(kù)亞霍加縣普查小區(qū)的重心激活 ArcToolbox > 選擇 Data Management Tools > Features > Fe

42、ature To Point 。在彈出的對(duì) 話框中，選擇 cuyautm 作為輸入要素，將輸出要素命名為 cuya_pt ，并選中 Inside 的選 項(xiàng)，于是得到普查小區(qū)重心的 shapefile 文件 cuya_pt 。2. 識(shí)別區(qū)市邊界內(nèi)的普查小區(qū)重心右鍵單擊目標(biāo)圖層 cuya_pt ，選擇 Joins and Relates > Join 。在對(duì)話框中，選擇 Join data from another layer based on spatial location ，設(shè)置源圖層( source layer )為 clevbnd polygon ， 連接選項(xiàng)為 it falls

43、inside ，將輸出結(jié)果命名為 tmp1 。圖 1.5 為空間連接對(duì)話框示意圖。打開(kāi) tmp1 的屬性表可以看到，對(duì)于城市邊界之內(nèi)的普查小區(qū)， clevbnd_id= 1，而邊界之外的普查小區(qū)， clevbnd_id= 0 。143. 將普查小區(qū)重心信息添加到多邊形圖層添加 cuyautm 圖層，右鍵單擊圖層并選擇 Joins and Relates > Join > 選擇 Join attributes from a table，用 tmp1 作為源數(shù)據(jù)表， STFID 為連接關(guān)鍵詞(目標(biāo)圖層 cuyautm 和源圖層 tmp1 都以此為關(guān)鍵詞)。4. 提取市區(qū)內(nèi)的普查小區(qū)打開(kāi)

44、 cuyautm 的屬性表 > 單擊 Options 標(biāo)簽 > 選擇 Select by Attributes > 設(shè)置選擇標(biāo) 準(zhǔn) tmp1.clevbnd_id = 1 。所有在城市之內(nèi)的多邊形將被選中并加亮。右鍵單擊圖層 cuyautm 并選擇 Data > Export Data ，需要注意的是，最頂端Export 一覽中為 Selectedfeatures，將輸出結(jié)果命名為 clevtrt 。所得 shapefile 文件 clevtrt 為克里夫蘭市的所有普 查小區(qū)。上面我們用到了空間連接。正如第 1.3 節(jié)中介紹的那樣，我們也可以用空間查詢Selection

45、 >( Selection by Location )，或地圖疊加工具( ArcToolbox > Analysis Tools > Overlay >Identity )來(lái)完成上述任務(wù)。例如，用空間查詢的方法為：在主菜單中，點(diǎn)擊Selection by Location > 使用查詢條件 select features from cuyautm that have their center in clevbnd polygon> 將所得結(jié)果輸出為 clevtrt 的 shapefile 文件。1.4.2 識(shí)別鄰接多邊形空間分析中，定義多邊形鄰接矩陣是一項(xiàng)非

46、常重要的任務(wù)。例如在第九章基于面單元的空 間聚類(lèi)和空間回歸時(shí)，我們用鄰接矩陣來(lái)定義空間權(quán)重,從面分析表征空間自相關(guān)。多邊形之間的鄰接包括兩種：（ 1）R 鄰接，也稱(chēng)邊鄰接，即兩個(gè)多邊形有一段共同的邊界；（2）Q 鄰接，也稱(chēng)廣義鄰接，指兩個(gè)多邊形有共同點(diǎn)或邊時(shí)即為鄰接多邊形（Cliff and Ord, 1973 ）。對(duì)于 R 鄰接，我們可以用 ArcInfo Workstation 里面的 PALINFO 命令來(lái)得到多邊形鄰接矩陣。下 面的案例分析使用 Q 鄰接來(lái)定義鄰接多邊形，依次來(lái)演示第1.3 節(jié)中介紹的空間分析工具。這里我們以一個(gè)普查小區(qū)為例，尋找完的鄰接多邊形。與第 1.4.1 節(jié)中第

47、 4 步類(lèi)似，從 clevtrt 圖層中選取 TRACTID = 1038' 的普查小區(qū)，并將其輸出成名為 zonei 的 shapefile 文件。下面將基于 Q 鄰接來(lái)尋找 zonei 的鄰接普查小區(qū)。圖 1.6 顯示了抽樣小區(qū)周 邊普查小區(qū)的 TRACTID 值?；?Q 鄰接，普查小區(qū) 1038 共有 6 個(gè)鄰接普查小區(qū)（ 1026、 1028、 1029、1035、1036 和 1039）。如果是基于 R 鄰接，則普查小區(qū) 1028 將不包括在鄰接 小區(qū)內(nèi)。下面是識(shí)別 zonei 的鄰接小區(qū)的具體過(guò)程。1. 生成緩沖普查小區(qū)在 ArcToolbox 中，選擇 Analysis

48、 Tools > Proximity > Buffer 。以一個(gè)較小的距離（如 30 米）在 zonei 周邊做一個(gè)緩沖區(qū)，將結(jié)果輸出名為 zonei_buff 的文件。緩沖區(qū)的距離必須 大于模糊容差（對(duì)于這里的 zonei，模糊容差大約為 1 米），但也要足夠小，以避免超出直接 相鄰多邊形的范圍。2. 從研究區(qū)中剪切緩沖區(qū)在 ArcToolbox 中，選擇 Analysis Tools > Extract > Clip 。選擇 clevtrt作為輸入要素， zonei_buff 作為剪切要素，將輸出結(jié)果命名為 zonei_clip3. 提取鄰接多邊形在 ArcTool

49、box 中，選擇 Analysis Tools > Overlay > Erase 。選取 zonei_clip 作為輸入要 素， zonei 作為刪除要素，將輸出結(jié)果命名為 zonejs 。所得的 shapefile 文件 zonejs 包含 了 zonei 基于 Q 鄰接的所有鄰接普查小區(qū)。圖 1.7 演示了操作過(guò)程。圖層 zonei 只包含一個(gè)普查小區(qū)即1038 。所得緩沖區(qū)zonei_buff 只包含一個(gè)多邊形（緩沖區(qū)里面的普查小區(qū) 1038 在圖 1.7 里示出來(lái)只是為了進(jìn) 行對(duì)比）。用緩沖區(qū)剪切研究區(qū)得到 zonei_clip ，包括普查小區(qū) 1038 一共七個(gè)多邊形。

50、 通過(guò)刪除操作可以從 zonei_clip 中去掉原來(lái)的普查小區(qū) zonei 。所得結(jié)果 zonejs 包含 6 個(gè)普查小區(qū)（圖 1.7 中給出了各自的 TRACTID 值）。生成研究區(qū)的多邊形鄰接矩陣需要對(duì)所有普查小區(qū)循環(huán)進(jìn)行上述操作。我們?cè)诠獗P(pán)中提供了一個(gè) AML （ Arc Micro Language）程序 Queen_Cont.aml 來(lái)實(shí)現(xiàn)這個(gè)任務(wù)，該程序是基于 （沈清 1994） 的工作開(kāi)發(fā)的。最后，我們可以用 ArcCatalog 刪除不需要的文件來(lái)節(jié)省存儲(chǔ)空間，我們主要保留 cuyautm 、 cuya_pt 和 clevtrt 三個(gè)文件即可，后續(xù)章節(jié)中將用到這三個(gè)文件。1.

51、5 小結(jié)16總結(jié)本章介紹的 GIS 及空間分析技巧主要包括如下幾點(diǎn)：1. ArcGIS 中的空間數(shù)據(jù)格式及格式轉(zhuǎn)換；2. 地圖投影及投影變換；3. 屬性數(shù)據(jù)管理（創(chuàng)建、編輯、刪除數(shù)據(jù)庫(kù)及數(shù)據(jù)列）；4. 屬性連接（包括屬性數(shù)據(jù)與空間數(shù)據(jù)的連接）；5. 繪制屬性地圖；6. 空間連接；7. 屬性與空間查詢；8. 地圖疊加操作（剪切、緩沖區(qū)、相交、合并、刪除）。2517其他重要的操作還包括：從公共資源中查找空間及屬性數(shù)據(jù)，計(jì)算 shapefile 文件中多邊形的面積，從多邊形圖層中提取重心來(lái)創(chuàng)建一個(gè)點(diǎn)圖層。后續(xù)章節(jié)將用到上述技巧。本章也介紹了一些基本概念如數(shù)據(jù)表之間的各種關(guān)系（一對(duì)一、多對(duì)一、一對(duì)多、

52、多對(duì)多），各種空間連接，空間查詢、空間連接與地圖疊加之間的區(qū)別。讀者如果想進(jìn)行更多的基于 GIS 的制圖練習(xí)，可以下載某個(gè)自己熟悉地區(qū)的普查數(shù)據(jù)和TIGER 數(shù)據(jù)，然后繪制該地區(qū)的一些人口統(tǒng)計(jì)屬性（人口、種族、年齡、性別等）和社會(huì)經(jīng)濟(jì)變量（收入、貧困狀況、教育水平、家庭結(jié)構(gòu)、住房屬性等） 。附錄 1 用 ArcGIS 輸入 /輸出 ASCII 文件小的 ASCII （文本）格式的數(shù)據(jù)文件很容易轉(zhuǎn)換成 ArcGIS 可以識(shí)別的格式。一般來(lái)說(shuō)， 我們可以用 Microsoft Excel 來(lái)打開(kāi)文件，在最頂端添加一行來(lái)作為列名，并將其保存為 CSV 文件（逗號(hào)分隔的文本文件）。 CSV 文件可以用

53、 ArcGIS 讀取。一個(gè)比較大的問(wèn)題是，按這種 方式處理時(shí)，沒(méi)法控制每列數(shù)據(jù)的類(lèi)型及格式。例如，在 ASCII 文件中，所有普查小區(qū)編碼或 STFID 值看起來(lái)都是數(shù)值型的（ numeric ）。從 ASCII 文件轉(zhuǎn)換過(guò)來(lái)的 CSV 文件自動(dòng)地定義各 列為數(shù)值型，從而使它很難與從 TIGER 文件提取的 GIS 圖層進(jìn)行連接，因?yàn)樵?TIGER 圖層 中，這些列通常定義為字符型（ character）。ASCII 文件。下面介紹如何用怎樣在 ArcGIS 中輸出 ASCII 文件呢 ? 在 ArcMap 中，我們可以打開(kāi)一個(gè)數(shù)據(jù)表（可能是 空間要素的屬性數(shù)據(jù)表，也可能是一張獨(dú)立的表），選擇

54、 Options > Export ，將其輸出成 dBase 格式的表格。然后，我們可以用 Microsoft Excel 打開(kāi) dBase 文件，再將其保存為 ASCII 文件。 另外一種方法是，使用 ArcToolbox > Spatial Statistics Tools > Utilities > Export Feature Attribute to Ascii ，將要素坐標(biāo)及所選屬性值輸出成空格、逗號(hào)或分號(hào)分隔的 ASCII 文本文件。但局限 的是，我們每次只能輸出一個(gè)變量（及其坐標(biāo)）。Microsoft Access 也常常用于屬性數(shù)據(jù)管理，即輸入、輸出2

55、9ArcInfo Workstationk 來(lái)輸入 / 輸出大型 ASCII 文本文件。1. 將 ASCII 文件轉(zhuǎn)換為 INFO 文件a. 在 ArcCatalog 中，打開(kāi)將用于保存新文件的目錄或工作空間（workspace），選擇File > New > INFO table 來(lái)創(chuàng)建一個(gè)新的 INFO 文件（比方說(shuō)可以命名為 ninfo ），定義所有的列（名稱(chēng)、數(shù)據(jù)類(lèi)型、數(shù)據(jù)格式）。b. 啟動(dòng) ArcInfo 命令行界面：點(diǎn)擊電腦的開(kāi)始菜單 > 程序 > ArcGIS > ArcInfoWorkstation > Arc 。 輸 入 w 以 切 換 到 相 應(yīng) 的 工 作 目 錄 （ 例 如 ： w COLUMNAR 選項(xiàng)， UNLOAD 命令輸出的結(jié)果為逗號(hào)分隔的 ASCII 文件，字符型數(shù) 據(jù)列將自動(dòng)在各數(shù)據(jù)上添加單引號(hào)（ ）'。c:Quant_GISproj1），輸入 tables以激活 TABLES 模塊。c. 在 TABLES 模塊下，輸入 select ninfo來(lái)選擇剛剛新建的 INFO 文件。d. 假設(shè) tfile 是需要轉(zhuǎn)換的文本文件。輸入 add from tfile來(lái)添加所有數(shù)據(jù)。轉(zhuǎn)換后的數(shù)據(jù)可以通過(guò) list 命令來(lái)查看。e. 輸入 quit 以退出 TABLES 模塊，再輸入 quit以退出 ArcI