我有一个来自 sentinel1 数据的投影 tiff 文件（s1a-ew-grd-hh-20180208t071218-20180208t071323-020512-023164-001.tiff）（我们称之为 test.tiff），其投影为 WKT 格式：
PROJCS["WGS 84",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Polar_Stereographic"],PARAMETER["latitude_of_origin",70],PARAMETER["central_meridian",-45],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]
这应该相当于北极极地立体（epsg 3995）。我的目标是通过插入一些 GCP 并扭曲它来稍微调整图像。为了仔细检查这个方法，我从一个简单的案例开始，我从投影图像 (test.tiff) 中获取角坐标：
Size is 17959, 18039
Corner Coordinates:
Upper Left ( 457598.467, -240151.021) ( 17d18'33.06"E, 85d13'55.33"N)
Lower Left ( 457598.467,-1142101.021) ( 23d 9'56.93"W, 78d40'39.65"N)
Upper Right ( 1355548.467, -240151.021) ( 34d57'13.08"E, 77d20'28.88"N)
Lower Right ( 1355548.467,-1142101.021) ( 4d53' 4.53"E, 73d44'33.76"N)
Center ( 906573.467, -691126.021) ( 7d40'47.73"E, 79d30'17.56"N)
我据此设置GCP。因此：
gdal_translate -gcp 0.0 0.0 457598.467 -240151.021 -gcp 0.0 18039.0 1355548.467 -240151.021 -gcp 17959.0 18039.0 1355548.467 -1142101.021 -gcp 17959.0 0.0 457598.467 -1142101.021 -a_srs WKT_PROJECTION test.tiff gcp.tiff
然后我使用 -tps 选项扭曲图像以使用 GCP 信息强制扭曲。
gdalwarp -tps gcp.tiff gcp_warped.tiff
人们会期望 test.tiff 和 gcp_warped.tiff 文件应该重合，但不幸的是情况并非如此：
在上述过程中有什么我想念的吗？还是有其他解决方法？