Methods, apparatuses, and computer program products for determining a modulation transfer function of an imaging system
View Patent ↗Methods, apparatuses, and computer program products are provided for determining a modulation transfer function of an imaging system. A method may include accessing an image of a phantom having a substantially circular shaped feature captured by the imaging system. The method may further include detecting the circular shaped feature within the image. The method may additionally include defining at least one line extending from a point within the detected circular shaped feature to a point outside of the circular shaped feature. The method may also include determining an edge spread function based at least in part on the defined at least one line. The method may further include determining the modulation transfer function of the imaging system based at least in part on the determined edge spread function. Corresponding apparatuses and computer program products are also provided.
1. A method for determining a modulation transfer function of an imaging system, the method comprising:
accessing an image of a phantom having a circular shaped feature captured by the imaging system;
detecting the circular shaped feature within the image;
defining at least one line extending from a point within the detected circular shaped feature to a point outside of the circular shaped feature;
determining, by a processor, an edge spread function based at least in part on the defined at least one line by at least:
determining an edge spread function for each defined line; and
determining an average edge spread function for the determined edge spread functions; and
determining the modulation transfer function of the imaging system based at least in part on the determined average edge spread function.
2. The method of claim 1 , wherein defining the at least one line comprises defining at least one radial line extending from a center point of the circular shaped feature and extending outward for a length greater than a radius of the circular shaped feature.
3. The method of claim 1 , wherein determining the edge spread function comprises determining the edge spread function based at least in part on pixel values along the at least one line.
4. The method of claim 1 , wherein in an instance in which a plurality of lines are defined and at least one line does not originate from a center point of the circular shaped feature, determining the edge spread function comprises:
determining a center point for each determined edge spread function based at least in part on an average value of pixels inside the circular shaped feature and on an average value of pixels in a region of the image outside of the circular shaped feature;
using the determined center points to align the edge spread functions; and
determining an average edge spread function based on the aligned edge spread functions.
5. The method of claim 1 , wherein determining the modulation transfer function of the imaging system based at least in part on the determined edge spread function comprises:
deriving a line spread function by calculating a derivative of the edge spread function; and
determining the modulation transfer function by determining the modulus of the discrete Fourier transform of the line spread function.
6. The method of claim 1 , wherein detecting the circular shaped feature within the image comprises using a Hough transform to detect the circular shaped feature.
7. The method of claim 1 , wherein the phantom comprises an American College of Radiology Computed Tomography Accreditation phantom.
8. The method of claim 1 , wherein the imaging system comprises a medical imaging system.
9. An apparatus for determining a modulation transfer function of an imaging system, the apparatus comprising at least one processor, wherein the at least one processor is configured to cause the apparatus to at least:
access an image of a phantom having a circular shaped feature captured by the imaging system;
detect the circular shaped feature within the image;
define at least one line extending from a point within the detected circular shaped feature to a point outside of the circular shaped feature;
determine an edge spread function based at least in part on the defined at least one line;
determine an edge spread function for each defined line;
determine an average edge spread function for the determined edge spread functions; and
determine the modulation transfer function of the imaging system based at least in part on the determined average edge spread function.
10. The apparatus of claim 9 , wherein the at least one processor is further configured to cause the apparatus to define the at least one line at least in part by defining at least one radial line extending from a center point of the circular shaped feature and extending outward for a length greater than a radius of the circular shaped feature.
11. The apparatus of claim 9 , wherein the at least one processor is further configured to cause the apparatus to determine the edge spread function based at least in part on pixel values along the at least one line.
12. The apparatus of claim 9 , wherein the at least one processor is further configured to cause the apparatus to:
determine a center point for each determined edge spread function based at least in part on an average value of pixels inside the circular shaped feature and on an average value of pixels in a region of the image outside of the circular shaped feature;
use the determined center points to align the edge spread functions; and
determine an average edge spread function based on the aligned edge spread functions.
13. The apparatus of claim 9 , wherein the at least one processor is further configured to cause the apparatus to determine the modulation transfer function of the imaging system at least in part by:
deriving a line spread function by calculating a derivative of the edge spread function; and
determining the modulation transfer function by determining the modulus of the discrete Fourier transform of the line spread function.
14. The apparatus of claim 9 , wherein the at least one processor is further configured to cause the apparatus to use a Hough transform to detect the circular shaped feature.
15. The apparatus of claim 9 , wherein the phantom comprises an American College of Radiology Computed Tomography Accreditation phantom.
16. The apparatus of claim 9 , wherein the imaging system comprises a medical imaging system.
17. A computer program product for determining a modulation transfer function of an imaging system, the computer program product comprising at least one non-transitory computer-readable storage medium having computer-readable program instructions stored therein, the computer-readable program instructions comprising:
program instructions configured to access an image of a phantom having a circular shaped feature captured by the imaging system;
program instructions configured to detect the circular shaped feature within the image;
program instructions configured to define at least one line extending from a point within the detected circular shaped feature to a point outside of the circular shaped feature;
program instructions configured to determine an edge spread function based at least in part on the defined at least one line;
program instructions configured to determine an edge spread function for each defined line;
program instructions configured to determine an average edge spread function for the determined edge spread functions; and
program instructions configured to determine the modulation transfer function of the imaging system based at least in part on the determined average edge spread function.
18. The computer program product of claim 17 , wherein the program instructions configured to determine the edge spread function comprise:
program instructions configured to determine an average edge spread function for the determined edge spread functions
wherein the program instructions configured to determine the modulation transfer function comprise program instructions configured to determine the modulation transfer function based at least in part on the average edge spread function.