给定一定的LineString和点p
from shapely.ops import nearest_points
from shapely.geometry import Point
p = Point(51.21745162000732, 4.41871738126533)
linestring = LineString([(51.2176008, 4.4177154), (51.21758, 4.4178548), (51.2175729, 4.4179023), (51.21745162000732, 4.41871738126533)])
与p最近的点是通过以下方法计算的：
n_p
浏览 0提问于2020-05-15得票数 1
回答已采纳
当我尝试在remote experitest its not working search filed中执行相同的脚本时，我有一个搜索字段，它在桌面版本的chrome中工作 1)我正在传递来自json的数据，它在桌面版本和更低版本的chrome中工作得很好。2)当我试图从chrome版本91远程执行它时，它通过存储而不是从json中获取数据来获取最近的数据。3)我尝试了等待命令，但它仍然不起作用 if (Region().contains("IN")) {
setText(this.gbnlPostcodeinputfield, Location);
浏览 9提问于2021-08-21得票数 0
我想使用osmnx从海地提取建筑层信息。我在不同的城市工作，而不是
import osmnx as ox
place_name = "Port-au-Prince"
buildings = ox.buildings_from_place(place_name)
我得到以下错误
TypeError: Geometry must be a shapely Polygon or MultiPolygon
浏览 2提问于2019-04-07得票数 0
1回答
我需要知道城市的纬度和经度。(使用xcode) google有自己的api地理编码，但它不想返回最近的城市(我只需要最近的城市)，它返回cyrilic的地名。如何解决我的问题？
浏览 3提问于2013-01-21得票数 0
回答已采纳
我有一个很长的线串，我用shapelies closest_points在线串上得到了两个点。这似乎工作得很好，但是调用带有直线和最近点之一的.split()作为args会返回原始对象，没有拆分。为什么会这样呢？ shapely_track = wkt.loads(act.track.wkt)
shapely_seg = wkt.loads(seg.track.wkt)
startpoint, endpoint = shapely_seg.boundary
segment_start, point = nearest_points(shapely_tra
浏览 13提问于2021-04-08得票数 0
2回答
我正在尝试加入两个空间数据集。第一个包含点和第二个多边形。
然而，其中一些点在多边形之外。
是否有一种简单的方法将这些点连接/截取到最近的多边形边界，而不是最近的多边形质心？
目前，我正在加入最近的多边形质心，但这并不产生我正在寻找的结果。
浏览 7提问于2020-11-18得票数 1
1回答
我正在使用CLLocationManager来尝试查找用户当前所在的地址。如果没有可视化的地图视图，我如何将类返回的坐标转换为实际的地址/业务？或者是否有一个不同的类来返回此数据？谢谢。
浏览 2提问于2011-07-23得票数 0
回答已采纳
我将.csv文件作为数据帧读取，如下所示：
import pandas as pd
df = pd.read_csv('myFile.csv')
df.head()
BoroName geometry
0 Brooklyn MULTIPOLYGON (((-73.97604935657381 40.63127590...
1 Queens MULTIPOLYGON (((-73.80379022888098 40.77561011...
2 Queens MULTIPOLYGON (((-73.8610972440186 40
浏览 0提问于2020-04-09得票数 14
1回答
我一直试图用形状来寻找线和多边形的交集，但是我遇到了一些浮点数的问题。
示例代码：
polygon = [(4.0, -2.0), (5.0, -2.0), (4.0, -3.0), (3.0, -3.0), (4.0, -2.0)]
shapely_poly = shapely.geometry.Polygon(polygon)
line = [(4.0, -2.0000000000000004), (2.0, -1.1102230246251565e-15)]
shapely_line = shapely.geometry.LineString(line)
intersection_l
浏览 2提问于2013-01-21得票数 12
回答已采纳
3回答
我希望使用Geopandas / Shapely实现与ArcPy 相当的功能。我对Geopandas和Shapely非常陌生，并且已经开发出了一种有效的方法，但我想知道是否有更有效的方法来做到这一点。
我有两个点文件数据集-人口普查块质心和餐馆。我正在寻找，对于每个人口普查区的质心，到它最近的餐厅的距离。对于同一家餐厅在多个街区是最近的餐厅，没有任何限制。
对我来说，这变得有点复杂的原因是因为是按元素计算的，基于索引进行匹配。因此，我的一般方法是将Restaurants文件转换为多点文件，然后将块文件的索引设置为相同的值。那么所有的街区质心和餐馆都有相同的索引值。
import pandas
浏览 1提问于2015-06-10得票数 15
回答已采纳
我试图找出直线字符串的起始点和用来寻找多边形最近距离的点之间的直线串的长度。
因此，我使用了下面的代码来获得行字符串和一些多边形之间的最小距离。
gdf['MinDistToTrack'] = gdf.geometry.apply(lambda l: min(rail_or.distance(l)))
我还想得到从行字符串开始到上面代码所使用的点的距离。
现在，我得到了包含具有值'MinDistToTrack‘的多边形的数据文件(我现在已经有了)，但也包含了一个值MinDistToTrack
所以，假设从线串开始到多边形，沿着线串的路径有22米，那么这就是我想与'
浏览 6提问于2022-03-02得票数 0
我正在使用cartopy来绘制地图，它以前工作得很好。但是最近，当我在JupyterLab中使用这个包时，我得到了以下错误：
'AttributeError: 'LGEOS360' object has no attribute 'GEOSBufferWithParams'.
我最近更新了我的Mac OS系统到Big Sur。我不确定这个错误是否是由升级引起的。如果是这样，我应该怎么做？
我已经尝试卸载和安装cartopy和Shapely。但这两种方法都没有奏效。
浏览 1提问于2020-11-27得票数 2
我得找出这点是否属于某一地区。信息由geojson文件提供。我通过以下方式创建了DataFrame：
choro = pd.concat([pd.DataFrame({'name':block['properties']['name'], 'area': shape(block['geometry'])}) for block in geojson['features']], ignore_index=True)
该代码适用于警告：
ShapelyDeprecationWarning:对多部分几何
浏览 2提问于2021-10-27得票数 1
1回答
我有两对lat/lon (以十进制表示)以及它们的半径(以米表示)。我想要实现的是找出这两个出口之间是否存在一个交集(当然，这在这里是不成立的，但我们的计划是在许多其他数据点上尝试这个算法)。为了检查这一点，我使用了Shapely的intersects()函数。然而，我的问题是，我应该如何处理不同的单位？我是否应该先做某种变换\投影(对于lat\lon和radius都是相同的单位)？
48.180759,11.518950,19.0
47.180759,10.518950,10.0
编辑：
我在这里找到了这个库()，这似乎很有帮助。但是，我不能100%确定我是否正确地应用了它。有什么想法吗？
浏览 3提问于2014-12-11得票数 3
回答已采纳
1回答
我第一次处理空间数据。我必须比较两个有lat和长细节的数据格式。我已经将两者转换为GeoPandas数据格式。
import pandas as pd
from pandas import DataFrame
import geopandas as gpd
from neighbors import nearest_neighbor
df = pd.DataFrame([[1973,22.525158,88.330775],[1976,72.85136,19.10840],[898,91.78523,26.15012]],columns=['id', 'lat
浏览 5提问于2021-03-04得票数 2
回答已采纳
1回答
我有一个由lat和lon给出的点，我想通过最小欧几里德距离找到离该点最近的边。例如 import osmnx as ox
track = [(40.7052, -74.0069)]
fig, ax = ox.plot_graph(G, show=False, close=False)
for pairs in track:
ax.scatter(pairs[1], pairs[0], c='red')
plt.show()
ox.distance.get_nearest_edge(G, track, return_geom=True, return_dist=Tru
浏览 7提问于2020-09-28得票数 0
回答已采纳
1回答
我尝试在Python中将文本文件中的数据(用纬度、经度和花粉通量值组织)绘制为栅格。我在https://autogis-site.readthedocs.io/en/latest/notebooks/L5/02_interactive-map-folium.html上使用Choropleth的代码来尝试显示数据。我的GeoPandas geodataframe具有点几何图形；但是，本教程中的点的几何图形看起来已经是多多边形，我假设它们是网格中的正方形。如何将我的数据(假设每个纬度/经度点都是网格中像素的中心)转换为网格geopandas (geodataframe)数据？我将使用的投影是La
浏览 114提问于2021-07-13得票数 0
回答已采纳
1回答
我试着在多边形接触的条件下创建它们的区域。在我的示例中，我有一个包含382个多边形的示例数据集，这些多边形需要组合在一起(但整个数据集包含6355个多边形)。(我想展示一张照片，但我没有足够的名气来这么做。)
我想做这种蛮力，但当然这需要很长时间，而且不是最优的。
def groupBuildings(blds):
# blds is a list with shapely polygons
groups = []
for bld in blds:
group = []
group.append(bld)
for othe
浏览 2提问于2017-02-18得票数 1
1回答
哪个GEO实现将用于数百万个点
、、、、
我试图找出哪个GEO的实现，以找到最近的点基础上的长/拉点到某个点。我将有数百万，如果不是数十亿不同的纬度/经度点，将需要比较。我一直在寻找许多不同的实现来完成我需要完成的工作。我研究过Postgis (看起来它很受欢迎，性能也很好)，Neo4J (图形数据库是一个新概念，我不知道它们是如何实现的)、AWSDynamoDBGeo散列(Scales非常好，但只有库是用node.js编写的，我希望用node.js编写一个库)等等，但无法确定哪种方法性能最好。我纯粹是在观察性能，而不是大量的特性。我所需要做的就是将一个点与所有点进行比较，找到最近的点(读操作)，并且能够快速地更改数据库中的一个点(写
浏览 0提问于2014-03-24得票数 5
回答已采纳
我在用火药库，我在想怎么做。
获取所有用户的纬度和逻辑，并将其保存到firebaseDB
在firebaseDB中对离我最近的用户排序
获取数据并显示最近的用户
但是我在文档中找不到数据库中的sotring，有什么方法可以在firebaseDB中对最近的用户排序吗？还是其他方式？
浏览 0提问于2018-05-18得票数 0
回答已采纳
4回答
形状和几何位置不精确的点多边形
、、、、
我正在用matplotlib和shapely测试点内多边形函数。
这里是一个包含百慕大三角形多边形的。
Google映射了的点多边形函数，清楚地显示出testingPoint和testingPoint2位于多边形内部，这是一个正确的结果。
如果我在matplotlib和shapely中测试这两个点，那么只有point2通过测试。
In [1]: from matplotlib.path import Path
In [2]: p = Path([[25.774252, -80.190262], [18.466465, -66.118292], [32.321384, -64.75737]])
浏览 6提问于2014-01-24得票数 9
回答已采纳
1回答
也许有人知道些什么，因为我找不到任何对我有意义的东西。我有一个数据集位置(经度，经度)，我想将它们捕捉到最近的道路上，并计算它们之间的距离。 到目前为止，我发现了OSM，但是我找不到一个关于如何通过python使用API的工作示例。 如果你们中的任何人能帮上忙，我非常感谢你们的细微细节。 在此期间，我会尝试自己找出答案，如果成功的话，我会发布答案(找不到任何类似的问题，所以它可能会在未来对某人有所帮助)
浏览 22提问于2020-09-28得票数 0
我正在试着解决一个问题，那就是如何通过动画字体来切换点击。 但是，为了避免这个问题，我在我的主题文件夹中创建了一个文件，并将其命名为Menu Bars。在此文件中，我放置了以下代码： <!DOCTYPE html>
<html>
<head>
<meta name="viewport"
content="width=device-
width, initial-scale=1">
<style>
.container {
display: inline-block;
cursor:
浏览 17提问于2020-04-08得票数 1
1回答
我想了解Shapely是如何工作的。
我可以用下面的代码画一条简单的线：
import matplotlib.pyplot as plt
A = Point(0,0)
B = Point(1,1)
AB = LineString([A,B])
plt.plot(AB)
然而，当我改变坐标时：
A = Point(1,0)
B = Point(3,4)
AB = LineString([A,B])
plt.plot(AB)
Shapely决定画两条线，这是我不理解的行为。
使用Shapely 1.7.0
浏览 10提问于2020-02-08得票数 1
回答已采纳
我正在尝试根据python中的一个形状文件来裁剪线条。我有一个可以工作的脚本，但它非常慢。有没有更快的方法来做到这一点？一个先决条件是我必须用python来做这件事。 from shapely.geometry import LineString, shape, Point
from shapely.affinity import rotate
from itertools import chain
from shapely.ops import cascaded_union
from fiona import open
import numpy as np
# open the shp-
浏览 98提问于2021-09-05得票数 0
回答已采纳
我有一个国家和坐标(country_gdf)列的数据格式。我有一个包含城市、纬度和经度(cities_gpd)列的单独数据格式。
我想在国家数据中添加一个新的列，用于统计一个国家内的城市数量。但是，我不知道如何处理坐标
我的样本数据如下；
import pandas as pd
import geopandas as gpd
from shapely import wkt
# Sample Country dataframe
country = pd.DataFrame(
{'Country': ['Argentina', 'Brazil
浏览 8提问于2022-07-20得票数 2
回答已采纳
1回答
我正在开发一个使用mongoose和MongoDB进行存储的node.js web应用程序。这是个地理定位应用。所以有时我必须找到存储在数据库中最近的10个位置。我有这个方法来计算这个距离：
var haversine = function(p1, p2) {
var R = 6371; // km
var dLat = toRad(p2.lat-p1.lat);
var dLon = toRad(p2.lon-p1.lon);
var lat1 = toRad(p1.lat);
var lat2 = toRad(p2.lat);
var a = Math.sin(dLat/2) * Ma
浏览 2提问于2012-08-14得票数 1
回答已采纳
我正在处理一些xarray数据数组，它们在特定的纬度/经度坐标上有一些数据。
对于数据数组1 (da1)中的每个lat/long坐标对，我希望在数据数组2 (da2)中找到最近的lat/long对。
遵循，一个似乎有效的解决方案是：
lats = xr.DataArray(da1.latitude.data, dims='z') #'z' is an arbitrary name placeholder
lons = xr.DataArray(da1.longitude.data, dims='z')
data = da2.sel(latitud
浏览 1提问于2022-01-26得票数 1
回答已采纳
1回答
如何从位置列表中得到最近的区域？
、、、、
我有一个API，它返回一个城市内与该地区的天气有关的不同区域的列表。我想根据我现在的位置得到最近的区域。
API返回
区域
纬度
经度
天气
如何根据这些数据找到最近的区域？
浏览 2提问于2012-08-19得票数 0
回答已采纳
1回答
我正在考虑使用Titan创建一个可伸缩的地理空间数据存储(我认为是R树)。在文档中，有一个GeoShape查询，文档中说泰坦可以使用Lucene或ElasticSearch处理geo数据。但是，这似乎非常慢，因为在cassandra中遍历节点本质上是在cassandra中执行连接查询，这是一个非常糟糕的主意。我想我可能误解了数据表示。
我读了，但我还是不太明白。如果所有的边缘都存储在卡桑德拉行中，那么土卫六仍然必须在顶点表上“加入”。解决这一问题的一种方法是使列值等于edge属性数据，然后您可以将顶点数据和边缘数据整齐地打包到行中。但是，当您想要执行比一个节点更深的查询时，就会发生故障，我们又
浏览 5提问于2014-03-15得票数 6
回答已采纳
1回答
我正在使用Presto 329，并试图通过替换有效的每行内部查询的交叉连接来重写以下交叉连接以提高效率： SELECT
a.store_id,
min(great_circle_distance(a.latitude, a.longitude,
CAST(b.latitude AS double), CAST(b.longitude AS double))) AS distance
FROM store_location a
CROSS JOIN places_of_interest b
GROUP BY a.store_id; 作为解释，store_loca
浏览 17提问于2021-01-28得票数 0
1回答
道路交叉路口建设中的误区
、、、
我试着把原始的道路数据分割成几段。首先，我创建了road_segments table.Then，我使用shapely将路况表的内容转换成一系列连接的路段。
对于每一条路，我首先将道路的名称和几何图形载入内存。
这里的代码工作正常。
我得到了道路的LineString几何，我想分裂在每一点，它接触或交叉另一条道路。为此，我试图为这条道路建立一个十字路口清单。
crossroads = []
cursor.execute("SELECT ST_AsText(centerline) FROM ROADS" +
"WHERE ST_T
浏览 0提问于2016-05-25得票数 1
回答已采纳
1回答
我正在使用OSMnx软件包来解决以下任务：-地图上有一个由经度和纬度定义的点X-我们需要检测包含该点X的多边形，并受相邻道路的约束-所以基本上点X在多边形内部，相邻道路将是该多边形的边界。 到目前为止，我只成功地在地图上绘制了图形的可视化，并找到了与点X最近的边/节点。 在附加的图像中，我已经突出显示了我想要提取的红色区域。 ?
浏览 20提问于2019-12-11得票数 2
回答已采纳
1回答
我想我不能问谷歌是如何做反向地理编码的，因为除了谷歌没有人会知道这一点。但是，像OSM (OpenStreetMap)这样的开源软件如何将GPS坐标解析成一个位置呢？
我知道在OSM (和Google Maps)中，每个地址都被认为是一个分配了GPS坐标的地点(或点)。反向地理编码实际上返回的不是具有相同坐标的精确点，而是最近的点。那么，这就是全部内容了吗？遍历数据库中的每个已知点并返回储藏点？或者是复杂得多？另外，它如何知道每个点属于哪个城市/国家？
提前谢谢。
浏览 0提问于2014-02-07得票数 0
误差
TypeError: 'LineString' object is not iterable
我试图在这个.shp文件中找到多边形的右下角。此.shp文件为正方形，但其他文件可能为矩形/三角形。
我想使用Geopandas，显然，我可以使用read_file()方法获得这个结果。我对SHPs很陌生，我有.shx，.dpf文件，但是当我在这个方法中只输入.shp时，我无法遍历多边形坐标。
下面是我的代码--我想在一个变量中捕获右下角，当前all_cords将捕获所有这些变量，所以我需要找到一种方法来克服这个错误，然后捕获右下角，
import geopandas as gp
浏览 25提问于2022-05-03得票数 0
2回答
我有一个应用程序托管在2个虚拟机在不同的地区(例如美国和澳大利亚)。由于数据驻留限制，用户的数据将位于这两个区域中的任何一个区域。
假设一个用户的数据在AU服务器上，正在前往加拿大，基于地理/延迟的路由将用户路由到最近的美国服务器。
我的应用程序可以在内部以某种方式响应，将此用户重定向到AU服务器。如何让所有后续请求直接进入AU服务器，而不是进入US服务器？
浏览 21提问于2021-10-08得票数 0
我有两个不同长度的数组。一个包含形状多边形，另一个包含形状点。我希望为两个数组中的每个元素组合运行a_polygon.contains(a_point) shapely函数。
我认为是构建包含行中所有可能组合的两列矩阵的理想中间步骤。但是，当输入数据很大时，‘cartersian(数组)’函数中的循环可能会影响性能。
我尝试广播其中一个数组，然后应用shapely函数：
Polygons_array[:,newaxis].contains(Points_array)
但这当然是行不通的。我知道最近发布的地质公园图书馆，但它不是我的阁楼安装的选择。
浏览 0提问于2014-07-24得票数 4
回答已采纳
我需要一个0.1°x0.1°分辨率的全球网格场。每个像素应该包含像素所在的国家/地区的名称。
我尝试使用Python语言中的reverse_geocode包生成字段。问题是，当您进入国际水域的某个位置时，包中提供的是最近国家的名称，而不是简单的“国际水域”(参见下面的示例代码)。我也考虑过使用google，但我想由于超过600万像素还有待确定，商业服务可能不是最佳选择。
import reverse_geocode as rg
# coordinates in atlantic ocean
lat_water = 35
lon_water = -29
coordinates = [(lat
浏览 31提问于2021-02-24得票数 0
回答已采纳
我是一个用python制作/绘制地图的新手，我一直在尝试按照这篇博文来生成一个世界地图()。我在这里遇到了一些问题：
安装底图(用于在地理投影上绘制数据的Matplotlib扩展)。
从mpl_toolkits.basemap导入底图回溯(最近一次调用)：文件"geos_demo.py"，第1行，从mpl_toolkits.basemap导入底图ImportError:没有命名为mpl_toolkits.basemapInstall的模块，但出现以下错误1然后我进入geos.py文件，将geos.dll的所有引用更改为geos_c.dll，如此处的链接所建议的：However，
浏览 263提问于2012-09-25得票数 25
我们已经使用Google Geocoding API很长时间了，没有出现任何问题。就在最近，由于各种原因，我们大幅减少了通话次数。你可以从这张30天的图表(蓝色)中看到：
但奇怪的是，正如我们所做的那样，延迟增加了三倍：
有什么想法可以解释为什么会这样吗？
浏览 1提问于2018-05-31得票数 0
1回答
我有两个数据帧，在第一个中，我有线串，在第二个中，我有很多点。我想找一个离线串最近的点。我试过一些东西，但我猜它不起作用。我该怎么做呢？
下面是我的代码：
def find_keybyvalue (value,dict):
for s,d in dict.items():
if value==d:
return s
liste={}
lis={}
for i in range(first.count()):
for n in range(second.count()):
point = second.geometry
浏览 11提问于2021-09-06得票数 0
我有一个系统，我喜欢根据纬度和经度数据找到离某人最近的人。我知道如何处理Mongodb中的位置数据，以及如何按半径等为每个人查找最近的人，但我不知道如果系统中有大量的用户，并且在postgresql中查询位置数据是否有可能使我的系统崩溃？哪个dbms很适合处理位置数据，Postgresql还是Mongodb？我很感激你的帮助。
非常感谢
浏览 0提问于2018-01-25得票数 2
回答已采纳
1回答
我正在尝试构建一个基本的平均堆栈应用程序，并且对于如何在前端显示存储在mongodb中的信息感到非常困惑。创建数据库查询并在视图中显示该查询的建立过程是什么。服务器应该发出db请求并以角度显示信息吗？在所有的教程中，我一直在研究mongodb和角/前端之间的交互，似乎非常模糊。
目前(我正在设置一个网球梯子应用程序)，我只想通过查询数据库中的信息来显示索引页面上的球员列表。这是正确的做法吗？
对于这个含糊不清的问题，我表示歉意，但我真的很难理解这是如何实现的。我最近一直使用Rails作为比较(虽然显然是非常不同的)，在这种情况下，您可以使用ActiveRecord查询数据库并根据查询结果定义一
浏览 1提问于2015-01-03得票数 2
回答已采纳
在我的移动应用程序中，我使用以下代码获得了用户的地理位置。我也有附近城市的地理位置(纬度，经度)。我可以将城市的地理位置数据放在数组或SQLite数据库中。如何使用jQuery或JavaScript获取和显示用户位置附近的3个最近城市。我没有使用谷歌地图。
$(document).ready(function() {
var Geo={};
if (navigator.geolocation) {
navigator.geolocation.getCurrentPosition(success, error);
}
//Get the lati
浏览 2提问于2016-03-10得票数 0
我们使用一个库来检查某些随机点是否存储在形状文件中的某些禁止区域中。
with fiona.open(path) as source:
geometry = get_exclusive_item(source[0])
geom = shapely.geometry.shape(geometry['geometry'])
def check(lat, lng):
point = shapely.geometry.Point(lng, lat)
return not geom.contains(point)
但是，最新的电话geom.contains(
浏览 4提问于2016-03-29得票数 0
回答已采纳
我有一个代码来执行一组三角形的合并，在python中使用函数cascaded_union()。
print cascaded_union(triangle)
当我打印它时，它会给出以下输出：
POLYGON ((211 55, 0 0, 150 204, 185 216, 215 159, 224 132, 211 55))
如何将这些点存储在具有x和y坐标的二维数组中？
浏览 0提问于2016-02-28得票数 1
回答已采纳
我的问题是:如何确定Aside和Bside的的哪一边是任意的LineString的“左”和“右”的边，然后再把这个几何图形分开呢？
为了解决这个问题，“左”和“右”被定义为当从一个节点“行走”到另一个节点时，一个LineString拆分器的左手和右边。
我创建了这个函数，用于将任意的shapely几何(非集合)分成两面--“左”和“右”：
import shapely.geometry as geo
import shapely.ops as ops
def splitLR(geom, splitter):
"""Split a geometry into a
浏览 3提问于2018-05-17得票数 2
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
