Кратчайшее расстояние от точки до контура полигона
Добавлено: 02 июл 2015, 18:24
Как найти кратчайшее расстояние от точки до контура замкнутого полигона?
Точка находится вне контура полигона.
Точка находится вне контура полигона.
Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
https://gis-lab.info/forum/
nickleb писал(а):Как найти кратчайшее расстояние от точки до контура замкнутого полигона?
Точка находится вне контура полигона.
Код: Выделить всё
import com.vividsolutions.jts.operation.distance.DistanceOp;
...
double distance = DistanceOp.distance(polygon, point);Код: Выделить всё
import com.vividsolutions.jts.operation.distance.DistanceOp;
...
double distance = DistanceOp.distance(polygon, point);Код: Выделить всё
>>> from shapely.geometry import Point, Polygon
>>> polygon = Polygon([(30, 10), (30, 30), (10, 30)]) # Implicitly closed
>>> point = Point(10, 20)
>>> point.distance(polygon)
7.0710678118654755
Как правильно подсказывают, это JTS.nickleb писал(а):Спасибо, Sergey. Это geolatte? Это я просто по-Googl'ил.... Если я правильно предполагаю... А как это попробовать использовать под Wind'ами?
Вставить в код программы в нужном месте.nickleb писал(а):А как это попробовать использовать под Wind'ами?
Да, Sergey, уже понял по Вашим комментариям и ссылкам дорогих и уважаемых gis-lab-овцев, что ничего - справиться можно. А задача конкретная такая: льдина дрейфует в океане и несёт на себе массово-балансовый буй и его ежечасные обсервации доступны по http. Требуется по этим обсервациям:Sergey Astakhov писал(а):Как правильно подсказывают, это JTS.nickleb писал(а):Спасибо, Sergey. Это geolatte? Это я просто по-Googl'ил.... Если я правильно предполагаю... А как это попробовать использовать под Wind'ами?
[quote=
Вставить в код программы в нужном месте.
Может вы для начала расскажите саму задачу? А то библиотек много всяких разных, под кучу языков.
Код: Выделить всё
library(plotKML)
library(sp)
temp <- tempfile()
download.file("http://imb.erdc.dren.mil/irid_data/2015D_clean.csv",temp)
buoy <- read.csv(temp,header=T,sep=",")
unlink(temp)
#Какие колонки в файле?
names(buoy)
coordinates(buoy) <- ~Longitude..degrees.+Latitude..degrees.
proj4string(buoy) <- CRS("+proj=longlat +datum=WGS84")
plotKML(buoy, var.name="Date")