Коррекция ЦМР и удаление шумов

[kazakov.gis]

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

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

Вы определитесь для начала: это будет обработка, чистка или коррекция цифровой модели рельефа? Подходы к решению задач и инструменты будут различаться. Что именно вам нужно?
Стандартного "скользящего окна" обычно мало для подобных задач. Например, есть хороший инструмент для фильтрации ЦМР - mdenoise.
Ну и, как обычно, желателен пример данных.

[kazakov.gis]

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

Здравствуйте. Мне бы хотелось в целом узнать как можно улучшить ЦМР, какая последовательность в обработке и какие инструменты для этого существуют. Фильтры meshdenosing и другие типо (фокальной статистики и sample filter) я уже применил получается интересно. Конечная цель состоит в создание корректной цмр для дальнейшего морфометрического анализа.

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

Какая ЦМР является корректной для морфометрического анализа в вашем понимании? Какой именно морфометрический анализ? Чем не устраивают уже имеющиеся результаты?

[kazakov.gis]

Какая ЦМР является корректной для морфометрического анализа в вашем понимании? Какой именно морфометрический анализ? Чем не устраивают уже имеющиеся результаты?

Результаты сглаживания и фильтрации меня в принципе устраивают, я не могу понять как корректно удалить локальные ошибки (артефакты). Из выше мной прикрепленных рисунков на восточном склоне особенно явно виден один из примеров таких ошибок. Да и в целом хотелось бы знать какими еще инструментами можно улучшить ЦМР.

[andreygeo]

Результаты сглаживания и фильтрации меня в принципе устраивают, я не могу понять как корректно удалить локальные ошибки (артефакты). Из выше мной прикрепленных рисунков на восточном склоне особенно явно виден один из примеров таких ошибок. Да и в целом хотелось бы знать какими еще инструментами можно улучшить ЦМР.

Ошибки - это такие пропилы вертикальные? Это точно не что-то реальное? Мало ли, кто их рыб поймёт
А в каком у вас изначальном формате ЦМР? В принципе, можно попробовать фильтровать инструментами для обработки LAS. Тот же бесплатный CloudCompare относительно неплохо справляется с выбросами.

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

В принципе, можно попробовать фильтровать инструментами для обработки LAS. Тот же бесплатный CloudCompare относительно неплохо справляется с выбросами.

Если зашла речь о LAS, то можно было бы попробовать утилиты LAStools, там есть в том числе фильтрация.

[kazakov.gis]

В принципе, можно попробовать фильтровать инструментами для обработки LAS. Тот же бесплатный CloudCompare относительно неплохо справляется с выбросами.

Если зашла речь о LAS, то можно было бы попробовать утилиты LAStools, там есть в том числе фильтрация.

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

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

Первое, что приходит в голову — сделать растровую маску по такому выделяющемуся "объекту", вырезать по маске кусок из ЦМР, обработать теми же фильтрами, а потом сшить обратно с ЦМР.

[kazakov.gis]

Первое, что приходит в голову — сделать растровую маску по такому выделяющемуся "объекту", вырезать по маске кусок из ЦМР, обработать теми же фильтрами, а потом сшить обратно с ЦМР.

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

[pendduduk]

Я бы для начала попробовал пройтись медианой с большим окном (скажем 40х40), потом маленьким (4х4), вычесть одно из другого и используя пороговое значение отбить эти ямки и бугорки. Затем уже вместо бугорков и ямок вставить сглаженные, например гауссианом, значения из растра с большой медианой.

Все это легко реализовать на питоне в связке с numpy и scipy, примерно так

Код: Выделить всё

import scipy.ndimage as ndimage

med40 = ndimage.filters.median_filter(data,(40,40))
med4 = ndimage.filters.median_filter(data,(4,4))
dif = med40[:]-med4[:]
lo_mask = dif<t_low
hi_mask = dif>t_low
out = np.copy(data)
gauss = ndimage.filters.gaussian_filter(med40,25)
out[lo_mask]=gauss[lo_mask]
out[hi_mask]=gauss[hi_mask]

положите картинку с буграми попробую их убрать.

[kazakov.gis]

[quote="pendduduk"]Я бы для начала попробовал пройтись медианой с большим окном (скажем 40х40), потом маленьким (4х4), вычесть одно из другого и используя пороговое значение отбить эти ямки и бугорки. Затем уже вместо бугорков и ямок вставить сглаженные, например гауссианом, значения из растра с большой медианой.

К сожалению я не силен в питоне. Подскажите а какой инструмент отвечает за данный фильтр (фокальная статистика?). Так же я не могу понять как корректно проинтерполировать две перекрывающиеся поверхности чтобы не было вот такого эффекта???

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

я не могу понять как корректно проинтерполировать две перекрывающиеся поверхности чтобы не было вот такого эффекта???

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

[kazakov.gis]

Встречный вопрос: а как эти поверхности стали перекрывающимися? Если делать через растровые маски, этого не должно происходить.[/quote]
У меня есть два грида с небольшим перекрытием разных лет съемки на разные участки. Я их просто взял и проинтерполировал и в следствии чего получил такую бяку. Вот мне бы и хотелось узнать что делать чтобы такого не было. Растровые маски говорите.. Можно поподробнее????

[pendduduk]

К сожалению я не силен в питоне. Подскажите а какой инструмент отвечает за данный фильтр (фокальная статистика?). Так же я не могу понять как корректно проинтерполировать две перекрывающиеся поверхности чтобы не было вот такого эффекта???

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