Might be fixed by [GDALSieveFilter](https://gdal.org/api/gdal_alg.html#_CPPv415GDALSieveFilter15GDALRasterBandH15GDALRasterBandH15GDALRasterBandHiiPPc16GDALProgressFuncPv) like: ```python from osgeo import gdal Image = gdal.Open('SomeImageName.tif', 1) # open image in read-write mode Band = Image.GetRasterBand(1) gdal.SieveFilter(srcBand=Band, maskBand=None, dstBand=Band, threshold=100, connectedness=8, callback=gdal.TermProgress_nocb) del Image, Band # close the datasets. ```