本教程的講義和源碼都是取自Utah State University的openGIS課程

本人只是做點(diǎn)翻譯，寫寫學(xué)習(xí)體會(huì)而已，版權(quán)屬于原作者。

歡迎轉(zhuǎn)載，不過別忘了上面這段話。

==================================================

為什么用open source？

優(yōu)點(diǎn)

免費(fèi)，適合個(gè)人和小公司

強(qiáng)大的開發(fā)工具，找bug更容易

跨平臺(tái)，windows和linux都能用

拉風(fēng)！

缺點(diǎn)

沒有內(nèi)嵌地理處理器

用的人少

Open source RS/GIS模塊

OGR矢量庫(kù)：簡(jiǎn)單的矢量數(shù)據(jù)讀寫，是GDAL的一部分

GDAL地理空間數(shù)據(jù)抽象庫(kù)：

a)???????? 讀寫柵格數(shù)據(jù)

b)???????? ArcGIS也是基于GDAL開發(fā)的

c)???????? C++庫(kù)，但是可以用python調(diào)用

相關(guān)模塊

Numeric：高速的數(shù)組處理，對(duì)柵格數(shù)據(jù)尤其重要

NumPy：下一代的Numeric

導(dǎo)入庫(kù)：

import ogr

或者：

from osgeo import ogr

萬(wàn)能的方法是：

try:

from osgeo import ogr

except:

import ogr

要讀取某種類型的數(shù)據(jù)，必須要先載入數(shù)據(jù)驅(qū)動(dòng)，也就是初始化一個(gè)對(duì)象，讓它“知道”某種數(shù)據(jù)結(jié)構(gòu)。

import ogr

driver = ogr.GetDriverByName(‘ESRI Shapefile’)

數(shù)據(jù)驅(qū)動(dòng)driver的open()方法返回一個(gè)數(shù)據(jù)源對(duì)象

open(, )

其中update為0是只讀，為1是可寫

例如：

from osgeo import ogr

driver = ogr.GetDriverByName('ESRI Shapefile')

filename = 'C:/Users/gongwei/Documents/My eBooks/python_and_sage/GDAL python/test/ospy_data1/sites.shp'

dataSource = driver.Open(filename,0)

if dataSource is None:

print 'could not open'

sys.exit(1)

print 'done!'

注意filename一定要寫絕對(duì)路徑！

因?yàn)橐欢ㄒ媒^對(duì)路徑，為了簡(jiǎn)化代碼，經(jīng)常會(huì)使用到os.chdir()

讀取數(shù)據(jù)層

layer = dataSource.GetLayer(0)

一般ESRI的shapefile都是填0的，如果不填的話默認(rèn)也是0.

再看看這個(gè)數(shù)據(jù)層里面有幾個(gè)點(diǎn)呢？

n = layer.GetFeatureCount()

print 'feature count:', n

讀出上下左右邊界

extent = layer.GetExtent()

print 'extent:', extent

print 'ul:', extent[0], extent[3]

print 'lr:', extent[1], extent[2]

讀取某一要素feature(總算切入正題了)，這里讀取的是一個(gè)點(diǎn)

feat = layer.GetFeature(41)

fid = feat.GetField('id')

print fid

feat = layer.GetFeature(0)

fid = feat.GetField('id') #should be a different id

print fid

另外還有按順序讀取feature，循環(huán)遍歷所有的feature

feat = layer.GetNextFeature() #讀取下一個(gè)

while feat:

feat = layer.GetNextFeature()

later.ResetReading() #復(fù)位

提取feature的幾何形狀

geom = feat.GetGeometryRef()

geom.GetX()

geom.GetY()

print geom.

釋放內(nèi)存

feature.Destroy()

關(guān)閉數(shù)據(jù)源，相當(dāng)于文件系統(tǒng)操作中的關(guān)閉文件

dataSource.Destroy()

讀完了再說怎么寫

創(chuàng)建新文件

driver.CreateDataSource()

但是這個(gè)文件不能已經(jīng)存在了，否則會(huì)出錯(cuò)

創(chuàng)建新的layer

dataSource.CreateLayer(,CreateLayer(, geom_type=, [srs])

舉個(gè)例子：

ds2 = driver.CreateDataSource('test.shp')

layer2 = ds2.CreateLayer('test', geom_type=ogr.wkbPoint)

要?jiǎng)h除一個(gè)shp文件

driver.DeleteDataSource('test.shp')

要添加一個(gè)新字段，只能在layer里面加，而且還不能有數(shù)據(jù)

添加的字段如果是字符串，還要設(shè)定寬度

fieldDefn = ogr.FieldDefn('id', ogr.OFTString)

fieldDefn.SetWidth(4)

layer.CreateField(fieldDefn)

添加一個(gè)新的feature，首先得完成上一步，把字段field都添加齊了

然后從layer中讀取相應(yīng)的feature類型，并創(chuàng)建feature

featureDefn = layer.GetLayerDefn()

feature = ogr.Feature(featureDefn)

設(shè)定幾何形狀

feature.SetGeometry(point)

設(shè)定某字段的數(shù)值

feature.SetField('id', 23)

將feature寫入layer

layer.CreateFeature(feature)

總結(jié)

以上是生活随笔為你收集整理的python的ogr模块_GDAL python教程（1）——用OGR读写矢量数据的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。