Где взять карту?

[wowka1319]

Возникла проблема.
В первом сообщении в этой теме я приложил картинку интересующей меня области земной поверхности. Теперь у меня есть Geotiff этой области, который осталось перепроецировать в WGS84. Но я не знаю в какую проекцию, на приложенной картинке это явно не меркатор. Для справки картинка - это скрин из ArcMap. В этой программе я смог найти что это WGS84 (отсюда и мои первоначальные требования к перепроецированию), но вот параметры прямоугольной проекции я не нашел. Может кто знает как это узнать либо по виду скрина сможет определить параметры проекции?

[ericsson]

Это "аборт по фотографии".
Возьмите, почитайте вот это http://resources.arcgis.com/en/help/mai ... 000v000000 там должно быть описано, как посмотреть проекцию слоя.

[wowka1319]

Спасибо. Нашел: WGS_1984_UTM_Zone_54N. Выставляю это в QGIS и ошибка (ошибка программы - исключение, короче косяк разрабов).
Пробовал утилитой gdalwarp - выдает какую-то ошибку, но после нажатия на "продолжить" делает что нужно. НО качество УЖАСАЮЩЕЕ!!!
Что посоветуете?

[ericsson]

Опять аборт по фотографии.
Во-первых, куда вы это вставляете в QGIS? Там должен быть список и место для вставки строки параметров, а не названия.
Во-вторых, "какую-то ошибку" - это не разговор.
В-третьих, "качество УЖАСАЮЩЕЕ!!!" - это вообще одни эмоции.
Еще раз, конкретно, что именно делали (командные строки и т.п.), в чем проблема (текст ошибок, скриншоты..).

Для общего образования: http://spatialreference.org/ref/epsg/32654/ - ваша проекция. Ищите в QGIS по коду EPSG:32654 или вставляйте параметры отсюда http://spatialreference.org/ref/epsg/32654/proj4/

Еще я вижу одну проблему. Ваше покрытие, изображенное на картинке в первом сообщении, выходит далеко за пределы 54й шестиградусной зоны. Это потенциальный источник проблем с данными, выходящими за допустимую область.

[wowka1319]

Насчет корректности выставления проекции можете не волноваться. И я знаю, что EPSG:32654 - моя СК - ее я и указываю. Насчет ошибки gdalwarp - см. на картинке.

+ выхлоп из gdalwarp с УЖАСАЮЩИМ качеством. И это не эмоции, на мой взгляд полное отсутствие интерполяции как таковой в процессе перепроецирования (не покажу "до/после" ибо долго грузить, просто поверь наслово).

Насчет шестиградусных зон UTM я вкурсе, но такова установка в ArcMap на карту, на которую мне нужно накладывать снимки. Конечно показанная карта далеко не шестиградусный лоскуток, но тем не менее 54 зону все-таки в себя включает (о. Сахалин). Кстати на сколько большой становится погрешность "зональных" проекций за ее границами?

[ericsson]

Вот "отсутствие интерполяции" - это конкретный ответ, а "ужасающее качество" - это пустые эмоции. Мне ни к чему "верить на слово" - я и так вижу, что результат вас не устраивает. Я спрашивал о деталях не чтобы поужасаться, а чтобы понять, что именно не так. От эмоций в этом деле толку, очевидно, ноль. Теперь (еще и из командной строки) вижу.

Вы http://www.gdal.org/gdalwarp.html читали? Начиная с "-r resampling_method". Словечко resampling ничего не напоминает?
Еще, учитывая форму, к которой придет в итоге картинка, я бы порекомендовал подумать над "-dstalpha" или "-dstnodata value [value...]", но не обязательно, это по вкусу.

Заметьте, это все есть в документации.

Про погрешность - при таком охвате на краях, которые еще через зону находятся (на глаз), там будут искажения формы на километры. Это неплохо видно по кривой форме краев. В таких случаях обычно используют конические проекции совершенно не зря.