Estoy tratando de crear una imagen de mew dicom a partir de una matriz numpy de tamaño estándar (512 x 512 o 256 x 256). Parece que esto debería ser sencillo, y he adaptado mi código dehttp://code.google.com/p/pydicom/source/browse/source/dicom/examples/write_new.py, que parece ejecutar el mismo proceso, pero cuando guardo el archivo, no puedo verlo ni en RadiAnt ni en MicroDicom.

import dicom, dicom.UID
from dicom.dataset import Dataset, FileDataset

def write_dicom(pixel_array,filename):

    file_meta = Dataset()
    ds = FileDataset(filename, {},file_meta = file_meta,preamble="\0"*128)
    ds.PixelData = pixel_array.tostring()
    ds.save_as(filename)
    return

if __name__ == "__main__":
    import numpy as np
    pixel_array = np.tile(np.arange(256).reshape(16,16), (16,16)) * 4
    write_dicom(pixel_array,'pretty.dcm')