英文:
Find the length of a road in R
问题
以下是您提供的问题的翻译部分：
我有一个包含假想道路的笛卡尔坐标的数据框`df`。我想要找到这条道路的长度。下面是一个绘图：library(ggplot2)ggplot(df) +geom_point(aes(x, y2))[![enter image description here][1]][1]## 我尝试找到长度的方法：我尝试将`df`转换为一个sf对象，然后找到长度，但它没有按预期工作：library(sf)df_sf <- st_as_sf(df, coords = c("x", "y2"))> sf::st_length(df_sf)[1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0[36] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0[71] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0[106] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0是否有一种方法可以将`df`转换为线性sf对象，然后找到长度？还有其他建议吗？### 数据df <- structure(list(x = c(12158, 12158, 12158, 12158, 12158, 12158,12158, 12158, 12158, 12158, 12157.999023, 12158, 12158, 12156.886719,12154.849609, 12151.196289, 12147.6875, 12143.267578, 12137.9375,12131.700195, 12124.556641, 12116.507813, 12107.557617, 12097.708984,12086.962891, 12075.326172, 12062.798828, 12049.386719, 12035.092773,12019.921875, 12003.879883, 11986.969727, 11969.197266, 11950.568359,11931.087891, 11910.761719, 11889.597656, 11867.599609, 11844.777344,11821.134766, 11796.681641, 11771.423828, 11745.369141, 11718.525391,11690.901367, 11662.505859, 11633.345703, 11603.432617, 11572.773438,11541.37793, 11509.255859, 11476.417969, 11442.87207, 11408.630859,11373.703125, 11338.099609, 11301.831055, 11264.910156, 11227.345703,11189.151367, 11150.337891, 11110.916016, 11070.899414, 11030.299805,10989.128906, 10947.400391, 10905.125, 10862.318359, 10818.991211,10794.041016, 10746.004883, 10685.342773, 10493.634766, 10488.304688,9822.144531, 9590.665039, 9495.904297, 9393.224609, 9333.71875,9274.916016, 9216.833008, 9159.488281, 9102.900391, 9047.083984,8992.058594, 8937.839844, 8884.444336, 8831.888672, 8780.1875,8729.358398, 8679.416016, 8630.375, 8582.251953, 8535.05957,8488.813477, 8443.52832, 8399.21582, 8355.890625, 8313.566406,8272.255859, 8231.970703, 8192.724609, 8154.52832, 8117.394043,8081.33252, 8046.355469, 8012.473633, 7979.696289, 7948.03418,7917.49707, 7888.09375, 7859.833984, 7832.725586, 7806.776855,7781.996094, 7758.390625, 7735.967773, 7714.734375, 7694.696777,7675.861328, 7658.23291, 7641.818359, 7626.621094, 7612.646973,7599.899414, 7588.383301, 7578.101074, 7569.056152, 7561.251953,7554.689941, 7548.160156, 7544.127441, 7540.84082, 7538, 7538,7538), y2 = c(15255, 14520, 13200, 11880, 10560, 9240, 7920,6600, 5280, 3960, 2786.666992, 2640, 1916.935059, 1725.494141,<details><summary>英文:</summary>## QuestionI have a dataframe `df` which contains cartesian coordinates of a hypothetical road. I want to find the length of the road. Following is a plot:library(ggplot2)ggplot(df) +geom_point(aes(x, y2))[![enter image description here][1]][1]## My attempt to find the length:I tried converting `df` to an sf object and then find length, but it didn't work as expected:library(sf)df_sf <- st_as_sf(df, coords = c("x", "y2"))> sf::st_length(df_sf)[1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0[36] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0[71] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0[106] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0Is there a way I can convert the `df` to a polyline sf object and then find the length? Any other suggestions?### Datadf <- structure(list(x = c(12158, 12158, 12158, 12158, 12158, 12158,12158, 12158, 12158, 12158, 12157.999023, 12158, 12158, 12156.886719,12154.849609, 12151.196289, 12147.6875, 12143.267578, 12137.9375,12131.700195, 12124.556641, 12116.507813, 12107.557617, 12097.708984,12086.962891, 12075.326172, 12062.798828, 12049.386719, 12035.092773,12019.921875, 12003.879883, 11986.969727, 11969.197266, 11950.568359,11931.087891, 11910.761719, 11889.597656, 11867.599609, 11844.777344,11821.134766, 11796.681641, 11771.423828, 11745.369141, 11718.525391,11690.901367, 11662.505859, 11633.345703, 11603.432617, 11572.773438,11541.37793, 11509.255859, 11476.417969, 11442.87207, 11408.630859,11373.703125, 11338.099609, 11301.831055, 11264.910156, 11227.345703,11189.151367, 11150.337891, 11110.916016, 11070.899414, 11030.299805,10989.128906, 10947.400391, 10905.125, 10862.318359, 10818.991211,10794.041016, 10746.004883, 10685.342773, 10493.634766, 10488.304688,9822.144531, 9590.665039, 9495.904297, 9393.224609, 9333.71875,9274.916016, 9216.833008, 9159.488281, 9102.900391, 9047.083984,8992.058594, 8937.839844, 8884.444336, 8831.888672, 8780.1875,8729.358398, 8679.416016, 8630.375, 8582.251953, 8535.05957,8488.813477, 8443.52832, 8399.21582, 8355.890625, 8313.566406,8272.255859, 8231.970703, 8192.724609, 8154.52832, 8117.394043,8081.33252, 8046.355469, 8012.473633, 7979.696289, 7948.03418,7917.49707, 7888.09375, 7859.833984, 7832.725586, 7806.776855,7781.996094, 7758.390625, 7735.967773, 7714.734375, 7694.696777,7675.861328, 7658.23291, 7641.818359, 7626.621094, 7612.646973,7599.899414, 7588.383301, 7578.101074, 7569.056152, 7561.251953,7554.689941, 7548.160156, 7544.127441, 7540.84082, 7538, 7538,7538), y2 = c(15255, 14520, 13200, 11880, 10560, 9240, 7920,6600, 5280, 3960, 2786.666992, 2640, 1916.935059, 1725.494141,1646.221191, 1567.054199, 1514.8125, 1462.640137, 1410.552734,1358.566406, 1306.697266, 1254.959961, 1203.371582, 1151.947266,1100.702148, 1049.652832, 998.813965, 948.202148, 897.831543,847.718262, 797.877441, 748.32373, 699.073242, 650.140137, 601.539551,553.286621, 505.395508, 457.881348, 410.758301, 364.040527, 317.742676,271.878418, 226.462402, 181.507324, 137.02832, 93.037598, 49.549316,6.57666, -35.867188, -77.769531, -119.118164, -159.899414, -200.101074,-239.711914, -278.71875, -317.109863, -354.874023, -391.999023,-428.474609, -464.288574, -499.430664, -533.890137, -567.65625,-600.719238, -633.067871, -664.693848, -695.585938, -725.736328,-755.134277, -771.569336, -802.256348, -839.523926, -951.065918,-954.111816, -1334.774658, -1468.841064, -1526.266846, -1591.573242,-1631.31665, -1672.092529, -1713.888184, -1756.691406, -1800.489014,-1845.267334, -1891.012939, -1937.712158, -1985.350342, -2033.913086,-2083.385742, -2133.753174, -2185, -2237.110596, -2290.069336,-2343.859619, -2398.465576, -2453.870117, -2510.056641, -2567.007813,-2624.706787, -2683.135254, -2742.276123, -2802.11084, -2862.621338,-2923.789307, -2985.596191, -3048.022949, -3111.050537, -3174.659912,-3238.831543, -3303.546143, -3368.783691, -3434.524414, -3500.748535,-3567.435547, -3634.56543, -3702.117188, -3770.071045, -3838.405762,-3907.100586, -3976.134766, -4045.486816, -4115.13623, -4185.061523,-4255.241211, -4325.654297, -4396.279297, -4467.094238, -4538.078125,-4609.208496, -4680.464844, -4770.95166, -4848.422363, -4943.858398,-5270.765625, -5280, -7920)), row.names = c(NA, -136L), class = "data.frame")[1]: https://i.stack.imgur.com/3wb6i.png</details># 答案1**得分**: 4要在`sf`中执行这个操作，你可以将数据转换成线串（linestring），然后使用`st_length`函数：``` rsf::st_length(sf::st_linestring(as.matrix(df)))#> [1] 25077.29
不过需要注意的是，在这里你其实不需要使用sf，因为你可以通过基本的欧几里德距离计算得到相同的结果：
sum(sqrt(diff(df$x)^2 + diff(df$y2)^2))#> [1] 25077.29
英文:
To do this in sf, you can convert into a linestring and use st_length
sf::st_length(sf::st_linestring(as.matrix(df)))#> [1] 25077.29
Note though that you don't need sf at all here, as you get the same result using a basic Euclidean distance calculation:
sum(sqrt(diff(df$x)^2 + diff(df$y2)^2))#> [1] 25077.29
通过集体智慧和协作来改善编程学习和解决问题的方式。致力于成为全球开发者共同参与的知识库，让每个人都能够通过互相帮助和分享经验来进步。
评论