Пара вопросов от дилетанта про ЛАНДСАТ

[AN1440]

Вот: получил по заказу ландсатовский снимок (интересовали области, погоревшие в процессе лесного пожара).

НО:

Во-первых снимок запроектирован не в ту проекцию, которой я обычно предпочитаю пользоваться. (УТМ-83 вместо желаемого "Альбертса"). То есть, работать я с им смогу, но будет несколько неудобно. Работаю я в АркВъю-330, и другого пока не предполагаю.

Во-вторых - запроектировано как-то не очень... Ошибка позиционирования явно доходит до 330 метров. Вездеходная дорога пошла прямо по фарватеру реки, а кордон оказался на противоположном берегу. Несмертельно, но неприятно. Есть ли возможность это как-то исправить?

В третьих - хотелось бы знать, когда это сфотографировано. Снимок называется "09052002" это есть это что - 9 мая или 5 сентября?

Ну и, в-четвёртых - "112_026" - по этому обозначению, надо полагать, можно определить местоположение снимка на земной поверхности. Как это делается? Ну, для этого-то конкретного снимка местоположение я знаю, а - в общем случае?

Понимаю, вопросы дурацкие, но ответы хотелось бы найти. И желательно, не на английском языке...

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

Если бы вы сообщили полное название снимка быо бы проще.

Ошибка позиционирования явно доходит до 330 метров. Вездеходная дорога пошла прямо по фарватеру реки, а кордон оказался на противоположном берегу. Несмертельно, но неприятно. Есть ли возможность это как-то исправить?

Допривязать - Arcview Image Analysis -Align, например или ERDAS и т.д.

В третьих - хотелось бы знать, когда это сфотографировано. Снимок называется "09052002" это есть это что - 9 мая или 5 сентября?

Если нельзя определить по косвенным признакам, например фенологии, то нужно смотреть по источнику снимка, а вообще назвать файл можно как угодно.

Ну и, в-четвёртых - "112_026" - по этому обозначению, надо полагать, можно определить местоположение снимка на земной поверхности. Как это делается? Ну, для этого-то конкретного снимка местоположение я знаю, а - в общем случае?

Это ряд и колонка в системе WRS, подробнее
http://gis-lab.info/qa/l7pathrowmap.html

[Clever]

В третьих - хотелось бы знать, когда это сфотографировано. Снимок называется "09052002" это есть это что - 9 мая или 5 сентября?

дд/мм/гггг обычно ландсаты так обозначают
ну или гггг/мм/дд

[geologic]

в американском формате месяц идет вперед. Например, снимок 1999 года на вашу территорию называется p112r026_7t19990922.742.preview.jpg, что, разумеется, обозначает 22-е сентября, а не наоборот

бесплатный Ландсат можно найти вот тут
ftp://ftp.glcf.umiacs.umd.edu/glcf/Land ... -111/WRS2/, дальше все разложено по Path/Row

Возможно ваш снимок просто "съехал" на 300 метров, у нас такое было с космоснимками не раз. Дело в том, что поставщики снимков считают свою задачу выполненной после орторектификации и установления проекции, их аргументация такова: "параметры сдвига вам легче будет 'лапой' подобрать" по местности" В этом случае не обязательно править изображение сильнодействующими средствами, достаточно просто сдвинуть его обратно, например, как описано тут для ArcView
http://geofaq.narod.ru/2wld.htm

[AN1440]

А что - эти снимки "ландсатовцы" поставляют уже готовыми, привязанными? Или "пере"поставщики перепривязывают самостоятельно?
Можно ли его перепривязать и перепроектировать?
-------------
{легче будет 'лапой' подобрать" по местности" } - Это реально можно сделать?

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

Дело в том, что поставщики снимков считают свою задачу выполненной после орторектификации и установления проекции, их аргументация такова: "параметры сдвига вам легче будет 'лапой' подобрать" по местности"

Странно сформулировано, не согласен.

Cсылка на полный архив GLCF:
ftp://ftp.glcf.umiacs.umd.edu/glcf/Landsat/

А что - эти снимки "ландсатовцы" поставляют уже готовыми, привязанными? Или "пере"поставщики перепривязывают самостоятельно?
Можно ли его перепривязать и перепроектировать?

Да, нет, да (GLCF = "перепоставщики")

[geologic]

каюсь, грешен... Но зато есть повод разъяснить смысл, который от формулировки мало зависит. В качестве поставщиков я не имел в виду сайт GLCF - это скорее публикация снимков... Поставщики на нашей территории это Сканекс, например и другие. Можно соглашаться или нет, но именно так мы снимки и получали от весьма уважаемых поставщиков, с упоминанием, что это ПРАВИЛО, а не исключение, что они ВСЕГДА так делают и могут для нас постараться, но точек мол нет на территорию, уточнение - отдельная работа, да и зачем вам проблемы с секретностью и т.п. Возмущаться же тут не было никакого смысла - работа по привязке по точкам, (точнее, по уточнению привязки - это все же ее последняя стадия, а не единственная) если разобраться, действительно трудоемкая если делать ее с гарантией точности. Тогда как посадить по одной точке и проверить по другой плоский снимок самим несложно. Имея базу, конечно.

Первый такой "двинутый" ландсат нас удивил, базы у нас под него не было, но выкрутились... двигали "лапой" прямо в процессе работы кому как вздумается, порождая творческие варианты Появилась база - снимок сел на место опять же простым сдвигом.

Предваряя другие возражения сразу скажу больше - высокоточные снимки также не могут служить сами по себе координатной базой, даже если поставщик уверяет, что они точно посажены. Последняя стадия всегда проверка/коррекция по ВАШЕЙ базе. Яркий пример - снимки гугл на Москву, при разрешении 1 м ошибки посадки 20-25 м. Но снимки орторектифицированы, проекция (теперь известна, коррекция - дело техники. Вот что имелось в виду.

2AN1440: Упомянутые снимки "ландсатовцы" поставляют привязанными, то есть со встроенной гис-геопривязкой. Мало того, для каждого PATH/ROW известен фиксированный "квадрат" привязки (можно скачать шейп-файл), и в особом файле метаданных документировано описание всех угловых точек, проекции, разрешения и т.п. - на случай, если ваш софт не понимает геотифф, например. Однако это не значит, что на реальной местности вершина горы совпадет с координатами тригопункта, например - некая ошибка все равно имеется (не 300 метров, конечно, но есть).

Кстати, там же в метафайле есть и точное время проведения съемки, потому что, как справедливо заметил Sim, имя файла еще ничего само по себе не гарантирует - только подсказка.

[AN1440]

Спасибо большое!
Только я понял так: "Привязка" снимка заложена в самом файле снимка. Ведь к ему не прикадывается никакой дополнительный файл.
То есть, для правки надо поменять с десяток байтов файла. Это реально сделать вручную? Надо знать, где нужные байты находятся... А менять байты в файлах - это мое любимое занятие...
-------------------------------------------
Кстати: много раз пользовался программой "ИмагеВорп", которая привязывает растры для АркъВъю. Странно - она почему-то округляет полученные результаты до величин, гораздо более грубых, чем это самое АркъВъю может воспринять.
Полученный world-файл приходится всегда "доправлять" ручками до нужной точности.

[geologic]

Привязка в геотиф действительно "спрятана внутри", но и в малобюджетном варианте есть способы ее оттуда выковырить. Как - описано, например, здесь: http://geofaq.narod.ru/2spatial_ps.htm#8
Обратно ее впихивать необязательно - можно использовать world-файл. Кстати, если вы читали вдумчиво, то наверно сообразили, что world-файл сделать легко по данным meta-описания (если его у вас нет, подойдет наверно любое для вашего path/row). Еще совет - из-за одного файла даже и так мучиться не надо, загрузите файл в какой-нибудь бесплатный гис-вьювер, например ECW-reader, там можно всю геопривязку вычитать из Info.

Мое мнение об ImageWarp изложено по первой моей ссылке - суета эти все полиномы, поверьте. Ограничьте преобразования растра сдвигом/вращением/масштабированьем, и будет вам щастя Для проекций же и выпучиваний лучше применять продвинутый софт.