R에서 sp 객체 투영

다른 CRS (대부분 WGS84 위도 / 경도)에 여러 가지 모양 파일이있어 일반적인 투영 (예 : Albers Equal Area Conic)으로 변환하고 싶지만 문제가 좋아지면 다른 질문을 선택하는 데 도움을 요청할 수 있습니다 -한정된).

R에서 공간 통계 자료를 작성하는 데 몇 달을 보냈지 만 5 년 전이었습니다. 내 인생 에서 한 투영에서 다른 투영으로 sp객체 (예 :) 를 변환하는 방법을 기억할 수 없습니다 SpatialPolygonsDataFrame.

예제 코드 :

P4S.latlon <- CRS("+proj=longlat +datum=WGS84")
hrr.shp <- readShapePoly("HRR_Bdry"), verbose=TRUE, proj4string=P4S.latlon) 
# Shapefile available at 
#   http://www.dartmouthatlas.org/downloads/geography/hrr_bdry.zip 
#   but you must rename all the filenames to have the same 
#   capitalization for it to work in R

이제 SpatialPolygonsDataFrame적절한 프로젝션 정보가 있지만 원하는 프로젝션으로 변환하고 싶습니다. 나는 이것에 대해 다소 직관적이지 않은 이름의 기능이 있다고 생각하지만, 그것이 무엇인지 기억할 수 없습니다.

CRS를 변경하는 것이 아니라 좌표를 일치시키기 위해 ( "재 투영", "변환"등) 변경하고 싶습니다.

편집하다

이 shapefile에 대해 멕시코에 성가신 AK / HI 제외 :

library(taRifx.geo)
hrr.shp <- 
  subset(hrr.shp, !(grepl( "AK-" , hrr.shp@data$HRRCITY ) |
                                     grepl( "HI-" , hrr.shp@data$HRRCITY )) )
proj4string(hrr.shp) <- P4S.latlon

spTransform()rgdal 에서 메소드를 사용할 수 있습니다. 예제를 사용하여 캔자스 (26978)의 경우 오브젝트를 NAD83으로 변환 할 수 있습니다.

library(rgdal)
library(maptools)

P4S.latlon <- CRS("+proj=longlat +datum=WGS84")
hrr.shp <- readShapePoly("HRR_Bdry", verbose=TRUE, proj4string=P4S.latlon)
plot(hrr.shp)

hrr.shp.2 <- spTransform(hrr.shp, CRS("+init=epsg:26978"))
plot(hrr.shp.2)

새 투영에 저장하려면 :

writePolyShape(hrr.shp.2, "HRR_Bdry_NAD83")

편집 : 또는 @Spacedman의 제안에 따라 (CRS 정보로 .prj 파일을 작성) :

writeOGR(hrr.shp.2, dsn = getwd(), layer = "HRR_Bdry_NAD83", driver="ESRI Shapefile")

어떤 CRS를 투사할지 확실하지 않은 경우 다음 게시물을 참조하십시오.

• R에서 shapefile을 읽기 위해 proj4string에 올바른 값을 선택 하시겠습니까?

데이터에 CRS가 없을 때 CRS를 정의 / 할당하려면 다음을 참조하십시오.

• R에서 셰이프 파일이 없을 때 CRS를 셰이프 파일에 할당?

sf-package가 도입 되었으므로 ( 여기서 소품 sf1 , sf2 , sf3 , sf4 및 마이그레이션 안내서를 참조하십시오 ) st_transform()벡터 데이터를 다시 투영 하는 데 사용할 수 있습니다 .

require(sf)

hrr_sf = st_read('HRR_Bdry.shp', stringsAsFactors = FALSE,
    crs = 4326) # has +proj=longlat +datum=WGS84
plot(hrr_sf)

hrr_sf2 = st_transform(hrr_sf, "+init=epsg:26978") # 1st option sp::CRS() not working/ needed
hrr_sf2 = st_transform(hrr_sf, 26978) # 2nd option - EPSG code as an integer
plot(hrr_sf2)

# don't think about doing this:
hrr_sf3 = st_read('HRR_Bdry.shp', stringsAsFactors = FALSE,
    crs = 26978)

# Output layer
st_write(hrr_sf2, dsn = getwd(), layer = "HRR_Bdry_NAD83", driver = "ESRI Shapefile")

sf는 앞으로 sp를 대체 할 것이며 sp와 비교할 때 단순성과 속도로 인해 몇 가지 장점을 가지고 있습니다.