diff --git a/nibabies/interfaces/conftest.py b/nibabies/interfaces/conftest.py
index fee46075..ead38398 100644
--- a/nibabies/interfaces/conftest.py
+++ b/nibabies/interfaces/conftest.py
@@ -10,6 +10,8 @@
     'epi.nii',
     'T1w.nii',
     'func_to_struct.mat',
+    'atlas.nii',
+    'label_list.txt',
     'sub-01_run-01_echo-1_bold.nii.gz',
     'sub-01_run-01_echo-2_bold.nii.gz',
     'sub-01_run-01_echo-3_bold.nii.gz',
diff --git a/nibabies/interfaces/workbench.py b/nibabies/interfaces/workbench.py
index ba8a15b4..25334bbe 100644
--- a/nibabies/interfaces/workbench.py
+++ b/nibabies/interfaces/workbench.py
@@ -1,7 +1,152 @@
-from nipype.interfaces.base import CommandLineInputSpec, File, traits, TraitedSpec
+from nipype.interfaces.base import CommandLineInputSpec, File, traits, TraitedSpec, Str
 from nipype.interfaces.workbench.base import WBCommand
 
 
+class CiftiCreateDenseTimeseriesInputSpec(CommandLineInputSpec):
+    out_file = File(
+        name_source=["in_file"],
+        name_template="%s.dtseries.nii",
+        keep_extension=False,
+        argstr="%s",
+        position=0,
+        desc="The output CIFTI file",
+    )
+    in_file = File(
+        exists=True,
+        mandatory=True,
+        argstr="-volume %s",
+        position=1,
+        desc="volume file containing all voxel data for all volume structures",
+    )
+    structure_label_volume = File(
+        exists=True,
+        mandatory=True,
+        argstr="%s",
+        position=2,
+        desc="label volume file containing labels for cifti structures",
+    )
+    left_metric = File(
+        exists=True,
+        argstr="-left-metric %s",
+        position=3,
+        desc="metric file for left surface",
+    )
+    roi_left = File(
+        exists=True,
+        argstr="-roi-left %s",
+        position=4,
+        requires=['left_metric'],
+        desc="ROI (as metric file) of vertices to use from left surface",
+    )
+    right_metric = File(
+        exists=True,
+        argstr="-right-metric %s",
+        position=5,
+        desc="metric file for right surface",
+    )
+    roi_right = File(
+        exists=True,
+        argstr="-roi-right %s",
+        position=6,
+        requires=['right_metric'],
+        desc="ROI (as metric file) of vertices to use from right surface",
+    )
+    cerebellum_metric = File(
+        exists=True,
+        argstr="-cerebellum-metric %s",
+        position=7,
+        desc="metric file for cerebellum",
+    )
+    roi_cerebellum = File(
+        exists=True,
+        argstr="-roi-cerebellum %s",
+        position=8,
+        requires=['cerebellum_metric'],
+        desc="ROI (as metric file) of vertices to use from cerebellum",
+    )
+    timestep = traits.Float(
+        1.0,
+        usedefault=True,
+        argstr="-timestep %g",
+        desc="the timestep, in seconds",
+    )
+    timestart = traits.Float(
+        0.0,
+        usedefault=True,
+        argstr="-timestart %g",
+        desc="the time at the first frame, in seconds",
+    )
+    unit = traits.Enum(
+        "SECOND", "HERTZ", "METER", "RADIAN",
+        usedefault=True,
+        argstr="-unit %s",
+        desc="use a unit other than time")
+
+
+class CiftiCreateDenseTimeseriesOutputSpec(TraitedSpec):
+    out_file = File(exists=True, desc="CIFTI dense timeseries file")
+
+
+class CiftiCreateDenseTimeseries(WBCommand):
+    """
+    Create a CIFTI dense timeseries.
+
+    All input files must have the same number of columns/subvolumes.  Only
+    the specified components will be in the output cifti.  At least one
+    component must be specified.
+
+    See -volume-label-import and -volume-help for format details of label
+    volume files.  The structure-label-volume should have some of the label
+    names from this list, all other label names will be ignored:
+
+    CORTEX_LEFT
+    CORTEX_RIGHT
+    CEREBELLUM
+    ACCUMBENS_LEFT
+    ACCUMBENS_RIGHT
+    ALL_GREY_MATTER
+    ALL_WHITE_MATTER
+    AMYGDALA_LEFT
+    AMYGDALA_RIGHT
+    BRAIN_STEM
+    CAUDATE_LEFT
+    CAUDATE_RIGHT
+    CEREBELLAR_WHITE_MATTER_LEFT
+    CEREBELLAR_WHITE_MATTER_RIGHT
+    CEREBELLUM_LEFT
+    CEREBELLUM_RIGHT
+    CEREBRAL_WHITE_MATTER_LEFT
+    CEREBRAL_WHITE_MATTER_RIGHT
+    CORTEX
+    DIENCEPHALON_VENTRAL_LEFT
+    DIENCEPHALON_VENTRAL_RIGHT
+    HIPPOCAMPUS_LEFT
+    HIPPOCAMPUS_RIGHT
+    INVALID
+    OTHER
+    OTHER_GREY_MATTER
+    OTHER_WHITE_MATTER
+    PALLIDUM_LEFT
+    PALLIDUM_RIGHT
+    PUTAMEN_LEFT
+    PUTAMEN_RIGHT
+    THALAMUS_LEFT
+    THALAMUS_RIGHT
+
+    >>> from nibabies.interfaces.workbench import CiftiCreateDenseTimeseries
+    >>> createdts = CiftiCreateDenseTimeseries()
+    >>> createdts.inputs.in_file = data_dir /'functional.nii'
+    >>> createdts.inputs.structure_label_volume = data_dir /'atlas.nii'
+    >>> createdts.cmdline  #doctest: +ELLIPSIS +NORMALIZE_WHITESPACE
+    'wb_command -cifti-create-dense-timeseries functional.dtseries.nii \
+    -volume .../functional.nii .../atlas.nii -timestart 0 -timestep 1 -unit SECOND'
+    """
+
+    input_spec = CiftiCreateDenseTimeseriesInputSpec
+    output_spec = CiftiCreateDenseTimeseriesOutputSpec
+    _cmd = "wb_command -cifti-create-dense-timeseries"
+
+
 class CiftiDilateInputSpec(CommandLineInputSpec):
     in_file = File(
         exists=True,
@@ -132,3 +277,297 @@ class CiftiDilate(WBCommand):
     input_spec = CiftiDilateInputSpec
     output_spec = CiftiDilateOutputSpec
     _cmd = "wb_command -cifti-dilate"
+
+
+class VolumeAffineResampleInputSpec(CommandLineInputSpec):
+    in_file = File(
+        exists=True,
+        mandatory=True,
+        argstr="%s",
+        position=0,
+        desc="volume to resample",
+    )
+    volume_space = File(
+        exists=True,
+        mandatory=True,
+        argstr="%s",
+        position=1,
+        desc="a volume file in the volume space you want for the output",
+    )
+    method = traits.Enum(
+        "CUBIC", "ENCLOSING_VOXEL", "TRILINEAR",
+        mandatory=True,
+        argstr="%s",
+        position=2,
+        desc="The resampling method. The recommended methods are CUBIC "
+             "(cubic spline) for most data, and ENCLOSING_VOXEL for label data.",
+    )
+    out_file = File(
+        name_source=["in_file"],
+        name_template="resampled_%s.nii.gz",
+        keep_extension=True,
+        argstr="%s",
+        position=3,
+        desc="the output volume",
+    )
+    affine = File(
+        exists=True,
+        mandatory=True,
+        argstr="-affine %s",
+        position=4,
+        desc="the affine file to apply",
+    )
+    flirt = traits.Bool(
+        argstr="-flirt %s %s",
+        position=5,
+        desc="Set ``True`` if ``affine`` is a FLIRT affine.",
+    )
+    flirt_source_volume = File(
+        exists=True,
+        desc="the source volume used when generating the affine; defaults to in_file",
+        requires=['flirt'],
+    )
+    flirt_target_volume = File(
+        exists=True,
+        desc="the target volume used when generating the affine; defaults to volume_space",
+        requires=['flirt'],
+    )
+
+
+class VolumeAffineResampleOutputSpec(TraitedSpec):
+    out_file = File(exists=True, desc="the output volume")
+
+
+class VolumeAffineResample(WBCommand):
+    """
+    Resample a volume file with an affine transformation.
+
+    >>> from nibabies.interfaces.workbench import VolumeAffineResample
+    >>> resample = VolumeAffineResample()
+    >>> resample.inputs.in_file = data_dir /'functional.nii'
+    >>> resample.inputs.volume_space = data_dir /'anatomical.nii'
+    >>> resample.inputs.method = 'CUBIC'
+    >>> resample.inputs.affine = data_dir / 'func_to_struct.mat'
+    >>> resample.cmdline  #doctest: +ELLIPSIS +NORMALIZE_WHITESPACE
+    'wb_command -volume-resample .../functional.nii .../anatomical.nii CUBIC \
+    resampled_functional.nii.gz -affine .../func_to_struct.mat'
+
+    If the affine was generated with FLIRT, this should be indicated.
+    By default, the interface will use the ``in_file`` and ``volume_space``
+    for references.
+
+    >>> resample.inputs.flirt = True
+    >>> resample.cmdline  #doctest: +ELLIPSIS +NORMALIZE_WHITESPACE
+    'wb_command -volume-resample .../functional.nii .../anatomical.nii CUBIC \
+    resampled_functional.nii.gz -affine .../func_to_struct.mat \
+    -flirt .../functional.nii .../anatomical.nii'
+
+    However, if other volumes were used to calculate the affine, they can
+    be provided:
+
+    >>> resample.inputs.flirt_source_volume = data_dir / 'epi.nii'
+    >>> resample.inputs.flirt_target_volume = data_dir /'T1w.nii'
+    >>> resample.cmdline  #doctest: +ELLIPSIS +NORMALIZE_WHITESPACE
+    'wb_command -volume-resample .../functional.nii .../anatomical.nii CUBIC \
+    resampled_functional.nii.gz -affine .../func_to_struct.mat \
+    -flirt .../epi.nii .../T1w.nii'
+    """
+
+    input_spec = VolumeAffineResampleInputSpec
+    output_spec = VolumeAffineResampleOutputSpec
+    _cmd = "wb_command -volume-resample"
+
+    def _format_arg(self, opt, spec, val):
+        if opt == "flirt" and val:
+            val = (
+                self.inputs.flirt_source_volume or self.inputs.in_file,
+                self.inputs.flirt_target_volume or self.inputs.volume_space,
+            )
+        return super()._format_arg(opt, spec, val)
+
+
+class VolumeAllLabelsToROIsInputSpec(CommandLineInputSpec):
+    in_file = File(
+        exists=True,
+        mandatory=True,
+        argstr="%s",
+        position=0,
+        desc="the input volume label file",
+    )
+    label_map = traits.Either(
+        traits.Int, Str,
+        mandatory=True,
+        argstr="%s",
+        position=1,
+        desc="the number or name of the label map to use",
+    )
+    out_file = File(
+        name_source=["in_file"],
+        name_template="%s_rois.nii.gz",
+        argstr="%s",
+        position=2,
+        desc="the output volume",
+    )
+
+
+class VolumeAllLabelsToROIsOutputSpec(TraitedSpec):
+    out_file = File(exists=True, desc="the output volume")
+
+
+class VolumeAllLabelsToROIs(WBCommand):
+    """
+    Make ROIs from all labels in a volume frame
+
+    The output volume has a frame for each label in the specified input
+    frame, other than the ??? label, each of which contains an ROI of all
+    voxels that are set to the corresponding label.
+
+    >>> from nibabies.interfaces.workbench import VolumeAllLabelsToROIs
+    >>> rois = VolumeAllLabelsToROIs()
+    >>> rois.inputs.in_file = data_dir / 'atlas.nii'
+    >>> rois.inputs.label_map = 1
+    >>> rois.cmdline  #doctest: +ELLIPSIS
+    'wb_command -volume-all-labels-to-rois .../atlas.nii 1 atlas_rois.nii.gz'
+    """
+
+    input_spec = VolumeAllLabelsToROIsInputSpec
+    output_spec = VolumeAllLabelsToROIsOutputSpec
+    _cmd = "wb_command -volume-all-labels-to-rois"
+
+
+class VolumeLabelExportTableInputSpec(CommandLineInputSpec):
+    in_file = File(
+        exists=True,
+        mandatory=True,
+        argstr="%s",
+        position=0,
+        desc="the input volume label file",
+    )
+    label_map = traits.Either(
+        traits.Int, Str,
+        mandatory=True,
+        argstr="%s",
+        position=1,
+        desc="the number or name of the label map to use",
+    )
+    out_file = File(
+        name_source=["in_file"],
+        name_template="%s_labels.txt",
+        argstr="%s",
+        position=2,
+        desc="the output text file",
+    )
+
+
+class VolumeLabelExportTableOutputSpec(TraitedSpec):
+    out_file = File(exists=True, desc="the output text file")
+
+
+class VolumeLabelExportTable(WBCommand):
+    """
+    Export label table from volume as text
+
+    Takes the label table from the volume label map, and writes it to a text
+    format matching what is expected by -volume-label-import.
+
+    >>> from nibabies.interfaces.workbench import VolumeLabelExportTable
+    >>> label_export = VolumeLabelExportTable()
+    >>> label_export.inputs.in_file = data_dir / 'atlas.nii'
+    >>> label_export.inputs.label_map = 1
+    >>> label_export.cmdline  #doctest: +ELLIPSIS
+    'wb_command -volume-label-export-table .../atlas.nii 1 atlas_labels.txt'
+    """
+
+    input_spec = VolumeLabelExportTableInputSpec
+    output_spec = VolumeLabelExportTableOutputSpec
+    _cmd = "wb_command -volume-label-export-table"
+
+
+class VolumeLabelImportInputSpec(CommandLineInputSpec):
+    in_file = File(
+        exists=True,
+        mandatory=True,
+        argstr="%s",
+        position=0,
+        desc="the input volume file",
+    )
+    label_list_file = File(
+        exists=True,
+        mandatory=True,
+        argstr="%s",
+        position=1,
+        desc="text file containing the values and names for labels",
+    )
+    out_file = File(
+        name_source=["in_file"],
+        name_template="%s_labels.nii.gz",
+        argstr="%s",
+        position=2,
+        desc="the output workbench label volume",
+    )
+    discard_others = traits.Bool(
+        argstr="-discard-others",
+        desc="set any voxels with values not mentioned in the label list to the ??? label",
+    )
+    unlabeled_values = traits.Int(
+        0,
+        usedefault=True,
+        argstr="-unlabeled-value %d",
+        desc="the value that will be interpreted as unlabeled",
+    )
+    subvolume = traits.Either(
+        traits.Int, Str,
+        argstr="-subvolume %s",
+        desc="select a single subvolume to import (number or name)",
+    )
+    drop_unused_labels = traits.Bool(
+        argstr="-drop-unused-labels",
+        desc="remove any unused label values from the label table",
+    )
+
+
+class VolumeLabelImportOutputSpec(TraitedSpec):
+    out_file = File(exists=True, desc="the output workbench label volume")
+
+
+class VolumeLabelImport(WBCommand):
+    """
+    Import a label volume to workbench format
+
+    Creates a label volume from an integer-valued volume file.  The label
+    name and color information is stored in the volume header in a nifti
+    extension, with a similar format as in caret5, see -volume-help.  You may
+    specify the empty string (use "") for <label-list-file>, which will be
+    treated as if it is an empty file.  The label list file must have the
+    following format (2 lines per label):
+
+    <labelname>
+    <key> <red> <green> <blue> <alpha>
+    ...
+
+    Label names are specified on a separate line from their value and color,
+    in order to let label names contain spaces.  Whitespace is trimmed from
+    both ends of the label name, but is kept if it is in the middle of a
+    label.  Do not specify the "unlabeled" key in the file, it is assumed
+    that 0 means not labeled unless -unlabeled-value is specified.  The value
+    of <key> specifies what value in the imported file should be used as this
+    label.  The values of <red>, <green>, <blue> and <alpha> must be integers
+    from 0 to 255, and will specify the color the label is drawn as (alpha of
+    255 means fully opaque, which is probably what you want).
+
+    By default, it will create new label names with names like LABEL_5 for
+    any values encountered that are not mentioned in the list file, specify
+    -discard-others to instead set these values to the "unlabeled" key.
+
+    >>> from nibabies.interfaces.workbench import VolumeLabelImport
+    >>> label_import = VolumeLabelImport()
+    >>> label_import.inputs.in_file = data_dir / 'atlas.nii'
+    >>> label_import.inputs.label_list_file = data_dir / 'label_list.txt'
+    >>> label_import.cmdline  #doctest: +ELLIPSIS +NORMALIZE_WHITESPACE
+    'wb_command -volume-label-import .../atlas.nii .../label_list.txt \
+    atlas_labels.nii.gz -unlabeled-value 0'
+    """
+
+    input_spec = VolumeLabelImportInputSpec
+    output_spec = VolumeLabelImportOutputSpec
+    _cmd = "wb_command -volume-label-import"