Guest User

Untitled

a guest
Jul 20th, 2018
99
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.12 KB | None | 0 0
  1. from osgeo import gdal, gdalconst, osr
  2.  
  3. def get_raster_info(src_raster_uri):
  4.  
  5. dataset = gdal.Open(src_raster_uri, gdalconst.GA_ReadOnly)
  6.  
  7. projection_wkt = dataset.GetProjection()
  8. srs = osr.SpatialReference(wkt = projection_wkt)
  9.  
  10. out_raster_info = {"Description":dataset.GetDescription(),\
  11. "DriverName":dataset.GetDriver().ShortName,\
  12. "ProjectionName":srs.GetAttrValue("LOCAL_CS", 0),\
  13. "ProjectionEPSG":int(srs.GetAttrValue("LOCAL_CS|AUTHORITY", 1)),\
  14. "DistanceUnit":srs.GetAttrValue("LOCAL_CS|UNIT", 0),\
  15. "PixelNumber_X,Y,All":(dataset.RasterXSize, dataset.RasterYSize, dataset.RasterXSize * dataset.RasterYSize ),\
  16. "RasterOrigin_X,Y":(dataset.GetGeoTransform()[0], dataset.GetGeoTransform()[3]),\
  17. "OnePixelSize_X,Y":(dataset.GetGeoTransform()[1], dataset.GetGeoTransform()[5]),\
  18. "Rotation":dataset.GetGeoTransform()[2],\
  19. "Metadata":dataset.GetMetadata()
  20. }
  21.  
  22. return out_raster_info
Add Comment
Please, Sign In to add comment