投影一个点
在上一章中,我们处理了埃菲尔铁塔的位置。我们再次提供了投影坐标系中的坐标,便于您计算距离等操作。现在回到这个地标,并用地理坐标来表示它的位置:48°51′29.6″N,2°17′40.2″E。或者用十进制度:纬度 48.8584,经度 2.2945。
Shapely 的几何对象本身没有 CRS 的概念,因此不能直接转换到另一个 CRS。为此,我们将使用 GeoPandas 将埃菲尔铁塔的点位置转换到另一种 CRS。做法是把这个单点放入 GeoSeries,使用 to_crs() 方法,然后再取出该点。
GeoPandas 已经导入。
本练习是课程的一部分
在 Python 中处理地理空间数据
练习说明
- 使用埃菲尔铁塔的坐标创建一个 shapely 点对象,并将其赋值给变量
eiffel_tower。 - 创建一个 GeoSeries(命名为
s_eiffel_tower),其唯一元素为该点,并将 CRS 指定为 EPSG:4326。 - 将
s_eiffel_tower转换为 EPSG:2154,并将结果命名为s_eiffel_tower_projected
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Construct a Point object for the Eiffel Tower
from shapely.geometry import Point
eiffel_tower = ____
# Put the point in a GeoSeries with the correct CRS
s_eiffel_tower = geopandas.GeoSeries([____], crs={'init': '____'})
# Convert to other CRS
s_eiffel_tower_projected = s_eiffel_tower.____
# Print the projected point
print(s_eiffel_tower_projected)