Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from osgeo import gdal, gdalconst, osr
- def get_raster_info(src_raster_uri):
- dataset = gdal.Open(src_raster_uri, gdalconst.GA_ReadOnly)
- projection_wkt = dataset.GetProjection()
- srs = osr.SpatialReference(wkt = projection_wkt)
- out_raster_info = {"Description":dataset.GetDescription(),\
- "DriverName":dataset.GetDriver().ShortName,\
- "ProjectionName":srs.GetAttrValue("LOCAL_CS", 0),\
- "ProjectionEPSG":int(srs.GetAttrValue("LOCAL_CS|AUTHORITY", 1)),\
- "DistanceUnit":srs.GetAttrValue("LOCAL_CS|UNIT", 0),\
- "PixelNumber_X,Y,All":(dataset.RasterXSize, dataset.RasterYSize, dataset.RasterXSize * dataset.RasterYSize ),\
- "RasterOrigin_X,Y":(dataset.GetGeoTransform()[0], dataset.GetGeoTransform()[3]),\
- "OnePixelSize_X,Y":(dataset.GetGeoTransform()[1], dataset.GetGeoTransform()[5]),\
- "Rotation":dataset.GetGeoTransform()[2],\
- "Metadata":dataset.GetMetadata()
- }
- return out_raster_info
Add Comment
Please, Sign In to add comment