Пулкова 1942 в ENVI

[AndreyL]

Чего-то я не понимаю.

Возникла необходимость привязать 25-ку. Программа – ENVI 4.1. Ставлю проекцию GK Zone18 (Pulkova 1942), а какой ставить датум? В списке датумов нет Красовского, по дефолту ставится <none>.

Для проверки беру космоснимок (Landsat7, UTM 48N, WGS84), проектирую его в GK Zone18 (Pulkova 1942) датум <none>, рисую координатную сетку – визуально сетка карты не совпадает с сеткой на снимке метров на 200. Как быть?

[Максим Дубинин]

В установочной директории ENVI есть такой файл, его надо отредактировать, добавить новый datum.

D:\RSI\IDL61\products\ENVI41\map_proj\datum.txt

визуально сетка карты не совпадает с сеткой на снимке метров на 200

Похоже на то, что карта в СК Пулково 1942, а снимок в WGS 1984, сдивг порядка 200 м в этом случае нормален. С другой стороны, ошибка привязки ландсат - до 250 м по спецификации. И что такое "сетка на снимке"?

[AndreyL]

В файле datum.txt нашел «AFG, Krassovsky, -43, -163, 45». Это то, что нужно, или не совсем? И за что отвечают три числа после названия эллипсоида? Для WGS-84, например, все три числа нули. И как добавить новый datum в этот файл – какие числа ставить?

При попытке добавить новую проекцию, импортируя ее из «GK Zone18 (Pulkova 1942)» появляется следующее:
Эллипсоид Krassovsky (не датум, а эллипсоид)
Проекция поперечная Меркатора
Восточный и северный сдвиги, параллель и меридиан, масштаб – вроде бы все похоже.

Теперь пытаюсь пользоваться конвертором координат:
Исходная проекция UTM 48N WGS84, новая проекция GK Zone18 (Pulkova 1942) datum <none>. Координаты (широта, долгота) в исходной и в новой ОДИНАКОВЫЕ!!! Т.е. получается, что у этих проекций одинаковый датум, и datum <none> (по дефолту во всех GK (Pulkova 1942)) никакого отношения к эллипсоиду Крассовского не имеет. Или я не прав?

что такое "сетка на снимке"?

Эту сетку я рисую в ENVI после того, как снимок из UTM спроектирован в GK (Pulkova 1942), так что расхождения в проекциях быть не должно – карта и снимок находятся в одинаковой проекции, а именно GK Zone18 (Pulkova 1942)

ошибка привязки ландсат - до 250 м по спецификации.

В соответствии с данными из ftp://ftp.glcf.umiacs.umd.edu/glcf/Landsat/README/ :
Position Accuracy 50m (CI 90%), это для TM (bands 7, pixel size 28.5 м)
Для черно-белого (pixel size 14.25 м) не нашел, но есть подозрение, что цифра должна быть той же. К тому же точки GPS на снимок ложатся с точностью метров 30 (причем ошибка систематическая, но это тема для другого раздела форума).

Но где то я читал, что на советских картах координатная сетка может быть смещена, чтобы враги не догадались. Тогда на мой, крайне некомпетентный, взгляд решение проблемы:
а) задание правильного датума
б) коррекция сдвигов

Как это правильно сделать?

[Максим Дубинин]

И за что отвечают три числа после названия эллипсоида?

По поводу цифр - прочитайте статьи по настройке GPS на этом же сайте, там довольно много материалов на эту тему. Если коротко - то надо ставить такие цифры для той СК, в которой сделана ваша сетка. Попробуйте 2 набора цифр для начала - ГОСТовский и NIMA.

карта и снимок находятся в одинаковой проекции

Проекция может и одна, но этого недостаточно. Если бы это было так, то предыдущего вопроса не было бы. UTM и GK - это названия проекций, а не систем координат. Исходно Landsat в WGS84, а карты как я понял в Pulkovo 1942, проекция не так важна в данном случае, важно правильно осуществить переход от одной СК к другой.

Position Accuracy 50m (CI 90%)

Ссылка которую вы привели - это ссылка на производный продукт на базе Landsat - Geocover, по его спецификации действительно такая ошибка. Для исходного Landsat'a - 250 м.

[AndreyL]

По поводу цифр - прочитайте статьи по настройке GPS на этом же сайте, там довольно много материалов на эту тему. Если коротко - то надо ставить такие цифры для той СК, в которой сделана ваша сетка. Попробуйте 2 набора цифр для начала - ГОСТовский и NIMA.

Да! Только что это за цифры (в файле, три штуки на датум)? Если это смещения центра эллипсоида, то в каком они порядке? Если это не смещения центра эллипсоида, то что это? Есть ли эта информация в документации ENVI, если есть, то где? Я пока не нашел – наверное плохо искал.

[AndreyL]

Извиняюсь, нашел. userguid.pdf -> Appendix B: ENVI File Formats -> ENVI Datum File (datum.txt)