서로 다른 CRS/투영에서 데이터 병합하기

모든 공간 객체에는 좌표 참조 체계(CRS)가 연결되어 있어요. 일반적으로 데이터 가져오기 시 설정되며, 공간 파일에서 직접 읽어 옵니다. neighborhoods 와 nyc_tracts 의 좌표계 정보도 이렇게 얻었어요.

sp 와 raster 패키지에는 객체의 CRS를 반환하는 proj4string() 함수가 있어요.

서로 다른 CRS를 사용하는 공간 데이터를 함께 다루는 것은, 한 데이터셋은 마일, 다른 데이터셋은 킬로미터 단위를 쓰는 상황과 비슷해요. 같은 것을 측정하지만 숫자는 바로 비교할 수 없죠.

이제 두 개의 폴리곤 객체를 살펴보겠습니다.

neighborhoods 에 proj4string() 을 호출하고, 이어서 nyc_tracts 에도 호출하세요. 두 문자열이 서로 다른지 확인하세요.
neighborhoods 의 coordinates() 에 head() 를 적용해 보고, nyc_tracts 에도 반복하세요. 문제가 보이나요? nyc_tracts 의 x 좌표는 대략 -70 주변인데, neighborhoods 는 약 1,000,000 부근이에요!
먼저 neighborhoods 를 그린 뒤, nyc_tracts 를 col = "red", add = TRUE 로 위에 겹쳐 그리세요.