comp.lang.idl-pvwave archive: archive » Re: Transferring image pixel values to a DICOM file

Home » Public Forums » archive » Re: Transferring image pixel values to a DICOM file

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: Transferring image pixel values to a DICOM file [message #76978 is a reply to message #76976]

Wed, 20 July 2011 03:37

dplatten
Messages: 32
Registered: December 2007

Member

Hi Vivek,

I used to use another method to create computed tomography DICOM files containing my own image data. I took an existing DICOM image, opened it, and then wrote my new image data over the existing image data in the file, and then saved the file. It was a nasty way to do it, but it worked. The CT images I was working on always had 512x512 pixels. An example of the code is below. There must be a more elegant way of doing this!

Regards,

David

PRO modifyDICOM

; create an array to hold the generated image
x_dim = 512
y_dim = 512
new_data = INTARR(x_dim,y_dim)
new_data[*,*] = 0

; generate some test data to write to the image
; a diagonal cross
FOR x = 0, (x_dim-1) DO BEGIN
new_data[x,x] = 1000
new_data[(x_dim-1)-x,x] = 1000
ENDFOR

; Or you could load in a tif file with your image data:
;image_array = read_tiff('c:\TEMP\StraightEdgeReallyBlurred.tif', channels=1)
;new_data = fix(image_array)

; Open an existing DICOM file to modify. This must have the
; same pixel dimensions as the one you want to create.
dicomFile = DIALOG_PICKFILE(/MUST_EXIST)

; open the file for writing
OPENU, dcmfile, dicomFile, /GET_LUN, /Append
POINT_LUN, -dcmfile, a

;position file pointer 512*512 from end of file
POINT_LUN, dcmfile, a-(512l*512*2)

; replace the pixel data with the generated data
WRITEU, dcmfile, new_data

; close the file and free the lun
CLOSE, dcmfile
FREE_LUN, dcmfile
END

Report message to a moderator

[Message index]

		Re: Transferring image pixel values to a DICOM file By: Vivek on Wed, 20 July 2011 07:52
		Re: Transferring image pixel values to a DICOM file By: dplatten on Wed, 20 July 2011 03:37
		Re: Transferring image pixel values to a DICOM file By: Vivek on Tue, 19 July 2011 07:18
		Re: Transferring image pixel values to a DICOM file By: Vivek on Tue, 19 July 2011 07:06
		Re: Transferring image pixel values to a DICOM file By: dplatten on Tue, 19 July 2011 01:11

Previous Topic:	multiple plots of different size on same page
Next Topic:	Shell startup file not recognized in IDL 8.0 Workbench

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Oct 11 02:27:22 PDT 2025

Total time taken to generate the page: 0.96148 seconds