本站消息
python gdal根据坐标（经纬度）创建空的矢量
发布于2020-02-26 11:58 阅读(1590) 评论(0) 点赞(22) 收藏(5)
前面的博文已经出现过这段代码了，这里单独拿出来吧，方便使用，坐标只需要知道左上角和右下角就行
# -*- coding: utf-8 -*-
import gdal
from osgeo import ogr, osr
shp_path = "./mask_use.shp"
driver = ogr.GetDriverByName("ESRI Shapefile")
data_source = driver.CreateDataSource(shp_path)
srs = osr.SpatialReference()
srs.ImportFromEPSG(4326) #这是WGS84,想用什么自己去搜下对应的编码就行了
layer = data_source.CreateLayer("polygon", srs, ogr.wkbPolygon)
feature = ogr.Feature(layer.GetLayerDefn())
wa = 116.122741699
ha = 40.080871582
wa1 = 116.122741699
ha1 = 40.168762207
wa2 = 116.251831055
ha2 = 40.168762207
wa3 = 116.251831055
ha3 = 40.080871582
wkt = "POLYGON((" + str(wa)+ " " +str(ha)+ "," + str(wa1) + " " + str(ha1) + "," + str(wa2)+ " " +str(ha2)+ "," + str(wa3)+ " " +str(ha3) + "))"
point = ogr.CreateGeometryFromWkt(wkt)
feature.SetGeometry(point)
layer.CreateFeature(feature)
feature = None
data_source = None
- 1
- 2
- 3
- 4
- 5
- 6
- 7
- 8
- 9
- 10
- 11
- 12
- 13
- 14
- 15
- 16
- 17
- 18
- 19
- 20
- 21
- 22
- 23
- 24
- 25
所属网站分类: 技术文章 > 博客
作者：goodbody
链接：https://www.pythonheidong.com/blog/article/234103/9d0dbe0c6433c0e2c472/
来源：python黑洞网
任何形式的转载都请注明出处,如有侵权 一经发现 必将追究其法律责任
22
0
收藏该文
昵称:
评论内容：(最多支持255个字符)
---无人问津也好，技不如人也罢，你都要试着安静下来，去做自己该做的事，而不是让内心的烦躁、焦虑，坏掉你本来就不多的热情和定力
站长公众号(new)
更多>
pdf(new)
更多>
问答(new)
更多>
游戏(new)
更多>