# -*- coding: utf-8 -*-
print(u"python与开源QGis课题研究组")
#print("汉字")
#&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;
#创建矢量数据文件
#&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;&＃43;
try:
from osgeo import ogr
except:
import ogr
driver &＃61; ogr.GetDriverByName("ESRI Shapefile")
import os
&＃39;&＃39;&＃39;
ds &＃61; driver.CreateDataSource("xx_tesp.shp")
layer &＃61; ds.CreateLayer(&＃39;test&＃39;,geom_type &＃61; ogr.wkbPoint)
fieldDefn &＃61; ogr.FieldDefn(&＃39;id&＃39;,ogr.OFTString)
fieldDefn.SetWidth(4)
layer.CreateField(fieldDefn)
featureDefn &＃61; layer.GetLayerDefn()
print(featureDefn)
feature &＃61; ogr.Feature(featureDefn)
#设定几何形状
point &＃61; ogr.Geometry(ogr.wkbPoint)
point.SetPoint(0,123,123)
#设定字段数值
feature.SetField(&＃39;id&＃39;,23)
#将要素写入到图层
layer.CreateFeature(feature)
ds.Destroy()
import os
out_shp &＃61; "xx_tesp.shp"
dir(out_shp)
if os.path.exists(out_shp):
driver.DeleteDataSource(out_shp)
dir(out_shp)
point &＃61; ogr.Geometry(ogr.wkbPoint)
print(point)
point.AddPoint(10,20)
print(point)
point.AddPoint(30,20)
print(point)
line &＃61; ogr.Geometry(ogr.wkbLineString)
print(line)
line.AddPoint(10,10)
print(line)
line.AddPoint(20,20)
print(line)
line.SetPoint(2,30,30)
print(line)
print(line.GetPointCount())
print(line.GetX(0))
print(line.GetX(1))
print(line.GetX(3))
ring &＃61; ogr.Geometry(ogr.wkbLinearRing)
ring.AddPoint(0,0)
ring.AddPoint(100,0)
ring.AddPoint(100,100)
ring.AddPoint(0,100)
ring.CloseRings()
print(ring)
print(type(ring))
outring &＃61; ogr.Geometry(ogr.wkbLinearRing)
outring.AddPoint(0,0)
outring.AddPoint(100,0)
outring.AddPoint(100,100)
outring.AddPoint(0,100)
outring.AddPoint(0,0)
inring &＃61; ogr.Geometry(ogr.wkbLinearRing)
inring.AddPoint(25,25)
inring.AddPoint(75,25)
inring.AddPoint(75,75)
inring.AddPoint(25,75)
inring.CloseRings()
polygon &＃61; ogr.Geometry(ogr.wkbPolygon)
polygon.AddGeometry(outring)
polygon.AddGeometry(inring)
print(polygon.GetGeometryCount())
#outring2 &＃61; polygon.GetGeometryRef(0)
#inring2 &＃61; polygon.GetGeometryRef(1)
#print(outring2)
#print(inring2)
for ringx in range(polygon.GetGeometryCount()):
print(polygon.GetGeometryRef(ringx))
multipoint &＃61; ogr.Geometry(ogr.wkbMultiPoint)
point &＃61; ogr.Geometry(ogr.wkbPoint)
point.AddPoint(10,10)
multipoint.AddGeometry(point)
print(multipoint)
point.AddPoint(20,20)
multipoint.AddGeometry(point)
print(multipoint)
mp &＃61; multipoint.GetGeometryCount()
print(mp)
for mpx in range(multipoint.GetGeometryCount()):
print(multipoint.GetGeometryRef(mpx))
extfile &＃61; &＃39;xx_data_pt.shp&＃39;
if os.access(extfile,os.F_OK):
driver.DeleteDataSource(extfile)
#1、创建数据源
newds &＃61; driver.CreateDataSource(extfile)
#print(dir(newds))
#print(dir(newds.CreateLayer))
#2、创建数据源图层
lyrn &＃61; newds.CreateLayer(&＃39;point&＃39;,None,ogr.wkbPoint)
#3、定义图层字段&＃xff0c;添加图层字段
fieldcnstr &＃61; ogr.FieldDefn("idx",ogr.OFTInteger)
lyrn.CreateField(fieldcnstr)
fieldf &＃61; ogr.FieldDefn("namex",ogr.OFTString)
lyrn.CreateField(fieldf)
point_coors_arr &＃61; [[1,0],[2,0],[3,0],[4,0]]
for idxx,point_coors in enumerate(point_coors_arr):
#print(type(idxx))
#print(type(point_coors))
wkt &＃61; &＃39;POINT (%f %f)&＃39; % (point_coors[0],point_coors[1])
#print(wkt)
geom &＃61; ogr.CreateGeometryFromWkt(wkt)
feat &＃61; ogr.Feature(lyrn.GetLayerDefn())
feat.SetField(&＃39;idx&＃39;,idxx)
feat.SetField(&＃39;namex&＃39;,&＃39;ID{0}&＃39;.format(idxx))
feat.SetGeometry(geom)
lyrn.CreateFeature(feat)
#print(lyrn.GetLayerDefn())
#print(&＃39;idx:%i&＃39; % (idxx))
#print("namex:%s" % &＃39;ID{0}&＃39;.format(idxx))
#x &＃61; &＃39;ID{0}&＃39;.format(idxx)
#print(x)
"""
wkt &＃61; &＃39;POINT (%f %f)&＃39; % (point_coors[0],point_coors[1])
geom &＃61; ogr.CreateGeometryFromWkt(wkt)
feat &＃61; ogr.Feature(lyrn.GetLayerDefn())
feat.SetField(&＃39;idx&＃39;,idxx)
feat.SetField(&＃39;namex&＃39;,&＃39;ID{0}&＃39;.format(idxx))
feat.SetGeometry(geom)
lyrn.CreateFeature(feat)
"""
newds.Destroy()
extfile &＃61; &＃39;xx_data_line.shp&＃39;
driver &＃61; ogr.GetDriverByName("ESRI Shapefile")
if os.access(extfile,os.F_OK):
driver.DeleteDataSource(extfile)
newds &＃61; driver.CreateDataSource(extfile)
lyrn &＃61; newds.CreateLayer(&＃39;line&＃39;,None,ogr.wkbLineString)
#字段定义创建
fieldcnstr &＃61; ogr.FieldDefn("id",ogr.OFTInteger)
fieldf &＃61; ogr.FieldDefn("name",ogr.OFTString)
lyrn.CreateField(fieldcnstr)
lyrn.CreateField(fieldf)
point_coors_arr &＃61; [[0,0,1,2,3,-2,6,0]]
#print(point_coors_arr)
#print(type(point_coors_arr))
#print(len(point_coors_arr))
for idx,point_coors in enumerate(point_coors_arr):
#wkt &＃61; &＃39;LINESTRING(%f %f,%f %f,%f %f,%f %f)&＃39; % (point_coors[len(point_coors_arr) - len(point_coors_arr)],point_coors[len(point_coors_arr) - len(point_coors_arr) &＃43; 1],)
wkt &＃61; &＃39;LINESTRING(%f %f,%f %f,%f %f,%f %f)&＃39; % (point_coors[0],point_coors[1],point_coors[2],point_coors[3],point_coors[4],point_coors[5]
,point_coors[6],point_coors[7])
print(wkt)
geom &＃61; ogr.CreateGeometryFromWkt(wkt)
feat &＃61; ogr.Feature(lyrn.GetLayerDefn())
feat.SetField(&＃39;id&＃39;,idx)
feat.SetField(&＃39;name&＃39;,&＃39;line_one&＃39;)
feat.SetGeometry(geom)
lyrn.CreateFeature(feat)
newds.Destroy()
extfile &＃61; &＃39;xx_data_polygon.shp&＃39;
driver &＃61; ogr.GetDriverByName(&＃39;ESRI Shapefile&＃39;)
if os.access(extfile,os.F_OK):
driver.DeleteDataSource(extfile)
#创建数据源文件
newds &＃61; driver.CreateDataSource(extfile)
#创建几何形状图层
lyrn &＃61; newds.CreateLayer(&＃39;polygon&＃39;,None,ogr.wkbPolygon)
#定义字段
fieldcnstr &＃61; ogr.FieldDefn(&＃39;id&＃39;,ogr.OFTInteger)
fieldf &＃61; ogr.FieldDefn(&＃39;name&＃39;,ogr.OFTString)
#图层添加字段
lyrn.CreateField(fieldcnstr)
lyrn.CreateField(fieldf)
wkt_polygon_1 &＃61; &＃39;POLYGON((2 1,12 1,12 4,2 4,1 2))&＃39;
wkt_polygon_2 &＃61; &＃39;POLYGON((4 1,8 1,8 3,4 3,3 1))&＃39;
wkt_polygon_3 &＃61; &＃39;POLYGON((8 4,10 4, 10 5,8 5,6 4))&＃39;
#print(type(wkt_polygon_1))
point_coors_arr &＃61; [wkt_polygon_1,wkt_polygon_2,wkt_polygon_3]
#print(type(point_coors_arr))
for idx,point_coors in enumerate(point_coors_arr):
#print(idx)
#print(point_coors)
wkt &＃61; point_coors
#使用wkt创建几何图形
geom &＃61; ogr.CreateGeometryFromWkt(wkt)
#获取图层要素 ogr.GetLayerDefn , ogr.Feature()
feat &＃61; ogr.Feature(lyrn.GetLayerDefn())
feat.SetField(&＃39;id&＃39;,idx)
print(&＃39;poly_{idx}&＃39;.format(idx &＃61; idx))
#feat.SetField(&＃39;name&＃39;,&＃39;poly_{idx}&＃39;.format(idx &＃61; idx))
feat.SetField(&＃39;name&＃39;,&＃39;poly_{idx}&＃39;.format(idx &＃61; idx))
feat.SetGeometry(geom)
lyrn.CreateFeature(feat)
newds.Destroy()
from osgeo import ogr
import os,math
inshp &＃61; "xx_data_polygon.shp"
ds &＃61; ogr.Open(inshp) #打开shp源文件
driver &＃61; ogr.GetDriverByName("ESRI Shapefile")
outputfile &＃61; "xx_data_polygon_copy.shp" #输出shp
if os.access(outputfile,os.F_OK):
driver.DeleteDataSource(outputfile)
pt_cp &＃61; driver.CopyDataSource(ds,outputfile)
pt_cp.Release()
from osgeo import ogr
import os,math
inshp &＃61; &＃39;xx_data_polygon.shp&＃39;
ds &＃61; ogr.Open(inshp)
driver &＃61; ogr.GetDriverByName("ESRI Shapefile")
outputfile &＃61; &＃39;cp_polygon.shp&＃39;
if os.access(outputfile,os.F_OK):
driver.DeleteDataSource(outputfile)
pt_cp &＃61; driver.CopyDataSource(ds,outputfile)
pt_cp.Release()
outputfile &＃61; &＃39;cp2.shp&＃39;
if os.access(outputfile,os.F_OK):
driver.DeleteDataSource(outputfile)
newds &＃61; driver.CreateDataSource(outputfile)
layer &＃61; ds.GetLayer()
pt_layer &＃61; newds.CopyLayer(layer,&＃39;xxx&＃39;)
#newds.Destroy()
outputfile &＃61; &＃39;cp3.shp&＃39;
if os.access(outputfile,os.F_OK):
driver.DeleteDataSource(outputfile)
newds &＃61; driver.CreateDataSource(outputfile)
layernew &＃61; newds.CreateLayer(&＃39;worldcopy&＃39;,None,ogr.wkbLineString)
layer &＃61; ds.GetLayer()
feature &＃61; layer.GetNextFeature()
if feature is not None:
layernew.CreateFeature(feature)
feature &＃61; layer.GetNextFeature()
newds.Destroy()
driver &＃61; ogr.GetDriverByName("ESRI Shapefile")
inshp &＃61; &＃39;xx_data_polygon.shp&＃39;
ds &＃61; ogr.Open(inshp)
outf &＃61; &＃39;cxp.shp&＃39;
if os.access(outf,os.F_OK):
driver.DeleteDataSource(outf)
print("exists")
newds &＃61; driver.CreateDataSource(outf)
layernew &＃61; newds.CreateLayer(&＃39;worldcopy&＃39;,None,ogr.wkbLineString)
layer &＃61; ds.GetLayer()
#print(layer)
feature &＃61; layer.GetNextFeature()
print(feature)
while feature is not None:
layernew.CreateFeature(feature)
feature &＃61; layer.GetNextFeature()
newds.Destroy()
&＃39;&＃39;&＃39;
from osgeo import ogr
ds &＃61; ogr.Open(&＃39;convChk.shp&＃39;)
layer &＃61; ds.GetLayer(0)
spatialRef &＃61; layer.GetSpatialRef()
print(spatialRef)
京公网安备 11010802041100号