R语言循环

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

从一固定点引出六条不同长度的线段，构成一个不规则的多边形，求这个不规则多边形的图心。以下为图例、数据及代码。   
代码：
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
结果：
         X坐标1.431466
     Y坐标1.913244这个代码指针对一组数据使用，但是手中有许多组这样的数据，想加一个循环，能够对各列进行循环操作，求得图心的X、Y坐标。求各位帮助！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：不规则多边形 不同长度 多边形 不规则 固定点