геопривязка данных Modis

[jj_kan]

Уважаемый Максим,
я так понимаю, что форум, в том числе создан и для того, что бы можно было решать какие то сложные, может быть не для Вас, как давно работающего в этой области, но сложные для тех, кто впервые с этим сталкивается, вопросы.

Я и пытался выяснить, что можно сделать в этом случае.

Все равно спасибо, за полезную информацию.
Многое для себя прояснил.

[Александр Мурый]

Уважаемый amuriy извините все что Вы показали это пишется в командной строке?
Т.Е. выделяю файл, далее ввожу тот скрипт, что Вы прислали и потом для перевода в зону пулково последнюю строчку, правильно?
Еще вопрос для этих целей (или вместо Workstation) в ArcGis 9.3 был вставлен модуль Command Line, а в ArcGis 10 эту роль играет Piton?

Для того, что я написал, ArcGIS не нужен вообще. Поставьте GDAL (см. статью здесь на сайте) --> откройте командную строку Windows --> перейдите в папку с данными --> введите команду (та, что в одну строку, после слов "надо выполнить") --> в папке появится перепроецированный tif-растр в СК "Pulkovo42 / GK_нужная_зона"

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

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

Вы все верно понимаете. Просто попытайтесь четко отвечать на наводящие вопросы и не скачите с факта на факт.

[jj_kan]

amuriy еще раз огромное спасибо!

Максим, постараюсь придерживаться Ваших рекомендаций.


Еще раз спасибо всем.

[gamm]

Судя по косвенным признакам, товарищу передали стандартный продукт MODIS (типа MCD12Q1.005), который зачем-то перевели из исходного формата в TIFF и лишили проекции (обычно такое наблюдается при имитации бурной деятельности). Уж больно вид малой Родины знакомый ... ему нужно найти оригинал, и cпокойно его загрузить через MRT

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

не, не стандартный, результат классификации в рамках проекта Nelda
да, проекции лишили так или иначе
так что первым делом надо пробовать SIN/Sphere

[gamm]

н
так что первым делом надо пробовать SIN/Sphere

The tiles are in a MODIS sinusoidal projection with a reference Earth radius of 6371007.181 meters, a central meridian at 0 meters, and no false easting or northing. The pixel size is exactly 463.312716525 meters.
www.fsl.orst.edu/nelda/nelc/docs/nelc_d ... ion_v2.pdf

[jj_kan]

Максим и gamm все нормально, перепроецирование произведено.
Действительно, у TIFF не было проекции.
Пришлось задействовать несколько источников описания.
Сейчас разбираюсь с точностью наложения Modis на Landsat. Т.к. разрешение разное то есть смещение.
Интересно, есть ли материал какой нибудь по поводу этого вопроса.

Заранее благодарен.

[gamm]

Сейчас разбираюсь с точностью наложения Modis на Landsat. Т.к. разрешение разное то есть смещение.
Интересно, есть ли материал какой нибудь по поводу этого вопроса.

если внутри программы имеет место округление размера пикселя до сантиметров, то будет набегать ошибка порядка 10км (или используется тип float, в этом случае ошибка будет меньше, 1-3 пикселя Модис). Если причина в этом, то ошибка по горизонтали будет примерно в 2 раза меньше, чем по вертикали, и Модис "съедет" влево-вниз относительно Ландсата.

[jj_kan]

Да gamm Вы правы, большая по величине ошибка именно по вертикали, но как бы видно что MOdis "съехал" не влево-вниз, а наоборот вправо-вверх.
А это можно исправить, что бы не по контрольным точкам, а более точным методом?

Еще раз спасибо.

[gamm]

Да gamm Вы правы, большая по величине ошибка именно по вертикали, но как бы видно что MOdis "съехал" не влево-вниз, а наоборот вправо-вверх.
А это можно исправить, что бы не по контрольным точкам, а более точным методом?

щас попробовали загрузить в Сканекс IP (с указанием проекции на сфероиде), и положили сверху вектор Волги (1:200,000) - все вроде легло без смещения. А у вас какие смещение ( в метрах)?

[jj_kan]

около 700 м

[gamm]

около 700 м

тогда нужно смотреть, по 200-ке понять такие сдвиги сложно, они сопоставимы с размером пикселя Модис. Если выложите небольшой кусок Ландсата в Geotiff (внутри границ Модис), могу попробовать наложить его на Модис, и послать картинку.

[jj_kan]

конечно gamm посылаю.
если можно, поподробнее что Вы делаете, что бы исправить эту погрешность.
Заранее благодарен.

[gamm]

если можно, поподробнее что Вы делаете, что бы исправить эту погрешность.

я не уверен, что там есть погрешность - возможно, погрешность классификации.

картинка с наложенным фрагментом береговой линии прилагается

picture.jpg (545.77 КБ) 11574 просмотра