Main Content

intensityFeatures

Radiomics intensity features

Since R2023b

collapse all in page

    Syntax

    I = intensityFeatures(R)
    I = intensityFeatures(R,Name=Value)

    Description

    example

    I = intensityFeatures(R) computes the radiomics intensity features I for the radiomics object R.

    I = intensityFeatures(R,Name=Value) specifies additional options using one or more optional name-value arguments.

    Examples

    collapse all

    Open Live Script

    Import a computed tomography (CT) image volume and the corresponding ROI mask volume from the IBSI validation data set [1][2][3] as medicalVolume objects.

    unzip("CTImageMaskNIfTI.zip")
data = medicalVolume("CT_image.nii.gz");
roi = medicalVolume("CT_mask.nii.gz");

    Visualize a slice of the CT image volume and the corresponding ROI.

    figure
imshowpair(data.Voxels(:,:,20),roi.Voxels(:,:,20),"montage")

    Create a radiomics object, using the CT image volume and ROI mask volume, with default preprocessing options.

    R = radiomics(data,roi)
    R = 
  radiomics with properties:

                             Data: [1x1 medicalVolume]
                         ROILabel: [1x1 medicalVolume]
                         Resample: 1
                        Resegment: 1
                       Discretize: 1
                    DiscretizeIVH: 1
            ResampledVoxelSpacing: 1
               DataResampleMethod: 'linear'
               MaskResampleMethod: 'linear'
              ResegmentationRange: []
                  ExcludeOutliers: 1
       DiscreteBinSizeOrBinNumber: []
                   DiscreteMethod: 'FixedBinNumber'
    DiscreteIVHBinSizeOrBinNumber: []
                DiscreteIVHMethod: 'FixedBinNumber'

    Compute the local intensity features of the ROI in the 2D-resampled CT image volume.

    I = intensityFeatures(R,Type="LocalIntensity",SubType="2D")
    I=1×3 table
    LabelID    LocalIntensityPeak2D    GlobalIntensityPeak2D
    _______    ____________________    _____________________

      "1"             71.689                  78.084

    [1] Vallières, Martin, Carolyn R. Freeman, Sonia R. Skamene, and Issam El Naqa. “A Radiomics Model from Joint FDG-PET and MRI Texture Features for the Prediction of Lung Metastases in Soft-Tissue Sarcomas of the Extremities.” The Cancer Imaging Archive, 2015. https://doi.org/10.7937/K9/TCIA.2015.7GO2GSKS.

    [2] Vallières, M, C R Freeman, S R Skamene, and I El Naqa. “A Radiomics Model from Joint FDG-PET and MRI Texture Features for the Prediction of Lung Metastases in Soft-Tissue Sarcomas of the Extremities.” Physics in Medicine and Biology 60, no. 14 (July 7, 2015): 5471–96. https://doi.org/10.1088/0031-9155/60/14/5471.

    [3] Clark, Kenneth, Bruce Vendt, Kirk Smith, John Freymann, Justin Kirby, Paul Koppel, Stephen Moore, et al. “The Cancer Imaging Archive (TCIA): Maintaining and Operating a Public Information Repository.” Journal of Digital Imaging 26, no. 6 (December 2013): 1045–57. https://doi.org/10.1007/s10278-013-9622-7.

    Input Arguments

    collapse all

    Data and ROI for feature computation, specified as a radiomics object. The radiomics object R contains details of the preprocessed data and region of interest (ROI) from which to compute the features.

    Name-Value Arguments

    Specify optional pairs of arguments as Name1=Value1,...,NameN=ValueN, where Name is the argument name and Value is the corresponding value. Name-value arguments must appear after other arguments, but the order of the pairs does not matter.

    Example: intensityFeatures(R,Type="LocalIntensity",SubType="2D") computes the local intensity features on 2D resampled data.

    Category of intensity features to compute, specified as one of these options.

    • "LocalIntensity"

    • "IntensityBasedStatistics"

    • "IntensityHistogram"

    • "IntensityVolumeHistogram"

    • "all"

    If you specify "all", the function computes every category of intensity features. For more information on which specific intensity features each category includes, see IBSI Standard and Radiomics Function Feature Correspondences.

    Data Types: char | string

    Resampling from which to compute intensity features, specified as one of these options.

    • "2D" — Computes features from the 2D resampled volume.

    • "3D" — Computes features from the 3D resampled volume.

    • "all" — Computes features for both options.

    When you 2D-resample the volume, the function makes the voxel spacing along the x- and y-dimensions isotropic, but the voxel spacing along the z-dimension is the same as in the input volume. When you 3D-resample the volume, the function makes the voxel spacing along all three spatial dimensions isotropic.

    Data Types: char | string

    Output Arguments

    collapse all

    Intensity features, returned as a table. The first column in I is LabelID. The subsequent columns are the intensity features. For more details on which intensity features are computed in each Type and SubType, see IBSI Standard and Radiomics Function Feature Correspondences.

    Version History

    Introduced in R2023b

    See Also

    | |

    Topics

    External Websites