An Improved SIFT Underwater Image Stitching Method

: Underwater image stitching is a technique employed to seamlessly merge images with overlapping regions, creating a coherent underwater panorama. In recent years, extensive research efforts have been devoted to advancing image stitching methodologies for both terrestrial and underwater applications. However, existing image stitching methods, which do not utilize detector information, heavily rely on matching feature pairs and tend to underperform in situations where underwater images contain regions with blurred feature textures. To address this challenge, we present an improved scale-invariant feature transform (SIFT) underwater image stitching method. This method enables the stitching of underwater images with arbitrarily acquired images featuring blurred feature contours and that do not require any detector information. Speciﬁcally, we perform a coarse feature extraction between the reference and training images, and then we acquire the target image and perform an accurate feature extraction between the reference and target images. In the ﬁnal stage, we propose an improved fade-in and fade-out fusion method to obtain a panoramic underwater image. The experimental results show that our proposed method demonstrates enhanced robustness, particularly in scenarios where detecting feature points is challenging, when compared to traditional SIFT methods. Additionally, our method achieves higher matching accuracy and produces higher-quality results in the stitching of underwater images.


Introduction
The process of stitching underwater images finds extensive applications in the realms of undersea exploration, marine archaeology, marine biological research, and the visualization of wrecks.Underwater conditions are characterized by the attenuation and backscattering of light, which restricts the field of view that cameras can capture.Consequently, it becomes unfeasible to encompass a larger area within a single frame.Nonetheless, panoramic imaging offers a solution by amalgamating images from diverse viewpoints into a single composite image that spans a wider area [1].This approach furnishes a more comprehensive reservoir of image data for subsequent research endeavors.
In recent years, significant advancements have been made in the fields of image stitching, mosaic construction, and full three-dimensional (3D) reconstruction.However, it is important to emphasize that these methods are primarily designed for terrestrial scenarios and do not directly address the unique challenges encountered in underwater imaging [2].Underwater image stitching presents distinct challenges compared to its land-based counterpart.The first major challenge lies in obtaining clear underwater images, given the complex and dynamic underwater environment.Underwater images are susceptible to blurring caused by water flow and the presence of microorganisms in the aquatic medium [3].The second challenge is associated with the complexity of underwater subjects.Underwater elements such as marine life, rocks, and wrecks often lack distinct texture features, resulting in blurry outlines that are difficult to extract.Additionally, marine organisms like seaweed, oysters, sea snails, earthworms, and sea grass can attach to Appl.Sci.2023, 13, 12251 2 of 22 shipwrecks, further obscuring their outlines and features.This complicates the accuracy of feature matching.Furthermore, the movement of marine life due to water currents can introduce alignment errors and ghosting artifacts in underwater images.The third challenge pertains to the prevalence of flat areas in underwater images, especially on extensive ocean floors [4].These regions exhibit similar textures around feature points, making it challenging to extract effective feature points.The fourth challenge is related to the assumptions underlying mosaic construction techniques, which often assume that the images are captured by a purely rotating camera, or that the scene is planar [5].These assumptions imply the absence of parallax caused by camera movement and a lack of three-dimensional effects.However, underwater images are subject to light attenuation and backscattering, which can result in significant distortions in the mosaic.These distortions undermine the applicability of these assumptions in underwater contexts.
To address the challenges of underwater image stitching, in recent years, many researchers have conducted work on underwater image stitching and video stitching.Underwater mosaicking has predominantly been explored within the realm of vision-based navigation and maintaining a stable position near the seafloor [6].However, these approaches are typically tailored to specific conditions and are unable to seamlessly stitch together any set of acquired images without external guidance [7].When it comes to underwater images captured by divers, they often lack the regularity of fixed time or distance intervals, and sensor data are often unavailable.Consequently, to facilitate the stitching of arbitrary underwater images, many researchers have turned to image registration methods like SIFT and SURF.These methods rely solely on image information and do not necessitate the use of detection devices or impose any special requirements specific to underwater images.Nonetheless, it is important to note that these image stitching methods heavily hinge on the information derived from matching pairs of distinctive image features.This reliance places significant demands on both the quantity and accuracy of these feature matches [8].Insufficient feature matches can pose a considerable challenge, especially in the context of underwater environments.Scale-invariant feature transform (SIFT) and speeded-up robust features (SURF) stand out as potent tools in terms of their ability to handle image transformation and deformation invariance.Lowe [9] proposed a method for extracting unique invariant features from images, which are invaluable for reliable matching across different objects or scenes.These features maintain their invariance to image scaling and rotation, and remarkably, SIFT remains robust even in the presence of variations in image illumination or affine distortion.SIFT's capacity to generate a substantial number of feature vectors, even from a small number of objects, renders it particularly suited for precise matching within large-scale feature databases.Bay et al. [10] introduced the SURF algorithm in response to the slow computational speed of SIFT.However, SURF exhibits lower stability and matching accuracy compared to SIFT, making it less effective in the context of stitching underwater images, particularly in regions with blurred features.In the deep sea, where underwater light rapidly diminishes due to attenuation, natural illumination becomes scarce beyond just a few tens of meters in depth.Consequently, most optical imaging in the deep sea relies on artificial lighting, often resulting in uneven illumination.Therefore, the use of SIFT for image registration in underwater image stitching is highly appropriate.In regions with rich and stable texture features, SIFT excels at extracting a sufficient number of features, enabling the generation of accurate matching pairs.However, underwater stitching frequently encounters challenges in flat areas with limited texture, unclear target contours in captured images, and the presence of underwater microorganisms, all of which can further complicate image quality.Solely relying on SIFT for feature extraction in such cases may result in an insufficient feature set and an inadequate number of accurate matching pairs, ultimately leading to alignment errors and ghosting artifacts in the stitching process.In contrast to existing stitching techniques, our goal is to introduce an underwater image stitching method that is capable of seamlessly merging arbitrarily acquired underwater images.This method remains effective even when the images exhibit blurred feature contours, and it does not depend on detector-related information.In this paper, we present an enhanced SIFT underwater image stitching method designed to address images with blurred feature contours without relying on detector information.
The overview of our approach is illustrated in Figure 1 and consists of three stages: the underwater image preprocessing stage, the SIFT-based coarse matching stage, and the SIFT-based precise matching stage.In the first stage, we employ the dynamic threshold white balance method and the contrast limited adaptive histogram equalization (CLAHE) algorithm to adjust colors and enhance brightness.In the second stage, we perform coarse feature extraction using the SIFT algorithm on both the reference image (image A) and the training image (image B).The training image represents the image to be transformed into the target image.Then, we match the features in both images using the nearest neighbor ratio algorithm and the random sample consensus (RANSAC) algorithm.Following this matching process, we rotate the training image using the homography matrix and correct the coordinate origin.Finally, any redundant blank areas are removed to obtain the target image (image C).In the third stage, we conduct accurate feature extraction based on the SIFT algorithm for both the reference image (image A) and the target image (image C).This precise feature extraction is followed by matching the features in both images, culminating in an image-paste-type fusion stitching.Ultimately, the panoramic image is created using the enhanced fade-in and fade-out method for seamless fusion.We employ a weighted smoothing algorithm to achieve the stitching of the target image.The weights assigned to the two images intended for merging undergo a gradual transition from 0 to 1 in either the horizontal or vertical direction, thereby achieving a seamless transition in the fused image.
In our experiments, we present compelling evidence demonstrating that our proposed method consistently yields excellent results, attesting to its effectiveness in the domain of underwater image stitching.
Appl.Sci.2023, 13, x FOR PEER REVIEW 3 of 24 when the images exhibit blurred feature contours, and it does not depend on detectorrelated information.In this paper, we present an enhanced SIFT underwater image stitching method designed to address images with blurred feature contours without relying on detector information.The overview of our approach is illustrated in Figure 1 and consists of three stages: the underwater image preprocessing stage, the SIFT-based coarse matching stage, and the SIFT-based precise matching stage.In the first stage, we employ the dynamic threshold white balance method and the contrast limited adaptive histogram equalization (CLAHE) algorithm to adjust colors and enhance brightness.In the second stage, we perform coarse feature extraction using the SIFT algorithm on both the reference image (image A) and the training image (image B).The training image represents the image to be transformed into the target image.Then, we match the features in both images using the nearest neighbor ratio algorithm and the random sample consensus (RANSAC) algorithm.Following this matching process, we rotate the training image using the homography matrix and correct the coordinate origin.Finally, any redundant blank areas are removed to obtain the target image (image C).In the third stage, we conduct accurate feature extraction based on the SIFT algorithm for both the reference image (image A) and the target image (image C).This precise feature extraction is followed by matching the features in both images, culminating in an image-paste-type fusion stitching.Ultimately, the panoramic image is created using the enhanced fade-in and fade-out method for seamless fusion.We employ a weighted smoothing algorithm to achieve the stitching of the target image.The weights assigned to the two images intended for merging undergo a gradual transition from 0 to 1 in either the horizontal or vertical direction, thereby achieving a seamless transition in the fused image.In our experiments, we present compelling evidence demonstrating that our proposed method consistently yields excellent results, attesting to its effectiveness in the domain of underwater image stitching.The primary contribution of this paper is the development of an enhanced SIFT underwater image stitching method strategically designed to minimize the ghosting artifacts commonly associated with stitching.Our approach excels even in areas characterized by blurred texture features, where traditional feature extraction methods struggle.By significantly increasing the number of accurate matching pairs, we enhance the matching The primary contribution of this paper is the development of an enhanced SIFT underwater image stitching method strategically designed to minimize the ghosting artifacts commonly associated with stitching.Our approach excels even in areas characterized by blurred texture features, where traditional feature extraction methods struggle.By significantly increasing the number of accurate matching pairs, we enhance the matching accuracy, improve the alignment of feature points between images, reduce stitching errors, and ultimately elevate the quality of underwater panoramic images.In a series of ten un-derwater image stitching experiments, our proposed method outperformed the traditional SIFT algorithm, increasing the non-reference underwater image quality measure (UIQM) by an average of 6.04%.Furthermore, when compared to references [11][12][13], our method enhances the UIQM index by 4.73%, 9.01%, and 11.06%, respectively.
The remainder of this paper is structured as follows: We describe the related works in Section 2. In Section 3, we describe our proposed improved SIFT underwater image stitching method.Section 4 describes the experiments.Finally, we present our conclusions in Section 5.

Related Work
In this section, we briefly review image stitching methods, both on land and underwater.While research on underwater image stitching has been relatively limited in recent years, notable advancements have been achieved in land-based image stitching methods.In land-based image stitching, a spectrum of approaches has been put forward to enhance local alignment.This includes the computation of multi-content-aware local warping as opposed to a singular global warping strategy.Moreover, some methodologies advocate for the use of seam-driven image stitching strategies aimed at mitigating artifacts resulting from image fusion [14,15].
Lowe et al. [16] proposed the AutoStitch image stitching algorithm, which leverages a global single-response matrix to establish correspondences across all images and generate a panoramic result through multiband fusion.This approach enabled the automated stitching of multiple images.However, a critical constraint of the AutoStitch algorithm is its reliance on the assumption of a consistent scene depth across overlapping image areas.Deviations from this assumption can lead to pronounced ghosting and alignment errors.To mitigate this limitation, researchers have proposed employing more robust matrices or multiple single-response matrices in place of a single global one.Lin et al. [17] proposed the smoothly varying affine (SVA) method, which incorporates smoothly varying affine stitching fields instead of global affine transformations.This technique exhibits superior local deformation and alignment capabilities and can also handle parallax to some extent.Gao et al. [18] proposed the dual-homography warping (DHW) approach, which divides the scene into a background plane and a foreground plane.Two distinct single-response matrices are employed, each aligned to the background and foreground, facilitating seamless stitching of a broader range of realistic scenes.Zaragoza et al. [11] proposed As-Projective-As-Possible (APAP).Firstly, SIFT features are extracted from two images, and erroneous features are removed using RANSAC.Then, an efficient computational method, moving DLT, is given to estimate the global homography matrix.When two images exhibit variations beyond mere rotation, or when they depict nonplanar scenes, moving DLT employs a location-dependent homography to warp each matching point.In the APAP algorithm, the image is divided into dense grids according to a certain pattern.Higher weights are provided for closer feature points, and lower weights are provided for more distant feature points.Each grid is aligned with a local single-response matrix, which is calculated from the global single-response matrix and the grid weights.Finally, the image stitching is performed using the local homography matrix.The APAP algorithm can align multiple images better via local weight constraint adjustment.Based on this, many researchers have been inspired to propose the use of grid optimization to achieve image stitching.However, since APAP performs affine transformation in non-overlapping regions of the image as well, it can cause severe visual distortion.Chang et al. [12] proposed shape-preserving half-projective (SPHP), which is based on the shape-preserving method of image scaling from the perspective of shape correction.The algorithm performs distorted alignment with single-strain transformations in overlapping regions and maintains the perspective of the image in non-overlapping regions with similarity transformations.At the same time, a smooth transition is formed between the overlapping and non-overlapping regions.The SPHP algorithm can reduce the perspective distortion in non-overlapping regions.The authors also proposed combining SPHP with APAP to make the proposed warping more robust.Lin et al. [13] proposed adaptive As-Natural-As-Possible (AANAP).This method uses a smooth stitching field over the entire target image while taking into account all local transformation variations.It is less dependent on the choice of parameters and can automatically compute the appropriate global similarity transform.The AANAP algorithm can also mitigate the perspective distortion in non-overlapping regions by linearizing the homography and slowly changing it into a global similarity transform.
Zhang et al. [19] developed an effective method to find a local alignment for optimal image stitching.This local stitching method uses a hybrid alignment model that simultaneously employs homography and content-preserving warping.Jing et al. [20] proposed a parallax-tolerant image stitching method based on robust elastic warping.They constructed an analysis warping function to eliminate parallax errors.Then, the input images are warped over the meshed image plane according to the calculated deformations.However, the effectiveness of these land-based methods in stitching images with blurry underwater conditions remains to be demonstrated.
Due to the complexity of the underwater environment, underwater image stitching techniques pose significant challenges.For more than thirty years, researchers have been studying to create automatic mosaics for underwater applications.For underwater mosaicking, some studies focus on vision-based navigation, while others focus on building large mosaics [21].Some algorithms can achieve image stitching but require assistance from external detectors or have specific requirements for underwater image acquisition.Fortunately, SIFT-based stitching methods can stitch images using only image information, but they require a high number and accuracy of feature matches [3].Insufficient feature matching can lead to stitching errors.
One of the earliest computer-assisted underwater image stitching systems was proposed by Haywood [22] in 1986.This is a manual system that uses camera positions to directly obtain the position relationship and coordinates of adjacent frames.The motion parameters between frames can be computed directly from the camera positions, without requiring any registration.However, this system is not automatic or real-time.
Negahdaripour et al. [6] developed a method for estimating camera motion in real time directly from the luminance variations of the image sequence.In [23], direct estimates of motion were successfully applied in underwater mosaic schemes.To achieve video mosaics along unconstrained vehicle paths, Fleischer et al. [24] utilized optimal estimation theory and smoother-follower techniques to identify the errors and reduce the image alignment errors propagating through the image chain.In [25], Fleischer et al. extended the theory of smoother-follower estimation to achieve real-time creation and online improvement of unconstrained underwater video mosaics.However, the image estimates depend on the sensor data and the system geometry.Richard et al. [7] used real-time image registration to achieve real-time creation of underwater video mosaics.Relative offsets between a pair of images can be used to determine when a new image should be acquired and where it belongs in the mosaic.In order to obtain images suitable for creating mosaics, the images are acquired snapping at a fixed time interval, snapping at a fixed position interval, or snapping at a fixed visual interval.Gracias et al. [26] proposed an integrated approach to solve the problem of pose estimation and mosaic construction in underwater mosaicking.This method uses a simple image motion model to compute the image motion and focus on the use of long image sequences with temporal distance overlays, such as those generated by circular trajectories or zigzag scanning patterns.However, this approach is based on the assumption that the dominant motion existing between consecutive frames is translation, and it does not account for lens distortion.
However, these approaches are unable to construct underwater mosaics without detector information or on arbitrarily obtained underwater images.Therefore, researchers have conducted further research on underwater mosaics of arbitrary images by using only the intrinsic information of the images.
Leone et al. [22] introduced a fully automated feature-based method for creating seamless, panoramic 2D mosaics from a sequence of images.Their texture-based approach, combined with statistical filtering, ensures accurate correspondences between adjacent frames and facilitates the removal of outliers.Raut et al. [8] proposed an enhancement to the SIFT algorithm to improve the automatic recognition of underwater images.Their method incorporates a Gabor filter as a pre-filter and utilizes the Hausdorff distance to compute the distance between key points, resulting in the detection of more key points and improved matching.Pizarro et al. [27] presented an approach for large-area underwater mosaicking, addressing challenges related to low overlap, unstructured motion, and varying lighting conditions in underwater images.Their method involves three key steps: determining image homography by estimating and compensating for radial distortion, topology estimation through feature-based pairwise image registration using a multiscale Harris interest point detector and a Zernike moment-based feature descriptor, and global alignment of all images based on the initial alignment obtained from the pairwise estimation.However, this approach heavily relies on effective feature detection and matching phases.
As far as our knowledge extends, there remain significant challenges in detecting a sufficient number of accurate feature pairs when the image feature contours are blurred.Insufficient feature matches can lead to errors in the mosaic.

Materials and Methods
To perform the stitching of underwater images with blurred feature contours and without relying on detector information, we present an enhanced stitching method tailored for SIFT underwater images.As illustrated in Figure 1, our approach encompasses three distinct stages, each serving a specific purpose.In this section, we commence by detailing the preprocessing steps applied to underwater images.Subsequently, we provide an overview of the traditional SIFT method, a widely employed technique in image stitching.Finally, we introduce our proprietary stitching algorithm.

Underwater Image Pre-Processing
When light traverses through water, a portion of it is reflected, while the remainder penetrates the water and is selectively absorbed.Different wavelengths of light experience varying degrees of absorption in water, with red light being significantly attenuated compared to blue-green light.Consequently, underwater images predominantly exhibit blue-green hues and tend to grow darker with increasing depth.Moreover, underwater images are susceptible to blurring and haziness due to the scattering of light by water, which can impair the accuracy of feature point extraction during the stitching process.Therefore, preprocessing of underwater images plays a crucial role, enabling the extraction of a larger number of feature points with enhanced precision.
Underwater image preprocessing can be categorized into two main types: image enhancement methods, and image restoration techniques.While image restoration methods necessitate prior knowledge of the environmental conditions, image enhancement methods offer a more straightforward approach that does not rely on such specific information.As such, this paper employs image enhancement methods.Specifically, dynamic threshold white balance and CLAHE are utilized for the preprocessing of underwater images.The CLAHE method [28] enhances local contrast in images by applying a threshold to limit contrast amplification within localized regions, addressing issues such as noise overamplification that may occur in the adaptive histogram equalization (AHE) method [29].The dynamic threshold white balance method [30], rooted in the YCbCr color space, offers image enhancement capabilities that enable adjustments to image tones and the restoration of accurate color representation.This ensures faithful rendering of the true colors of the captured subject in the images obtained by the camera.
In this paper, we improve the quality of underwater images by optimizing their color and contrast.The dynamic threshold white balance method is used to adjust the color, while the CLAHE algorithm enhances the contrast and brightness.We set the contrast limit parameter of CLAHE to 2 and use a 4 × 4 grid size to block the image.
We selected a set of captured images as the preprocessed dataset.The results are shown in Figure 2. As can be seen, there was a significant improvement in image quality.Moreover, we used different methods to preprocess the images and perform image feature point matching.We compared our method with the integrated color model (ICM) [31], underwater light attenuation prior (ULAP) [32], underwater dark channel prior (UDCP) [33], and image blurriness and light absorption (IBLA) [34].The results of the underwater image quality assessment are summarized in Table 1.The accuracy rate of interest point matching is improved in the preprocessed images, with more features correctly matched.Specifically, our method improves the accuracy rate from 44.76% in the original images to 59.06% in the preprocessed images.Compared with other algorithms, our algorithm has higher accuracy in feature matching.Overall, our proposed method effectively enhances the detection of interest points in underwater images.
restoration of accurate color representation.This ensures faithful rendering of the true colors of the captured subject in the images obtained by the camera.
In this paper, we improve the quality of underwater images by optimizing their color and contrast.The dynamic threshold white balance method is used to adjust the color, while the CLAHE algorithm enhances the contrast and brightness.We set the contrast limit parameter of CLAHE to 2 and use a 4 × 4 grid size to block the image.
We selected a set of captured images as the preprocessed dataset.The results are shown in Figure 2. As can be seen, there was a significant improvement in image quality.Moreover, we used different methods to preprocess the images and perform image feature point matching.We compared our method with the integrated color model (ICM) [31], underwater light attenuation prior (ULAP) [32], underwater dark channel prior (UDCP) [33], and image blurriness and light absorption (IBLA) [34].The results of the underwater image quality assessment are summarized in Table 1.The accuracy rate of interest point matching is improved in the preprocessed images, with more features correctly matched.Specifically, our method improves the accuracy rate from 44.76% in the original images to 59.06% in the preprocessed images.Compared with other algorithms, our algorithm has higher accuracy in feature matching.Overall, our proposed method effectively enhances the detection of interest points in underwater images.As described in the literature [9], the SIFT algorithm utilizes a 128-dimensional feature vector to describe feature points in an image with scale, rotation, and translation  As described in the literature [9], the SIFT algorithm utilizes a 128-dimensional feature vector to describe feature points in an image with scale, rotation, and translation invariance.The matching process is based on the similarity between these feature point descriptors.
Firstly, SIFT detects preliminary extrema points through the construction of a Gaussian difference pyramid.It then filters low-contrast and edge-response points to obtain robust key points.Subsequently, SIFT assigns an orientation to each key point, calculates the gradient magnitude and orientation, and generates orientation gradient histograms.The radius of the image area to be calculated is then determined, and the coordinate axis is rotated to align with the main direction of the key point.An 8 × 8 window surrounding the key point is selected as the feature descriptor, and a feature point descriptor is generated [9].Finally, SIFT performs feature point matching.
In the feature point matching, the reference image (image A) is denoted by I and the training image (image B) is denoted by I .The expression is as follows: where x = x y T and x = x y T denote the pixel coordinates of the reference image and the pixel coordinates of the training image, respectively, while H ∈ R 3×3 represents the homography matrix.The flush coordinates of the target image and the reference image pixel points are expressed as follows: The homography matrix has a degree of freedom of 8. Therefore, we assume that h 9 = 1, and the above equation can be organized as follows: Organize Equations ( 3) and ( 4) into the following form: where x 1 y 1 T and x 1 y 1 T indicate the coordinates of the first set of matched pairs.
When we obtain the coordinates of the 4 sets of matched pairs, we can solve the above system of equations to obtain the homography matrix.

The Proposed Improved SIFT Method
An improved SIFT method is proposed to address the drawback of only performing feature extraction once in the traditional SIFT algorithm.The reference image (image A) is denoted by I and the target image (image C) is denoted by I .The expression is as follows: where x = x y T and x = x y T denote the pixel coordinates in the reference image and the pixel coordinates in the target image, respectively, H ∈ R 3×3 represents the homography matrix from the training image (image B) to the reference image (image A), and H ∈ R 3×3 represents the homography matrix from the target image (image C) to the training image (image B).Equation (7) shows the projection transformation of the target image (image C) to the reference image (image A).Again, the above equation can be organized as follows: During the stitching phase, we undertake a two-step approach involving coarse feature extraction and precise feature extraction, both based on the SIFT algorithm.Initially, coarse feature extraction is conducted on the reference image (image A) and the training image (image B).Subsequently, we employ the k-nearest neighbor algorithm (KNN) to gauge the similarity of key points between these two images.The KNN method utilizes the Euclidean distance as the measure for assessing feature point similarity.It performs feature point matching by systematically traversing all key points in the two images slated for concatenation.The Euclidean distance between any two corresponding feature points within these two images to be joined (D) can be defined as follows: where F i Q and F j T represent feature points in the reference and training images, respectively, K represents the dimensions of the feature points, and N is the number of feature points in the training image.Next, we iterate through all of the feature points in the reference image, identifying those within the training image that possess the closest Euclidean distances t 1 and the second-closest Euclidean distances t 2 .After applying this selection process to the feature points, we can establish potential matching pairs.The calculation formula is provided below: where V f represents a predetermined threshold, typically within the range of 0.4 to 0.8.In this paper, we set V f to 0.75.For our underwater images, this threshold was thoughtfully selected to maintain a specific quantity of feature matching pairs while ensuring a high level of matching accuracy.If result is below 0, it signifies that the distance ratio is below the defined threshold, indicating a successful match of feature pairs.Conversely, if result exceeds 0, it implies that the distance ratio surpasses the specified threshold, and subsequently, the removal of the potential matching pair.
To further improve the accuracy of feature alignment and stitching, the RANSAC algorithm is employed.Through multiple iterations, it fits the optimal mathematical model to eliminate matching pairs with significant errors.Firstly, a subset of elements is selected from the set of coarsely matched feature points and designated as the inlier group.Within the inlier group, each element is incorporated into the model to compute the model parameters.The calculation formula of the model is shown in Formula (2).Secondly, the remaining points from the feature point set, excluding those in the inlier group, are inserted into the parameter model obtained in the previous steps.The distance from the model is computed for each element point, and points that fall within the allowable error range are marked as inlier elements, while those outside this range are designated as outliers.The error threshold is set to 4. Utilizing all of the element points contained within the inlier group at this stage, we repeat steps 1 and 2 N times, where the expression for the number of iterations, N, is as follows: where p represents the confidence level, typically in the range of 0.95 to 0.99, φ represents the proportion of inlier element points in the dataset, and m represents the total number of feature points in the dataset.The best-fitting model, containing the most inlier element points during the iterative process, is chosen as the optimal model.The element points within this inlier group at this stage are considered to be the final feature point matches, and the corresponding homography matrix H is calculated.
Once the homography matrix H for image A and image B is obtained, multiply this matrix with image B and conduct a coordinate system correction.The correction aligns two of the four vertices of the training image with the X and Y axes, ensuring that all coordinate points are positive.This enhancement significantly accelerates the computational speed of subsequent stitching.Finally, the superfluous blank areas resulting from the rotation are eliminated to yield the target image (image C).The proposed SIFTbased coarse matching stage is shown in Algorithm 1.The SIFT-based coarse matching stage employs the next input parameters-the reference image (image A) re f erenceImg and the training image (image B) trainImg-and employs the next output parameter: the target image (image C) targetImg; kpsRe f erence and kpsTrain represent the key points detected in image A and image B, respectively; re f erenceDescriptors and trainDescriptors represent the feature descriptors for key points detected in image A and image B, respectively; matches represents the matching feature point pairs obtained using the KNN method; H represents the estimated perspective transformation matrix; h q and w q represent the height and width of image A, respectively; origin represents the boundary corners of the original image; pts represents the coordinates of these corner points; dst represents the positions of these corner points on the target image after perspective transformation.Input: re f erenceImg, trainImg Output: targetImg kpsRe f erence, re f erenceDescriptors kpsTrain, trainDescriptors matches ← matchKeyPointsKNN H ← getHomography h q , w q , _ ← shape (re f erenceImg) origin ← array [0, 0], 0, h q −1], w q -1, h q −1], w q -1, 0]]) pts ← origin.reshape (-1,1,2) dst ← perspectiveTransform (pts , H). squeeze () if min (dst ) < 0 then if min (dst [: , 0]) < min (dst [: ,1]) then dst←dst− min (dst In the second phase, precise feature extraction is carried out using the SIFT algorithm on both the target image (image C) and the reference image (image A).The SIFT algorithm is independently applied to extract features within the overlapping regions of these two images.During feature extraction, a bidirectional k-nearest neighbors (KNN) matching algorithm is utilized for feature point matching.This approach is particularly effective in capturing more precise feature points, especially in flat areas of the image.To elucidate the process, an initial key point is selected from image 1, the two closest feature points in image 2 to this key point in image 1 are identified, and their positions are recorded.If the distance ratio between the nearest feature point and the second-nearest feature point falls below the matching threshold of 0.75, the nearest feature point is retained.A similar procedure is then applied to determine the two closest feature points in image 1 for each feature point in image 2 among the nearest feature points in image 2.
To create a realistic large-scale underwater scene image, this paper introduces an image-pasting fusion stitching method.Firstly, the feature matching pairs are ranked in descending order according to their matching scores.Then, the logarithmic values of the matching feature points derived from the final filter are compared against a predefined threshold.If the values are below this threshold, the top 10% of the feature matching pairs with the highest scores are selected to determine the relative position of the two underwater images, denoted as x sub and y sub .The equations of the relative position correction symbols for the two images are as follows: where x sub is the relative position correction symbol in the horizontal direction, y sub is the relative position correction symbol in the vertical direction, (x A , y A ) and (x B , y B ) are the coordinates of matching feature points for the two images, respectively, and n represents the number of feature matching points between two images.If the number of feature pairs resulting from the final matching process exceeds the predefined threshold, Formulas ( 13) and ( 14) are employed to determine the relative positions of the features in all matched pairs.The threshold value is chosen from the range of 50 to 100, depending on the image size.Through experiments conducted on the images presented in this article, it was observed that setting the threshold below 50 leads to a significant increase in matching errors.Conversely, if the threshold exceeds 100, it results in reduced computational efficiency.Therefore, in this study, the threshold was set to 60.
The determination of relative positions is based on the polarity of the relative position correction symbol.When x sub is positive, it indicates that the image to be matched aligns to the left of the reference image, while the reference image is to the right.Conversely, when y sub is positive, it indicates that the image to be matched aligns to the top, while the reference image is below.Once the positions of the two images are established, the relative positions can be utilized for image stitching.
To enhance the quality of the stitched images, we employed an improved fusion algorithm in this study.This approach involves calculating the weighted sum of pixel values within the overlapping regions of two images, taking into account the distance between these regions and the two images.These weights play a critical role in seamlessly blending the overlapping areas, thus reducing the ghosting effects and visible stitch marks in the final image.The calculation formula for the weighted sum is as follows: where f is the pixel value of the fused image, (x, y) is the coordinate of the concatenated image, f 1 and f 2 denote the pixel values of the two original images, respectively, c 1 and c 2 are weighting factors, and we have After completing the stitching process, the identification of the stitching direction involves the determination of the overlapping area between two images.The fusion algorithm is divided into horizontal and vertical fusion categories, depending on the observed horizontal and vertical stitching characteristics during the image stitching process.The dimensions of the overlapping area, specifically its width and height, are crucial in approximating the stitching direction of the image.When the width of the overlapping part exceeds the height, the two images are categorized as being horizontally stitched, and the horizontal fusion method is applied.Conversely, if the height of the overlapping part is greater than the width, the two images are considered to be vertically stitched, and the vertical fusion method is employed.This algorithm enhances the stitching accuracy by considering the direction of stitching between the images.
In the horizontal fusion method, f 1 (x, y) and f 2 (x, y) denote the pixel values of the original image on the left and right sides of the overlapping area, respectively, while c 1 and c 2 are calculated as shown below: ) where x i indicates the horizontal coordinate of the current pixel point, x l indicates the left boundary of the overlapping area, x r indicates the right boundary of the overlapping area, and x r = x l .
In the vertical fusion method, f 1 (x, y) and f 2 (x, y) denote the pixel values of the original image at the top and bottom of the overlapping area, respectively; c 1 and c 2 are calculated as shown below: where y i indicates the vertical coordinate of the current pixel point, y u indicates the upper boundary of the overlapping area, y d indicates the lower boundary of the overlapping area, and y d = y u .

Results
In this section, we present a comprehensive evaluation of our method by performing both qualitative and quantitative comparisons with other stitching algorithms.Our algorithm simulations were implemented using Python.We commenced with a qualitative assessment, comparing our improved SIFT-based underwater image stitching technique with several established stitching algorithms.Specifically, we focused on feature matching performance when contrasting our method with the traditional SIFT approach.Additionally, we employed the UIQM metric to assess the performance of our method and of other stitching algorithms in stitching underwater images.To further substantiate the stitching performance, we conducted experimental assessments involving multiple underwater images.

Qualitative Comparisons
In this paper, we conducted a comparative analysis by selecting several representative image stitching algorithms.Specifically, we categorized the images into two groups (underwater images and land images) based on their content.Subsequently, we compared our method with four other approaches: traditional SIFT, APAP, SPHP, and AANAP.To facilitate alignment comparison across these stitching methods, we employed a straightforward pixel-weighted average fusion technique.To maintain conciseness, we present the results for the Temple [11] and Rock datasets only.Figures 3 and 4 depict the stitching results for the Rock and Temple images, respectively.The Rock dataset was obtained from the internet, while the Temple dataset was provided by the authors of [18].As observed in the figures, the results of the other methods exhibited varying degrees of ghosting artifacts within the overlapping regions.In contrast, our method demonstrated superior performance in eliminating artifacts from these overlapping areas.If we focus our attention on the regions enclosed by the red ellipses in the figures, it becomes evident that our results offer a more natural appearance.
Appl.Sci.2023, 13, x FOR PEER REVIEW 14 of 24 the internet, while the Temple dataset was provided by the authors of [18].As observed in the figures, the results of the other methods exhibited varying degrees of ghosting artifacts within the overlapping regions.In contrast, our method demonstrated superior performance in eliminating artifacts from these overlapping areas.If we focus our attention on the regions enclosed by the red ellipses in the figures, it becomes evident that our results offer a more natural appearance.Figure 3 illustrates the results for the Rock dataset, with each row representin output of a different method.The sequence of results is as follows: traditional SIFT rithm, APAP, AANAP, SPHP, and our method.The images' intricate textural fea particularly in underwater settings, present a substantial challenge for image stitc We highlight two aspects for each result: red boxes indicate parallax errors in the ove ping regions where the feature textures are clear, and blue boxes indicate parallax e in the overlapping regions where the feature textures are blurred.The outcomes o traditional SIFT method exhibit low alignment accuracy, with ghosting effects appe at the edges of the rocks, leading to a significant loss of detail information.
The results in the second row indicate that APAP performs well in regions with tively rich texture features but faces challenges in aligning areas with blurred de Moving to the third row, we can observe the outcomes of the AANAP method, w exhibit noticeable splicing artifacts and alignment errors at the junction of overlap and non-overlapping regions.In the fourth row, we present the results of the S method.SPHP demonstrates effective alignment in regions with abundant texture tures but encounters difficulties in regions characterized by blurred details.
The results of the Temple dataset are presented in Figure 4, with each row c sponding to a different method.The order of the results is as follows: the traditional algorithm, APAP, AANAP, SPHP, and our method.Two aspects of each result have highlighted, with blue boxes indicating parallax errors in the overlapping region an Figure 3 illustrates the results for the Rock dataset, with each row representing the output of a different method.The sequence of results is as follows: traditional SIFT algorithm, APAP, AANAP, SPHP, and our method.The images' intricate textural features, particularly in underwater settings, present a substantial challenge for image stitching.We highlight two aspects for each result: red boxes indicate parallax errors in the overlapping regions where the feature textures are clear, and blue boxes indicate parallax errors in the overlapping regions where the feature textures are blurred.The outcomes of the traditional SIFT method exhibit low alignment accuracy, with ghosting effects appearing at the edges of the rocks, leading to a significant loss of detail information.
The results in the second row indicate that APAP performs well in regions with relatively rich texture features but faces challenges in aligning areas with blurred details.Moving to the third row, we can observe the outcomes of the AANAP method, which exhibit noticeable splicing artifacts and alignment errors at the junction of overlapping and non-overlapping regions.In the fourth row, we present the results of the SPHP method.SPHP demonstrates effective alignment in regions with abundant texture features but encounters difficulties in regions characterized by blurred details.
The results of the Temple dataset are presented in Figure 4, with each row corresponding to a different method.The order of the results is as follows: the traditional SIFT algorithm, APAP, AANAP, SPHP, and our method.Two aspects of each result have been high-lighted, with blue boxes indicating parallax errors in the overlapping region and red boxes showing alignment errors at the junction of the overlapping and non-overlapping regions.
The traditional SIFT algorithm exhibits poor alignment in the overlapping areas, as evident in both the stone lion and roof areas, which suffer from ghosting issues.The second row displays the results of the APAP algorithm, which effectively aligns the overlapping areas but still exhibits some ghosting problems.In the third row, the results of AANAP reveal a parallax error in the overlapping area.In addition to ghosting, errors in alignment are noticeable, particularly in the houses at the overlapping edges of the two images.The fourth row illustrates the outcomes of the SPHP method.As described in the SPHP paper, this method prioritizes reducing perspective distortion rather than alignment accuracy.Parallax errors are evident in the overlapping area, particularly noticeable in the angles of the power poles, antennas, and houses.Our method demonstrates remarkable robustness when compared to the traditional SIFT algorithm, especially in scenarios where the textures surrounding the features exhibit high similarity.Figure 5 presents a visual comparison of the matching outcomes between our method and the conventional SIFT algorithm.The colored lines represent feature matching between two images, while the red boxes represent better stitching images obtained by our algorithm compared to traditional SIFT algorithms.Image A refers to the reference image and image B refers to the training image.As evident in the figure, the traditional SIFT algorithm struggled to obtain a satisfactory number of feature matches in the experimental image.Out of the 133 detected matching pairs, only 73 were accurate, resulting in an accuracy rate of 54.89%.This inaccuracy leads to improper alignment and significant ghosting artifacts in the traditional SIFT algorithm's results.Similar to the traditional SIFT algorithm, the outcomes of other stitching methods are also suboptimal.In contrast, our algorithm utilizes the projected transformed image for precise feature extraction and matching.This approach resulted in the identification of 85 correct matching pairs out of 139 detected, corresponding to an accuracy rate of 61.15%.Consequently, our algorithm excels in detecting a greater number of feature pairs and achieving higher matching accuracy, even in cases where feature detection is challenging in underwater images.The stitching results further underscore the superior alignment performance of our algorithm.
Figure 6 illustrates the comparative matching results between our method and the conventional SIFT algorithm using a different set of images.The colored lines represent feature matching between two images, while the red boxes represent better stitching images obtained by our algorithm compared to traditional SIFT algorithms.Image A refers to the reference image and image B refers to the training image.As depicted in Figure 6, the traditional SIFT algorithm successfully identified 277 correct matches out of 430 pairs, achieving an accuracy rate of 64.42%.In contrast, our algorithm identified 304 correct matches out of 434 pairs, yielding an accuracy rate of 70.05%.The results highlighted within the red box emphasize the limitations of the traditional SIFT algorithm, which excels in matching regions with distinct texture features but falters in regions with less prominent textures.Conversely, our algorithm exhibits the capability to accurately detect feature matches in regions with both prominent feature points and areas where the image contours are less distinct, leading to an enhanced quantity and accuracy of matches.This experiment underscores the heightened robustness and improved matching performance of our approach, stemming from our precise feature extraction and matching methodology.

Comparisons with Image Stitching Algorithms
UIQM is grounded in the framework proposed by Karen Panetta in 2015 [35], which takes into consideration three vital metrics for assessing underwater image quality: colorfulness measure (UICM), sharpness index (UISM), and contrast measure (UIConM).UIQM is computed as a linear combination of these three metrics.The higher the UIQM value, the more favorable the color balance, sharpness, and contrast of the image.UIQM serves as a comprehensive standard for evaluating the quality of underwater images, with a specific focus on their degradation characteristics.Widely recognized and authoritative in the realm of underwater image quality assessment, UIQM played a pivotal role in our study.Consequently, we conducted ten experiments on an image dataset and employed the UIQM metric to compare the performance of our proposed stitching method with that of several other stitching algorithms.
The formula for calculating the UIQM is shown below: Based on the results of [35], we took the scale factors as c 1 = 0.0282, c 2 = 0.2953, and c 3 = 3.5753.The results are shown in Table 2. Generally speaking, the value range of UIQM was between 0 and 3.For images with good imaging quality and clarity, the value of UIQM was generally above 1.5.As shown in Table 2, except for image 2, the UIQM value was relatively low due to the blurriness and poor quality of the image itself.Our method showed significant improvements in UIQM metrics compared to the traditional SIFT algorithm, with all datasets having higher scores and an average improvement of 6.04%.This demonstrates that our algorithm enhances the quality and sharpness of underwater panoramic images.Additionally, our method outperformed other algorithms such as APAP, AANAP, and SPHP, with improvements in UIQM of 4.73%, 9.01%, and 11.06%, respectively.
Appl.Sci.2023, 13, x FOR PEER REVIEW 18 of 24 Figure 6 illustrates the comparative matching results between our method and the conventional SIFT algorithm using a different set of images.The colored lines represent feature matching between two images, while the red boxes represent better stitching images obtained by our algorithm compared to traditional SIFT algorithms.Image A refers to the reference image and image B refers to the training image.As depicted in Figure 6, the traditional SIFT algorithm successfully identified 277 correct matches out of 430 pairs, achieving an accuracy rate of 64.42%.In contrast, our algorithm identified 304 correct matches out of 434 pairs, yielding an accuracy rate of 70.05%.The results highlighted within the red box emphasize the limitations of the traditional SIFT algorithm, which excels in matching regions with distinct texture features but falters in regions with less prominent textures.Conversely, our algorithm exhibits the capability to accurately detect feature matches in regions with both prominent feature points and areas where the image contours are less distinct, leading to an enhanced quantity and accuracy of matches.This experiment underscores the heightened robustness and improved matching performance of our approach, stemming from our precise feature extraction and matching methodology.

Comparisons with Image Stitching Algorithms
UIQM is grounded in the framework proposed by Karen Panetta in 2015 [35], which takes into consideration three vital metrics for assessing underwater image quality: colorfulness measure (UICM), sharpness index (UISM), and contrast measure (UIConM).UIQM is computed as a linear combination of these three metrics.The higher the UIQM value, the more favorable the color balance, sharpness, and contrast of the image.UIQM serves as a comprehensive standard for evaluating the quality of underwater images, with a specific focus on their degradation characteristics.Widely recognized and authoritative in the realm of underwater image quality assessment, UIQM played a pivotal role in our study.Consequently, we conducted ten experiments on an image dataset and employed the UIQM metric to compare the performance of our proposed stitching method with that of several other stitching algorithms.
The formula for calculating the UIQM is shown below: Based on the results of [35], we took the scale factors as  2. Generally speaking, the value range of UIQM

Stitching of Multiple Images
To demonstrate the effectiveness of our method in stitching multiple images, we present a panorama example consisting of multiple input images.In order to achieve better visual effects, our proposed fade-in and fade-out fusion algorithm was employed to handle the stitching gaps during the process.Figure 7 showcases a large underwater shipwreck scene.As shown in Figure 7, when compared with SIFT, APAP, and AANAP, our approach accomplishes the seamless stitching of the underwater shipwreck images while preserving the integrity of the image content and presenting a natural panorama.Parallax errors and perspective distortions are imperceptible in the resulting image.Even in the presence of a weak feature texture and blurred surface contours of the shipwreck caused by numerous attachments, our method effectively tackles the demanding challenge of feature extraction and matching.This showcases the resilience and comprehensiveness of our approach in managing the fusion of multiple images.

Discussion
The effectiveness of image stitching heavily depends on the image quality and the number of matched features.The proposed method maintains alignment accuracy and shows robustness in all of the abovementioned challenging examples.It performs better in aligning relatively texture-rich regions and can achieve effective image stitching in regions with blurred details.
The traditional SIFT algorithm requires that the targets to be stitched are in or approximately in the same plane.When there are occlusion and angle differences, ghosting of the stitched image occurs.APAP uses a local homography matrix for the calculation, which improves the performance of the stitching compared with the traditional SIFT algorithm approach.However, there are still some alignment errors.Moreover, the local model of APAP is vulnerable to the number and distribution of matching features.AA-NAP reduces the perspective error based on APAP, but it brings certain alignment error problems.AANAP linearizes the homology in the non-overlapping regions and uses global similarity transformation with automatic estimation in the overlapping regions.Smooth interpolation is used for extrapolation between the overlapping and non-overlapping regions.In this process, it is able to achieve a slow change between the overlapping and non-overlapping regions.However, the compensation effect of smooth interpolation is limited when there is already a large alignment error between the overlapping and nonoverlapping regions.The current parameter selection of the SPHP algorithm does not take into account the image content (e.g., line features), which can lead to the algorithm failing to reduce the distortion of the line.

Conclusions
In this paper, we introduce an enhanced SIFT method designed for the stitching of underwater images characterized by smeared feature contours.Our method stands out As shown in Figure 7, when compared with SIFT, APAP, and AANAP, our approach accomplishes the seamless stitching of the underwater shipwreck images while preserving the integrity of the image content and presenting a natural panorama.Parallax errors and perspective distortions are imperceptible in the resulting image.Even in the presence of a weak feature texture and blurred surface contours of the shipwreck caused by numerous attachments, our method effectively tackles the demanding challenge of feature extraction and matching.This showcases the resilience and comprehensiveness of our approach in managing the fusion of multiple images.

Discussion
The effectiveness of image stitching heavily depends on the image quality and the number of matched features.The proposed method maintains alignment accuracy and shows robustness in all of the abovementioned challenging examples.It performs better in aligning relatively texture-rich regions and can achieve effective image stitching in regions with blurred details.
The traditional SIFT algorithm requires that the targets to be stitched are in or approximately in the same plane.When there are occlusion and angle differences, ghosting of the stitched image occurs.APAP uses a local homography matrix for the calculation, which improves the performance of the stitching compared with the traditional SIFT algorithm approach.However, there are still some alignment errors.Moreover, the local model of APAP is vulnerable to the number and distribution of matching features.AANAP reduces the perspective error based on APAP, but it brings certain alignment error problems.AANAP linearizes the homology in the non-overlapping regions and uses global similarity transformation with automatic estimation in the overlapping regions.Smooth interpolation is used for extrapolation between the overlapping and non-overlapping regions.In this process, it is able to achieve a slow change between the overlapping and non-overlapping regions.However, the compensation effect of smooth interpolation is limited when there is already a large alignment error between the overlapping and non-overlapping regions.The current parameter selection of the SPHP algorithm does not take into account the image content (e.g., line features), which can lead to the algorithm failing to reduce the distortion of the line.

Conclusions
In this paper, we introduce an enhanced SIFT method designed for the stitching of underwater images characterized by smeared feature contours.Our method stands out because it does not rely on detector information for feature extraction.Even in regions with blurred texture features, where traditional feature extraction methods falter, our approach significantly augments the number of accurately matched pairs, thereby elevating the precision of feature matching.Our method incorporates two critical phases for the stitching process, each rooted in the SIFT algorithm: coarse feature extraction and precise feature extraction.To enhance the image contrast and quality, we integrate a dynamic threshold white balance algorithm and CLAHE.Additionally, we have refined the fade-in and fadeout technique for generating panoramic images.A series of experiments were conducted to validate the effectiveness of our algorithm.Our results underscore the superiority of our method over the conventional SIFT algorithm.It effectively mitigated ghosting issues in stitched images and enhanced the UIQM measure by an average of 6.04% across ten underwater image stitching experiments.When compared to alternative algorithms such as APAP, SPHP, and AANAP, our method outperformed them, boosting the UIQM index by 4.73%, 9.01%, and 11.06%, respectively.Furthermore, we demonstrated the efficacy of our method in seamlessly stitching multiple images, preserving the integrity of the image content in expansive scenes without significant ghosting artifacts.Our future research endeavors will concentrate on compensating for parallax in cases of substantial movements and addressing the influence of moving objects within the scene on stitching quality.

Figure 1 .
Figure 1.Various stages in the improved SIFT underwater image stitching technique.

Figure 1 .
Figure 1.Various stages in the improved SIFT underwater image stitching technique.

Figure 2 .
Figure 2. (a,b) The original underwater images; (c,d) the preprocessed underwater images generated by our method.

Figure 2 .
Figure 2. (a,b) The original underwater images; (c,d) the preprocessed underwater images generated by our method.

Figure 5 .
Figure 5.An image feature matching case: (a) Image inputs to be stitched.(b) Our feature matching and the traditional SIFT algorithm.(c) Stitching results of the traditional SIFT algorithm and ours.Figure 5.An image feature matching case: (a) Image inputs to be stitched.(b) Our feature matching and the traditional SIFT algorithm.(c) Stitching results of the traditional SIFT algorithm and ours.

Figure 5 .
Figure 5.An image feature matching case: (a) Image inputs to be stitched.(b) Our feature matching and the traditional SIFT algorithm.(c) Stitching results of the traditional SIFT algorithm and ours.Figure 5.An image feature matching case: (a) Image inputs to be stitched.(b) Our feature matching and the traditional SIFT algorithm.(c) Stitching results of the traditional SIFT algorithm and ours.

Figure 6 .
Figure 6.An image feature matching case: (a) Image inputs to be stitched.(b) The feature matching of image B. (c) Our feature matching and the traditional SIFT algorithm.

.
The results are shown in Table

Figure 6 .
Figure 6.An image feature matching case: (a) Image inputs to be stitched.(b) The feature matching of image B. (c) Our feature matching and the traditional SIFT algorithm.

Figure 7 .
Figure 7. Stitching results of multiple images in perspective projection on the Ship image pair: (a-c) Image inputs to be stitched.(d) SIFT.(e) APAP, As-Projective-As-Possible warps.(f) AANAP, adaptive As-Natural-As-Possible;.(g) Our method.

Figure 7 .
Figure 7. Stitching results of multiple images in perspective projection on the Ship image pair: (a-c) Image inputs to be stitched.(d) SIFT.(e) APAP, As-Projective-As-Possible warps.(f) AANAP, adaptive As-Natural-As-Possible;.(g) Our method.

Table 1 .
Underwater image quality evaluations.

Table 1 .
Underwater image quality evaluations.

Table 2 .
UIQM index of different stitching methods.