首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >问答首页 >找出构成邻里区域的每个点的纬度和经度?

找出构成邻里区域的每个点的纬度和经度?
EN

Stack Overflow用户
提问于 2019-03-24 08:45:06
回答 2查看 508关注 0票数 0

我试图找到沿特定邻域边界的每个点的纬度和经度。例如,我尝试使用geopy -当我输入一个特定的邻域时,它会返回一对坐标。我希望收到一个坐标列表,该列表将勾勒出指定邻域的轮廓。

具体地说,我正在努力获得纽约市曼哈顿的社区。

谢谢。

代码语言:javascript
运行
复制
from geopy.geocoders import Nominatim
geolocator = Nominatim()

location = geolocator.geocode('Gramercy NYC')
print(location)

输出:位置(美国纽约市纽约市曼哈顿曼哈顿社区委员会6号,纽约市,纽约市,(40.7355189,-73.9840794,0.0))

我正在得到一组坐标,我想要多个坐标。

EN

回答 2

Stack Overflow用户

发布于 2019-03-25 01:16:01

geopy的Nominatim类支持返回完整的几何数据,而不是单个点:

代码语言:javascript
运行
复制
from geopy.geocoders import Nominatim
geolocator = Nominatim(user_agent="specify_your_app_name_here")

location = geolocator.geocode('Gramercy NYC', geometry='wkt')
geometry = location.raw['geotext']

geometry将包含以下内容:

代码语言:javascript
运行
复制
'POLYGON((-73.9901251 40.7377992,-73.9869574 40.736466,-73.9887732 40.7339641,-73.9825505 40.7313605,-73.9785138 40.7368725,-73.9847526 40.7395063,-73.9856806 40.7382199,-73.9873061 40.7389028,-73.9877245 40.7383154,-73.9897603 40.7391731,-73.9901251 40.7377992))'

几何图形支持的格式为wktsvgkmlgeojson。有关更多详细信息,请参阅the Nominatim.geocode docs

您可能还想考虑先在https://nominatim.openstreetmap.org/页面上尝试查询,它提供了一个方便的web界面,可以显示这些几何图形。

票数 1
EN

Stack Overflow用户

发布于 2019-03-24 09:44:47

我不知道geopy能提供这一点。

你真正需要的是你感兴趣的地理位置的shapefile。美国人口普查局为不同类型的地理位置提供了其中的几个:Cartographic Boundary Shapefiles。由于您的问题与纽约市有关,我们可以使用NYU提供的shapefile。

使用纽约大学的GeoJSON文件:

代码语言:javascript
运行
复制
import geopandas
df = geopandas.read_file('nyu-2451-34561-geojson.json')

df[df.ntaname == 'Gramercy'].geometry
#54    (POLYGON ((-73.97849845639804 40.7367909565254...
#Name: geometry, dtype: object

这些几何图形是Shapely MultiPolygons的,因此您可以访问一些有用的属性,例如,您可以将所有内容都估计为矩形,而不是使用多重多边形:

代码语言:javascript
运行
复制
df[df.ntaname == 'Gramercy'].geometry.bounds
#         minx      miny       maxx       maxy
#54 -73.989967  40.73135 -73.978052  40.743325

# Or get the full boundary:
df[df.ntaname == 'Gramercy'].geometry.boundary
#54    (LINESTRING (-73.97849845639804 40.73679095652...
#dtype: object

如果在jupyter笔记本上工作,你可以对它的形状有一个感觉:

代码语言:javascript
运行
复制
df[df.ntaname == 'Gramercy'].geometry.item()

票数 0
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/55319747

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档