Pointe-in-polygone pointu et matplotlib non précis avec la géolocalisation

Je viens de faire cela pour tester si les points sont réellement à l'intérieur du triangle: xxx

Maintenant, lorsque vous utilisez Google Maps et que le polygone est mappé sur des coordonnées sphériques, le triangle est déformé, une chose à garder à l'esprit.

Quoi qu'il en soit, tracer vos données avec KML dans la Terre de Gokele montrent également le point situé à l'extérieur du triangle ?! xxx

même apparence que dans l'image MATPLOTLIB, Le point 1 est glaciaire à l'extérieur du triangle, lorsqu'il est tracé dans des coordonnées de 2D euclidiennes. Pour les calculs géométriques dans les coordonnées géométriques, vérifiez la console de Python QGIS ou des outils GDAL / OGR. Ou vous utiliseriez l'API Google Maps, comme dans l'exemple, qui est lié sur Page , où le sujet 2D-géométries vs géométries géodésiques est Coverd.

N'oubliez pas: le monde n'est pas plat! em> Si la projection de Google Maps 'est la réponse souhaitée, vous devez projeter les coordonnées géographiques sur Mercator sphérique pour obtenir un ensemble différent de coordonnées X et Y. pyproj peut vous aider, assurez-vous d'inverser vos axes de coordonnées (c.-à-d.: x, Y ou longitude, latitude).

import pyproj
from shapely.geometry import Polygon, Point
from shapely.ops import transform
from functools import partial

project = partial(
    pyproj.transform,
    pyproj.Proj(init='epsg:4326'),
    pyproj.Proj('+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs'))

poly = Polygon(([-80.190262, 25.774252], [-66.118292, 18.466465], [-64.75737, 32.321384]))
p1 = Point(-76.728515625, 27.254629577800088)

# Old answer, using long/lat coordinates
poly.contains(p1)  # False
poly.distance(p1)  # 0.01085626429747994 degrees

# Translate to spherical Mercator or Google projection
poly_g = transform(project, poly)
p1_g = transform(project, p1)

poly_g.contains(p1_g)  # True
poly_g.distance(p1_g)  # 0.0 meters