Форумы GIS-Lab.info

На входе - растр lat/long. Нужно получить меркатор epsg:3857.
Сначала через gdal_translate -gcp... привязываю исходный растр.
Потом через gdalwarp получаю растр в меркаторе.
Всё работает, только если в исходном растре нету 180 меридиана, т.е. перехода из E в W (из конца проекции в ее начало).
Разделить исходные данные на W и E не получится, ведь 180 меридиан проходит в исходном растре проходит под углом.
Как с этим бороться? Как мне получить два (W и E) корректных растра в меркаторе?

На входе - растр lat/long.

ведь 180 меридиан проходит в исходном растре проходит под углом

Это как?

Если под углом, то это не широта/долгота, а какая-то географическая проекция...

1. А если есть 180 меридиан, что происходит?
2. В какой проекции исходный растр?
3. Нельзя ли как промежуточный этап использовать проекцию Меркатора с центральным меридианом, отличным от 0?
4. Вы уверены, что после трансформации получится стоящий результат?

Если у вас исходный растр "в lat lon", то он должен выглядеть вот так, и на нем 180-й меридиан не проходит под углом:

Если он так не выглядит, то у вас проблема не с перепроецированием.