我试图找到沿特定邻域边界的每个点的纬度和经度。例如,我尝试使用geopy -当我输入一个特定的邻域时,它会返回一对坐标。我希望收到一个坐标列表,该列表将勾勒出指定邻域的轮廓。
具体地说,我正在努力获得纽约市曼哈顿的社区。
谢谢。
from geopy.geocoders import Nominatim
geolocator = Nominatim()
location = geolocator.geocode('Gramercy NYC')
print(location)
输出:位置(美国纽约市纽约市曼哈顿曼哈顿社区委员会6号,纽约市,纽约市,(40.7355189,-73.9840794,0.0))
我正在得到一组坐标,我想要多个坐标。
发布于 2019-03-24 17:16:01
geopy的Nominatim
类支持返回完整的几何数据,而不是单个点:
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
将包含以下内容:
'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))'
几何图形支持的格式为wkt
、svg
、kml
和geojson
。有关更多详细信息,请参阅the Nominatim.geocode docs。
您可能还想考虑先在https://nominatim.openstreetmap.org/页面上尝试查询,它提供了一个方便的web界面,可以显示这些几何图形。
发布于 2019-03-24 01:44:47
我不知道geopy
能提供这一点。
你真正需要的是你感兴趣的地理位置的shapefile。美国人口普查局为不同类型的地理位置提供了其中的几个:Cartographic Boundary Shapefiles。由于您的问题与纽约市有关,我们可以使用NYU提供的shapefile。
使用纽约大学的GeoJSON文件:
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的,因此您可以访问一些有用的属性,例如,您可以将所有内容都估计为矩形,而不是使用多重多边形:
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笔记本上工作,你可以对它的形状有一个感觉:
df[df.ntaname == 'Gramercy'].geometry.item()
https://stackoverflow.com/questions/55319747
复制相似问题