python gdal根據座標（經緯度）創建空的矢量

前面的博文已經出現過這段代碼了，這裏單獨拿出來吧，方便使用，座標只需要知道左上角和右下角就行

# -*- 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