gdal.ComputeProximity始终返回0

Question

我尝试创建一个光栅层，使其与具有目标值的单元格之间的距离最小，因此我使用gdal.ComputeProximity（）

但我必须做些错事，因为我的输出光栅只包含0个值，我不知道哪里是我的错误

先谢谢你

我的数据特征： uint8光栅，994列，1320行 没有数据值是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"])