marcnol / pyHiM

Multiplexed DNA-FISH data analysis pipeline
GNU Lesser General Public License v3.0
4 stars 2 forks source link

Bug (localize_3d) - error when pyHiM try to segment reference RT image #201

Open XDevos opened 8 months ago

XDevos commented 8 months ago

The problem should come from ~image_3d_0 = reinterpolate_z(...) inside segment_sources_3d_file method.

parameters.json:

"barcode": {
    "segmentedObjects": {
        "3Dmethod": "thresholding",
        "reducePlanes": false
    }
}

terminal:

>>>Processing roi:[005] cycle:[RT27]<<<
$ File:scan_001_RT27_005_ROI_converted_decon_ch01.tif
$ Reduced Z-planes from 60 to 30
$ z_range used = 0-30
> Removing inhomogeneous background from 30 planes using 1 worker...
100%|█████████████████████████████████████████| 30/30 [00:01<00:00, 16.99it/s]
$ Running reference fiducial cycle: no shift applied!
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
> Segmenting 30 planes using 1 worker...
  0%|   
Traceback (most recent call last):
  File "/home/xdevos/Repositories/repo-marcnol/pyHiM/src/pyHiM.py", line 161, in <module>
    main()
  File "/home/xdevos/Repositories/repo-marcnol/pyHiM/src/pyHiM.py", line 116, in main
    pipe.segment_sources_3d(current_param, label, datam.processed_roi)
  File "/home/xdevos/Repositories/repo-marcnol/pyHiM/src/core/function_caller.py", line 275, in segment_sources_3d
    _segment_sources_3d.segment_sources_3d()
  File "/home/xdevos/Repositories/repo-marcnol/pyHiM/src/imageProcessing/segmentSources3D.py", line 580, in segment_sources_3d
    self.segment_sources_3d_in_folder()
  File "/home/xdevos/Repositories/repo-marcnol/pyHiM/src/imageProcessing/segmentSources3D.py", line 506, in segment_sources_3d_in_folder
    output_tables.append(self.segment_sources_3d_file(filename_to_process))
  File "/home/xdevos/Repositories/repo-marcnol/pyHiM/src/imageProcessing/segmentSources3D.py", line 314, in segment_sources_3d_file
    _, segmented_image_3d = self._segment_3d_volumes(image_3d_aligned)
  File "/home/xdevos/Repositories/repo-marcnol/pyHiM/src/imageProcessing/segmentSources3D.py", line 222, in _segment_3d_volumes
    binary, segmented_image_3d = _segment_3d_volumes_by_thresholding(
  File "/home/xdevos/Repositories/repo-marcnol/pyHiM/src/imageProcessing/segmentMasks.py", line 1052, in _segment_3d_volumes_by_thresholding
    image_2d_segmented = _segment_2d_image_by_thresholding(
  File "/home/xdevos/Repositories/repo-marcnol/pyHiM/src/imageProcessing/segmentMasks.py", line 937, in _segment_2d_image_by_thresholding
    if segm.nlabels <= 0:
AttributeError: 'NoneType' object has no attribute 'nlabels'