Chiếu một Point
Ở chương trước, bạn đã làm việc với vị trí của Tháp Eiffel. Một lần nữa, chúng tôi cung cấp cho bạn các tọa độ trong một hệ tọa độ phẳng (projected), để bạn có thể, chẳng hạn, tính khoảng cách. Giờ hãy quay lại địa danh mang tính biểu tượng này và biểu diễn vị trí của nó bằng tọa độ địa lý: 48°51′29.6″B, 2°17′40.2″Đ. Hoặc ở dạng thập phân: vĩ độ 48.8584 và kinh độ 2.2945.
Đối tượng hình học của Shapely không có khái niệm CRS, vì vậy không thể chuyển trực tiếp sang CRS khác. Do đó, chúng ta sẽ dùng GeoPandas để biến đổi điểm vị trí Tháp Eiffel sang một CRS khác. Ta sẽ đặt điểm đơn lẻ này vào một GeoSeries, dùng phương thức to_crs(), rồi trích xuất lại điểm.
GeoPandas đã được import sẵn.
Bài tập này là một phần của khóa học
Làm việc với Dữ liệu Không gian địa lý trong Python
Hướng dẫn bài tập
- Tạo một đối tượng point của shapely với tọa độ của Tháp Eiffel và gán vào biến
eiffel_tower. - Tạo một GeoSeries (gọi là
s_eiffel_tower) với Tháp Eiffel là phần tử duy nhất và chỉ định CRS là EPSG:4326. - Chuyển
s_eiffel_towersang EPSG:2154 và đặt kết quả làs_eiffel_tower_projected
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# 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)