多亏了卡蒂娅,她完成了这项工作!
要计算三维空间中给定点列表的半径,我使用以下代码: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)
再次感谢你!代码很好!在