서로 다른 CRS/투영에서 데이터 병합하기
모든 공간 객체에는 좌표 참조 체계(CRS)가 연결되어 있어요. 일반적으로 데이터 가져오기 시 설정되며, 공간 파일에서 직접 읽어 옵니다. neighborhoods 와 nyc_tracts 의 좌표계 정보도 이렇게 얻었어요.
sp 와 raster 패키지에는 객체의 CRS를 반환하는 proj4string() 함수가 있어요.
서로 다른 CRS를 사용하는 공간 데이터를 함께 다루는 것은, 한 데이터셋은 마일, 다른 데이터셋은 킬로미터 단위를 쓰는 상황과 비슷해요. 같은 것을 측정하지만 숫자는 바로 비교할 수 없죠.
이제 두 개의 폴리곤 객체를 살펴보겠습니다.
이 연습은 강의의 일부입니다
R로 지리공간 데이터 시각화하기
연습 안내
neighborhoods에proj4string()을 호출하고, 이어서nyc_tracts에도 호출하세요. 두 문자열이 서로 다른지 확인하세요.neighborhoods의coordinates()에head()를 적용해 보고,nyc_tracts에도 반복하세요. 문제가 보이나요?nyc_tracts의 x 좌표는 대략 -70 주변인데,neighborhoods는 약 1,000,000 부근이에요!- 먼저
neighborhoods를 그린 뒤,nyc_tracts를col = "red",add = TRUE로 위에 겹쳐 그리세요.
실습형 인터랙티브 연습
이 예제를 이 샘플 코드를 완성하여 풀어보세요.
library(sp)
# proj4string() on nyc_tracts and neighborhoods
# coordinates() on nyc_tracts and neighborhoods
# plot() neighborhoods and nyc_tracts