當桌面gis遇上python的簡單介紹

• 1、有哪些 GIS+Python 的開發經驗值得分享
• 2、arcgis 怎麼用python做統計分析
• 3、arcgis中打開python後，如何編程？
• 4、關於arcgis的python腳本編程， shape文件出png圖片問題
• 5、arcgis python怎麼運行

GISer入門指南電子書第一季 分享版.pptx免費下載

鏈接:

提取碼: wttn 

ArcGIS入門級教程,詳細介紹了ArcGIS產品族及成員產品 1、GIS 的概念和需求 2、什麼是ArcGIS 3、Geodatabase中的GIS數據概念 4、桌面GIS產品:ArcView, ArcEditor和ArcInfo等

方法/步驟

1

打開arcmap軟體，點擊工具欄上的ArcToolbox工具箱，在工具箱裡面選擇要進行數據處理的工具（工具的選擇視數據處理的內容而定），本文選擇了「Define Projection」定義投影的工具進行說明。

雙擊此工具，在彈出的工具窗口點擊右下角的「Show Help」按鈕。

2

點擊「Tool Help」按鈕，此時就打開了軟體的幫助文檔並定位到了此工具（define projection）對應的頁面，往下拉動滾動條，定位到此工具的python腳本的內容：

Syntax 介紹了此工具腳本的語法

Code Sample 介紹了腳本使用的例子

拷貝腳本使用例子裡面的代碼。

3

在軟體的主界面點擊「Geoprocessing」菜單–點擊「Python」菜單項，即出現內嵌在arcigis中的python工具窗口。

python工具窗口的左邊是代碼編輯區域，右邊是信息提示區域。

4

把第2步拷貝的代碼語句粘貼到python工具窗口的左邊區域，修改一下工具腳本的參數（視具體工具和數據而定），連續點擊兩下回車鍵。

5

等待腳本程序的執行完成，完成之後在桌面右下角會彈出提示信息，同時pythong工具窗口的右邊區域也會有具體的工具腳本執行信息。

6

執行完成之後，在arcmap載入處理過的數據，右鍵查看屬性，可以看到通過此工具腳本（define projection）進行投影定義的數據已經具有了投影參數。

此功能和Python的編輯模式是一樣的：

1.

編寫一行，運行一行；

2.

最大的的好處是，可以很方便的調用ArcGIS的介面，處理GIS數據；

3.

當然也可以使用ArcGIS

Toolboxes：

1、點類型的轉柵格：

PointToRaster example 1 (Python 

window)

Converts point features to a raster dataset.

import arcpy

from arcpy import env

env.workspace = “c:/data”

arcpy.PointToRaster_conversion(“ca_ozone_pts.shp”, “ELEVATION”, 

                                “c:/output/ca_elev”, “MAXIMUM”, “”, 2000)

PointToRaster example 2 

(stand-alone script)

Converts point features to a raster dataset.

# Name: PointToRaster_Ex_02.py

# Description: Converts point features to a raster dataset.

# Requirements: ArcInfo

# Import system modules

import arcpy

from arcpy import env

# Set environment settings

env.workspace = “C:/data”

# Set local variables

inFeatures = “ca_ozone_pts.shp”

valField = “ELEVATION”

outRaster = “c:/output/ca_elev02”

assignmentType = “MAXIMUM”

priorityField = “”

cellSize = 2000

# Execute PointToRaster

arcpy.PointToRaster_conversion(inFeatures, valField, outRaster, 

                               assignmentType, priorityField, cellSize)

2、面類型的轉柵格：

PolygonToRaster example 1 

(Python window)

Converts polygon features to a raster dataset.

import arcpy

from arcpy import env

env.workspace = “c:/data”

arcpy.PolygonToRaster_conversion(“ca_counties.shp”, “NAME”, 

                                 “c:/output/ca_counties.img”, 

                                 “MAXIMUM_AREA”, “MALES”, 0.25)

PolygonToRaster example 2 

(stand-alone script)

Converts polygon features to a raster dataset.

# Name: PolygonToRaster_Ex_02.py

# Description: Converts polygon features to a raster dataset.

# Requirements: ArcInfo

# Import system modules

import arcpy

from arcpy import env

# Set environment settings

env.workspace = “C:/data”

# Set local variables

inFeatures = “ca_counties.shp”

valField = “NAME”

outRaster = “c:/output/ca_counties”

assignmentType = “MAXIMUM_AREA”

priorityField = “MALES”

cellSize = 0.5

# Execute PolygonToRaster

arcpy.PolygonToRaster_conversion(inFeatures, valField, outRaster, 

                                 assignmentType, priorityField, cellSize)

3、線類型轉柵格

PolylineToRaster example 1 

(Python window)

Converts polyline features to a raster dataset.

import arcpy

from arcpy import env

env.workspace = “c:/data”

arcpy.PolylineToRaster_conversion(“roads.shp”, “CLASS”, “c:/output/roads.img”,

                                   “MAXIMUM_COMBINED_LENGTH”, “LENGTH”, 30) 

PolylineToRaster example 2 

(stand-alone script)

Converts polyline features to a raster dataset.

# Name: PolylineToRaster_Ex_02.py

# Description: Converts polyline features to a raster dataset.

# Requirements: ArcInfo

# Import system modules

import arcpy

from arcpy import env

# Set environment settings

env.workspace = “C:/data”

# Set local variables

inFeatures = “roads.shp”

valField = “CLASS”

outRaster = “c:/output/roads.tif”

assignmentType = “MAXIMUM_COMBINED_LENGTH”

priorityField = “LENGTH”

cellSize = 30

# Execute PolylineToRaster

arcpy.PolylineToRaster_conversion(inFeatures, valField, outRaster, 

                                  assignmentType, priorityField, cellSize)

我只知道桌面端，在地理處理菜單欄里（geoprocessing)打開Python，將語句拷貝進去，回車就可以執行，左邊欄里會有狀態顯示.