Gabriele:

New implementation:

Precision: 0.9928, Recall: 0.9133, F1: 0.9514
Total true vertices: 150
Matched vertices: 137
Missed true vertices: 13
Total fitted vertices: 138
Fake vertices: 1
Mean signed residuals (x,y,z): [0.00038846 0.00010329 0.00066991]
Std  signed residuals (x,y,z): [0.00853282 0.00977771 0.0052135 ]
Mean absolute residuals (x,y,z): [0.0043629  0.00525042 0.00273546]
Std  absolute residuals (x,y,z): [0.00733391 0.00823679 0.00448273]

Missed true vertices:
  ROF         x         y         z  nContributors
 11.0 -0.000259  0.022460  -3.69227           10.0
 13.0 -0.000585  0.002898   8.07200           11.0
 30.0 -0.024436 -0.009381   7.03565            1.0
 35.0 -0.023708 -0.009344  -5.51584            1.0
 48.0 -0.010583 -0.008603  -3.05986           14.0
 75.0  0.011355 -0.002970   1.90074            1.0
 90.0  0.002929 -0.018336   2.66971            1.0
102.0  0.001209  0.020929  -4.71790            1.0
123.0 -0.001262 -0.002698   8.83149           22.0
128.0  0.006334 -0.001214  -5.15652           10.0
158.0 -0.022409 -0.003296  -7.26120            3.0
161.0  0.002261  0.001900 -10.54620            1.0
164.0 -0.008004  0.000373  -5.64786            8.0

Old implementation:

Precision: 1.0000, Recall: 0.9067, F1: 0.9510
Total true vertices: 150
Matched vertices: 136
Missed true vertices: 14
Total fitted vertices: 136
Fake vertices: 0
Mean signed residuals (x,y,z): [-7.04813868e-05  1.47534449e-04  3.28789706e-04]
Std  signed residuals (x,y,z): [0.00157584 0.00351721 0.00311071]
Mean absolute residuals (x,y,z): [0.00092318 0.00105691 0.00067765]
Std  absolute residuals (x,y,z): [0.00127659 0.00335669 0.00305332]

Missed true vertices:
  ROF         x         y         z  nContributors
 11.0 -0.000259  0.022460  -3.69227           10.0
 30.0 -0.024436 -0.009381   7.03565            1.0
 35.0 -0.023708 -0.009344  -5.51584            1.0
 48.0 -0.010583 -0.008603  -3.05986           14.0
 49.0 -0.011022  0.000661  -3.79249           39.0
 75.0  0.011355 -0.002970   1.90074            1.0
 90.0  0.002929 -0.018336   2.66971            1.0
 97.0  0.011611  0.000815   0.00732           28.0
102.0  0.001209  0.020929  -4.71790            1.0
123.0 -0.001262 -0.002698   8.83149           22.0
128.0  0.006334 -0.001214  -5.15652           10.0
158.0 -0.022409 -0.003296  -7.26120            3.0
161.0  0.002261  0.001900 -10.54620            1.0
164.0 -0.008004  0.000373  -5.64786            8.0

Next steps:

• Use different trackleting and improve fitting to reduce residuals
• Attempt another approach suggested by Felix still using Ray Tracing techniques
  • Move from sampling of line inside voxels to voxel traversing