Question

我正在尝试使用指定的测量包转换一些数据，但我没有成功。

我的数据：

Long        Lat
62ᵒ36.080   58ᵒ52.940
61ᵒ28.020   54ᵒ59.940
62ᵒ07.571   56ᵒ48.873
62ᵒ04.929   57ᵒ33.605
63ᵒ01.419   60ᵒ30.349
63ᵒ09.555   61ᵒ29.199
63ᵒ43.499   61ᵒ23.590
64ᵒ34.175   62ᵒ30.304
63ᵒ16.342   59ᵒ16.437
60ᵒ55.090   54ᵒ49.269
61ᵒ28.013   54ᵒ59.928
62ᵒ07.868   56ᵒ48.040
62ᵒ04.719   57ᵒ32.120
62ᵒ36.083   58ᵒ51.766
63ᵒ01.644   60ᵒ30.714
64ᵒ33.897   62ᵒ30.772
63ᵒ43.604   61ᵒ23.426
63ᵒ09.288   61ᵒ29.888
63ᵒ16.722   59ᵒ16.204

我正在尝试：

library(measurements)

library(readxl)

coord = read.table('coord_converter.txt', header = T, stringsAsFactors = F) 

# change the degree symbol to a space
lat = gsub('°','', coord$Lat)
long = gsub('°','', coord$Long)

# convert from decimal minutes to decimal degrees
lat = measurements::conv_unit(lat, from = 'deg_dec_min', to = 'dec_deg')
long = measurements::conv_unit(long, from = 'deg_dec_min', to = 'dec_deg')

我在倒数第二行得到的是：

警告讯息：

在split(as.numeric(unlist(strsplit(x, " "))) * c(3600, 60), f = rep(1:length(x)，：NAs introduced by coercion
在as.numeric(unlist(strsplit(x, " "))) * c(3600, 60)：longer object length is not a multiple of shorter object length
在split.default(as.numeric(unlist(strsplit(x, " "))) * c(3600,：data length is not a multiple of split variable

有人可以指出我的错误或提出如何进行的建议吗？

谢谢！

Answer 1

我认为这里的问题是gsub调用后，measurements::conv_unit所要求的度数和分钟不是以空格分隔的。

例如，这样可以正常工作（对于这个可重复的例子，我也将“ᵒ”改为“°”）：

library(measurements)

#read your data
txt <- 
"Long        Lat
62°36.080   58°52.940
61°28.020   54°59.940
62°07.571   56°48.873
62°04.929   57°33.605
63°01.419   60°30.349
63°09.555   61°29.199
63°43.499   61°23.590
64°34.175   62°30.304
63°16.342   59°16.437
60°55.090   54°49.269
61°28.013   54°59.928
62°07.868   56°48.040
62°04.719   57°32.120
62°36.083   58°51.766
63°01.644   60°30.714
64°33.897   62°30.772
63°43.604   61°23.426
63°09.288   61°29.888
63°16.722   59°16.204"

coord <- read.table(text = foo, header = TRUE, stringsAsFactors = F)

# change the degree symbol to a space
lat = gsub('°',' ', coord$Lat)
long = gsub('°',' ', coord$Long)

# convert from decimal minutes to decimal degrees
lat = measurements::conv_unit(lat, from = 'deg_dec_min', to = 'dec_deg')
long = measurements::conv_unit(long, from = 'deg_dec_min', to = 'dec_deg')

...产量

> cbind(long, lat)
      long               lat               
 [1,] "62.6013333333333" "58.8823333333333"
 [2,] "61.467"           "54.999"          
 [3,] "62.1261833333333" "56.81455"        
 [4,] "62.08215"         "57.5600833333333"
 [5,] "63.02365"         "60.5058166666667"
 [6,] "63.15925"         "61.48665"        
 [7,] "63.7249833333333" "61.3931666666667"
 [8,] "64.5695833333333" "62.5050666666667"
 [9,] "63.2723666666667" "59.27395"        
[10,] "60.9181666666667" "54.82115"        
[11,] "61.4668833333333" "54.9988"         
[12,] "62.1311333333333" "56.8006666666667"
[13,] "62.07865"         "57.5353333333333"
[14,] "62.6013833333333" "58.8627666666667"
[15,] "63.0274"          "60.5119"         
[16,] "64.56495"         "62.5128666666667"
[17,] "63.7267333333333" "61.3904333333333"
[18,] "63.1548"          "61.4981333333333"
[19,] "63.2787"          "59.2700666666667"

使用测量包转换地理坐标

1 个答案: