Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Turbulence-resistant free space optical communication via chaotic block-matching and 3D filtering

Open Access Open Access

Abstract

In this paper, we propose a chaotic block-matching and three-dimensional (C-BM3D) filtering algorithm to remove the noise and enhance the security in the turbulent channel of free space optical (FSO) communication. We experimentally demonstrate the performance of C-BM3D by comparing it with chaotic non-local means filtering (C-NLM), chaotic Gaussian filtering and chaotic Median filtering based on Log-normal and Gamma-Gamma turbulence models. The results show that the peak signal-to-noise ratios (PSNRs) of C-BM3D in the weak turbulence under Log-normal and Gamma-Gamma models are up to 96.2956 and 93.2853, respectively. The C-BM3D also achieves superior image similarity in Log-normal turbulent channel, with its structural similarity index measures (SSIMs) nearly equal to 1. Additionally, the signal-to-noise ratio (SNR) of C-BM3D ranks the highest, and its bit error rate (BER) improves by at least 15 dB compared to that of the other three algorithms. The experimental results indicate that the C-BM3D can be a good candidate for the next generation of FSO communication in security and turbulence resistance.

© 2024 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Free space optical (FSO) communication has been widely applied to low earth orbit (LEO) satellite, underwater and unmanned aerial vehicle (UAV) communications because of its license-free, easy deployment, high capacity, etc. For example, C. Li et al. studied a wavelength-division-multiplexing (WDM) four-level pulse amplitude modulation (PAM4) FSO-underwater wireless optical communication (UWOC) integrated system based on 100 Gb/s channel capacity [1]. Since FSO systems do not require license fees, FSO has been viewed as an alternative technology to conventional radiofrequency (RF) systems [24] in geostationary satellites [57]. To further improve the FSO communication performance, relay-assisted FSO system has been proposed in [8], which outperforms the pure FSO system. FSO links [9] achieve robust and low-latency communications between autonomous underwater vehicle (AUV) and UAV. Authors in [10] also combined UAV with dual-hop RF/FSO relay system to achieve extensive and reliable communications in space-air-ground integrated networks (SAGIN).

However, turbulence is one of the main factors that adversely influences the performance of FSO communication. If a laser transmits through turbulence, intensity fluctuation and optical power attenuation exist, which affects communication performance in turn [11]. There have been many studies about transmitting general information in FSO channels [1215], which achieved great progress in resisting turbulence. However, there is little research on FSO image transmission. The previous experiments of traditional Gaussian filtering [16] and Median filtering [17] verified that although they have good denoising effects, they are at the expense of communication quality in FSO transmission. The newly proposed non-local means (NLM) filtering [18] can make full use of redundant information and maintain the details of the image, but the noise of similar blocks increases when searching for weights. Based on the above reduction methods, block-matching and three-dimensional (BM3D) filtering is proposed [19] to combine the spatial and frequency domain algorithms, which achieves better visual effects and higher similarity. The latest studies [2023] show that, compared with traditional filtering algorithms, two-step iterative shrink/threshold segmentation and enhanced Lagrange algorithm, the introduction of the BM3D algorithm not only improves the image reconstruction quality but also enhances the anti-noise capability of the technique.

Due to the exposure of optical signals in FSO communication, the transmitted signals are prone to be eavesdropped and accessed by illegal attackers [24]. The chaotic technique has been widely used in image encryption [25] and security enhancement [26,27] of fiber optic communication thanks to its high sensitivity to initial values, strong unpredictability, etc. For instance, WDM signals can achieve secure transmission on 50 km standard single-mode optical fiber using the private chaotic phase in disarray [28]. Chaotic compressive sensing (CS) successfully enhances the orthogonal frequency division multiplexing-passive optical network (OFDM-PON) transmission security [29]. Chaotic three-dimensional (3D) constellation scrambling could effectively enhance the security and transmission performance of coherent optical orthogonal frequency-division multiplexing (CO-OFDM) systems [30].

In this paper, we combine logistic chaotic mapping with the BM3D filtering algorithm to resist the turbulence and enhance the security of FSO communication. Specifically, the image is encrypted at the transmitter by digital chaos and decrypted at the receiver. At the receiver, the image with turbulence will be filtered by the C-BM3D and other typical filtering algorithms. To clarify the superiority of C-BM3D, we will compare the bit error rate (BER), structural similarity index measure (SSIM) and peak signal-to-noise ratio (PSNR) performance of different filtering algorithms in different turbulent conditions.

2. Principle

2.1 Atmospheric turbulence in FSO channels

Among common statistical models for simulating turbulent conditions, the Log-normal and Gamma-Gamma models are widely used to study turbulence-resistant techniques [31]. The probability distribution function (PDF) of Log-normal is expressed as [32]:

$$p(I) = \frac{1}{{\sqrt {2\pi \sigma _l^2} }}\frac{1}{I}\textrm{exp} ( - \frac{{{{(In(\frac{1}{{{I_0}}} + \frac{{\sigma _l^2}}{2}))}^2}}}{{2\sigma _l^2}}),$$
where I is the irradiance intensity in turbulent medium, I0 is the irradiance intensity in free space, σ2 l is log irradiance variance that represents the strength of atmospheric turbulence. The PDF of Gamma-Gamma is expressed as:
$${f_{{\gamma _b}}}({\gamma _b}) = \frac{{\alpha {\beta ^{(\alpha + \beta )/2}}{{\bar{\gamma }}_b}^{ - (\alpha + \beta )/4}}}{{\Gamma (\alpha )\Gamma (\beta )}}{\gamma _b}^{ - (\alpha + \beta )/4 - 1} \times {K_{\alpha - \beta }}[2\sqrt {\alpha \beta \sqrt {\frac{{{\gamma _b}}}{{{{\bar{\gamma }}_b}}}} } ],$$
where Γ(.) represents the Gamma function, and Kn (.) is the modified Bessel function of the second kind of order n. The ${\bar{\gamma }_b}$ is average signal-to-noise ratio (SNR) per bit, given by ${\bar{\gamma }_b} = E[{A^2}]{E_b}/{N_0}$. The α and β are given by (3) and (4), which are shaping parameters that are related to the effective number of large-scale and small-scale eddies in the scattering process, respectively.
$$\alpha = {[\textrm{exp} (\frac{{0.49\sigma _R^2}}{{{{(1 + 1.11\sigma _R^{12/5})}^{5/6}}}}) - 1]^{ - 1}},$$
$$\beta = {[\textrm{exp} (\frac{{0.51\sigma _R^2}}{{{{(1 + 0.69\sigma _R^{12/5})}^{5/6}}}}) - 1]^{ - 1}},$$
where σR is the Rytov variance, namely the intensity of turbulence.

2.2 Logistic chaotic mapping

Logistic mapping is a typical nonlinear iterative equation, as shown in the equation:

$${x_{k + 1}} = \mu \times {x_k}(1 - {x_k}),k = 0,1, \cdots \cdots ,n,$$
where μ is the control parameters, μ∈ (0, 3.5699456], xk is the chaotic sequence, x0∈ (0, 1), and k is the iteration time step.

One main characteristic of chaotic systems is the high initial sensitivity. When μ = 0.190903252535798, x1 = 0.699076722656686, x2 = 0.699076722656687, the initial sensitivity of the logistic chaotic mapping is shown in Fig. 1. After certain iterations, the mixed chaotic sequence splits into two sequences with different translations. As the number of iterations increases, the values of the two sequences show bigger difference. It shows that the chaotic system still maintains the initial value sensitivity of logistic chaotic mapping.

 figure: Fig. 1.

Fig. 1. Initial value sensitivity.

Download Full Size | PDF

2.3 BM3D algorithm

BM3D algorithm [33] can find similar blocks to reference blocks of two-dimensional (2D) images. It stacks 2D-image similar blocks into 3D groups, which can be performed by collaborative filtering. We estimate the recovered image by aggregating the processed results and original image blocks. The overall process of the algorithm is shown in Fig. 2, which can be divided into two steps: Hard-thresholding filtering and Wiener filtering.

 figure: Fig. 2.

Fig. 2. Execution process of BM3D and the reference block is marked with “R”.

Download Full Size | PDF

We use Zx to represent a square block of N1×N1 on the noisy image z with the top-left coordinate x. xX, X is the image domain. Zs is a 3D block, and the top-left coordinates of all 2D blocks are defined by the set S. To represent the reference block in each 3D block, we add the R (reference block) identifier under the corresponding coordinate, namely xR. Similarly, we use Yx to represent a block on a noiseless image y. Because Step1 and Step2 are the same structure, the meanings of some parameters are also the same. In order to distinguish the parameters in the two steps, we add the “ht” and “wie” marks, respectively. For the basic estimation image of Step1, we use ${\hat{y}^{basic}}$ to represent it, and the final estimation image obtained by Step2 is ${\hat{y}^{final}}$.

  • 1. Basic estimates: Hard-thresholding filtering

    (1) Estimates by block: For each block of the noisy image:

    a) Grouping. After finding similar blocks to the reference block, they are altogether stacked into a 3D array. When the noise variance is large or the block is relatively small, it is not accurate to find similar blocks based on the noisy image z. Therefore, the two blocks on the noisy image z are separated by 2D normalized transformation. Then, the coefficients whose amplitude is less than a certain threshold are set to zero, so that the matching error of the two blocks can be expressed as the mean square error of these coefficients.

    $$d({Z_{{x_R}}},{Z_x}) = \frac{{||\gamma ^{\prime}(\tau _{2D}^{ht}({Z_{{x_R}}})) - \gamma ^{\prime}(\tau _{2D}^{ht}({Z_x}))||_2^2}}{{{{(N_1^{ht})}^2}}},$$
    where $\gamma ^{\prime}$ is a hard threshold operator with threshold ${\lambda _{2D}}\sigma $, ||.||22 is the 2-norm, $N_1^{ht}$ is the block size used in basic estimates, $\tau _{2D}^{ht}$ is the 2D normalized transformation.

    After calculating the matching error between the current reference block and all other blocks according to Eq. (6), we only keep those blocks whose error is less than a certain threshold and get the corresponding coordinate set:

    $$S_{{x_R}}^{ht} = \{ x \in X:d({Z_{{x_R}}},{Z_x}) \le \tau _{match}^{ht}\} ,$$
    where $\tau _{match}^{ht}$ is the maximum error threshold for determining the similarity of two blocks, which is generally set according to experience. The current reference block will be judged to be a similar block because its matching error is 0, so there is at least one element in the coordinate set. By stacking all the similar blocks together, we get a 3D array ${Z_{S_{xR}^{ht}}}$ of shape $N_1^{ht} \times N_1^{ht} \times |S_{{x_R}}^{ht}|$, where $|S_{{x_R}}^{ht}|$ represents the number of elements in the collection.

    b) Collaborative hard thresholding. After obtaining the 3D array corresponding to the reference block, we can carry out 3D collaborative transformation and filtering, which can be formally expressed as (8). Firstly, the 3D normalized transformation is applied to the above groups. After obtaining the non-zero component from the weight coefficient by hard threshold, the inverse 3D normalized transformation is performed to get the estimation of all grouped blocks. Finally, they are returned to their original positions.

    $$\hat{Y}_{S_{{x_R}}^{ht}}^{ht} = \tau {_{3D}^{ht}{^{ - 1}}}(\gamma (\tau _{3D}^{ht}({Z_{S_{{x_R}}^{ht}}}))),$$
    where $\gamma$ is a hard threshold operator with threshold ${\lambda _{3D}}\sigma $, ${Z_{S_{{x_R}}^{ht}}}$ is the formed 3D groups, $\tau _{3D}^{ht}$ is the 3D normalized transformation and $\tau {_{3D}^{ht}{^{ - 1}}}$ is the inverse 3D normalized transformation.

    (2) Aggregation: The superposition weights are calculated by the number of non-zero components. Finally, the estimated image is obtained by dividing the superimposed blocks by the weights of each block.

    $${\hat{y}^{basic}}(x) = \frac{{\sum\limits_{{x_R} \in X} {\sum\limits_{{x_m} \in S_{{x_R}}^{ht}} {\omega _{{x_R}}^{ht}\hat{Y}_{{x_m}}^{ht,{x_R}}(x)} } }}{{\sum\limits_{{x_R} \in X} {\sum\limits_{{x_m} \in S_{{x_R}}^{ht}} {\omega _{{x_R}}^{ht}{\chi _{{x_m}}}(x)} } }},\forall x \in X.$$

    We assume that each block has been filled to the same size as the original image by zeros and ${\chi _{{x_m}}}:X \to \{ 0,1\}$ is used to determine whether a pixel x is on block xm.

  • 2. Final estimates: Wiener Filtering

    (1) Estimates by block: For each block of the obtained image from basic estimates:

    a) Grouping. Within the basic estimates range, block matching (BM) is used to find the location of all blocks that are similar to the reference block, namely R. In this section, we find similar blocks by (10) based on the basic estimation image. Because the current noise is relatively small, it is not necessary to do transformation and hard-thresholding like (6).

    $$S_{{x_R}}^{wie} = \left\{ {x \in X:\frac{{||\hat{Y}_{{x_R}}^{basic} - \hat{Y}_x^{basic}||_2^2}}{{{{(N_1^{wie})}^2}}} < \tau_{match}^{wie}} \right\},$$
    where $N_1^{wie}$ is the block size used in Wiener filtering, $\tau _{match}^{wie}$ is the maximum d-distance for which two blocks are considered similar, $\hat{Y}_{{x_R}}^{basic}$ and $\hat{Y}_x^{basic}$ are the basic estimate blocks located at xR and x, respectively. Now there are two groups formed, one from the noisy image and another one from the basic estimates: ${\hat{Y}^{wie}}_{S_{xR}^{wie}}$ and ${Z_{S_{xR}^{wie}}}$. ${\hat{Y}^{wie}}_{S_{xR}^{wie}}$ is obtained by stacking together the basic estimate blocks $\hat{Y}_{x \in S_{{x_R}}^{wie}}^{basic}$, and ${Z_{S_{xR}^{wie}}}$ is obtained by stacking together the noisy blocks ${Z_{x \in S_{{x_R}}^{wie}}}$.

    b) Collaborative Wiener filtering. The 3D normalized transformation is applied to the above two groups. The Wiener filtering is performed on the noisy one. Next, the inverse 3D normalized transformation is applied to the filtered Wiener shrinkage coefficients to estimate all blocks. Finally, these block estimates are returned to their original positions. The estimates are shown as:

    $${\hat{Y}^{wie}}_{S_{{x_R}}^{wie}} = \tau _{3D}^{wi{e^{ - 1}}}({W_{S_{{x_R}}^{wie}}}\tau _{3D}^{wie}({Z_{S_{{x_R}}^{wie}}})),$$
    where ${W_{S_{{x_R}}^{wie}}}$ is the Wiener shrinkage coefficient, ${Z_{S_{{x_R}}^{wie}}}$ is the noisy blocks in final estimates.

    (2) Aggregation: The final estimate ${\hat{y}^{final}}(x)$ of the real image is calculated by aggregating all the obtained local estimates using the weighted average.

    $${\hat{y}^{final}}(x) = \frac{{\sum\limits_{{x_R} \in X} {\sum\limits_{{x_m} \in S_{{x_R}}^{wie}} {\omega _{{x_R}}^{wie}\hat{Y}_{{x_m}}^{wie,{x_R}}(x)} } }}{{\sum\limits_{{x_R} \in X} {\sum\limits_{{x_m} \in S_{{x_R}}^{wie}} {\omega _{{x_R}}^{wie}{\chi _{{x_m}}}(x)} } }},\forall x \in X.$$

3. Simulation process and analysis

In this experiment, we aim to simulate and verify the denoising effect and communication quality of the C-BM3D algorithm by comparing it with the C-NLM algorithm, chaotic Gaussian filtering and chaotic Median filtering under Log-normal and Gamma-Gamma turbulence models. Specifically, the matrix size of the experimental image is 256 × 256. Firstly, the image is encrypted by logistic chaotic mapping. The logistic map is used to generate an array of sequence of random numbers. After obtaining random sequence, the image is cluttered and diffused at the bit plane. Finally, the bitstream of 1*65536 is got by XOR operation on the random image and the random sequence by bit, which is transmitted through the FSO channel. The decryption process is the reverse of encryption. At the receiver, the image recovered from the received bitstream will be filtered by C-BM3D, C-NLM, chaotic Gaussian filtering and chaotic Median filtering. The experimental FSO communication system is shown in Fig. 3, and the simulation parameters are listed in Table 1.

 figure: Fig. 3.

Fig. 3. The experimental FSO communication system.

Download Full Size | PDF

Tables Icon

Table 1. The simulation parameters

Additionally, we calculate the SNR and BER to evaluate the communication quality after the four algorithms. And we calculate the PSNR and SSIM of the recovered image and original image to judge the image quality. The formulas are as:

$$SNR = 10\ast {\log _{10}}(Ps/Pn)$$
where Ps represents the average power of the signal and Pn represents the average power of the noise.
$${P_b} = {I_e}/I,$$
where Ie represents the number of bits accepted in error, and I represents the total number of transmitted bits.
$$PSNR(f,g) = 10\lg \frac{{{{({2^n} - 1)}^2}}}{{\sum\limits_i {\sum\limits_j {\frac{{{{(f(i,j) - g(i,j))}^2}}}{{MN}}} } }},$$
where f(i, j) is the original image, g(i, j) is the restored image. i and j are the image pixels. M × N is the image size.
$$SSIM(f,g) = \frac{{(2{\mu _f}{\mu _g} + {c_1})(2{\sigma _{fg}} + {c_2})}}{{(\mu _f^2 + \mu _g^2 + {c_1})(\sigma _f^2 + \sigma _g^2 + {c_2})}},$$
where μf is the mean value of f and μg is the mean value of g. σ2f is the variance of f and σ2 g is the variance of g. σfg is the covariance of f and g. c1 = (k1 L)2 and c1 = (k1 L)2 are two constants to avoid division by zero. The closer the SSIM value is to 1, the higher similarity between the original image and the recovered image.

3.1 Image quality under turbulence of different intensity

In this section, we aim to compare the PSNR and SSIM of filtering algorithms to estimate the image quality under Log-normal and Gamma-Gamma turbulence models by combining Opti-System with MATLAB. The experiment will be carried out in three different intensities of turbulence in each model: weak, moderate and strong. The values of attenuation (a) are set at 5 dB, 12 dB, 15 dB, respectively, and the values of index refraction structure (c2n) are set at 8 × 10−15m−2/3, 3 × 10−14m−2/3, 1 × 10−13m−2/3, respectively. The results of Log-normal and Gamma-Gamma turbulence models are illustrated in Table 2 and Table 3, respectively. From the data in the two tables below, the PSNRs of C-BM3D in weak turbulence under Log-normal and Gamma-Gamma models are up to 96.2956 and 93.2853, respectively, which are much higher than that of the other three algorithms. Furthermore, C-BM3D and chaotic Gaussian filtering achieve superior image similarity in the Log-normal turbulent channel, with their SSIMs nearly equal to 1. However, it is also apparent that the PSNRs of C-BM3D in strong turbulence under the Gamma-Gamma model present dramatic decrease, even as low as 26.0561. Under strong turbulence, the noise is so dense that the similar blocks may contain some noise. The denoising effect decreases with the estimation error increases. For chaotic Median filtering and C-NLM, their values do not experience obvious change, which have been kept at lower numbers than C-BM3D and chaotic Gaussian filtering. Overall, these results indicate that the C-BM3D has good image denoising and reconstruction performance, with its higher PSNRs and SSIMs in Log-normal and Gamma-Gamma turbulent channels.

Tables Icon

Table 2. Comparison in Log-normal Turbulent Channel

Tables Icon

Table 3. Comparison in Gamma-Gamma Turbulent Channel

3.2 Communication quality under turbulence of different intensity

In this section, we demonstrate the SNR distribution of filtering algorithms to estimate their communication quality under the Gamma-Gamma turbulence model (Since C-BM3D has already shown superior PSNR performance in the Log-normal turbulent channel from the former section). The experiment will be carried out in four different intensities of turbulence, as shown in Fig. 4. The values of a are set at 5 dB, 10 dB, 12 dB and 15 dB, respectively, and the values of c2n are set at 8 × 10−15m−2/3, 8 × 10−14m−2/3, 8 × 10−13m−2/3 and 8 × 10−12 m−2/3, respectively. It is obvious from Fig. 4 that the SNR of C-BM3D has ranked the highest, even peaking at over 100, while the SNRs of other algorithms have kept at less than 40. Thus, C-BM3D has greater advantages in improving communication quality in FSO communication. However, it can also be seen that C-BM3D shows inferior stability when it comes to strong turbulence under the Gamma-Gamma model.

 figure: Fig. 4.

Fig. 4. The SNR distribution of C-NLM, C-BM3D, chaotic Gaussian filtering, chaotic Median filtering under different Gamma-Gamma turbulent conditions: (a) a = 5 dB/km, $c_n^2 = 8 \times {10^{ - 15}}{m^{ - 2/3}}$; (b) a = 10 dB/km, $c_n^2 = 8 \times {10^{ - 14}}{m^{ - 2/3}}$; (c) a = 12 dB/km, $c_n^2 = 8 \times {10^{ - 13}}{m^{ - 2/3}}$; (d) a = 15 dB/km, $c_n^2 = 8 \times {10^{ - 12}}{m^{ - 2/3}}$.

Download Full Size | PDF

To further verify the performance of C-BM3D and logistic chaotic mapping, we calculate the BER distribution of the BM3D, NLM, Gaussian filtering, Median filtering and FSO channel with and without chaos, as shown in Fig. 5. Specifically, Fig. 5(a) and Fig. 5(b) demonstrate the results under Log-normal and Gamma-Gamma turbulence models, respectively. It can be seen from Fig. 5 that the BER difference between with and without chaos of every algorithm is little, indicating that logistic chaotic mapping has no adverse effect on the FSO communication quality but also improves communication security. Importantly, the BER values of the BM3D algorithm and FSO channel are almost equal, indicating that the BM3D algorithm improves the image quality without sacrificing the communication quality. From Fig. 5(a), the BM3D algorithm improves at least 24 dB BER performance than the other three algorithms in Log-normal turbulent condition. The superiority of the BM3D algorithm is also verified under the Gamma-Gamma turbulence model. From Fig. 5(b), the BER value of the BM3D algorithm remains unchanged at approximately 0.014 in weak turbulence of the Gamma-Gamma model, improving by at least 15 dB compared with other algorithms. Thus, the BM3D algorithm improves the image quality and communication quality greatly, which can be a good candidate for the next generation of FSO communication in security and turbulence resistance.

 figure: Fig. 5.

Fig. 5. The BER distribution of C-NLM, C-BM3D, chaotic Gaussian filtering, chaotic Median filtering under Log-normal (a) and Gamma-Gamma (b) turbulent conditions.

Download Full Size | PDF

4. Conclusion

In this paper, we compare the communication performance of C-BM3D, C-NLM, chaotic Gaussian filtering and chaotic Median filtering under Log-normal and Gamma-Gamma turbulence models in the FSO system. The results show that the PSNRs of C-BM3D in the weak turbulence under Log-normal and Gamma-Gamma models are up to 96.2956 and 93.2853, respectively. The C-BM3D also achieves superior image similarity in the Log-normal turbulent channel, with its SSIMs nearly equal to 1. Additionally, the SNR of C-BM3D has ranked the highest, and its BER has remained at the lowest among the four algorithms in FSO communication, improving at least 15 dB BER performance compared with other algorithms. Overall, C-BM3D outperforms other methods since it not only recovers the image with superior quality from the chaotic system of the FSO communication system but also enhances communication quality and security. However, one problem is that the C-BM3D lacks stability in strong turbulence, which needs to be further improved and studied in the future.

Funding

National Natural Science Foundation of China (62025105, 62301097, U21B2005); China Postdoctoral Science Foundation (2023T160775); Natural Science Foundation of Chongqing (CSTB2022NSCQ-MSX0542); Chongqing Municipal Education Commission (KJQN202200610).

Disclosures

The authors declare no conflicts of interest.

Data availability

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

References

1. C. Y. Li, X. H. Huang, H. H. Lu, et al., “A WDM PAM4 FSO–UWOC integrated system with a channel capacity of 100 Gb/s,” J. Lightwave Technol. 38(7), 1766–1776 (2020). [CrossRef]  

2. Z. Hu, C. Chen, Z. Zhang, et al., “Secure Cooperative Transmission for Mixed RF/FSO Spectrum Sharing Networks,” IEEE Trans. Commun. 68(5), 3010–3023 (2020). [CrossRef]  

3. S. Malik, P. Saxena, and Y. H. Chung, “Performance analysis of a UAV-based IRS-assisted hybrid RF/FSO link with pointing and phase shift errors,” J. Opt. Commun. Netw. 14(4), 303–315 (2022). [CrossRef]  

4. Z. Zhang, Q. Sun, M. López-Benítez, et al., “Performance Analysis of Dual-Hop RF/FSO Relaying Systems With Imperfect CSI,” IEEE Trans. Veh. Technol. 71(5), 4965–4976 (2022). [CrossRef]  

5. Y. Ma, T. Lv, G. Pan, et al., “On secure uplink transmission in hybrid RF-FSO cooperative satellite-aerial-terrestrial networks,” IEEE Trans. Commun. 70(12), 8244–8257 (2022). [CrossRef]  

6. Q. Huang, M. Lin, W. P. Zhu, et al., “Uplink Massive Access in Mixed RF/FSO Satellite-Aerial-Terrestrial Networks,” IEEE Trans. Commun. 69(4), 2413–2426 (2021). [CrossRef]  

7. E. Zedini, A. Kammoun, and M. S. Alouini, “Performance of multibeam very high throughput satellite systems based on FSO feeder links with HPA nonlinearity,” IEEE Trans. Commun. 19(9), 5908–5923 (2020). [CrossRef]  

8. G. Xu, Q. Zhang, Z. Song, et al., “Relay-Assisted Deep Space Optical Communication System Over Coronal Fading Channels,” IEEE Trans. Aerosp. Electron. Syst. 59(6), 8297–8312 (2023). [CrossRef]  

9. P. Agheli, H. Beyranvand, and M. J. Emadi, “UAV-assisted underwater sensor networks using RF and optical wireless links,” J. Lightwave Technol. 39(22), 7070–7082 (2021). [CrossRef]  

10. L. Qu, G. Xu, Z. Zeng, et al., “UAV-Assisted RF/FSO Relay System for Space-Air-Ground Integrated Network: A Performance Analysis,” IEEE Trans. Wireless Commun. 21(8), 6211–6225 (2022). [CrossRef]  

11. M. K. Miao and X. F. Li, “Performance analysis of FSO systems over a Lognormal-Rician turbulence channel with generalized pointing errors,” J. Lightwave Technol. 40(13), 4206–4216 (2022). [CrossRef]  

12. T. W. Wu, W. Zeng, Y. J. Liu, et al., “Secure turbulence-resistant coherent free-space optical communications via chaotic region-optimized probabilistic constellation shaping,” Opt. Lett. 48(3), 684–687 (2023). [CrossRef]  

13. L. B. Stotts and L. C. Andrews, “Adaptive optics model characterizing turbulence mitigation for free space optical communications link budgets,” Opt. Express 29(13), 20307–20321 (2021). [CrossRef]  

14. Y. T. Huang, H. Z. Huang, H. S. Chen, et al., “Free-space optics communications employing elliptical-aperture multimode diversity reception under anisotropic turbulence,” J. Lightwave Technol. 40(5), 1502–1508 (2022). [CrossRef]  

15. C. Qiu, W. Q. Ma, F. Wang, et al., “Atmospheric turbulence resistant heterodyne coherent receiver of few-mode fiber,” Opt. Express 30(16), 28312–28324 (2022). [CrossRef]  

16. A. K. Singh, “Major development under Gaussian filtering since unscented Kalman filter,” IEEE/CAA J. Autom. Sinica 7(5), 1308–1325 (2020). [CrossRef]  

17. Q. Wang, K. Zhao, X. Yi, et al., “Low-cost OFDR distributed sensing based on optical fiber with enhanced rayleigh backscattering profiles and Median filtering,” in 27th International Conference on OFS (Optica Publishing Group, 2022), paper Th4.48.

18. M. F. Wu, Y. F. Chen, P. B. Zhu, et al., “NLM parameter optimization for ϕ -OTDR signal,” J. Lightwave Technol. 40(17), 6045–6051 (2022). [CrossRef]  

19. Y. H. Do, Y. Cho, S. H. Kang, et al., “Optimization of block-matching and 3D filtering (BM3D) algorithm in brain SPECT imaging using fan beam collimator: Phantom study,” Nucl. Eng. Technol. 54(9), 3403–3414 (2022). [CrossRef]  

20. M. Pan, P. D. Hua, Z. Y. Ding, et al., “Long distance distributed strain sensing in OFDR by BM3D-SAPCA image denoising,” J. Lightwave Technol. 40(24), 7952–7960 (2022). [CrossRef]  

21. J. Y. Yao, Y. H. Yao, F. Y. Cao, et al., “Total variation and block-matching 3D filtering-based image reconstruction for single-shot compressed ultrafast photography,” Opt. Lasers Eng. 139(29), 106475 (2021). [CrossRef]  

22. P. Zdankowski, M. Trusiak, D. McGloin, et al., “Quasi-noise-free stimulated emission depletion microscopy imaging of thick samples using adaptive optics and block-matching 3D filtering,” in Imaging and Applied Optics (Optica Publishing Group, 2019), paper ITh1C.4.

23. J. A. Liang, Y. F. Guo, and B. Liu, “BM3D-based denoising method for color polarization filter array,” Opt. Express 30(12), 22107–22122 (2022). [CrossRef]  

24. A. Yang, Z. Ma, C. Zhang, et al., “Review on application progress of federated learning model and security hazard protection,” Digit. Commun. Netw. 9(1), 146–158 (2023). [CrossRef]  

25. E. Francisca, R. Hurtado-Perez, M. Hernandez-Romo, et al., “Digital image encryption using chaotic logistic mapping,” in Frontiers in Optics (Optica Publishing Group, 2017), paper JTu2A.100.

26. J. Tian, L. Z. Yang, C. Qin, et al., “Refractive index sensing based on chaotic correlation fiber loop ring down system using tapered fiber,” IEEE Sensors J. 20(8), 4215–4220 (2020). [CrossRef]  

27. W. J. Zeng, C. F. Zhang, X. S. Liang, et al., “Chaotic phase noise-like encryption based on geometric shaping for coherent data center interconnections,” Opt. Express 32(2), 1595–1608 (2024). [CrossRef]  

28. Z. Yang, L. L. Yi, J. X. Ke, et al., “Chaotic optical communication over 1000 km transmission by coherent detection,” J. Lightwave Technol. 38(17), 4648–4655 (2020). [CrossRef]  

29. T. W. Wu, C. F. Zhang, Y. H. Chen, et al., “Compressive sensing chaotic encryption algorithms for OFDM-PON data transmission,” Opt. Express 29(3), 3669–3684 (2021). [CrossRef]  

30. Y. Q. Zhang, N. Jiang, A. K. Zhao, et al., “Security enhancement in coherent OFDM optical transmission with chaotic three-dimensional constellation scrambling,” J. Lightwave Technol. 40(12), 3749–3760 (2022). [CrossRef]  

31. Z. Y. Hu, Z. Z. Chen, Y. M. Li, et al., “Adaptive Transceiver Design for High-Capacity Multi-Modal Free-Space Optical Communications With Commercial Devices and Atmospheric Turbulence,” J. Lightwave Technol. 41(11), 3397–3406 (2023). [CrossRef]  

32. S. Malik and P. K. Sahu, “Assessment of the FSO communication system using adaptive and MIMO MPPM with pointing errors and an atmospheric turbulence channel,” Appl. Opt. 60(6), 1719–1728 (2021). [CrossRef]  

33. C. A. N. Santos, D. L. N. Martins, and N. D. A. Mascarenhas, “Ultrasound Image Despeckling Using Stochastic Distance-Based BM3D,” IEEE Trans. on Image Process. 26(6), 2632–2643 (2017). [CrossRef]  

Data availability

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. Initial value sensitivity.
Fig. 2.
Fig. 2. Execution process of BM3D and the reference block is marked with “R”.
Fig. 3.
Fig. 3. The experimental FSO communication system.
Fig. 4.
Fig. 4. The SNR distribution of C-NLM, C-BM3D, chaotic Gaussian filtering, chaotic Median filtering under different Gamma-Gamma turbulent conditions: (a) a = 5 dB/km, $c_n^2 = 8 \times {10^{ - 15}}{m^{ - 2/3}}$; (b) a = 10 dB/km, $c_n^2 = 8 \times {10^{ - 14}}{m^{ - 2/3}}$; (c) a = 12 dB/km, $c_n^2 = 8 \times {10^{ - 13}}{m^{ - 2/3}}$; (d) a = 15 dB/km, $c_n^2 = 8 \times {10^{ - 12}}{m^{ - 2/3}}$.
Fig. 5.
Fig. 5. The BER distribution of C-NLM, C-BM3D, chaotic Gaussian filtering, chaotic Median filtering under Log-normal (a) and Gamma-Gamma (b) turbulent conditions.

Tables (3)

Tables Icon

Table 1. The simulation parameters

Tables Icon

Table 2. Comparison in Log-normal Turbulent Channel

Tables Icon

Table 3. Comparison in Gamma-Gamma Turbulent Channel

Equations (16)

Equations on this page are rendered with MathJax. Learn more.

p ( I ) = 1 2 π σ l 2 1 I exp ( ( I n ( 1 I 0 + σ l 2 2 ) ) 2 2 σ l 2 ) ,
f γ b ( γ b ) = α β ( α + β ) / 2 γ ¯ b ( α + β ) / 4 Γ ( α ) Γ ( β ) γ b ( α + β ) / 4 1 × K α β [ 2 α β γ b γ ¯ b ] ,
α = [ exp ( 0.49 σ R 2 ( 1 + 1.11 σ R 12 / 5 ) 5 / 6 ) 1 ] 1 ,
β = [ exp ( 0.51 σ R 2 ( 1 + 0.69 σ R 12 / 5 ) 5 / 6 ) 1 ] 1 ,
x k + 1 = μ × x k ( 1 x k ) , k = 0 , 1 , , n ,
d ( Z x R , Z x ) = | | γ ( τ 2 D h t ( Z x R ) ) γ ( τ 2 D h t ( Z x ) ) | | 2 2 ( N 1 h t ) 2 ,
S x R h t = { x X : d ( Z x R , Z x ) τ m a t c h h t } ,
Y ^ S x R h t h t = τ 3 D h t 1 ( γ ( τ 3 D h t ( Z S x R h t ) ) ) ,
y ^ b a s i c ( x ) = x R X x m S x R h t ω x R h t Y ^ x m h t , x R ( x ) x R X x m S x R h t ω x R h t χ x m ( x ) , x X .
S x R w i e = { x X : | | Y ^ x R b a s i c Y ^ x b a s i c | | 2 2 ( N 1 w i e ) 2 < τ m a t c h w i e } ,
Y ^ w i e S x R w i e = τ 3 D w i e 1 ( W S x R w i e τ 3 D w i e ( Z S x R w i e ) ) ,
y ^ f i n a l ( x ) = x R X x m S x R w i e ω x R w i e Y ^ x m w i e , x R ( x ) x R X x m S x R w i e ω x R w i e χ x m ( x ) , x X .
S N R = 10 log 10 ( P s / P n )
P b = I e / I ,
P S N R ( f , g ) = 10 lg ( 2 n 1 ) 2 i j ( f ( i , j ) g ( i , j ) ) 2 M N ,
S S I M ( f , g ) = ( 2 μ f μ g + c 1 ) ( 2 σ f g + c 2 ) ( μ f 2 + μ g 2 + c 1 ) ( σ f 2 + σ g 2 + c 2 ) ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.