Discussion :

[pâlelune]

Bonjour,

Dans R, je dispose de données géographiques de différents points inscrits dans des départements. Le souci est qu'un même point peut être à cheval entre 2 ou 3 départements.

J'ai besoin d'effectuer la correspondance entre chaque point pour un département donné pour pouvoir utiliser les données spécifiques.

Je cherche à allouer à chaque point les valeurs du département avec lequel la superposition est la plus grande.

Pourriez-vous me conseiller une approche sous R

Je vous remercie d'avance.

[ced]

Bonjour,

Comment un point peut-il être à cheval sur plusieurs départements ?
Il doit y avoir des superpositions dans votre couche de départements, ce qui semble plutôt curieux...

Par ailleurs, qu'entendez-vous par :

Je cherche à allouer à chaque point les valeurs du département avec lequel la superposition est la plus grande.

Un point est ou n'est pas dans un département (intersection avec un polygone). Il n'y a pas de notion de taille là-dedans, un point n'ayant pas de surface...

Ou alors, je n'ai rien compris.

[pâlelune]

Bonjour,

Comment un point peut-il être à cheval sur plusieurs départements ?
Il doit y avoir des superpositions dans votre couche de départements, ce qui semble plutôt curieux...

Par ailleurs, qu'entendez-vous par :

Un point est ou n'est pas dans un département (intersection avec un polygone). Il n'y a pas de notion de taille là-dedans, un point n'ayant pas de surface...

Ou alors, je n'ai rien compris.

Bonsoir Ced,

Merci de votre réponse.

En fait, il ne s'agit d'un point en tant que tel mais plutôt de mailles de 12X12 km définies sur l'ensemble de la France (désolé de n'avoir pas été très explicite, R me cause de sacrés traumatismes ).

Par ailleurs, je dispose de données de prix et je cherche à allouer à chaque maille les valeurs correspondantes à un département spécifique. J'ai pu voir que certaines de mes mailles étaient à cheval entre 2 ou 3 départements voire 4 plus rarement.

Dans ce cas de figure, je voudrais donc attribuer les données sur les prix de commodités à la maille dont la surface dans un département est plus importante par rapport à un autre ou aux autres.

Si vous pouviez me conseillez une approche, Merci par avance.

[Sengar]

Dans ce cas de figure, je voudrais donc attribuer les données sur les prix de commodités à la maille dont la surface dans un département est plus importante par rapport à un autre ou aux autres.

Salut,

Pourquoi ne pas tout simplement utiliser le milieu de la maille plutôt que de se compliquer la vie?

[pâlelune]

Salut,

Pourquoi ne pas tout simplement utiliser le milieu de la maille plutôt que de se compliquer la vie?

Bonjour SENGAR,

Je serai bien partant pour ne pas me compliquer la vie mais étant débutant sous si vous pouviez me conseiller des commandes en particulier.

Je me permets de vous exposer ma démarche.

Après avoir importé les données géographiques sous R via le package rgdal et la fonction readOGR, je les ai spécifié en tableau comportant les longitudes et latitudes des mailles, les mailles, leur identifiant et le département relié aux mailles tel que:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
> head(DataMdep)
    ID MAILLE CodeNUTS   LONG     LAT
7  120    241    FR302 559388 2659441
8  121    242    FR302 571388 2659441
9  122    243    FR301 583388 2659441
10 122    243    FR302 583388 2659441
11 123    244    FR301 595388 2659441
12 124    245    FR301 607388 2659441

Le souci comme vous pouvez avec la maille 243 c'est qu'une maille peut se trouver entre 2 et 3 départements voire 4 comme suit:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
 
 
 > table(table(DataMdep$MAILLE))
 
   1    2    3    4    5
 
2310 1256  196    7    2

Je voudrais aboutir à associer à chaque maille un département (soit en considérant le département incorporant la superficie de la maille la plus importante ou par moyenne).

Je vous remercie par avance de votre aide.[/code]

[Sengar]

Ah d'accord.

Du coup le problème c'est même pas la maille de 12*12km qui chevauche plusieurs departements.

C'est plutôt tes données qui ne permettent pas d'avoir des coordonnées assez fines pour distinguer les départements à leurs frontières.

Pas d'espoir de calculer des surfaces d'appartenance à chacun des departements de la maille dans ces conditions.
A moins que tu ais une autre base de données pour les départements qui definit plus précisément les frontières?

[pâlelune]

Ah d'accord.

Du coup le problème c'est même pas la maille de 12*12km qui chevauche plusieurs departements.

C'est plutôt tes données qui ne permettent pas d'avoir des coordonnées assez fines pour distinguer les départements à leurs frontières.

Pas d'espoir de calculer des surfaces d'appartenance à chacun des departements de la maille dans ces conditions.
A moins que tu ais une autre base de données pour les départements qui definit plus précisément les frontières?

Bonjour Sengar,

Merci de ta réponse.

Je disposais surtout d'un shapefile des départements de la France. Cependant, dans une autre base de données comportant mes données sur les prix de biens agricoles dans les différents départements, je dispose des coordonnées géographiques des différents départements. Je ne sais pas s'il est possible dans ce cadre de déterminer les surfaces de mailles ?