计算三角形区域 - R.

时间:2018-01-15 20:23:02

标签: r

我试图用固定点(Cx,Cy)计算一对运行点(X1,Y1,... Xn,Yn)的三角形区域,但没有成功。有人可以告诉我,问题是什么,所以我可以尝试解决它?

脚本:

library(ggplot2)

nElem <- 100
xData <- as.data.frame(seq(1,nElem,5))
yData <- as.data.frame(seq(5,nElem,5))
xyDATA<- cbind(xData,yData)
colnames(xyDATA) <- c("xCoord","yCoord")

Cx <- 10
Cy <- 1

ggplot(xyDATA) + geom_point(aes(x = xCoord, y = yCoord)) + geom_point(aes(x = Cx, y = Cy),colour="red",size=4)

for(i in 1:19) 
{
  Ax <- xyDATA[i,1]
  Ay <- xyDATA[i,2]
  Bx <- xyDATA[i+1,1]
  By <- xyDATA[i+1,2]


  s <-  abs(0.5*((Ax*(By-Cy))+(Bx*(Cy-Ay))+(Cx*(Ay-By))))
  # print(Ax)
  # print(Ay)
  # print(Bx)
  # print(By)
  print(s)
}

1 个答案:

答案 0 :(得分:0)

如果您没有看到绘制的点图,则应修改ggplot行,如:

    p <- ggplot(xyDATA) + geom_point(aes(x = xCoord, y = yCoord)) +
 geom_point(aes(x = Cx, y = Cy),colour="red",size=4)

    print(p)

图形对象的print()方法生成实际显示。

现在图像显示为: enter image description here