Le portail francophone de la géomatique

Bonjour à tous,

je souhaite faire différentes carto d'évolution de variables climatiques (températures moyennes, cumul des pluies, etc.) à l'échelle de territoires, en partant des données de projection mises à disposition sur DRIAS.

Le site met à dispo les informations par points, à partir d'un ensemble de « 143*134 points, numérotés de 0 à 19161 (origine en 41.4 N / 4.1 W, résolution 8 km) », appelée grille Safran.

À l'échelle d'un territoire, lorsque je récupère un petit nombre de points, je peux assez facilement recréer une grille vecteur autour de ces points en créant une grille de 8 km, en décalant les xmin, xmax, ymin, ymax de 4 000 pour que les points soient centrés. Cette grille, en y faisant une jonction spatiale avec le point à l'intérieur de chaque polygone, me permet une représentation graphique des évolutions à partir d'une maille de 8km*8km.
Seulement, je vois bien que le décalage ne doit pas être exactement de 8 km, car même sur un territoire limité (disons un gros EPCI), je sens un léger décalage de la grille par endroit.
Ce n'est pas gênant, puisque vu la maille et le degré de précision des projections, on est plus là pour donner un ordre d'idée des évolutions qu'une analyse locale précise.

Cela dit, pour gagner du temps et ne pas recréer ma grille à chaque fois autour d'un territoire, j'ai pensé à récupérer tous les points de France et créer ma grille, créer une carte de France pour chaque indicateur qui m'intéresse, et je n'aurai plus qu'à faire une extraction à chaque fois à partir du territoire qui m'intéresse.

Sauf que là, le petit décalage se sent, puisque nulle part ou presque le point est un centroïde, et la grille se décale partout.

J'imagine qu'il doit y avoir une notion de lien avec la projection (L93) et le fait que la terre est ronde ? Des trucs que je ne maitrise pas du tout en fait Parce que si je calcule la distance entre deux points au nord-ouest de ma carte avec la formule du théorème de Pythagore, et que je fais la même chose avec deux points au sud-est, je n'ai pas la même distance (respectivement 7 994,991518 mètres et 7 991,387593 mètres).

De là, comment puis-je assez facilement recréer une grille de polygones carrés autour de chaque point qui ne se décale pas et qui s'assure à chaque fois de mettre la limite de son polygone à mi chemin avec le point suivant, et ainsi la rendre exploitable ?
À moins que la grille elle-même soit téléchargeable, mais je doute, je n'ai rien trouvé sur le site, et j'imagine qu'en mettant à disposition la grille de points, ils considèrent que cela suffit évidemment.

Pour info, pour chaque point de la grille, j'ai les infos suivantes :
idPoint    E_LambIIe    N_Lamb2e    E_Lamb93    N_Lamb93    Longitude      Latitude      Altitude

Merci par avance !

Groflo

[edit : d'ailleurs pour ceux qui voudrait faire un test, voir ce lien qui donne accès à la grille de points.]


[edit 2 : une alternative que je viens de trouver serait de changer le symbole des points directement, sans recréer de couches, en modifiant leur forme de point par une forme de carré de 8 000 mètres de côté. S'il y a ainsi de micro décalage dans les chevauchements, ils sent presque invisibles pour l'essentiel et résolvent le problème du décalage grandissant vis-à-vis des points, puisque c'est chaque point qui est directement une zone de 8km*8km !
Et en voyant l'image, je me rends compte que la grille est un peu « penchée », ce qui doit expliquer le décalage quand je crée une grille vecteur qui elle, crée une grille « droite », dans l'axe est-ouest / nord-sud.
Cela dit, si certains ont tout de même une solution vraiment propre à me proposer, je prends ! Merci]

Dernière modification par Groflo (Mon 22 October 2018 17:38)