gdal.ComputeProximity始终返回0

Question

我尝试创建距离目标值具有最小距离的栅格图层，因此我使用gdal.ComputeProximity（）。

但是我必须做错什么，因为我的输出栅格仅包含0个值，并且看不到我的错误在哪里。

谢谢。

我的数据的特征： 具有994列和1320行的uint8栅格 no_data值为0 espg：2154（公制）

from osgeo import gdal

raster1 = r'/path/to/my/raster/os_1951.tif'

r1 = gdal.Open(raster1)
r1band=r1.GetRasterBand(1)

r2 = r'/path/to/my/raster/dist13_1951_t2.tif'
drv = gdal.GetDriverByName('GTiff')
r2 = drv.Create( r2,
                r1.RasterXSize, r1.RasterYSize, 1,
                gdal.GetDataTypeByName('Float32'))
r2.SetGeoTransform(r1.GetGeoTransform())
r2.SetProjection(r1.GetProjectionRef())

r2band = r2.GetRasterBand(1)
gdal.ComputeProximity(r1band,r2band,["VALUES=13","DISTUNITS=GEO"])