from IPython.display import Image
Image(url='http://www.portailsig.org/sites/default/files/zeropoint_logo.png')

import os
os.environ['GISBASE'] = '/Applications/GRASS-6.4.app/Contents/MacOS'
os.environ['PYTHONPATH']="/Applications/GRASS-6.4.app/Contents/MacOS/etc/python/"
os.environ['LD_LIBRARY_PATH']="/Applications/GRASS-6.4.app/Contents/MacOS/lib"
os.environ['GIS_LOCK'] = '$$'
os.environ['GISRC'] = '/Users/martinlaloux/.grassrc6'
import grass.script as grass


gisbase = os.environ['GISBASE']
gisdb="/Users/martinlaloux/grassdata"
location="geol"
mapset="test_python"
import grass.script.setup as gsetup
gsetup.init(gisbase, gisdb, location, mapset)

debligne =grass.read_command("v.to.db", flags="p", map="testgrass", type="line", option="start", units="meters" , quiet=True)



debx=float(debligne.split("|")[1])
deby=float(debligne.split("|")[2])
debz=float(debligne.split("|")[3])

finligne=grass.read_command("v.to.db", flags="p", map="testgrass",type="line", option="end", units="meters", quiet=True)

finx=float(finligne.split("|")[1])
finy=float(finligne.split("|")[2])
finz=float(finligne.split("|")[3])
 

from shapely.geometry import Point

Point1 = Point(debx,deby)

Point2 = Point(finx,finy)

from shapely.geometry import LineString

line = LineString([Point1,Point2])

%matplotlib inline

print line

import matplotlib.pyplot as plt 
fig = plt.figure() 
ax = fig.gca() 
x, y = line.xy
ax.plot(x, y, color='b')

load_ext rpy2.ipython

%%R -w 480 -h 480 -u px
library(ggplot2)
dat <- data.frame(x = rnorm(2), y = rnorm(2), 
                  lab = sample(c('deb', 'fin'), 10, replace = TRUE))
x <- ggplot(dat, aes(x = x, y = y, color = lab)) + geom_point() + geom_line()
print(x)

%load_ext oct2py.ipython

%octave -i x x

%octave -i y y

%%octave -s 480,480 -f svg
plot(x,y);

%Rpush x y

%R lm(y~x)$coef

Xr = x - x.mean(); Yr = y - y.mean()
pente = (Yr*Yr).sum() / (Xr**2).sum()
interception = x.mean() - x.mean() * pente
(interception, pente)

%R resid(lm(y~x)); coef(lm(x~y))

b = %R a=resid(lm(y~x))

print b

v1 = %R plot(x,y); print(summary(lm(y~x))); vv=mean(x)*mean(y)
print 'v1 est égal à :', v1
v2 = %R mean(x)*mean(y)
print 'v2 est égal à:', v2