# Tckedit or tckgen for tract dissection
**URL:** https://community.mrtrix.org/t/tckedit-or-tckgen-for-tract-dissection/1654
**Category:** Uncategorized
**Created:** [May 22, 2018, 12:55pm UTC](https://community.mrtrix.org/t/tckedit-or-tckgen-for-tract-dissection/1654 "2018-05-22T12:55:31Z")
**Posts on this page:** 4
**Page:** 1
### Author:  [@anege](https://community.mrtrix.org/u/anege)
#### Post date: [May 22, 2018, 12:55pm UTC](https://community.mrtrix.org/t/tckedit-or-tckgen-for-tract-dissection/1654/1 "2018-05-22T12:55:32Z")
Dear experts,
I am trying to dissect fibers passing through two ROIs.
In order to do so, I am using
> tckedit WB\_tractogram.tck -include ROI1 -include ROI2 OUTPUT.tck -force
where the WholeBrain (WB) tractogram is a 20-Million-fiber whole-brain tractogram.
• My first question is **whether using the 20M-fiber WB tractogram is a good approach or I should decrease the number of fibers of the tractogram from 20M to 2M** and then run tckedit. I opted for the 20M-fiber tractogram since I get more fibers between the two ROIs, but I am not sure whether this is right. Since the threshold to consider that a tract exists at an individual level is usually set at 10 fibers, I am concerned that my results (the existence/not existence of that tract) are biased by the number of fibers of my WB tractogram. In other words, it is easier to find 10 fibers that pass through my ROIs with a 20M WB tractogram than with a 2M tractogram.
• My second question is **whether I should use tckedit (with the WB tractogram) or tckgen (just considering my ROIs)**.
I apologize if my questions are too basic and thank you very much for the help you provide everyday.
Best
---
### Author:  [@ThijsDhollander](https://community.mrtrix.org/u/ThijsDhollander)
#### Post date: [May 23, 2018, 8:01am UTC](https://community.mrtrix.org/t/tckedit-or-tckgen-for-tract-dissection/1654/2 "2018-05-23T08:01:11Z")
Hi @anege,
I would use “true” (debatable) targeted tractography, i.e. `tckgen` directly every single day without a doubt for that scenario, especially if your aim is mostly “_dissecting_” (i.e. segmentation, potentially to create a ROI for further quantification, but not even per se). You’ll have much faster access and control of parameters which are only available during the `tckgen` step itself, such as step size, curvature constraints, even fancy things like the iFOD2 power parameter, etc… A lot of targeted tractography scenarios involve some level of supervision (which, I reckon, is what ultimately enables them to be reliable), and with that, a bit of trial and error. The best results are those that are semi-automated (i.e. semi-driven by the data), but also consciously inspected by someone with anatomical knowledge of what structure they are after. And the latter bit may lead to feedback to improve the parameters, both direct streamline tracking parameters (e.g. curvature and the like) as well as ROI induced constraints (inclusion, exclusion, masks, etc…).
Cheers,
Thijs
---
### Author:  [@LiuYuchen](https://community.mrtrix.org/u/LiuYuchen)
#### Post date: [November 13, 2019, 11:33am UTC](https://community.mrtrix.org/t/tckedit-or-tckgen-for-tract-dissection/1654/3 "2019-11-13T11:33:59Z")
If it be possible to generate a .nii file include the voxels of the inclusion ROI which the tracts passing through?
---
### Author:  [@ThijsDhollander](https://community.mrtrix.org/u/ThijsDhollander)
#### Post date: [November 14, 2019, 12:32am UTC](https://community.mrtrix.org/t/tckedit-or-tckgen-for-tract-dissection/1654/4 "2019-11-14T00:32:36Z")
For future reference, @LiuYuchen’s question was answered over here:
> [@Estimating tracts for a specific region](https://community.mrtrix.org/t/estimating-tracts-for-a-specific-region/885/4):
>
> If I understand what you mean correctly, this should do it:
>
> ```auto
> tckmap your_tracks.tck -template your_inclusion_ROI.mif - | mrcalc - your_inclusion_ROI.mif -mult final_map.mif
>
> ```
>
> Essentially: use `tckmap` to create an image of streamline counts for your tracks in ROI space, then multiply that by your ROI image itself. Is this what you were after?