Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
https://gis-lab.info/forum/
Здравствуйте. Мне бы хотелось в целом узнать как можно улучшить ЦМР, какая последовательность в обработке и какие инструменты для этого существуют. Фильтры meshdenosing и другие типо (фокальной статистики и sample filter) я уже применил получается интересно. Конечная цель состоит в создание корректной цмр для дальнейшего морфометрического анализа.Александр Мурый писал(а):Вы определитесь для начала: это будет обработка, чистка или коррекция цифровой модели рельефа? Подходы к решению задач и инструменты будут различаться. .
Результаты сглаживания и фильтрации меня в принципе устраивают, я не могу понять как корректно удалить локальные ошибки (артефакты). Из выше мной прикрепленных рисунков на восточном склоне особенно явно виден один из примеров таких ошибок. Да и в целом хотелось бы знать какими еще инструментами можно улучшить ЦМР.Александр Мурый писал(а):Какая ЦМР является корректной для морфометрического анализа в вашем понимании? Какой именно морфометрический анализ? Чем не устраивают уже имеющиеся результаты?
Ошибки - это такие пропилы вертикальные? Это точно не что-то реальное? Мало ли, кто их рыб поймётkazakov.gis писал(а): Результаты сглаживания и фильтрации меня в принципе устраивают, я не могу понять как корректно удалить локальные ошибки (артефакты). Из выше мной прикрепленных рисунков на восточном склоне особенно явно виден один из примеров таких ошибок. Да и в целом хотелось бы знать какими еще инструментами можно улучшить ЦМР.
Если зашла речь о LAS, то можно было бы попробовать утилиты LAStools, там есть в том числе фильтрация.andreygeo писал(а):В принципе, можно попробовать фильтровать инструментами для обработки LAS. Тот же бесплатный CloudCompare относительно неплохо справляется с выбросами.
Спасибо вам за интересный софт, будет время обязательно поюзаю. А подскажите в ручном режиме можно как-то отредактировать ошибки вот такого рода.Александр Мурый писал(а):Если зашла речь о LAS, то можно было бы попробовать утилиты LAStools, там есть в том числе фильтрация.andreygeo писал(а):В принципе, можно попробовать фильтровать инструментами для обработки LAS. Тот же бесплатный CloudCompare относительно неплохо справляется с выбросами.
Похоже только это и остается. А вот еще вопрос каким инструментом корректно слить два грида с перекрытием в один, чтобы в области перекрытия посчитались новые значения и по ним постройлся новый гридАлександр Мурый писал(а):Первое, что приходит в голову — сделать растровую маску по такому выделяющемуся "объекту", вырезать по маске кусок из ЦМР, обработать теми же фильтрами, а потом сшить обратно с ЦМР.
Код: Выделить всё
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 писал(а): я не могу понять как корректно проинтерполировать две перекрывающиеся поверхности чтобы не было вот такого эффекта???
По первому вопросу, не совсем понял... если вы имеете в виду где искать в меню медианный фильтр, то не подскажу т.к. не представляю с каким ПО вы работаете.kazakov.gis писал(а):К сожалению я не силен в питоне. Подскажите а какой инструмент отвечает за данный фильтр (фокальная статистика?). Так же я не могу понять как корректно проинтерполировать две перекрывающиеся поверхности чтобы не было вот такого эффекта???