The synthesis of patient data of a certain medical image modality by applying an image processing pipeline starting from other modality is receiving a lot of interest recently, as it allows to save acquisition time and sometimes avoid radiation to the patient. An example of this is the creation of computerized tomography volumes from magnetic resonance imaging data, which can be useful for several applications such as electromagnetic simulations, cranial morphometry and attenuation correction in PET/MR systems. We present a fast patch-based algorithm for this purpose, implemented using graphics processing unit computing techniques and gaining up to x 15.9 of speedup against a multicore CPU solution and up to about x 75 against a single core CPU solution.