Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- def CreateGeoTiff(Name, Array, driver, NDV,
- xsize, ysize, GeoT, Projection, DataType):
- if DataType == 'Float32':
- DataType = gdal.GDT_Float32
- NewFileName = Name+'.tif'
- # Set nans to the original No Data Value
- Array[np.isnan(Array)] = NDV
- # Set up the dataset
- DataSet = driver.Create( NewFileName, xsize, ysize, 1, DataType)
- # the '1' is for band 1.
- DataSet.SetGeoTransform(GeoT)
- DataSet.SetProjection( Projection.ExportToWkt() )
- # Write the array
- DataSet.GetRasterBand(1).WriteArray( Array )
- DataSet.GetRasterBand(1).SetNoDataValue(NDV)
- return NewFileName
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement