# set values for output image
# image size
imageExtent = (int(extX), int(extY), int(extZ), 1, 1, 1)
outImage.setImageExtent(imageExtent)
# intensity range
outImage.setMinVoxelValue(float(wCenter)-float(wWidth))
outImage.setMaxVoxelValue(float(wCenter)+float(wWidth))
# spacing
outImage.setVoxelSize(int(voxelSpacingXY[0]), int(voxelSpacingXY[1]), int(voxelSpacingZ))
# orientation and spacing
oX = [float(imageOrientationPatient[0]),float(imageOrientationPatient[1]),float(imageOrientationPatient[2])]
oY = [float(imageOrientationPatient[3]),float(imageOrientationPatient[4]),float(imageOrientationPatient[5])]
oZ = cross(oX, oY)
s = matrix([[float(voxelSpacingXY[0]), 0, 0, 0], [0, float(voxelSpacingXY[1]), 0, 0], [0, 0, float(voxelSpacingZ), 0], [0, 0, 0, 1]])
r = matrix([[oX[0],oY[0],oZ[0],0], [oX[1],oY[1],oZ[1], 0], [oX[2],oY[2],oZ[2],0], [0,0,0,1]])
t = matrix([[1., 0, 0, float(imagePositionPatient[0])], [0, 1, 0, float(imagePositionPatient[1])], [0, 0, 1, float(imagePositionPatient[2])], [0, 0, 0, 1]])

wm = t * r * s

outImage.setVoxelToWorldMatrix(wm)