[gdal-dev] Issues in Raster Resampling

[Even Rouault]

Mike,

if you look at the nodata value set on 2015_daytime_mean_clip.tif, it is 
at 3.39999995214436425e+38, and not at nan.

If you add the srcNodata=math.nan argument to gdal.Warp(), this will fix 
your issue.

Even

Le 23/04/2022 à 06:40, Chao Li a écrit :
> Dear Gdal Developer,
>
> Thank you for reading my email.
>
> I used gdal.Wrap() to resample from a high resolution to a lower. 
> However, my raster has no value (nan). So, when I set resampleAlg, the 
> larger grids with nan(s) will become nan.
>
> Here is my reprex in Python:
>
> from osgeo import gdal
>
> ### resample and reproject
> ### data are from MOD11A2 and MYD11A2, and have been converted into 
> annual mean values
> raster_rprj = gdal.Warp("./2015_daytime_mean_re.tif",
>                         "./2015_daytime_mean_clip.tif", dstSRS = 
> "EPSG:4326",
>                         xRes = 0.008, yRes = 0.008, resampleAlg = 
> "average")
>
> raster_rprj = None
>
> I hope it runs as the function np.nanmean()
>
> Thank you for your time.
>
> Best regards,
>
> Mike Li
>
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
http://www.spatialys.com
My software is free, but my time generally not.