Well, I would like to say thanks Wolverine for all the effort you are putting into pulling this all together.
I'm not sure if I speak for the others in this thread, but I'm finding the challenge of working through the difficulty of getting the coords reliably together quite fun. It definitely adds a new dimension to the game for me, and ED wouldn't be ED without a bit of a challenge
I believe the result that comes back from Snuble's trilateration function is rounded to 5 decimal places.
To be honest I'm not sure why I used 6 decimal places on the averaged value.
I can't say I really understand all this 1/32 talk. Are you saying it would be better to use 3 decimal places from both the trilateration function and the final averaged output?
I meant I would check a few more in game. As there is a slight overlap with Harbringer's data I thought I would compare the two. I don't fancy doing any recalculations, already spending more time in PyCharm than ED, which is not good!
function trilateration3d($p1,$p2,$p3,$p4){
$ex = vector_unit(vector_diff($p2, $p1));
$i = vector_dot_product($ex, vector_diff($p3, $p1));
$ey = vector_unit(vector_diff(vector_diff($p3, $p1), vector_multiply($ex, $i)));
$ez = vector_cross($ex,$ey);
$d = vector_length($p2, $p1);
$r1 = $p1[3]; $r2 = $p2[3]; $r3 = $p3[3]; $r4 = $p4[3];
if($d - $r1 >= $r2 || $r2 >= $d + $r1){
return array();
}
$j = vector_dot_product($ey, vector_diff($p3, $p1));
$x = (($r1*$r1) - ($r2*$r2) + ($d*$d)) / (2*$d);
$y = ((($r1*$r1) - ($r3*$r3) + ($i*$i) + ($j*$j)) / (2*$j)) - (($i*$x) / $j);
$z = $r1*$r1 - $x*$x - $y*$y;
if($z < 0){
return array();
}
$z1 = sqrt($z);
$z2 = $z1 * -1;
$result1 = $p1;
$result1 = vector_sum($result1, vector_multiply($ex, $x));
$result1 = vector_sum($result1, vector_multiply($ey, $y));
$result1 = vector_sum($result1, vector_multiply($ez, $z1));
$result2 = $p1;
$result2 = vector_sum($result2, vector_multiply($ex, $x));
$result2 = vector_sum($result2, vector_multiply($ey, $y));
$result2 = vector_sum($result2, vector_multiply($ez, $z2));
$r1 = vector_length($p4, $result1);
$r2 = vector_length($p4, $result2);
$t1 = $r1 - $r4;
$t2 = $r2 - $r4;
$coords = array();
if(abs($t1) < abs($t2)){
$coords = array((round(($result1[0]*32),0)/32), (round(($result1[1]*32),0)/32), (round(($result1[2]*32),0)/32));
}
else{
$coords = array((round(($result2[0]*32),0)/32), (round(($result2[1]*32),0)/32), (round(($result2[2]*32),0)/32));
}
return $coords;
}
>trade.py local --ly 10 --dp 6 hermitage
Local systems to HERMITAGE within 10.0 ly.
------------------------------------------
4.291161 LHS 3281
5.093846 WOLF 654
7.196801 LP 275-68
7.709129 TILIAN
9.480758 LHS 3262
9.773859 ZETA HERCULIS
9.917876 LHS 6309function trilateration3d($p1,$p2,$p3,$p4){
$ex = vector_unit(vector_diff($p2, $p1));
$i = vector_dot_product($ex, vector_diff($p3, $p1));
$ey = vector_unit(vector_diff(vector_diff($p3, $p1), vector_multiply($ex, $i)));
$ez = vector_cross($ex,$ey);
$d = vector_length($p2, $p1);
$r1 = $p1[3]; $r2 = $p2[3]; $r3 = $p3[3]; $r4 = $p4[3];
if($d - $r1 >= $r2 || $r2 >= $d + $r1){
return array();
}
$j = vector_dot_product($ey, vector_diff($p3, $p1));
$x = (($r1*$r1) - ($r2*$r2) + ($d*$d)) / (2*$d);
$y = ((($r1*$r1) - ($r3*$r3) + ($i*$i) + ($j*$j)) / (2*$j)) - (($i*$x) / $j);
$z = $r1*$r1 - $x*$x - $y*$y;
if($z < 0){
return array();
}
$z1 = sqrt($z);
$z2 = $z1 * -1;
$result1 = $p1;
$result1 = vector_sum($result1, vector_multiply($ex, $x));
$result1 = vector_sum($result1, vector_multiply($ey, $y));
$result1 = vector_sum($result1, vector_multiply($ez, $z1));
$result2 = $p1;
$result2 = vector_sum($result2, vector_multiply($ex, $x));
$result2 = vector_sum($result2, vector_multiply($ey, $y));
$result2 = vector_sum($result2, vector_multiply($ez, $z2));
$r1 = vector_length($p4, $result1);
$r2 = vector_length($p4, $result2);
$t1 = $r1 - $r4;
$t2 = $r2 - $r4;
$coords = array();
if(abs($t1) < abs($t2)){
$result1[0]+=(1/64);
$result1[0]*=32;
$result1[0]=floor($result1[0]);
$result1[0]/=32;
$result1[1]+=(1/64);
$result1[1]*=32;
$result1[1]=floor($result1[1]);
$result1[1]/=32;
$result1[2]+=(1/64);
$result1[2]*=32;
$result1[2]=floor($result1[2]);
$result1[2]/=32;
$coords = array($result1[0], $result1[1], $result1[2]);
}
else{
$result2[0]+=(1/64);
$result2[0]*=32;
$result2[0]=floor($result2[0]);
$result2[0]/=32;
$result2[1]+=(1/64);
$result2[1]*=32;
$result2[1]=floor($result2[1]);
$result2[1]/=32;
$result2[2]+=(1/64);
$result2[2]*=32;
$result2[2]=floor($result2[2]);
$result2[2]/=32;
$coords = array($result2[0], $result2[1], $result2[2]);
}
return $coords;
}I've looked over my code and made an error.
After fixing the error your alternate method gives the same result as the previous method.
Does that mean the coordinates you posted earlier need recalculating?I've looked over my code and made an error.
After fixing the error your alternate method gives the same result as the previous method.