Форумы GIS-Lab.info

Здравствуйте уважаемые форумчане. Не хочу плодить темы, спрошу тут, тем более вопрос по теме.

В данной теме я полный нуб, так что сильно не пинайте за корявые вопросы)

Задача: Создать 3Д рельеф для квадрата E65N072 (далее будет использоваться в симуляторе x-plane).
С данного сайта скачал ASTGTM2_N65E072_dem, утилитами скомпилировал меш. На выходе получилось http://pixs.ru/showimage/RRRRRRSRSR_5974113_9304803.png Видно что есть впадины и скалы чего не должно быть. С помощью GM открыл dem где и обнаружил эти самые огрехи. Но что дальше делать ума не приложу.
Подскажите, как и чем я могу отредактировать сетку. Пока задача стоит просто убрать скалы и впадины, особая точность рельефа не требуется. Спасибо.

Убрать такие неровности можно путем фильтрации/сглаживания данных различными методами (сплайн и т.п.), во всяком случае когда-то в Surfer так делал...

Vasilidze писал(а): Пока задача стоит просто убрать скалы и впадины, особая точность рельефа не требуется. Спасибо.

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

Vasilidze писал(а):Видно что есть впадины и скалы чего не должно быть.

Астэр - такой Астэр...
В ГИС есть фильтрация для сглаживания артефактов. В ArcGIS за это отвечает инструмент Focal Statistics с параметром neighborhood "circle" и радиусом в пикселях (обычно 3-5, но можно и больше сглаживать). В итоге получите - более сглаженный рельеф, если "ям" мало, то можна делать руками.
Советую использовать http://srtm.csi.cgiar.org/ который так же можно передескритизировать в 30м ячейку, если нужно.

dime1 писал(а): Советую использовать http://srtm.csi.cgiar.org/ который так же можно передескритизировать в 30м ячейку, если нужно.

Квадрат E65N072 лежит за пределами покрытия данными SRTM.

oops
Было бы интересно сделать инструмент, которому задаётся некий порог, ниже которого значения заполняются наибольшим значением из окружающих ячеек.
Допустим, реально на растре рельеф от 50 до 250м, мы инструменту задаём <=50 и заполняем эту яму наибольшим значением по контуру (в 1 пиксель шириной) ямы.

Легко "подчистить" ASTER GDEM до приличного состояния лично у меня не выходило пока ни разу. Тем более с районах с перепадом высот в несколько десятков метров, как в квадрате E65N072.
Один из вариантов — использовать альтернативные источники ЦМР, например, GMRT, ACE2.

dime1 писал(а): Было бы интересно сделать инструмент, которому задаётся некий порог, ниже которого значения заполняются наибольшим значением из окружающих ячеек.

В GRASS это можно сделать сочетанием модулей r.mapcalc и r.neighbors.

dime1 писал(а): Допустим, реально на растре рельеф от 50 до 250м, мы инструменту задаём <=50 и заполняем эту яму наибольшим значением по контуру (в 1 пиксель шириной) ямы.

А вот здесь, видимо, придётся попрограммировать. Мне вот быстро не приходит в голову, как сделать стандартными инструментами ГИС "наибольшее значение по контуру".

Александр Мурый писал(а):А вот здесь, видимо, придётся попрограммировать. Мне вот быстро не приходит в голову, как сделать стандартными инструментами ГИС "наибольшее значение по контуру".

Как и много другое. Просто высказал идею, интересно было послушать отзывы. А то она либо: уже кем-то реализована (в рамках "промышленной" подготовки ЦМР) или я что-то упустил, как с покрытием srtm)

Александр Мурый писал(а):А вот здесь, видимо, придётся попрограммировать. Мне вот быстро не приходит в голову, как сделать стандартными инструментами ГИС "наибольшее значение по контуру".

можно сделать и проще, и лучше:
1) создаем маску "плохих мест"
2) все данные "хороших мест" используем для интерполяции (нужен быстрый интерполятор, из известных мне годится только MBA, который есть, например, в R - или использовать TIN или RBF, если есть - будет хуже, но сгодится.). Интерполируем на ту же сетку.
3) впечатываем результат интерполяции в "плохие места". Поскольку симулятор, то точность не важна - добавляем дисперсии для красоты (или фракталами, или случайными полями, посчитав вариаограмму в окрестностях)

P.S. учитывая специфику GDEM, разумнее попробовать "уравнять" среднее значение "плохого места" и его "хорошего" окружения, рельеф там обычно правильный, но смещенный.

gamm писал(а): 1) создаем маску "плохих мест"

А вот с этого момента поподробнее, пож-та. Какие методы выделения "плохих мест" лучше подходят в случае ASTER GDEM?