Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

FLAIR image with corresponding bval/bvec files #449

Closed
dlevitas opened this issue Nov 17, 2020 · 4 comments
Closed

FLAIR image with corresponding bval/bvec files #449

dlevitas opened this issue Nov 17, 2020 · 4 comments

Comments

@dlevitas
Copy link

I believe this in relation to #393, but wanted to clarify. I've got what appears to be a FLAIR acquisition (though only 2D) with bval and bvec files (all values are zero). My thought is that this is some form of localizer and should be disregarded, but want to double check before proceeding. I've shared the specific dicoms through Google Drive.

Thanks,

Dan

@neurolabusc
Copy link
Collaborator

Your sample DICOMs include the private tag 0043,1039) which suggests it does not have a b-weighted. dcm2niix infers that this might be a b=0 series used for spatial undistortion with TOPUP. In your case, this is a false positive, as this is not a diffusion sequence despite the presence of 0043,1039. I would suggest you just ignore the bvec/bval file created. For the developmental branch I have changed the heuristic for GE data so the image must both have a b-value associated with it as well as being an EPI sequence (0018,9018). This will avoid your exemplar of a false alarm, though risks missing true diffusion sequences where 0018,9018 is not provided.

@mr-jaemin do you have a preferred method for determining if a GE DICOM series is from a diffusion sequence (including series where all volumes are b=0)?

@mr-jaemin
Copy link
Collaborator

You could use the private tag "Vas collapse flag" (0043, 1030) SS.

For DTI, the values should be 14 (T2) or 16 (DTI including b0 with gradient (0, 0, 0). For DWI, the values should be 14 (T2) or one of the following values depending on diffusion gradient direction applied:

// Diffusion image type values for the diffusion direction annotation field (0x0043,0x1030)
DiffusionRightLeftDicomValue = 3;
DiffusionAnteriorPosteriorDicomValue = 4;
DiffusionSuperiorInferiorDicomValue = 5;
DiffusionObliqueRightLeftDicomValue = 9;
DiffusionObliqueAnteriorPosteriorDicomValue = 10;
DiffusionObliqueSuperiorInferiorDicomValue = 11;
DiffusionT2DicomValue = 14;
DiffusionCombinedDicomValue = 15;
DiffusionDtiDicomValue = 16;
DiffusionDirection1DicomValue = 43;
DiffusionDirection2DicomValue = 44;
DiffusionDirection3DicomValue = 45;
DiffusionDirection4DicomValue = 46;

DTI_ADC = 20
DTI_FA = 18
DTI_T2W = 19

Alternately, For EPI, you can re-use this logic in the sliceTimingGE:
(This criteria doesn't apply to PROFELLER Diffusion (DWI, not DTI)
// Diffusion
else if ( (d->epiVersionGE == 2) || (d->internalepiVersionGE == 2) || (strstr(ioptGE,"DIFF") != NULL) )

@dlevitas
Copy link
Author

Great, thanks. Should I also disregard the FLAIR acquisition, since it appears to be a 2D slice?

@dlevitas
Copy link
Author

I recant my earlier statement -- I realize now that the FLAIR is anisotropic. All good here, thanks!

yarikoptic added a commit to neurodebian/dcm2niix that referenced this issue Apr 6, 2021
* commit 'v1.0.20200427-97-g0587941': (22 commits)
  Overflow for Siemens data [missing protocol name] (rordenlab#466)
  MacOS notarization, minor tweaks
  Increase details for series stacking, enhance merge override mode (rordenlab#463)
  Bump version date
  Retain Philips Scaling Values for images where scaling does not vary but volumes must be separated (rordenlab#461)
  Support huge PAR/REC files (rordenlab#460)
  Prevent dti4D overflow
  Fix PAR/REC ADC detection (rordenlab#462)
  DICOM field map calibrated in Hz given name _fieldmaphz (rordenlab#455)
  PAR/REC field map calibrated in Hz given name _fieldmaphz (rordenlab#455)
  Ignore LocationsInAcquisition (0020,1002) if number of slices is not evenly divisible with this value (rordenlab#459)
  Convert DICOM series where bit allocated (0028,0100) varies across slices (rordenlab#458)
  Clear RepetitionTimeExcitation (rordenlab#457)
  leak (rordenlab#454)
  Single file mode memory allocation (rordenlab#454)
  When -n is specified, only save BIDS for requested series (rordenlab#453)
  Use 1st study time if multiple times provided.
  Apple. M1. use C++ not. C
  Fix CMake for Apple Silicon (note similar change required for CloudFlare zlib)
  Only report b-value for GE data with 0043,1039 if EPI (0018,9018) (rordenlab#449)
  ...
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants