在 R 中找到一条道路的长度。

huangapple go评论65阅读模式
英文:

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>

## Question

I 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&#39;t work as expected:  

    library(sf)
    df_sf &lt;- st_as_sf(df, coords = c(&quot;x&quot;, &quot;y2&quot;))
    
    &gt; 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 0


Is there a way I can convert the `df` to a polyline sf object and then find the length? Any other suggestions?

### Data

    df &lt;- 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 = &quot;data.frame&quot;)

  [1]: https://i.stack.imgur.com/3wb6i.png

</details>


# 答案1
**得分**: 4

要在`sf`中执行这个操作,你可以将数据转换成线串(linestring),然后使用`st_length`函数:
``` r
sf::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)))
#&gt; [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))
#&gt; [1] 25077.29

huangapple
  • 本文由 发表于 2023年2月23日 22:11:23
  • 转载请务必保留本文链接:https://go.coder-hub.com/75545954.html
匿名

发表评论

匿名网友

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen:

确定