Forum

Author Topic: How is the total marker error (pix) being calculated  (Read 5657 times)

colin.spexigeo

  • Newbie
  • *
  • Posts: 13
    • View Profile
How is the total marker error (pix) being calculated
« on: February 07, 2020, 01:11:31 AM »
Hi there,

I'm interested in how the "Total" value of marker errors is being calculated. I had assumed that it was an average (mean) of all camera projection errors, but this does not seem to be the case.

For example, Agisoft reports the following:

Errors
-dYUxAmaeU0.143
-YDOHdBJSz0.458
-YOo6Otjv50.629
CYzKGPFuwI0.437
f6n4TJCcSD0.054
FqIO_ToeyF0.397
Fy4qsbY3ft0.661
ladYl1TVsF0.497
nziQ02gyo50.322
pfZjrVEIQA0.518
rFtTrwWuiN1.332
tQGuj6cN-y0.585
VcznBnOJaW0.358
XDxBg3CLSh0.135
Total0.554 pix

However, the *average* values of all these errors is actually 0.466.

Does anybody know how the "total" values is being calculated?

Thanks in advance,
Colin

Paulo

  • Hero Member
  • *****
  • Posts: 1557
    • View Profile
Re: How is the total marker error (pix) being calculated
« Reply #1 on: February 07, 2020, 03:43:46 AM »
The total error is the RMS of  marker errors

Marker           Error
dYUxAmaeU   0.143
-YDOHdBJSz   0.458
-YOo6Otjv5   0.629
CYzKGPFuwI   0.437
f6n4TJCcSD   0.054
FqIO_ToeyF   0.397
Fy4qsbY3ft   0.661
ladYl1TVsF   0.497
nziQ02gyo5   0.322
pfZjrVEIQA   0.518
rFtTrwWuiN   1.332
tQGuj6cN-y   0.585
VcznBnOJaW   0.358
XDxBg3CLSh   0.135
RMS           0.554  SQRT(sum.squares(B2:B15)/14)

« Last Edit: February 07, 2020, 03:53:45 AM by Paulo »
Best Regards,
Paul Pelletier,
Surveyor

colin.spexigeo

  • Newbie
  • *
  • Posts: 13
    • View Profile
Re: How is the total marker error (pix) being calculated
« Reply #2 on: February 07, 2020, 04:38:52 AM »
Thank you Paulo!

P.S. Off-hand, you wouldn't happen to know an easy way to calculate the "Error (m)" value, would you? I imagine you'd have to project the marker points/error offsets through the cameras (maybe projected onto a plane that's facing the sensor) and run an RMS on the cartesian-space error values. I'm trying to calculate this information and surface it to the user outside of the program.

Paulo

  • Hero Member
  • *****
  • Posts: 1557
    • View Profile
Re: How is the total marker error (pix) being calculated
« Reply #3 on: February 07, 2020, 07:06:23 AM »
Assuming the markers have reference locations, then individual error for each marker is calculated as:

Code: [Select]
for marker in chunk.markers:
      source = marker.reference.location
      estim = chunk.crs.project(chunk.transform.matrix.mulp(marker.position))
      error = estim - source
      total = error.norm()
      print(marker.label, error.x, error.y, error.z, total)

and total error again is calculated as RMS of each marker error...Example:

Code: [Select]
Marker Error.X(m)  Error.Y(m) Error.Z(m)   Error(m) Error(pix)
PA03 0.077        0.047 0.017        0.092 0.120
PA08 0.041       -0.068 -0.020        0.082 0.302
PA12 -0.052        0.031 -0.005        0.061 0.254
PA18 -0.053        0.019 -0.005        0.056 0.086
PA01 -0.016       -0.041 -0.004        0.044 0.228
PA19 0.006       -0.022 0.017        0.029 0.094
PA16 -0.017       -0.006 -0.005        0.019 0.135
PA11 0.007        0.014 0.010        0.019 0.117
              RMS      0.056    0.184

« Last Edit: February 07, 2020, 07:11:50 AM by Paulo »
Best Regards,
Paul Pelletier,
Surveyor