Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
  • 03:57

    skoudoro on master

    BF: Fixes #2010, by using a mor… Merge pull request #2011 from a… (compare)

  • 03:57
    skoudoro closed #2011
  • 03:57
    skoudoro closed #2010
  • 03:57
    skoudoro commented #2011
  • 03:56
    skoudoro labeled #2011
  • 03:56
    skoudoro milestoned #2011
  • 01:37
    arokem opened #2011
  • Dec 05 23:38
    Garyfallidis assigned #2010
  • Dec 05 23:38
    Garyfallidis opened #2010
  • Dec 05 20:38
    codecov[bot] commented #1965
  • Dec 05 20:38
    pep8speaks commented #1965
  • Dec 05 20:38
    codecov[bot] commented #1965
  • Dec 05 20:38
    Garyfallidis synchronize #1965
  • Dec 05 20:38

    Garyfallidis on horizon_updates_b1

    DOC: improved documentation for… (compare)

  • Dec 03 14:40
    skoudoro commented #1959
  • Dec 03 14:05
    RafaelNH commented #2009
  • Dec 03 10:44
    RafaelNH commented #2009
  • Dec 03 10:42
    RafaelNH commented #2009
  • Dec 03 10:02
    RafaelNH commented #2009
  • Dec 02 23:58
    skoudoro commented #1965
Jon Haitz Legarreta Gorroño
@jhlegarreta

@robertl1996

Does anybody know where the source code for this denoise library can be found: https://dipy.org/documentation/1.0.0./reference/dipy.denoise/. Particularly the nl means filters

The source code is found here: https://github.com/nipy/dipy/blob/master/dipy/denoise/

Jon Haitz Legarreta Gorroño
@jhlegarreta

@zlf123-marker

I have a look, but I want to display my TRK file and t1 data, I don't know what to change?Can you help me @arokem

DIPY examples/tutorials use data that is fetched directly from a few sites that host the data of interest. If you want to adapt the examples to use your own data, then, instead of using the fetch_* and read_ * commands, you will need to read your local files.

In order to read your t1, you will need to use NiBabel as in
https://nipy.org/nibabel/nibabel_images.html

something like

import nibabel as nib

t1_img =  nib.load(your_t1_filename)
t1 = t1_img.get_data()
fig = plt.figure()
plt.imshow(t1, cmap="gray")

Similarly, to show your TRK file, you can do something similar to this

from dipy.io.streamline import load_tractogram
from fury import window, actor

tractogram = load_tractogram(your_trk_filename, your_t1_filename)
streamlines = tractogram.streamlines
renderer = window.Renderer()
stream_actor = actor.line(streamlines)
renderer.add(stream_actor)
window.record(renderer)

If you want to display your streamlines and your T1 at the same time, you will need to add both actors to the same scene, adding to the above something like

t1_actor = actor.slicer(t1)
renderer.add(t1_actor)
Jakob Wasserthal
@wasserth
Hi, is there a function in dipy to resample streamlines to make all points have a specific distance (in mm) from each other? So longer streamlines would have more points than shorter streamlines.
Wade Weber
@wadeatUT
Hi all. When I convert my Freesurfer apar+aseg back to native space to use as my WM label file in DIPY, it won't let me reslice it. I used Freesurfer's recomendations to go from Tal back to native space and to convert to NIFTI. When I load that file into DIPY and use DIPY's reslice command (from 0.8mm to 2mm), it looks like it's processing, but the shape and zooms haven't changed. Is there something special that I have to do with data converted from Freesurfer? Thanks!
zlf123-marker
@zlf123-marker
@jhlegarreta Thank you very much
Shreyas Fadnavis
@ShreyasFadnavis
Hi @wadeatUT, Can you share the piece of code where you are reslicing the data.. I don't think anything special is needed.. it would be helpful to see what the exact issue is!
Wade Weber
@wadeatUT
@ShreyasFadnavis Thanks for the help. Here's a screenshot of my Jupyter Notebook:
reslice.png
@ShreyasFadnavis Also, if I change "new_zooms" to (3., 3., 3.), the new vox size comes out as (1.2, 1.2, 1.2) instead of 3mm iso? Weird.
Shreyas Fadnavis
@ShreyasFadnavis
@wadeatUT I guess this line is missing: zooms = orig_t1FS_img.header.get_zooms()[:3] and should solve the issue!
Shreyas Fadnavis
@ShreyasFadnavis
dipy_workshop.png
Hi everyone! The DIPY Workshop 2020 Early Bird registration ends in 3 days! To all those interested, visit: https://workshop.dipy.org/ | for any questions mail us: dipy@python.org and for updates, follow us on Twitter: @dipymri !! Looking forward to see you in March at Indiana University Bloomington!
Eleftherios Garyfallidis
@Garyfallidis
:star2: :star: :rocket: :rocket: :rocket: :rocket: :rocket: :rocket: :star: :star2:
Eleftherios Garyfallidis
@Garyfallidis
Thank you for sharing @ShreyasFadnavis
Serge Koudoro
@skoudoro
:+1: :+1: :+1:
zlf123-marker
@zlf123-marker
image.png
@jhlegarreta Traceback (most recent call last):
File "C:/Users/20190401-1/Desktop/image-registration/viz_t1_dti.py", line 55, in <module>
tractogram = load_tractogram('AF_left.trk',example_filename)
File "D:\Anaconda3\lib\site-packages\dipy\io\streamline.py", line 129, in load_tractogram
if not is_header_compatible(filename, reference):
File "D:\Anaconda3\lib\site-packages\dipy\io\utils.py", line 233, in is_header_compatible
reference_2)
File "D:\Anaconda3\lib\site-packages\dipy\io\utils.py", line 195, in get_reference_info
voxel_order = ''.join(nib.aff2axcodes(affine))
TypeError: sequence item 0: expected str instance, NoneType found
@jhlegarreta It's running like an error
Francois Rheault
@frheault
@zlf123-marker is your affine all zeros?
Scott Trinkle
@scott-trinkle
quick visualization question. is there a way to get different segments of individual streamlines to have different colors? right now the default "directional" colormap seems to give the whole streamline the same color based on the overall direction, but I would like to see how the direction changes across long streamlines
zlf123-marker
@zlf123-marker
image.png
image.png
@frheault
not all zeroes
zlf123-marker
@zlf123-marker
Hello everyone, I want to know how to display TRK file and t1 file at the same time?
Everybody, I want to know how to display my TRK file and t1 file at the same time?
Eleftherios Garyfallidis
@Garyfallidis
@scott-trinkle yes you can. You can set the colors for every point of the streamline.
@zlf123-marker are your t1 and TRK in the same coordinate system? If yes then see https://dipy.org/documentation/1.0.0./examples_built/viz_advanced/#example-viz-advanced and dipy_horizon. If not then you have to do registration first.
zlf123-marker
@zlf123-marker
Ok, thank you very much
paras sethi
@dreadysharp_twitter
Hey guys.. I wanted to read about how blockwise 3d denoising has been implemented ..
But I see 2 methods nlmeans and non_local_means different
Can you help which one is the actual implementation
Serge Koudoro
@skoudoro
Hi @dreadysharp_twitter, the actual one is nlmeans. We plan to remove the non_local_means and improve memory management for nlmeans.
paras sethi
@dreadysharp_twitter
I need to know the one which uses block averaging
I am thinking about working on it for a college project . I feel like the nlmeans doesn't have block averaging abilities as mentioned in the documentation. Can you mention the implementation which uses block averaging for 3d denoising
Boojum
@intein
Hello. No eddy current correction at the moment?
Eleftherios Garyfallidis
@Garyfallidis
@Boojum you can do eddy current correction by affinely registering all non b0 volumes to the b0. We are planning to have a tutorial and interface about that but in the meantime I can send you an example script.
Boojum
@intein
@Garyfallidis Sounds great. If it's not a problem for you, please send me the script
maruti07
@maruti07
Hi, i donot have a reverse phase encoded b0. I do have a spin echo feild maps but its echospacing and EPI factor do not match the diffusion parameters. I wanted to know if i can use them for distortion correction?
paras sethi
@dreadysharp_twitter
I'm sorry to ask again.. but can someone kindly tell which of the nlmeans implementation works with block averaging . It will be very helpful
Eleftherios Garyfallidis
@Garyfallidis
it does
Serge Koudoro
@skoudoro
both
zlf123-marker
@zlf123-marker
Excuse me, what should I do if I want to display t1 in volume data instead of slices?
Eleftherios Garyfallidis
@Garyfallidis
Not sure what you mean @zlf123-marker can you elaborate?
Eleftherios Garyfallidis
@Garyfallidis
@intein can you please send me an email so I can forward to you the script?
zlf123-marker
@zlf123-marker
@Garyfallidis That is to say I want to display t1 data as a whole rather than as a slice of x,y, and z, what should I do?
zlf123-marker
@zlf123-marker
May I ask how the dipy library displays NRRD data?
paras sethi
@dreadysharp_twitter
I would also need to know how can we display the data in volume instead of just some slices
My email Id is psethi432@gmail.com
anm-ranjan
@anm-ranjan
@Garyfallidis Could you provide me with the eddy current correction script as well? My email ID is anm.ranjan@gmail.com. It will be of great help!