Форумы GIS-Lab.info

rhot писал(а):Не могли бы вы выложить данные в более свободном формате, чем детище Рошала?

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

АлексЮстасу писал(а): Чтобы можно было оценить и сравнить результаты работы разных программ, разных алгоритмов.

АлексЮстасу писал(а): ЦМР нерегулярная, получена из данных воздушного лазерного сканирования.

АлексЮстасу писал(а): На выходе - эта же ЦМР (в смысле - то же число точек, положение в плане точек может изменяться, но только за счет поправок за шумы, если они вдруг вносятся во все три координаты, а не только в Z), но с удаленными шумами.

rhot писал(а):Вот например. Первые 3 из списка вполне подойдут.

Александр Мурый писал(а):Задумка хорошая, как раз для статьи на ГИС-Лаб

Александр Мурый писал(а):Вы хотите сравнивать алгоритмы или ПО (это разные задачи)?

Александр Мурый писал(а):

АлексЮстасу писал(а): ЦМР нерегулярная, получена из данных воздушного лазерного сканирования.

Т.е. это "облако точек", а не ЦМР.Прежде чем очищать облако точек, не надо ли отделить собственно рельеф от объектов на нём? Или здесь это уже сделано? Для выделения из данных лидарной съёмки собственно ЦМР есть спец. алгоритмы, в т.ч., в открытых ГИС (GRASS, SAGA). Кроме этого, есть спец. утилиты (например, libLAS, оно же LAStools.

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

library(MBA)

p<-read.table("place_left_3.xyz",header=FALSE)
names(p)<-c("X","Y","Z")

plot(p[,1:2],pch=20,cex=0.3,asp=1)

tmp_011_08<-mba.surf(p,100,100,n=1,m=1,h=6,extend=TRUE)

x.pos<-tmp_011_08$xyz.est$x
y.pos<-tmp_011_08$xyz.est$y
z.mtr<-tmp_011_08$xyz.est$z
z.lim<-range(as.double(z.mtr))

image(x.pos,y.pos,z.mtr,col=gray((126:32)/128),asp=1)
contour(x.pos,y.pos,z.mtr,levels=seq(floor(z.lim[1]),ceiling(z.lim[2]),by=0.5),col=" black",add=TRUE,asp=1)

gamm писал(а): картинку приложить не удалось - и bmp, и emf запрещено.

gamm писал(а):начнем с того, что у вас нет ЦМР

gamm писал(а):для вашего случая придумано два подхода: традиционный - цифровая фильтрация (например, на основе локального Фурье и прочих Кальманов), и новый - быстрая аппроксимаций В-сплайнами MBA - multilevel B-spline approximation, она есть в R, пакет MBA. Последняя хороша тем, что работает с линейной скоростью, так что не боится миллионов точек.

gamm писал(а):смастерил за пару минут ...

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

library(MBA)

p<-read.table("place_left_3.xyz",header=FALSE)
names(p)<-c("X","Y","Z")

plot(p[,1:2],pch=20,cex=0.3,asp=1)

tmp_011_08<-mba.surf(p,100,100,n=1,m=1,h=6,extend=TRUE)

x.pos<-tmp_011_08$xyz.est$x
y.pos<-tmp_011_08$xyz.est$y
z.mtr<-tmp_011_08$xyz.est$z
z.lim<-range(as.double(z.mtr))

image(x.pos,y.pos,z.mtr,col=gray((126:32)/128),asp=1)
contour(x.pos,y.pos,z.mtr,levels=seq(floor(z.lim[1]),ceiling(z.lim[2]),by=0.5),col=" black",add=TRUE,asp=1)

gamm писал(а):картинку приложить не удалось - и bmp, и emf запрещено.

Что такое "R"?

Александр Мурый писал(а):А почему только bmp и emf? PNG ещё никто не отменял

Variant писал(а):пользуясь случаем, хочу упомянуть здесь нестандартный подход, показавшийся мне довольно любопытным. В отличие от обычного сглаживания помеха здесь не даёт практически никакого вклада в итоговый результат:
http://www.miigaik.ru/nauka/dissertacio ... 9-7732.doc
не знаю, правда, удастся ли каким-то образом учесть его в сравнении, разве что попытаться выйти на контакт с автором

АлексЮстасу писал(а):На нерегулярную ЦМР множество точек поверхности не тянет?

Что такое "R"?