如何获得椭圆的半轴长度？在R [英] How to obtain the lengths of semi axes of an ellipse? in R

问题描述

我有这组x和y坐标：

  x <-c（1.798805,2.402390,2.000000,3.000000， 1.000000）
 y <-c（0.3130147,0.4739707,0.2000000,0.8000000,0.1000000）
 as.matrix（cbind（x，y）） - > d 
  

我想计算包含这组点的椭球，我使用函数 ellipsoidhull（） $ b

 > $ c $在包集群中，我得到：椭圆体（d）
'椭圆体'2维：`
 center =（2.00108 0.36696）;平方ave.radius d ^ 2 = 2` 
和形状矩阵= 
x 0.66590 0.233106 $ b $ 0.23311 0.095482 
因此，面积= 0.60406 
  

然而，从这些结果我可以得出这个椭圆的半长轴的长度并不明显。

有什么想法？

非常感谢您提前。

Tina 。

解决方案

半轴的平方是形状矩阵的
特征值，乘以平均半径的平方

pre $ x < - c（1.798805,2.402390,2.000000,3.000000,1.000000）
y< - c（ 0.3130147,0.4739707,0.2000000,0.8000000,0.1000000）
d < - cbind（x，y）
library（cluster）
r < - ellipsoidhull（d）
plot（x， y，asp = 1，xlim = c（0,4））
lines（predict（r））
e <-str（eigen（r $ cov）$ values）
a< - sqrt（r $ d2）* e [1]＃半长轴
b < - sqrt（r $ d2）* e [2]＃半短轴
theta < - seq（ 0 ，2 * pi，length = 200）
lines（r $ loc [1] + a * cos（theta），r $ loc [2] + a * sin（theta））
lines（r $ loc [1] + b * cos（theta），r $ loc [2] + b * sin（theta））


I have this set of x and y coordinates:

x<-c(1.798805,2.402390,2.000000,3.000000,1.000000)
y<-c(0.3130147,0.4739707,0.2000000,0.8000000,0.1000000)
as.matrix(cbind(x,y))->d

and I want to calculate the ellipsoid that contains this set of points, I use the function ellipsoidhull() in the package "cluster", and I get:

> ellipsoidhull(d)
'ellipsoid' in 2 dimensions:`
 center = ( 2.00108 0.36696 ); squared ave.radius d^2 =  2`
and shape matrix =
x 0.66590 0.233106
y 0.23311 0.095482
  hence, area  =  0.60406

However it's not obvious to me how I can get from these results, the lengths of the semi-major axes of this ellipse.

Any idea?

Thank you very much in advance.

Tina.

解决方案

The square of the semi-axes are the eigenvalues of the shape matrix, times the average squared radius.

x <- c(1.798805,2.402390,2.000000,3.000000,1.000000)
y <- c(0.3130147,0.4739707,0.2000000,0.8000000,0.1000000)
d <- cbind( x, y )
library(cluster)
r <- ellipsoidhull(d)
plot( x, y, asp=1, xlim=c(0,4) )
lines( predict(r) )
e <- sqrt(eigen(r$cov)$values)
a <- sqrt(r$d2) * e[1]  # semi-major axis
b <- sqrt(r$d2) * e[2]  # semi-minor axis
theta <- seq(0, 2*pi, length=200)
lines( r$loc[1] + a * cos(theta), r$loc[2] + a * sin(theta) )
lines( r$loc[1] + b * cos(theta), r$loc[2] + b * sin(theta) )

这篇关于如何获得椭圆的半轴长度？在R的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！