ฉันมีรูปร่างหลายแบบใน CRS ที่แตกต่างกัน (ส่วนใหญ่เป็น WGS84 lat / lon) ที่ฉันต้องการแปลงเป็นรูปแบบทั่วไป (น่าจะเป็น Albers Equal Area Conic แต่ฉันอาจขอความช่วยเหลือในการเลือกคำถามอื่นเมื่อปัญหาของฉันดีขึ้น -defined)
ฉันใช้เวลาสองสามเดือนในการทำสถิติเชิงพื้นที่ใน 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 แต่เพื่อเปลี่ยนพิกัดให้ตรงกัน ("reproject", "transform" ฯลฯ )
แก้ไข
ไม่รวม 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