gis - GDAL warp 太多点无法转换。它适用于相似的图像和来源
问题描述
我正在尝试将 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,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 文件。
推荐阅读
- python - Aws Lambda + Django + django-storages + collectstatic 不一致的行为
- javascript - sql while循环.net通过隐藏变量解析值到javascript函数以创建图表
- php - 具有超过 1 个 IN 函数的安全 pdo mysql 语句
- java - PowerMock 在不使用 @PrepareForTest 注释的情况下为测试准备一个类
- python - 如何在按钮中制作图像?
- c# - 将 json 读取到字符串时遇到问题
- python - 如何使用 Python 提取 1000 条推文?
- amazon-web-services - AWS CLI(aws rds) 未知选项:--cloudwatch-logs-export-configuration
- python - While循环减去
- mysql - TIMEOUT in Laravel