Creating georeference geotiff image

Question

I have tried to create a geotiff image after extracting texture feature of a grayscale geotiff image using the below codes.

#Creating texture images
Texture_Contrast = Image.fromarray(Contrast)

#image coordinates corner
(left_x_min,bottom_y_min,right_x_max,top_y_max) = im_geo.bounds

# Obtain the projection system
CRS = im_geo.crs

# set geotransform
nx = width
ny = height
xres = (right_x_max - left_x_min) / float(nx)
yres = (top_y_max - bottom_y_min) / float(ny)
geotransform = (left_x_min, xres, 0, top_y_max , 0, -yres)

dst_format = 'Gtiff'
dst_datatype = gdal.GDT_Byte
dst_options = ['COMPRESS=LZW']
dst_file = 'myGEOTIFF.tif'
dst_nbands = 1

driver = gdal.GetDriverByName(dst_format)
dst_ds = driver.Create(dst_file, nx, ny, dst_nbands, dst_datatype, dst_options)
dst_ds.GetRasterBand(1).WriteArray(Texture_Contrast)

dst_ds.SetGeoTransform(geotransform)    # specify coords
srs = osr.SpatialReference()            # establish encoding
srs.ImportFromEPSG(32752)                # WGS84 lat/long
dst_ds.SetProjection(srs.ExportToWkt()) # export coords to file

dst_ds.GetRasterBand(1).WriteArray(Texture_Contrast)   # write r-band to the raster

dst_ds.FlushCache()                     # write to disk
dst_ds = None

However, after running the WriteArray(Texture_Contrast) the following error was arise.

  Input In [47] in <module>
    out_band.WriteArray(Texture_Contrast)

  File ~\miniconda3\envs\spyder-env\lib\site-packages\osgeo\gdal.py:4090 in WriteArray
    return gdal_array.BandWriteArray(self, array, xoff, yoff,

  File ~\miniconda3\envs\spyder-env\lib\site-packages\osgeo\gdal_array.py:509 in BandWriteArray
    if array is None or len(array.shape) != 2:

  File ~\miniconda3\envs\spyder-env\lib\site-packages\PIL\Image.py:546 in __getattr__
    raise AttributeError(name)

AttributeError: shape

Could anyone give an enlightenment and the solution of this? Thank you in advance.

Answer 1

Are the shapes (a,b) and nx, ny of the array you extracted the same? nx and ny must be the same shape as the array you created.

or

You said you wrote a gray image, but the result of Texture_Contrast.shape might be (x, y, z). If so, you need to select a channel and make it 2dim.

The error message seems to be the case below. Please confirm it.