GDALWARP如何输出ISEA（二十面体斯奈德等积）投影栅格？

Question

使用 gdalwarp 生成 ISEA“二十面体斯奈德等面积”GEOTIFF/光栅。
使用 PROJ/GDAL 无法生成二十面体投影吗？

使用 gdalwarp 生成 ISEA“二十面体斯奈德等面积”GEOTIFF/光栅。
我正在尝试使用 gdal 工具将 Plate-Carree 转换为 ISEA 并返回。

输入 GEOTIFF 为“ NE1_50M_SR_W.tif ”
直接来自 Natural Earth Downloads 的“The Natural Earth 快速入门套件”。

我找不到满足此错误消息的输入格式或参数排列。

gdalwarp  -te -180 -90 +180 +90     -t_srs "+proj=isea +elipse=sphere" NE1_50M_SR_W.tif out-isea.tif
ERROR 1: Too many points (10201 out of 10201) failed to transform, unable to compute output bounds.

gdalinfo NE1_50M_SR_W.tif 
Driver: GTiff/GeoTIFF
Files: NE1_50M_SR_W.tif
Size is 10800, 5400
Coordinate System is:
GEOGCRS["WGS 84",
    DATUM["World Geodetic System 1984",
        ELLIPSOID["WGS 84",6378137,298.257223563,
            LENGTHUNIT["metre",1]]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["degree",0.0174532925199433]],
    CS[ellipsoidal,2],
        AXIS["geodetic latitude (Lat)",north,
            ORDER[1],
            ANGLEUNIT["degree",0.0174532925199433]],
        AXIS["geodetic longitude (Lon)",east,
            ORDER[2],
            ANGLEUNIT["degree",0.0174532925199433]],
    ID["EPSG",4326]]
Data axis to CRS axis mapping: 2,1
Origin = (-179.999999999999972,90.000000000000000)
Pixel Size = (0.033333333333330,-0.033333333333330)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_DATETIME=2014:10:18 09:32:38
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_SOFTWARE=Adobe Photoshop CC 2014 (Macintosh)
  TIFFTAG_XRESOLUTION=342.85699
  TIFFTAG_YRESOLUTION=342.85699
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (-180.0000000,  90.0000000) (180d 0' 0.00"W, 90d 0' 0.00"N)
Lower Left  (-180.0000000, -90.0000000) (180d 0' 0.00"W, 90d 0' 0.00"S)
Upper Right ( 180.0000000,  90.0000000) (180d 0' 0.00"E, 90d 0' 0.00"N)
Lower Right ( 180.0000000, -90.0000000) (180d 0' 0.00"E, 90d 0' 0.00"S)
Center      (  -0.0000000,   0.0000000) (  0d 0' 0.00"W,  0d 0' 0.00"N)
Band 1 Block=10800x1 Type=Byte, ColorInterp=Red
Band 2 Block=10800x1 Type=Byte, ColorInterp=Green
Band 3 Block=10800x1 Type=Byte, ColorInterp=Blue
molly@morfydd:~/Documents/earthgeo/gdal/scripttest$ 

Answer 1

再次尝试使用新的 GDAL 和 新的 PROJ 实现，版本 9.5 或更高版本。它说“可用形式：正向和反向，球形”。

逆是一个多年的问题，最终由 Jerome St-Louis 在 2024 年 8 月的拉取请求中实现。