Skip to content

Image registration with the QuadDIRECT optimization algorithm

License

Notifications You must be signed in to change notification settings

HolyLab/RegisterQD.jl

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

76 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

RegisterQD

CIcodecov

RegisterQD performs image registration using the global optimization routine QuadDIRECT. Unlike many other registration packages, this is not "greedy" descent based on an initial guess---it attempts to find the globally-optimal alignment of your images.

To use this package, users must choose between using either the CPU or the GPU. For CPU processing, you must manually load the RegisterMismatch package: using RegisterMismatch, RegisterQD. For GPU processing, you should instead load the RegisterMismatchCuda package: using RegisterMismatchCuda, RegisterQD. Note that loading both mismatch packages in the same session will cause method conflicts. Both mismatch packages are registered in the publicly-available HolyLabRegistry, and users are advised to add that registry. In the current absense of Github resources for GPU code, "gpu_test.jl" should be run on your personal machine as required.

This package exports the following registration functions:

  • qd_translate: register images by shifting one with respect to another (translations only)
  • qd_rigid: register images using rotations and translations
  • qd_affine: register images using arbitrary affine transformations

In general, using more degrees of freedom allows you to solve harder optimization problems, but also makes it harder to find the global optimum. Your best strategy is to permit no more degrees of freedom than needed to solve the problem.

See the help on these functions for details about how to call them.

Another important feature of this package is that it supports images that were sampled anisotropically. This is particularly common for three-dimensional biomedical imaging, where MRI and optical microscopy typically have one axis sampled at lower resolution. A rotation (from a rigid transformation) in physical space needs to be modified before applying it to an anisotropically-sampled image; see arrayscale and getSD for more information.

NOTE: see NEWS.md for information about a recent breaking change.