/Inter-Modality-Image-Registration

Deep Learning based Inter-Modality Image Registration Supervised by Intra-Modality Similarity

Primary LanguageJupyter Notebook

Deep Learning based Inter-Modality Image Registration Supervised by Intra-Modality Similarity

The aim of the project is to reveal the geometric accuracy of tomography images and the soft tissue resolution quality of MRI images with deep learning-based image registration method.

Principles of Radiography

Radiology is the use of x-rays and other imaging methods in medicine for diagnosis and treatment.

COMPUTED TOMOGRAPHY ( CT )

By taking many two-dimensional X-ray images of an object from different angles, it is tried to obtain a three-dimensional image of the internal structure of that object.

  • Creates a detailed and layered picture by scanning the relevant area or areas of the body with X-r
  • Gives the closest results geometrically.
  • Helpful in applications such as drug dosage for disease treatment.

image

MAGNETIC RESONANCE ( MRI )

Magnetic Resonance Imaging (MRI) is a technique that uses magnetic fields and radio waves to create detailed images of organs and tissues in the body

  • High-resolution images of the inside of the body are obtained to diagnose various problems.
  • Most widely used method for disease diagnosis.

image

DIFFERENCE BETWEEN MRI AND CT

CT
image

  • High geometric accuracy.
  • Low image resolution.
    MR
    image

PROJECT DEVELOPMENT ENVIRONMENT

Examining the development environment, Anaconda3, Spyder IDE used for the project because Spyder provides more convenience for the data processing, such as data visualization, ITK, Variable Tracking. Keras, Tensorflow libraries and deep learning modules will be used for the image registration. PyQt library used for the desktop application.

BASIC IMAGE REGISTRATION

There are four basic steps incluiding Interpolation, Metric, Optimization, Transformation in the project.

image

Interpolation

  • When a point is mapped from one space to another with a transformation, it is generally mapped to a non-grid location. Therefore, interpolation is necessary to evaluate the image density at the mapped location.
  • Transformation maps point from the still image coordinate system to the moving image coordinate system.

image

Interpolation Types

Linear Interpolation : A method that assumes the density varies linearly between grid positions.

B-Spline Interpolation : A method for comparing image density using the B-spline function.

Windowed Sinc Interpolation : A method that performs interpolation based on Fourier analysis.

Windowed Sinc Interpolation used in the project

Metric

  • This is the most critical part of the registration process.

  • The choice of the metric algorithm is directly determined by the problem to be solved.

  • We used the histogram-based mutual information method because it is highly recommended in medical registrations.

image

Mutual Information Metric Graph

image

Optimization

  • The basic input to an optimizer is a cost function or metric object.

  • The metric is set using optimization algorithm.

  • Once the optimization has finished, the final parameters can be obtained.

  • Gradient Descent method was used in the Project.

Transformation

  • This stage is where the conversion is done between the two images.

  • We performed your process using the Euler 3D method, which is more recommended for geometric transformation in the literature.

image

Registration Result

We have a poor result in terms of geometrically smooth image quality in our tomography image. In MR, we have a very good geometric output in terms of image quality and a poor output in terms of size. So how are things in our fusion image? Let's focus especially on the bone area. Geometrically it is the same as in the tomography view, and in terms of image quality it is the same as in the MRI image. So an image that meets our expectations.

image

Alpha Blending

It combines the pixels in the two images and allows the images we have joined to be seen overlapping.

PURPLE : Purple Areas are images obtained as a result of tomography.
GREEN : Green Areas are images obtained as a result of MRI.

image

Landmark Validation

image

image

DEEP LEARNING METHODS!

Registration Without Any Feature Extraction ( Classic Deep Learning )

  • 15 patient data with jpeg format was read.
  • Output layer activation function is sigmoid applied.
  • Occuracy %30

image

Canny Edge Feature Extraction

  • Sk-Image library used for feature extraction.
  • Suitable for framing the Picture.

image

Sobel Edge Detection Feature Extraction

  • Sk-Image library used for feature extraction.
  • More advanced version of the high-pass filter.

image

AUTOENCODERS

  • The purpose of the autoencoder is to learn a symbolic vector that represents the data.

image

ARRANGEMENT THE NEURAL NETWORK’S HIDDEN LAYERS

  • Output activation function changed as ReLU.
  • More hidden layer added.
  • Accuracy %94

image

HIGH PASS & LOW PASS FEATURE EXTRACTION METHOD!

High-Pass : It was not used because of Sobel Edge Detection. Low-Pass : It was used to help improve accuracy in MR and CT images. The result of the study was yield.

Superpixel Segmentation Benchmark Work

image

image

image

image

Result

  • Mean absolute error decrease to 2.72.
  • Accuracy increased by 3.7 after superpixel and low-pass filter feature extractions.
  • Accuracy %97.7

image