포인트 투영하기
이전 장에서 에펠탑의 위치를 다뤘습니다. 그때는 거리를 계산하는 등 작업을 할 수 있도록 투영 좌표계에서 좌표를 제공했죠. 이제 이 상징적인 랜드마크로 돌아가, 지리 좌표(geographical coordinates)로 표현해 볼게요: 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)