I think much of what you’re trying to do here can probably be done using tckedit: if any streamlines are selected when specifying an inclusion ROI, then these streamlines at least intersect with that ROI. You can also use the -ends_only option if you’re interested about whether the streamlines end in those regions.
That’s not something diffusion MRI can tell you. Diffusion is an inherently symmetric process, there’s no way to tell afferent from efferent…
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?