Load, open and display 2dseq file (MRI, bruker)

조회 수: 33 (최근 30일)
viet le
viet le 2017년 6월 27일
댓글: Maqbool 2024년 1월 16일
Hello,
I have an MR image data which acquired at Bruker 4.7 T MRI scanner. (Bruker file attached). This code by Matlab can not run and I did not know how to fix? really appreciate your help.
Thank you.
num_dif = 1;
dimension = [128 128 100]; % [x y z] the dimension of the 2dseq file, suppose average=1
voxel_size = [0.1 0.1 0.1]; % [x y z] pixel spacing in mm
f1 = fopen('2dseq','r');
A=fread(f1,'int16');
A=reshape(A,dimension(1),dimension(2),dimension(3),num_dif);
% slope is a dimension(3) x num_dif matrix that defines the scaling parameters
% input the slope here
slope=[81.3139111191806 81.3139111191806 81.3139111191806 81.3139111191806
...
]
slope = slope./min(min(slope));
slope = reshape(slope',dimension(3),num_dif);
for i = 1:dimension(3)
for j = 1:num_dif
A(:,:,i,j) = A(:,:,i,j)/slope(i,j);
end
end
% image0 ~ imageN, where N is the number of diffusion gradient encoding
image0 = reshape(A(:,:,:,1),prod(dimension),1);
  댓글 수: 1
Oscar Jalnefjord
Oscar Jalnefjord 2017년 6월 28일
There is only data from a single slice in the attached file 2dseq. if you set
dimension = [128 128 100]
and set slope to a scalar value your code runs

댓글을 달려면 로그인하십시오.

답변 (2개)

Khieu
Khieu 2018년 4월 16일
편집: Khieu 2018년 4월 16일

1. You should take a look at parameters files: acqp, method. (text file) these files contain all image parameters that suitable for you. 2. I found that your image is 2D based on: - method file: ##$PVM_Matrix=( 2 ) 128 128 So your image data size should be: 128x128 (complex or real you should check on method file or reco file)

A=fread(f1,'int16');
% if complex data, do three following lines 
RealI = A(1:2:end);
ImgI = A(2:2:end);
cmplxImage = complex(RealI, ImgI);

Then you just reshape that variable to correct image dimension...

Good luck

ps: your image is T1-weighted image, not diffusion image, so you do not need a variable like num_dif (number of diffusions)

  댓글 수: 1
Maqbool
Maqbool 2024년 1월 16일
This only for one real and imaginary lines, what about I needs all of them and to do FFT

댓글을 달려면 로그인하십시오.


Cecil Yen
Cecil Yen 2018년 10월 29일
I wrote a Matlab function to read Bruker's 2dseq file for my own purposes years ago. It will phrase the header file and load the binary image for you. I use some tricks to speed up the reading, which shall be the fastest Matlab implementation this far. I uploaded it to Matlab Central such that everyone can use it.
https://www.mathworks.com/matlabcentral/fileexchange/69177-read_2dseq-quickly-reads-bruker-s-2dseq-mri-images
  댓글 수: 2
Avigdor
Avigdor 2018년 10월 29일
Hi,
There seems to be an issue with the find_newline function and the read_2dseq does not run.
Amaan Khan
Amaan Khan 2023년 7월 25일
Hi, is there any documentation for this by any chance?

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Convert Image Type에 대해 자세히 알아보기

태그

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!