楼主: RRRUP
672 1

[程序分享] R语言循环 [推广有奖]

  • 0关注
  • 0粉丝

等待验证会员

大专生

10%

还不是VIP/贵宾

-

威望
0
论坛币
2 个
通用积分
1.0004
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
157 点
帖子
14
精华
0
在线时间
62 小时
注册时间
2018-12-20
最后登录
2020-6-6

相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
求助各位!现在手中有一组代码是计算不规则多边形图心的,但是他只针对一组数据,现在有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

二维码

扫码加我 拉你入群

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

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

关键词:不规则多边形 多边形图 多边形 不规则 组数据

帮你顶顶帖子让更多人看到哈~~

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-5-23 00:42