R语言如何计算三角形周长

请问各位，已知二维平面上三个点的坐标，如下，我用了2种方法计算出三角形周长，
问题1：为什么两方法算出来不一样，perimeter函数是怎么算的？
问题2：如何用R快速实现三角形周长的计算？


prob11 <- -180
prob12 <- -20
prob21 <- -140
prob22 <- 55
prob31 <- 10
prob32 <- 0
#method1
dist12 <- sqrt((prob22-prob12)^2+(prob21-prob11)^2)
dist13 <- sqrt((prob32-prob12)^2+(prob31-prob11)^2)
dist23 <- sqrt((prob22-prob32)^2+(prob21-prob31)^2)
(dist123 <- dist12+dist13+dist23) #435.8152
#method2
library(geosphere)
xy <- rbind(c(prob11,prob12), c(prob21,prob22), c(prob31, prob32))
perimeter(xy)  #40005149

矩阵的生成按你给的顺序则如下，结果是一个3X2矩阵，一样的。

1. Tri1 <- matrix(c(-180,-20,-140,55,10,0), 3, byrow = T)

复制代码

第一种方法：计算的是直角坐标平面上的三角形周长，两点间距离是在坐标系上的距离。
第二种方法：用的是球面三角学的包，虽然你所给出的参数与上面一样，但在perimeter()函数中，是作为经纬度参与计算的，默认的a和f参数是地球的，而结果是地球表面上的三角周长，单位是米。

至于有没有更简单的计算周长的方法，第一种不是最简单最直接的吗？

llb_321 发表于 2022-10-20 14:51
第一种方法：计算的是直角坐标平面上的三角形周长，两点间距离是在坐标系上的距离。
第二种方法：用的是球 ...

谢谢回复！第二种方法明白，我的主要问题是第一种方法R有函数实现吗？（除了自己把那几行放到自编函数里）是不是有能提升一点速度的算法？在多次运算中还是有用的。

mengqinqing 发表于 2022-10-20 21:19
谢谢回复！第二种方法明白，我的主要问题是第一种方法R有函数实现吗？（除了自己把那几行放到自编函数里） ...

其实不用自编函数，试试下面的代码

1. Tri<-matrix(c(-180,-140,10,-20,55,0),3)
   
2. sum(sqrt(dist(Tri[,1])^2+dist(Tri[,2])^2))

复制代码

llb_321 发表于 2022-10-21 10:07
其实不用自编函数，试试下面的代码

谢谢回复！

mengqinqing 发表于 2022-10-24 11:59
谢谢回复！

呵呵，我之前说的代码其实也搞复杂了，你再试试这个，最简单的：

1. Tri <- matrix(c(-180,-20,-140,55,10,0), 3, byrow = T)
   
2. sum(dist(Tri))

复制代码