<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD with MathML3 v1.1d2 20140930//EN" "JATS-journalpublishing1-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="1.1d2" xml:lang="en">
  <front>
    <journal-meta>
      <journal-id journal-id-type="nlm-ta">JIAP</journal-id>
      <journal-id journal-id-type="publisher-id">IECE</journal-id>
      <journal-title-group>
        <journal-title>IECE Journal of Image Analysis and Processing</journal-title>
      </journal-title-group>
      <issn pub-type="ppub" publication-format="print">pending</issn>
      <issn pub-type="epub" publication-format="electronic">pending</issn>
      <publisher>
        <publisher-name>Institute of Emerging and Computer Engineering Inc</publisher-name>
        <publisher-loc>522 W RIVERSIDE AVE STE N, SPOKANE, WA, 99201-0508, UNITED STATES</publisher-loc>
      </publisher>
    </journal-meta>
    <article-meta>
      <article-id pub-id-type="doi">10.62762/JIAP.2024.764051</article-id>
      <article-categories>
        <subj-group subj-group-type="heading">
          <subject>Research Article</subject>
        </subj-group>
      </article-categories>
      <title-group>
        <article-title>High-Quality Multi-Focus Image Fusion: A Comparative Analysis of DCT-Based Approaches with Their Variants</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <name>
            <surname>Osama</surname>
            <given-names>Muhammad</given-names>
          </name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <name>
            <surname>Khan</surname>
            <given-names>Sarwar Shah</given-names>
          </name>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <contrib-id contrib-id-type="orcid">https://orcid.org/0009-0001-2556-7068</contrib-id>
          <name>
            <surname>Khan</surname>
            <given-names>Sajid</given-names>
          </name>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <name>
            <surname>Ahmad</surname>
            <given-names>Shehzad</given-names>
          </name>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <contrib-id contrib-id-type="orcid">https://orcid.org/0000-0002-4216-1246</contrib-id>
          <name>
            <surname>Mehmood</surname>
            <given-names>Gulzar</given-names>
          </name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <contrib-id contrib-id-type="orcid">https://orcid.org/0000-0003-4556-5606</contrib-id>
          <name>
            <surname>Ali</surname>
            <given-names>Ijaz</given-names>
          </name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
      </contrib-group>
      <author-notes>
        <corresp id="cor2">Corresponding Author: Sarwar Shah Khan. Email: <email>sskhan0092@gmail.com</email></corresp>
      </author-notes>
      <pub-date date-type="pub" pub-type="epub" publication-format="online">
        <day>14</day>
        <month>3</month>
        <year>2025</year>
      </pub-date>
      <volume>1</volume>
      <issue>1</issue>
      <fpage>27</fpage>
      <lpage>35</lpage>
      <history>
        <date date-type="received">
          <day>21</day>
          <month>11</month>
          <year>2024</year>
        </date>
        <date date-type="accepted">
          <day>24</day>
          <month>12</month>
          <year>2024</year>
        </date>
      </history>
      <permissions>
        <copyright-statement>© 2025 by the Authors. Published by Institute of Emerging and Computer Engineers. This is an open access article under the CC BY license (https://creativecommons.org/licenses/by/4.0/).</copyright-statement>
        <copyright-year>2025</copyright-year>
        <copyright-holder>Institute of Emerging and Computer Engineering Inc</copyright-holder>
        <license xlink:href="https://creativecommons.org/licenses/by/4.0/">
        <license-p>This work is licensed under a <ext-link ext-link-type="uri" xlink:type="simple" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution 4.0 International License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
        </license>
      </permissions>
      <self-uri xlink:href="https://www.iece.org/article/abs/jiap.2024.764051">This article is available from https://www.iece.org/article/abs/jiap.2024.764051</self-uri>
      <abstract>
        <p>Image fusion, especially in the context of multi-focus image fusion, plays a crucial role in digital image processing by enhancing the clarity and detail of visual content through the combination of multiple source images. Traditional spatial domain methods often suffer from issues like spectral distortion and low contrast, which has led researchers to explore techniques in the frequency domain, such as the Discrete Cosine Transform (DCT). DCT-based methods are particularly valued for their computational efficiency, making them a strong alternative, especially in applications like image compression and fusion. This study focuses on DCT-based approaches, including variants that incorporate Singular Value Decomposition (SVD) and a combination of Correlation Coefficient with Energy-Correlation (Corr_Eng), both with and without Consistency Verification (CV). Extensive testing on multi-focus image datasets revealed that the DCT + SVD + CV method consistently shows better results in both qualitative and quantitative assessments. This indicates that integrating DCT+SVD+CV provides a powerful approach for achieving effective and efficient image fusion.</p>
      </abstract>
      <kwd-group kwd-group-type="author" xml:lang="en">
        <kwd>multi-focused</kwd>
        <kwd>image fusion</kwd>
        <kwd>discrete cosine transform</kwd>
        <kwd>spatial domain and frequency domain approaches</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <p>
      <fig id="fig1"/>
    </p>
    <sec id="S1">
      <label>1.</label>
      <title>Introduction</title>
      <p id="S1.p1">Over the past two decades, image fusion has become one of the most significant research fields in digital image processing, leading to the development and implementation of numerous approaches aimed at enhancing accuracy. This process involves merging two source images into a single, resultant image that conveys more meaningful and informative content than any individual source image could provide on its own. In some cases, the required information can only be revealed by merging a couple of images. High-quality visuals are essential in many such as security, computer vision [<xref rid="ref001" ref-type="bibr">1</xref>], medical, military [<xref rid="ref002" ref-type="bibr">2</xref>], remote sensing [<xref rid="ref003" ref-type="bibr">3</xref>], navigation guidance for pilots, and weather forecasting [<xref rid="ref004" ref-type="bibr">4</xref>].</p>
      <p id="S1.p2">Multi-focus imaging is one of the essential types of image fusion that has seen extensive research interest over the last few decades. Optical lenses are limited by their depth of field, meaning that only objects at a specific distance from the lens will appear sharp and in focus. As a result, in any given image, only one object will be in focus, while another object at a different distance from the lens will be out of focus and, hence blurred. Several factors contribute to the extent of this blurring, including the distance from the object, the focal length, the number of lenses used, and the distance between the lens and the sensor plane [<xref rid="ref005" ref-type="bibr">5</xref>].</p>
      <p id="S1.p3">Multi-focus image fusion necessitates a variety of traditional and advanced approaches to produce a more informative resultant image. These approaches can be broadly categorized into a couple of classes such as spatial domain and frequency domain approaches. In the spatial domain, to deal with the image in its original form, meaning that the pixel values are directly manipulated based on the scene. On the other hand, frequency domain approaches focus on the rate at which pixel values change in the spatial domain.</p>
      <p id="S1.p4">Spatial domain approaches include averaging methods, Principal Component Analysis (PCA), simple maximum and minimum methods [<xref rid="ref006" ref-type="bibr">6</xref>], and Intensity Hue Saturation (IHS) [<xref rid="ref007" ref-type="bibr">7</xref>]. However, these approaches often produce subpar results due to spectral distortions, leading to low-contrast images with less information [<xref rid="ref004" ref-type="bibr">4</xref>]. Additionally, the spatial domain doesn’t provide enough robustness and perceivable [<xref rid="ref008" ref-type="bibr">8</xref>]. While, frequency domain approaches encompass techniques like pyramid transform, Discrete Wavelet Transform (DWT), Stationary Wavelet Transform (SWT) [<xref rid="ref009" ref-type="bibr">9</xref>], and Discrete Cosine Transform (DCT) [<xref rid="ref010" ref-type="bibr">10</xref>]. The benefits of the frequency domain are minimal processing complexity, ease of view, the ability to manipulate the frequency composition and special transformed domain property is easily applicable [<xref rid="ref008" ref-type="bibr">8</xref>].</p>
      <p>
        <table-wrap id="S1.tab1">
          <label>Table 1</label>
          <caption>
            <p>Advantages and disadvantages of different multi-focus image fusion.</p>
          </caption>
          <table>
            <tr>
              <td style="border-top: 1px solid black;" align="center">
                <bold>Fusion Method</bold>
              </td>
              <td style="border-top: 1px solid black;" align="center">
                <bold>Advantages</bold>
              </td>
              <td style="border-top: 1px solid black;" align="center">
                <bold>Disadvantages</bold>
              </td>
            </tr>
            <tr>
              <td style="border-top: 1px solid black;" align="justify">
                <p>[c]DWT</p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I1">
                      <list-item id="S1.I1.i1">
                        <p id="S1.I1.i1.p1">DWT is an effective multi-focus image fusion approach.</p>
                      </list-item>
                      <list-item id="S1.I1.i2">
                        <p id="S1.I1.i2.p1">DWT allows the decomposition of an image into various frequency sub-bands, enabling detailed analysis at different scales.</p>
                      </list-item>
                      <list-item id="S1.I1.i3">
                        <p id="S1.I1.i3.p1">The DWT decreases the spectral distortion in an image.</p>
                      </list-item>
                      <list-item id="S1.I1.i4">
                        <p id="S1.I1.i4.p1">It is particularly effective in representing image edges and textures.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I2">
                      <list-item id="S1.I2.i1">
                        <p id="S1.I2.i1.p1">Artifacts may be introduced during the decomposition and reconstruction process.</p>
                      </list-item>
                      <list-item id="S1.I2.i2">
                        <p id="S1.I2.i2.p1">Due to DWT’s sensitivity to changes in the input image, the fused image may be inaccurate and misaligned.</p>
                      </list-item>
                      <list-item id="S1.I2.i3">
                        <p id="S1.I2.i3.p1">DWT preserves only the vertical and horizontal properties.</p>
                      </list-item>
                      <list-item id="S1.I2.i4">
                        <p id="S1.I2.i4.p1">It experiences anomalous ringing, which reduces the resolution of the resulting image.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
            </tr>
            <tr>
              <td style="border-top: 1px solid black;" align="justify">
                <p>[c]SWT</p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I3">
                      <list-item id="S1.I3.i1">
                        <p id="S1.I3.i1.p1">SWT produces more accurate fusion results than DWT because it is shift-invariant.</p>
                      </list-item>
                      <list-item id="S1.I3.i2">
                        <p id="S1.I3.i2.p1">Regarding multi-focus picture fusion, SWT performs better at maintaining image edges and details information.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I4">
                      <list-item id="S1.I4.i1">
                        <p id="S1.I4.i1.p1">In real-time applications, the redundancy in SWT can be undesirable since it increases memory usage and computational costs.</p>
                      </list-item>
                      <list-item id="S1.I4.i2">
                        <p id="S1.I4.i2.p1">The redundancy can sometimes lead to overfitting.</p>
                      </list-item>
                      <list-item id="S1.I4.i3">
                        <p id="S1.I4.i3.p1">SWT is less efficient.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
            </tr>
            <tr>
              <td style="border-top: 1px solid black;" align="justify">
                <p>[c]DCT</p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I5">
                      <list-item id="S1.I5.i1">
                        <p id="S1.I5.i1.p1">DCT efficiently compacts most of the signal’s energy into a few low-frequency components.</p>
                      </list-item>
                      <list-item id="S1.I5.i2">
                        <p id="S1.I5.i2.p1">DCT is computationally less expensive.</p>
                      </list-item>
                      <list-item id="S1.I5.i3">
                        <p id="S1.I5.i3.p1">It is faster and easier to implement.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I6">
                      <list-item id="S1.I6.i1">
                        <p id="S1.I6.i1.p1">In block-based DCT, block artifacts may occur.</p>
                      </list-item>
                      <list-item id="S1.I6.i2">
                        <p id="S1.I6.i2.p1">It’s possible that DCT is less successful in maintaining high-frequency elements.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
            </tr>
            <tr>
              <td style="border-top: 1px solid black;" align="justify">
                <p>[c]DCT + SVD</p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I7">
                      <list-item id="S1.I7.i1">
                        <p id="S1.I7.i1.p1">By combining the advantages of both methods—the energy compaction of DCT and the major feature extraction of SVD—detail preservation is improved.</p>
                      </list-item>
                      <list-item id="S1.I7.i2">
                        <p id="S1.I7.i2.p1">This approach is resistant to noise and small image distortions.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I8">
                      <list-item id="S1.I8.i1">
                        <p id="S1.I8.i1.p1">Adding SVD to DCT increases the computational complexity.</p>
                      </list-item>
                      <list-item id="S1.I8.i2">
                        <p id="S1.I8.i2.p1">The performance heavily depends on the selection of singular values.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
            </tr>
            <tr>
              <td style="border-top: 1px solid black;" align="justify">
                <p>[c] DCT + SVD +  Correlation  Coefficient (CV)</p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I9">
                      <list-item id="S1.I9.i1">
                        <p id="S1.I9.i1.p1">The correlation coefficient helps to ensure that the fused image maintains the structural similarity to the original images.</p>
                      </list-item>
                      <list-item id="S1.I9.i2">
                        <p id="S1.I9.i2.p1">This combination is effective in fusing images with complementary information.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I10">
                      <list-item id="S1.I10.i1">
                        <p id="S1.I10.i1.p1">The fusion process is made more difficult by the requirement to optimize several parameters.</p>
                      </list-item>
                      <list-item id="S1.I10.i2">
                        <p id="S1.I10.i2.p1">As with any multi-phase procedure, overfitting is a possibility.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
            </tr>
            <tr>
              <td style="border-top: 1px solid black;" align="justify">
                <p>[c] DCT + Correlation  Coefficient and  Energy-Correlation</p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I11">
                      <list-item id="S1.I11.i1">
                        <p id="S1.I11.i1.p1">By balancing the correlation and energy information between images, this approach creates a fused image that preserves structural details as well as energy properties.</p>
                      </list-item>
                      <list-item id="S1.I11.i2">
                        <p id="S1.I11.i2.p1">The technique is flexible for a range of fusion tasks.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
              <td style="border-top: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I12">
                      <list-item id="S1.I12.i1">
                        <p id="S1.I12.i1.p1">The overall complexity is increased when energy and correlation metrics are combined with DCT.</p>
                      </list-item>
                      <list-item id="S1.I12.i2">
                        <p id="S1.I12.i2.p1">This technique could need a lot of processing power.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
            </tr>
            <tr>
              <td style="border-top: 1px solid black;border-bottom: 1px solid black;" rowspan="2" align="justify">
                <p>[c] DCT + Correlation  Coefficient and  Energy-Correlation + CV</p>
              </td>
              <td style="border-top: 1px solid black;border-bottom: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I13">
                      <list-item id="S1.I13.i1">
                        <p id="S1.I13.i1.p1">This method maximizes the retention of structural and statistical information from the input images.</p>
                      </list-item>
                      <list-item id="S1.I13.i2">
                        <p id="S1.I13.i2.p1">The fused image is likely to be of high accuracy, maintaining the essential features and details of the original images.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
              <td style="border-top: 1px solid black;border-bottom: 1px solid black;" align="justify">
                <p>
                  <p>
                    <list list-type="bullet" id="S1.I14">
                      <list-item id="S1.I14.i1">
                        <p id="S1.I14.i1.p1">The integration of multiple techniques leads to a significant increase in computational requirements.</p>
                      </list-item>
                      <list-item id="S1.I14.i2">
                        <p id="S1.I14.i2.p1">The complexity of integrating and tuning multiple techniques can make the implementation challenging.</p>
                      </list-item>
                    </list>
                  </p>
                </p>
              </td>
            </tr>
          </table>
          <sec id="S2">
            <label>2.</label>
            <title>Multi-focus Image fusion Approaches</title>
            <p id="S2.p1">Various approaches have been developed for both spatial and frequency domains in the area of multi-focus image fusion. Additionally, the frequency domain offers more advantages over the spatial domain [<xref rid="ref008" ref-type="bibr">8</xref>]. Hence, this study mainly focuses on the DCT approach. Its purpose is to highlight the maximum effective approaches by analyzing their characteristics and quality, as well as presenting experimental results obtained from image sets.</p>
            <sec id="S2.SS1">
              <label>2.1</label>
              <title>Discrete Cosine Transform</title>
              <p id="S2.SS1.p1">The DCT approach facilitates the transformation from the spatial domain to the frequency domain, making it possible to extract detailed and outline information from an image based on pixel frequencies. DCT is an effective approach for handling frequencies, offering a fast and straightforward solution by utilizing only cosine functions for the transformation. The Inverse Discrete Cosine Transform (IDCT) can then be used to reconstruct the original pixel values from the frequencies derived through DCT [<xref rid="ref020" ref-type="bibr">20</xref>]. A finite sequence of data points is significantly represented by DCT as the sum of cosine functions that oscillate at different frequencies [<xref rid="ref021" ref-type="bibr">21</xref>]. The process of DCT evaluation is as follows.</p>
              <p id="S2.SS1.p2">The two-dimensional DCT transform of an N×N (usually 8×8) block of an image <inline-formula><mml:math alttext="x" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> (m,n) and the inverse DCT (IDCT) are defined in Eq.(1) and (3) respectively:</p>
              <p>
                <disp-formula-group id="S4.EGx1">
                  <disp-formula id="S2.Ex1">
                    <mml:math alttext="\displaystyle d(k,l)=\frac{2a(k)a(l)}{N}\sum_{m=0}^{N-1}\sum_{n=0}^{N-1}x(m,n)%&#10;\cos\left[\frac{(2m+1)\pi k}{2n}\right]" display="inline">
                      <mml:mrow>
                        <mml:mrow>
                          <mml:mi>d</mml:mi>
                          <mml:mo>⁢</mml:mo>
                          <mml:mrow>
                            <mml:mo stretchy="false">(</mml:mo>
                            <mml:mi>k</mml:mi>
                            <mml:mo>,</mml:mo>
                            <mml:mi>l</mml:mi>
                            <mml:mo stretchy="false">)</mml:mo>
                          </mml:mrow>
                        </mml:mrow>
                        <mml:mo>=</mml:mo>
                        <mml:mrow>
                          <mml:mstyle displaystyle="true">
                            <mml:mfrac>
                              <mml:mrow>
                                <mml:mn>2</mml:mn>
                                <mml:mo>⁢</mml:mo>
                                <mml:mi>a</mml:mi>
                                <mml:mo>⁢</mml:mo>
                                <mml:mrow>
                                  <mml:mo stretchy="false">(</mml:mo>
                                  <mml:mi>k</mml:mi>
                                  <mml:mo stretchy="false">)</mml:mo>
                                </mml:mrow>
                                <mml:mo>⁢</mml:mo>
                                <mml:mi>a</mml:mi>
                                <mml:mo>⁢</mml:mo>
                                <mml:mrow>
                                  <mml:mo stretchy="false">(</mml:mo>
                                  <mml:mi>l</mml:mi>
                                  <mml:mo stretchy="false">)</mml:mo>
                                </mml:mrow>
                              </mml:mrow>
                              <mml:mi>N</mml:mi>
                            </mml:mfrac>
                          </mml:mstyle>
                          <mml:mo>⁢</mml:mo>
                          <mml:mrow>
                            <mml:mstyle displaystyle="true">
                              <mml:munderover>
                                <mml:mo movablelimits="false">∑</mml:mo>
                                <mml:mrow>
                                  <mml:mi>m</mml:mi>
                                  <mml:mo>=</mml:mo>
                                  <mml:mn>0</mml:mn>
                                </mml:mrow>
                                <mml:mrow>
                                  <mml:mi>N</mml:mi>
                                  <mml:mo>−</mml:mo>
                                  <mml:mn>1</mml:mn>
                                </mml:mrow>
                              </mml:munderover>
                            </mml:mstyle>
                            <mml:mrow>
                              <mml:mstyle displaystyle="true">
                                <mml:munderover>
                                  <mml:mo movablelimits="false">∑</mml:mo>
                                  <mml:mrow>
                                    <mml:mi>n</mml:mi>
                                    <mml:mo>=</mml:mo>
                                    <mml:mn>0</mml:mn>
                                  </mml:mrow>
                                  <mml:mrow>
                                    <mml:mi>N</mml:mi>
                                    <mml:mo>−</mml:mo>
                                    <mml:mn>1</mml:mn>
                                  </mml:mrow>
                                </mml:munderover>
                              </mml:mstyle>
                              <mml:mrow>
                                <mml:mi>x</mml:mi>
                                <mml:mo>⁢</mml:mo>
                                <mml:mrow>
                                  <mml:mo stretchy="false">(</mml:mo>
                                  <mml:mi>m</mml:mi>
                                  <mml:mo>,</mml:mo>
                                  <mml:mi>n</mml:mi>
                                  <mml:mo stretchy="false">)</mml:mo>
                                </mml:mrow>
                                <mml:mo lspace="0.167em">⁢</mml:mo>
                                <mml:mrow>
                                  <mml:mi>cos</mml:mi>
                                  <mml:mo>⁡</mml:mo>
                                  <mml:mrow>
                                    <mml:mo>[</mml:mo>
                                    <mml:mstyle displaystyle="true">
                                      <mml:mfrac>
                                        <mml:mrow>
                                          <mml:mrow>
                                            <mml:mo stretchy="false">(</mml:mo>
                                            <mml:mrow>
                                              <mml:mrow>
                                                <mml:mn>2</mml:mn>
                                                <mml:mo>⁢</mml:mo>
                                                <mml:mi>m</mml:mi>
                                              </mml:mrow>
                                              <mml:mo>+</mml:mo>
                                              <mml:mn>1</mml:mn>
                                            </mml:mrow>
                                            <mml:mo stretchy="false">)</mml:mo>
                                          </mml:mrow>
                                          <mml:mo>⁢</mml:mo>
                                          <mml:mi>π</mml:mi>
                                          <mml:mo>⁢</mml:mo>
                                          <mml:mi>k</mml:mi>
                                        </mml:mrow>
                                        <mml:mrow>
                                          <mml:mn>2</mml:mn>
                                          <mml:mo>⁢</mml:mo>
                                          <mml:mi>n</mml:mi>
                                        </mml:mrow>
                                      </mml:mfrac>
                                    </mml:mstyle>
                                    <mml:mo>]</mml:mo>
                                  </mml:mrow>
                                </mml:mrow>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:mrow>
                      </mml:mrow>
                    </mml:math>
                  </disp-formula>
                  <disp-formula id="S2.E1">
                    <mml:math alttext="\displaystyle\quad\times\cos\left[\frac{(2n+1)\pi l}{2n}\right]" display="inline">
                      <mml:mrow>
                        <mml:mi/>
                        <mml:mo lspace="1.222em" rspace="0.222em">×</mml:mo>
                        <mml:mrow>
                          <mml:mi>cos</mml:mi>
                          <mml:mo>⁡</mml:mo>
                          <mml:mrow>
                            <mml:mo>[</mml:mo>
                            <mml:mstyle displaystyle="true">
                              <mml:mfrac>
                                <mml:mrow>
                                  <mml:mrow>
                                    <mml:mo stretchy="false">(</mml:mo>
                                    <mml:mrow>
                                      <mml:mrow>
                                        <mml:mn>2</mml:mn>
                                        <mml:mo>⁢</mml:mo>
                                        <mml:mi>n</mml:mi>
                                      </mml:mrow>
                                      <mml:mo>+</mml:mo>
                                      <mml:mn>1</mml:mn>
                                    </mml:mrow>
                                    <mml:mo stretchy="false">)</mml:mo>
                                  </mml:mrow>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mi>π</mml:mi>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mi>l</mml:mi>
                                </mml:mrow>
                                <mml:mrow>
                                  <mml:mn>2</mml:mn>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mi>n</mml:mi>
                                </mml:mrow>
                              </mml:mfrac>
                            </mml:mstyle>
                            <mml:mo>]</mml:mo>
                          </mml:mrow>
                        </mml:mrow>
                      </mml:mrow>
                    </mml:math>
                  </disp-formula>
                </disp-formula-group>
              </p>
              <p>where <inline-formula><mml:math alttext="k,l=0,1,\ldots,N-1" display="inline"><mml:mrow><mml:mrow><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>l</mml:mi></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi mathvariant="normal">…</mml:mi><mml:mo>,</mml:mo><mml:mrow><mml:mi>N</mml:mi><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:mrow></mml:mrow></mml:math></inline-formula> and</p>
              <p>
                <disp-formula id="S2.E2">
                  <mml:math alttext="a(k)=\begin{cases}\frac{1}{\sqrt{2}},&amp;\text{if }k=0\\&#10;1,&amp;\text{otherwise}\end{cases}" display="block">
                    <mml:mrow>
                      <mml:mrow>
                        <mml:mi>a</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:mrow>
                          <mml:mo stretchy="false">(</mml:mo>
                          <mml:mi>k</mml:mi>
                          <mml:mo stretchy="false">)</mml:mo>
                        </mml:mrow>
                      </mml:mrow>
                      <mml:mo>=</mml:mo>
                      <mml:mrow>
                        <mml:mo>{</mml:mo>
                        <mml:mtable columnspacing="5pt" displaystyle="true" rowspacing="0pt">
                          <mml:mtr>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mrow>
                                <mml:mstyle displaystyle="false">
                                  <mml:mfrac>
                                    <mml:mn>1</mml:mn>
                                    <mml:msqrt>
                                      <mml:mn>2</mml:mn>
                                    </mml:msqrt>
                                  </mml:mfrac>
                                </mml:mstyle>
                                <mml:mo>,</mml:mo>
                              </mml:mrow>
                            </mml:mtd>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mrow>
                                <mml:mrow>
                                  <mml:mtext>if </mml:mtext>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mi>k</mml:mi>
                                </mml:mrow>
                                <mml:mo>=</mml:mo>
                                <mml:mn>0</mml:mn>
                              </mml:mrow>
                            </mml:mtd>
                          </mml:mtr>
                          <mml:mtr>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mrow>
                                <mml:mn>1</mml:mn>
                                <mml:mo>,</mml:mo>
                              </mml:mrow>
                            </mml:mtd>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mtext>otherwise</mml:mtext>
                            </mml:mtd>
                          </mml:mtr>
                        </mml:mtable>
                      </mml:mrow>
                    </mml:mrow>
                  </mml:math>
                </disp-formula>
              </p>
              <p>
                <disp-formula-group id="S2.E3">
                  <disp-formula id="S2.E3X">
                    <mml:math alttext="\displaystyle x(m,n)=\sum_{k=0}^{N-1}\sum_{j=0}^{N-1}\frac{2a(k)a(l)}{N}\times&#10;d%&#10;(k,l)" display="inline">
                      <mml:mrow>
                        <mml:mrow>
                          <mml:mi>x</mml:mi>
                          <mml:mo>⁢</mml:mo>
                          <mml:mrow>
                            <mml:mo stretchy="false">(</mml:mo>
                            <mml:mi>m</mml:mi>
                            <mml:mo>,</mml:mo>
                            <mml:mi>n</mml:mi>
                            <mml:mo stretchy="false">)</mml:mo>
                          </mml:mrow>
                        </mml:mrow>
                        <mml:mo>=</mml:mo>
                        <mml:mrow>
                          <mml:mstyle displaystyle="true">
                            <mml:munderover>
                              <mml:mo movablelimits="false">∑</mml:mo>
                              <mml:mrow>
                                <mml:mi>k</mml:mi>
                                <mml:mo>=</mml:mo>
                                <mml:mn>0</mml:mn>
                              </mml:mrow>
                              <mml:mrow>
                                <mml:mi>N</mml:mi>
                                <mml:mo>−</mml:mo>
                                <mml:mn>1</mml:mn>
                              </mml:mrow>
                            </mml:munderover>
                          </mml:mstyle>
                          <mml:mrow>
                            <mml:mstyle displaystyle="true">
                              <mml:munderover>
                                <mml:mo movablelimits="false">∑</mml:mo>
                                <mml:mrow>
                                  <mml:mi>j</mml:mi>
                                  <mml:mo>=</mml:mo>
                                  <mml:mn>0</mml:mn>
                                </mml:mrow>
                                <mml:mrow>
                                  <mml:mi>N</mml:mi>
                                  <mml:mo>−</mml:mo>
                                  <mml:mn>1</mml:mn>
                                </mml:mrow>
                              </mml:munderover>
                            </mml:mstyle>
                            <mml:mrow>
                              <mml:mrow>
                                <mml:mstyle displaystyle="true">
                                  <mml:mfrac>
                                    <mml:mrow>
                                      <mml:mn>2</mml:mn>
                                      <mml:mo>⁢</mml:mo>
                                      <mml:mi>a</mml:mi>
                                      <mml:mo>⁢</mml:mo>
                                      <mml:mrow>
                                        <mml:mo stretchy="false">(</mml:mo>
                                        <mml:mi>k</mml:mi>
                                        <mml:mo stretchy="false">)</mml:mo>
                                      </mml:mrow>
                                      <mml:mo>⁢</mml:mo>
                                      <mml:mi>a</mml:mi>
                                      <mml:mo>⁢</mml:mo>
                                      <mml:mrow>
                                        <mml:mo stretchy="false">(</mml:mo>
                                        <mml:mi>l</mml:mi>
                                        <mml:mo stretchy="false">)</mml:mo>
                                      </mml:mrow>
                                    </mml:mrow>
                                    <mml:mi>N</mml:mi>
                                  </mml:mfrac>
                                </mml:mstyle>
                                <mml:mo lspace="0.222em" rspace="0.222em">×</mml:mo>
                                <mml:mi>d</mml:mi>
                              </mml:mrow>
                              <mml:mo>⁢</mml:mo>
                              <mml:mrow>
                                <mml:mo stretchy="false">(</mml:mo>
                                <mml:mi>k</mml:mi>
                                <mml:mo>,</mml:mo>
                                <mml:mi>l</mml:mi>
                                <mml:mo stretchy="false">)</mml:mo>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:mrow>
                      </mml:mrow>
                    </mml:math>
                  </disp-formula>
                  <disp-formula id="S2.E3Xa">
                    <mml:math alttext="\displaystyle\times\cos\left[\frac{(2m+1)\pi k}{2n}\right]" display="inline">
                      <mml:mrow>
                        <mml:mi/>
                        <mml:mo lspace="0.222em" rspace="0.222em">×</mml:mo>
                        <mml:mrow>
                          <mml:mi>cos</mml:mi>
                          <mml:mo>⁡</mml:mo>
                          <mml:mrow>
                            <mml:mo>[</mml:mo>
                            <mml:mstyle displaystyle="true">
                              <mml:mfrac>
                                <mml:mrow>
                                  <mml:mrow>
                                    <mml:mo stretchy="false">(</mml:mo>
                                    <mml:mrow>
                                      <mml:mrow>
                                        <mml:mn>2</mml:mn>
                                        <mml:mo>⁢</mml:mo>
                                        <mml:mi>m</mml:mi>
                                      </mml:mrow>
                                      <mml:mo>+</mml:mo>
                                      <mml:mn>1</mml:mn>
                                    </mml:mrow>
                                    <mml:mo stretchy="false">)</mml:mo>
                                  </mml:mrow>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mi>π</mml:mi>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mi>k</mml:mi>
                                </mml:mrow>
                                <mml:mrow>
                                  <mml:mn>2</mml:mn>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mi>n</mml:mi>
                                </mml:mrow>
                              </mml:mfrac>
                            </mml:mstyle>
                            <mml:mo>]</mml:mo>
                          </mml:mrow>
                        </mml:mrow>
                      </mml:mrow>
                    </mml:math>
                  </disp-formula>
                  <disp-formula id="S2.E3Xb">
                    <mml:math alttext="\displaystyle\times\cos\left[\frac{(2n+1)\pi k}{2n}\right]" display="inline">
                      <mml:mrow>
                        <mml:mi/>
                        <mml:mo lspace="0.222em" rspace="0.222em">×</mml:mo>
                        <mml:mrow>
                          <mml:mi>cos</mml:mi>
                          <mml:mo>⁡</mml:mo>
                          <mml:mrow>
                            <mml:mo>[</mml:mo>
                            <mml:mstyle displaystyle="true">
                              <mml:mfrac>
                                <mml:mrow>
                                  <mml:mrow>
                                    <mml:mo stretchy="false">(</mml:mo>
                                    <mml:mrow>
                                      <mml:mrow>
                                        <mml:mn>2</mml:mn>
                                        <mml:mo>⁢</mml:mo>
                                        <mml:mi>n</mml:mi>
                                      </mml:mrow>
                                      <mml:mo>+</mml:mo>
                                      <mml:mn>1</mml:mn>
                                    </mml:mrow>
                                    <mml:mo stretchy="false">)</mml:mo>
                                  </mml:mrow>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mi>π</mml:mi>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mi>k</mml:mi>
                                </mml:mrow>
                                <mml:mrow>
                                  <mml:mn>2</mml:mn>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mi>n</mml:mi>
                                </mml:mrow>
                              </mml:mfrac>
                            </mml:mstyle>
                            <mml:mo>]</mml:mo>
                          </mml:mrow>
                        </mml:mrow>
                      </mml:mrow>
                    </mml:math>
                  </disp-formula>
                </disp-formula-group>
              </p>
              <p id="S2.SS1.p5">where <inline-formula><mml:math alttext="m,n=0,1,\ldots,N-1" display="inline"><mml:mrow><mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mo>,</mml:mo><mml:mi>n</mml:mi></mml:mrow><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mo>,</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mi mathvariant="normal">…</mml:mi><mml:mo>,</mml:mo><mml:mrow><mml:mi>N</mml:mi><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:mrow></mml:mrow></mml:math></inline-formula>.</p>
              <p id="S2.SS1.p6">In Eq.(1), <inline-formula><mml:math alttext="d(0,0)" display="inline"><mml:mrow><mml:mi>d</mml:mi><mml:mo>⁢</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>0</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> is the DC coefficient, which is the coefficient with zero frequency in both dimensions, and the other <inline-formula><mml:math alttext="d(k,l)" display="inline"><mml:mrow><mml:mi>d</mml:mi><mml:mo>⁢</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>l</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula>s are the AC coefficients, which are the remaining coefficients with non-zero frequency of the block.</p>
            </sec>
            <sec id="S2.SS2">
              <label>2.2</label>
              <title>Singular Value Decomposition (SVD) in DCT Domain</title>
              <p id="S2.SS2.p1">SVD is a mathematical technique used to factorize a matrix into three component matrices. In the context of the DCT domain for image fusion, SVD helps in combining the important features of multiple images into a single, more informative image [<xref rid="ref013" ref-type="bibr">13</xref>]. Decompose the DCT-transformed images Eq.(1) using SVD.</p>
              <p>
                <disp-formula id="S2.E4">
                  <mml:math alttext="A=U\Sigma V^{T}" display="block">
                    <mml:mrow>
                      <mml:mi>A</mml:mi>
                      <mml:mo>=</mml:mo>
                      <mml:mrow>
                        <mml:mi>U</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:mi mathvariant="normal">Σ</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:msup>
                          <mml:mi>V</mml:mi>
                          <mml:mi>T</mml:mi>
                        </mml:msup>
                      </mml:mrow>
                    </mml:mrow>
                  </mml:math>
                </disp-formula>
              </p>
              <p id="S2.SS2.p3">where <inline-formula><mml:math alttext="A" display="inline"><mml:mi>A</mml:mi></mml:math></inline-formula> is the DCT-transformed matrix, <inline-formula><mml:math alttext="U" display="inline"><mml:mi>U</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math alttext="V" display="inline"><mml:mi>V</mml:mi></mml:math></inline-formula> are orthogonal matrices, and <inline-formula><mml:math alttext="\Sigma" display="inline"><mml:mi mathvariant="normal">Σ</mml:mi></mml:math></inline-formula> is a diagonal matrix containing the singular values. Combine the singular values from the corresponding matrices of the input images using a fusion rule, such as averaging or maximum selection.</p>
              <p>
                <disp-formula id="S2.E5">
                  <mml:math alttext="\Sigma_{f}=\frac{\Sigma_{1}+\Sigma_{2}}{2}" display="block">
                    <mml:mrow>
                      <mml:msub>
                        <mml:mi mathvariant="normal">Σ</mml:mi>
                        <mml:mi>f</mml:mi>
                      </mml:msub>
                      <mml:mo>=</mml:mo>
                      <mml:mfrac>
                        <mml:mrow>
                          <mml:msub>
                            <mml:mi mathvariant="normal">Σ</mml:mi>
                            <mml:mn>1</mml:mn>
                          </mml:msub>
                          <mml:mo>+</mml:mo>
                          <mml:msub>
                            <mml:mi mathvariant="normal">Σ</mml:mi>
                            <mml:mn>2</mml:mn>
                          </mml:msub>
                        </mml:mrow>
                        <mml:mn>2</mml:mn>
                      </mml:mfrac>
                    </mml:mrow>
                  </mml:math>
                </disp-formula>
              </p>
              <p id="S2.SS2.p5">where <inline-formula><mml:math alttext="\Sigma_{f}" display="inline"><mml:msub><mml:mi mathvariant="normal">Σ</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:math></inline-formula> is the resultant singular value matrix, <inline-formula><mml:math alttext="\Sigma_{1}" display="inline"><mml:msub><mml:mi mathvariant="normal">Σ</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:math></inline-formula> and <inline-formula><mml:math alttext="\Sigma_{2}" display="inline"><mml:msub><mml:mi mathvariant="normal">Σ</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:math></inline-formula> are the singular value matrices of the reference image. Apply inverse SVD to the resultant matrices and then use IDCT to transform using Eq.(3) back to the spatial domain [<xref rid="ref014" ref-type="bibr">14</xref>].</p>
              <p>
                <disp-formula id="S2.E6">
                  <mml:math alttext="A_{f}=U\Sigma_{f}V^{T}" display="block">
                    <mml:mrow>
                      <mml:msub>
                        <mml:mi>A</mml:mi>
                        <mml:mi>f</mml:mi>
                      </mml:msub>
                      <mml:mo>=</mml:mo>
                      <mml:mrow>
                        <mml:mi>U</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:msub>
                          <mml:mi mathvariant="normal">Σ</mml:mi>
                          <mml:mi>f</mml:mi>
                        </mml:msub>
                        <mml:mo>⁢</mml:mo>
                        <mml:msup>
                          <mml:mi>V</mml:mi>
                          <mml:mi>T</mml:mi>
                        </mml:msup>
                      </mml:mrow>
                    </mml:mrow>
                  </mml:math>
                </disp-formula>
              </p>
            </sec>
            <sec id="S2.SS3">
              <label>2.3</label>
              <title>Singular Value Decomposition (SVD) + Consistency Verification (CV) in DCT domain</title>
              <p id="S2.SS3.p1">SVD + CV are applied in the DCT domain for multi-focus image fusion. This method enhances the fusion process by ensuring that the merged image retains the highly relevant information from both input images. Using SVD in the DCT domain helps in effectively capturing the frequency information, while CV ensures that the resultant image is consistent and free from inconsistencies or artifacts. Ensure that the fused singular values using Eq.(4) are consistent with the input images using Eq.(3) by verifying the consistency of the coefficients:</p>
              <p>
                <disp-formula id="S2.E7">
                  <mml:math alttext="CV(i,j)=\begin{cases}1&amp;\text{if }|F_{1}(i,j)-F_{2}(i,j)|&lt;T\\&#10;0&amp;\text{otherwise}\end{cases}" display="block">
                    <mml:mrow>
                      <mml:mrow>
                        <mml:mi>C</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:mi>V</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:mrow>
                          <mml:mo stretchy="false">(</mml:mo>
                          <mml:mi>i</mml:mi>
                          <mml:mo>,</mml:mo>
                          <mml:mi>j</mml:mi>
                          <mml:mo stretchy="false">)</mml:mo>
                        </mml:mrow>
                      </mml:mrow>
                      <mml:mo>=</mml:mo>
                      <mml:mrow>
                        <mml:mo>{</mml:mo>
                        <mml:mtable columnspacing="5pt" displaystyle="true" rowspacing="0pt">
                          <mml:mtr>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mn>1</mml:mn>
                            </mml:mtd>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mrow>
                                <mml:mrow>
                                  <mml:mtext>if </mml:mtext>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mrow>
                                    <mml:mo stretchy="false">|</mml:mo>
                                    <mml:mrow>
                                      <mml:mrow>
                                        <mml:msub>
                                          <mml:mi>F</mml:mi>
                                          <mml:mn>1</mml:mn>
                                        </mml:msub>
                                        <mml:mo>⁢</mml:mo>
                                        <mml:mrow>
                                          <mml:mo stretchy="false">(</mml:mo>
                                          <mml:mi>i</mml:mi>
                                          <mml:mo>,</mml:mo>
                                          <mml:mi>j</mml:mi>
                                          <mml:mo stretchy="false">)</mml:mo>
                                        </mml:mrow>
                                      </mml:mrow>
                                      <mml:mo>−</mml:mo>
                                      <mml:mrow>
                                        <mml:msub>
                                          <mml:mi>F</mml:mi>
                                          <mml:mn>2</mml:mn>
                                        </mml:msub>
                                        <mml:mo>⁢</mml:mo>
                                        <mml:mrow>
                                          <mml:mo stretchy="false">(</mml:mo>
                                          <mml:mi>i</mml:mi>
                                          <mml:mo>,</mml:mo>
                                          <mml:mi>j</mml:mi>
                                          <mml:mo stretchy="false">)</mml:mo>
                                        </mml:mrow>
                                      </mml:mrow>
                                    </mml:mrow>
                                    <mml:mo stretchy="false">|</mml:mo>
                                  </mml:mrow>
                                </mml:mrow>
                                <mml:mo>&lt;</mml:mo>
                                <mml:mi>T</mml:mi>
                              </mml:mrow>
                            </mml:mtd>
                          </mml:mtr>
                          <mml:mtr>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mn>0</mml:mn>
                            </mml:mtd>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mtext>otherwise</mml:mtext>
                            </mml:mtd>
                          </mml:mtr>
                        </mml:mtable>
                      </mml:mrow>
                    </mml:mrow>
                  </mml:math>
                </disp-formula>
              </p>
              <p id="S2.SS3.p2">where, <inline-formula><mml:math alttext="T" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> is a threshold, <inline-formula><mml:math alttext="F_{1}(i,j)" display="inline"><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>⁢</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math alttext="F_{2}(i,j)" display="inline"><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>⁢</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> are the DCT coefficients of the input images. If the difference is within the threshold, the value is kept; otherwise, it is discarded [<xref rid="ref013" ref-type="bibr">13</xref>].</p>
            </sec>
            <sec id="S2.SS4">
              <label>2.4</label>
              <title>Correlation Coefficient and Energy-Correlation in the DCT domain</title>
              <p id="S2.SS4.p1">The Correlation Coefficient and Energy-Correlation Coefficient are computed in the DCT domain for an image fusion scenario. These metrics are essential in evaluating the quality of the resultant image, particularly in maintaining the correlation and energy consistency between the original and resultant images. Correlation coefficients are used to measure how effectively the fused image retains the information from the underlying photos, ensuring the quality of the resulting image. The Energy-Correlation Coefficient specifically measures how well the energy of the image (related to its contrast and brightness) is preserved during the fusion process.</p>
              <p id="S2.SS4.p2">The Correlation Coefficient between the original image <inline-formula><mml:math alttext="A" display="inline"><mml:mi>A</mml:mi></mml:math></inline-formula> and the fused image <inline-formula><mml:math alttext="F" display="inline"><mml:mi>F</mml:mi></mml:math></inline-formula> is computing using:</p>
              <p>
                <disp-formula id="S2.E8">
                  <mml:math alttext="CC=\frac{\sum\limits_{i=1}^{N}\sum\limits_{j=1}^{M}(A(i,j)-\overline{A})(F(i,j%&#10;)-\overline{F})}{\sqrt{\sum\limits_{i=1}^{N}\sum\limits_{j=1}^{M}(A(i,j)-%&#10;\overline{A})^{2}}\sqrt{\sum\limits_{i=1}^{N}\sum\limits_{j=1}^{M}(F(i,j)-%&#10;\overline{F})^{2}}}" display="block">
                    <mml:mrow>
                      <mml:mrow>
                        <mml:mi>C</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:mi>C</mml:mi>
                      </mml:mrow>
                      <mml:mo>=</mml:mo>
                      <mml:mfrac>
                        <mml:mrow>
                          <mml:munderover>
                            <mml:mo movablelimits="false">∑</mml:mo>
                            <mml:mrow>
                              <mml:mi>i</mml:mi>
                              <mml:mo>=</mml:mo>
                              <mml:mn>1</mml:mn>
                            </mml:mrow>
                            <mml:mi>N</mml:mi>
                          </mml:munderover>
                          <mml:mrow>
                            <mml:munderover>
                              <mml:mo lspace="0.167em" movablelimits="false" rspace="0em">∑</mml:mo>
                              <mml:mrow>
                                <mml:mi>j</mml:mi>
                                <mml:mo>=</mml:mo>
                                <mml:mn>1</mml:mn>
                              </mml:mrow>
                              <mml:mi>M</mml:mi>
                            </mml:munderover>
                            <mml:mrow>
                              <mml:mrow>
                                <mml:mo stretchy="false">(</mml:mo>
                                <mml:mrow>
                                  <mml:mrow>
                                    <mml:mi>A</mml:mi>
                                    <mml:mo>⁢</mml:mo>
                                    <mml:mrow>
                                      <mml:mo stretchy="false">(</mml:mo>
                                      <mml:mi>i</mml:mi>
                                      <mml:mo>,</mml:mo>
                                      <mml:mi>j</mml:mi>
                                      <mml:mo stretchy="false">)</mml:mo>
                                    </mml:mrow>
                                  </mml:mrow>
                                  <mml:mo>−</mml:mo>
                                  <mml:mover accent="true">
                                    <mml:mi>A</mml:mi>
                                    <mml:mo>¯</mml:mo>
                                  </mml:mover>
                                </mml:mrow>
                                <mml:mo stretchy="false">)</mml:mo>
                              </mml:mrow>
                              <mml:mo>⁢</mml:mo>
                              <mml:mrow>
                                <mml:mo stretchy="false">(</mml:mo>
                                <mml:mrow>
                                  <mml:mrow>
                                    <mml:mi>F</mml:mi>
                                    <mml:mo>⁢</mml:mo>
                                    <mml:mrow>
                                      <mml:mo stretchy="false">(</mml:mo>
                                      <mml:mi>i</mml:mi>
                                      <mml:mo>,</mml:mo>
                                      <mml:mi>j</mml:mi>
                                      <mml:mo stretchy="false">)</mml:mo>
                                    </mml:mrow>
                                  </mml:mrow>
                                  <mml:mo>−</mml:mo>
                                  <mml:mover accent="true">
                                    <mml:mi>F</mml:mi>
                                    <mml:mo>¯</mml:mo>
                                  </mml:mover>
                                </mml:mrow>
                                <mml:mo stretchy="false">)</mml:mo>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:mrow>
                        <mml:mrow>
                          <mml:msqrt>
                            <mml:mrow>
                              <mml:munderover>
                                <mml:mo movablelimits="false">∑</mml:mo>
                                <mml:mrow>
                                  <mml:mi>i</mml:mi>
                                  <mml:mo>=</mml:mo>
                                  <mml:mn>1</mml:mn>
                                </mml:mrow>
                                <mml:mi>N</mml:mi>
                              </mml:munderover>
                              <mml:mrow>
                                <mml:munderover>
                                  <mml:mo lspace="0.167em" movablelimits="false" rspace="0em">∑</mml:mo>
                                  <mml:mrow>
                                    <mml:mi>j</mml:mi>
                                    <mml:mo>=</mml:mo>
                                    <mml:mn>1</mml:mn>
                                  </mml:mrow>
                                  <mml:mi>M</mml:mi>
                                </mml:munderover>
                                <mml:msup>
                                  <mml:mrow>
                                    <mml:mo stretchy="false">(</mml:mo>
                                    <mml:mrow>
                                      <mml:mrow>
                                        <mml:mi>A</mml:mi>
                                        <mml:mo>⁢</mml:mo>
                                        <mml:mrow>
                                          <mml:mo stretchy="false">(</mml:mo>
                                          <mml:mi>i</mml:mi>
                                          <mml:mo>,</mml:mo>
                                          <mml:mi>j</mml:mi>
                                          <mml:mo stretchy="false">)</mml:mo>
                                        </mml:mrow>
                                      </mml:mrow>
                                      <mml:mo>−</mml:mo>
                                      <mml:mover accent="true">
                                        <mml:mi>A</mml:mi>
                                        <mml:mo>¯</mml:mo>
                                      </mml:mover>
                                    </mml:mrow>
                                    <mml:mo stretchy="false">)</mml:mo>
                                  </mml:mrow>
                                  <mml:mn>2</mml:mn>
                                </mml:msup>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:msqrt>
                          <mml:mo>⁢</mml:mo>
                          <mml:msqrt>
                            <mml:mrow>
                              <mml:munderover>
                                <mml:mo movablelimits="false">∑</mml:mo>
                                <mml:mrow>
                                  <mml:mi>i</mml:mi>
                                  <mml:mo>=</mml:mo>
                                  <mml:mn>1</mml:mn>
                                </mml:mrow>
                                <mml:mi>N</mml:mi>
                              </mml:munderover>
                              <mml:mrow>
                                <mml:munderover>
                                  <mml:mo lspace="0.167em" movablelimits="false" rspace="0em">∑</mml:mo>
                                  <mml:mrow>
                                    <mml:mi>j</mml:mi>
                                    <mml:mo>=</mml:mo>
                                    <mml:mn>1</mml:mn>
                                  </mml:mrow>
                                  <mml:mi>M</mml:mi>
                                </mml:munderover>
                                <mml:msup>
                                  <mml:mrow>
                                    <mml:mo stretchy="false">(</mml:mo>
                                    <mml:mrow>
                                      <mml:mrow>
                                        <mml:mi>F</mml:mi>
                                        <mml:mo>⁢</mml:mo>
                                        <mml:mrow>
                                          <mml:mo stretchy="false">(</mml:mo>
                                          <mml:mi>i</mml:mi>
                                          <mml:mo>,</mml:mo>
                                          <mml:mi>j</mml:mi>
                                          <mml:mo stretchy="false">)</mml:mo>
                                        </mml:mrow>
                                      </mml:mrow>
                                      <mml:mo>−</mml:mo>
                                      <mml:mover accent="true">
                                        <mml:mi>F</mml:mi>
                                        <mml:mo>¯</mml:mo>
                                      </mml:mover>
                                    </mml:mrow>
                                    <mml:mo stretchy="false">)</mml:mo>
                                  </mml:mrow>
                                  <mml:mn>2</mml:mn>
                                </mml:msup>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:msqrt>
                        </mml:mrow>
                      </mml:mfrac>
                    </mml:mrow>
                  </mml:math>
                </disp-formula>
              </p>
              <p id="S2.SS4.p3">where <inline-formula><mml:math alttext="\overline{A}" display="inline"><mml:mover accent="true"><mml:mi>A</mml:mi><mml:mo>¯</mml:mo></mml:mover></mml:math></inline-formula> and <inline-formula><mml:math alttext="\overline{F}" display="inline"><mml:mover accent="true"><mml:mi>F</mml:mi><mml:mo>¯</mml:mo></mml:mover></mml:math></inline-formula> are the mean values of the original and resultant images, respectively.</p>
              <p id="S2.SS4.p4">The energy of an image in the DCT domain can be represented as:</p>
              <p>
                <disp-formula id="S2.E9">
                  <mml:math alttext="E_{A}=\sum_{u=0}^{N-1}\sum_{v=0}^{M-1}|F_{A}(u,v)|^{2}" display="block">
                    <mml:mrow>
                      <mml:msub>
                        <mml:mi>E</mml:mi>
                        <mml:mi>A</mml:mi>
                      </mml:msub>
                      <mml:mo rspace="0.111em">=</mml:mo>
                      <mml:mrow>
                        <mml:munderover>
                          <mml:mo movablelimits="false" rspace="0em">∑</mml:mo>
                          <mml:mrow>
                            <mml:mi>u</mml:mi>
                            <mml:mo>=</mml:mo>
                            <mml:mn>0</mml:mn>
                          </mml:mrow>
                          <mml:mrow>
                            <mml:mi>N</mml:mi>
                            <mml:mo>−</mml:mo>
                            <mml:mn>1</mml:mn>
                          </mml:mrow>
                        </mml:munderover>
                        <mml:mrow>
                          <mml:munderover>
                            <mml:mo movablelimits="false" rspace="0em">∑</mml:mo>
                            <mml:mrow>
                              <mml:mi>v</mml:mi>
                              <mml:mo>=</mml:mo>
                              <mml:mn>0</mml:mn>
                            </mml:mrow>
                            <mml:mrow>
                              <mml:mi>M</mml:mi>
                              <mml:mo>−</mml:mo>
                              <mml:mn>1</mml:mn>
                            </mml:mrow>
                          </mml:munderover>
                          <mml:msup>
                            <mml:mrow>
                              <mml:mo stretchy="false">|</mml:mo>
                              <mml:mrow>
                                <mml:msub>
                                  <mml:mi>F</mml:mi>
                                  <mml:mi>A</mml:mi>
                                </mml:msub>
                                <mml:mo>⁢</mml:mo>
                                <mml:mrow>
                                  <mml:mo stretchy="false">(</mml:mo>
                                  <mml:mi>u</mml:mi>
                                  <mml:mo>,</mml:mo>
                                  <mml:mi>v</mml:mi>
                                  <mml:mo stretchy="false">)</mml:mo>
                                </mml:mrow>
                              </mml:mrow>
                              <mml:mo stretchy="false">|</mml:mo>
                            </mml:mrow>
                            <mml:mn>2</mml:mn>
                          </mml:msup>
                        </mml:mrow>
                      </mml:mrow>
                    </mml:mrow>
                  </mml:math>
                </disp-formula>
              </p>
              <p id="S2.SS4.p5">where <inline-formula><mml:math alttext="E_{A}" display="inline"><mml:msub><mml:mi>E</mml:mi><mml:mi>A</mml:mi></mml:msub></mml:math></inline-formula> is the energy of the image <inline-formula><mml:math alttext="A" display="inline"><mml:mi>A</mml:mi></mml:math></inline-formula>, and <inline-formula><mml:math alttext="F_{A}(u,v)" display="inline"><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:mo>⁢</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>u</mml:mi><mml:mo>,</mml:mo><mml:mi>v</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> are the DCT coefficients.</p>
              <p id="S2.SS4.p6">The Energy-Correlation Coefficient among the original image <inline-formula><mml:math alttext="A" display="inline"><mml:mi>A</mml:mi></mml:math></inline-formula> and the resultant image <inline-formula><mml:math alttext="F" display="inline"><mml:mi>F</mml:mi></mml:math></inline-formula> is computing using:</p>
              <p>
                <disp-formula id="S2.E10">
                  <mml:math alttext="ECC=\frac{\sum\limits_{u=0}^{N-1}\sum\limits_{\nu=0}^{M-1}(|F_{A}(u,\nu)|\cdot%&#10;|F_{F}(u,\nu)|)}{\sqrt{\sum\limits_{u=0}^{N-1}\sum\limits_{\nu=0}^{M-1}|F_{A}(%&#10;u,\nu)|^{2}}\sqrt{\sum\limits_{u=0}^{N-1}\sum\limits_{\nu=0}^{M-1}|F_{F}(u,\nu%&#10;)|^{2}}}" display="inline">
                    <mml:mrow>
                      <mml:mrow>
                        <mml:mi>E</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:mi>C</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:mi>C</mml:mi>
                      </mml:mrow>
                      <mml:mo>=</mml:mo>
                      <mml:mfrac>
                        <mml:mrow>
                          <mml:mstyle displaystyle="false">
                            <mml:munderover>
                              <mml:mo movablelimits="false">∑</mml:mo>
                              <mml:mrow>
                                <mml:mi>u</mml:mi>
                                <mml:mo>=</mml:mo>
                                <mml:mn>0</mml:mn>
                              </mml:mrow>
                              <mml:mrow>
                                <mml:mi>N</mml:mi>
                                <mml:mo>−</mml:mo>
                                <mml:mn>1</mml:mn>
                              </mml:mrow>
                            </mml:munderover>
                          </mml:mstyle>
                          <mml:mrow>
                            <mml:mstyle displaystyle="false">
                              <mml:munderover>
                                <mml:mo movablelimits="false">∑</mml:mo>
                                <mml:mrow>
                                  <mml:mi>ν</mml:mi>
                                  <mml:mo>=</mml:mo>
                                  <mml:mn>0</mml:mn>
                                </mml:mrow>
                                <mml:mrow>
                                  <mml:mi>M</mml:mi>
                                  <mml:mo>−</mml:mo>
                                  <mml:mn>1</mml:mn>
                                </mml:mrow>
                              </mml:munderover>
                            </mml:mstyle>
                            <mml:mrow>
                              <mml:mo stretchy="false">(</mml:mo>
                              <mml:mrow>
                                <mml:mrow>
                                  <mml:mo stretchy="false">|</mml:mo>
                                  <mml:mrow>
                                    <mml:msub>
                                      <mml:mi>F</mml:mi>
                                      <mml:mi>A</mml:mi>
                                    </mml:msub>
                                    <mml:mo>⁢</mml:mo>
                                    <mml:mrow>
                                      <mml:mo stretchy="false">(</mml:mo>
                                      <mml:mi>u</mml:mi>
                                      <mml:mo>,</mml:mo>
                                      <mml:mi>ν</mml:mi>
                                      <mml:mo stretchy="false">)</mml:mo>
                                    </mml:mrow>
                                  </mml:mrow>
                                  <mml:mo rspace="0.055em" stretchy="false">|</mml:mo>
                                </mml:mrow>
                                <mml:mo rspace="0.222em">⋅</mml:mo>
                                <mml:mrow>
                                  <mml:mo stretchy="false">|</mml:mo>
                                  <mml:mrow>
                                    <mml:msub>
                                      <mml:mi>F</mml:mi>
                                      <mml:mi>F</mml:mi>
                                    </mml:msub>
                                    <mml:mo>⁢</mml:mo>
                                    <mml:mrow>
                                      <mml:mo stretchy="false">(</mml:mo>
                                      <mml:mi>u</mml:mi>
                                      <mml:mo>,</mml:mo>
                                      <mml:mi>ν</mml:mi>
                                      <mml:mo stretchy="false">)</mml:mo>
                                    </mml:mrow>
                                  </mml:mrow>
                                  <mml:mo stretchy="false">|</mml:mo>
                                </mml:mrow>
                              </mml:mrow>
                              <mml:mo stretchy="false">)</mml:mo>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:mrow>
                        <mml:mrow>
                          <mml:msqrt>
                            <mml:mrow>
                              <mml:mstyle displaystyle="false">
                                <mml:munderover>
                                  <mml:mo movablelimits="false">∑</mml:mo>
                                  <mml:mrow>
                                    <mml:mi>u</mml:mi>
                                    <mml:mo>=</mml:mo>
                                    <mml:mn>0</mml:mn>
                                  </mml:mrow>
                                  <mml:mrow>
                                    <mml:mi>N</mml:mi>
                                    <mml:mo>−</mml:mo>
                                    <mml:mn>1</mml:mn>
                                  </mml:mrow>
                                </mml:munderover>
                              </mml:mstyle>
                              <mml:mrow>
                                <mml:mstyle displaystyle="false">
                                  <mml:munderover>
                                    <mml:mo movablelimits="false">∑</mml:mo>
                                    <mml:mrow>
                                      <mml:mi>ν</mml:mi>
                                      <mml:mo>=</mml:mo>
                                      <mml:mn>0</mml:mn>
                                    </mml:mrow>
                                    <mml:mrow>
                                      <mml:mi>M</mml:mi>
                                      <mml:mo>−</mml:mo>
                                      <mml:mn>1</mml:mn>
                                    </mml:mrow>
                                  </mml:munderover>
                                </mml:mstyle>
                                <mml:msup>
                                  <mml:mrow>
                                    <mml:mo stretchy="false">|</mml:mo>
                                    <mml:mrow>
                                      <mml:msub>
                                        <mml:mi>F</mml:mi>
                                        <mml:mi>A</mml:mi>
                                      </mml:msub>
                                      <mml:mo>⁢</mml:mo>
                                      <mml:mrow>
                                        <mml:mo stretchy="false">(</mml:mo>
                                        <mml:mi>u</mml:mi>
                                        <mml:mo>,</mml:mo>
                                        <mml:mi>ν</mml:mi>
                                        <mml:mo stretchy="false">)</mml:mo>
                                      </mml:mrow>
                                    </mml:mrow>
                                    <mml:mo stretchy="false">|</mml:mo>
                                  </mml:mrow>
                                  <mml:mn>2</mml:mn>
                                </mml:msup>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:msqrt>
                          <mml:mo>⁢</mml:mo>
                          <mml:msqrt>
                            <mml:mrow>
                              <mml:mstyle displaystyle="false">
                                <mml:munderover>
                                  <mml:mo movablelimits="false">∑</mml:mo>
                                  <mml:mrow>
                                    <mml:mi>u</mml:mi>
                                    <mml:mo>=</mml:mo>
                                    <mml:mn>0</mml:mn>
                                  </mml:mrow>
                                  <mml:mrow>
                                    <mml:mi>N</mml:mi>
                                    <mml:mo>−</mml:mo>
                                    <mml:mn>1</mml:mn>
                                  </mml:mrow>
                                </mml:munderover>
                              </mml:mstyle>
                              <mml:mrow>
                                <mml:mstyle displaystyle="false">
                                  <mml:munderover>
                                    <mml:mo movablelimits="false">∑</mml:mo>
                                    <mml:mrow>
                                      <mml:mi>ν</mml:mi>
                                      <mml:mo>=</mml:mo>
                                      <mml:mn>0</mml:mn>
                                    </mml:mrow>
                                    <mml:mrow>
                                      <mml:mi>M</mml:mi>
                                      <mml:mo>−</mml:mo>
                                      <mml:mn>1</mml:mn>
                                    </mml:mrow>
                                  </mml:munderover>
                                </mml:mstyle>
                                <mml:msup>
                                  <mml:mrow>
                                    <mml:mo stretchy="false">|</mml:mo>
                                    <mml:mrow>
                                      <mml:msub>
                                        <mml:mi>F</mml:mi>
                                        <mml:mi>F</mml:mi>
                                      </mml:msub>
                                      <mml:mo>⁢</mml:mo>
                                      <mml:mrow>
                                        <mml:mo stretchy="false">(</mml:mo>
                                        <mml:mi>u</mml:mi>
                                        <mml:mo>,</mml:mo>
                                        <mml:mi>ν</mml:mi>
                                        <mml:mo stretchy="false">)</mml:mo>
                                      </mml:mrow>
                                    </mml:mrow>
                                    <mml:mo stretchy="false">|</mml:mo>
                                  </mml:mrow>
                                  <mml:mn>2</mml:mn>
                                </mml:msup>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:msqrt>
                        </mml:mrow>
                      </mml:mfrac>
                    </mml:mrow>
                  </mml:math>
                </disp-formula>
              </p>
              <p id="S2.SS4.p7">where <inline-formula><mml:math alttext="F_{A}(u,v)" display="inline"><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>A</mml:mi></mml:msub><mml:mo>⁢</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>u</mml:mi><mml:mo>,</mml:mo><mml:mi>v</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math alttext="F_{F}(u,v)" display="inline"><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi>F</mml:mi></mml:msub><mml:mo>⁢</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>u</mml:mi><mml:mo>,</mml:mo><mml:mi>v</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> are the DCT coefficients of the original and resultant images, respectively.</p>
            </sec>
            <sec id="S2.SS5">
              <label>2.5</label>
              <title>Correlation Coefficient and Energy-Correlation + Consistency Verification in DCT domain</title>
              <p id="S2.SS5.p1">For multi-focus image fusion in the DCT domain, the Correlation Coefficient, Energy-Correlation Coefficient, and Consistency Verification are employed. These techniques work together to ensure that the fused image maintains a high level of quality, retaining important information and energy consistency from the original images while ensuring that the fusion process produces reliable and accurate results. The DCT domain is particularly useful in this context because it allows the image’s frequency components to be manipulated, making it easier to merge images effectively. In order to ensure that crucial information is retained, the Correlation Coefficient calculates how well the fused image corresponds with the original images. The Energy-Correlation Coefficient evaluates how well the image’s energy (related to its contrast and brightness) is retained. Finally, consistency verification makes sure that the integrity of the fused image is maintained by preventing artifacts or inconsistencies from being introduced throughout the fusion process.</p>
              <p>
                <disp-formula id="S2.E11">
                  <mml:math alttext="CV(i,j)=\begin{cases}1&amp;\text{if }|F_{1}(i,j)-F_{2}(i,j)|&lt;T\\&#10;0&amp;\text{otherwise}\end{cases}" display="block">
                    <mml:mrow>
                      <mml:mrow>
                        <mml:mi>C</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:mi>V</mml:mi>
                        <mml:mo>⁢</mml:mo>
                        <mml:mrow>
                          <mml:mo stretchy="false">(</mml:mo>
                          <mml:mi>i</mml:mi>
                          <mml:mo>,</mml:mo>
                          <mml:mi>j</mml:mi>
                          <mml:mo stretchy="false">)</mml:mo>
                        </mml:mrow>
                      </mml:mrow>
                      <mml:mo>=</mml:mo>
                      <mml:mrow>
                        <mml:mo>{</mml:mo>
                        <mml:mtable columnspacing="5pt" displaystyle="true" rowspacing="0pt">
                          <mml:mtr>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mn>1</mml:mn>
                            </mml:mtd>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mrow>
                                <mml:mrow>
                                  <mml:mtext>if </mml:mtext>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mrow>
                                    <mml:mo stretchy="false">|</mml:mo>
                                    <mml:mrow>
                                      <mml:mrow>
                                        <mml:msub>
                                          <mml:mi>F</mml:mi>
                                          <mml:mn>1</mml:mn>
                                        </mml:msub>
                                        <mml:mo>⁢</mml:mo>
                                        <mml:mrow>
                                          <mml:mo stretchy="false">(</mml:mo>
                                          <mml:mi>i</mml:mi>
                                          <mml:mo>,</mml:mo>
                                          <mml:mi>j</mml:mi>
                                          <mml:mo stretchy="false">)</mml:mo>
                                        </mml:mrow>
                                      </mml:mrow>
                                      <mml:mo>−</mml:mo>
                                      <mml:mrow>
                                        <mml:msub>
                                          <mml:mi>F</mml:mi>
                                          <mml:mn>2</mml:mn>
                                        </mml:msub>
                                        <mml:mo>⁢</mml:mo>
                                        <mml:mrow>
                                          <mml:mo stretchy="false">(</mml:mo>
                                          <mml:mi>i</mml:mi>
                                          <mml:mo>,</mml:mo>
                                          <mml:mi>j</mml:mi>
                                          <mml:mo stretchy="false">)</mml:mo>
                                        </mml:mrow>
                                      </mml:mrow>
                                    </mml:mrow>
                                    <mml:mo stretchy="false">|</mml:mo>
                                  </mml:mrow>
                                </mml:mrow>
                                <mml:mo>&lt;</mml:mo>
                                <mml:mi>T</mml:mi>
                              </mml:mrow>
                            </mml:mtd>
                          </mml:mtr>
                          <mml:mtr>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mn>0</mml:mn>
                            </mml:mtd>
                            <mml:mtd class="ltx_align_left" columnalign="left">
                              <mml:mtext>otherwise</mml:mtext>
                            </mml:mtd>
                          </mml:mtr>
                        </mml:mtable>
                      </mml:mrow>
                    </mml:mrow>
                  </mml:math>
                </disp-formula>
              </p>
              <p id="S2.SS5.p3">where, <inline-formula><mml:math alttext="T" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> is a threshold, <inline-formula><mml:math alttext="F_{1}(i,j)" display="inline"><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>⁢</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math alttext="F_{2}(i,j)" display="inline"><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>⁢</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> are the DCT coefficients of the input images. If the difference is within the threshold, the value is kept; otherwise, it is discarded [<xref rid="ref013" ref-type="bibr">13</xref>].</p>
              <sec id="S3">
                <label>3.</label>
                <title>Experimentations</title>
                <sec id="S3.SS1">
                  <label>3.1</label>
                  <title>Performance Metrics</title>
                  <p id="S3.SS1.p1"><bold>Entropy</bold> is one of the commonly used and important metrics to evaluate the information content of the resultant image. The higher values mean good results.</p>
                  <p>
                    <disp-formula id="S3.E1">
                      <mml:math alttext="E=-\sum_{k=0}^{G-1}S_{k}\log S_{k}" display="block">
                        <mml:mrow>
                          <mml:mi>E</mml:mi>
                          <mml:mo>=</mml:mo>
                          <mml:mrow>
                            <mml:mo>−</mml:mo>
                            <mml:mrow>
                              <mml:munderover>
                                <mml:mo movablelimits="false">∑</mml:mo>
                                <mml:mrow>
                                  <mml:mi>k</mml:mi>
                                  <mml:mo>=</mml:mo>
                                  <mml:mn>0</mml:mn>
                                </mml:mrow>
                                <mml:mrow>
                                  <mml:mi>G</mml:mi>
                                  <mml:mo>−</mml:mo>
                                  <mml:mn>1</mml:mn>
                                </mml:mrow>
                              </mml:munderover>
                              <mml:mrow>
                                <mml:msub>
                                  <mml:mi>S</mml:mi>
                                  <mml:mi>k</mml:mi>
                                </mml:msub>
                                <mml:mo lspace="0.167em">⁢</mml:mo>
                                <mml:mrow>
                                  <mml:mi>log</mml:mi>
                                  <mml:mo lspace="0.167em">⁡</mml:mo>
                                  <mml:msub>
                                    <mml:mi>S</mml:mi>
                                    <mml:mi>k</mml:mi>
                                  </mml:msub>
                                </mml:mrow>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:mrow>
                      </mml:math>
                    </disp-formula>
                  </p>
                  <p id="S3.SS1.p3"><bold>Correlation Coefficient</bold> (Corr) is an important measure that presents the correlation and also computes the similarity of spectral features between the reference and resultant images. The best value is close to the positive one, which depicts that the reference and resultant images are similar and if the image is dissimilar then the value is closer to zero [<xref rid="ref016" ref-type="bibr">16</xref>].</p>
                  <p>
                    <disp-formula-group id="S3.E2">
                      <disp-formula id="S3.E2X">
                        <mml:math alttext="\displaystyle Corr=\frac{2C_{zp}}{C_{z}+C_{p}}" display="inline">
                          <mml:mrow>
                            <mml:mrow>
                              <mml:mi>C</mml:mi>
                              <mml:mo>⁢</mml:mo>
                              <mml:mi>o</mml:mi>
                              <mml:mo>⁢</mml:mo>
                              <mml:mi>r</mml:mi>
                              <mml:mo>⁢</mml:mo>
                              <mml:mi>r</mml:mi>
                            </mml:mrow>
                            <mml:mo>=</mml:mo>
                            <mml:mstyle displaystyle="true">
                              <mml:mfrac>
                                <mml:mrow>
                                  <mml:mn>2</mml:mn>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:msub>
                                    <mml:mi>C</mml:mi>
                                    <mml:mrow>
                                      <mml:mi>z</mml:mi>
                                      <mml:mo>⁢</mml:mo>
                                      <mml:mi>p</mml:mi>
                                    </mml:mrow>
                                  </mml:msub>
                                </mml:mrow>
                                <mml:mrow>
                                  <mml:msub>
                                    <mml:mi>C</mml:mi>
                                    <mml:mi>z</mml:mi>
                                  </mml:msub>
                                  <mml:mo>+</mml:mo>
                                  <mml:msub>
                                    <mml:mi>C</mml:mi>
                                    <mml:mi>p</mml:mi>
                                  </mml:msub>
                                </mml:mrow>
                              </mml:mfrac>
                            </mml:mstyle>
                          </mml:mrow>
                        </mml:math>
                      </disp-formula>
                      <disp-formula id="S3.E2Xa">
                        <mml:math alttext="\displaystyle C_{zp}=\sum_{a=1}^{M}\sum_{b=1}^{N}I_{z}(a,b)*I_{p}(a,b)" display="inline">
                          <mml:mrow>
                            <mml:msub>
                              <mml:mi>C</mml:mi>
                              <mml:mrow>
                                <mml:mi>z</mml:mi>
                                <mml:mo>⁢</mml:mo>
                                <mml:mi>p</mml:mi>
                              </mml:mrow>
                            </mml:msub>
                            <mml:mo>=</mml:mo>
                            <mml:mrow>
                              <mml:mstyle displaystyle="true">
                                <mml:munderover>
                                  <mml:mo movablelimits="false">∑</mml:mo>
                                  <mml:mrow>
                                    <mml:mi>a</mml:mi>
                                    <mml:mo>=</mml:mo>
                                    <mml:mn>1</mml:mn>
                                  </mml:mrow>
                                  <mml:mi>M</mml:mi>
                                </mml:munderover>
                              </mml:mstyle>
                              <mml:mrow>
                                <mml:mstyle displaystyle="true">
                                  <mml:munderover>
                                    <mml:mo movablelimits="false">∑</mml:mo>
                                    <mml:mrow>
                                      <mml:mi>b</mml:mi>
                                      <mml:mo>=</mml:mo>
                                      <mml:mn>1</mml:mn>
                                    </mml:mrow>
                                    <mml:mi>N</mml:mi>
                                  </mml:munderover>
                                </mml:mstyle>
                                <mml:mrow>
                                  <mml:mrow>
                                    <mml:mrow>
                                      <mml:msub>
                                        <mml:mi>I</mml:mi>
                                        <mml:mi>z</mml:mi>
                                      </mml:msub>
                                      <mml:mo>⁢</mml:mo>
                                      <mml:mrow>
                                        <mml:mo stretchy="false">(</mml:mo>
                                        <mml:mi>a</mml:mi>
                                        <mml:mo>,</mml:mo>
                                        <mml:mi>b</mml:mi>
                                        <mml:mo rspace="0.055em" stretchy="false">)</mml:mo>
                                      </mml:mrow>
                                    </mml:mrow>
                                    <mml:mo rspace="0.222em">∗</mml:mo>
                                    <mml:msub>
                                      <mml:mi>I</mml:mi>
                                      <mml:mi>p</mml:mi>
                                    </mml:msub>
                                  </mml:mrow>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:mrow>
                                    <mml:mo stretchy="false">(</mml:mo>
                                    <mml:mi>a</mml:mi>
                                    <mml:mo>,</mml:mo>
                                    <mml:mi>b</mml:mi>
                                    <mml:mo stretchy="false">)</mml:mo>
                                  </mml:mrow>
                                </mml:mrow>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:math>
                      </disp-formula>
                      <disp-formula id="S3.E2Xb">
                        <mml:math alttext="\displaystyle C_{z}=\sum_{a=1}^{M}\sum_{b=1}^{N}I_{z}(a,b)^{2}" display="inline">
                          <mml:mrow>
                            <mml:msub>
                              <mml:mi>C</mml:mi>
                              <mml:mi>z</mml:mi>
                            </mml:msub>
                            <mml:mo>=</mml:mo>
                            <mml:mrow>
                              <mml:mstyle displaystyle="true">
                                <mml:munderover>
                                  <mml:mo movablelimits="false">∑</mml:mo>
                                  <mml:mrow>
                                    <mml:mi>a</mml:mi>
                                    <mml:mo>=</mml:mo>
                                    <mml:mn>1</mml:mn>
                                  </mml:mrow>
                                  <mml:mi>M</mml:mi>
                                </mml:munderover>
                              </mml:mstyle>
                              <mml:mrow>
                                <mml:mstyle displaystyle="true">
                                  <mml:munderover>
                                    <mml:mo movablelimits="false">∑</mml:mo>
                                    <mml:mrow>
                                      <mml:mi>b</mml:mi>
                                      <mml:mo>=</mml:mo>
                                      <mml:mn>1</mml:mn>
                                    </mml:mrow>
                                    <mml:mi>N</mml:mi>
                                  </mml:munderover>
                                </mml:mstyle>
                                <mml:mrow>
                                  <mml:msub>
                                    <mml:mi>I</mml:mi>
                                    <mml:mi>z</mml:mi>
                                  </mml:msub>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:msup>
                                    <mml:mrow>
                                      <mml:mo stretchy="false">(</mml:mo>
                                      <mml:mi>a</mml:mi>
                                      <mml:mo>,</mml:mo>
                                      <mml:mi>b</mml:mi>
                                      <mml:mo stretchy="false">)</mml:mo>
                                    </mml:mrow>
                                    <mml:mn>2</mml:mn>
                                  </mml:msup>
                                </mml:mrow>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:math>
                      </disp-formula>
                      <disp-formula id="S3.E2Xc">
                        <mml:math alttext="\displaystyle C_{p}=\sum_{a=1}^{M}\sum_{b=1}^{N}I_{p}(a,b)^{2}" display="inline">
                          <mml:mrow>
                            <mml:msub>
                              <mml:mi>C</mml:mi>
                              <mml:mi>p</mml:mi>
                            </mml:msub>
                            <mml:mo>=</mml:mo>
                            <mml:mrow>
                              <mml:mstyle displaystyle="true">
                                <mml:munderover>
                                  <mml:mo movablelimits="false">∑</mml:mo>
                                  <mml:mrow>
                                    <mml:mi>a</mml:mi>
                                    <mml:mo>=</mml:mo>
                                    <mml:mn>1</mml:mn>
                                  </mml:mrow>
                                  <mml:mi>M</mml:mi>
                                </mml:munderover>
                              </mml:mstyle>
                              <mml:mrow>
                                <mml:mstyle displaystyle="true">
                                  <mml:munderover>
                                    <mml:mo movablelimits="false">∑</mml:mo>
                                    <mml:mrow>
                                      <mml:mi>b</mml:mi>
                                      <mml:mo>=</mml:mo>
                                      <mml:mn>1</mml:mn>
                                    </mml:mrow>
                                    <mml:mi>N</mml:mi>
                                  </mml:munderover>
                                </mml:mstyle>
                                <mml:mrow>
                                  <mml:msub>
                                    <mml:mi>I</mml:mi>
                                    <mml:mi>p</mml:mi>
                                  </mml:msub>
                                  <mml:mo>⁢</mml:mo>
                                  <mml:msup>
                                    <mml:mrow>
                                      <mml:mo stretchy="false">(</mml:mo>
                                      <mml:mi>a</mml:mi>
                                      <mml:mo>,</mml:mo>
                                      <mml:mi>b</mml:mi>
                                      <mml:mo stretchy="false">)</mml:mo>
                                    </mml:mrow>
                                    <mml:mn>2</mml:mn>
                                  </mml:msup>
                                </mml:mrow>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:math>
                      </disp-formula>
                    </disp-formula-group>
                  </p>
                  <p id="S3.SS1.p5"><bold>Signal to Noise Ratio</bold> is a performance metric to use to find the ratio between information and noise of the resultant image. SNR higher values express that both the reference and resultant images are similar [<xref rid="ref017" ref-type="bibr">17</xref>].</p>
                  <p>
                    <disp-formula id="S3.E3">
                      <mml:math alttext="SNR=10\log_{10}\left(\frac{\sum\limits_{a=1}^{M}\sum\limits_{b=1}^{N}(I_{z}(a,%&#10;b))^{2}}{\sum\limits_{a=1}^{M}\sum\limits_{b=1}^{N}(I_{z}(a,b)-I_{p}(a,b))^{2}%&#10;}\right)" display="block">
                        <mml:mrow>
                          <mml:mrow>
                            <mml:mi>S</mml:mi>
                            <mml:mo>⁢</mml:mo>
                            <mml:mi>N</mml:mi>
                            <mml:mo>⁢</mml:mo>
                            <mml:mi>R</mml:mi>
                          </mml:mrow>
                          <mml:mo>=</mml:mo>
                          <mml:mrow>
                            <mml:mn>10</mml:mn>
                            <mml:mo lspace="0.167em">⁢</mml:mo>
                            <mml:mrow>
                              <mml:msub>
                                <mml:mi>log</mml:mi>
                                <mml:mn>10</mml:mn>
                              </mml:msub>
                              <mml:mo>⁡</mml:mo>
                              <mml:mrow>
                                <mml:mo>(</mml:mo>
                                <mml:mfrac>
                                  <mml:mrow>
                                    <mml:munderover>
                                      <mml:mo movablelimits="false">∑</mml:mo>
                                      <mml:mrow>
                                        <mml:mi>a</mml:mi>
                                        <mml:mo>=</mml:mo>
                                        <mml:mn>1</mml:mn>
                                      </mml:mrow>
                                      <mml:mi>M</mml:mi>
                                    </mml:munderover>
                                    <mml:mrow>
                                      <mml:munderover>
                                        <mml:mo lspace="0.167em" movablelimits="false" rspace="0em">∑</mml:mo>
                                        <mml:mrow>
                                          <mml:mi>b</mml:mi>
                                          <mml:mo>=</mml:mo>
                                          <mml:mn>1</mml:mn>
                                        </mml:mrow>
                                        <mml:mi>N</mml:mi>
                                      </mml:munderover>
                                      <mml:msup>
                                        <mml:mrow>
                                          <mml:mo stretchy="false">(</mml:mo>
                                          <mml:mrow>
                                            <mml:msub>
                                              <mml:mi>I</mml:mi>
                                              <mml:mi>z</mml:mi>
                                            </mml:msub>
                                            <mml:mo>⁢</mml:mo>
                                            <mml:mrow>
                                              <mml:mo stretchy="false">(</mml:mo>
                                              <mml:mi>a</mml:mi>
                                              <mml:mo>,</mml:mo>
                                              <mml:mi>b</mml:mi>
                                              <mml:mo stretchy="false">)</mml:mo>
                                            </mml:mrow>
                                          </mml:mrow>
                                          <mml:mo stretchy="false">)</mml:mo>
                                        </mml:mrow>
                                        <mml:mn>2</mml:mn>
                                      </mml:msup>
                                    </mml:mrow>
                                  </mml:mrow>
                                  <mml:mrow>
                                    <mml:munderover>
                                      <mml:mo movablelimits="false">∑</mml:mo>
                                      <mml:mrow>
                                        <mml:mi>a</mml:mi>
                                        <mml:mo>=</mml:mo>
                                        <mml:mn>1</mml:mn>
                                      </mml:mrow>
                                      <mml:mi>M</mml:mi>
                                    </mml:munderover>
                                    <mml:mrow>
                                      <mml:munderover>
                                        <mml:mo lspace="0.167em" movablelimits="false" rspace="0em">∑</mml:mo>
                                        <mml:mrow>
                                          <mml:mi>b</mml:mi>
                                          <mml:mo>=</mml:mo>
                                          <mml:mn>1</mml:mn>
                                        </mml:mrow>
                                        <mml:mi>N</mml:mi>
                                      </mml:munderover>
                                      <mml:msup>
                                        <mml:mrow>
                                          <mml:mo stretchy="false">(</mml:mo>
                                          <mml:mrow>
                                            <mml:mrow>
                                              <mml:msub>
                                                <mml:mi>I</mml:mi>
                                                <mml:mi>z</mml:mi>
                                              </mml:msub>
                                              <mml:mo>⁢</mml:mo>
                                              <mml:mrow>
                                                <mml:mo stretchy="false">(</mml:mo>
                                                <mml:mi>a</mml:mi>
                                                <mml:mo>,</mml:mo>
                                                <mml:mi>b</mml:mi>
                                                <mml:mo stretchy="false">)</mml:mo>
                                              </mml:mrow>
                                            </mml:mrow>
                                            <mml:mo>−</mml:mo>
                                            <mml:mrow>
                                              <mml:msub>
                                                <mml:mi>I</mml:mi>
                                                <mml:mi>p</mml:mi>
                                              </mml:msub>
                                              <mml:mo>⁢</mml:mo>
                                              <mml:mrow>
                                                <mml:mo stretchy="false">(</mml:mo>
                                                <mml:mi>a</mml:mi>
                                                <mml:mo>,</mml:mo>
                                                <mml:mi>b</mml:mi>
                                                <mml:mo stretchy="false">)</mml:mo>
                                              </mml:mrow>
                                            </mml:mrow>
                                          </mml:mrow>
                                          <mml:mo stretchy="false">)</mml:mo>
                                        </mml:mrow>
                                        <mml:mn>2</mml:mn>
                                      </mml:msup>
                                    </mml:mrow>
                                  </mml:mrow>
                                </mml:mfrac>
                                <mml:mo>)</mml:mo>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:mrow>
                      </mml:math>
                    </disp-formula>
                  </p>
                  <p id="S3.SS1.p6"><bold>Peak Signal to Noise Ratio</bold> is a widely used performance metric, which is calculated by the number of gray levels in the image divided by the corresponding pixels in the reference and the Resultant images. The higher values indicate that the Resultant and reference images are similar [<xref rid="ref018" ref-type="bibr">18</xref>, <xref rid="ref019" ref-type="bibr">19</xref>].</p>
                  <p>
                    <disp-formula id="S3.E4">
                      <mml:math alttext="PSNR=20\log_{10}\left(\frac{G^{2}}{\frac{1}{M\times N}\sum\limits_{a=1}^{M}%&#10;\sum\limits_{b=1}^{N}(I_{z}(a,b)-I_{p}(a,b))^{2}}\right)" display="inline">
                        <mml:mrow>
                          <mml:mrow>
                            <mml:mi>P</mml:mi>
                            <mml:mo>⁢</mml:mo>
                            <mml:mi>S</mml:mi>
                            <mml:mo>⁢</mml:mo>
                            <mml:mi>N</mml:mi>
                            <mml:mo>⁢</mml:mo>
                            <mml:mi>R</mml:mi>
                          </mml:mrow>
                          <mml:mo>=</mml:mo>
                          <mml:mrow>
                            <mml:mn>20</mml:mn>
                            <mml:mo lspace="0.167em">⁢</mml:mo>
                            <mml:mrow>
                              <mml:msub>
                                <mml:mi>log</mml:mi>
                                <mml:mn>10</mml:mn>
                              </mml:msub>
                              <mml:mo>⁡</mml:mo>
                              <mml:mrow>
                                <mml:mo>(</mml:mo>
                                <mml:mfrac>
                                  <mml:msup>
                                    <mml:mi>G</mml:mi>
                                    <mml:mn>2</mml:mn>
                                  </mml:msup>
                                  <mml:mrow>
                                    <mml:mfrac>
                                      <mml:mn>1</mml:mn>
                                      <mml:mrow>
                                        <mml:mi>M</mml:mi>
                                        <mml:mo lspace="0.222em" rspace="0.222em">×</mml:mo>
                                        <mml:mi>N</mml:mi>
                                      </mml:mrow>
                                    </mml:mfrac>
                                    <mml:mo>⁢</mml:mo>
                                    <mml:mrow>
                                      <mml:mstyle displaystyle="false">
                                        <mml:munderover>
                                          <mml:mo movablelimits="false">∑</mml:mo>
                                          <mml:mrow>
                                            <mml:mi>a</mml:mi>
                                            <mml:mo>=</mml:mo>
                                            <mml:mn>1</mml:mn>
                                          </mml:mrow>
                                          <mml:mi>M</mml:mi>
                                        </mml:munderover>
                                      </mml:mstyle>
                                      <mml:mrow>
                                        <mml:mstyle displaystyle="false">
                                          <mml:munderover>
                                            <mml:mo movablelimits="false">∑</mml:mo>
                                            <mml:mrow>
                                              <mml:mi>b</mml:mi>
                                              <mml:mo>=</mml:mo>
                                              <mml:mn>1</mml:mn>
                                            </mml:mrow>
                                            <mml:mi>N</mml:mi>
                                          </mml:munderover>
                                        </mml:mstyle>
                                        <mml:msup>
                                          <mml:mrow>
                                            <mml:mo stretchy="false">(</mml:mo>
                                            <mml:mrow>
                                              <mml:mrow>
                                                <mml:msub>
                                                  <mml:mi>I</mml:mi>
                                                  <mml:mi>z</mml:mi>
                                                </mml:msub>
                                                <mml:mo>⁢</mml:mo>
                                                <mml:mrow>
                                                  <mml:mo stretchy="false">(</mml:mo>
                                                  <mml:mi>a</mml:mi>
                                                  <mml:mo>,</mml:mo>
                                                  <mml:mi>b</mml:mi>
                                                  <mml:mo stretchy="false">)</mml:mo>
                                                </mml:mrow>
                                              </mml:mrow>
                                              <mml:mo>−</mml:mo>
                                              <mml:mrow>
                                                <mml:msub>
                                                  <mml:mi>I</mml:mi>
                                                  <mml:mi>p</mml:mi>
                                                </mml:msub>
                                                <mml:mo>⁢</mml:mo>
                                                <mml:mrow>
                                                  <mml:mo stretchy="false">(</mml:mo>
                                                  <mml:mi>a</mml:mi>
                                                  <mml:mo>,</mml:mo>
                                                  <mml:mi>b</mml:mi>
                                                  <mml:mo stretchy="false">)</mml:mo>
                                                </mml:mrow>
                                              </mml:mrow>
                                            </mml:mrow>
                                            <mml:mo stretchy="false">)</mml:mo>
                                          </mml:mrow>
                                          <mml:mn>2</mml:mn>
                                        </mml:msup>
                                      </mml:mrow>
                                    </mml:mrow>
                                  </mml:mrow>
                                </mml:mfrac>
                                <mml:mo>)</mml:mo>
                              </mml:mrow>
                            </mml:mrow>
                          </mml:mrow>
                        </mml:mrow>
                      </mml:math>
                    </disp-formula>
                  </p>
                  <p>
                    <fig id="F1">
                      <label>Figure 1.</label>
                      <caption>
                        <p>The fused images and error images of different frequency domain methods on disk image set.</p>
                      </caption>
                      <graphic xlink:href="fig1.jpg"/>
                    </fig>
                  </p>
                </sec>
                <sec id="S3.SS2">
                  <label>3.2</label>
                  <title>Results and discussion</title>
                  <p id="S3.SS2.p1">In this study, we conducted a comparative analysis of various image fusion techniques, including DWT, SWT, and DCT, along with several variations of DCT-based methods. The specific variations of DCT methods we explored include DCT + SVD, DCT + SVD + CV, DCT + Correlation and Energy (DCT + Corr_Eng), and DCT + Correlation and Energy with CV (DCT + Corr_Eng + CV). The effectiveness of these fusion techniques was rigorously assessed through three different types of performance metrics: qualitative error image (QEI) analysis, quantitative measures, and qualitative assessments. In this letter, the experimentations are performed on Test image; “Clocks”. The grayscale image set is provided by “Lytro multi-focus datasets” [<xref rid="ref015" ref-type="bibr">15</xref>].</p>
                  <p id="S3.SS2.p2">To evaluate the quantitative performance, we utilized four specific metrics: entropy, SNR (Signal-to-Noise Ratio), PSNR (Peak Signal-to-Noise Ratio), and correlation. These were selected for their effectiveness in capturing various dimensions of fusion quality. Experiments were conducted on widely recognized multi-focus image datasets, specifically the disk dataset, with each image in the set having a resolution of 520×520 pixels.</p>
                  <p id="S3.SS2.p3">For the image fusion process, we focused on fusing two images at a time, although the algorithms used are flexible enough to handle more than two multi-focus images. Additionally, we used the qualitative error image (QEI) technique to evaluate the fusion results. The QEI is essentially a difference image obtained by subtracting the resultant fused image from a reference image [<xref rid="ref010" ref-type="bibr">10</xref>]. The less visible the QEI, the closer the fused image is to the reference image, indicating better fusion quality. This method provides a clear visual indication of how well the fusion process has preserved the important features of the original images [<xref rid="ref004" ref-type="bibr">4</xref>].</p>
                  <p>
                    <table-wrap id="T2">
                      <label>Table 2</label>
                      <caption>
                        <p>The quantitative results of a Disk image set.</p>
                      </caption>
                      <table>
                        <tr>
                          <td style="border-top: 1px solid black;" align="center">
                            <bold>Metrics</bold>
                          </td>
                          <td style="border-top: 1px solid black;" align="center">
                            <bold>DWT</bold>
                          </td>
                          <td style="border-top: 1px solid black;" align="center">
                            <bold>SWT</bold>
                          </td>
                          <td style="border-top: 1px solid black;" align="center">
                            <bold>DCT</bold>
                          </td>
                          <td style="border-top: 1px solid black;" align="center">
                            <bold>DCT+Corr_eng</bold>
                          </td>
                          <td style="border-top: 1px solid black;" align="center">
                            <bold>DCT+Corr_eng+CV</bold>
                          </td>
                          <td style="border-top: 1px solid black;" align="center">
                            <bold>DCT+SDV</bold>
                          </td>
                          <td style="border-top: 1px solid black;" align="center">
                            <bold>DCT+SDV+CV</bold>
                          </td>
                        </tr>
                        <tr>
                          <td style="border-top: 1px solid black;" align="center">entropy</td>
                          <td style="border-top: 1px solid black;" align="center">0.1790</td>
                          <td style="border-top: 1px solid black;" align="center">0.1795</td>
                          <td style="border-top: 1px solid black;" align="center">0.1834</td>
                          <td style="border-top: 1px solid black;" align="center">
                            <bold>0.1838</bold>
                          </td>
                          <td style="border-top: 1px solid black;" align="center">0.1675</td>
                          <td style="border-top: 1px solid black;" align="center">0.1656</td>
                          <td style="border-top: 1px solid black;" align="center">0.1690</td>
                        </tr>
                        <tr>
                          <td align="center">SNR</td>
                          <td align="center">15.6230</td>
                          <td align="center">15.9438</td>
                          <td align="center">15.6593</td>
                          <td align="center">15.6855</td>
                          <td align="center">15.5727</td>
                          <td align="center">15.5455</td>
                          <td align="center">
                            <bold>15.7767</bold>
                          </td>
                        </tr>
                        <tr>
                          <td align="center">PSNR</td>
                          <td align="center">35.9750</td>
                          <td align="center">36.1354</td>
                          <td align="center">35.9931</td>
                          <td align="center">36.0062</td>
                          <td align="center">35.9498</td>
                          <td align="center">35.9363</td>
                          <td align="center">
                            <bold>36.0519</bold>
                          </td>
                        </tr>
                        <tr>
                          <td style="border-bottom: 1px solid black;" align="center">Correlation</td>
                          <td style="border-bottom: 1px solid black;" align="center">0.9882</td>
                          <td style="border-bottom: 1px solid black;" align="center">0.9891</td>
                          <td style="border-bottom: 1px solid black;" align="center">0.9883</td>
                          <td style="border-bottom: 1px solid black;" align="center">0.9884</td>
                          <td style="border-bottom: 1px solid black;" align="center">0.9881</td>
                          <td style="border-bottom: 1px solid black;" align="center">0.9880</td>
                          <td style="border-bottom: 1px solid black;" align="center">
                            <bold>0.9886</bold>
                          </td>
                        </tr>
                      </table>
                    </table-wrap>
                  </p>
                  <p id="S3.SS2.p4">The qualitative results for all the fusion methods are illustrated in Figure <xref ref-type="fig" rid="F1">1</xref>, where we showcase both the final fused images and the corresponding error (or difference) images for each technique. While at first glance, the fused images across all methods may appear quite similar, with only minor differences, the true distinctions become evident when examining the qualitative error images. Notably, the fusion results and error images generated by the extended DCT + SVD approach are particularly impressive. It is clear that the difference image resulting from the DCT + SVD + CV method is more informative compared to DCT + SVD alone.</p>
                  <p id="S3.SS2.p5">In our observations, when we closely examine and compare the error images in Figure <xref ref-type="fig" rid="F1">1</xref> for all the fusion methods, it is evident that the qualitative performance of the DCT + SVD + CV method stands out across all three datasets. This approach consistently produces superior results, highlighting the effectiveness of incorporating CV into the fusion process.</p>
                  <p id="S3.SS2.p6">The quantitative results, evaluated using four different metrics, reveal some intriguing insights. The statistical values across all the fusion methods are very close, indicating that each method performs well. However, as shown in Table <xref rid="T2" ref-type="table">2</xref>, the DCT + SVD method stands out as the best performer for the disk image set. Despite the strong performance across the board, when we combine both visual and statistical assessments, it becomes clear that the DCT + SVD + CV method is the most superior among all the fusion techniques. This combination consistently outperforms the others, demonstrating its effectiveness in producing higher-quality fused images.</p>
                  <sec id="S4">
                    <label>4.</label>
                    <title>Conclusion</title>
                    <p id="S4.p1">This study offers an in-depth examination of various image fusion techniques, with a particular focus on methods based on the Discrete Cosine Transform (DCT) and its advanced variations. It emphasizes the effectiveness of DCT in the frequency domain, especially for multi-focus image fusion when combined with Singular Value Decomposition (SVD) and Consistency Verification (CV). The results demonstrate that the DCT + SVD + CV method consistently outperforms others in both qualitative and quantitative assessments, making it the most effective for achieving high-quality image fusion. Experimental findings suggest that while techniques such as Discrete Wavelet Transform (DWT) and Stationary Wavelet Transform (SWT) have their advantages, DCT-based methods, particularly when enhanced with SVD and CV, offer a more balanced approach to preserving image details and structural integrity. The DCT + SVD + CV technique stands out for its ability to retain essential image information while reducing artifacts, leading to superior overall fusion quality.</p>
                    <p id="S4.p2">Future research could explore several avenues for improving DCT-based fusion techniques. One potential area is the optimization of threshold values in Consistency Verification to adaptively enhance fusion accuracy for different image types. Additionally, incorporating machine learning models with DCT-based methods could automate the fusion process, potentially yielding even better results. Further research could also investigate extending these fusion techniques to handle more than two input images and applying them to different image modalities, such as infrared and visible light fusion, to test their robustness and versatility in various scenarios.</p>
                    <p id="S4.p3">
                      <bold>Data Availability Statement</bold>
                    </p>
                    <p id="S4.p4">Data will be made available on request.</p>
                    <p id="S4.p5">his work was supported without any funding.</p>
                    <p id="S4.p6">he authors declare no conflicts of interest.</p>
                    <p id="S4.p7">
                      <bold>Ethical Approval and Consent to Participate</bold>
                    </p>
                    <p id="S4.p8">Not applicable.</p>
                  </sec>
                </sec>
              </sec>
            </sec>
          </sec>
          <table-wrap-foot>
            <p>Frequency domain methods like DWT are effective for extracting frequency information from images but fall short in providing detailed directional information. The DWT approach has limitations such as the presence of ringing artifacts, longer processing times, problems with shift variance and additive noise, and higher energy consumption [<xref rid="ref010" ref-type="bibr">10</xref>]. To resolve the DWT issues the SWT was introduced and this approach is a fully shift-invariant transform that eliminates the down-sampling step of the decimated approach by up-sampling the filters, placing zeros among the filter coefficients. This approach has a simpler architecture and provides higher time-frequency localization [<xref rid="ref011" ref-type="bibr">11</xref>].</p>
            <p>However, due to the aforementioned issues with frequency methods, researchers have been increasingly interested in using the DCT for multi-focus image fusion [<xref rid="ref010" ref-type="bibr">10</xref>]. DCT-based methods are particularly efficient for transmitting and archiving images encoded in the JPEG standard. In the compressed domain, these approaches have performed better, avoiding the need for the complex and time-consuming decoding and encoding operations required by spatial-based algorithms [<xref rid="ref012" ref-type="bibr">12</xref>]. As a result, DCT-based multi-focus image fusion algorithms consume importantly less energy and time [<xref rid="ref010" ref-type="bibr">10</xref>].</p>
            <p>In this article, the focus is on DCT-based approaches, which come in many variations. These variations are thoroughly analyzed to expand our understanding of the functionality of different DCT approaches. Moreover, the article provides a comparative analysis with advanced approaches, including DWT, SWT, and DCT-based variations, highlighting their respective advantages and disadvantages in Table 1. The purpose of this comparison is to help new readers with basic concepts and explore potential modifications for new approaches. The DCT-based methods examined include DCT + SVD, DCT + SVD + CV, DCT + Corr_Eng, and DCT + Corr_Eng + CV.</p>
          </table-wrap-foot>
        </table-wrap>
      </p>
    </sec>
  </body>
  <back>
    <ack>
      <title>Acknowledgments</title>
      <p id="ack.p1">T</p>
    </ack>
    <sec id="sec0100" sec-type="COI-statement">
      <title>Conflict of interest</title>
      <p>T</p>
    </sec>
    <ref-list>
      <title>References</title>
      <ref id="ref001">
        <label>[1]</label>
        <mixed-citation> Farid, M.S., Mahmood, A., &amp; Al-Maadeed, S.A. (2019). Multi-focus image fusion using content adaptive blurring. <italic>Information Fusion</italic>, 45, 96–112. [<uri>https://doi.org/10.1016/j.inffus.2018.02.001</uri>] </mixed-citation>
      </ref>
      <ref id="ref002">
        <label>[2]</label>
        <mixed-citation> Muller, A. C., &amp; Narayanan, S. (2009). Cognitively-engineered multisensor image fusion for military applications. <italic>Information Fusion</italic>, 10(2), 137-149. [<uri>https://doi.org/10.1016/j.inffus.2008.08.008</uri>] </mixed-citation>
      </ref>
      <ref id="ref003">
        <label>[3]</label>
        <mixed-citation> Wang, J., Lu, T., Huang, X., Zhang, R., &amp; Feng, X. (2024). Pan-sharpening via conditional invertible neural network. <italic>Information Fusion</italic>, 101, 101980. [<uri>https://doi.org/10.1016/j.inffus.2023.101980</uri>] </mixed-citation>
      </ref>
      <ref id="ref004">
        <label>[4]</label>
        <mixed-citation> Bovith, T., Nielsen, A., Hansen, L., Overgaard, S., &amp; Gill, R. (2006, July). Detecting weather radar clutter by information fusion with satellite images and numerical weather prediction model output. In <italic>2006 IEEE International Symposium on Geoscience and Remote Sensing</italic> (pp. 511-514). IEEE. [<uri>https://doi.org/10.1109/IGARSS.2006.136</uri>] </mixed-citation>
      </ref>
      <ref id="ref005">
        <label>[5]</label>
        <mixed-citation> Zhang, S., Shen, X., Lin, Z., Měch, R., Costeira, J. P., &amp; Moura, J. M. (2018). Learning to understand image blur. In <italic>Proceedings of the IEEE conference on computer vision and pattern recognition</italic> (pp. 6586-6595). </mixed-citation>
      </ref>
      <ref id="ref006">
        <label>[6]</label>
        <mixed-citation> Morris, C., &amp; Rajesh, R. S. (2014, December). A novel and improved Spatial domain fusion method using Simple—PCA techniques. In <italic>2014 International Conference on Communication and Network Technologies</italic> (pp. 90-94). IEEE. [<uri>https://doi.org/10.1109/CNT.2014.7062731</uri>] </mixed-citation>
      </ref>
      <ref id="ref007">
        <label>[7]</label>
        <mixed-citation> Singh, G., Khosla, A., &amp; Anwar, M. I. (2016, February). Spatial domain color image enhancement based on local processing. In <italic>2016 3rd International Conference on Signal Processing and Integrated Networks (SPIN)</italic> (pp. 265-269). IEEE. [<uri>https://doi.org/10.1109/SPIN.2016.7566701</uri>] </mixed-citation>
      </ref>
      <ref id="ref008">
        <label>[8]</label>
        <mixed-citation> Ackar, H., Abd Almisreb, A., &amp; Saleh, M.A. (2019). A review on image enhancement techniques. <italic>Southeast Europe Journal of Soft Computing</italic>, 8(1). </mixed-citation>
      </ref>
      <ref id="ref009">
        <label>[9]</label>
        <mixed-citation> Khan, S.S., Khan, M., Alharbi, Y., Haider, U., Ullah, K., &amp; Haider, S. (2021). Hybrid Sharpening Transformation Approach for Multifocus Image Fusion Using Medical and Nonmedical Images. <italic>Journal of Healthcare Engineering</italic>, 2021. [<uri>https://doi.org/10.1155/2021/7000991</uri>] </mixed-citation>
      </ref>
      <ref id="ref010">
        <label>[10]</label>
        <mixed-citation> Amin-Naji, M., &amp; Aghagolzadeh, A. (2018). Multi-Focus Image Fusion in DCT Domain using Variance and Energy of Laplacian and Correlation Coefficient for Visual Sensor Networks. <italic>Journal of AI and Data Mining</italic>, 6(2), 233–250. [<uri>https://doi.org/10.22044/JADM.2017.5169.1624</uri>] </mixed-citation>
      </ref>
      <ref id="ref011">
        <label>[11]</label>
        <mixed-citation> Gharbia, R., Hassanien, A.E., El-Baz, A.H., Elhoseny, M., &amp; Gunasekaran, M. (2018). Multi-spectral and panchromatic image fusion approach using stationary wavelet transform and swarm flower pollination optimization for remote sensing applications. <italic>Future Generation Computer Systems</italic>, 88(11), 501–511. [<uri>https://doi.org/10.1016/j.future.2018.06.010</uri>] </mixed-citation>
      </ref>
      <ref id="ref012">
        <label>[12]</label>
        <mixed-citation> Tang, J., Peli, E., &amp; Acton, S. (2003). Image enhancement using a contrast measure in the compressed domain. <italic>IEEE Signal Processing Letters</italic>, 10(10), 289–292. [<uri>https://doi.org/10.1109/LSP.2003.817178</uri>] </mixed-citation>
      </ref>
      <ref id="ref013">
        <label>[13]</label>
        <mixed-citation> Amin-Naji, M., Ranjbar-Noiey, P., &amp; Aghagolzadeh, A. (2017). Multi-focus image fusion using singular value decomposition in DCT domain. <italic>2017 10th Iranian Conference on Machine Vision and Image Processing (MVIP)</italic>, 45–51. </mixed-citation>
      </ref>
      <ref id="ref014">
        <label>[14]</label>
        <mixed-citation> Rajakumar, C., &amp; Satheeskumaran, S. (2022). Singular value decomposition and saliency-map based image fusion for visible and infrared images. <italic>International Journal of Image and Data Fusion</italic>, 13(1), 21–43. [<uri>https://doi.org/10.1080/19479832.2021.1976832</uri>] </mixed-citation>
      </ref>
      <ref id="ref015">
        <label>[15]</label>
        <mixed-citation> Nejati, Mansour. (2016). Lytro Multi-focus Image Dataset. [<uri>https://doi.org/10.13140/RG.2.1.4802.6647</uri>] </mixed-citation>
      </ref>
      <ref id="ref016">
        <label>[16]</label>
        <mixed-citation> Karunasingha, D. S. K. (2022). Root mean square error or mean absolute error? Use their ratio as well. <italic>Information Sciences</italic>, 585, 609-629. [<uri>https://doi.org/10.1016/j.ins.2021.11.036</uri>] </mixed-citation>
      </ref>
      <ref id="ref017">
        <label>[17]</label>
        <mixed-citation> Moushmi, S., Sowmya, V., &amp; Soman, K. P. (2016). Empirical wavelet transform for multifocus image fusion. In <italic>Proceedings of the International Conference on Soft Computing Systems: ICSCS 2015, Volume 1</italic> (pp. 257-263). Springer India. [<uri>https://doi.org/10.1007/978-81-322-2671-0_25</uri>] </mixed-citation>
      </ref>
      <ref id="ref018">
        <label>[18]</label>
        <mixed-citation> Shah, M., et al. (2023). Multi-Focus Image Fusion using Unsharp Masking with Discrete Cosine Transform, 1–5. </mixed-citation>
      </ref>
      <ref id="ref019">
        <label>[19]</label>
        <mixed-citation> Khan, S.S., Ran, Q., &amp; Khan, M. (2020). Image pan-sharpening using enhancement based approaches in remote sensing. <italic>Multimedia Tools and Applications</italic>, 79(43), 32791-32805. [<uri>https://doi.org/10.1007/s11042-020-09682-z</uri>] </mixed-citation>
      </ref>
      <ref id="ref020">
        <label>[20]</label>
        <mixed-citation> Kumar, S., &amp; B. K. (2013). Multifocus and multispectral image fusion based on pixel significance using discrete cosine harmonic wavelet transform. <italic>Signal, Image and Video Processing</italic>, 7, 1125-1143. [<uri>https://doi.org/10.1007/s11760-012-0368-3</uri>] </mixed-citation>
      </ref>
      <ref id="ref021">
        <label>[21]</label>
        <mixed-citation> Lee, J., Vijaykrishnan, N., Irwin, M. J., &amp; Radhakrishnan, R. (2004). Inverse discrete cosine transform architecture exploiting sparseness and symmetry properties. <italic>IEEE Workshop on Signal Processing Systems (SIPS)</italic>, 361–366. [<uri>https://doi.org/10.1109/SIPS.2004.1363061</uri>] </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>
