test.py
751 Bytes
# coding=utf-8
#author: 4N
#createtime: 2022/3/15
#email: nheweijun@sina.com
from osgeo import ogr
from osgeo.ogr import *
import math
import os
class util:
@classmethod
def envelop_2_polygon(cls,env):
ring = ogr.Geometry(ogr.wkbLinearRing)
ring.AddPoint(env[0], env[2])
ring.AddPoint(env[0], env[3])
ring.AddPoint(env[1], env[3])
ring.AddPoint(env[1], env[2])
ring.AddPoint(env[0], env[2])
# Create polygon
poly = ogr.Geometry(ogr.wkbPolygon)
poly.AddGeometry(ring)
return poly
p1= util.envelop_2_polygon([0,1,0,1])
p2 = util.envelop_2_polygon([2,3,2,3])
lll = p1.GetGeometryRef(0)
print(lll.GetPoints())