Scaling of intensity value of captured jpg image

I have a captured CT image (.jpg) of segittal plane. I want to scale it (.jpg) corresponding to .dcm image (512x512 uint16). The intensity range of .dcm image is from [0 to 3070]. Or is there any other way to convert it into dicom format with same intensity range?

Thanks for help...

Answer 1

[KALYAN ACHARJYA]

Is this?

input_image=imread('filename.jpg');
dicomwrite(input_image,'dicom_img_name.dcm');

Or

You want to map the RGB pixel values [0-255] ranges to [0 to 3070].

Good Luck!

[kritika joshi]

Yes I want RGB pixel values [0-255] ranges to [0 to 3070].

[KALYAN ACHARJYA]

I have tested with sample jpg image, it has succesfully save the image in dcm format, but still the pixel values in the range 0- 255.

im_data=imread('1-1.jpg');
dicomwrite(im_data,'dicom_img_name.dcm');
%Check the save dcm image
dicom_data=dicomread('dicom_img_name.dcm')

Mapping, Try

out_put_data=double(gray_image).*(3070/255);

after that you can dicomwrite the out_put_data as dicom image??

[kritika joshi]

Thanks for your efforts. I had tried this and getting a round off range of (0 to 2966).

[KALYAN ACHARJYA]

It's my pleasure, happy to know that, the problem has been solved

Answer 2

[darova]

I'm not familiar with DICOM images but looks simple:

I1 = imread('3.jpg');
I2 = rgb2gray(I1);
I3 = I2-min(I2);                    % make min value 0 
I3 = I3*3070/max(I3(:));            % make max value 3070
dicomwrite(I2,'test_scaled_3070.dcm')
imshow(I2)

[kritika joshi]

Thanks for your answer. I have tried this, but I3 = I3*3070/max(I3(:)); is giving a binary matrix.