问题描述
我正在尝试将 netCDF 文件从 GOES 完整磁盘转换为 geotiff,但在最后一步出错。
显然,从 NC 作品和 fulldisk.tif 生成 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,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。我需要全磁盘域中的中美洲和加勒比海地区。
感谢您的指导。
解决方法
刚换的
-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。