发布于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
作者:goodbody
链接:https://www.pythonheidong.com/blog/article/234103/9d0dbe0c6433c0e2c472/
来源:python黑洞网
任何形式的转载都请注明出处,如有侵权 一经发现 必将追究其法律责任
昵称:
评论内容:(最多支持255个字符)
---无人问津也好,技不如人也罢,你都要试着安静下来,去做自己该做的事,而不是让内心的烦躁、焦虑,坏掉你本来就不多的热情和定力
Copyright © 2018-2021 python黑洞网 All Rights Reserved 版权所有,并保留所有权利。 京ICP备18063182号-1
投诉与举报,广告合作请联系vgs_info@163.com或QQ3083709327
免责声明:网站文章均由用户上传,仅供读者学习交流使用,禁止用做商业用途。若文章涉及色情,反动,侵权等违法信息,请向我们举报,一经核实我们会立即删除!