Propagating camera calibration information after performing a geometric operation

• See also
Reference:
McalAlloc McalAssociate McalWarp MgenWarpParameter MimFlip MimResize MimRotate MimTranslate MimWarp
Availability
Not available in MIL-Lite


Available in MIL
Previous
Next

When you have a calibrated image and perform a geometric operation, the resulting image does not retain its camera calibration information. However, you can create an appropriately warped version of the camera calibration context and associate it with the uncalibrated geometrically transformed image.

McalWarp() allows you to apply a warping to the pixel-to-world mapping of a camera calibration context and use it to configure a specified camera calibration context. The newly configured context can then be associated with a similarly warped image. McalWarp() is meant to be used in conjunction with MimWarp(). As such, it takes the same warping parameter values used to transform the image. Therefore, McalWarp() can perform first-order polynomial warpings, perspective polynomial warpings, and custom warpings using LUTs.

If you want to apply a flip, resize, rotate, or translate operation on an image and wanted the destination image to be appropriately calibrated, you must perform the transformation of the image using MimWarp() instead of MimFlip(), MimResize(), MimRotate(), and MimTranslate(); you can then use McalWarp() followed by McalAssociate() to appropriately calibrate the transformed image. Note that, MimWarp() and McalWarp() use a matrix of coefficients (or LUTs) to describe its transformations. To generate the coefficients that represent the required transform, you can use MgenWarpParameter() to first generate the required transformation matrix (or LUTs) and subsequently call MimWarp() and McalWarp().

To perform a geometric operation on an image and preserve the camera calibration information, perform the following:

When warping the pixel-to-world mapping of a camera calibration context, the destination camera calibration context need not be in the same camera calibration mode as the source context. McalWarp() approximates the equivalent calibration when converting from one mode to the other.

The following code snippet shows how to perform a rotation using MimWarp(). It generates the required transformation matrices for the rotation using MgenWarpParameter(). It then shows how to warp the pixel-to-world mapping of the existing camera calibration context and propagate the new camera calibration information to the warped image using McalWarp() and then McalAssociate() respectively.

// Generate the warping matrix to rotate the source image +90 degrees around its center.
// This is the composition of a translation to the source's center, a rotation, and then a
// translation from the destination's center.
MIL_DOUBLE SrcCenterX = 0.5*(SrcSizeX-1);
MIL_DOUBLE SrcCenterY = 0.5*(SrcSizeY-1);
MIL_DOUBLE DstCenterX = 0.5*(DstSizeX-1);
MIL_DOUBLE DstCenterY = 0.5*(DstSizeY-1);
MIL_ID WarpMatrixArrayBufId = MbufAlloc2d(SysId, 3, 3, 32+M_FLOAT, M_ARRAY, M_NULL);
MgenWarpParameter(M_NULL, WarpMatrixArrayBufId, M_NULL, M_WARP_POLYNOMIAL,
   M_TRANSLATE, -SrcCenterX, -SrcCenterY);
MgenWarpParameter(WarpMatrixArrayBufId, WarpMatrixArrayBufId, M_NULL, M_WARP_POLYNOMIAL,
   M_ROTATE, 90.0, M_NULL);
MgenWarpParameter(WarpMatrixArrayBufId, WarpMatrixArrayBufId, M_NULL, M_WARP_POLYNOMIAL,
   M_TRANSLATE, DstCenterX, DstCenterY);

// Warp the image.
MimWarp(SrcImageBufId, DstImageBufId, WarpMatrixArrayBufId, M_NULL, M_WARP_POLYNOMIAL,
   M_BILINEAR+M_OVERSCAN_CLEAR);
// Same as:
// MimRotate(SrcImageBufId, DstImageBufId, 90.0, M_DEFAULT, M_DEFAULT, M_DEFAULT,
//    M_DEFAULT, M_BILINEAR+M_OVERSCAN_CLEAR);

// Warp the pixel-to-world mapping of the calibration context 
// and associate the context to the warped image.
MIL_ID DstContextCalId = McalAlloc(SysId, M_LINEAR_INTERPOLATION, M_DEFAULT, M_NULL);
McalWarp(SrcImageBufId, DstContextCalId, WarpMatrixArrayBufId, M_NULL, 0.0, 0.0,
   (MIL_DOUBLE)DstSizeX, (MIL_DOUBLE)DstSizeY, M_DEFAULT, M_DEFAULT, M_WARP_POLYNOMIAL,
   M_DEFAULT);

if (McalInquire(DstContextCalId, M_CALIBRATION_STATUS, M_NULL) == M_CALIBRATED)
   {
   McalAssociate(DstContextCalId, DstImageBufId, M_DEFAULT);
   }
else
   {
   MosPrintf(MIL_TEXT("Calibration warping has failed.\n"));
   }