지리 좌표가도 분 초 형식이지만 다음과 같은 구분 기호가없는 CSV 장소가 있습니다.
df <- data.frame(name = c("farm_1", "farm_2", "seabrook_1", "rocks_road"),
lat = c(425319.3, 425317, 425317.1, 425323.3),
long = c(705045.5, 705101.1, 705145.4, 705219.8))
name long lat
farm_1 425319.3 705045.5
farm_2 425317 705101.1
seabrook_1 425317.1 705145.4
rocks_road 425323.3 705219.8
다음과 같이 지리 좌표가도 분 분 형식으로 된 또 다른 CSV 장소가 있습니다.
df_2 <- data.frame(name = c("exeter_road", "hampton_hill", "portsmouth_ave", "pebble_ln"),
GPS_cordinates_DMM = c("N42 58.855 W70 56.473", "N42 58.666 W70 54.981",
"N42 56.579 W70 52.550", "N42 55.949 W70 53.631"))
name GPS_cordinates_DMM
exeter_road N42 58.855 W70 56.473
hampton_hill N42 58.666 W70 54.981
portsmouth_ave N42 56.579 W70 52.550
pebble_ln N42 55.949 W70 53.631
각 데이터 프레임의 좌표를 구문 분석하여 십진법 위도와 경도로 변환하고 싶습니다. 예를 들어 첫 번째 데이터 프레임은 다음과 같습니다.
df_dec <- data.frame(name = c("farm_1", "farm_2", "seabrook_1", "rocks_road"),
latitude = c(42.88869444, 42.88805556, 42.88808333, 42.88980556),
longitude = c(70.84597222, 70.85030556, 70.86261111, 70.87216667))
name latitude longitude
farm_1 42.88869 70.84597
farm_2 42.88806 70.85031
seabrook_1 42.88808 70.86261
rocks_road 42.88981 70.87217
두 번째 데이터 프레임은 다음과 같습니다.
df_2_dec <- df_2 <- data.frame(name = c("exeter_road", "hampton_hill", "portsmouth_ave", "pebble_ln"),
latitude = c(42.98091667, 42.97776667, 42.94298333, 42.93248333),
longitude = c(70.94121667, 70.91635, 70.87583333, 70.89385))
name latitude longitude
exeter_road 42.98092 70.94122
hampton_hill 42.97777 70.91635
portsmouth_ave. 42.94298 70.87583
pebble_ln 42.93248 70.89385
그런 다음 결국 결합하고 매핑 / 분석 할 수 있습니다.
이러한 좌표 유형을 구문 분석하고 변환 할 수있는 패키지 또는 기능이 있습니까?
그렇지 않다면 첫 번째 데이터 세트의 두 번째 행 위도에 소수가없는 것과 같은 문제를 처리 할 수있는 강력하고 강력한 문서를 작성하는 것이 좋습니다.
를 사용 substr하면 위치에 따라 문자열에서도, 분, 초의 숫자 값을 긁어 내고 ( substring종료 위치가 필요하지 않음) 숫자로 바꾸고 계산할 수 있습니다.
f1 <- function(x) (as.numeric(substr(x, 1, 2))*60^2 + as.numeric(substr(x, 3, 4))*60 +
as.numeric(substring(x, 5)))/60^2
res1 <- data.frame(name=df$name, lapply(df[-1], f1))
res1
# name lat long
# 1 farm_1 42.88869 70.84597
# 2 farm_2 42.88806 70.85031
# 3 seabrook_1 42.88808 70.86261
# 4 rocks_road 42.88981 70.87217
두 번째 표본은 N, S, E, W에서 분할 할 수 있습니다. strsplit기본적으로 첫 번째 표본 과 동일한 작업을 수행합니다.
tmp <- as.data.frame(
gsub("\\D", "", do.call(rbind, strsplit(df_2$GPS_cordinates_DMM, "[NSEW]"))[,-1]))
f2 <- function(x) as.numeric(substr(x, 1, 2)) +
as.numeric(substring(x, 3))/1e3/60
res2 <- data.frame(name=df_2$name, setNames(lapply(tmp, f2), c("lat", "lon")))
res2
# name lat lon
# 1 exeter_road 42.98092 70.94122
# 2 hampton_hill 42.97777 70.91635
# 3 portsmouth_ave 42.94298 70.87583
# 4 pebble_ln 42.93248 70.89385
이 기사는 인터넷에서 수집됩니다. 재 인쇄 할 때 출처를 알려주십시오.
침해가 발생한 경우 연락 주시기 바랍니다[email protected] 삭제
몇 마디 만하겠습니다