aspace_SDE.rsx

##Ellipse=group
##Layer=vector
##Field=Field Layer
##SDE_shape = output vector
##weighted=boolean True.
# weighted: pondéré par le champs Field ou non pondéré
# le fichier SDE_Output.txt est nécessaire pour la commande,même s'il ne sert à rien ici
library(aspace)
coords<-coordinates(Layer)
if(weighted) {
calc_sde(id=1, filename="/Users/nouvportable/SDE_Output.txt"", centre.xy=NULL, calccentre=TRUE, weighted=TRUE, weights=Layer[[Field]], points=coords, verbose=FALSE)
}else {
calc_sde(id=1, filename="/Users/nouvportable/SDE_Output.txt", centre.xy=NULL, calccentre=TRUE, weighted=FALSE, weights=NULL, points=coords, verbose=FALSE)
}
aa = c(sdeloc[2],sdeloc[3])
sp = SpatialPolygons( list(  Polygons(list(Polygon(aa)), 1)))
res = SpatialPolygonsDataFrame(sp, sdeatt)
proj4string(res) <- Layer@proj4string
SDE_shape= res