Overcoming Periodic Stripe Noise in Infrared Linear Array Images: The Fourier-Assisted Correlative Denoising Method

Infrared linear array detectors frequently experience vertical, low-frequency, and periodic stripe noise during imaging, stemming from electro-mechanical interference. Unlike conventional periodic disturbances, this interference showcases long periodicities and is uniquely columnar in orientation. Its presence, especially within the low-frequency domain, renders conventional filtering techniques ineffective and, at times, detrimental to image quality. Addressing this challenge, we introduce Fourier-Assisted Correlative Denoising (FACD), a correlation-centric denoising approach tailored for such unique interference patterns. This mechanism begins with the capture of a pure background image, inclusive of periodic noise, during the non-uniform correction phase of the infrared detector. Leveraging the noise’s frequency domain attributes, we extract a one-dimensional single-cycle noise signal. The infrared image is subsequently segmented into parts, and using the detected noise periodicity, the one-dimensional signals for each segment are computed. By leveraging the correlation between these signals and the benchmark one-dimensional noise pattern, we ascertain the noise profile within each segment. This profile is then employed for spatial domain denoising across the entire image frame. Empirical assessments confirm that the FACD outperforms contemporary denoising techniques by augmenting the peak signal-to-noise ratio by approximately 2.5 dB, underscoring its superior robustness. Furthermore, in light of its specificity to this noise model, FACD rapidly denoises high-resolution real infrared linear array scans, thus meeting the stringent real-time and resolution imperatives of advanced infrared linear array scanning apparatuses.


Introduction
As advancements continue in large-scale focal plane arrays and infrared linear array detector techniques [1,2], the 360-degree rotational scanning technique has increasingly caught the eye for its applicability to wide-angle imaging scenarios.Infrared linear array detectors have been effectively deployed in diverse fields, encompassing aerial target surveillance [3][4][5], medical diagnostics [6], facial recognition [7], industrial inspections [8], and aerospace applications [9].Parallel to these developments, there has been a burgeoning interest in the research community toward refining denoising techniques for images obtained through infrared linear array detectors [10,11].In imagery, disturbances related to infrared linear arrays primarily fall into stripe noises [12,13] and high-frequency mixing noises.The former can be further divided into those resulting from detector non-uniformities [14] and those arising from electromagnetic interferences with the detector [15].
Stripe noise, originating from detector inconsistencies, frequently demonstrates a directional preference, manifesting horizontally or vertically.This noise distribution does not adhere to a consistent pattern and typically corresponds to the detector's photosensitive element's response to illumination.Addressing the challenges posed by image stripes Sensors 2023, 23, 8716 2 of 20 due to detector non-uniformity, Li Fangzhou et al. [11] proposed an enhanced NL-means model capable of concurrently eliminating stripe noise and Gaussian noise within images.In similar endeavors, Wang E et al. [16,17] amalgamated guided filtering with wavelet transformations, aiming to extract stripe noise from images within the frequency domain.Periodic stripe noise, attributed to electromagnetic disturbances, frequently presents with short intervals and high frequencies, while its direction generally does not appear horizontal or vertical.It is challenging to denoise this noise directly from the spatial domain.A prevailing method encompasses transitioning the image to the frequency domain, facilitating more efficient noise mitigation.In this vein, Hamd M H et al. [18] proposed a filter anchored in the noise signal's peak value, designed to target and neutralize the periodic noise crest within the frequency domain.In a parallel approach, Yadav V P et al. [19] outlined a denoising methodology predicated on a localized threshold.
Infrared linear array detectors, when employed in 360-degree rotational scanning, require instantaneous imaging due to their rapid acquisition rates.A single pass can produce tens of thousands of image columns.This imaging modality introduces a peculiar kind of periodic stripe noise attributable to electromagnetic interference.Unlike noises from other devices, this noise showcases a protracted period and a lower frequency, aligning conspicuously with the direct current (DC) component within the frequency domain.
Notably, deploying denoising methodologies as suggested by Hamd M H et al. [18] and Yadav V P et al. [19], which primarily target peak noise in the frequency domain, might compromise the image's low-frequency constituents, leading to marked image degradation.Additionally, the periodic noise in infrared linear array rotational scanning imagery is elongated, with gradual pixel transitions across adjacent columns, lacking any distinct singular stripes.As a result, denoising strategies, as posited by Li Fangzhou et al. [11] and Wang E et al. [16,17] encounter reduced efficacy.The enhanced NL-means paradigm presented by Li Fangzhou et al. [11] embeds a sliding search operation.Simultaneously, the methodologies proposed by Wang E et al. [16,17] intertwine guided filtering, the MHE algorithm, and TV-L1 regularization-factors that can potentially constrain the achievement of the real-time imaging prerequisites of infrared linear array systems.With the advancements in deep learning, denoising algorithms based on this framework have shown significant effectiveness [20,21].However, these methodologies' substantial dependency on curated image datasets is a significant limitation.The intrinsically high resolution of infrared linear array images complicates the task of compiling a comprehensive dataset, presenting formidable challenges.The self-supervised denoising method pioneered by Jaakko Lehtinen et al. [22], along with the approach by Qingliang Jiao et al. [23] designed to minimize reliance on pristine training data are particularly significant.While both methodologies cater to application requirements, the neural network model's consistent input size necessitates size preprocessing for the THz spectrum, amplifying the complexity.Furthermore, the substantial computational demand of these techniques poses challenges for real-time processing of high-resolution images.
Considering the abovementioned challenges, this paper proposes an infrared linear array image denoising algorithm based on signal correlation.The foundational framework of the denoising algorithm hinges on a correlation-centric approach.Specifically, our proposed algorithm, FACD (Fourier-Assisted Correlative Denoising), operates by leveraging the intrinsic frequency domain attributes of periodic noise.At its core, FACD begins by capturing a background image saturated with the said noise.It then delves into the frequency spectrum to extract a one-dimensional, single-cycle noise signal.By segmenting the infrared image and matching the detected noise periodicity, the algorithm harnesses the correlation between these segmented signals and the established one-dimensional noise pattern.This correlation-driven methodology enables FACD to precisely ascertain the noise profile within each segment of the image, facilitating effective spatial domain denoising.The primary advantage of the FACD approach lies in its precision and specificity to target the unique interference patterns presented by electromechanical interferences in infrared linear array detectors.Experimental results show this algorithm proficiently eliminates Sensors 2023, 23, 8716 3 of 20 periodic noise in high-resolution 360-degree rotational scanning images, satisfying the real-time imperatives of infrared linear array detectors.

Characteristics of the Periodic Noise
Periodic noise is a prevalent noise type in images, typically emanating from interference during image capture or subsequent transmission.Traditionally, this noise is typified by its high frequency, manifesting as luminous sectors within the image's two-dimensional Fourier transform domain, as delineated in Figure 1.An example would be a newspaper image sampled under the Moiré pattern.spatial domain denoising.The primary advantage of the FACD approach lie precision and specificity to target the unique interference patterns presen electromechanical interferences in infrared linear array detectors.Experimental show this algorithm proficiently eliminates periodic noise in high-resolution 360 rotational scanning images, satisfying the real-time imperatives of infrared linea detectors.

Characteristics of the Periodic Noise
Periodic noise is a prevalent noise type in images, typically emanatin interference during image capture or subsequent transmission.Traditionally, this typified by its high frequency, manifesting as luminous sectors within the image dimensional Fourier transform domain, as delineated in Figure 1.An example w a newspaper image sampled under the Moiré pattern.This paper concentrates on a distinct type of noise, divergent from trad periodic disturbances.Infrared linear array detectors capture the images under s during their comprehensive 360-degree rotational scanning.Electromec disturbances during this procedure introduce noise with a lower frequency.In the domain of the image, this noise manifests as subtle, gradually varying stripes.Conv within the image's two-dimensional Fourier transform domain, it is observed as lu sectors proximate to the DC component.The noise image of the infrared detector o through blackbody is shown in Figure 2a, and its two-dimensional Fourier fre domain image is presented in Figure 2b.Upon juxtaposing the frequency domain images of periodic noise depicted in 1 and 2, one discerns a clear contrast.Traditional periodic noise, when viewed frequency domain, is characterized by multiple bright regions distanced from This paper concentrates on a distinct type of noise, divergent from traditional periodic disturbances.Infrared linear array detectors capture the images under scrutiny during their comprehensive 360-degree rotational scanning.Electromechanical disturbances during this procedure introduce noise with a lower frequency.In the spatial domain of the image, this noise manifests as subtle, gradually varying stripes.Conversely, within the image's two-dimensional Fourier transform domain, it is observed as luminous sectors proximate to the DC component.The noise image of the infrared detector obtained through blackbody is shown in Figure 2a, and its two-dimensional Fourier frequency domain image is presented in Figure 2b.
electromechanical interferences in infrared linear array detectors.Experimental show this algorithm proficiently eliminates periodic noise in high-resolution 360rotational scanning images, satisfying the real-time imperatives of infrared linea detectors.

Characteristics of the Periodic Noise
Periodic noise is a prevalent noise type in images, typically emanating interference during image capture or subsequent transmission.Traditionally, this n typified by its high frequency, manifesting as luminous sectors within the image dimensional Fourier transform domain, as delineated in Figure 1.An example wo a newspaper image sampled under the Moiré pattern.This paper concentrates on a distinct type of noise, divergent from trad periodic disturbances.Infrared linear array detectors capture the images under sc during their comprehensive 360-degree rotational scanning.Electromech disturbances during this procedure introduce noise with a lower frequency.In the domain of the image, this noise manifests as subtle, gradually varying stripes.Conv within the image's two-dimensional Fourier transform domain, it is observed as lum sectors proximate to the DC component.The noise image of the infrared detector ob through blackbody is shown in Figure 2a, and its two-dimensional Fourier freq domain image is presented in Figure 2b.Upon juxtaposing the frequency domain images of periodic noise depicted in F 1 and 2, one discerns a clear contrast.Traditional periodic noise, when viewed frequency domain, is characterized by multiple bright regions distanced from t Employing a standard Gaussian notch filter to address this periodic noise can unintentionally weaken the image's inherent low-frequency elements, negatively impacting the denoising outcome.To better address the low-frequency periodic noise for denoising, analyzing the infrared noise image and introducing some prior information is essential.Examining the pure background noise image sourced from the blackbody, it is evident that the noise primarily appears as vertical stripes with varying brightness.Consequently, a comprehensive statistical analysis of the noise is conducted in both horizontal and vertical orientations, as delineated in Figure 3.To better address the low-frequency periodic noise for denoising, analyzing the infrared noise image and introducing some prior information is essential.Examining the pure background noise image sourced from the blackbody, it is evident that the noise primarily appears as vertical stripes with varying brightness.Consequently, a comprehensive statistical analysis of the noise is conducted in both horizontal and vertical orientations, as delineated in Figure 3.As illustrated in Figure 3, the noise distinctly displays a periodic characteristic along the image's columnar axis, whereas along the row axis, it resembles stochastic Gaussian noise.Comprehensive assessments using the blackbody substantiate that the amplitude of this periodic noise remains unaffected by variations in the blackbody's temperature.This inference posits that the periodic noise prevalent in infrared detectors originates exclusively from electromechanical interference, devoid of influence from other external variables.Every time the infrared linear array scanner is powered on, it necessitates a nonuniformity adjustment, obtaining the images of periodic noise under a pure background.Subsequently, the noise period could be determined, and denoising could be performed in the spatial domain based on signal correlation.

Noise Period Calculation Based on Fourier Transform
Normally, infrared detectors need to be calibrated to ensure consistent response from the detection unit before each use.In this study, the non-uniformity correction deployed for the infrared detector is based on the single-point method.Given that the infrared images acquired during this correction phase primarily reflect the pure background, they remain relatively insulated from external scene perturbations.Such images exhibit a pronounced presence of periodic noise.Therefore, the algorithm employs the Fourier transform to analyze the images obtained during non-uniformity correction, aiming to discern the noise periodicity attributable to electromechanical interference.In this example, upon examination of the periodic noise within the pure background imagery, it was ascertained that the noise periodicity for the selected infrared detector approximates 850 columns.Furthermore, recognizing the high resolution inherent to infrared linear array scanning images and the need for multiple noise periods in Fourier transform As illustrated in Figure 3, the noise distinctly displays a periodic characteristic along the image's columnar axis, whereas along the row axis, it resembles stochastic Gaussian noise.Comprehensive assessments using the blackbody substantiate that the amplitude of this periodic noise remains unaffected by variations in the blackbody's temperature.This inference posits that the periodic noise prevalent in infrared detectors originates exclusively from electromechanical interference, devoid of influence from other external variables.Every time the infrared linear array scanner is powered on, it necessitates a non-uniformity adjustment, obtaining the images of periodic noise under a pure background.Subsequently, the noise period could be determined, and denoising could be performed in the spatial domain based on signal correlation.

Noise Period Calculation Based on Fourier Transform
Normally, infrared detectors need to be calibrated to ensure consistent response from the detection unit before each use.In this study, the non-uniformity correction deployed for the infrared detector is based on the single-point method.Given that the infrared images acquired during this correction phase primarily reflect the pure background, they remain relatively insulated from external scene perturbations.Such images exhibit a pronounced presence of periodic noise.Therefore, the algorithm employs the Fourier transform to analyze the images obtained during non-uniformity correction, aiming to discern the noise periodicity attributable to electromechanical interference.In this example, upon examination of the periodic noise within the pure background imagery, it was ascertained that the noise periodicity for the selected infrared detector approximates 850 columns.Furthermore, recognizing the high resolution inherent to infrared linear array scanning images and the need for multiple noise periods in Fourier transform calculations, this section crops the line scanning images, adjusting their resolution to 1024 × 4096.
Since the periodic noise of infrared linear array images primarily shows in the column direction, the clipped image is first converted into a one-dimensional signal in the column direction, as shown by the following formula: where I c is the image used for non-uniform correction, i, the number of columns, obtains i ∈ [1, N col ], where N col = 4096; N row = 1024 is the number of rows of the image.Finally, I c1 is the transformed one-dimensional signal sequence, as shown in Figure 3a.
In the next step, Fourier transform is applied to I c , as shown in the following formula: In this formula, X(m) is the first half of the sequence after the Fourier transform, excluding the DC component, where m ∈ 1, 1 2 N col .Using the modulus values of X(m), the maximum value of X(m) is calculated, then record the current subscript m MAX is recorded, as shown in the following equation: In this context, the function max(•) retrieves the maximum value, while abs(•) performs the modulus operation on complex numbers.The noise period's estimated value, denoted as P e , is elucidated in the subsequent equation: Since the number of columns of the cropped image N col may not necessarily be an integer multiple of the noise period, P e is the estimated value of the noise signal period.The closer the number of columns of the cropped image is to an integer multiple of the noise period and the larger that number, the closer P e is to the actual noise period.
In addition, other noises that are present in the collected signal can also influence the estimated value P e of the electromechanical noise period.Using P e as a constraining parameter and leveraging the periodicity of I c1 , our objective is to accurately determine the noise period, denoted P.
Initially, the median value, 'Mid', of the sequence I c1 is calculated as indicated in the formula below: In this context, min(•) represents the function used to determine the minimum value, while max(•) signifies the function to ascertain the maximum value.Subsequently, we find points in the sequence I c1 that are near 'Mid', as demonstrated in the following equation: Within this framework, P array denotes a set of points in I c1 proximate to the median 'Mid'.The function Find(•) is used to identify values in I c1 that satisfy designated conditions.In this equation, V 1 acts as a range threshold parameter.Given that the infrared image of this device is articulated in 8-bit grayscale, this threshold is pragmatically set at 2. The objective is to discern values in the I c1 series that fluctuate between the [Mid − 2, Mid + 2] range.With the noise period estimate serving as a benchmark restraint, the precise value P's computation is illustrated in the subsequent equation: T max = mean P array Find P array > P array (1) + P e − V 2 & P array < P array (1) T min = mean P array Find P array > P array (1) − V 2 & P array < P array (1) + V 2 (8) In this discussion, P represents the definitive noise period value, while P array (1) extracts the first value from the P array series.The mean(•) function procures the average value, and the Find(•) function is employed to pinpoint values within P array that resonate with the pre-established criteria.As evident from Figure 3, which delineates the noise image's columnar orientation statistical depiction, the electromechanical noise periodicity is roughly around 800 columns.The range threshold V 2 in the formula is manually set to 25 with the intention to constrain the range for P so that P ∈ [700, 900], ensuring its proximity to the actual noise period.

Fourier-Assisted Correlative Denoising (FACD) Mechanism
The periodic noise in the infrared columnar image caused by electromechanical interference is a type of additive noise, and its image signal model is shown in Equation (10).
In this context, v i represents the original infrared imagery filled with noise disruptions, while u i signifies the ideal noise-free infrared imagery.s i delineates the periodic noise due to electromechanical disturbances, and t i captures other noise modalities.The aim of the algorithm in this subsection is to utilize the correlation of periodic signals to eliminate the periodic additive noise on the infrared image.As observed from Section 2.2, even with our computational derivation of the noise image's periodic magnitude P via Fourier transformation, direct denoising of this noise image remains challenging.This challenge arises when the linear array detector acquires images in a scanning fashion.Even if multiple images represent the same scene, the moments of capture might vary, leading to potential differences in the phase offset of noise stripes, as depicted in Figure 4.

𝑃 = 𝑇 − 𝑇
In this discussion,  represents the definitive noise period value, while extracts the first value from the  series.The  (•) function proc average value, and the  (•) function is employed to pinpoint values with that resonate with the pre-established criteria.As evident from Figure 3, which d the noise image's columnar orientation statistical depiction, the electromechani periodicity is roughly around 800 columns.The range threshold  in the fo manually set to 25 with the intention to constrain the range for  so that  ∈ [7 ensuring its proximity to the actual noise period.

Fourier-Assisted Correlative Denoising (FACD) Mechanism
The periodic noise in the infrared columnar image caused by electrom interference is a type of additive noise, and its image signal model is shown in (10).

𝑣 = 𝑢 + 𝑠 + 𝑡
In this context,  represents the original infrared imagery filled wi disruptions, while  signifies the ideal noise-free infrared imagery. delin periodic noise due to electromechanical disturbances, and  captures oth modalities.The aim of the algorithm in this subsection is to utilize the corre periodic signals to eliminate the periodic additive noise on the infrared im observed from Section 2.2, even with our computational derivation of the noise periodic magnitude  via Fourier transformation, direct denoising of this noi remains challenging.This challenge arises when the linear array detector acquire in a scanning fashion.Even if multiple images represent the same scene, the mo capture might vary, leading to potential differences in the phase offset of noise s depicted in Figure 4.   Based on the analysis of the linear array detector's periodic noise in Section 2.1, the periodic noise possesses a robust one-dimensional periodic characteristic in the column direction.It is feasible to exploit the correlation in the noise signal's columnar direction to Sensors 2023, 23, 8716 7 of 20 determine the phase of the periodic noise in the image.The formulas for auto-correlation and cross-correlation of the signal are provided in the succeeding equations.
Let x(n) and y(n) be one-dimensional sequence signals with ∈ [1, N], where N sig- nifies the length of the sequence and sign(•) is the signum function.During the nonuniformity correction stage of the detector, a periodic noise image is secured against a pure background, as exemplified in Figure 2a.The average value for each column of the image is determined as outlined in Equation (13).
In this context, n(j) refers to the one-dimensional sequence depicted in Figure 3a, which is recognized as the multi-cycle noise present in the image's columnar orientation.The term I n (i, j) signifies a pixel within the noise image when set against a pure background.Here, with ∈ where C represents the total number of columns in the image and R indicates its rows.From the sequence n(j), we extract a contiguous sequence with length paralleling the noise cycle P, termed n(j), where j is within the interval j ∈ 1, P .
To address the infrared image affected by periodic noise, it is segmented into individual blocks.Each block has a column size matching the noise cycle P and row size approximating 1/4 of the complete frame's row count.This modular decomposition is graphically represented in Figure 5.To address the infrared image affected by periodic noise, it is se individual blocks.Each block has a column size matching the noise cycle  approximating 1/4 of the complete frame's row count.This modular dec graphically represented in Figure 5.The subsequent phase involves computing the one-dimensional sequence in the column direction for each subdivided block, as represented by the following equation:

𝑅 (𝑙
In this equation: 1. R denotes the row count of the image I(i, j);

2.
n stands for the count of segmentation blocks in the row direction and n ∈ [0, 3]; 3.
m stands for the count of segmentation blocks in the column direction; m ∈ 0, f loor C/ P − 1 , where the function f loor(•) rounds a number down to its nearest integer, C signifies the column count of the image.
Use Formula ( 14) to calculate the offset m p at which the cross-correlation between each small block's one-dimensional column direction sequence I p (j) and n (j) reaches its maximum.At this time, the length of the sequence N = P, and m p is transformed as shown in the following formula: Sort all the m p values to obtain the sequence m p (i Here, N m is the number of small image blocks as shown in Figure 5.Take the middle portion of the sorted m p sequence offsets and compute their mean to obtain the corresponding one-dimensional noise sequence for the small image block, as shown in the following formula: where mean(•) is the average function and k ∈ 1, P .Use n p (k) to denoise the image I(i, j), as shown in the equation: I(i, j) − n p mod j, P , mod j, P = 0 where I u (i, j) is the image after denoising using correlation, mod(•) is the modulo function, and P is the noise period obtained using the Fourier Transform in Section 2.2.

Overall Algorithmic Procedure
As illustrated in Figure 6, the comprehensive algorithmic procedure is delineated systematically as follows: • Step 1: Single-Cycle Noise Sequence Extraction n (j) 1.
During the non-uniformity correction process of the infrared camera, collect the cyclical noise image under a pure background; 2.
Implementing the Fourier Transform-centric periodic solution algorithm, as introduced in Section 2.2, to identify the noise period P for the observed noise signal; 3.
Employing Formula (13) and the ascertained noise period P, derive the singlecycle noise sequence n (j).

•
Step 2: Image Data Partition and Processing Begin by subdividing the procured images into specific blocks, as illustrated in Figure 5. Subsequently, for each image block, compute the one-dimensional sequence I p (j) utilizing Formula (14).

•
Step 3: Cross-Correlation Function and Offset Determination For each one-dimensional image sequence I p (j), its cross-correlation function with the noise sequence n (j) is calculated using Formula (12).Record the offset m p when each cross-correlation function reaches its maximum value.Thereafter, a transformative Equation ( 15) is applied to the offset m p .

•
Step 4: Sequence Ordering and Offset Refinement After the previous steps, the offset series m p is systematically arranged to form the sequence m p (i), i ∈ [1, N m ].Herein, N m represents the total number of image blocks, as Sensors 2023, 23, 8716 9 of 20 explicated in Figure 5. Eliminate both the top quartile of the larger sequence values and the bottom quartile of the smaller ones.Following this, Equation ( 16) is employed to produce the one-dimensional noise sequence n p (j), constrained by j ∈ 1, P .

•
Step 5: Image Denoising Collaboratively, with the sequence n p (j), Equation ( 17) is leveraged to execute a denoising process on the image.
Sensors 2023, 23, 8716 9 of 21 explicated in Figure 5. Eliminate both the top quartile of the larger sequence values and the bottom quartile of the smaller ones.Following this, Equation ( 16) is employed to produce the one-dimensional noise sequence  (), constrained by  ∈ 1,  .

•
Step 5: Image Denoising Collaboratively, with the sequence  () , Equation ( 17) is leveraged to execute a denoising process on the image.

Investigation into the Consistency of Periodic Noise in Infrared Linear Array Detectors
The preceding section articulated that the Fourier-Assisted Correlative Denoising (FACD) mechanism requires an image characterized by periodic noise against a pure background.We conducted a comprehensive experiment to assess the applicability and robustness of the proposed mechanism in real-world scenarios.For the empirical study, the experiment was designed to meticulously record and analyze the fluctuations in noise periodicity and amplitude over successive time intervals following the initial power-up of the detector without undergoing any rebooting processes.
We employed a 1024-pixel long-wave infrared linear array detector, which boasts a specialized shutter mechanism.This shutter is purposefully activated when a nonuniformity correction command is initiated, facilitating the precise non-uniformity correction procedure.The images obtained at this critical juncture form the foundation of our data source.By applying histogram statistical methods, we discerned the columnar distribution associated with the noise periodicity within the current image framework.

Investigation into the Consistency of Periodic Noise in Infrared Linear Array Detectors
The preceding section articulated that the Fourier-Assisted Correlative Denoising (FACD) mechanism requires an image characterized by periodic noise against a pure background.We conducted a comprehensive experiment to assess the applicability and robustness of the proposed mechanism in real-world scenarios.For the empirical study, the experiment was designed to meticulously record and analyze the fluctuations in noise periodicity and amplitude over successive time intervals following the initial power-up of the detector without undergoing any rebooting processes.
We employed a 1024-pixel long-wave infrared linear array detector, which boasts a specialized shutter mechanism.This shutter is purposefully activated when a non-uniformity correction command is initiated, facilitating the precise non-uniformity correction procedure.The images obtained at this critical juncture form the foundation of our data source.By applying histogram statistical methods, we discerned the columnar distribution associated with the noise periodicity within the current image framework.Data were harvested across varied operational intervals, culminating in four distinct datasets.Each of these sets encapsulated an average of 55 noise cycles.The resultant histograms are illustrated in the subsequent figure.Within this graphical representation, the x-axis demarcates the column count corresponding to the noise periodicity, while the y-axis quantifies the cycle counts spanning diverse noise periods.
Figure 7 provides a sequence of histograms documenting the performance of the infrared detector at different intervals.Specifically, Figure 7a showcases the data shortly after the initial power-on and subsequent cooling.The following histograms, Figure 7b-d, were captured at 6 h intervals, marking 6, 12, and 18 h of continuous operation.A consistent observation across these figures is the stable noise periodicity, consistently measuring 847 columns.
Sensors 2023, 23, 8716 10 of 21 Data were harvested across varied operational intervals, culminating in four distinct datasets.Each of these sets encapsulated an average of 55 noise cycles.The resultant histograms are illustrated in the subsequent figure.Within this graphical representation, the x-axis demarcates the column count corresponding to the noise periodicity, while the y-axis quantifies the cycle counts spanning diverse noise periods.
Figure 7 provides a sequence of histograms documenting the performance of the infrared detector at different intervals.Specifically, Figure 7a showcases the data shortly after the initial power-on and subsequent cooling.The following histograms, Figures 7bd, were captured at 6 h intervals, marking 6, 12, and 18 h of continuous operation.A consistent observation across these figures is the stable noise periodicity, consistently measuring 847 columns.For a deeper insight into noise amplitude fluctuations, we compiled and analyzed data from 220 noise periods across the four intervals.Figure 8 presents this analysis, where the x-axis demarcates experimental batches of 22 periodic signals and the y-axis displays amplitude variations in pixel values.The noise amplitude across all periods demonstrated remarkable stability, fluctuating minimally between −30 and 30 pixels.Drawing from these observations, the infrared linear array detector for the experiment, once powered up, maintains consistent noise characteristics in terms of period and amplitude.This consistency underscores the effectiveness of the FACD algorithm.Initialized with a single-period noise signal extracted at power-up, the FACD algorithm can consistently denoise images without needing recalibration during extended operation periods.For a deeper insight into noise amplitude fluctuations, we compiled and analyzed data from 220 noise periods across the four intervals.Figure 8 presents this analysis, where the x-axis demarcates experimental batches of 22 periodic signals and the y-axis displays amplitude variations in pixel values.The noise amplitude across all periods demonstrated remarkable stability, fluctuating minimally between −30 and 30 pixels.
datasets.Each of these sets encapsulated an average of 55 noise cycles.histograms are illustrated in the subsequent figure.Within this graphical r the x-axis demarcates the column count corresponding to the noise periodi y-axis quantifies the cycle counts spanning diverse noise periods.
Figure 7 provides a sequence of histograms documenting the perfo infrared detector at different intervals.Specifically, Figure 7a showcases th after the initial power-on and subsequent cooling.The following histogram d, were captured at 6 h intervals, marking 6, 12, and 18 h of continuous consistent observation across these figures is the stable noise periodicity measuring 847 columns.Drawing from these observations, the infrared linear array det experiment, once powered up, maintains consistent noise characteristic period and amplitude.This consistency underscores the effectiveness algorithm.Initialized with a single-period noise signal extracted at poweralgorithm can consistently denoise images without needing recalibration du Drawing from these observations, the infrared linear array detector for the experiment, once powered up, maintains consistent noise characteristics in terms of period and amplitude.This consistency underscores the effectiveness of the FACD algorithm.Initialized with a single-period noise signal extracted at power-up, the FACD algorithm can consistently denoise images without needing recalibration during extended operation periods.

Impact of Image Block Width
As depicted in Figure 5, this approach employs a strategy of partitioning the image into blocks to enhance the denoising results of the correlation-based algorithm.The length of each image block corresponds to the signal period deduced from the algorithm in Section 2.2, while the width is manually set to a quarter of the entire frame's row count.During real-world capture by the camera, the impact of scene elements within the image on inter-signal correlation cannot be disregarded.When the width of an image block is set smaller, more blocks can be derived from an infrared linear array image of the exact resolution.Since fewer pixels exist within a smaller block, the calculation time for the one-dimensional column-direction signal of the block image is reduced.However, the signal within the block is largely determined by the scene elements, making the influence of the image scene on the correlation algorithm more prominent.
Conversely, when one opts for larger block widths, the number of blocks extracted from an infrared linear array image of the exact dimensions is reduced.The calculation time for the column-direction signal increases.Given that the correlation between the image scene and periodic noise is relatively weak, the influence of the image scene on the correlation algorithm is considerably reduced when the block's width is larger.
In this subsection, we conduct an experiment utilizing 20 infrared line-scanning images, each with a resolution of 1024 × 4096.A pure background periodic noise image of identical resolution is acquired through a blackbody.We create a simulated infrared image with noise by superimposing this periodic noise image onto the clean infrared line-scanning image.Using the algorithm described in Section 2.2 and the pure background noise image, we identify a signal period covering 847 columns.From these data, we derive a single-cycle, one-dimensional noise signal sequence.For the batch of 20 noisy infrared images, we randomly pick image segments sized 1024 × 847 and segment them further into blocks.The subsequent denoising process is executed using the approach described in Section 3.2.The relationship between the block width of the images and their corresponding denoising results is illustrated in Figure 9.

Impact of Image Block Width
As depicted in Figure 5, this approach employs a strategy of partitioni into blocks to enhance the denoising results of the correlation-based algorithm of each image block corresponds to the signal period deduced from the Section 2.2, while the width is manually set to a quarter of the entire frame' During real-world capture by the camera, the impact of scene elements with on inter-signal correlation cannot be disregarded.When the width of an imag smaller, more blocks can be derived from an infrared linear array image resolution.Since fewer pixels exist within a smaller block, the calculation tim dimensional column-direction signal of the block image is reduced.Howev within the block is largely determined by the scene elements, making the infl image scene on the correlation algorithm more prominent.
Conversely, when one opts for larger block widths, the number of blo from an infrared linear array image of the exact dimensions is reduced.Th time for the column-direction signal increases.Given that the correlation image scene and periodic noise is relatively weak, the influence of the image correlation algorithm is considerably reduced when the block's width is larg In this subsection, we conduct an experiment utilizing 20 infrared l images, each with a resolution of 1024 × 4096.A pure background periodic n identical resolution is acquired through a blackbody.We create a simula image with noise by superimposing this periodic noise image onto the clean scanning image.Using the algorithm described in Section 2.2 and the pure noise image, we identify a signal period covering 847 columns.From these da a single-cycle, one-dimensional noise signal sequence.For the batch of 20 n images, we randomly pick image segments sized 1024 × 847 and segment into blocks.The subsequent denoising process is executed using the approa in Section 3.2.The relationship between the block width of the image corresponding denoising results is illustrated in Figure 9.   Figure 9 displays the average PSNR value across 20 image frames.With an image block width of 1, every row of a 1024 × 847 resolution image is correlated with a singular noise cycle.Even though there are 1024 such image blocks, the scene details notably affect signal correlation, resulting in a post-denoising PSNR average of 30.95 dB.As the image block width expands to four, the number of smaller blocks within the 1024 × 847 image diminishes.However, the scene details' impact on signal correlation decreases, holding the denoised image's PSNR steady at around 32 dB, indicating improved denoising.Given that the infrared images used in this study have 1024 rows, for simplification, further experiments segmented the authentic images into blocks, each sized 256 × 847.

Robustness Experiment
In the context of image acquisition with infrared linear array detectors, the content of the captured scenes significantly influences the denoising effectiveness of the algorithm presented in this paper.For images that primarily feature background elements, the fluctuations in the one-dimensional signal, as elaborated in Section 2.3, are primarily dictated by inherent periodic noise.Thus, the algorithm leverages signal correlation to achieve denoising results.However, when the captured image encompasses intricate background details, the changes in the one-dimensional signal will culminate in both the scene elements and the periodic noise.Given the unpredictable and variable nature of scene content, it can be inferred that such content operates independently from the periodic noise signal.
To systematically replicate the interplay of scene content on the denoising efficiency of the algorithm, this section incorporates a one-dimensional signal characterized by a Gaussian distribution centered at zero and exhibiting diverse variances.This signal is superimposed onto the one-dimensional signal permeated with periodic noise, computed per the methodology above.The underlying logic remains lucid: an augmentation in the variance of the Gaussian distribution signifies increased scene complexity, as illustratively presented in Figure 10.

Robustness Experiment
In the context of image acquisition with infrared linear array detectors, the content of the captured scenes significantly influences the denoising effectiveness of the algorithm presented in this paper.For images that primarily feature background elements, the fluctuations in the one-dimensional signal, as elaborated in Section 2.3, are primarily dictated by inherent periodic noise.Thus, the algorithm leverages signal correlation to achieve denoising results.However, when the captured image encompasses intricate background details, the changes in the one-dimensional signal will culminate in both the scene elements and the periodic noise.Given the unpredictable and variable nature of scene content, it can be inferred that such content operates independently from the periodic noise signal.
To systematically replicate the interplay of scene content on the denoising efficiency of the algorithm, this section incorporates a one-dimensional signal characterized by a Gaussian distribution centered at zero and exhibiting diverse variances.This signal is superimposed onto the one-dimensional signal permeated with periodic noise, computed per the methodology above.The underlying logic remains lucid: an augmentation in the variance of the Gaussian distribution signifies increased scene complexity, as illustratively presented in Figure 10.In Figure 10a, the depicted curve elucidates a one-dimensional signal derived from a chosen image block extracted from an infrared linear array image infused with periodic noise.This periodic noise-pervaded image was synthesized by overlaying a noise-free image with a pure background noise image.
Figure 10b through Figure 10h, respectively, display the outcomes after adding onedimensional Gaussian noise sequences with zero mean and varying variances to the signal in Figure 10a.In this context, σ denotes the standard deviation of the Gaussian noise.
To verify the influence of the one-dimensional Gaussian noise on the denoising efficacy of the presented algorithm, an experimental setup, encapsulated within this section, leveraged 20 frames of immaculate real infrared linear array rotational scanning images, each boasting a resolution of 1024 × 4096.By overlaying these images with the pure background periodic noise sourced from the blackbody, we achieved 20 frames of In Figure 10a, the depicted curve elucidates a one-dimensional signal derived from a chosen image block extracted from an infrared linear array image infused with periodic noise.This periodic noise-pervaded image was synthesized by overlaying a noise-free image with a pure background noise image.
Figure 10b through Figure 10h, respectively, display the outcomes after adding onedimensional Gaussian noise sequences with zero mean and varying variances to the signal in Figure 10a.In this context, σ denotes the standard deviation of the Gaussian noise.To verify the influence of the one-dimensional Gaussian noise on the denoising efficacy of the presented algorithm, an experimental setup, encapsulated within this section, leveraged 20 frames of immaculate real infrared linear array rotational scanning images, each boasting a resolution of 1024 × 4096.By overlaying these images with the pure background periodic noise sourced from the blackbody, we achieved 20 frames of infrared linear array images replete with periodic noise and identical resolution parameters.
As the denoising mechanism of our algorithm took course, diverse variances of Gaussian noise sequences were added to the one-dimensional signal extracted from the periodic noise-laden linear array images.Table 1 lucidly tabulates the aggregate outcomes, underscoring the influence of the one-dimensional Gaussian noise on the algorithm's efficacy.In the results presented in Table 1, it is evident that even with an increase in the standard deviation, σ, of the overlaid one-dimensional Gaussian noise continuously increasing, the PSNR (peak signal-to-noise ratio) of the denoised image remains relatively stable.It underscores the robustness of the proposed algorithm.

Simulation Experiment
To substantiate the denoising efficacy of the presented algorithm, we undertook simulation experiments utilizing noise-free infrared images.We sourced 50 infrared frames from the CAMEL dataset [24] and incorporated one pure background infrared linear array image, exhibiting periodic noise, captured via a blackbody.This noise image, with a resolution of 1024 × 4096, exhibits a noise periodicity spanning roughly 850 columns.Given that the CAMEL dataset [24] contains infrared images of 480 × 640 resolution and acknowledges the additive characteristic of infrared linear array periodic noise, it was imperative to integrate this noise into our selected dataset.The infrared images were resized to 1024 × 1365 from their original 480 × 640 by employing interpolation techniques.The columnar dimension of the noise image was then constricted to produce a resolution of 1024 × 1365 and an inferred noise periodicity of about 284 columns.By computing the mean of the noise image, each pixel value was reduced by this mean and then added to the infrared set, forming a simulated dataset with additive noise.
We then applied two periodic noise removal algorithms, Two-Circle Filter [18] and Local Threshold [19], and three improved infrared image stripe noise denoising methods: the enhanced NL-means algorithm [11], the Wavelet Gradient-Equalization algorithm [16], and the Wavelet TV-L1 algorithm [17].Due to the noise's low-frequency characteristics, the Two-Circle Filter [18] was adjusted to circumvent the image's DC component by configuring the filter radius d to 2.
The enhanced NL-means algorithm [11] operated with neighborhood search blocks sized at 3 × 3, governed by a parameter γ = 500.Both the Wavelet Gradient-Equalization [16] and Wavelet TV-L1 [17] algorithms employed the Haar wavelet basis with a horizontal guiding filter window of 1 × 8 and a vertical one of 10 × 1, and a guided filter parameter ε = 0.16.The Wavelet Gradient-Equalization algorithm [16] set the MHE algorithm parameter n to 2, indicating that it processes two adjacent columns of the pixel.In the proposed method,   As discerned from Figure 11, it is evident that our algorithm provides the best denoising effect.The denoised image is virtually indistinguishable from the noise-free image, devoid of discernible noise stripes.However, denoised images obtained from the Two-Circle Filter [18] and Local Threshold [19] still exhibit periodic noise stripes.The low frequency of infrared linear array periodic noise primarily affects the image's low-frequency components near the DC, unlike other high-frequency periodic noises, which manifest clear peak noise frequencies in their spectra.Solely suppressing noise peak frequencies in the spectrum is not enough for effective denoising.Due to the extended noise period and minimal changes in adjacent column pixels, there are no noticeable stripe noises, making the enhanced NL-means [11], Wavelet Gradient-Equalization [16], and Wavelet TV-L1 [17] less effective.Our proposed strategy integrates a priori knowledge of the intrinsic periodic noise in the background.Through meticulous analysis of the periodic attributes of the noise, our algorithm studies the correlation between noise patterns and image content and then denoises images in the spatial domain.The results are superior.For a quantitative assessment, the pristine noise-free image served as the benchmark against which the PSNR of the denoised image was ascertained, with the detailed outcomes delineated in Table 2, followed by a line chart (Figure 12) to further present the trends.denoising effect.The denoised image is virtually indistinguishable from image, devoid of discernible noise stripes.However, denoised images obta Two-Circle Filter [18] and Local Threshold [19] still exhibit periodic noise st frequency of infrared linear array periodic noise primarily affects the frequency components near the DC, unlike other high-frequency periodic manifest clear peak noise frequencies in their spectra.Solely suppressin frequencies in the spectrum is not enough for effective denoising.Due to noise period and minimal changes in adjacent column pixels, there are no no noises, making the enhanced NL-means [11], Wavelet Gradient-Equaliza Wavelet TV-L1 [17] less effective.Our proposed strategy integrates a priori the intrinsic periodic noise in the background.Through meticulous an periodic attributes of the noise, our algorithm studies the correlation b patterns and image content and then denoises images in the spatial domai are superior.For a quantitative assessment, the pristine noise-free image benchmark against which the PSNR of the denoised image was ascertai detailed outcomes delineated in  Table 2 evaluates the average PSNR across the 50 denoised infrared method outperforms other algorithms.The PSNR of images denoised us Circle algorithm [18] was the lowest.It can be attributed to the proximity of l periodic noise peaks in the frequency domain to the DC component.
Despite calibrating the filter radius r to a value of 2, the intrinsic ope Table 2 evaluates the average PSNR across the 50 denoised infrared frames.Our method outperforms other algorithms.The PSNR of images denoised using the Two-Circle algorithm [18] was the lowest.It can be attributed to the proximity of low-frequency periodic noise peaks in the frequency domain to the DC component.
Despite calibrating the filter radius r to a value of 2, the intrinsic operations of the denoising algorithm might have inadvertently perturbed the DC component within the frequency domain, consequently altering the global mean value of the image frame.Such a phenomenon correlates with the visually diminished luminosity observed in images processed using the Two-Circle algorithm, as depicted in Figure 11.
Table 3 delineates the various denoising algorithms' SSIM (structural similarity) indexes.It is evident that our proposed method consistently achieves superior results in preserving structural integrity.While the Wavelet TV-L1 approach and other traditional algorithms deliver commendable performance, the SSIM index for our proposed method is consistently closer to unity, underscoring its ability to retain the inherent structure of the image post-denoising.This superior structural preservation, as highlighted in our algorithm, complements the PSNR results presented earlier and further consolidates the efficacy of our method, as vividly depicted in the comparison of the datasets.

Real Data Experiment
To validate the effectiveness of the proposed algorithm in the context of authentic high-resolution infrared linear array images, this section employs two real infrared linear array scanned images with dimensions of 1024 × 4096 containing periodic noise for the denoising experiment, along with one pure background noise image obtained during the equipment's non-uniformity correction process.For the Wavelet Gradient-Equalization algorithm [16] and the Wavelet TV-L1 algorithm [17], the vertical guiding filter window dimensions were modified to 100 × 1, optimizing their capacity to extract stripe noise in high-resolution imagery.Other algorithm parameters remain consistent with those in Section 3.4.The denoising results are graphically represented in Figure 13.
In Figure 13, the proposed algorithm effectively removes periodic noise while preserving the details of the image, especially the intricate features within the scene.On the other hand, the Two-Circle Filter [18] and Local Threshold [19] algorithms significantly blur the detailed information about the clouds in the sky.Such distortion can be traced back to the confluence of low-frequency elements, such as the clouds mentioned above, and the frequency domain characteristics of the periodic noise.This overlap poses a distinct challenge to denoising endeavors, mainly when conducted directly within the frequency domain, further compounded by increasing image resolution.
Furthermore, in the real infrared linear array-scanned images, the noise period is approximately 850 columns, much higher than that in the simulated images.Consequently, the change in pixel values between adjacent image columns is even more gradual.Such slow variation exacerbates the inadequacy of the improved NL-means algorithm [11], Wavelet Gradient-Equalization algorithm [16], and Wavelet TV-L1 algorithm [17] when confronted with this specific type of noise.Their denoising efficacy under these circumstances is less than ideal.In Figure 13, the proposed algorithm effectively removes periodic noise while preserving the details of the image, especially the intricate features within the scene.On the other hand, the Two-Circle Filter [18] and Local Threshold [19] algorithms significantly blur the detailed information about the clouds in the sky.Such distortion can be traced back to the confluence of low-frequency elements, such as the clouds mentioned above, and the frequency domain characteristics of the periodic noise.This overlap poses a distinct challenge to denoising endeavors, mainly when conducted directly within the frequency domain, further compounded by increasing image resolution.
Furthermore, in the real infrared linear array-scanned images, the noise period is approximately 850 columns, much higher than that in the simulated images.Consequently, the change in pixel values between adjacent image columns is even more gradual.Such slow variation exacerbates the inadequacy of the improved NL-means algorithm [11], Wavelet Gradient-Equalization algorithm [16], and Wavelet TV-L1

Real-Time Performance of the Algorithm
Given the fast imaging capabilities and high resolution provided by infrared linear array detectors, it is imperative for associated algorithms to maintain optimal real-time performance.In this section, we assess the denoising time for each algorithm, utilizing authentic infrared linear array images (resolution 1024 × 4096) afflicted with periodic noise.The specifics of these algorithms remain consistent with Section 3.5.The computational environment comprises an Intel i7-9750H processor with 16 GB of memory, operating on the Matlab2020b platform.
Table 4 delineates the average runtime for each algorithm to denoise an individual image frame across ten trials.Notably, the enhanced NL-means algorithm [11] demands computational time for proximity search around every pixel, thereby prolonging its duration.The Wavelet Gradient-Equalization [16] and Wavelet TV-L1 [17] algorithms necessitate MHE and L1 regularization processes post-wavelet decomposition on sub-images (resolution 512 × 2048).This requirement incrementally extends their processing durations, making them less compatible with the real-time demands of infrared linear array detectors.Contrastingly, our proposed algorithm capitalizes on the unique periodic noise variation in the columnar direction.By leveraging periodic signal correlation, it facilitates direct spatial domain image denoising, ensuring minimal time consumption.As the image resolution escalates, there is an evident proportional increase in the runtime across all algorithms.This increment is particularly pronounced in methods like the improved NL-means, which exhibit a substantial time surge even with modest increments in resolution.However, a compelling observation is the minimal runtime escalation of our proposed FACD method, which manifests only slight increments even when the image resolution multiplies.This indicates the scalability and robustness of our approach, making it particularly adept for applications that might necessitate varied image resolutions.Furthermore, when juxtaposed with other algorithms, the FACD's time increment remains negligible, reinforcing its prowess in ensuring consistent and rapid denoising irrespective of resolution expansions.
Our algorithm emerges as a potent tool within infrared linear array imaging systems.According to the characteristics of this system, it first needs to close the baffle to perform the non-uniformity correction every time it powers up.Collect a frame of the image as the pure background image at this moment.Subsequent analysis yields noise periodicity and amplitude.The baffle is then reopened, and leveraging our algorithm's methodology, the noise's correlation with authentic images is ascertained, facilitating spatial domain denoising.After completing the above steps, the image denoising process is simply arithmetic operations on every pixel after receiving a column image.This streamlined approach ensures brevity in execution, addressing long-standing challenges posed by periodic noise due to electromechanical disturbances within infrared linear array imaging systems.

Discussion
Within the confines of the experiments outlined in the preceding section, our proposed algorithm exhibited marked proficiency in attenuating low-frequency noise, while simultaneously upholding commendable real-time operational efficiency.When one broaches the subject of high-resolution, real-time image processing systems, it becomes pivotal to note that the periodicity of columns in infrared linear array detectors might shrink to approximately 30 microseconds.Such systems predominantly find utility in the domains of target extraction and recognition.Given the prodigious generation rate, surpassing 30,000 image columns every second, resorting to a full-frame buffering approach for image elaboration would entail a substantial memory overhead, jeopardizing the real-time assurance.Conventionally, systems of this ilk adopt a staggered processing paradigm, buffering at intervals of 128 or 256 columns.In other words, with every caching instance, the system extracts targets within the image to achieve real-time target output.Thus, the advantages of column-by-column noise reduction are evident.Our algorithm can be deployed during the embedded image preprocessing phase.After performing rapid noise reduction on each column, the image is passed to the next stage.The resultant computational footprint is minimal, both in terms of resource consumption and parameter storage overhead.Such efficiency earmarks ample computational latitude for downstream algorithms, rendering our method an ideal fit for embedded systems that are relatively resource-constrained.
However, our algorithm requires a preliminary step-namely, collecting pure noise images against a blackbody or a baffle.It may not resonate well with imaging systems without a baffle or a blackbody.Although one could obtain relatively clean background images by tilting the detector toward the sky, this method has limitations (for instance, when the working environment is indoors or when the sky has complex cloud layers).Future endeavors could contemplate refining this preliminary step, striving for simultaneous acquisition of periodic noise frequencies while not compromising the real-time operation.While there is potential to dovetail neural network frameworks into our model, a salient point of contention remains.As underscored earlier, extant deep convolutional neural networks (CNNs) tailored for noisy image mitigation predominantly fixate on entire frame denoising.Their deployment complexities, especially when juxtaposed with real-time processing demands, stand as impediments.Additionally, the conspicuous absence of high-resolution training datasets further hinders the optimization and effectiveness of such neural networks.

Conclusions
This research offers a novel solution to the prevalent challenge of low-frequency periodic signal noise induced by electromechanical interferences in infrared linear array detectors.We introduced an efficient denoising algorithm grounded in signal correlation principles, focusing on extracting and utilizing noise signals from pure backgrounds to achieve comprehensive denoising.The robustness and efficiency of our approach were validated through extensive simulations and real imagery tests, demonstrating notable enhancements in PSNR and superior denoising capabilities compared to existing methods.Notably, our algorithm maintained high performance even under substantial perturbations and noise variations, evidencing its adaptability to different noise conditions.Furthermore, the algorithm exhibited remarkable real-time efficiency in processing high-resolution images, showcasing its practical applicability and readiness for deployment in the embedded image preprocessing phase, ensuring minimal computational footprint and resource consumption.However, it is crucial to note that our algorithm necessitates a preliminary step of collecting pure noise images, which may pose limitations in certain working environments or system configurations.Future endeavors could focus on refining this preliminary step and exploring the integration of neural network frameworks, addressing the deployment complexities associated with real-time processing and the lack of high-resolution training datasets.

Figure 1 .
Figure 1.(a) Traditional periodic noise.(b) The spectrum of the traditional periodic noise.

Figure 1 .
Figure 1.(a) Traditional periodic noise.(b) The spectrum of the traditional periodic noise.

Figure 1 .
Figure 1.(a) Traditional periodic noise.(b) The spectrum of the traditional periodic noise.

Figure 2 .
Figure 2. (a) Noise image obtained through blackbody.(b) Two-dimensional Fourier frequency domain image.Upon juxtaposing the frequency domain images of periodic noise depicted in Figures1 and 2, one discerns a clear contrast.Traditional periodic noise, when viewed in the frequency domain, is characterized by multiple bright regions distanced from the DC component.It indicates that such noise primarily comprises high-intensity, high-frequency components, significantly affecting

Figure 3 .
Figure 3. Analysis of infrared linear array scanning images: (a) Statistics of pixel mean values for each column of images.(b) Statistics of pixel mean values for each row of images.

Figure 3 .
Figure 3. Analysis of infrared linear array scanning images: (a) Statistics of pixel mean values for each column of images.(b) Statistics of pixel mean values for each row of images.

Figure 4 .
Figure 4.The same scene contains noise of different phases.(a) The image with noise in o (b) The image with the same noise but different phase

Figure 4 .
Figure 4.The same scene contains noise of different phases.(a) The image with noise in one phase (b) The image with the same noise but different phase.
) = 1  ()( −  − min(( − ), 0) × )  () = 1  ()( −  − (( − ), 0) × ) Let () and () be one-dimensional sequence signals with ∈ [1 signifies the length of the sequence and  (•) is the signum function.D uniformity correction stage of the detector, a periodic noise image is sec pure background, as exemplified in Figure 2a.The average value for each image is determined as outlined in Equation (13).() = 1   (, ) In this context, () refers to the one-dimensional sequence depicted which is recognized as the multi-cycle noise present in the image's column The term  (, ) signifies a pixel within the noise image when set a background.Here, with ∈ [1, ] ,  ∈ [1, ] , where  represents the to columns in the image and R indicates its rows.From the sequence () contiguous sequence with length paralleling the noise cycle  , termed ( within the interval  ∈ 1,  .

Figure 5 .
Figure 5.The process of image block segmentation.

Figure 5 .
Figure 5.The process of image block segmentation.

Figure 6 .
Figure 6.Flowchart of the proposed denoising algorithm.The flowchart diagram illustrates that the first step involves employing Fourier Transform to calculate the periodicity of the absorbed signals.The upper section of the diagram shows the extraction of a single-period signal from images containing objects, represented as a one-dimensional curve.Meanwhile, the lower section represents the extraction of single-period signals with different phases using a black body.Subsequently, noise reduction is achieved in the spatial domain through the utilization of correlation.

Figure 6 .
Figure 6.Flowchart of the proposed denoising algorithm.The flowchart diagram illustrates that the first step involves employing Fourier Transform to calculate the periodicity of the absorbed signals.The upper section of the diagram shows the extraction of a single-period signal from images containing objects, represented as a one-dimensional curve.Meanwhile, the lower section represents the extraction of single-period signals with different phases using a black body.Subsequently, noise reduction is achieved in the spatial domain through the utilization of correlation.

Figure 7 .
Figure 7. Histograms for noise periods against cycle counts captured in 0 (a), 6 (b), 12 (c), and 18 (d) hour(s) after the initial power-on of an infrared linear array detector.In the histogram, the yellow bars represent the counts of signals at different periods.

Figure 8 .
Figure 8. Noise amplitude fluctuation.It shows the amplitude fluctuations of the signal typically range between approximately −30 and 30.
of colum ns for sig nal p eriod

Figure 7 .
Figure 7. Histograms for noise periods against cycle counts captured in 0 (a), 6 (b), 12 (c), and 18 (d) hour(s) after the initial power-on of an infrared linear array detector.In the histogram, the yellow bars represent the counts of signals at different periods.

Figure 7 .
Figure 7. Histograms for noise periods against cycle counts captured in 0 (a), 6 (b), 1 hour(s) after the initial power-on of an infrared linear array detector.In the histog bars represent the counts of signals at different periods.

Figure 8
presents this a the x-axis demarcates experimental batches of 22 periodic signals and the y amplitude variations in pixel values.The noise amplitude across all periods remarkable stability, fluctuating minimally between −30 and 30 pixels.

Figure 8 .
Figure 8. Noise amplitude fluctuation.It shows the amplitude fluctuations of the si range between approximately −30 and 30.

Figure 8 .
Figure 8. Noise amplitude fluctuation.It shows the amplitude fluctuations of the signal typically range between approximately −30 and 30.

Figure 9 .
Figure 9.The influence of image block width on algorithm denoising.

Figure 9
Figure 9 displays the average PSNR value across 20 image frames.W block width of 1, every row of a 1024 × 847 resolution image is correlated wi noise cycle.Even though there are 1024 such image blocks, the scene details n signal correlation, resulting in a post-denoising PSNR average of 30.95 dB.A block width expands to four, the number of smaller blocks within the 1024 diminishes.However, the scene details' impact on signal correlation decrea the denoised image's PSNR steady at around 32 dB, indicating improved den

Figure 9 .
Figure 9.The influence of image block width on algorithm denoising.

Figure 10 .
Figure 10.Simulating image background complexity using Gaussian noise with zero mean and diverse variances.With an increase in the variance of added noise, there is a noticeable amplification in the signal curve's fluctuations within the images.This trend leads to a progressive loss of clarity in the distinctive features of the original signal, resulting in an overall blurring effect.

Figure 10 .
Figure 10.Simulating image background complexity using Gaussian noise with zero mean and diverse variances.With an increase in the variance of added noise, there is a noticeable amplification in the signal curve's fluctuations within the images.This trend leads to a progressive loss of clarity in the distinctive features of the original signal, resulting in an overall blurring effect.
were selected in Step 2. Five images were handpicked from the dataset to qualitatively analyze the denoising effects of different algorithms, as shown in Figure 11.Sensors 2023, 23, 8716 14 of 21 handpicked from the dataset to qualitatively analyze the denoising effects of different algorithms, as shown in Figure 11.

Figure 13 .
Figure13.Denoising effect comparison on real infrared images.Both "Denoised by Two-Circle" and "Denoised by Local Threshold" methods, when applied in the frequency domain, led to a significant deterioration of the low-frequency information in the images.Conversely, "Denoised by Improved NL-means", "Denoised by Wavelet and Gradient Equalization", and "Denoised by Wavelet and TV-L1" methods exhibited limitations in effectively eliminating this low-frequency periodic noise.In contrast, our proposed algorithm demonstrated a notable capability to successfully remove this type of noise, providing a more effective denoising solution.

Figure 13 .
Figure13.Denoising effect comparison on real infrared images.Both "Denoised by Two-Circle" and "Denoised by Local Threshold" methods, when applied in the frequency domain, led to a significant deterioration of the low-frequency information in the images.Conversely, "Denoised by Improved NL-means", "Denoised by Wavelet and Gradient Equalization", and "Denoised by Wavelet and TV-L1" methods exhibited limitations in effectively eliminating this low-frequency periodic noise.In contrast, our proposed algorithm demonstrated a notable capability to successfully remove this type of noise, providing a more effective denoising solution.
It indicates that such noise primarily comprises high-intensity, highfrequency components, significantly affecting image quality.Conversely, the infrared linear array frequency domain image scarcely showcases distinct luminous areas, revealing only sporadic bright spots proximate to the DC component.It underscores the predominance of noise characterized by reduced intensity and frequency.Employing a standard Gaussian notch filter to address this periodic noise can unintentionally weaken the image's inherent low-frequency elements, negatively impacting the denoising outcome.

Table 1 .
Effect of Gaussian noise on algorithmic denoising.

Table 2 .
Denoising results of different algorithms.

Table 2
, followed by a line chart (Figure present the trends.

Table 2 .
Denoising results of different algorithms.Line chart for Table2: denoising results of different algorithms.

Table 3 .
SSIM index of different algorithms.

Table 4 .
Average run time of different algorithms.