为什么从正交到长纬度的 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
当地图包含极纬度时,将正交 (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