Camera model representation in memory and on disk
Python interface
Interfacing with camera models is done in Python with the mrcal.cameramodel
class. This describes one camera; a calibrated set of cameras is represented
by multiple objects. Each object always contains
- The
intrinsics
: lens parameters describing one of the lens models - The
extrinsics
: a pose of the camera in space. This pose is represented as a transformation between some common reference coordinate system and the camera coordinate system. The specific meaning of the reference coordinate system is arbitrary, but all the cameras in a calibrated set must be defined in respect to the one common reference.
Each camera model object may also contain:
- The
optimization_inputs
: all the data used to compute the model initially. Used for the uncertainty computations and any after-the-fact analysis. - The
valid_intrinsics_region
: a contour in the imager where the projection behavior is "reliable". This is usually derived from the uncertainty plot, and used as a shorthand. It isn't as informative as the uncertainty plot, but such a valid-intrinsics contour is often convenient to have and to visualize.
C interface
The C API uses the mrcal_cameramodel_t
structure to represent a model. This
contains just the bare minimum:
- intrinsics (
mrcal_lensmodel_t lensmodel
,double intrinsics[0]
) - extrinsics (
double rt_cam_ref[6]
) - imager size (
unsigned int imagersize[2]
)
Note that the intrinsics data has size 0 because the size of this array depends on the specific lens model being used, and is unknown at compile time.
So it is an error to define this on the stack. Do not do this:
void f(void) { mrcal_cameramodel_t model; // ERROR }
If you need to define a known-at-compile-time model on the stack you can use the lensmodel-specific cameramodel types:
void f(void) { mrcal_cameramodel_LENSMODEL_OPENCV8_t model; // OK }
This only exists for models that have a constant number of parameters; notably
there is no mrcal_cameramodel_LENSMODEL_SPLINED_STEREOGRAPHIC_t
. When reading
a model from disk, mrcal dynamically allocates the right amount of memory, and
returns a mrcal_cameramodel_t*
.
The C API has a simple interface for reading/writing .cameramodel
data:
mrcal_cameramodel_t* mrcal_read_cameramodel_string(const char* string, int len); mrcal_cameramodel_t* mrcal_read_cameramodel_file (const char* filename); void mrcal_free_cameramodel(mrcal_cameramodel_t** cameramodel); bool mrcal_write_cameramodel_file(const char* filename, const mrcal_cameramodel_t* cameramodel);
File formats
Two different file formats are available:
.cameramodel
: the mrcal-native format. This is a plain text representation of a Pythondict
describing all the fields. This is the preferred format, and the only format supported by the C API.cahvor
: the legacy format available for compatibility with existing JPL tools. If you don't need to interoperate with tools that require this format, there's little reason to use it. This format cannot store splined models oroptimization_inputs
, so it cannot be used for the uncertainty computations. Not supported in C
Note that both file formats can describe lenses using the LENSMODEL_CAHVOR
model.
The mrcal.cameramodel
class will intelligently pick the correct file format
based on the data (if reading) and the filename (if writing). The
mrcal-to-cahvor
and mrcal-from-cahvor
tools can be used to convert between
the two file formats.
Sample usages
See the API documentation for usage details.
Grafting two models
A trivial example to
- read two models from disk
- recombine into a joint model that uses the lens parameters from one model with geometry from the other
- write to disk
import mrcal model_for_intrinsics = mrcal.cameramodel('model0.cameramodel') model_for_extrinsics = mrcal.cameramodel('model1.cameramodel') model_joint = mrcal.cameramodel( model_for_intrinsics ) extrinsics = model_for_extrinsics.extrinsics_rt_fromref() model_joint.extrinsics_rt_fromref(extrinsics) model_joint.write('model-joint.cameramodel')
This is the basic operation of the mrcal-graft-models
tool.
Re-optimizing a model
To re-optimize a model from its optimization_inputs
:
import mrcal m = mrcal.cameramodel('camera.cameramodel') optimization_inputs = m.optimization_inputs() mrcal.optimize(**optimization_inputs) model_reoptimized = \ mrcal.cameramodel( optimization_inputs = m.optimization_inputs(), icam_intrinsics = m.icam_intrinsics() )
Here we asked mrcal to re-optimize the data used to compute the given model originally. We didn't make any changes to the inputs, and we should already have an optimal solution, so this re-optimized model would be the same as the initial one. But we could tweak optimization problem before reoptimizing, and this would give us an nice way to observe the effects of those changes. We can add input noise or change the lens model or regularization terms or anything else.