Hi Peter
If you don’t know the calibration you minimize
\sum |z-K.uncalibrate(pose.transform(point))|^2
If you already know the calibration, you can use
\sum |K.calibrate(z)-project(pose.transform(point))|^2
where K.calibrate needs to be done only once. Saves a ton of compute
Look for calibrate method of whatever calibration model you think of using.
FD