def transform_geom(wkt, srid_in, srid_out):
proj_in = SpatialReference(int(srid_in))
proj_out = SpatialReference(int(srid_out))
ogr = OGRGeometry(wkt)
if hasattr(ogr, 'srs'):
ogr.srs = proj_in
else:
ogr.set_srs(proj_in)
ogr.transform_to(proj_out)
return ogr.wkt
def transform_geom(wkt, srid_in, srid_out):
proj_in = SpatialReference(int(srid_in))
proj_out = SpatialReference(int(srid_out))
ogr = OGRGeometry(wkt)
if hasattr(ogr, 'srs'):
ogr.srs = proj_in
else:
ogr.set_srs(proj_in)
ogr.transform_to(proj_out)
return ogr.wkt
def get_extent_from_text(points, srid_in, srid_out):
"""Transform an extent from srid_in to srid_out."""
proj_in = SpatialReference(srid_in)
proj_out = SpatialReference(srid_out)
if srid_out == 900913:
if int(float(points[0])) == -180:
points[0] = -179
if int(float(points[1])) == -90:
points[1] = -89
if int(float(points[2])) == 180:
points[2] = 179
if int(float(points[3])) == 90:
points[3] = 89
wkt = 'POINT(%f %f)' % (float(points[0]), float(points[1]))
wkt2 = 'POINT(%f %f)' % (float(points[2]), float(points[3]))
ogr = OGRGeometry(wkt)
ogr2 = OGRGeometry(wkt2)
if hasattr(ogr, 'srs'):
ogr.srs = proj_in
ogr2.srs = proj_in
else:
ogr.set_srs(proj_in)
ogr2.set_srs(proj_in)
ogr.transform_to(proj_out)
ogr2.transform_to(proj_out)
wkt = ogr.wkt
wkt2 = ogr2.wkt
mins = wkt.replace('POINT (', '').replace(')', '').split(' ')
maxs = wkt2.replace('POINT (', '').replace(')', '').split(' ')
mins.append(maxs[0])
mins.append(maxs[1])
return mins
def get_extent_from_text(points, srid_in, srid_out):
"""Transform an extent from srid_in to srid_out."""
proj_in = SpatialReference(srid_in)
proj_out = SpatialReference(srid_out)
if srid_out == 900913:
if int(float(points[0])) == -180:
points[0] = -179
if int(float(points[1])) == -90:
points[1] = -89
if int(float(points[2])) == 180:
points[2] = 179
if int(float(points[3])) == 90:
points[3] = 89
wkt = 'POINT(%f %f)' % (float(points[0]), float(points[1]))
wkt2 = 'POINT(%f %f)' % (float(points[2]), float(points[3]))
ogr = OGRGeometry(wkt)
ogr2 = OGRGeometry(wkt2)
if hasattr(ogr, 'srs'):
ogr.srs = proj_in
ogr2.srs = proj_in
else:
ogr.set_srs(proj_in)
ogr2.set_srs(proj_in)
ogr.transform_to(proj_out)
ogr2.transform_to(proj_out)
wkt = ogr.wkt
wkt2 = ogr2.wkt
mins = wkt.replace('POINT (', '').replace(')', '').split(' ')
maxs = wkt2.replace('POINT (', '').replace(')', '').split(' ')
mins.append(maxs[0])
mins.append(maxs[1])
return mins
