Le portail francophone de la géomatique

Bonjour,

J'essaie depuis plusieurs semaines de "territorialiser" des points sur une carte de la France, je m'explique :

Objectif : Avoir une carte de la France par canton ou département avec le nombre de jours anormalement chauds (comme en page 194 ici : https://www.ecologie.gouv.fr/sites/defa … CC_WEB.pdf)

Je commence par télécharger une liste de points sur drias-climat.fr/qui indique l'écart à la normale Ecart de nombre de jours anormalement chauds (jour(s)) ainsi que la Latitude;Longitude; dans un fichier texte (WGS 84)

je l'importe comme fichier texte délimité (pas de soucis), comme suivant[img]https://i.ibb.co/KLbVTRZ/capture-tete-d-limit.png[/img] , me voila avec une vague de points en forme de France.

J'apporte ensuite un fond de la carte de la France au format .shp, SCR : Lambert 93. La superposition avec la couche précédente "colle" bien. Le niveau de la maille de cette carte étant trop fin (certaines zones n'ont pas de "points chaleurs") je regroupe par cantons ma carte de la France.

Et depuis c moment, impossible de joindre par localisation (résumé) la moyenne(ou le max) des valeurs des points de ma couche anormalement chauds ... ou en passant par joindre par localisation (non résumé) comme ceci

Les tests que j'ai pu réaliser  :
- joindre par localisation avec ma Carte France avant regroupement ( prend du temps à calculer, fonctionne assez bien mais il y a des polygones sans graduations car pas de points "chaleurs" associés)
- joindre par localisation en intersec OU / contient / ou à l'intérieur et avec plusieurs types de jointures (chevauchement,  premier élément, plusieurs éléments distincts) comme ici https://ibb.co/tqgjKt1
-joindre par localisation plus proche (en ajoutant la colonne ATnx à ma carte France_regroupé par canton)
- joindre par localisation (résumé) avec médian ou moyenne ou max (ibid qu'au dessus pour les entités)
-passer par un champ virtuel et la commande overlay_nearest ou overleay_intersec comme suivant " array_to_string(( overlay_nearest('indicesALADIN63_CNRM-CM5_vagues chaleur RCP4.5 H1',ATX35))) " (la gradation en sortie est trop faible)

et les conclusions que j'en tire :
- Pas besoin de faire une reprojection de WGS 84 vers Lambert 93 ou vice versa pour utiliser la jonction par localisation
- utiliser l'intersec pour les prédicats géométriques
- utiliser le champ virtuel avec un overleay_intersec et combiner avec un max pour n'avoir qu'un seul 1 (avec la valeur la plus élevé) pour gradué le polygone

Je suis assez novice sur QGIS et j'essaie par tâtonnements, mais je suis tombé sur trop compliqué pour moi. J'aurais souhaité avoir vos conseils pour surmonter cet obstacle.

Merci pour votre temps, n'hésitez pas à me demander plus de conseils ou de capture d'écrans pour que je puisse éclairer mon propos
Baptiste