Medical Imaging Interaction Toolkit  2016.11.0
Medical Imaging Interaction Toolkit
The Random Parcellation View

Table of Contents

QmitkRandomParcellationIcon.png
Icon of the Random Parcellation View

Summary

This view is used to create a random parcellation of a segmented image. It was designed to parcellate grey matter to acquire random nodes for network analysis. This document describes how to use this view.

Overview

When parcellating an image you can either do a parcellation and use the generated parcellation. In which case it will have as many parcellation as seed points were specified. Instead you can choose to merge some of the original parcels untill either none below a threshold size remain, or until a specified minimum number of parcels has been reached. A further option is to "Just Merge Small Parcels". This will result in only neighbouring parcels below a threshold size being merged. Using this option may result in final parcels below the given threshold size if they are not adjacent to another parcel below the threshold size.

In all there are four possibilities to get a random parcellation

  • Without merging
  • Merging with respect to the number of nodes afterwards
  • Merging with respect to the number of voxels of the smallest parcel
  • Merging with respect to the number of voxels of the smallest parcel, where just small parcels are merged

Usage

First you have to get a segmentation of the target structure. You can use the Segmentation View to create one. Select it in the Data Manager. Now you can choose which of the possible ways of parcellation you would like. See Overview for these possibilities. Start the parcellation process using the "Select Random Nodes" button. You will get a new image in the Data Manager. Every voxel of a given parcel will have the same value (0 for background).

Methods

This Plug-In uses a kind of region growing algorithm. Random seed voxels are set on the segmentation and each is an own region at the beginning. Then voxels are added iteratively to each smallest region, first just the 6-connected ones and later the 26-connected ones and isolated voxels. Specific constraints have to be fulfilled before a voxel is added to a region. The algorithm stops when all voxels are contained in the regions. For merging it is checked which regions are neighbors of the smallest parcel and then a suitable one, according to a cost function, is merged to this parcel. When the merging algorithm stops depends on which stop-mechanism is chosen.