Сшивка DEM

[remotesens]

Всем доброго времени дня. Столкнулся с такой проблемой: при сшивке нескольких фрагментов DEM в одну мозаику (неважно - SRTM или ASTER GDEM) видимо в результате интерполяции растр сползает на один - пол- пикселя в зависимости от применяемого софта.
Когда мозаикируешь просто снимки это не суть важно, но если сшиваешь DEM - это серьезно. Если потом делать анализ по этому DEM-у или использовать его для подложки во время ортотрансформирования это может снижать качество результатов. Ведь абсолютная точность по горизонтали того же SRTM- 30 метров, а если учитывать, что размер пиксела 90 метров сдвиг на один пиксел это критически много.
Короче говоря, буду очень благодарен, если у кого есть софт или подход, чтобы сшивать DEM-ы без результирующего сдвига.

[andreygeo]

Всем доброго времени дня. Столкнулся с такой проблемой: при сшивке нескольких фрагментов DEM в одну мозаику (неважно - SRTM или ASTER GDEM) видимо в результате интерполяции растр сползает на один - пол- пикселя в зависимости от применяемого софта.

Дурацкий вопрос:..а как это проявляется? Сколько сшиваю в GM никогда не замечал, может быть, конечно, не обращал внимания.
А какие точности (масштаб) вы хотите вытянуть при ортотрансформировании с использование SRTM?

[remotesens]

Хороший метод ответа - сначала обозвать вопрос дурацким, а потом заявить, что не обращал внимания и поэтому не в курсе. Попробуйте в GM даже просто пересохранить растр - увидите сдвиг.
Если у кого есть ответ на конкретный вопрос, буду благодарен за ответ.

[andreygeo]

Хороший метод ответа - сначала обозвать вопрос дурацким, а потом заявить, что не обращал внимания и поэтому не в курсе. Попробуйте в GM даже просто пересохранить растр - увидите сдвиг.
Если у кого есть ответ на конкретный вопрос, буду благодарен за ответ.

Да я ж собсна свой вопрос и обозвал дурацким Как проявляется этот сдвиг? Вот такой дурацкий вопрос. В GM сшил, пересохранил растр...не могу понять, сдивнулся или нет. Туплюс...может принтскрин выложите?

[ericsson]

Для начала, вопрос: вы в курсе, что для разных форматов высотных данных указанные в метаданных координаты могут соответствовать в разных случаях разным местам (углам или центрам угловых пикселей)? Если да, то вы знаете точно, как должно быть в случае каждого из используемых вами форматов?

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

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

Как уже сказали, неплохо было бы увидеть скриншот этих самых "сползаний". Сшивал SRTM и ASTER GDEM много раз (с помощью gdal_merge, GRASS, QGIS), при этом не наблюдал подобных проблем.
Вы сшиваете растры в исходной СК, не перепроецированными?

[ericsson]

Сдвиг такого рода не является чем-то невероятным - собственно, у GDAL с этим проблемы были - в багтрекере это отражено, и у GlobalMapper - тоже, довольно давно. Вопрос, скорее, в том, каким образом контролируется результат, какие именно данные объединяются и чем.
Потому что пока не упомянуты ни конкретные данные, ни конкретный софт.

[remotesens]

Всем доброго времени дня. Столкнулся с такой проблемой: при сшивке нескольких фрагментов DEM в одну мозаику (неважно - SRTM или ASTER GDEM) видимо в результате интерполяции растр сползает на один - пол- пикселя в зависимости от применяемого софта.

Дурацкий вопрос:..а как это проявляется? Сколько сшиваю в GM никогда не замечал, может быть, конечно, не обращал внимания.
А какие точности (масштаб) вы хотите вытянуть при ортотрансформировании с использование SRTM?

А понял, насчет дурацкого вопроса. Извиняюсь тогда.
Выражается в сдвиге...
Global Mapper - полпикселя
ENVI - 1 пиксел
PCI Geomatica - полпикселя
INPHO Orthovista - полпикселя.
Между прочим сдвиг возникает даже просто при пересохранении в файла, например в Globalmapper, в других не знаю. Видимо это связано с особенностями расчета положения пикселов (решение наверное какое-нибудь вероятностное, но я тут некомпетентен - говорить не буду).
В архиве по ссылке картинки.
SRTM - до мозаикирования в GM, AFTER - после (никакого преобразования систем координат не выполнялось - все в родной географии). Взят краешек с водой, чтоб ярче было видно, растры без свертки, красная линия - краешек, очерченный по первичному SRTM-у.
Цель: понятно, что в общераспространенных ЦММ до фига погрешностей и кроме этой, но не хотелось бы еще искусственно ухудшать их при мозаикировании сдвигом на 45/90 метров.

[lam]

Экспорт в формат аски файла (*.asc) и правите в шапке координаты XY на величину смещения полпискеля, затем обратно импорт в нужную вам программу.

[remotesens]

Спасибо, попробуем. О результате сообщу.

[ericsson]

http://gis-lab.info/qa/asciigrid.html там есть полезное вам замечание про corner/center

[remotesens]

http://gis-lab.info/qa/asciigrid.html там есть полезное вам замечание про corner/center

Тоже спасибо за информацию. Будем изучать.