This site uses cookies to assist with navigation and your ability to provide feedback. Cookie Policy

Radiomic features

First Order Features - Histogram

Histogram calculation

To build a histogram, it is necessary to determine a bin width ("bin" parameter). The indices derived from the histogram will depend on this bin width parameter.

This dependence, similar to that found in texture index calculations, is often overlooked in publications.

In LIFEx, with the absolute model the histogram is built a number of bins equal to that entered by the user in the "number of grey level" and "size of bin" fields of the resampling menu.

In LIFEx, with the relative model the histogram is built only with "number of grey level" fields of the resampling menu that entered by the user and min and max are extracted values of each ROI.

Read more: First Order Features - Histogram

First Order Features - Shape

SHAPE_Sphericity is how spherical a Volume of Interest is. Sphericity is equal to 1 for a perfect sphere.


where \(V\) and \(A\) correspond to the volume and the surface of the Volume Of Interest based on the Delaunay triangulation.

Read more: First Order Features - Shape

Grey-Level Zone Length Matrix (GLZLM)

The grey-level zone length matrix (GLZLM) provides information on the size of homogeneous zones for each grey-level in 3 dimensions. It is also named Grey Level Size Zone Matrix (GLSZM).

From this matrix, 11 texture indices can be computed. Element \((i,j)\) of GLZLM corresponds to the number of homogeneous zones of \(j\) voxels with the intensity \(i\) in an image and is called \(GLZLM(i,j)\) thereafter.

Read more: Grey-Level Zone Length Matrix (GLZLM)

Grey-Level Run Length Matrix (GLRLM)

The grey-level run length matrix (GLRLM) gives the size of homogeneous runs for each grey level. This matrix is computed for the 13 different directions in 3D (4 in 2D) and for each of the 11 texture indices derived from this matrix, the 3D value is the average over the 13 directions in 3D (4 in 2D). The element \((i,j)\) of GLRLM corresponds to the number of homogeneous runs of \(j\) voxels with intensity \(i\) in an image and is called \(GLRLM(i,j)\) thereafter.

Read more: Grey-Level Run Length Matrix (GLRLM)

Neighborhood Grey-Level Different Matrix (NGLDM)

The neighborhood grey-level different matrix (NGLDM) corresponds to the difference of grey-level between one voxel and its 26 neighbours in 3 dimensions (8 in 2D). Three texture indices can be computed from this matrix. An element \((i,1)\) of NGLDM corresponds to the probability of occurrence of level \(i\) and an element \((i,2)\) is equal to:
NGLDM(i,2)= \sum_{p}\sum_{q} \left\lbrace
|\overline{M}(p,q)-i| \mbox{~if $I(p,q)=i$} \\
0 \mbox{~else}

where \(\overline{M}(p,q)\) is the average of intensities over the 26 neighbour voxels of voxel \((p,q)\).

Read more: Neighborhood Grey-Level Different Matrix (NGLDM)

Grey Level Co-occurrence Matrix (GLCM)

The grey level co-occurrence matrix (GLCM) [Haralick] takes into account the arrangements of pairs of voxels to calculate textural indices. The GLCM is calculated from 13 different directions in 3D with a \(\delta\)-voxel distance (\(\|\overrightarrow{d\|}\)) relationship between neighboured voxels. The index value is the average of the index over the 13 directions in space (X, Y, Z). Six textural indices can be computed from this matrix. An entry \((i,j)\) of GLCM for one direction is equal to:

GLCM_{\Delta x,\Delta y}(i,j)= \frac{1}{Pairs_{ROI}}\sum_{p=1}^{N-\Delta x}\sum_{q=1}^{M-\Delta y}
\mbox{1 if ($I(p,q)=i$, $I(p+\Delta x, q+\Delta y)=j$) } \\
\mbox{ and $I(p,q), I(p+\Delta x, q+\Delta y) \in ROI$ } \\
\mbox{0 otherwise}

where \(I(p,q)\) corresponds to voxel \((p,q)\) in an image (\(I\)) of size \(N * M\). The vector \(\overrightarrow{d}=(\Delta x,\Delta y)\) covers the 4 directions (D1, D2, D3, D4, in 2D space or 13 directions (D1, D2, ..., D13, in 3D space and \(Pairs_{ROI}\)

Read more: Grey Level Co-occurrence Matrix (GLCM)

Login Form