目錄：

一、處理點要素

1.創建編輯單點

?2.在一個要素中創建多點

二、處理線要素

1.創建和編輯單條線

2.創建和編輯多線為一個要素

三、處理面要素

1.創建和編輯環行線和單個多邊形

2.創建多個多邊形為一個要素

---

一、處理點要素

1.創建編輯單點

? ? ? ? 對于創建單點的操作，我們需要OGR庫中的幾個函數。ogr.Geometry(ogr.__)表示用來創建一個幾何對象要素，括號內__參數如下：

? ? ? ? 對于創建單點對象要素的操作，我們這里對于輸入ogr.wkbPoint確認幾何類型為2D的點。在點創建完成后可以使用Add.Point()來創建一個新的單點，同時可以采用.GetX/Y()獲取他的X/Y值，或直接打印該點查看。此外可以用SetPoint（point,x,y,[z]）來編輯這個點，其中point代表對象中點的幾何索引FID(如：0、1、2....，如果只有一個點，則僅有唯一FID為0)，X/Y/Z分別代表他的三個坐標。

實例代碼如下：

# 導入模板 from osgeo import ogr import os from ospybook.vectorplotter import VectorPlotter # 創建點 new_point = ogr.Geometry(ogr.wkbPoint) new_point.AddPoint(100,20) # 查看創建點信息 x,y=new_point.GetX(),new_point.GetY() print(x,y,new_point) # 打印到圖像上 vp=VectorPlotter(False) vp.plot(new_point,'rs') vp.draw() # 編輯第一個點，FID索引為0 new_point.SetPoint(0,10,20) print(new_point)

?2.在一個要素中創建多點

? ? ? ? 前面我們使用過單點要素的創建，但有時我們有時讓一個變量（要素）包含多個點。對于該操作我們至少需要創建兩個幾何對象。至少需要一個點對象和一個多點對象來容納。

# 導入模板 from osgeo import ogr import os from ospybook.vectorplotter import VectorPlotter # 創建單點和多點要素 new_points=ogr.Geometry(ogr.wkbMultiPoint) new_point=ogr.Geometry(ogr.wkbPoint) # 多次將單點導入單點 new_point.AddPoint(10,20) new_points.AddGeometry(new_point) new_point.AddPoint(20,30) new_points.AddGeometry(new_point) new_point.AddPoint(50,10) new_points.AddGeometry(new_point) # 將多點打印在圖像上 vp=VectorPlotter(False) vp.plot(new_points,'bo') vp.draw()

二、處理線要素

1.創建和編輯單條線

? ? ? ? 線要素是通過直線將多個頂點或點相連。創建線的第一步就是創建一個空的幾何對象，然后添加頂點（按順序）。

# 導入模板 from osgeo import ogr import os from ospybook.vectorplotter import VectorPlotter # 第一部分創建單線 # 創建單條線 new_line=ogr.Geometry(ogr.wkbLineString) new_line.AddPoint(5,6) new_line.AddPoint(7,9) new_line.AddPoint(10,14) new_line.AddPoint(1,5) # 打印在圖像中 vp=VectorPlotter(False) vp.plot(new_line,'b-') vp.draw() print(new_line) # 第二部分修改點 # 更改第四個頂點位置并打印 new_line.SetPoint(3,12,16) vp.plot(new_line,'y--') vp.draw() print(new_line.GetPoints()) # 第三部分在兩個點之間增加一個點 # 添加一個點_將頂點轉為元組 line_list=new_line.GetPoints() line_list[2:2]=[(30,40)] print(line_list) # 添加一個點_元組導成頂點 add_line=ogr.Geometry(ogr.wkbLineString) for i in line_list:add_line.AddPoint(*i) vp.plot(add_line,'r:') vp.draw() print(add_line)

? ? ? ? 在代碼第一個部分中，首先我們利用ogr.Geometry(ogr.wkbLineString)創建了一個線要素，然后挨個添加頂點（在線要素中，不會頂替上一個點的信息，這一點與點不同），增加的頂點會自己連成線，然后輸出打印。

? ? ? ? 在第二部分修改點坐標中，與點要素類似，通過SetPoint(索引FID,X,Y)來修改點的信息。

? ? ? ? 第三部分則實現了在兩點之間增加新點的操作。首先通過 對象.GetPoints()來獲取一個新的元組列表，用list[a:a]=[坐標元組]，將新的坐標值以元組形式傳入列表中。然后創建一個新的線要素，利用for循環將列表的元組轉化為坐標加入到新線要素中（注意：for循環時使用*將頂點炸開為兩個參數，而一個元組整體是不能應用于AddPoint（）函數的），并打印出圖。

2.創建和編輯多線為一個要素

? ? ? ? 對于多線的操作與點類似，分別建立不同的單線要素，再將這些要素添加到多線要素中。

# 導入模板 from osgeo import ogr import os from ospybook.vectorplotter import VectorPlotter import ospybook as pb # 創建單條線3組 new_line1=ogr.Geometry(ogr.wkbLineString) new_line1.AddPoint(5,6) new_line1.AddPoint(7,9) new_line1.AddPoint(10,14) new_line1.AddPoint(1,5) new_line2=ogr.Geometry(ogr.wkbLineString) new_line2.AddPoint(8,6) new_line2.AddPoint(7,4) new_line2.AddPoint(20,14) new_line3=ogr.Geometry(ogr.wkbLineString) new_line3.AddPoint(15,16) new_line3.AddPoint(4,9) new_line3.AddPoint(10,14) new_line3.AddPoint(11,5) # 創立多線要素，并將單線添加到里面 line=ogr.Geometry(ogr.wkbMultiLineString) line.AddGeometry(new_line1) line.AddGeometry(new_line2) line.AddGeometry(new_line3) # 結果顯示 vp=VectorPlotter(False) vp.plot(line,'b-') vp.draw() print(line)

三、處理面要素

1.創建和編輯環行線和單個多邊形

? ? ? ? 簡單來看，多邊形是由線構成，我們把這些線稱為環形線，因此當我們想構建一個多邊形就要先構建圍城多邊形的環形線。

? ? ? ? 開始構建環形線的操作與線相同，逐個添加頂點并形成線。這里注意的時，環形線需要閉合，在構建環形線時確定第一個點和最后一個點的坐標相同或環形或多邊形上調用CloseRings()函數使其閉合。??

? ? ? ? 環形線構成后添加到新建的面要素中即可完成對單面的創建。

# 導入模板 from osgeo import ogr import os from ospybook.vectorplotter import VectorPlotter import ospybook as pb # 創建環形線 new_ring=ogr.Geometry(ogr.wkbLinearRing) new_ring.AddPoint(10,14) new_ring.AddPoint(15,14) new_ring.AddPoint(15,22) new_ring.AddPoint(10,22) # 創建單個多邊形 new_polygon=ogr.Geometry(ogr.wkbPolygon) new_polygon.AddGeometry(new_ring) # 閉合多邊形，如果在環形線添加最后一個點也第一個點坐標相同可不進行此操作 new_polygon.CloseRings() # 繪制圖形 vp=VectorPlotter(False) vp.plot(new_polygon,'b') vp.draw()

? ? ? ? 對于多邊形的編輯我們首先要改變他的環形線，再對多邊形進行編輯。

# 導入模板 from osgeo import ogr import os from ospybook.vectorplotter import VectorPlotter import ospybook as pb # 此操作同創建單個多邊形 # 創建環形線 new_ring=ogr.Geometry(ogr.wkbLinearRing) new_ring.AddPoint(10,14) new_ring.AddPoint(15,14) new_ring.AddPoint(15,22) new_ring.AddPoint(10,22) # 創建單個多邊形 new_polygon=ogr.Geometry(ogr.wkbPolygon) new_polygon.AddGeometry(new_ring) # 閉合多邊形，如果在環形線添加最后一個點也第一個點坐標相同可不進行此操作 new_polygon.CloseRings() # 繪制圖形 vp=VectorPlotter(False) vp.plot(new_polygon,'b') vp.draw() # 編輯多邊形 # 改變某個點位置改變多邊形 ring=new_polygon.GetGeometryRef(0) ring.SetPoint(0,9,14) # 通過列表轉換，用兩個頂點代替第三個頂點 list=ring.GetPoints() list[2:3]=((16,14),(16,23)) for i in range(len(list)):ring.SetPoint(i,*list[i]) polygon=ogr.Geometry(ogr.wkbPolygon) polygon.AddGeometry(ring) vp.plot(polygon,'b') vp.draw()

2.創建多個多邊形為一個要素

? ? ? ? 創建多個多邊形的操作與線相同，分別創建不同的多邊形進行加入。

# 導入模板 from osgeo import ogr import os from ospybook.vectorplotter import VectorPlotter import ospybook as pb # 創建兩個多邊形 ring1=ogr.Geometry(ogr.wkbLinearRing) ring1.AddPoint(10,10) ring1.AddPoint(10,15) ring1.AddPoint(15,15) polygon1=ogr.Geometry(ogr.wkbPolygon) polygon1.AddGeometry(ring1) ring2=ogr.Geometry(ogr.wkbLinearRing) ring2.AddPoint(10,12) ring2.AddPoint(15,5) ring2.AddPoint(10,4) polygon2=ogr.Geometry(ogr.wkbPolygon) polygon2.AddGeometry(ring2) # 創建一個要素多個多邊形 polygon=ogr.Geometry(ogr.wkbMultiPolygon) polygon.AddGeometry(polygon1) polygon.AddGeometry(polygon2) polygon.CloseRings() vp=VectorPlotter(False) vp.plot(polygon,fill=False,ec='b') vp.draw()

總結

以上是生活随笔為你收集整理的OGR几何要素处理（创建编辑点线面）——Python地理数据处理学习分享的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。