全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
966 1
2018-12-20
求助各位!现在手中有一组代码是计算不规则多边形图心的,但是他只针对一组数据,现在有1000组这样的数据,如何加个循环能一次求出这些数据的图心?
options(digits = 3)
# data1 include distance and angle
data1$x <- data1$距离*sin(data1$角度/180*pi)
data1$y <- data1$距离*cos(data1$角度/180*pi)
#
# data1$x <- data1$X
# data1$y <- data1$Y
# data <- data1
order(data1$角度,decreasing = T)
data <- data1[order(data1$角度,decreasing = F),]
n <- nrow(data)
z1 <- numeric(n-2);z2 <- numeric(n-2);ss <- numeric(n-2)
for (i in 1:(n-2)) {
  d1 <- data[c(1,i+1,i+2),c(3,4)]
  zx <- mean(d1$x);zy <- mean(d1$y)
  s <- 1/2*(d1$x[1]*d1$y[2]+d1$x[2]*d1$y[3]+d1$x[3]*d1$y[1]-d1$x[1]*d1$y[3]-d1$x[2]*d1$y[1]-d1$x[3]*d1$y[2])
  z1 <- zx
  z2 <- zy
  ss <- s
}
zzx <- sum(z1*ss)/sum(ss)
zzy <- sum(z2*ss)/sum(ss)
zhongxin <- c(zzx,zzy)
zhongxin
QQ截图20181220141633.jpg

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2018-12-30 15:52:55
帮你顶顶帖子让更多人看到哈~~
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

分享

扫码加好友,拉您进群