我正在使用gdalwarp在日期线附近(即180°,又称antimeridian)操纵SRTM磁贴。SRTM磁贴与子午线的重叠非常小(1/2像素)。您可以使用gdalinfo查看此内容:
gdalinfo S16W180.hgt
Driver: SRTMHGT/SRTMHGT File Format
Files: S16W180.hgt
Size is 1201, 1201
[...]
Lower Left (-180.0004167, -16.0004167) (180d 0' 1.50"W, 16d 0' 1.50"S)
Upper Right (-178.9995833, -14.9995833) (178d59'58.50"W, 14d59'58.50"S)
[...]
因此,源跨越了日期线很小的一部分。
这会导致gdalwarp出现问题,最终会产生巨大的全球性产出。
gdalwarp -t_srs "epsg:900913" S16W180.hgt test.tif
gdalinfo test.tif
Driver: GTiff/GeoTIFF
Files: test.tif
Size is 1703, 5
[...]
Lower Left (-20037508.330,-1806798.473) (180d 0' 0.00"W, 16d 7'13.00"S)
Upper Right (20032839.451,-1689152.120) (179d57'29.01"E, 15d 5'45.84"S)
请注意,经度几乎覆盖了整个地球,而且行数出乎意料的小(5)
这是gdalwarp中的错误吗?如果没有,传递给gdalwarp以获得合理输出的正确选择是什么?
dds.cr.usgs.gov/srtm/version2_1/SRTM3/Australia/S16W180.hgt.zip,以备您尝试。
—
重力风暴2012年
添加SOURCE_EXTRA参数,请参见code.google.com/p/maptiler/issues/detail?id=6-试试gdalwarp -t_srs epsg:900913 -wo SOURCE_EXTRA = 120 S16W180.hgt test.tif
—
Mapperz
也许使用-te参数作为“目标范围”,或者首先使用gdal_translate和a_ullr来覆盖现有范围,或者使用-projwin来缩小范围内想要的位,以固定范围
—
mdsumner