Hi, It seems to work for me with GDAL 3.6.0dev on Windows.
gdalwarp -t_srs EPSG:4326 -r average -ot float64 -multi MOD17A2H_GPP_2020_001.tif mod64m.tif Creating output file that is 2969P x 2572L. Processing MOD17A2H_GPP_2020_001.tif [1/1] : 0Using internal nodata values (e.g. 32767) for image MOD17A2H_GPP_2020_001.tif. Copying nodata values from source MOD17A2H_GPP_2020_001.tif to destination mod64m.tif. ...10...20...30...40...50...60...70...80...90...100 - done. gdalinfo mod64m.tif Driver: GTiff/GeoTIFF Files: mod64m.tif Size is 2969, 2572 Coordinate System is: GEOGCRS["WGS 84", … Band 1 Block=2969x1 Type=Float64, ColorInterp=Gray NoData Value=32767 I don’t know why nodata is lost for you. Have you already tried to set the output nodata explicitly with -dstnodata 32767 as a workaround? -Jukka Rahkonen- Lähettäjä: Chao Li <[email protected]> Lähetetty: maanantai 27. kesäkuuta 2022 10.50 Vastaanottaja: Rahkonen Jukka <[email protected]> Aihe: Re: [gdal-dev] gdalwarp Resampling Input Raster Nodata Setting Dear Ms./Mr. Jukka Thank you for your reply. I try to delete the -srcnodata 32767, but it still doesn't work. I have attach a part of raster, in this email. Thanks for your help Best regards Rahkonen Jukka <[email protected]<mailto:[email protected]>> 于2022年6月27日周一 14:52写道： Hi, The right syntax is average -srcnodata 32767 but you should not need to use it at all because GDAL recognizes nodata from your input automatically as you can see from the gdalinfo report. Could you clip and share a small sample from your source data so that some nodata is included? You can use gdal_translate with -srswin or -projwin for that. -Jukka Rahkonen- Lähettäjä: gdal-dev <[email protected]<mailto:[email protected]>> Puolesta Chao Li Lähetetty: maanantai 27. kesäkuuta 2022 4.09 Vastaanottaja: [email protected]<mailto:[email protected]> Aihe: [gdal-dev] gdalwarp Resampling Input Raster Nodata Setting Dear GDAL experts, I have a problem using gdalwarp. I want to reproject a raster and change its resolution by the average method. This is the gdalinfo of the raster [q70176a@ito-2 ~]$ gdalinfo DP12/11_RawData/02_GGP/Gross_PP_8Days_500m_v6/GPP/MOD17A2H_GPP_2016_001.tif Driver: GTiff/GeoTIFF Files: DP12/11_RawData/02_GGP/Gross_PP_8Days_500m_v6/GPP/MOD17A2H_GPP_2016_001.tif Size is 86400, 31200 Coordinate System is: PROJCRS["unknown", BASEGEOGCRS["unknown", DATUM["unknown", ELLIPSOID["unknown",6371007.181,0, LENGTHUNIT["metre",1, ID["EPSG",9001]]]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433, ID["EPSG",9122]]]], CONVERSION["Sinusoidal", METHOD["Sinusoidal"], PARAMETER["Longitude of natural origin",0, ANGLEUNIT["degree",0.0174532925199433], ID["EPSG",8802]], PARAMETER["False easting",0, LENGTHUNIT["metre",1], ID["EPSG",8806]], PARAMETER["False northing",0, LENGTHUNIT["metre",1], ID["EPSG",8807]]], CS[Cartesian,2], AXIS["easting",east, ORDER[1], LENGTHUNIT["metre",1, ID["EPSG",9001]]], AXIS["northing",north, ORDER[2], LENGTHUNIT["metre",1, ID["EPSG",9001]]]] Data axis to CRS axis mapping: 1,2 Origin = (-20015109.353999998420477,7783653.637667000293732) Pixel Size = (463.312716527771215,-463.312716527787586) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left (-20015109.354, 7783653.638) (166d17' 5.25"W, 70d 0' 0.00"N) Lower Left (-20015109.354,-6671703.118) ( 0d 0' 0.00"E, 60d 0' 0.00"S) Upper Right (20015109.354, 7783653.638) (166d17' 5.25"E, 70d 0' 0.00"N) Lower Right (20015109.354,-6671703.118) ( 0d 0' 0.00"W, 60d 0' 0.00"S) Center ( -0.000, 555975.260) ( 0d 0' 0.00"W, 5d 0' 0.00"N) Band 1 Block=86400x1 Type=Int16, ColorInterp=Gray NoData Value=32767 This is the gdalwarp commend: [q70176a@ito-2 ~]$ gdalwarp -t_srs EPSG:4326 -te -180 -90 180 90 -tr 0.1 0.1 -ot Float64 -r average -srcnodata "32767" -multi -overwrite DP12/11_RawData/02_GGP/Gross_PP_8Days_500m_v6/GPP/MOD17A2H_GPP_2016_001.tif DP12/11_RawData/02_GGP/Gross_PP_8Days_500m_v6/GPP_01_test/test2.tif I also tried [q70176a@ito-2 ~]$ gdalwarp -t_srs EPSG:4326 -te -180 -90 180 90 -tr 0.1 0.1 -ot Float64 -r average -srcnodata 32767 -multi -overwrite DP12/11_RawData/02_GGP/Gross_PP_8Days_500m_v6/GPP/MOD17A2H_GPP_2016_001.tif DP12/11_RawData/02_GGP/Gross_PP_8Days_500m_v6/GPP_01_test/test2.tif However, the 32767 is not set to be nodata, and be averaged in the output raster. How can I set this? Additionally, I am using a HPC with multinodes. With 8 nodes, I have 288. If I set -multi and -wo NUM_THREADS=val/ALL_CPUS, will all that 288 run together? Since there nothing returns, I really do not know. Best regards. -- ------------------------------------------------------------------------------------ Mike Li Department of Urban and Environmental Engineering Graduate School of Engineering Kyushu University 744 Motooka, Nishi-ku, Fukuoka 819-0395 Japan Tel: 090-8304-8953 E-mail: [email protected]<mailto:[email protected]> ------------------------------------------------------------------------------------- 李潮（リ チョウ） 九州大学 大学院 工学府 都市環境システム工学専攻 都市工学研究室 ------------------------------------------------------------------------------------- -- ------------------------------------------------------------------------------------ Mike Li Department of Urban and Environmental Engineering Graduate School of Engineering Kyushu University 744 Motooka, Nishi-ku, Fukuoka 819-0395 Japan Tel: 090-8304-8953 E-mail: [email protected]<mailto:[email protected]> ------------------------------------------------------------------------------------- 李潮（リ チョウ） 九州大学 大学院 工学府 都市環境システム工学専攻 都市工学研究室 -------------------------------------------------------------------------------------
_______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
