Yes Alexey,
i would like to save in a TPs file, the X, Y, Z, and Error X, Error Y, Error Z for each point of course related to crs projected system as in
Coordenadas de puntos de enlace con sus errores RMS en X,Y,Z, error de reproyeccion promedio y numero de proyecciones
Sistema de coordenadas WGS 84 / UTM zone 16N + GGM2010 geoid height
Index X Y Z Error_X(cm) Error_Y(cm) Error_Z(cm) Error_total(cm) num_obs errrep(pix)
0 236376.294 2325813.957 8.117 2.649 3.897 5.733 7.421 2 0.508
Right now I am using this piece of code to transform error components from internal CS to projected CRS:
M = chunk.transform.matrix
s = chunk.transform.scale
points = chunk.point_cloud.points
npoints = len(points)
....
for point_id in range(npoints):
if not points[point_id].valid:
continue
if chunk.crs:
V = M * (points[point_id].coord)
V.size = 3
#print (V)
X = chunk.crs.project(V)
m = chunk.crs.localframe(V)
R = m.rotation() * M.rotation()
rmsXYZ = s * Metashape.Vector([math.sqrt(points[point_id].cov[0,0]), math.sqrt(points[point_id].cov[1,1]), math.sqrt(points[point_id].cov[2,2])])
rmsXYZ = R * rmsXYZ
rmsXYZ = Metashape.Vector([abs(rmsXYZ.x),abs(rmsXYZ.y),abs(rmsXYZ.z)])
rmsXYZ = 100 * rmsXYZ
I do appreciate your input,