<div dir="ltr">Dear Gdal Developer,<div><br></div><div>Thank you for reading my email.<br><div><br></div><div><div>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.</div><div><br></div><div>Here is my reprex in Python:</div><div><br></div><div>from osgeo import gdal</div><div><br></div><div>### resample and reproject</div><div>### data are from MOD11A2 and MYD11A2, and have been converted into annual mean values<br>raster_rprj = gdal.Warp("./2015_daytime_mean_re.tif",<br>                        "./2015_daytime_mean_clip.tif", dstSRS = "EPSG:4326",<br>                        xRes = 0.008, yRes = 0.008, resampleAlg = "average")</div><div><br></div><div>raster_rprj = None</div><div><br></div><div>I hope it runs as the function np.nanmean()</div><div><br></div><div>Thank you for your time.</div><div><br></div><div>Best regards,</div><div><br></div><div>Mike Li</div></div></div><div><br></div></div>