开始使用免费开始使用

投影一个点

在上一章中,我们处理了埃菲尔铁塔的位置。我们再次提供了投影坐标系中的坐标,便于您计算距离等操作。现在回到这个地标,并用地理坐标来表示它的位置: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)
编辑并运行代码