Как лечиь ASTER DEM

[Григорий]

При использовании ASTER DEM столкнулся с артефактами (значения сильно отличаются от реальных).Размер артефактов менее 10 пикселей. Есть ли скрипты, опции, расширения, которые позволяют редактировать растр (лучше grid) попиксельно?

[Boris]

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

[Григорий]

выгружать ничего никуда не надо. Потом привязывать заново устанешь. Можно это дело в ENVI легко организовать. Вопрос теперь в другом. Где взять правильные значения пикселов, т.е. реальные высоты вместо заоблачных значений? Как Вы решали этот вопрос ?

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

Григорий, редактировать растры (в т.ч., и попиксельно, вы можете почти в любой растровой ГИС -- GRASS, SAGA, ILWIS и т.д.) Но сомнительно, что таким редактированием вы сильно улучшите исходные данные. Хотя локально подровнять, конечно, можно. ЦМР вам нужны для работы/анализа или чисто для отображения рельефа?

"заоблачные значения" -- это просто ошибки ?.. вы бы примеры привели с картинками

[Григорий]

Мне нужны как можно более точные значения высот, причем для каждого пикселя.
Я использую это для анализа. Из ЦМР я собираюсь извлечь максимальные и минимальные высоты ледников (они красной границей показаны). Пример выглядит так:

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

[КО]: При расчётах имейте в виду заявленную вертикальную точность этих ЦМР -- около 20 м (по личным наблюдениям 10-15 м в горных районах).

Почитал файл README.pdf, который идёт с каждым архивом ASTER GDEM.
В растровом QA-файле (ASTGMT_...._num.tif) даны либо отрицательные значения (источник данных, использованный для замены ошибочных значений в dem), либо положительные (число исходных снимков для создания конечной dem).
Можете привести значения "сбойных" высот и соответствующие значения из QA-растра?

Как правило (по крайней мере, на терр. Евразии до 60 град. с.ш.), для замены "плохих" пикселей исп-ся SRTM v.3.
Видимо, в вашем случае значения не заменились как надо. Тогда это можно сделать самостоятельно (скачать SRTM v.3 на этот район; сделать растровую "маску" или ещё как-то убрать аномальные значения исходной ЦМР; пропатчить её значениями из SRTM).

[Григорий]

Какое пространственное разрешение у SRTM? У меня с моими данными на участок с артефактом пришелся 1 (один!) пиксел. Сбойные высоты в ASTER DEM составляют 55 с лишним тысяч метров или определяются как no data.

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

Какое пространственное разрешение у SRTM?

У SRTM 3'' горизонт. разрешение = 90м, но иногда (для определённых целей) можно переинтерполировать на более высокое разрешение (вплоть до 30м). Примеры статей, где это используется, я видел в сети.

+ SRTM более выверен, там меньше ошибок, чем в сыром пока ещё ASTER GDEM, поэтому для высотных характеристик чего-либо можно вполне уверенно применять.
А если объекты меньше 90м, то, наверное, есть другие пути уточнения (топокарты, ЦМР по космоснимкам, GPS, ..??)

[Григорий]

К сожалению, на мою территорию исследований SRTM3 нет, слишком далеко на севере. Придется цифровать топокарту:(

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

Может, просто тем или иным способом удалить ошибочные пиксели? А оставшиеся "дырки" залечить интерполяцией..