DICOM Image Processing

I have a DICOM image that I am trying to calculate and plot the pre-sampled Modulation Transfer Function of. I have the paper that shows me what method to use but I am not sure how to code it.

Please contact me for more details.