python中三个点是什么意思_python或R的点列表中三个点的半径

多亏了卡蒂娅,她完成了这项工作!

要计算三维空间中给定点列表的半径,我使用以下代码:rm(list=ls())

setwd("...")

## ask Katia @ stackoverflow

# distance between 2 points

eucl.dist

# area of triangle

area

v1 = x1 - x2

v2 = x1 - x3

0.5* sqrt( (v1[2]*v2[3] - v1[3]*v2[2])**2 +

(v1[3]*v2[1] - v1[1]*v2[3])**2 +

(v1[1]*v2[2] - v1[2]*v2[1])**2)

}

# function to calculate curvature (optional for this problem)

curvature

4 * area(x1, x2, x3) / (eucl.dist(x1,x2)*eucl.dist(x2,x3)*eucl.dist(x3,x1))

}

# function to calculate radius

radius

0.25 * (eucl.dist(x1,x2)*eucl.dist(x2,x3)*eucl.dist(x3,x1)) / area(x1, x2, x3)

}

# reading data

path3D

colnames(path3D)

# calculate vector with radii

rads

rads

#add header

colnames(rads)

#considering, that for the first value no radius could be calculated

rads

#...also for the last one

rads

#merge the data

path3D

#write to file

write.table(path3D, "....txt", sep=";",row.names=FALSE)

再次感谢你!代码很好!在

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值