signal intensity value in MRI dicom file

11 views (last 30 days)
Can Anybody help me to extract matlab code to extract signal intensity value in MRI dicom file with the corresponding pixel coordination? here i have attached dicom file Brain http://www.osirix-viewer.com/datasets/DATA/BRAINIX.zip

Accepted Answer

harjeet singh
harjeet singh on 28 Dec 2015
hello Ranga use this to exract signal intensity
X = uint8(dicomread('IM-0001-0001'));
figure(1)
imshow(X)
pixel_co=[256 200];
pixel_vale=X(pixel_co(1),pixel_co(2))
  6 Comments
Walter Roberson
Walter Roberson on 29 Dec 2015
If you get 0 every time then your data might have been negative and you should probably not have converted to uint8... as I cautioned could be a problem. Remove the uint8() part and try again and see what the values are.
Image Analyst
Image Analyst on 29 Dec 2015
If you remove uint8() and find that it is negative, you will probably have to use [] in imshow():
imshow(X, []);

Sign in to comment.

More Answers (3)

Image Analyst
Image Analyst on 29 Dec 2015
Try calling impixelinfo() to interactively view the value as you mouse around over the image.
grayImage = dicomread(fullFileName);
imshow(grayImage);
colorbar;
hp = impixelinfo();
Or you can call imtool().

Ranga Tudugala
Ranga Tudugala on 30 Dec 2015
Okay Thanks...Guys..
I have extracted pixel value which I needed but now im struggling. I need to select every pixel which having particular value (range) in the image. any body ?
  7 Comments
Ranga Tudugala
Ranga Tudugala on 31 Dec 2015
Articulacy It CAN.. that the researchers Do... what they seems cnt.. I think you are genius than me. you can do it easily.. I think i have found the best resource person.. for that..
here I hv attached the paper "MRI-based treatment planning with pseudo CT generated through atlas registration." http://www.ncbi.nlm.nih.gov/pubmed/24784377
Image Analyst
Image Analyst on 31 Dec 2015
Well Star could answer better than me, but I'm reading "A pseudo CT...[snip]...was generated by deforming atlas CT images...[snip]. The required deformation for atlas CT images was derived from a nonlinear registration of conjugated atlas MR images"
So it sounds like they're aligning/warping CT images with the formula for warping based on MRI images, or what I called multi-modality registration. Maybe that's what you meant. I was thinking that you wanted to take MRI images and convert them to have it look just as if you'd taken a CT of the patient instead of an MRI. I think you can align CT or MRI volumes in a variety of ways, but I don't think you can get one from the other. Otherwise hospitals would not have two separate multi-million dollar instruments when they could get by with one.
Star, what do you think?

Sign in to comment.


Ranga Tudugala
Ranga Tudugala on 31 Dec 2015
Edited: Walter Roberson on 31 Dec 2015
Im getting this error after this codes
movingRegisteredAffineWithIC = imregister(moving,fixed,'affine',optimizer,metric,...
'InitialTransformation',tformSimilarity);
This is the error
Error using imregtform>parseInputs (line 240)
Parameters must be a character.
Error in imregtform (line 119)
parsedInputs = parseInputs(varargin{:});
Error in imregister (line 115)
tform = imregtform(varargin{:});
can somebody tell what happened ?

Community Treasure Hunt

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

Start Hunting!