load_ext rpy2.ipython
The rpy2.ipython extension is already loaded. To reload it, use: %reload_ext rpy2.ipython
import numpy as np
from pandas import DataFrame, Series
import fiona
import matplotlib.pyplot as plt
c = fiona.open('/Users/martinlaloux/rs3D/rs3D.shp')
x,y = zip(*[i['geometry']['coordinates'] for i in c])
z = [i['properties'][u'elev'] for i in c]
x = np.array(x)
y = np.array(y)
z = np.array(z)
%Rpush x y z
%%R
library(gstat) # Loads gstat library
library(sp)
df = data.frame(x, y, z)
names(df)
[1] "x" "y" "z"
%%R
coordinates(df) = ~x+y
%%R
walk.var2 <- variogram(z ~ x+y,data=df,width=10,cutoff=534)
%%R
plot(walk.var2,cex=1.3,pch=16,col=1,xlab="Distance",ylab="Semivariance",main="Variogram for z, Lag Spacing = 10")
%%R
walk.var2 <- variogram(z ~ x+y,data=df,width=10)
%%R
plot(walk.var2,cex=1.3,pch=16,col=1,xlab="Distance",ylab="Semivariance",main="Variogram for z, Lag Spacing = 10")
%%R
model1.out <- fit.variogram(walk.var2,vgm(70000,"Sph",40,20000),fit.method=2)
%%R
plot(walk.var2, model=model1.out,xlab="Distance",ylab="Semivariance",main="Variogram for V, Lag Spacing = 10", pch=16,col=1,cex=1.3,lwd=2)
%%R
model1.out <- fit.variogram(walk.var2,vgm(70000,"Mat",40,20000),fit.method=2)
plot(walk.var2, model=model1.out,xlab="Distance",ylab="Semivariance",main="Variogram for V, Lag Spacing = 10", pch=16,col=1,cex=1.3,lwd=2)
%%R
model1.out <- fit.variogram(walk.var2,vgm(70000,"Gau",40,20000),fit.method=2)
plot(walk.var2, model=model1.out,xlab="Distance",ylab="Semivariance",main="Variogram for V, Lag Spacing = 10", pch=16,col=1,cex=1.3,lwd=2)
%%R
model1.out <- fit.variogram(walk.var2,vgm(70000,"Exp",40,20000),fit.method=2)
plot(walk.var2, model=model1.out,xlab="Distance",ylab="Semivariance",main="Variogram for V, Lag Spacing = 10", pch=16,col=1,cex=1.3,lwd=2)