Mapping/Rescaling intensities of 3d images

I am working with 3D CT images with HU units of 0-4017. I would like to map the values to 800-1200 HU. I Have been searching a lot to find a solution for mapping 3D matrices/images, but no success unfortunately!
Any help is very much appreciated…
Thank you