mosamatic2 2.0.4__py3-none-any.whl → 2.0.5__py3-none-any.whl

mosamatic2/cli.py

@@ -5,6 +5,7 @@ from mosamatic2.commands import (
     segmentmusclefatl3tensorflow,
     createpngsfromsegmentations,
     dicom2nifti,
+    selectslicefromscans,
 )
 from mosamatic2.core.utils import show_doc_command
 
@@ -31,4 +32,5 @@ main.add_command(rescaledicomimages.rescaledicomimages)
 main.add_command(segmentmusclefatl3tensorflow.segmentmusclefatl3tensorflow)
 main.add_command(createpngsfromsegmentations.createpngsfromsegmentations)
 main.add_command(dicom2nifti.dicom2nifti)
+main.add_command(selectslicefromscans.selectslicefromscans)
 main.add_command(show_doc_command(main)) # Special command to show long description for command

mosamatic2/commands/selectslicefromscans.py

@@ -0,0 +1,65 @@
+import click
+
+from mosamatic2.core.tasks import SelectSliceFromScansTask
+
+
+@click.command(help='Selects specific slice from CT scans')
+@click.option(
+    '--scans', 
+    required=True, 
+    type=click.Path(exists=True), 
+    help='Directory with scans (each patient separate subdirectory)',
+)
+@click.option(
+    '--output', 
+    required=True, 
+    type=click.Path(), 
+    help='Output directory'
+)
+@click.option(
+    '--vertebra', 
+    required=True,
+    help='Vertebral level for selecting slice (default: "L3")'
+)
+@click.option(
+    '--overwrite', 
+    type=click.BOOL, 
+    default=False, 
+    help='Overwrite [true|false]'
+)
+def selectslicefromscans(scans, vertebra, output, overwrite):
+    """
+    Selects specific slice from list of CT scans
+    
+    Parameters
+    ----------
+    --scans : str
+        Directory to scans. Each patient's scan should be in a separate
+        subdirectory. For example:
+
+        /scans
+        /scans/patient1
+        /scans/patient1/file1.dcm
+        /scans/patient1/file2.dcm
+        ...
+        /scans/patient2
+        ...
+
+    --output : str
+        Path to output directory where selected slices will be placed. Each
+        slice's file name will be the same as the scan directory name, so in
+        the example above that would be "patient1", "patient2", etc.
+
+    --vertebra : str
+        Vertebral level where to take slice [L3|T4]
+
+    --overwrite : bool
+        Overwrite contents output directory [true|false]
+    """
+    task = SelectSliceFromScansTask(
+        inputs={'scans': scans},
+        params={'vertebra': vertebra},
+        output=output,
+        overwrite=overwrite,
+    )
+    task.run()

mosamatic2/core/data/multinumpyimage.py

@@ -0,0 +1,26 @@
+import os
+from mosamatic2.core.managers.logmanager import LogManager
+from mosamatic2.core.data.filedata import FileData
+from mosamatic2.core.data.numpyimage import NumPyImage
+
+LOG = LogManager()
+
+
+class MultiNumPyImage(FileData):
+    def __init__(self):
+        super(MultiNumPyImage, self).__init__()
+        self._images = []
+
+    def images(self):
+        return self._images
+    
+    def load(self):
+        if self.path():
+            for f in os.listdir(self.path()):
+                f_path = os.path.join(self.path(), f)
+                image = NumPyImage()
+                image.set_path(f_path)
+                if image.load():
+                    self._images.append(image)
+            return True
+        return False

mosamatic2/core/data/numpyimage.py

@@ -0,0 +1,13 @@
+from mosamatic2.core.data.filedata import FileData
+from mosamatic2.core.utils import (
+    is_numpy,
+    load_numpy_array,
+)
+
+class NumPyImage(FileData):
+    def load(self):
+        if self.path():
+            if is_numpy(self.path()):
+                self.set_object(load_numpy_array(self.path()))
+                return True
+        return False

mosamatic2/core/tasks/__init__.py

@@ -2,4 +2,5 @@ from mosamatic2.core.tasks.rescaledicomimagestask.rescaledicomimagestask import
 from mosamatic2.core.tasks.segmentmusclefatl3tensorflowtask.segmentmusclefatl3tensorflowtask import SegmentMuscleFatL3TensorFlowTask
 from mosamatic2.core.tasks.calculatescorestask.calculatescorestask import CalculateScoresTask
 from mosamatic2.core.tasks.createpngsfromsegmentationstask.createpngsfromsegmentationstask import CreatePngsFromSegmentationsTask
-from mosamatic2.core.tasks.dicom2niftitask.dicom2niftitask import Dicom2NiftiTask
+from mosamatic2.core.tasks.dicom2niftitask.dicom2niftitask import Dicom2NiftiTask
+from mosamatic2.core.tasks.selectslicefromscanstask.selectslicefromscanstask import SelectSliceFromScansTask

mosamatic2/core/tasks/calculatescorestask/calculatescorestask.py

@@ -4,10 +4,8 @@ import pandas as pd
 from mosamatic2.core.tasks.task import Task
 from mosamatic2.core.managers.logmanager import LogManager
 from mosamatic2.core.data.multidicomimage import MultiDicomImage
+from mosamatic2.core.data.numpyimage import NumPyImage
 from mosamatic2.core.utils import (
-    is_dicom, 
-    load_dicom,
-    is_jpeg2000_compressed,
     get_pixels_from_dicom_object,
     calculate_area,
     calculate_mean_radiation_attenuation,
@@ -41,13 +39,11 @@ class CalculateScoresTask(Task):
                 if file_type == '.seg.npy':
                     f_seg_name = f_seg_name.removesuffix(file_type)
                     if f_seg_name == f_img_name:
-                        # img_seg_pairs.append((f_img_path, f_seg_path))
                         img_seg_pairs.append((image, f_seg_path))
                 elif file_type == '.tag':
                     f_seg_name = f_seg_name.removesuffix(file_type).removesuffix('.dcm')
                     f_img_name = f_img_name.removesuffix('.dcm')
                     if f_seg_name == f_img_name:
-                        # img_seg_pairs.append((f_img_path, f_seg_path))
                         img_seg_pairs.append((image, f_seg_path))
                 else:
                     raise RuntimeError('Unknown file type')
@@ -78,7 +74,11 @@ class CalculateScoresTask(Task):
 
     def load_segmentation(self, f, file_type='npy'):
         if file_type == 'npy':
-            return np.load(f)
+            segmentation = NumPyImage()
+            segmentation.set_path(f)
+            if segmentation.load():
+                return segmentation.object()
+            LOG.error(f'Could not load segmentation file {f}')
         if file_type == 'tag':
             pixels = get_pixels_from_tag_file(f)
             try:

mosamatic2/core/tasks/segmentmusclefatl3tensorflowtask/segmentmusclefatl3tensorflowtask.py

@@ -28,7 +28,7 @@ class SegmentMuscleFatL3TensorFlowTask(Task):
     def __init__(self, inputs, params, output, overwrite=True):
         super(SegmentMuscleFatL3TensorFlowTask, self).__init__(inputs, params, output, overwrite)
 
-    def load_images(self):
+    def load_images(self):        
         image_data = MultiDicomImage()
         image_data.set_path(self.input('images'))
         if image_data.load():
@@ -55,7 +55,6 @@ class SegmentMuscleFatL3TensorFlowTask(Task):
                     import tensorflow as tf
                     tfLoaded = True
                 with tempfile.TemporaryDirectory() as model_dir_unzipped:
-                # model_dir_unzipped = os.path.join(os.path.split(f_path)[0], 'model_unzipped')
                     os.makedirs(model_dir_unzipped, exist_ok=True)
                     with zipfile.ZipFile(f_path) as zipObj:
                         zipObj.extractall(path=model_dir_unzipped)
@@ -65,7 +64,6 @@ class SegmentMuscleFatL3TensorFlowTask(Task):
                     import tensorflow as tf
                     tfLoaded = True
                 with tempfile.TemporaryDirectory() as contour_model_dir_unzipped:
-                # contour_model_dir_unzipped = os.path.join(os.path.split(f_path)[0], 'contour_model_unzipped')
                     os.makedirs(contour_model_dir_unzipped, exist_ok=True)
                     with zipfile.ZipFile(f_path) as zipObj:
                         zipObj.extractall(path=contour_model_dir_unzipped)

mosamatic2/core/tasks/selectslicefromscanstask/selectslicefromscanstask.py

@@ -0,0 +1,109 @@
+import os
+import math
+import tempfile
+import shutil
+import nibabel as nib
+import numpy as np
+from totalsegmentator.python_api import totalsegmentator
+from mosamatic2.core.tasks.task import Task
+from mosamatic2.core.managers.logmanager import LogManager
+from mosamatic2.core.utils import load_dicom
+
+LOG = LogManager()
+
+TOTAL_SEGMENTATOR_OUTPUT_DIR = os.path.join(tempfile.gettempdir(), 'total_segmentator_output')
+TOTAL_SEGMENTATOR_TASK = 'total'
+Z_DELTA_OFFSETS = {
+    'vertebrae_L3': 0.333,
+    'vertebrae_T4': 0.5,
+}
+
+
+class SelectSliceFromScansTask(Task):
+    INPUTS = ['scans']
+    PARAMS = ['vertebra']
+
+    def __init__(self, inputs, params, output, overwrite):
+        super(SelectSliceFromScansTask, self).__init__(inputs, params, output, overwrite)
+
+    def load_scan_dirs(self):
+        scan_dirs = []
+        for d in os.listdir(self.input('scans')):
+            scan_dir = os.path.join(self.input('scans'), d)
+            if os.path.isdir(scan_dir):
+                scan_dirs.append(scan_dir)
+        return scan_dirs
+    
+    def extract_masks(self, scan_dir):
+        os.makedirs(TOTAL_SEGMENTATOR_OUTPUT_DIR, exist_ok=True)
+        totalsegmentator(input=scan_dir, output=TOTAL_SEGMENTATOR_OUTPUT_DIR, fast=True)
+        # os.system(f'TotalSegmentator -i {scan_dir} -o {TOTAL_SEGMENTATOR_OUTPUT_DIR} --fast')
+
+    def delete_total_segmentator_output(self):
+        if os.path.exists(TOTAL_SEGMENTATOR_OUTPUT_DIR):
+            shutil.rmtree(TOTAL_SEGMENTATOR_OUTPUT_DIR)
+
+    def get_z_delta_offset_for_mask(self, mask_name):
+        if mask_name not in Z_DELTA_OFFSETS.keys():
+            return None
+        return Z_DELTA_OFFSETS[mask_name]
+
+    def find_slice(self, scan_dir, vertebra):
+        if vertebra == 'L3':
+            vertebral_level = 'vertebrae_L3'
+        elif vertebra == 'T4':
+            vertebral_level = 'vertebrae_T4'
+        else:
+            raise RuntimeError(f'Unknown vertbra {vertebra}. Options are "L3" and "T4"')
+        # Find Z-positions DICOM images
+        z_positions = {}
+        for f in os.listdir(scan_dir):
+            f_path = os.path.join(scan_dir, f)
+            p = load_dicom(f_path, stop_before_pixels=True)
+            if p is not None:
+                z_positions[p.ImagePositionPatient[2]] = f_path
+        # Find Z-position L3 image
+        mask_file = os.path.join(TOTAL_SEGMENTATOR_OUTPUT_DIR, f'{vertebral_level}.nii.gz')
+        mask_obj = nib.load(mask_file)
+        mask = mask_obj.get_fdata()
+        affine_transform = mask_obj.affine
+        indexes = np.array(np.where(mask == 1))
+        index_min = indexes.min(axis=1)
+        index_max = indexes.max(axis=1)
+        world_min = nib.affines.apply_affine(affine_transform, index_min)
+        world_max = nib.affines.apply_affine(affine_transform, index_max)
+        z_direction = affine_transform[:3, 2][2]
+        z_sign = math.copysign(1, z_direction)
+        z_delta_offset = self.get_z_delta_offset_for_mask(vertebral_level)
+        if z_delta_offset is None:
+            return None
+        z_delta = 0.333 * abs(world_max[2] - world_min[2]) # This needs to be vertebra-specific perhaps
+        z_l3 = world_max[2] - z_sign * z_delta
+        # Find closest L3 image in DICOM set
+        positions = sorted(z_positions.keys())
+        closest_file = None
+        for z1, z2 in zip(positions[:-1], positions[1:]):
+            if min(z1, z2) <= z_l3 <= max(z1, z2):
+                closest_z = min(z_positions.keys(), key=lambda z: abs(z - z_l3))
+                closest_file = z_positions[closest_z]
+                LOG.info(f'Closest image: {closest_file}')
+                break
+        return closest_file
+
+    def run(self):
+        scan_dirs = self.load_scan_dirs()
+        vertebra = self.param('vertebra')
+        nr_steps = len(scan_dirs)
+        for step in range(nr_steps):
+            scan_dir = scan_dirs[step]
+            scan_name = os.path.split(scan_dir)[1]
+            self.extract_masks(scan_dir)
+            file_path = self.find_slice(scan_dir, vertebra)
+            if file_path is not None:
+                extension = '' if file_path.endswith('.dcm') else '.dcm'
+                target_file_path = os.path.join(self.output(), vertebra + '_' + scan_name + extension)
+                shutil.copyfile(file_path, target_file_path)
+            else:
+                LOG.error(f'Could not find slice for vertebral level: {vertebra}')
+            self.delete_total_segmentator_output()
+            self.set_progress(step, nr_steps)

mosamatic2/core/utils.py

@@ -119,10 +119,18 @@ def is_numpy_array(value):
     return isinstance(value, np.array)
 
 
-def load_numpy_array(f):
+def is_numpy(f):
     try:
+        np.load(f)
+        return True
+    except:
+        return False
+
+
+def load_numpy_array(f):
+    if is_numpy(f):
         return np.load(f)
-    except Exception as e:
+    else:
         return None
 
 

mosamatic2/server.py

@@ -6,6 +6,7 @@ from mosamatic2.core.tasks import SegmentMuscleFatL3TensorFlowTask
 from mosamatic2.core.tasks import CalculateScoresTask
 from mosamatic2.core.tasks import CreatePngsFromSegmentationsTask
 from mosamatic2.core.tasks import Dicom2NiftiTask
+from mosamatic2.core.tasks import SelectSliceFromScansTask
 
 app = Flask(__name__)
 
@@ -105,6 +106,22 @@ def run_dicom2nifti():
     return 'PASSED'
 
 
+@app.route('/selectslicefromscans')
+def run_selectslicefromscans():
+    scans = request.args.get('scans')
+    vertebra = request.args.get('vertebra')
+    output = request.args.get('output')
+    overwrite = request.args.get('overwrite', default=True, type=bool)
+    task = SelectSliceFromScansTask(
+        inputs={'scans': scans},
+        params={'vertebra': vertebra},
+        output=output,
+        overwrite=overwrite,
+    )
+    task.run()
+    return 'PASSED'
+
+
 def main():
     parser = argparse.ArgumentParser()
     parser.add_argument('--port', type=int, default=constants.MOSAMATIC2_SERVER_PORT)

mosamatic2/ui/mainwindow.py

@@ -19,6 +19,7 @@ from mosamatic2.ui.widgets.panels.tasks.segmentmusclefatl3tensorflowtaskpanel im
 from mosamatic2.ui.widgets.panels.tasks.createpngsfromsegmentationstaskpanel import CreatePngsFromSegmentationsTaskPanel
 from mosamatic2.ui.widgets.panels.tasks.calculatescorestaskpanel import CalculateScoresTaskPanel
 from mosamatic2.ui.widgets.panels.tasks.dicom2niftitaskpanel import Dicom2NiftiTaskPanel
+from mosamatic2.ui.widgets.panels.tasks.selectslicefromscanstaskpanel import SelectSliceFromScansTaskPanel
 from mosamatic2.ui.widgets.panels.pipelines.defaultpipelinepanel import DefaultPipelinePanel
 
 LOG = LogManager()
@@ -36,6 +37,7 @@ class MainWindow(QMainWindow):
         self._create_pngs_from_segmentations_task_panel = None
         self._calculate_scores_task_panel = None
         self._dicom2nifti_task_panel = None
+        self._select_slice_from_scans_task_panel = None
         self._default_pipeline_panel = None
         self.init_window()
 
@@ -75,12 +77,15 @@ class MainWindow(QMainWindow):
         create_pngs_from_segmentations_task_action.triggered.connect(self.handle_create_pngs_from_segmentations_task_action)
         dicom2nifti_task_action = QAction('Dicom2NiftiTask', self)
         dicom2nifti_task_action.triggered.connect(self.handle_dicom2nifti_task_action)
+        select_slice_from_scans_task_action = QAction('SelectSliceFromScansTask', self)
+        select_slice_from_scans_task_action.triggered.connect(self.handle_select_slice_from_scans_task_action)
         tasks_menu = self.menuBar().addMenu('Tasks')
         tasks_menu.addAction(rescale_dicom_images_task_action)
         tasks_menu.addAction(segment_muscle_fat_l3_tensorflow_task_action)
         tasks_menu.addAction(calculate_scores_task_action)
         tasks_menu.addAction(create_pngs_from_segmentations_task_action)
         tasks_menu.addAction(dicom2nifti_task_action)
+        tasks_menu.addAction(select_slice_from_scans_task_action)
 
     def init_pipelines_menu(self):
         default_pipeline_action = QAction('DefaultPipeline', self)
@@ -106,6 +111,7 @@ class MainWindow(QMainWindow):
             self._main_panel.add_panel(self.create_pngs_from_segmentations_task_panel(), 'createpngsfromsegmentationstaskpanel')
             self._main_panel.add_panel(self.calculate_scores_task_panel(), 'calculatescorestaskpanel')
             self._main_panel.add_panel(self.dicom2nifti_task_panel(), 'dicom2niftitaskpanel')
+            self._main_panel.add_panel(self.select_slice_from_scans_task_panel(), 'selectslicefromscanstaskpanel')
             self._main_panel.add_panel(self.default_pipeline_panel(), 'defaultpipelinepanel')
             self._main_panel.select_panel('defaultpipelinepanel')
         return self._main_panel
@@ -142,6 +148,11 @@ class MainWindow(QMainWindow):
         if not self._dicom2nifti_task_panel:
             self._dicom2nifti_task_panel = Dicom2NiftiTaskPanel()
         return self._dicom2nifti_task_panel
+
+    def select_slice_from_scans_task_panel(self):
+        if not self._select_slice_from_scans_task_panel:
+            self._select_slice_from_scans_task_panel = SelectSliceFromScansTaskPanel()
+        return self._select_slice_from_scans_task_panel
     
     def default_pipeline_panel(self):
         if not self._default_pipeline_panel:
@@ -170,6 +181,9 @@ class MainWindow(QMainWindow):
     def handle_dicom2nifti_task_action(self):
         self.main_panel().select_panel('dicom2niftitaskpanel')
 
+    def handle_select_slice_from_scans_task_action(self):
+        self.main_panel().select_panel('selectslicefromscanstaskpanel')
+
     def handle_default_pipeline_action(self):
         self.main_panel().select_panel('defaultpipelinepanel')
 
@@ -184,6 +198,7 @@ class MainWindow(QMainWindow):
         self.create_pngs_from_segmentations_task_panel().save_inputs_and_parameters()
         self.calculate_scores_task_panel().save_inputs_and_parameters()
         self.dicom2nifti_task_panel().save_inputs_and_parameters()
+        self.select_slice_from_scans_task_panel().save_inputs_and_parameters()
         self.default_pipeline_panel().save_inputs_and_parameters()
         return super().closeEvent(event)
 

mosamatic2/ui/resources/VERSION

@@ -1 +1 @@
-2.0.4 
+2.0.5 

mosamatic2/ui/widgets/panels/tasks/{selectslicefromscantaskpanel.py → selectslicefromscanstaskpanel.py}

@@ -3,7 +3,7 @@ import os
 from PySide6.QtWidgets import (
     QLineEdit,
     QCheckBox,
-    QSpinBox,
+    QComboBox,
     QHBoxLayout,
     QVBoxLayout,
     QFormLayout,
@@ -16,26 +16,26 @@ from PySide6.QtCore import (
     Slot,
 )
 
-from mosamaticdesktop.core.utils.logmanager import LogManager
-from mosamaticdesktop.ui.panels.taskpanel import TaskPanel
-from mosamaticdesktop.ui.settings import Settings
-from mosamaticdesktop.ui.utils import is_macos
-from mosamaticdesktop.ui.worker import Worker
-
-from mosamatic.tasks import SelectSliceFromScanTask
+from mosamatic2.core.managers.logmanager import LogManager
+from mosamatic2.ui.widgets.panels.taskpanel import TaskPanel
+from mosamatic2.ui.settings import Settings
+from mosamatic2.ui.utils import is_macos
+from mosamatic2.ui.worker import Worker
+from mosamatic2.core.tasks import SelectSliceFromScansTask
 
 LOG = LogManager()
 
-PANEL_TITLE = 'Automatically select L3 or T4 slice from full CT scan'
-PANEL_NAME = 'selectslicefromscantaskpanel'
+PANEL_TITLE = 'Automatically select L3 slice from full CT scan'
+PANEL_NAME = 'selectslicefromscanstaskpanel'
 
 
-class SelectSliceFromScanTaskPanel(TaskPanel):
+class SelectSliceFromScansTaskPanel(TaskPanel):
     def __init__(self):
-        super(SelectSliceFromScanTaskPanel, self).__init__()
+        super(SelectSliceFromScansTaskPanel, self).__init__()
         self.set_title(PANEL_TITLE)
         self._scans_dir_line_edit = None
         self._scans_dir_select_button = None
+        self._vertebra_combobox = None
         self._output_dir_line_edit = None
         self._output_dir_select_button = None
         self._overwrite_checkbox = None
@@ -58,6 +58,13 @@ class SelectSliceFromScanTaskPanel(TaskPanel):
             self._scans_dir_select_button.clicked.connect(self.handle_scans_dir_select_button)
         return self._scans_dir_select_button
     
+    def vertebra_combobox(self):
+        if not self._vertebra_combobox:
+            self._vertebra_combobox = QComboBox()
+            self._vertebra_combobox.addItems(['L3'])
+            self._vertebra_combobox.setCurrentText(self.settings().get(f'{PANEL_NAME}/vertebra'))
+        return self._vertebra_combobox
+    
     def output_dir_line_edit(self):
         if not self._output_dir_line_edit:
             self._output_dir_line_edit = QLineEdit(self.settings().get(f'{PANEL_NAME}/output_dir'))
@@ -101,6 +108,7 @@ class SelectSliceFromScanTaskPanel(TaskPanel):
         output_dir_layout.addWidget(self.output_dir_line_edit())
         output_dir_layout.addWidget(self.output_dir_select_button())
         self.form_layout().addRow('Scans directory', scans_dir_layout)
+        # self.form_layout().addRow('Vertebra', self.vertebra_combobox())
         self.form_layout().addRow('Output directory', output_dir_layout)
         self.form_layout().addRow('Overwrite', self.overwrite_checkbox())
         layout = QVBoxLayout()
@@ -134,11 +142,11 @@ class SelectSliceFromScanTaskPanel(TaskPanel):
             LOG.info('Running task...')
             self.run_task_button().setEnabled(False)
             self.save_inputs_and_parameters()
-            self._task = SelectSliceFromScanTask(
-                self.scans_dir_line_edit().text(), 
-                self.output_dir_line_edit().text(), 
-                'vertebrae_L3',
-                self.overwrite_checkbox().isChecked()
+            self._task = SelectSliceFromScansTask(
+                inputs={'scans': self.scans_dir_line_edit().text()},
+                params={'vertebra': 'L3'},
+                output=self.output_dir_line_edit().text(),
+                overwrite=self.overwrite_checkbox().isChecked(),
             )
             self._worker = Worker(self._task)
             self._thread = QThread()
@@ -180,5 +188,6 @@ class SelectSliceFromScanTaskPanel(TaskPanel):
     
     def save_inputs_and_parameters(self):
         self.settings().set(f'{PANEL_NAME}/scans_dir', self.scans_dir_line_edit().text())
+        self.settings().set(f'{PANEL_NAME}/vertebra', self.vertebra_combobox().currentText())
         self.settings().set(f'{PANEL_NAME}/output_dir', self.output_dir_line_edit().text())
         self.settings().set(f'{PANEL_NAME}/overwrite', self.overwrite_checkbox().isChecked())

{mosamatic2-2.0.4.dist-info → mosamatic2-2.0.5.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: mosamatic2
-Version: 2.0.4
+Version: 2.0.5
 Summary: 
 Author: Ralph Brecheisen
 Author-email: r.brecheisen@maastrichtuniversity.nl

{mosamatic2-2.0.4.dist-info → mosamatic2-2.0.5.dist-info}/RECORD

@@ -1,13 +1,14 @@
 models.py,sha256=Kx6oWKt7IpTTxrhBDrX61X-ZX12J7yPkJFuhVDsDHoQ,8807
 mosamatic2/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 mosamatic2/app.py,sha256=RIUa5tvMYFcmEII4xZPLZZdx9dXWqBvwkxkl_R97Jkw,860
-mosamatic2/cli.py,sha256=bVd2zf7CYpVIc-zRweWgRO7jOPZi68JhnnDBuq-L43E,1201
+mosamatic2/cli.py,sha256=RTogXBXiYW_EDTmnX7TXx2T22Yth_1AwLFiJP2IUDMc,1289
 mosamatic2/commands/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 mosamatic2/commands/calculatescores.py,sha256=Lb8Q8L2yq7Tt6VBJ6_lltRuldrev_pac6fcgF-xzZyE,1984
 mosamatic2/commands/createpngsfromsegmentations.py,sha256=uUAQJVTqOkBCfENzi21RBNYvf6_nuesx1MeR3j_-7dM,1682
 mosamatic2/commands/dicom2nifti.py,sha256=4hMvglxdOid7p7P_Jc-Tudsf93JRMlYaQsEQctLq2d4,1015
 mosamatic2/commands/rescaledicomimages.py,sha256=25QdCzB5s0sRwkTb3o5zco2bIwy6LttNf7i97kGBDYQ,1280
 mosamatic2/commands/segmentmusclefatl3tensorflow.py,sha256=CdScmA_EQicaN4GY5bBUOYwfhDPqy9om2sxY3WrtmM0,1424
+mosamatic2/commands/selectslicefromscans.py,sha256=3398PM2uBcxF6wpb0-c-Itp_qxoAxBf0SE2nDDI3Ct4,1715
 mosamatic2/constants.py,sha256=MVYMwO-x2jQSN37o3zNkRseVQ1nYRA3mLv3v_Q0Mlds,1284
 mosamatic2/core/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 mosamatic2/core/data/__init__.py,sha256=j9iGqUTJlGF0N0gPrzzpe_Dhv0Bj9c6FdQ1g7U-_j2g,298
@@ -16,6 +17,8 @@ mosamatic2/core/data/dicomimageseries.py,sha256=OZkNi15crL8nEA-PGYsM0k9NMi2mMHRv
 mosamatic2/core/data/dixonseries.py,sha256=kq9fy65MSM2XwiScqp7b3rQ09JmpyGwbG6ldZsuPRrM,516
 mosamatic2/core/data/filedata.py,sha256=hCnpizGqOpxzIADJkDS2_NSmKVLL1u49TYjSJE5UXQo,515
 mosamatic2/core/data/multidicomimage.py,sha256=cdd0H4Dq49h7NLKBx51_h_HZVnH7-reu48PY8m6tXwU,1034
+mosamatic2/core/data/multinumpyimage.py,sha256=MU5_ec00vx9_JhkibgkcgyGFfaj69STy85ifa8MNsSc,753
+mosamatic2/core/data/numpyimage.py,sha256=ZbJ9JRAVOgvXzTZ6kscEuGZJqKT-MPdt42IAlmh3iXo,359
 mosamatic2/core/managers/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 mosamatic2/core/managers/logmanager.py,sha256=NEaXvhl0aILjBbK710GaWanVuuNvB51HpHhE5rgYvng,1391
 mosamatic2/core/managers/logmanagerlistener.py,sha256=Gaig07yjBnyQq9I8sN85olTEeDCDyCFQnEJdwzvmgvc,99
@@ -23,9 +26,9 @@ mosamatic2/core/pipelines/__init__.py,sha256=MXoylu34AlcyudJpjuLa4B7K9gZmxLnZGDa
 mosamatic2/core/pipelines/defaultpipeline.py,sha256=EUblRhLZXhslz9bdBmx8v9yjpgOdsHviqQZ5elK8gNQ,2901
 mosamatic2/core/pipelines/pipeline.py,sha256=mRxKXLKwgKDpc8R9mCI6gDKGJ2lKVxRQ__Sf0Mfn_Qc,384
 mosamatic2/core/singleton.py,sha256=FV0k_LlOCmFhlWN6gf1c2x7YXWyd8-7DsIMvOKrI6NY,224
-mosamatic2/core/tasks/__init__.py,sha256=3bWBH2Q5X44_k7s3wei7AHm86MLEqoySyxSN7UKuqRs,545
+mosamatic2/core/tasks/__init__.py,sha256=az_IatVtpT7M2AiOl7ez6EW_oR7_n9GtC1nL0_mYwZk,655
 mosamatic2/core/tasks/calculatescorestask/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-mosamatic2/core/tasks/calculatescorestask/calculatescorestask.py,sha256=1CUJK9IsePBbEuF3-0LnZp1sWMxPDiJaNy0a-8DaTUE,6576
+mosamatic2/core/tasks/calculatescorestask/calculatescorestask.py,sha256=jEbogrEnEtMbmfQ8KOGSUdGMhdeX69Bv01-82faRIlY,6617
 mosamatic2/core/tasks/createpngsfromsegmentationstask/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 mosamatic2/core/tasks/createpngsfromsegmentationstask/createpngsfromsegmentationstask.py,sha256=1UpOsp1CH668BQ0g4tALou_tFgisC306VcvqOKSDuTo,1884
 mosamatic2/core/tasks/dicom2niftitask/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
@@ -34,17 +37,19 @@ mosamatic2/core/tasks/rescaledicomimagestask/__init__.py,sha256=47DEQpj8HBSa-_TI
 mosamatic2/core/tasks/rescaledicomimagestask/rescaledicomimagestask.py,sha256=vGSpMJoXFtE-IHGxTEO9DMkerMJcfG5r9tXgtvkxm6Y,3053
 mosamatic2/core/tasks/segmentmusclefatl3tensorflowtask/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 mosamatic2/core/tasks/segmentmusclefatl3tensorflowtask/paramloader.py,sha256=VxTCOYK_1VRAG83P-ulm0LPvqXI-0iT5BCr0Rdr7MWg,900
-mosamatic2/core/tasks/segmentmusclefatl3tensorflowtask/segmentmusclefatl3tensorflowtask.py,sha256=bDRG9phYSF4o4mFsdJhyJ613CqOixFPTd5yL2gBFSPY,5351
+mosamatic2/core/tasks/segmentmusclefatl3tensorflowtask/segmentmusclefatl3tensorflowtask.py,sha256=E3GHzGtl-v_qIdLjSy3cbDnf8vfWK6Bax3RAZcKCqXg,5149
+mosamatic2/core/tasks/selectslicefromscanstask/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+mosamatic2/core/tasks/selectslicefromscanstask/selectslicefromscanstask.py,sha256=1a_QHcIdS4VB39cNEBniaaDzUx3k8DcR0wCtdr8G1d0,4586
 mosamatic2/core/tasks/task.py,sha256=APPnid6dpSGkPuDqU1vm2RIMR5vkpvbP1CPHUMjympg,1691
-mosamatic2/core/utils.py,sha256=RHBcvrUshese6NG_OigO4ZGGR3iCOmmeDJjgPnLpaBo,10753
-mosamatic2/server.py,sha256=UDALqk_TZ5Vz7ZOjzhV2-7-kHIPM-_OKXFV7lf2R1t0,3617
+mosamatic2/core/utils.py,sha256=zh44FNOWxNKuZ4FcM7VIfMvdlzOr4AA8_PJ1r-6_83k,10855
+mosamatic2/server.py,sha256=jiUc7-aFaDMSY2hg6BhMC9C2FU4V2OrMBSoiIAP8Zcs,4158
 mosamatic2/ui/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-mosamatic2/ui/mainwindow.py,sha256=efmZAC5vYEdCg5bw1ixrcACzYN2xe15xyrNwnXersyI,10100
+mosamatic2/ui/mainwindow.py,sha256=k5yxHTjwvpwXXWx94sJ45q1BcH-06B9I7Qnqfd1Zbr8,11119
 mosamatic2/ui/resources/icons/mosamatic2.icns,sha256=OfhC-diJTIgaNMOezxKKilGsY7mRkaGdU5dGr0MOjIA,2994125
 mosamatic2/ui/resources/icons/mosamatic2.ico,sha256=ySD3RYluHK3pgS0Eas7eKrVk_AskdLQ4qs_IT-wNhq4,12229
 mosamatic2/ui/resources/icons/spinner.gif,sha256=rvaac6GUZauHSPFSOLWr0RmLfjmtZih2Q8knQ2WP3Po,16240
 mosamatic2/ui/resources/images/body-composition.jpg,sha256=KD-BudbXwThB4lJOZZN-ad5-TZRaaZ5cKTH0Ar1TOZs,21227
-mosamatic2/ui/resources/VERSION,sha256=CxHeq1Bk1fry3nkuFpVDAvnw27-hr4X2onfG3ZD1M6U,8
+mosamatic2/ui/resources/VERSION,sha256=I2Qe9unc-Hx4iKTROuxBK26Iy4oCUcMjfIY6KzAtKAY,8
 mosamatic2/ui/settings.py,sha256=YEVHYJIfNsqMO3v1pjzgh7Pih9GGoUX7S9s8S-sBNUk,2121
 mosamatic2/ui/utils.py,sha256=6bbPIrh4RJ_yhQKNZrgPbL4XeUEogjIjbk_e5c3QS5g,853
 mosamatic2/ui/widgets/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
@@ -65,10 +70,10 @@ mosamatic2/ui/widgets/panels/tasks/createpngsfromsegmentationstaskpanel.py,sha25
 mosamatic2/ui/widgets/panels/tasks/dicom2niftitaskpanel.py,sha256=mb5LxDuhkwA08LcHGigul7xmef_wTGrmr-DZ72O_joc,7474
 mosamatic2/ui/widgets/panels/tasks/rescaledicomimagestaskpanel.py,sha256=_K-XM4zvzq_AMugLqN4p87iEsfX594qYSVnfjNlmAfU,7591
 mosamatic2/ui/widgets/panels/tasks/segmentmusclefatl3tensorflowtaskpanel.py,sha256=pLEpUXx_vi6PfOykAwp95m5yLRrIvdF_XaupHBOW31U,9624
-mosamatic2/ui/widgets/panels/tasks/selectslicefromscantaskpanel.py,sha256=Q5txmTPssKEQpy_yc-iKCN22K5h3B_7WzPP885uP79Q,7506
+mosamatic2/ui/widgets/panels/tasks/selectslicefromscanstaskpanel.py,sha256=ysIoX_IHTPBfEhJQkdM54S3ZL24IAeA8A_ldzgr8BwU,8071
 mosamatic2/ui/widgets/splashscreen.py,sha256=MS-OczOWfwwEQNQd-JWe9_Mh57css0cSQgbu973rwQo,4056
 mosamatic2/ui/worker.py,sha256=v7e3gq7MUudgpB1BJW-P7j5wurzu6-HG5m7I6WHgJp0,699
-mosamatic2-2.0.4.dist-info/entry_points.txt,sha256=MCUpKkgbej1clgp8EqlLQGs0BIKwGPcBPiVWLfGz9Gw,126
-mosamatic2-2.0.4.dist-info/METADATA,sha256=89x_i4kNHoZwtq3qBwy5QPOHC-ljsHJsntQDCh5km50,1402
-mosamatic2-2.0.4.dist-info/WHEEL,sha256=b4K_helf-jlQoXBBETfwnf4B04YC67LOev0jo4fX5m8,88
-mosamatic2-2.0.4.dist-info/RECORD,,
+mosamatic2-2.0.5.dist-info/entry_points.txt,sha256=MCUpKkgbej1clgp8EqlLQGs0BIKwGPcBPiVWLfGz9Gw,126
+mosamatic2-2.0.5.dist-info/METADATA,sha256=ODtnfcQqYDrCk6jNK9J1vr7eV3XgUi1yT33ctWPzQ_E,1402
+mosamatic2-2.0.5.dist-info/WHEEL,sha256=b4K_helf-jlQoXBBETfwnf4B04YC67LOev0jo4fX5m8,88
+mosamatic2-2.0.5.dist-info/RECORD,,