Method for image intensity correction using extrapolation and adaptive smoothing

A method of image intensity correction, comprising:

• a) receiving original image data for an original image, wherein the image data for the original image comprises a plurality of magnitude values corresponding to a plurality of pixels in the original image;
• b) creating an intensity correction template, wherein the intensity correction template comprises a plurality of correction values corresponding to at least a portion of the plurality of pixels, wherein creating the intensity correction template comprises:
• i) determining hole pixels and image support pixels based on the plurality of magnitude values corresponding to the plurality of pixels in the original image, wherein each of the pixels of the plurality of pixels is determined to be a hole pixel or an image support pixel,
• ii) setting the correction values of the image support pixels in the intensity correction template to the magnitude values of the corresponding pixels of the original image;
• and
• iii) setting the correction values for at least a portion of the hole pixels in the intensity correction template near the image support pixels, wherein each of the correction values for the hole pixels in the intensity correction template near the image support pixels are related to the correction values of the image support pixels with respect to which the hole pixel is near;
• c) generating an intensity correction map, wherein generating an intensity correction map comprises altering the correction values of the image support pixels of the intensity correction template based on the values of the image support pixels of the intensity correction template and the values of the hole pixels of the intensity correction template near the image support pixel;
• and
• d) producing corrected image data for a corrected image, wherein the corrected image data for the corrected image comprises a plurality of corrected magnitude values corresponding to the image support pixels, wherein the corrected magnitude value of each image support pixel is proportional to the magnitude value of the image support pixel in the original image divided by the correction value of the image support pixel in the intensity correction map.