28 Nov
2022
28 Nov
'22
11:56 a.m.
Message posté par : Thais781
----------------------------------------
Bonjour,
Cf. l'image ci-joint.
J'ai un polygon (en bleu) et un point (celui en orange)
Et je cherche a calculer les points de contact les plus proches entre ce point et le
polygon
Pour le moment, je trace un cercle a partir du point orange et je fais évoluer le diamètre
jusqu'a qu'il touche (intercepts) le polygone et j'ai donc le point de
contact.
Cette solution est un peu ... bourin ... peut être qu'il existe une solution plus
élégante pour obtenir ces point (les verts dans le dessin)
N'hésitez pas si vous avez une idée ;-)
Thais
----------------------------------------
Ce message est accompagné de fichiers, pour les télécharger, suivre le lien ci-dessous.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=356620#p356620
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
28 Nov
28 Nov
1:37 p.m.
Message posté par : Marion Robin
----------------------------------------
Bonjour,
Il y a la fonction Shapely nearest_points qui fait ça. Je viens de tester et ça a
l'air de plutôt bien fonctionner :-)
import geopandas as gpd
from shapely.geometry import Polygon, Point
from shapely.ops import nearest_points
poly = Polygon([(4.338074,50.848677), (4.344961,50.833264), (4.366227,50.840809),
(4.367945,50.852455), (4.346693,50.858306)])
poly = gpd.GeoSeries(poly, crs='EPSG:3395')
point = gpd.GeoSeries([Point(4.382617,50.811948)], crs='EPSG:3395')
nearest_edge, _ = nearest_points(poly[0], point[0])
nearest_edge = gpd.GeoSeries(nearest_edge, crs='EPSG:3395')
ax = poly.plot()
point.plot(ax=ax, color='red')
nearest_edge.plot(ax=ax, color='green')
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=356626#p356626
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
29 Nov
29 Nov
4:04 p.m.
Message posté par : Thais781
----------------------------------------
Merci Marion
C'est effectivement beaucoup plus propre comme ca ;-)
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=356691#p356691
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
1047
days inactive
1048
days old
2 comments
1 participants
participants (1)
-
geolibre＠georezo.net