從座標集精確定位中心點

例如從組意見：

point xcoord ycoord accuracy time 
1  x1  y1  2000  t1 
2  x2  y2  2000  t2 
3  x3  y3  2000  t3 
n  xn  yn  2000  tn

Techically所有這些點都與他們的準確性半徑的圓中心點，我們知道，我們是在時間ATLEAST一個圓內。

我想計算圓心相交區域的中心點和精度。

這是一個粗略的例子。綠線是試圖代表準確性 crude example picture

來源

2013-10-11 warbaque

+1非常有趣的問題......這是我在C++

全局表和變量的方法來存儲所有需要：

const int N=5;  // number of measurement circles 
struct _circle { double x,y,r; _circle() { x=0; y=0; r=0.0; } } circle[N],circle0,circle1;

circle[] - 測量位置
circle0 - 原始實際位置（未知）
circle1 - 計算更精確的位置

這裏是我的隨機測量（複製您的測量值代替）：

int i; 
double r,a,x0,x1,y0,y1; 
// set real position 
circle0.x=50.0; 
circle0.y=50.0; 
circle0.r=25.0; 
// set meassured positions 
Randomize(); 
for (i=0;i<N;i++) 
    { 
    r=circle0.r*(0.2+(0.3*double(Random(101))*0.01)); 
    a= 2.0*M_PI*   double(Random(101))*0.01; 
    circle[i]=circle0; 
    circle[i].x+=r*cos(a); 
    circle[i].y+=r*sin(a); 
    }

這是如何計算平均作風位置：

// compute more accurate position (average style) 
circle1.x=0.0; 
circle1.y=0.0; 
circle1.r=circle0.r; 
for (i=0;i<N;i++) 
    { 
    circle1.x+=circle[i].x; 
    circle1.y+=circle[i].y; 
    } 
circle1.x/=N; 
circle1.y/=N; 
for (i=0;i<N;i++) 
    { 
    r=circle1.x-circle[i].x; r*=r; 
    a=circle1.y-circle[i].y; a*=a; 
    r=circle[i].r-sqrt(a+r); 
// if (circle1.r>r) circle1.r=r; // best accuracy 
    if (circle1.r<r) circle1.r=r; // worse accuracy 
    }

在兩個最後的ifs中選擇了你想要的精度半徑...
po sition是基於平均所有實測

這是如何計算幾何風格位置：

// compute more accurate position (geometry style) 
x0=circle[i].x; x1=x0; 
y0=circle[i].y; y1=y0; 
for (i=0;i<N;i++) 
    { 
    a=circle[i].x; if (x0>a) x0=a; if (x1<a) x1=a; 
    a=circle[i].y; if (y0>a) y0=a; if (y1<a) y1=a; 
    } 
circle1.x=0.5*(x0+x1); x1-=x0; x1*=x1; 
circle1.y=0.5*(y0+y1); y1-=y0; y1*=y1; 
circle1.r=1.0*sqrt(x1+y1);

位置佔地面積

的中心，這裏有一些預覽我的代碼輸出：

accuracy overlap

來源

2013-12-28 17:00:34 Spektre

從座標集精確定位中心點

回答

相關問題