Measurement of Contrast and Spatial Resolution for the Photothermal Imaging Method

: Analyzing the quality of images generated from an imaging method is essential for determining the limits and applicability of that method. This study analyzed the quality of images resulting from a photothermal imaging method by applying the line spread function and the modulation transfer function to the spatial resolution and contrast, on the basis of certain parameters of the photothermal imaging method for a copper-resin double-layered structure. The parameters are the ratio of the ﬁrst-layer thickness to the thermal di ﬀ usion length ( L f / L d ) and the ratio of the pump-beam radius to the thermal di ﬀ usion length ( R b / L d ). The phase delay proﬁle (edge response function, ERF) of the subsurface structure derived from the photothermal imaging method becomes dimensionless upon division by the thermal di ﬀ usion length; as the ratio L f / L d increases, the spatial resolution and contrast increase.


Introduction
Recently, the photothermal imaging method has been studied as an industrial tomographic tool. It acquires morphological information from the inside of an object by means of the photothermal effect of a temporally modulated heat source [1][2][3][4]. The photothermal effect is a physical phenomenon that occurs when a material absorbs light energy and converts it to thermal energy [5,6]. As shown in Figure 1a, when a material is heated by a temporally modulated pump beam, the temperature of both the air layer on the material surface together and the material itself increases as heat is conducted from the material to the air. In addition, the refractive index of the air layer changes with temperature because of the density change. It is possible to achieve a modulated temperature change through the refraction of the probe beam that passes through the air layer. The physical phenomena that result from the photothermal effect are based on the thermal properties of the target material to be measured. Hence, the phase delay between the period of the refraction of the probe beam and that of the heat source (light energy) must be analyzed in order to acquire information on the thermal properties of the specimen [7][8][9][10][11]. This method can be applied to imaging. As illustrated in Figure 1b, when the surface of a material is scanned by the pump beam and probe beam, the phase delay changes as the thermal properties are changed by the distribution of internal materials. This photothermal imaging method has been studied for more than ten years by McDonald, Wetsel, and Friedrich [12][13][14], since Murphy and Aamodt published a study on how to acquire morphological information on a subsurface internal structure from the phase delay [4]. The related studies to the photothermal imaging method have been focused on the possibility of detecting subsurface structures using photothermal signals (phase delay). Previous studies have shown that the photothermal imaging method is experimentally feasible. However, there have been attempts to study the quality of images acquired from the photothermal imaging method. Therefore, the limitation of measurement in photothermal imaging has not been investigated or reported.
The applicability of the developed imaging method depends on the image quality for the image of the internal structure acquired from the imaging method. Image quality is represented by spatial resolution, contrast, sharpness, and other metrics, and the scope of imaging method applicability depends on these factors [15][16][17].
Image quality varies with the type of imaging method, specifically owing to the configuration of the components of the relevant equipment, the characteristics of the object to be measured, and the measuring system operator′s skills. Analyzing how such factors affect image quality not only requires comprehensive and complicated tasks, but is also of great importance in evaluating a given imaging method [13][14][15].
The image quality analysis techniques and processes are well established. However, the quality of images acquired from the photothermal imaging method has not been quantitatively analyzed. Moreover, no one has investigated whether the image quality analysis techniques can be applied to photothermal imaging. The present study was carried out to analyze image quality by focusing on the contrast and spatial resolution of an imaging method that utilizes the photothermal effect. In addition, the factors that influence the acquired images from the photothermal imaging method are investigated, expecting to provide experimental measurement conditions for the photothermal imaging method.

Contrast
Contrast is a measure of the difference between the magnitude of signals from an internal structure and those from the surrounding or background of the internal structure. In general, the purpose of an imaging method is to examine the inside of an object to be measured and to detect any This photothermal imaging method has been studied for more than ten years by McDonald, Wetsel, and Friedrich [12][13][14], since Murphy and Aamodt published a study on how to acquire morphological information on a subsurface internal structure from the phase delay [4]. The related studies to the photothermal imaging method have been focused on the possibility of detecting subsurface structures using photothermal signals (phase delay). Previous studies have shown that the photothermal imaging method is experimentally feasible. However, there have been attempts to study the quality of images acquired from the photothermal imaging method. Therefore, the limitation of measurement in photothermal imaging has not been investigated or reported.
The applicability of the developed imaging method depends on the image quality for the image of the internal structure acquired from the imaging method. Image quality is represented by spatial resolution, contrast, sharpness, and other metrics, and the scope of imaging method applicability depends on these factors [15][16][17].
Image quality varies with the type of imaging method, specifically owing to the configuration of the components of the relevant equipment, the characteristics of the object to be measured, and the measuring system operator s skills. Analyzing how such factors affect image quality not only requires comprehensive and complicated tasks, but is also of great importance in evaluating a given imaging method [13][14][15].
The image quality analysis techniques and processes are well established. However, the quality of images acquired from the photothermal imaging method has not been quantitatively analyzed. Moreover, no one has investigated whether the image quality analysis techniques can be applied to photothermal imaging. The present study was carried out to analyze image quality by focusing on the contrast and spatial resolution of an imaging method that utilizes the photothermal effect. In addition, the factors that influence the acquired images from the photothermal imaging method are investigated, expecting to provide experimental measurement conditions for the photothermal imaging method.

Contrast
Contrast is a measure of the difference between the magnitude of signals from an internal structure and those from the surrounding or background of the internal structure. In general, the purpose of an imaging method is to examine the inside of an object to be measured and to detect any structural Appl. Sci. 2019, 9, 1996 3 of 15 differences from the intended design or any abnormal defects. In order to easily detect such defects, imaging methods with high contrast are preferred [18].
This section focuses on the modulation index that quantifies contrast in reference to periodic signals of constant amplitude, as shown in Figure 2.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 15 structural differences from the intended design or any abnormal defects. In order to easily detect such defects, imaging methods with high contrast are preferred [18]. This section focuses on the modulation index that quantifies contrast in reference to periodic signals of constant amplitude, as shown in Figure 2. The modulation index ( ) of a periodic signal ( ) whose maximum and minimum amplitudes of signal are and , respectively, can be defined as Equation (1) [18,19], where it is the ratio of the average value of this signal (( + )/2) to the size of the amplitude (( − )/2) of the periodic signal ( ).
In the case of an imaging method that utilizes the photothermal effect, as with other imaging methods, the contrast changes depending on the size of the object to be measured (in this study, the size of the internal structure). The size of an object to be measured can be quantified on the basis of the spatial frequency, which can be derived by applying the Fourier transform to a spatial structure in the space. In summary, an imaging method′s contrast can be quantified from the modulation index, and the size of a structure to be measured can be quantified from the spatial frequency [15,17]. The relation between the spatial frequency and modulation index is called the modulation transfer function (MTF). The MTF of an imaging method can be calculated by either of the following methods: A For line pairs of varied sizes, the modulation index of a spatial frequency corresponding to each line pair is calculated [20,21]. B A Fourier transform is applied to the line spread function (LSF) derived from line pairs of the same size [20,21].
The first method (Method A), which uses line pairs of varied sizes, is the most typical way to derive the MTF, but it is disadvantageous in that a large number of specimens is required in order to derive the entire MTF profile. For this reason, comparing the MTF in various conditions drastically increases the number of numerical analyses or experiments.
The second method (Method B), which utilizes the LSF, can be implemented in a relatively simple manner [21]. On the assumption that there are infinite specimens with the internal structure shown in Figure 3a in the x-and y-directions, 1-D scanning on the specimen surface using the photothermal imaging method will produce the phase delay change results shown in Figure 3b. The sigmoidal phase difference signal derived near the boundary shown in this figure is called the edge response function (ERF). As the ERF is differentiated, the LSF as in Figure 3c can be obtained.
(2) The modulation index (m f req ) of a periodic signal f (x) whose maximum and minimum amplitudes of signal are A max and A min , respectively, can be defined as Equation (1) [18,19], where it is the ratio of the average value of this signal ((A max + A min )/2) to the size of the amplitude ((A max − A min )/2) of the periodic signal f (x).
In the case of an imaging method that utilizes the photothermal effect, as with other imaging methods, the contrast changes depending on the size of the object to be measured (in this study, the size of the internal structure). The size of an object to be measured can be quantified on the basis of the spatial frequency, which can be derived by applying the Fourier transform to a spatial structure in the space. In summary, an imaging method s contrast can be quantified from the modulation index, and the size of a structure to be measured can be quantified from the spatial frequency [15,17]. The relation between the spatial frequency and modulation index is called the modulation transfer function (MTF). The MTF of an imaging method can be calculated by either of the following methods: A For line pairs of varied sizes, the modulation index of a spatial frequency corresponding to each line pair is calculated [20,21]. B A Fourier transform is applied to the line spread function (LSF) derived from line pairs of the same size [20,21].
The first method (Method A), which uses line pairs of varied sizes, is the most typical way to derive the MTF, but it is disadvantageous in that a large number of specimens is required in order to derive the entire MTF profile. For this reason, comparing the MTF in various conditions drastically increases the number of numerical analyses or experiments.
The second method (Method B), which utilizes the LSF, can be implemented in a relatively simple manner [21]. On the assumption that there are infinite specimens with the internal structure shown in Figure 3a in the xand y-directions, 1-D scanning on the specimen surface using the photothermal imaging method will produce the phase delay change results shown in Figure 3b. The sigmoidal phase difference signal derived near the boundary shown in this figure is called the edge response function (ERF). As the ERF is differentiated, the LSF as in Figure 3c can be obtained. Physically and mathematically, the LSF and MTF are in the relation of a Fourier transform as in Equation (2) [18]. By normalizing the results of applying the Fourier transform to the LSF, the MTF of the imaging system can be derived as shown in Figure 3d. Physically and mathematically, the LSF and MTF are in the relation of a Fourier transform as in Equation (2) [18]. By normalizing the results of applying the Fourier transform to the LSF, the MTF of the imaging system can be derived as shown in Figure 3d. These two methods of deriving the MTF (Methods A and B) have been found in a number of studies to have the same mathematical and physical meanings [15,20]. The present study derived the MTF by using the relatively simpler Method B and verified contrast changes depending on the spatial frequency. These two methods of deriving the MTF (Methods A and B) have been found in a number of studies to have the same mathematical and physical meanings [15,20]. The present study derived the MTF by using the relatively simpler Method B and verified contrast changes depending on the spatial frequency.

Spatial Resolution
The spatial resolution also presents a quantified value of image quality. The spatial resolution can be viewed as the ability of a measuring system to accurately distinguish two different objects at a certain distance from each other [15,17,18].
The determination of spatial resolution begins with the LSF. The LSF can be viewed as a response of the imaging system to line impulses (a linear boundary of the internal structure). Hence, the spatial resolution is determined from the form and size of the LSF [13,15,16]. Many studies have clarified that the spatial resolution can be quantified through the full width half maximum (FWHM) of the LSF of the imaging system [15,22].

Numerical and Experimental Analysis
This study numerically analyzed the effect of parameters of the measuring system (namely, the pump beam modulation frequency and pump beam radius) on the contrast and spatial resolution, specifically in relation to imaging a copper-resin double-layered structure that contains slits. This study also experimentally validated the representative result of numerical analysis.
The specimens of the copper-resin double-layered structure that were used in the numerical analysis and experiment are illustrated in Figure 4. The properties of the copper and resin are presented in Table 1. Additionally, Figure 5 presents the optical arrangement of the specimens and laser applied to the numerical analysis and experiment.

Spatial Resolution
The spatial resolution also presents a quantified value of image quality. The spatial resolution can be viewed as the ability of a measuring system to accurately distinguish two different objects at a certain distance from each other [15,17,18].
The determination of spatial resolution begins with the LSF. The LSF can be viewed as a response of the imaging system to line impulses (a linear boundary of the internal structure). Hence, the spatial resolution is determined from the form and size of the LSF [13,15,16]. Many studies have clarified that the spatial resolution can be quantified through the full width half maximum (FWHM) of the LSF of the imaging system [15,22].

Numerical and Experimental Analysis
This study numerically analyzed the effect of parameters of the measuring system (namely, the pump beam modulation frequency and pump beam radius) on the contrast and spatial resolution, specifically in relation to imaging a copper-resin double-layered structure that contains slits. This study also experimentally validated the representative result of numerical analysis.
The specimens of the copper-resin double-layered structure that were used in the numerical analysis and experiment are illustrated in Figure 4. The properties of the copper and resin are presented in Table 1. Additionally, Figure 5 presents the optical arrangement of the specimens and laser applied to the numerical analysis and experiment.

Spatial Resolution
The spatial resolution also presents a quantified value of image quality. The spatial resolution can be viewed as the ability of a measuring system to accurately distinguish two different objects at a certain distance from each other [15,17,18].
The determination of spatial resolution begins with the LSF. The LSF can be viewed as a response of the imaging system to line impulses (a linear boundary of the internal structure). Hence, the spatial resolution is determined from the form and size of the LSF [13,15,16]. Many studies have clarified that the spatial resolution can be quantified through the full width half maximum (FWHM) of the LSF of the imaging system [15,22].

Numerical and Experimental Analysis
This study numerically analyzed the effect of parameters of the measuring system (namely, the pump beam modulation frequency and pump beam radius) on the contrast and spatial resolution, specifically in relation to imaging a copper-resin double-layered structure that contains slits. This study also experimentally validated the representative result of numerical analysis.
The specimens of the copper-resin double-layered structure that were used in the numerical analysis and experiment are illustrated in Figure 4. The properties of the copper and resin are presented in Table 1. Additionally, Figure 5 presents the optical arrangement of the specimens and laser applied to the numerical analysis and experiment.

Setting the Study Parameters
The phase delay signals derived by the photothermal imaging method are significantly affected by the thermal diffusion length (L d ), which is expressed as Equation (3) and presented in units of length (m).
Here, α is thermal diffusivity (m 2 /s) of the material to be measured, and f is modulation frequency (Hz) of the pump beam.
As shown in Figure 5, the length variables of the photothermal imaging method that affect the photothermal effect and phase delay include the pump-beam radius (R b ), the thickness of the first layer (here, the copper layer, L f ) of the copper-resin double layer, the probe beam's relative position (X 0 ), and the z-position of the probe beam (L h ). Theoretically, if the ratio of each factor to the thermal diffusion length is the same, the same phase delay profile should appear when they are nondimensionalized on the basis of the thermal diffusion length (L d ) [1,2,24]. As the experimental results in Figure 6 show, the values of R b /L d , L f /L d , X 0 /L d , and L h /L d are the same in three conditions of different thermal diffusion lengths depending on the pump-beam frequency, as listed in Table 2. Thus, they formed the same scanning profile.

Setting the Study Parameters
The phase delay signals derived by the photothermal imaging method are significantly affected by the thermal diffusion length ( ), which is expressed as Equation (3) and presented in units of length (m).
Here, α is thermal diffusivity (m 2 /s) of the material to be measured, and is modulation frequency (Hz) of the pump beam.
As shown in Figure 5, the length variables of the photothermal imaging method that affect the photothermal effect and phase delay include the pump-beam radius ( ), the thickness of the first layer (here, the copper layer, ) of the copper-resin double layer, the probe beam's relative position ( 0 ), and the z-position of the probe beam ( ℎ ). Theoretically, if the ratio of each factor to the thermal diffusion length is the same, the same phase delay profile should appear when they are nondimensionalized on the basis of the thermal diffusion length ( ) [1,2,24]. As the experimental results in Figure 6 show, the values of / , / , 0 / , and ℎ / are the same in three conditions of different thermal diffusion lengths depending on the pump-beam frequency, as listed in Table 2. Thus, they formed the same scanning profile.  Among these factors, the relative position and z-position of the probe-beam ( ) are determined when a photothermal phenomenon is detected or phase delay data are acquired. These do not directly affect the phenomena of the photothermal effect. Hence, the present study derives image quality change related to two parameters: the first (copper) layer thickness ( ) of the copper-resin doublelayer and the pump-beam radius ( ).

Numerical Analysis
In order to analyze the image quality of the photothermal imaging method, numerical analysis was conducted on specimens with the simple subsurface structure shown in Figure 4. The numerical  Among these factors, the relative position and z-position of the probe-beam (L f ) are determined when a photothermal phenomenon is detected or phase delay data are acquired. These do not directly affect the phenomena of the photothermal effect. Hence, the present study derives image quality change related to two parameters: the first (copper) layer thickness (L f ) of the copper-resin double-layer and the pump-beam radius (R b ).

Numerical Analysis
In order to analyze the image quality of the photothermal imaging method, numerical analysis was conducted on specimens with the simple subsurface structure shown in Figure 4. The numerical analysis model had the boundary conditions presented in Figure 7, and it was conducted on the basis of the model established by Kim et al. [3]. Here, the infinite condition means that the physical properties of the adjacent material are infinite [25]. The size and shape of the numerical model specimen are the same as those in Figure 4, and the thickness of the top and bottom air layers is set to 15 mm. Equation (4) is the governing equation for the three areas that consist of three layers: the top air layer, the specimens, and the bottom air layer: where the subscript i indicates each domain to be interpreted (0: top air layer, 1: sample, 2: bottom air layer), T i is the temperature of the layer, k i is the thermal conductivity of the layer, t is the time, and Q i indicates the modulated heat source. Equation (5) is the heat source equation for each layer.
Here, the light absorption coefficient (λ, cm −1 ) of copper is 623,100 cm −1 , a m 2 is the area of the pump beam, and x c and y c are the center coordinates of the pump beam, respectively, which are changed to implement laser scanning. y c is fixed at 0 mm, and x c is changed from 0 to 10 mm (0.1 mm intervals), as shown in Figure 4.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 7 of 15 analysis model had the boundary conditions presented in Figure 7, and it was conducted on the basis of the model established by Kim et al. [3]. Here, the infinite condition means that the physical properties of the adjacent material are infinite [25]. The size and shape of the numerical model specimen are the same as those in Figure 4, and the thickness of the top and bottom air layers is set to 15 mm. Equation (4) is the governing equation for the three areas that consist of three layers: the top air layer, the specimens, and the bottom air layer: where the subscript indicates each domain to be interpreted (0: top air layer, 1: sample, 2: bottom air layer), is the temperature of the layer, is the thermal conductivity of the layer, is the time, and indicates the modulated heat source. Equation (5) is the heat source equation for each layer. Here, the light absorption coefficient ( , cm −1 ) of copper is 623,100 cm −1 , (m 2 ) is the area of the pump beam, and and are the center coordinates of the pump beam, respectively, which are changed to implement laser scanning.
is fixed at 0 mm, and is changed from 0 to 10 mm (0.1 mm intervals), as shown in Figure 4. The specific parameters used in the numerical analysis are presented in Table 3. With the other parameters fixed, the numerical analysis was conducted with five different values of the first (copper) layer thickness ( ) and four different values of the pump-beam radius ( ).  The specific parameters used in the numerical analysis are presented in Table 3. With the other parameters fixed, the numerical analysis was conducted with five different values of the first (copper) layer thickness (L f ) and four different values of the pump-beam radius (R b ).  Figure 8 presents the general configuration of the experimental equipment, which consisted mainly of the pump beam, specimen, probe-beam, and data-acquisition configurations. The light source from the pump beam to the probe beam was modulated by a mechanical chopper (SRS 40), passed through a focus lens with a focal length of 1000 mm, and then radiated onto the specimen surface. The intensity of the laser used in this process was 10 W, and a Sprout-G 10-W model with a wavelength of 532 nm, which is a diode-pumped solid-state (DPSS) laser with a Gaussian profile, was used. The radius of the pump beam is determined by the distance from the maximum intensity to 1/e 2 in the Gaussian distribution [19], and that radius is adjusted by the focus lens. The specimen to be measured was installed at the specimen section. The specimen was mounted on the Newport optic mount and then fixed on an x-y-axis automatic transport system (x-axis: M-IMS300CCHA; y-axis: M-ILS150PP) consisting of two automatic transport units. The probe beam section was located near the specimen section. The probe beam was an He-Ne laser (Newport 1125 P) with a wavelength of 632.8 nm, an intensity of 5 mW, and a diameter of 500 µm. The modulated refraction of the probe-beam that passed through the air layer on the heated specimen was detected by the position sensitive diodes (PSD) sensor (C10443-01). Thereafter, the reference modulation frequency set by the mechanical chopper and the probe-beam response frequency detected by the PSD sensor were collected by the lock-in amplifier (AMETEK 7270 DSP), and the phase difference between the two signals was calculated. Figure 5 presents the arrangement and relative distance between the probe-beam and pump beam.  Figure 8 presents the general configuration of the experimental equipment, which consisted mainly of the pump beam, specimen, probe-beam, and data-acquisition configurations. The light source from the pump beam to the probe beam was modulated by a mechanical chopper (SRS 40), passed through a focus lens with a focal length of 1000 mm, and then radiated onto the specimen surface. The intensity of the laser used in this process was 10 W, and a Sprout-G 10-W model with a wavelength of 532 nm, which is a diode-pumped solid-state (DPSS) laser with a Gaussian profile, was used. The radius of the pump beam is determined by the distance from the maximum intensity to 1/e 2 in the Gaussian distribution [19], and that radius is adjusted by the focus lens. The specimen to be measured was installed at the specimen section. The specimen was mounted on the Newport optic mount and then fixed on an x-y-axis automatic transport system (x-axis: M-IMS300CCHA; yaxis: M-ILS150PP) consisting of two automatic transport units. The probe beam section was located near the specimen section. The probe beam was an He-Ne laser (Newport 1125 P) with a wavelength of 632.8 nm, an intensity of 5 mW, and a diameter of 500 µ m. The modulated refraction of the probebeam that passed through the air layer on the heated specimen was detected by the position sensitive diodes (PSD) sensor (C10443-01). Thereafter, the reference modulation frequency set by the mechanical chopper and the probe-beam response frequency detected by the PSD sensor were collected by the lock-in amplifier (AMETEK 7270 DSP), and the phase difference between the two signals was calculated. Figure 5 presents the arrangement and relative distance between the probebeam and pump beam. The numerical analysis results were verified by experimentally examining the image quality depending on the / change. Specifically, experiments were conducted on specimens with value of 0.25 mm, 0.50 mm, and 1.00 mm (with / values of 0.362, 0.725, and 1.449, respectively), and the pump-beam radius ( ) was fixed at 100 µ m. The other parameters were set as shown in Table 3.

Results of Numerical Analysis
As mentioned above, the 1-D scanning results were numerically analyzed while changing the first layer thickness ( ) and the pump-beam radius ( ). In order to clarify common tendencies in a number of results according to each parameter change, changes in image quality were examined mainly for two groups: one group in which was fixed at 100 µ m ( / = 0.145) and another group in which was fixed at 0.5 mm ( / = 0.725). Additionally, every graph presented in this section was nondimensionalized on the basis of the thermal diffusion length.  Table 3.

Results of Numerical Analysis
As mentioned above, the 1-D scanning results were numerically analyzed while changing the first layer thickness (L f ) and the pump-beam radius (R b ). In order to clarify common tendencies in a number of results according to each parameter change, changes in image quality were examined mainly for two groups: one group in which R b was fixed at 100 µm (R b /L d = 0.145) and another group in which L f was fixed at 0.5 mm (L f /L d = 0.725). Additionally, every graph presented in this section was nondimensionalized on the basis of the thermal diffusion length.

Image Quality
Change with L f /L d Figure 9 shows the phase delay signals (the ERF on the internal structure) derived by the photothermal imaging method with changed various of L f /L d . As L f /L d decreased (closer to the surface), the response to the internal structural change became more intense. As L f /L d increased (i.e., farther from the surface), the response weakened.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 9 of 15 Figure 9 shows the phase delay signals (the ERF on the internal structure) derived by the photothermal imaging method with changed various of / . As / decreased (closer to the surface), the response to the internal structural change became more intense. As / increased (i.e., farther from the surface), the response weakened. The LSF can be obtained by differentiating the ERF as shown in Figure 10. Thus, the spatial resolution from the LSF can be analyzed. As mentioned previously, the spatial resolution of the imaging system is the FWHM of the LSF. Figure 9 indicates that, as / increases (located deeper), the FWHM decreases. In other words, the spatial resolution is improved. The spatial resolutions for five different values of / are summarized in Table 4.  The contrast was verified from the MTF derived by applying the FFT to the LSF. Figure 11 shows the contrast according to the change of / . As in the case of spatial resolution, the contrast increased together with / . In the case of the MTF, it is difficult to compare the general tendency on the basis of a single value. Instead, Table 4 presents the dimensionless spatial frequency (the The LSF can be obtained by differentiating the ERF as shown in Figure 10. Thus, the spatial resolution from the LSF can be analyzed. As mentioned previously, the spatial resolution of the imaging system is the FWHM of the LSF. Figure 9 indicates that, as L f /L d increases (located deeper), the FWHM decreases. In other words, the spatial resolution is improved. The spatial resolutions for five different values of L f /L d are summarized in Table 4.

Image Quality Change with /
Appl. Sci. 2019, 9, x FOR PEER REVIEW 9 of 15 Figure 9 shows the phase delay signals (the ERF on the internal structure) derived by the photothermal imaging method with changed various of / . As / decreased (closer to the surface), the response to the internal structural change became more intense. As / increased (i.e., farther from the surface), the response weakened. The LSF can be obtained by differentiating the ERF as shown in Figure 10. Thus, the spatial resolution from the LSF can be analyzed. As mentioned previously, the spatial resolution of the imaging system is the FWHM of the LSF. Figure 9 indicates that, as / increases (located deeper), the FWHM decreases. In other words, the spatial resolution is improved. The spatial resolutions for five different values of / are summarized in Table 4.  The contrast was verified from the MTF derived by applying the FFT to the LSF. Figure 11 shows the contrast according to the change of / . As in the case of spatial resolution, the contrast increased together with / . In the case of the MTF, it is difficult to compare the general tendency on the basis of a single value. Instead, Table 4 presents the dimensionless spatial frequency (the  The contrast was verified from the MTF derived by applying the FFT to the LSF. Figure 11 shows the contrast according to the change of L f /L d . As in the case of spatial resolution, the contrast increased together with L f /L d . In the case of the MTF, it is difficult to compare the general tendency on the basis of a single value. Instead, Table 4 presents the dimensionless spatial frequency (the thermal diffusion length multiplied by the spatial frequency) at which the MTF became 0.5 and the contrast decreased by half.

Image Quality Change with /
Appl. Sci. 2019, 9, x FOR PEER REVIEW 10 of 15 thermal diffusion length multiplied by the spatial frequency) at which the MTF became 0.5 and the contrast decreased by half. Figure 11. Performance of the MTF according to changes in / .
The changes in the spatial resolution and contrast according to / revealed that, as / increased, the spatial resolution and contrast increased accordingly. A small / means that the subsurface structure (resin structure) within the range of the thermal wave has a large influence, and a large / means that this influence is small. Figure 12 depicts the physical phenomena depending on / . As shown in Figure 12, when a thermal wave of the same size exists in the same position, the subsurface structure can be detected when / is small but cannot be detected when / is large. For this reason, as / increases, the FWHM of the LSF decreases and the resolution increases. However, the photothermal imaging signal should come from the phase delay, which reflects the distribution of the thermal properties within the range of the thermal wave. Therefore, if / is large, the influence of the subsurface structure on phase delay becomes small, so the intensity of the ERF decreases. Therefore, as / reached a certain level, the intensity of the ERF decreased significantly. Hence, it is estimated that experimental detection would be impossible, according to the signal noise analysis by the photothermal imaging method suggested by Kim et al. [2]. The changes in the spatial resolution and contrast according to L f /L d revealed that, as L f /L d increased, the spatial resolution and contrast increased accordingly. A small L f /L d means that the subsurface structure (resin structure) within the range of the thermal wave has a large influence, and a large L f /L d means that this influence is small. Figure 12 depicts the physical phenomena depending on L f /L d . As shown in Figure 12, when a thermal wave of the same size exists in the same position, the subsurface structure can be detected when L f /L d is small but cannot be detected when L f /L d is large. For this reason, as L f /L d increases, the FWHM of the LSF decreases and the resolution increases. However, the photothermal imaging signal should come from the phase delay, which reflects the distribution of the thermal properties within the range of the thermal wave. Therefore, if L f /L d is large, the influence of the subsurface structure on phase delay becomes small, so the intensity of the ERF decreases.

Image Quality Change with /
Appl. Sci. 2019, 9, x FOR PEER REVIEW 10 of 15 thermal diffusion length multiplied by the spatial frequency) at which the MTF became 0.5 and the contrast decreased by half. Figure 11. Performance of the MTF according to changes in / .
The changes in the spatial resolution and contrast according to / revealed that, as / increased, the spatial resolution and contrast increased accordingly. A small / means that the subsurface structure (resin structure) within the range of the thermal wave has a large influence, and a large / means that this influence is small. Figure 12 depicts the physical phenomena depending on / . As shown in Figure 12, when a thermal wave of the same size exists in the same position, the subsurface structure can be detected when / is small but cannot be detected when / is large. For this reason, as / increases, the FWHM of the LSF decreases and the resolution increases. However, the photothermal imaging signal should come from the phase delay, which reflects the distribution of the thermal properties within the range of the thermal wave. Therefore, if / is large, the influence of the subsurface structure on phase delay becomes small, so the intensity of the ERF decreases. Therefore, as / reached a certain level, the intensity of the ERF decreased significantly. Hence, it is estimated that experimental detection would be impossible, according to the signal noise analysis by the photothermal imaging method suggested by Kim et al. [2]. Therefore, as L f /L d reached a certain level, the intensity of the ERF decreased significantly. Hence, it is estimated that experimental detection would be impossible, according to the signal noise analysis by the photothermal imaging method suggested by Kim et al. [2]. Figure 13 shows the ERF according to changes in R b /L d . Unlike the case of L f /L d , there was a minor change in the ERF intensity as R b /L d increased in the case of R b /L d . Because the ERF change was insignificant, the LSF also showed no major change. As a result, the spatial resolution change derived from the LSF, as shown in Figure 14, showed a difference of only 10% from the average. The MTF indicating the contrast in Figure 15 also showed the same tendency. In addition, the positions where the contrast decreased by half showed a difference of less than 10% from the average. The reason for this is that the diameter of the thermal wave in the radial direction cannot have influences on the detection of subsurface structure.  Numerical analysis was conducted for each of the 20 cases in which both / and / changed, as well as the cases in which each parameter was fixed, and the results are presented in Because the ERF change was insignificant, the LSF also showed no major change. As a result, the spatial resolution change derived from the LSF, as shown in Figure 14, showed a difference of only 10% from the average. The MTF indicating the contrast in Figure 15 also showed the same tendency. In addition, the positions where the contrast decreased by half showed a difference of less than 10% from the average. The reason for this is that the diameter of the thermal wave in the radial direction cannot have influences on the detection of subsurface structure. Because the ERF change was insignificant, the LSF also showed no major change. As a result, the spatial resolution change derived from the LSF, as shown in Figure 14, showed a difference of only 10% from the average. The MTF indicating the contrast in Figure 15 also showed the same tendency. In addition, the positions where the contrast decreased by half showed a difference of less than 10% from the average. The reason for this is that the diameter of the thermal wave in the radial direction cannot have influences on the detection of subsurface structure.  Numerical analysis was conducted for each of the 20 cases in which both / and / changed, as well as the cases in which each parameter was fixed, and the results are presented in Table 5. Likewise, the image quality change according to changes in / was insignificant. In Because the ERF change was insignificant, the LSF also showed no major change. As a result, the spatial resolution change derived from the LSF, as shown in Figure 14, showed a difference of only 10% from the average. The MTF indicating the contrast in Figure 15 also showed the same tendency. In addition, the positions where the contrast decreased by half showed a difference of less than 10% from the average. The reason for this is that the diameter of the thermal wave in the radial direction cannot have influences on the detection of subsurface structure.  Numerical analysis was conducted for each of the 20 cases in which both / and / changed, as well as the cases in which each parameter was fixed, and the results are presented in Table 5. Likewise, the image quality change according to changes in / was insignificant. In Numerical analysis was conducted for each of the 20 cases in which both L f /L d and R b /L d changed, as well as the cases in which each parameter was fixed, and the results are presented in Table 5. Likewise, the image quality change according to changes in R b /L d was insignificant. In contrast, for the image quality change according to changes in L f /L d , the spatial resolution increased as L f /L d increased. Table 5. Image quality depending on changes in L f and R b . This tendency suggests that the pump-beam radius has an insignificant effect on image quality and that an increase in L f /L d within an experimentally measurable range is essential to achieve good image quality. Figure 16 shows the ERF measured in the experiments in comparison with the numerical analysis results. It was found that the numerical analysis results corresponded well with the experimental results, with a determination coefficient (R 2 ) of 0.974. Noise was removed from the experimental results by means of the low-pass fast Fourier transform (FFT) filter proposed by Kim et al. The LSF (Figure 17) was then derived, and the MTF (Figure 18) was derived from the FFT in order to analyze the image quality. However, when the ratio of L f /L d was 1.449, it was difficult to obtain a smoothed ERF through the removal of noise. This is because the intensity of the noise is larger than intensity of the ERF, as shown in Figure 18. Therefore, it was impossible to carry out the analysis of image quality. For this reason, Table 6 does not indicate the spatial resolution and contrast for the L f /L d of 1.449. When the numerical analysis results were compared with the experimental results for the final spatial resolution and contrast according to each L f /L d value, the difference was within 8% in all cases. contrast, for the image quality change according to changes in / , the spatial resolution increased as / increased. This tendency suggests that the pump-beam radius has an insignificant effect on image quality and that an increase in / within an experimentally measurable range is essential to achieve good image quality. Figure 16 shows the ERF measured in the experiments in comparison with the numerical analysis results. It was found that the numerical analysis results corresponded well with the experimental results, with a determination coefficient ( 2 ) of 0.974. Noise was removed from the experimental results by means of the low-pass fast Fourier transform (FFT) filter proposed by Kim et al. The LSF (Figure 17) was then derived, and the MTF (Figure 18) was derived from the FFT in order to analyze the image quality. However, when the ratio of / was 1.449, it was difficult to obtain a smoothed ERF through the removal of noise. This is because the intensity of the noise is larger than intensity of the ERF, as shown in Figure 18. Therefore, it was impossible to carry out the analysis of image quality. For this reason, Table 6 does not indicate the spatial resolution and contrast for the / of 1.449. When the numerical analysis results were compared with the experimental results for the final spatial resolution and contrast according to each / value, the difference was within 8% in all cases.   For various L f / L d Figure 16. ERF comparison between the experimental and numerical analysis results.   1% -

Conclusions
In this study, the image quality changed according to different measuring conditions, which included the depth of a subsurface structure from the surface of the object (first layer thickness, ) and the pump-beam radius ( ). As image quality metrics, the spatial resolution and contrast were derived from the FWHM of the LSF and the MTF.
1. It was verified that the phase delay profile (edge response function, ERF) of the subsurface structure derived by the photothermal imaging method could be nondimensionalized in relation to the thermal diffusion length. Therefore, various studies on / and / can be conducted   1% -

Conclusions
In this study, the image quality changed according to different measuring conditions, which included the depth of a subsurface structure from the surface of the object (first layer thickness, ) and the pump-beam radius ( ). As image quality metrics, the spatial resolution and contrast were derived from the FWHM of the LSF and the MTF.
1. It was verified that the phase delay profile (edge response function, ERF) of the subsurface structure derived by the photothermal imaging method could be nondimensionalized in relation to the thermal diffusion length. Therefore, various studies on / and / can be conducted

Conclusions
In this study, the image quality changed according to different measuring conditions, which included the depth of a subsurface structure from the surface of the object (first layer thickness, L f ) and the pump-beam radius (R b ). As image quality metrics, the spatial resolution and contrast were derived from the FWHM of the LSF and the MTF.

1.
It was verified that the phase delay profile (edge response function, ERF) of the subsurface structure derived by the photothermal imaging method could be nondimensionalized in relation to the thermal diffusion length. Therefore, various studies on L f /L d and R b /L d can be conducted to quantitatively analyze the spatial resolution and contrast of photothermal imaging regardless of the absolute size and thermal properties of the material to be measured.

2.
Analyzing the spatial resolution and contrast under different measuring conditions revealed that, as the ratio of L f /L d increases, the spatial resolution and contrast increase as well. On the other hand, the ratio of R b /L d has an insignificant effect on the spatial resolution and contrast. 3.
The range of thermal waves that vary with L d has a significant effect on the spatial resolution and contrast of photothermal imaging. L d can be adjusted by the thermal diffusivity of the material to be measured and the modulation frequency of the pump beam as shown in Equation (3). Therefore, to increase the spatial resolution and contrast, it is advantageous for the photothermal imaging method to select L d , which is relatively short compared to L f .

4.
As the L f /L d increases, however, the influence of the subsurface structure within the range of the thermal wave decreases. This means that the phase delay change may not reflect the subsurface structure. Therefore, the photothermal imaging method in large L f /L d is not effective for the measurement. 5.
In conclusion, the image quality analysis of this paper is expected to provide a validation method for the results of photothermal imaging and effective measurement conditions for the photothermal imaging method.
These results suggest that the experimentally measurable spatial resolution and contrast can be improved if the signal-to-noise ratio of the photothermal imaging system can be increased. The noise of a photothermal signal can be controlled by the surface roughness of the specimen to be measured and the pump-beam intensity of the photothermal imaging system [2]. Therefore, in future research, image quality can be analyzed according to the magnitude of noise in the photothermal imaging method.