GDAL warp 太多点无法转换。它适用于相似的图像和来源

Posted 2023-03-23

【中文标题】GDAL warp 太多点无法转换。它适用于相似的图像和来源

【英文标题】：GDAL warp Too many points failed to transform. It works with similar image and source

【发布时间】：2021-07-18 15:11:54

【问题描述】：

我正在尝试将 netCDF 文件从 GOES Full 磁盘转换为 geotiff，但最后一步出错。

显然，从 NC 生成 tif 的过程有效，并且 fulldisk.tif 已生成但未进行地理参考，我需要将其覆盖在传单地图中。

gdal_translate -ot float32 -unscale -CO COMPRESS=deflate NETCDF:"fulldisk.nc":CMI fulldisk.tif
Input file size is 2500, 1500
0...10...20...30...40...50...60...70...80...90...100 - done.

tiff 可以在这里下载：https://nube1.on.gt/fulldisk.tif

但是当尝试使用以下过程为 Lat,Lon 重新投影时：

gdalwarp -t_srs EPSG:4326 -dstnodata -999.0 fulldisk.tif /var/www/nube1.on.gt/fulldisk_geo.tif

上面写着：

Processing fulldisk.tif [1/1] : 0Using internal nodata values (e.g. -1) for image fulldisk_geo.tif.
ERROR 1: Too many points (529 out of 529) failed to transform, unable to compute output bounds.
Warning 1: Unable to compute source region for output window 0,0,2500,1500, skipping.
...10...20...30...40...50...60...70...80...90...100 - done.

NC 文件是来自 AWS 的 GOES-16 you can download，gdalinfo 输出如下：

Driver: netCDF/Network Common Data Format
Files: fulldisk.nc
Size is 512, 512
Coordinate System is `'
Metadata:
  NC_GLOBAL#cdm_data_type=Image
  NC_GLOBAL#Conventions=CF-1.7
  NC_GLOBAL#dataset_name=OR_ABI-L2-CMIPF-M6C13_G16_s20211141750164_e20211141759483_c20211141759577.nc
  NC_GLOBAL#date_created=2021-04-24T17:59:57.7Z
  NC_GLOBAL#id=71a29d9f-02b2-44b5-81b2-73a3fec6a62e
  NC_GLOBAL#institution=DOC/NOAA/NESDIS > U.S. Department of Commerce, National Oceanic and Atmospheric Administration, National Environmental Satellite, Data, and Information Services
  NC_GLOBAL#instrument_ID=FM1
  NC_GLOBAL#instrument_type=GOES R Series Advanced Baseline Imager
  NC_GLOBAL#iso_series_metadata_id=8c9e8150-3692-11e3-aa6e-0800200c9a66
  NC_GLOBAL#keywords=SPECTRAL/ENGINEERING > INFRARED WAVELENGTHS > BRIGHTNESS TEMPERATURE
  NC_GLOBAL#keywords_vocabulary=NASA Global Change Master Directory (GCMD) Earth Science Keywords, Version 7.0.0.0.0
  NC_GLOBAL#license=Unclassified data.  Access is restricted to approved users only.
  NC_GLOBAL#Metadata_Conventions=Unidata Dataset Discovery v1.0
  NC_GLOBAL#naming_authority=gov.nesdis.noaa
  NC_GLOBAL#orbital_slot=GOES-East
  NC_GLOBAL#platform_ID=G16
  NC_GLOBAL#processing_level=National Aeronautics and Space Administration (NASA) L2
  NC_GLOBAL#production_data_source=Realtime
  NC_GLOBAL#production_environment=OE
  NC_GLOBAL#production_site=NSOF
  NC_GLOBAL#project=GOES
  NC_GLOBAL#scene_id=Full Disk
  NC_GLOBAL#spatial_resolution=2km at nadir
  NC_GLOBAL#standard_name_vocabulary=CF Standard Name Table (v35, 20 July 2016)
  NC_GLOBAL#summary=Single emissive band Cloud and Moisture Imagery Products are digital maps of clouds, moisture and atmospheric windows at IR bands.
  NC_GLOBAL#timeline_id=ABI Mode 6
  NC_GLOBAL#time_coverage_end=2021-04-24T17:59:48.3Z
  NC_GLOBAL#time_coverage_start=2021-04-24T17:50:16.4Z
  NC_GLOBAL#title=ABI L2 Cloud and Moisture Imagery
Subdatasets:
  SUBDATASET_1_NAME=NETCDF:"fulldisk.nc":CMI
  SUBDATASET_1_DESC=[5424x5424] toa_brightness_temperature (16-bit integer)
  SUBDATASET_2_NAME=NETCDF:"fulldisk.nc":DQF
  SUBDATASET_2_DESC=[5424x5424] status_flag (8-bit integer)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)

我使用来自同一站点的另一个非全盘映像测试了相同的过程，它按预期工作，https://nube1.on.gt/conus_geo.tif。我需要全盘域中的中美洲和加勒比地区。

感谢您的指导。

【参考方案1】：

刚刚更换

-dstnodata -999.0

与

-dstnodata -999

强制目的地无数据值为整数，得到

gdalwarp -t_srs EPSG:4326 -dstnodata -999 fulldisk.tif fulldisk_geo.tif
Processing input file fulldisk.tif.
Using internal nodata values (e.g. -1) for image fulldisk.tif.
0...10...20...30...40...50...60...70...80...90...100 - done.

没有错误并生成 168.1 MiB 文件。

【讨论】：

谢谢，由于某种原因，现在 gdalwarp 可以使用 -999 和 -999.0。再次感谢您的尝试。如果其他发现可能与该问题相关，将发布。