Форумы GIS-Lab.info

qssaka писал(а): Обязательно покапаю GRASS, но.. я уже копался в Geostatistical Analist, но так и не смог найти там скользящее окно. Надеюсь здесь найдется кто-нибудь, кто уже сталкивался с подобным типом задач.

qssaka писал(а):На следующем шаге программа должна мочь/уметь построить линии изоконцентраций по усредненным данным. и да здесь тоже желателен выбор каким образом интерполировать (ближайший neighbor, кригинг, etc).

#load data
p<-read.table("my_data.dat",header=TRUE)
# X Y Val1 Val2
#12345.98 321345.89 0.9834
dx<-100 # 100m
dy<-100 # 100m
X0<-min(p$X)-0.5*dx
Y0<-min(p$Y)-0.5*dy
iX<-as.integer((p$X-X0)/dx)
iY<-as.integer((p$Y-Y0)/dx)
# координаты лучше тоже усреднить, чтобы уменьшить ошибку
tmp<-aggregate(p,by=list(iX=iX,iY=iY),FUN=mean)
#save data
write.table(tmp,"new_data.dat",row.names=FALSE,sep="\t")

qssaka писал(а):ума не приложу где в серфере может находится Moving Window

Surfer Help писал(а):Moving Average
The Moving Average gridding method assigns values to grid nodes by averaging the data within the grid node's search ellipse

qssaka писал(а):По разным причинам проводилось неравномерное опробование (не по сети). При этом во множестве точек опробования выявлены рудные концентрации свинца. Значения концентраций в таких точках являются выбросами. => Таким образом, при интерполяции вокруг "рудных" точек , стоящих особнячком, возникают чудовищно большие аномалии. Я интерполировал в Surfer и в ArcMap Geostsat. analyst методами nearest neighbor и kriging (ordinary и simple).

qssaka писал(а):to juffin_h (to gamm) Да у кригинга есть эллипс в который он подбирает точки. Сейчас я почти на 100% уверен, что с этой бедой (выбросы --> несоразмерно большие аномалии) можно бороться подставляя в центр этого окна (эллипса, etc) не средннее, а медиану. Тогда окно должно быть достаточно большим, чтобы захватывать минимум 5 - 10 точек (потому как медиана по 2ум точка = среднее, не так?) Но, в этом случае окно должно двигаться с перекрытием, а не "пошагово"!, иначе вся площадь может быть покрыта 40-50 такими окнами, а это равносильно 40-50 точкам для дальнейшей интерполяции. Теперь бы алгоритм движения окна продумать...

juffin_h писал(а):В Surfer, среди методов интерполяции (Gridding Methods), есть метод Data Metrics. По сути, это не метод интерполяции, а инструмент для анализа набора исходных данных. На основе этого анализа делается выбор оптимального метода интерполяции.