Suite

GDistance pour déterminer le point le plus proche du polygone : la modification de la projection modifie-t-elle les résultats qualitativement ?

GDistance pour déterminer le point le plus proche du polygone : la modification de la projection modifie-t-elle les résultats qualitativement ?


j'utilisegDistancedugéospackage pour R pour trouver le point le plus proche d'un polygone. La projection des deux points et du polygone est+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0. Mon code ressemble à ceci :

min(gDistance(MonPolygone, MesPoints, byid=TRUE))

L'exécution de ceci donne un résultat et les avertissements suivants :

Messages d'avertissement : 1 : Dans RGEOSDistanceFunc(spgeom1, spgeom2, byid, "rgeos_distance") : L'objet spatial 1 n'est pas projeté ; GEOS attend des coordonnées planaires 2 : Dans RGEOSDistanceFunc(spgeom1, spgeom2, byid, "rgeos_distance") : l'objet spatial 2 n'est pas projeté ; GEOS attend des coordonnées planaires

Ma compréhension de ceci est quegDistancen'est pas satisfait de ma projection (ou de son absence). L'ampleur réelle du résultat n'est pas importante, juste que le point est en effet le point le plus proche du polygone. Ma question est la suivante : ce résultat changera-t-il si je modifie la projection ? Autrement dit, quelle que soit la valeur numérique réelle de la distance entre le polygone et le point, le même point sera-t-il identifié comme étant le plus proche si je modifie la projection ?


Plus de détails

Mes points sont des données qui ont une couverture mondiale sur les océans à une résolution de 1 degré. Mes polygones sont des régions de l'océan généralement de quelques degrés de longueur, mais généralement assez étroites. Néanmoins, étant donné que les points et les polygones proviennent de sources différentes, les deux ne s'alignent pas nécessairement - c'est pourquoi j'ai besoin de trouver les points les plus proches des polygones. Cela dit, il y aura toujours un point à moins de 1 degré d'un polygone. Mes résultats n'ont pas besoin d'être très précis (c'est-à-dire des estimations approximatives).


À la place degDistance(pour les coordonnées planes) Vous pouvez utiliserdist2Ligne(pour les coordonnées angulaires), à partir dugéosphèrepaquet. Bien qu'il s'appelledist2Ligne, la fonction fonctionne également pour les polygones (spatiaux)*.


Il ne veut pas que vous utilisiez des données non projetées car les polygones ne sont pas correctement représentés à moins que vous n'utilisiez une projection à surface égale.

Jetez un œil à Groenland vs Afrique sur une carte GCS et cela vous dit que le Groenland est plus grand que l'Afrique - même pas proche.

Ainsi, afin d'obtenir des mesures appropriées, assurez-vous d'appliquer à vos données une projection de surface égale qui convient à la surface de la Terre et à l'échelle sur laquelle vous travaillez.


Voir la vidéo: Constructing ANGLES WITH GIVEN SIZE: GeoGebra Beginner Exercise 6