Skip to content

How does ANTs handle qform and sform in NIFTI 1 images?

Jump to bottom
Philip Cook edited this page Oct 5, 2020 · 29 revisions

ANTs NIFTI I/O will change in the next release

Because ITK has changed its NIFTI-1 I/O conventions, ANTs binaries will also change. The source code as already been updated to reflect this, so the next release will contain the new method. This page will be updated accordingly. Briefly, the preference for qform has been removed, and sform is now used wherever possible.

ITK NIFTI I/O for ANTs 2.3.4 and earlier

Short answer

ANTs uses ITK image I/O, which uses a common framework to read many different image formats. ITK supports a rotation, but not a full affine transform, to translate coordinates from index space (voxel indices) to physical space. Therefore it can only use qform when reading a NIFTI-1 image. When writing NIFTI output, only the qform is written; the sform is replaced with zeros.

Long answer

The exact procedure for defining the transform to physical space is here: itkNiftiImageIO.cxx

In general, qform is used on input. The rotation comes from the quaternion, the scaling comes from the pixel dimensions, and the translation from the qoffset fields. This is NIFTI "method 2" as explained in the NIFTI-1 documentation. This is the procedure whenever the qform code is greater than 0, regardless of the sform code.

In rare cases, the qform code is 0. This should never be the case except for backwards compatibility, ie reading old ANALYZE files. If the qform code is 0, indicating no valid qform, the rotation and scaling are extracted from the sform, if it is present.

When writing NIFTI files, ITK encodes the rotation, translation and scaling in the qform. The qform code is set to 1 (NIFTI_XFORM_SCANNER_ANAT), and the sform code is set to 0. The loss of the sform is unavoidable because the ITK I/O code reads images into a common internal structure, independent of the individual image format on disk, and only supports a rigid transformation.

Precision errors in header transforms

ITK reads and writes a variety of file formats. Images are read into a generic ITK image object, and after processing the results are written out in the desired format. As a result, there can be precision errors in the output NIFTI header.

Potential interoperability problems with other software

  1. qform and sform codes may be interpreted differently. ITK uses qform whenever qform_code > 0.

  2. The sform will be overwritten in output images, and sform code will be set to 0. If the sform is not identical to the qform, information will be lost.

  3. Another issue: "Note that in cases where both an SForm and QForm are present, ITK gives precedence to the QForm while SPM/MRIcroGL/FSLeyes/Mango give precedence to the SForm. The rationale for the SPM approach is clear: QForms were designed to provide an idea of native space, while the SForm could provide the user an easy way to provide a starting estimate. When informed of this difference the ITK group elected to keep their method to ensure backward compatibility. To support both conventions, recent versions of SPM change both the SForm as well as QForm when the user changes the origin with the "Display" button. While the original intent of the SForm and QForm was to provide a method to encode both native scanner as well as some losslessly re-oriented coordinates, the differences in implementations supports SPM's method of using the same value for both to ensure consistent display and processing."

  4. In registration, the fixed image may have an sform, which will be ignored in the registration. The warped moving image will have only a qform. Thus the registration may work correctly but appear to be incorrect because of differences in how the NIFTI header is interpreted.

Clone this wiki locally