为什么 gdalwarp 从正交到长纬度会截断极点?这一定是可能的,否则南极洲就不会出现在盘子上

问题描述 投票:0回答:1

为什么从正交到长纬度的 gdalwarp 会截断极点? Lonlat 应该显示 -90 到 +90 它与 eqc 相同,不是吗?它们是 Plate-Carre 的同义词还是我弄错了?为什么会出现差距? 这一定是可能的,否则南极洲就不会出现在板块上。 我是否缺少插值参数?

将正交图像重新投影为等距柱状投影存在某种问题。 两极正在被截断。我缺少一个参数吗? “eqc”“longlat”“plate-carree”之间有区别吗? 这是gdal工具的根本缺陷吗?


gdal_translate 
-of Gtiff 
-a_ullr -1666666 1666666 1666666 -1666666  
-a_srs 
    "+proj=ortho  
    +ellps=WGS84 
    +datum=WGS84 
    +lat_0=85 
    +lon_0=45"  
kat.png 
QQQQ-kat-ortho.tif 



gdalwarp 
-t_srs 
    "+proj=longlat 
    +datum=WGS84 
    +ellps=WGS84 
    +units=m"   
QQQQ-kat-ortho.tif
QQQQ-kat-latlon.tif 

gdal_merge.py  
-o QQQQ-kat-eq-merged.tif     
NE1_50M_SR_W_tenth.tif 
QQQQ-kat-latlon.tif

convert QQQQ-kat-eq-merged.tif QQQQ-kat-eq-merged.png

convert QQQQ-kat-latlon.tif QQQQ-kat-latlon.png

syntax gis gdal map-projections proj
1个回答
0
投票

当地图包含极纬度时,将正交 (ortho) 或透视 (nsper) 图像重新投影为等距柱状投影 (eqc),

需要目标范围参数:“-te -180 -90 +180 +90”

或其他一些 -te 角集,明确地一直延伸到 +/-90 纬度,否则默认 xmax 设置为纬度 85。也许 gdal 确实期望墨卡托并剪掉极点。 ??这不会出现在任何地方的文档中。

示例:

gdalwarp -te -180 -90 +180 +90 -t_srs "+proj=eqc +datum=WGS84 +ellps=WGS84 +units=m" source_ORTHO.tif target_EQC.tif

© www.soinside.com 2019 - 2024. All rights reserved.