如何计算R中xy点的周长？

Question

我正在计算 xy 点的周长（包括按顺序排列的所有点），其中第一个点和最后一个点并不总是相等。我无法弄清楚为什么下面用于计算周长的代码不起作用。

理想情况下，我想要一种在第一个和最后一个 xy 点不同时进行标记的方法。

谢谢

  m
       x       y
 [1,]  606.3   95.4
 [2,]  612.4  178.7
 [3,]  610.2  222.6
 [4,]  610.2  222.8
 [5,]  625.8  249.8
 [6,]  625.8  250.1
 [7,]  633.9  268.9
 [8,]  668.7  272.2
 [9,]  693.7  222.6
[10,]  723.2  157.0
[11,]  738.6  109.9
[12,]  681.2   90.5
[13,]  606.3   95.4
[14,]  833.3  154.6
[15,]  753.7  267.5
[16,]  747.8  305.1
[17,]  773.8  354.7
[18,]  767.0  393.8
[19,]  763.0  442.0
[20,]  817.4  446.9
[21,]  817.6  446.9
[22,]  840.2  412.3
[23,]  892.1  317.7
[24,]  875.3  218.8
[25,]  833.3  154.6


library(geosphere)
perimeter(m)


**Error in perimeter(m) : could not find function "perimeter"**

Answer 1

这样就可以了..

样本数据

library( data.table )
m <- fread("x       y
  606.3   95.4
  612.4  178.7
  610.2  222.6
  610.2  222.8
  625.8  249.8
  625.8  250.1
  633.9  268.9
  668.7  272.2
  693.7  222.6
  723.2  157.0
  738.6  109.9
  681.2   90.5
  606.3   95.4
  833.3  154.6
  753.7  267.5
  747.8  305.1
  773.8  354.7
  767.0  393.8
  763.0  442.0
  817.4  446.9
  817.6  446.9
  840.2  412.3
  892.1  317.7
  875.3  218.8
  833.3  154.6")

代码

library( grDevices )

hull <- chull( m ) #create coordinates of convex hull
coords <- m[ c( hull, hull[1]), ]  # create a closed perimeter-polygon

#plot it
plot(m)        #plot points
lines(coords, col="red") #plot perimeter