У меня есть растровое изображение геотифа, которое имеет систему координат с долготой от 0 до 360. Горизонтальный центр изображения составляет 180 долготы. Смотрите изображение ниже:
Я хочу преобразовать его в EPSG: 4326 SRS с дальностью -180 180. И я хочу, чтобы центр изображения был на гринвичском меридиане (0). Я предполагаю, что этот srs очень широко используется. Я ожидаю, что результат будет выглядеть так:
Поэтому я использую команду gdalwarp для перепроектирования:
gdalwarp -s_srs '+proj=latlong +datum=WGS84 +pm=180dW' -t_srs EPSG:4326 test_col.tif test_4326.tif
Но я получаю только tiff с большими размерами (больше пикселей) и метаданными EPSG: 4326. Само изображение выглядит так же, как и исходное. Но я ожидаю, что это поменяет полушария.
Вопрос в том, как мне сделать изображение строго -180 180 EPSG: 4326 с центром в 0 долготы?
Это gdalinfo моего исходного файла:
Origin = (-0.102272598067084,89.946211604095552)
Pixel Size = (0.204545196134167,-0.204423208191126)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -0.1022726, 89.9462116) ( 0d 6' 8.18"W, 89d56'46.36"N)
Lower Left ( -0.1022726, -89.9462116) ( 0d 6' 8.18"W, 89d56'46.36"S)
Upper Right ( 359.897, 89.946) (359d53'50.18"E, 89d56'46.36"N)
Lower Right ( 359.897, -89.946) (359d53'50.18"E, 89d56'46.36"S)
Center ( 179.8975000, -0.0000000) (179d53'51.00"E, 0d 0' 0.00"S)
Это gdalinfo после gdalwarp
Origin = (-180.102727401932952,89.946211604095552)
Pixel Size = (0.091397622896436,-0.091420837939082)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left (-180.1027274, 89.9462116) (180d 6' 9.82"W, 89d56'46.36"N)
Lower Left (-180.1027274, -89.9699975) (180d 6' 9.82"W, 89d58'11.99"S)
Upper Right ( 179.8211116, 89.9462116) (179d49'16.00"E, 89d56'46.36"N)
Lower Right ( 179.8211116, -89.9699975) (179d49'16.00"E, 89d58'11.99"S)
Center ( -0.1408079, -0.0118929) ( 0d 8'26.91"W, 0d 0'42.81"S)
-tr xres yres
флаг?