Calcul de surface de polygone

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Cyrielle Jac
Messages : 55
Enregistré le : 13 Mar 2017, 08:30

Calcul de surface de polygone

Messagepar Cyrielle Jac » 21 Avr 2017, 10:53

Bonjour,

J'ai représenter deux polygones sur un plot avec les lignes :

plot(pcaxes[,1:2], type="n", ylim=c(-3.5,4.8), xlim=c(-2.3,6.2))

### Faire seulement deux polygones
## 1
tata<-toto[,4:5]
pts=chull(tata[,1],tata[,2])
around.points=cbind(tata[pts,1],tata[pts,2])
polygon(around.points[,1],around.points[,2],col=rgb(0.6,1,0))

## 2
tutu<-titi[,4:5]
pts=chull(tutu[,1],tutu[,2])
around.points=cbind(tutu[pts,1],tutu[pts,2])
polygon(around.points[,1],around.points[,2],col=rgb(0.6,0.5,1,0.4))


tutu et tata sont deux dataframe de forme :
Dim1 Dim2
Sp1 -1.25 0.02
Sp2 0.23 2.54
Sp3 0.43 1.68

Et je voudrais calculer la surface de chacun de ces deux polygones ainsi que la surface de leur interaction. Je suis tomber sur un précédent post à ce sujet mais je ne suis pas arriver à appliquer le script à mon script. quelqu'un a t-il une solution?

Bonne journée

Cyrielle

Serge Rapenne
Messages : 1426
Enregistré le : 20 Aoû 2007, 15:17
Contact :

Re: Calcul de surface de polygone

Messagepar Serge Rapenne » 21 Avr 2017, 13:57

Bonjour,

Il me semble que la solution la plus simple est de passer par les fonctions spatiales. cf par ex : https://stackoverflow.com/questions/997 ... erlap-in-r

Serge

Cyrielle Jac
Messages : 55
Enregistré le : 13 Mar 2017, 08:30

Re: Calcul de surface de polygone

Messagepar Cyrielle Jac » 24 Avr 2017, 09:18

Bonjour,

J'avais déjà regarder avec ces packages mais je dois mal les utiliser car rien ne fonctionne :
Error in identical(spgeom1@proj4string, spgeom2@proj4string) :
trying to get slot "proj4string" from an object of a basic class ("matrix") with no slots

Pouvez-vous m'aidez?

Bonne journée

cyrielle


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité

cron