wkt raster

e86c5ba7 · Mario Chirinos Colunga · 6403b201 · e86c5ba7
Commit e86c5ba7 authored 6 years ago by Mario Chirinos Colunga
Show whitespace changes
Inline Side-by-side

Showing with 1 addition and 1 deletion

rasterWkt.py geosentinel/rasterWkt.py +1 -1

No files found.
--- a/geosentinel/rasterWkt.py
+++ b/geosentinel/rasterWkt.py
@@ -37,7 +37,7 @@ def getPolygonArray(image, wkt_geom):


 	srs = osr.SpatialReference()
-	srs.ImportFromEPSG(4326)
+	srs.ImportFromWkt(image.GetProjection())
 	layer = data_source.CreateLayer("wkt", srs, ogr.wkbPolygon)
 	geometry = ogr.CreateGeometryFromWkt(wkt_geom)
 	feature = ogr.Feature(layer.GetLayerDefn())