Skip to Content

Nous adressons toutes nos pensées à la famille de notre ami Jérôme !

http://www.forumsig.org/showthread.php/43488-Disparition-de-Phoenix

aspace_SDE.rsx

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