ArcMap. Draft NAC georeferencing.

Описание: на сайте lroc.asu.edu
Archive - Data Products (EDRs, CDRs, RDRs) - Image Search
Выделяем NACL, NACR,ShowThumbnails(Yes)
задаем координаты на нужный район Coordinate Range или ищем конкретный снимок например M165645602:
жмем Search
кликаем на нужный снимок
например на
M165645602RC

Cлева в таблице нам нужны координаты (они внизу таблицы)
Upper right latitude26.05
Upper right longitude30.79
Lower right latitude25.57
Lower right longitude30.79
Lower left latitude25.57
Lower left longitude30.72
Upper left latitude26.05
Upper left longitude30.71
(здесь надо еще отследить вывернут ли снимок наизнанку или развернут на 180 градусов)

Может понадобиться параметры солнца
Sub solar azimuth190.27
Sub solar latitude-0.91
Sub solar longitude323.36

перекопируем их в notepad



Выкачиваем слева вверху снимок Download CDR PTIF (multi-resolution, lossy
compression) или CDR IMG-файл (-он больше но без искажений),

далее в ArcCatalog в Properties для скачанного файла
1) делаем пирамиду (Build Pyramid)
2)Задаем географическую проекцию (в списке SolarSystem-Earth-Moon)
GCS_Moon_2000
Angular Unit: Degree (0,017453292519943295)
Prime Meridian: Reference_Meridian (0,000000000000000000)
Datum: D_Moon_2000
Spheroid: Moon_2000_IAU_IAG
Semimajor Axis: 1737400,000000000000000000
Semiminor Axis: 1737400,000000000000000000
Inverse Flattening: 0,000000000000000000
3)Открываем ArcMap, включаем toolbar Georeferencing
4)Загружаем полученный IMG файл (проект автоматом принимает координатную систему 1го добавленного слоя)
5)В тулбаре ищем кнопку XY (Goto to XY)
возникаем подменю - копируем из notepad первую пару координат, нажимаем на лупу+, если все правильно жмем в этом подменю на синюю точку с XY - появляется вспомогательная точка,
6) Повторяем для оставшихся пар
7)В тулбаре Georeferencing смотрим что высвечен слой M165645602RC, начинаем для кажого угла процедуру привязки (кликаем на соседнюю кнопку, кликаем на угол, в тулбаре Goto to XY отыскиваем нужную точку или самостоятельно вбиваем координаты, Zoom to, максимально увеличиваем, кликаем 2-й раз - угол переехал. Тоже с другими углами.
8) В тулбаре Georeferencing жмем Update Georeferencing. Снимок привязан, закрываем проект.
10) Чтобы померить метры, км, нужно создать проекцию Projected Coordinate System новую проекцию на основе PlateeCaree и геоида Луны.
11) Открываем ArcMap чистый бланк, в свойствах Layers задаем Coordinate System, жмем на глобус Add Coordianate System-New PROJECTED System - Projection Name2: PlateCaree, Change: выбираем Moon, Name1 (название созданной проекции): Bla-Bla-Bla
12) Добавляем в получившийся проект слой M165645602RC и можем мерить километры, а в обющем меню ArcMap подправляем Customize-ArcMap Options-DataView-Use these units instead - Decimal Degrees-OK (внизу справа в углу, теперь отображаются текущие градусы вместо метров).