UKF Tractography Output Issues

Operating system: Windows 10
Slicer version: 4.10.2
Expected behavior: UKF Tractography
Actual behavior: Empty output file from UKF Tractography

Hi all,

I am trying to run a UKF Tractography on Windows 10, and am having problems. Namely, the program will run for somewhere between 4-20 hours, tell me ‘Status Completed’, but then output a completely empty bundle file.

It has worked once correctly, whereby after about 2 hours I had a complete output. The other times it never worked, despite me keeping all settings and parameters the same…

Here is the log file output from one of the times it did not work:

——————————————————————————————————

[DEBUG][Qt] 14.09.2020 13:12:55 [] (unknown:0) - UKF Tractography standard output:

Using the 2T simple model. Setting the default parameters accordingly:

“*”: set by user

“-”: default setting

  • stoppingFA: 0.12

  • seedingThreshold: 0.18

  • Qm: 0.001

  • Ql: 50

  • Rs: 0.02

  • stepLength: 0.3

  • recordLength: 2

  • stoppingThreshold: 0.1

  • seedsPerVoxel: 1

Found 28 cores on your system.

Running tractography with 28 thread(s).

Number of non-zero gradients: 64

Number of zero gradients: 1

Permuting the axis order to: 0 1 2 3

Resizing the data to: 64 220 220 144

Computing the baseline image

Dividing the signal by baseline image

Data normalization finished!

Using 2-tensor simple model.

Branching disabled

Using unconstrained filter

nmse_avg=0

[DEBUG][Qt] 14.09.2020 13:12:55 [] (unknown:0) - UKF Tractography completed without errors

[ERROR][VTK] 14.09.2020 13:12:55 [vtkCompositeDataPipeline (0000014846A160B0)] (D:\D\S\Slicer-4102-build\VTK\Common\ExecutionModel\vtkDemandDrivenPipeline.cxx:709) - Input port 0 of algorithm vtkExtractSelectedPolyDataIds(00000148AEB3D000) has 0 connections but is not optional.

[ERROR][VTK] 14.09.2020 13:12:55 [vtkCompositeDataPipeline (0000014846A160B0)] (D:\D\S\Slicer-4102-build\VTK\Common\ExecutionModel\vtkDemandDrivenPipeline.cxx:709) - Input port 0 of algorithm vtkExtractSelectedPolyDataIds(00000148AEB3D000) has 0 connections but is not optional.

[ERROR][VTK] 14.09.2020 13:12:55 [vtkCompositeDataPipeline (0000014846A160B0)] (D:\D\S\Slicer-4102-build\VTK\Common\ExecutionModel\vtkDemandDrivenPipeline.cxx:709) - Input port 0 of algorithm vtkExtractSelectedPolyDataIds(00000148AEB3D000) has 0 connections but is not optional.

[ERROR][VTK] 14.09.2020 13:13:33 [vtkXMLDataParser (00000148B0221620)] (D:\D\S\Slicer-4102-build\VTK\IO\XMLParser\vtkXMLDataParser.cxx:495) - Error reading beginning of compression header. Read 0 of 12 bytes.

[ERROR][VTK] 14.09.2020 13:13:33 [vtkXMLDataParser (00000148B0221620)] (D:\D\S\Slicer-4102-build\VTK\IO\XMLParser\vtkXMLDataParser.cxx:847) - ReadCompressionHeader failed. Aborting read.

[ERROR][VTK] 14.09.2020 13:13:33 [vtkXMLPolyDataReader (00000148ADD18430)] (D:\D\S\Slicer-4102-build\VTK\IO\XML\vtkXMLDataReader.cxx:415) - Cannot read point data array “FA1” from PointData in piece 0. The data array in the element may be too short.

———————————————————————————————————

Could someone tell me what is going wrong?

I really appreciate your help, thanks!

Not sure if anyone got back to you on this. A couple suggestions: use the latest preview version of 4.11. If you are still having trouble can you try reproducing this with one of the datasets in the sample data. Let us know how it goes.

Thank you for your response! I did the update, and tried to run it on the sample data, which worked well.

Comparing the files from the sample and my data, I noticed the issue was coming from the brain mask I had created, which was not of data type short. This must have been the problem, as now that I have changed it things seem to be running fine.

Thanks for everything!