
In order to generate high-quality single-photon emission computed tomography (SPECT) images under low-dose acquisition mode, a sinogram denoising method was studied for suppressing random oscillation and enhancing contrast in the projection domain. A conditional generative adversarial network with cross-domain regularization (CGAN-CDR) is proposed for low-dose SPECT sinogram restoration. The generator stepwise extracts multiscale sinusoidal features from a low-dose sinogram, which are then rebuilt into a restored sinogram. Long skip connections are introduced into the generator, so that the low-level features can be better shared and reused, and the spatial and angular sinogram information can be better recovered. A patch discriminator is employed to capture detailed sinusoidal features within sinogram patches; thereby, detailed features in local receptive fields can be effectively characterized. Meanwhile, a cross-domain regularization is developed in both the projection and image domains. Projection-domain regularization directly constrains the generator via penalizing the difference between generated and label sinograms. Image-domain regularization imposes a similarity constraint on the reconstructed images, which can ameliorate the issue of ill-posedness and serves as an indirect constraint on the generator. By adversarial learning, the CGAN-CDR model can achieve high-quality sinogram restoration. Finally, the preconditioned alternating projection algorithm with total variation regularization is adopted for image reconstruction. Extensive numerical experiments show that the proposed model exhibits good performance in low-dose sinogram restoration. From visual analysis, CGAN-CDR performs well in terms of noise and artifact suppression, contrast enhancement and structure preservation, particularly in low-contrast regions. From quantitative analysis, CGAN-CDR has obtained superior results in both global and local image quality metrics. From robustness analysis, CGAN-CDR can better recover the detailed bone structure of the reconstructed image for a higher-noise sinogram. This work demonstrates the feasibility and effectiveness of CGAN-CDR in low-dose SPECT sinogram restoration. CGAN-CDR can yield significant quality improvement in both projection and image domains, which enables potential applications of the proposed method in real low-dose study.
Citation: Si Li, Limei Peng, Fenghuan Li, Zengguo Liang. Low-dose sinogram restoration enabled by conditional GAN with cross-domain regularization in SPECT imaging[J]. Mathematical Biosciences and Engineering, 2023, 20(6): 9728-9758. doi: 10.3934/mbe.2023427
[1] | Robert Stephen Cantrell, Chris Cosner, Yuan Lou . Evolution of dispersal and the ideal free distribution. Mathematical Biosciences and Engineering, 2010, 7(1): 17-36. doi: 10.3934/mbe.2010.7.17 |
[2] | Minus van Baalen, Atsushi Yamauchi . Competition for resources may reinforce the evolution of altruism in spatially structured populations. Mathematical Biosciences and Engineering, 2019, 16(5): 3694-3717. doi: 10.3934/mbe.2019183 |
[3] | Jinyu Wei, Bin Liu . Coexistence in a competition-diffusion-advection system with equal amount of total resources. Mathematical Biosciences and Engineering, 2021, 18(4): 3543-3558. doi: 10.3934/mbe.2021178 |
[4] | Andrea Pugliese, Abba B. Gumel, Fabio A. Milner, Jorge X. Velasco-Hernandez . Sex-biased prevalence in infections with heterosexual, direct, and vector-mediated transmission: a theoretical analysis. Mathematical Biosciences and Engineering, 2018, 15(1): 125-140. doi: 10.3934/mbe.2018005 |
[5] | Nancy Azer, P. van den Driessche . Competition and Dispersal Delays in Patchy Environments. Mathematical Biosciences and Engineering, 2006, 3(2): 283-296. doi: 10.3934/mbe.2006.3.283 |
[6] | Zhilan Feng, Robert Swihart, Yingfei Yi, Huaiping Zhu . Coexistence in a metapopulation model with explicit local dynamics. Mathematical Biosciences and Engineering, 2004, 1(1): 131-145. doi: 10.3934/mbe.2004.1.131 |
[7] | Yang Kuang, Angela Peace, Hao Wang . Special issue: Resource explicit population models. Mathematical Biosciences and Engineering, 2019, 16(1): 538-540. doi: 10.3934/mbe.2019025 |
[8] | Bethan Morris, Lee Curtin, Andrea Hawkins-Daarud, Matthew E. Hubbard, Ruman Rahman, Stuart J. Smith, Dorothee Auer, Nhan L. Tran, Leland S. Hu, Jennifer M. Eschbacher, Kris A. Smith, Ashley Stokes, Kristin R. Swanson, Markus R. Owen . Identifying the spatial and temporal dynamics of molecularly-distinct glioblastoma sub-populations. Mathematical Biosciences and Engineering, 2020, 17(5): 4905-4941. doi: 10.3934/mbe.2020267 |
[9] | Azmy S. Ackleh, Shuhua Hu . Comparison between stochastic and deterministic selection-mutation models. Mathematical Biosciences and Engineering, 2007, 4(2): 133-157. doi: 10.3934/mbe.2007.4.133 |
[10] | Muntaser Safan . Mathematical analysis of an SIR respiratory infection model with sex and gender disparity: special reference to influenza A. Mathematical Biosciences and Engineering, 2019, 16(4): 2613-2649. doi: 10.3934/mbe.2019131 |
In order to generate high-quality single-photon emission computed tomography (SPECT) images under low-dose acquisition mode, a sinogram denoising method was studied for suppressing random oscillation and enhancing contrast in the projection domain. A conditional generative adversarial network with cross-domain regularization (CGAN-CDR) is proposed for low-dose SPECT sinogram restoration. The generator stepwise extracts multiscale sinusoidal features from a low-dose sinogram, which are then rebuilt into a restored sinogram. Long skip connections are introduced into the generator, so that the low-level features can be better shared and reused, and the spatial and angular sinogram information can be better recovered. A patch discriminator is employed to capture detailed sinusoidal features within sinogram patches; thereby, detailed features in local receptive fields can be effectively characterized. Meanwhile, a cross-domain regularization is developed in both the projection and image domains. Projection-domain regularization directly constrains the generator via penalizing the difference between generated and label sinograms. Image-domain regularization imposes a similarity constraint on the reconstructed images, which can ameliorate the issue of ill-posedness and serves as an indirect constraint on the generator. By adversarial learning, the CGAN-CDR model can achieve high-quality sinogram restoration. Finally, the preconditioned alternating projection algorithm with total variation regularization is adopted for image reconstruction. Extensive numerical experiments show that the proposed model exhibits good performance in low-dose sinogram restoration. From visual analysis, CGAN-CDR performs well in terms of noise and artifact suppression, contrast enhancement and structure preservation, particularly in low-contrast regions. From quantitative analysis, CGAN-CDR has obtained superior results in both global and local image quality metrics. From robustness analysis, CGAN-CDR can better recover the detailed bone structure of the reconstructed image for a higher-noise sinogram. This work demonstrates the feasibility and effectiveness of CGAN-CDR in low-dose SPECT sinogram restoration. CGAN-CDR can yield significant quality improvement in both projection and image domains, which enables potential applications of the proposed method in real low-dose study.
In recent years, multi-agent consensus control has been recognized as a vital element of distributed collaborative control for applications such as distributed computing, unmanned aerial vehicle formation, and intelligent transportation systems. Researchers have shown significant interest in this area, and a wide range of control mechanisms have been explored in the past few years. These control mechanisms include adaptive [1,2,3,4], fault-tolerant [5], impulse [6,7], and sliding mode [8] methods.
In practical systems, the stability of the system state is susceptible to disruption from unexpected factors, including nonlinear dynamics and system uncertainties. Existing research has primarily focused on continuous-time systems, where the state information of intelligent agents is continuously transmitted between nodes, leading to significant network usage and energy consumption. However, the development of event-triggered control solves this problem by avoiding constant communication. Earlier studies on event-triggered control that can be found have investigated centralized, distributed, and self-triggered event-triggered control techniques [9,10,11]. One researcher proposed an adaptive event-triggered control scheme for strongly connected networks that dynamically adjusted the triggering time interval on the basis of sampled data [12]. Another mechanism using a dynamic event-triggering mechanism was proposed to reduce communication resource wastage compared with traditional event-triggering mechanisms [13]. Others have assumed that system parameters such as the efficiency factor of the executor, external disturbances, and precursor control input signals are all unknown and introduced a fault-tolerant control to obtain sufficient conditions for consistent tracking [14]. However, these findings primarily investigate traditional event triggering mechanisms. In [15,16], researchers proposed a dynamic event-triggered mechanism, which can significantly reduce the number of triggers and conserve communication resources. researchers respectively proposed centralized and distributed dynamic event-triggered mechanisms in [15,16], while scholars suggested both centralized and distributed mechanisms in [17], verifying their superiority.
However, the rapid growth of network information technologies has also led to a rise in cyber attacks. Among them, denial-of-service (DoS) attacks are the most common, being relatively easy to execute in the attack space. These attacks usually target the controller or exhaust the resources of the target system directly, resulting in the system being unable to provide normal services or communication. In some cases, these attacks cause the system to crash. Therefore, countering DoS attacks has received significant research attention. Researchers studied the multi-agent systems(MASs) under DoS attacks in given attack frequency and upper bounds on attack duration [18,19]. Compared with linear systems, nonlinear MASs are more widely used in real life. Among these, a secure controller based on event triggering was proposed to solve the lead-following consensus problem of second-order nonlinear systems [20]. This is more common than linear systems. Another proposal was for an event-triggered adaptive fault-tolerant control strategy, which reduced the computational cost of heterogeneity [21]. For nonperiodic DoS attacks, the upper bounds of network attacks, actuator failures, attack duration, and frequency are obtained. Another method uses a security mechanism employing a prediction-based switching observer scheme to address the issue of invalidation in event-triggered mechanisms during attack intervals [22]. A novel framework for observer-based event-triggered containment control, taking into account the occurrence of DoS attacks, has also been introduced [23]. This framework establishes a resilient event-triggered controller, using a specially designed observer. The goal is to achieve consistent control of MASs in the presence of DoS attacks.
Based on these observations, we aim to explore the security consensus problem of nonlinear MAS with external disturbances under DoS attacks in this paper. Our contributions are as follows.
1) In this paper, a nonlinear system with external disturbances is considered, and the effects of the nonlinear dynamics and uncertainty of the system are eliminated by designing an adaptive scheme and state-feedback control gains by updating the laws of the adaptive parameters online.
2) Compared with [22,23], a dynamic variable is introduced to adjust the triggering instances under DoS attacks. Therefore, the event-triggered mechanism proposed in this paper is more flexible and can effectively save communication resources. In addition, continuous communication between agents is not required to determine whether a trigger condition satisfies the trigger condition.
Notation R is the set of real numbers, and RN×N is the set of N×N real matrix. ‖⋅‖ represents a Euclidean norm of vectors or matrices. The superscripts A−1 and AT represent the inverse and transpose of matrix A. λmax(A) is the maximum eigenvalue, and λmin(A) is the minimum eigenvalue of matrix A. D+(⋅) denotes the righthand derivative of a function, and ⨂ is Kronecker product. diag{A1,⋯,An} is the diagonal matrix. ∩ is the intersection of sets, and ∪ denotes the union of sets.
For a given MAS, the digraph G is (V,E), where V={1,2,⋯,N} is the set of nodes, and E⊆V×V represents the edge set of followers. The information exchange between each node can be described by the adjacency matrix A and the Laplacian matrix L. A=[aij]∈RN×N if agents i and j communicate with one another, aij=1; otherwise aij=0 and L=[lij]∈Rn×n where L=D−A, The degree matrix D=diag(di) with di=∑Ni=1aij. In this paper, we assume that the agents are linked by a balancing topology, i.e., aij=aji. If the agent i communicates with the leader, then bi=1; otherwise, bi=0.
For a leader-following system, the dynamics of the leader are described by the equation
˙x0=Ax0+f(t,x0(t)). | (2.1) |
The ith follower system is
˙xi=Axi+Bui(t)+f(t,xi(t))+wi. | (2.2) |
In the preceding, x(t)∈Rn are positions of the agent, ui(t)∈R is the control input, A∈Rn×n and B∈Rn×p are system matrices, f(x) is a nonlinear function, and wi is the uncertainty input satisfying
wi≤ςi1|u|+ςi2|x|+γi, | (2.3) |
where ςi1<1, ςi2 and γi are unknown constants.
Lemma 1. If the nonlinear function f(t,xi(t)),i=1,2,..., is continuously differentiable in a region S∈R2 and xi(t0)∈S, then for any xi(t0)∈S, the following formula is satisfied:
f(xi(t),t)−f(xi(t0),t)=∂f(⋅)∂xi×(xi(t)−xi(t0)), | (2.4) |
where f(⋅)=f(xi(t0))+Δ(xi(t)−xi(t0)), 0<Δ<1.
Assumption 1. If there is a continuously differentiable function f(t,xi(t)) and the highest order s∈{1,2,⋯N}, there exist bounded positive scalars δix, such that
|∂f(zi)∂xi|≤δixxis, | (2.5) |
where xis=|xi|s+|xi|s−1+⋯+1,s≥1. We also need some assumptions to ensure that the purpose is achieved.
Assumption 2. A,B can be stabilized, and the digraph G is strongly connected.
Next, we define the position errors ei(t):
ei(t)=N∑j=0,j≠iaij(xi(t)−xj(t))+bi(xi(t)−x0(t)). | (2.6) |
According to Definition 1, we have
˙ei(t)=N∑j=0,j≠iaij(˙xi(t)−˙xj(t))+bi(˙xi(t)−˙x0(t))=N∑j=0,j≠iaij(A(xi−xj)+B(ui−uj)+[f(t,xi(t))−f(t,xj(t))]+(wi−wj))+bi(A(xi−x0)+Bui)+wi+bi[f(t,xi(t))−f(t,x0(t))]=N∑j=0,j≠iAaij((xi−xj)+bi(xi−x0))+Baij(ui−uj)+aij∂f(zi)∂xi(xi(t)−xj(t0))+aij(wi−wj))+bi(Bui+wi))+bi[∂f(zi)∂xi(xi(t)−x0(t))]=N∑j=0,j≠iA(aij(xi−xj)+bi(xi−x0))+∂f(zi)∂xi[aij(xi(t)−xj(t0))+bi(xi(t)−x0(t))]+aij(B(ui−uj))+(wi−wj))+bi(Bui+wi)). | (2.7) |
We also have
˙e(t)=(L+L0)(Bu(t)+w(t))+(A+∂f(zi)∂x)e(t), | (2.8) |
where e(t)=[e1(t),e2(t),⋯,ei(t)], ∂f(zi)∂x=diagNi=1[∂f(zi)∂xi], L0=diagNi=1[bi]. L is defined as
lij={∑Nk=1,k≠iaik,j=i−aij,j≠i. | (2.9) |
, then Eq (2.8) can be expressed as
˙e(t)=(A+ΔA)e(t)+ˉL⊗(Bu(t)+w(t)), | (2.10) |
where A=[IN×N000], ˉL=[0L+L0], ΔA=[000∂f(zi)∂x].
Definition 1. MAS (2.1) and (2.2) are said to have consensus if each agent's position state in the system satisfies
limt→∞‖xi(t)−x0(t)‖=0,i=1,2,⋯,N. | (2.11) |
The distributed adaptive control input is
ui(t)=kiei(t)−ˉeiΨi(t), | (2.12) |
where ˉei=∑Nj=1ajieji,aii=ai0+∑Nj=1,j≠iaij. ki is the control gain determined based on the linear matrix inequality(LMI)
(A+ˉL⊗BK)TP+P(A+ˉL⊗BK)≤0, | (2.13) |
where P=[P00N×N0N×NIN×N]>0,K=diag[ki],andP0 is a positive definition matrix. Ψi(t) is defined as
Ψi={11−ςi1(ςi1|kiei|+ςi2|xis|+γi|ei|+2ˆδix|e|2xis|ei|2),ei≠00,ei=0 | (2.14) |
where ˆδix is the estimate of unknown parameters δix. The following describes the updated laws for the adaptive parameters:
dˆδidt=|ei|2xis. | (2.15) |
Since δi is an unknown constant, defined as ˜δi(t)=ˆδi(t)−δi, the adaptive error systems are described by
d˜δidt=dˆδidt. | (2.16) |
It follows from Eq (2.12) that
u(t)=ke(t)−ˉeΨ(t), | (2.17) |
where k(t)=[k1(t),k2(t),⋯]T and ˉe=[ˉe1,ˉe2⋯ˉen]. According to Eqs (2.11) and (2.15), we have
˙e(t)=(A+ˉL⊗BK)e(t)+ΔAe(t)−ˉL(BˉevΨ(t))−ˉLw(t). | (2.18) |
Next, we define the event trigger time series as {tjk} for the jth agent. Therefore, the next triggering time tik+1 for the ith agent can be expressed as
tik+1=inf{t>tik|Hi(t)≥0}. | (2.19) |
The function Hi(⋅) is given by
Hi(⋅)=−θiχi(t)+αi‖qi(t)‖2−ηi‖ei(t)‖2, | (2.20) |
where θi>0, αi,ηi∈Rn>0. qi(t) is defined as the measurement error according to Eq (2.6):
qi(t)=ei(tik)−ei(t). | (2.21) |
χi(t) satisfies
˙χi(t)=−βiχi(t)+ηi‖ei(t)‖2−αi‖qi(t)‖2, | (2.22) |
where βi>0, initial value χi(0)>0 could be randomly selected, and ti0=0.
Remark 1: The internal dynamic variable updates according to internal variables such as self-feedback, measurement error, and neighborhood error. In comparison with the conventional static triggering strategy [22,23], the dynamic event triggered control protocol we proposed can more effectively reduce network communication and save resources.
A DoS attack aims to block the communication channels so the targeted system cannot exchange information normally. Communication channels are not the only things affected by DoS attacks because the attacks can damage communication equipment along with hindering data transmission, measurement, and control channels simultaneously. DoS attacks are extremely destructive to the system, but their energy consumption requires attackers to replenish energy supplies after the attack is over, which takes time. Therefore, the time series can be split into two sections based on whether a DoS assault was launched. In the absence of the DoS attack, the system functions and communicates properly. However, in the presence of a DoS attack, communication is cut off, and the controller stops functioning. Here, we assume that the time interval of DoS attacks is {tm}m∈N, where tm is the moment of the DoS attack, and [tm,tm+Δm] is the mth DoS time interval, and Δm is the time duration of the mth attack. The DoS attack interval is the same for all multi-agents. Thus, the set time instants where communication is blocked (the interval of the DoS attack) are
Ξa(t0,t)={⋃m∈N[tm,tm+Δm]}∩[t0,t]. | (3.1) |
Similarly, the sequence of time intervals without attacks is given by
Ξs(t0,t)=[t0,t]∖Ξa(t0,t). | (3.2) |
Because of the recovery mechanism, the MAS cannot immediately restore communication after the end of a DoS attack, and due to the event-triggering mechanism, there is an upper bound for the time when the two events occur consecutively. We assume that they can exist at the same time. Therefore, the actual DoS attack lasts longer, and consequently, the mth DoS attack's actual time frame may be described as [tm,tm+ˉΔm]. The new time period of the DoS attack is as follows:
˜Ξa(t0,t)={⋃m∈N[tm,tm+ˉΔm]}∩[t0,t] | (3.3) |
˜Ξs(t0,t)=[t0,t]∖˜Ξa(t0,t) | (3.4) |
Assumption 2. Define na(t0,t) as the number of attacks in the period [t0,t], so the attack frequency Fa(t0,t)>0 satisfies
Fa(t0,t)=na(t0,t)t−t0. | (3.5) |
Assumption 3. Define Na(t0,t) as the total time interval of the DoS attack in the period [t0,t]. The constants T0≥0,F0≥0,0<1T1<1,0<1F1<1 are such that
|Ξa(t0,t)|≤Ξ0+t−t0T1and | (3.6) |
Na(t0,t)≤F0+t−t0F1, | (3.7) |
where 1T1 is the attack strength.
Lemma 1. Previous research considers Eq (2.1) and this DoS attack model under Assumptions 2 and 3 [18]. If the Lyapunov function V1(t),V2(t) satisfies
{˙V1(t)≤−l0V(t)+τ0t∈˜Ξs˙V2(t)≤l1V(t)+τ1t∈˜Ξa, | (3.8) |
where l0,l1,τ0,τ1 are positive constants. T1,F1 defined in Assumption 3 satisfies
1T1<l0−η∗l0+l11F1<η∗2lnκ+(l0+l1)ρ, | (3.9) |
where 0<η∗<l0 is the time to restore communication. ρ>0, κ≥1 is a constant satisfying
{κV2((tm+ˉΔm)−)−V1(tm+ˉΔm)≥0κV1(t−m)−V2(tm+1)≥0. | (3.10) |
Thus, we say that V(t) are bounded.
Remark 2: Lemma 1 gives an upper bound on DoS attack frequency and duration, ensuring that the Lyapunov function remains stable over the entire time span [18].
Remark 3: The DoS attack considered in this paper mainly attacks the communication channels between agents. Thus, when the DoS attack comes, there is no information interaction between neighboring agents, and the event-triggering control is not triggered. In addition, we consider a DETC. Compared with the traditional event-triggering control, we introduce a dynamic variable that uses communication resources more effectively. In the simulation section below, we compare our method with the traditional event-triggering mechanism.
In this section, we prove system stability. Our presentation has two sections: the stability study of the MAS (2.1) and (2.2) under a DoS attack and the proof of non-Zeno behavior.
Theorem 1. For the MAS (2.1) and (2.2) under DoS attacks, we consider Assumption 1 and the controller (2.12). If the LMI (2.13) satisfies (A+ˉL⊗BK)TP+P(A+ˉL⊗BK)≤ξiP, where ξi∈Rn=σiηi,σi>1, then a feasible solution exists and the MAS is said to achieve leader-following consensus.
Proof of Theorem 1. The system stability proof is also divided into two parts. The communication of the system is damaged under a DoS attack, but the system is not always in an impassable state. The proof is divided between DoS attacks and non-DoS attacks, as per the prior section. When there are non-DoS attacks in the system, we consider the Lyapunov function
W(t)=V(t)+N∑i=1χi(t)=eT(t)Pe(t)+N∑i=1κ−1i~δ2i+N∑i=1χi(t). | (4.1) |
It follows from Eqs (2.20)–(2.22) that
˙χi=−βiχi−θiχi. | (4.2) |
The preceding implies that
χi(t)≥χi(0)e−(βi+θi)t>0, | (4.3) |
which leads to W(t)>0.
The derivative of W is
˙W(t)=eT[(A+ˉL⊗BK)TP+P(A+ˉL⊗BK)]e+2eTPΔAe−2eTPˉLˉeΨ(t)+2eTPˉLw(t)+N∑i2κ−1i˙˜δi˜δi+N∑i=1˙χi(t). | (4.4) |
According to Eq (2.3), the condition in Assumption 1, and the control protocol in {Eq (2.12)}, we have
˙W(t)≤eT[(A+ˉL⊗BK)TP+P(A+ˉL⊗BK)]e+2N∑i=1|e|2∂f(zi)∂xi−2N∑i=1|e|2Ψi+2N∑i=1|e|(ςi1|u|+ςi2|x|+δi)+N∑i2κ−1i˙˜δi˜δi+N∑i=1˙χi(t)≤eT[(A+ˉL⊗BK)TP+P(A+ˉL⊗BK)]e+2N∑i=1δix|e|2xis−2N∑i=1|e|2(1−ςi1)Ψi+2N∑i=1|e|(ςi1|kiei|+ςi2|x|+δi)+N∑i2κ−1i˙˜δi˜δi+N∑i=1˙χi(t). | (4.5) |
Choosing Ψi(t) as in Eq (2.15), we obtain
˙W(t)≤eT[(A+ˉL⊗BK)TP+P(A+ˉL⊗BK)]e−2N∑i=1δix|e|2xis+N∑i2κ−1i˙˜δi˜δi+N∑i=1˙χi(t)≤eT[(A+ˉL⊗BK)TP+P(A+ˉL⊗BK)]e+N∑i=1˙χi(t). | (4.6) |
On the other hand
N∑i=1ξieTi(t)Pei(t)≤λmax(P)N∑i=1ξi‖ei(t)‖2. | (4.7) |
Based on the condition in Theorem 1, Eq (2.22), and Eq (4.6), we have
˙W(t)≤−N∑i=1(ξi−ηi)‖ei(t)‖2−N∑i=1αi‖qi(t)‖2−N∑i=1βiχi≤−N∑i=1ηi(σi−1)‖ei(t)‖2−N∑i=1βiχi, | (4.8) |
then
˙W(t)≤−(σM−1)N∑i=1ξi‖ei(t)‖2−N∑i=1βiχi≤−l0W(t)+τ0, | (4.9) |
where l0=min([(σM−1)/λmax(P)],1, βm)>0, σM=max[σi], βm=min[βi], τ0=∑Ni=1κ−1i~δ2i, l0,andτ0 are positive constants.
When there are DoS attacks in the system, then communication and control channel blockages exist. In this case, the control input becomes 0, ui(t)=0, so the Lyapunov function can be expressed as
V(t)=eT(t)Pe(t)+N∑i=1κ−1i~δ2i. | (4.10) |
Similar to (4.4), (4.10) can be written as
˙V(t)≤V(t)+2N∑i=1δix|e|2xis+2N∑i=1|e|(ςi2|x|+δi)≤l1V(t)+τ1, | (4.11) |
where l1=1, and τ1=2∑Ni=1δix|e|2xis+2∑Ni=1|e|(ςi2|x|+δi). According to the conditions of (3) and Assumption 1, we know that τ1 has an upper bound. From Lemma 1, we know that the system stabilizes in a limited time under a DoS attack. The proof is completed.
Next is the proof of no Zeno behavior. We assume that there is a positive constant T0 such that limk→∞tik=T0. Based on the property of limit, we know that for any ε0>0, there exists N(ε0) such that tik∈(T0−ε0,T0+ε0),∀k≥N(ε0). This means tiN(ε0+1)−tiN(ε0)<2ε0.
According to (4.11), W(t) gradually decreases to 0, Then ξmλmin(P)‖ei(t)‖2≤V(t)<W(t). Therefore, we have
‖ei(t)‖≤√W0ξmλmin(P)=ϖ0. | (4.12) |
Because ‖ei(t)‖ and ‖qi(t)‖ are bounded, the Dini derivative of ‖qi(t)‖ is
D+‖qi(t)‖≤‖˙qi(t)‖=‖−N∑j=1aij(˙xi(t)−˙xj(t))+bi(˙xi(t)−˙x0(t))‖≤‖A+ΔA‖‖ei(t)‖+‖ˉL‖‖B‖‖N∑j=1(ui(t))‖+‖ˉL‖‖wi(t)‖≤‖ˉA‖ϖ0+‖ˉL‖‖B‖M1+‖L‖M2=ˆW0, | (4.13) |
where ˉA=A+ΔA. According to Eqs (2.3), (2.12), and (2.14), we obtain ui(t). wi(t) has an upper bound, and M1,M2 is their upper bound.
Since only the trigger condition in {Eq (2.19)} is met and the event is triggered when ‖qi(t)‖ is reset to 0, then ‖qi(t)‖≥√ηiαi‖ei(t)‖2+θiαiχi≥√θiαiχi,ti−k,k=1,2⋯, which implies that
‖qi(ti−k)‖≥√θiαiχi(ti−k)=√θiαiχi(0)e−βi+θi2ti−k, | (4.14) |
then, we can obtain
tiN(ε0+1)−tiN(ε0)≥1ˆW0√θiαiχi(0)e−βi+θ12ti−N(ε0+1). | (4.15) |
If ε0>0 is a solution of
1ˆW0√θiαiχi(0)e−βi+θi2TO=2ε0eβi+θi2ε0, | (4.16) |
then
tiN(ε0+1)−tiN(ε0)≥1ˆW0√θiαiχi(0)e−βi+θi2(T0+ε0)=2ε0. | (4.17) |
As a result, the aforementioned assumption is false, concluding the evidence that the agent i does not have Zeno behavior.
To show the efficacy of the proposed control strategy, we present a simulation example in this section. Our simulation uses MASs composed of six agents as shown in Figure 1, where agent 1 is the leader, and others are followers. The system is
˙xi(t)=Axi(t)+Bui(t)+(−sin(xi(t))+1.5cos(2.5∗t))+wi |
The system parameters are set as
A=[0I3A1A2],B=[0I3] |
A1=[00002ϕ2000−ϕ2],A2=[02ϕ0−2ϕ00000] |
ςi1=0.1×√1+i2,ςi2=0.5×√1+i2 |
In this example, we consider the flight of an aircraft, ϕ=0.002 is the angular velocity of the aircraft, and I3 represents the identity matrix of 3 × 3.
αi=87.5,βi=0.004,θi=3.5 |
ηi=[0.210.1050.1050.210.210.105] |
Figures 2 and 3 show the response curves and consistency errors of the system state for all agents. They show that the followers' states converge toward those of the leader as time progresses. Figures 3 and 4 show the control input curves and event trigger time instant for all agents. There are four times DoS attacks, with T0=3,F0=4. The duration of the DoS attack is |Ξa(0,40)|=3.5. In Table 1, we can see that the dynamic event-triggered mechanism proposed in this paper has far fewer triggering instances in the same time than the other two literatures [22,23], which can effectively save communication resources. In addition, continuous communication between agents is not required to determine whether a trigger condition satisfies a trigger condition. Considering the static event-triggered control protocol, we have
tik+1=inf{t>tik|‖qi(t)‖2−ρ‖ei(t)‖2≥0} |
tik+1=inf{t>tik|‖qi(t)‖−ϱ‖ei(t)‖≥0} |
Agent i | 1 | 2 | 3 | 4 | 5 | 6 |
[22] | 260 | 1897 | 1914 | 1704 | 1635 | 1861 |
[23] | 1632 | 1899 | 1917 | 1694 | 1633 | 1869 |
Our DETC | 50 | 65 | 70 | 106 | 42 | 6 |
where ρ and ϱ are positive constants. Our DETC effectively reduces communication frequency.
In this paper, we propose a dynamic event-triggered adaptive control approach to address the leader-following consensus problem for nonlinear MASs experiencing DoS attacks. We have presented a distributed control strategy and adaptive update laws to ensure system stability in the presence of uncertainties. The Lyapunov stability theory is used to derive conditions for achieving consensus. The DoS attacks considered here mainly target the MASs' communication channels. In reality, there are other types, scales, and levels of DoS attacks. Formulating mathematical models of these other types of DoS attacks and solving these models is the direction of our future research.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This research was funded in part by National Natural Science Foundation of China under Grant 62273109, in part by the Guangdong Basic and Applied Basic Research Foundation under Grant 2023A1515010168, Grant 2019A1515010830, in part by the Key Special Foundation for General Universities in Guangdong Province under Grant 2022ZDZX1018, and in part by the Maoming Science and Technology Plan Foundation under Grant 2022S043.
The authors declare there is no conflict of interest.
[1] | H. A. Ziessman, J. P. O'Malley, J. H. Thrall, Nuclear Medicine: The Requisites, Elsevier, 2014. |
[2] |
P. Ritt, H. Vija, J. Hornegger, T. Kuwert, Absolute quantification in SPECT, Eur. J. Nucl. Med. Mol. Imaging, 38 (2011), 69–77. https://doi.org/10.1007/s00259-011-1770-8 doi: 10.1007/s00259-011-1770-8
![]() |
[3] |
X. Niu, Y. Yang, M. Jin, M. N. Wernick, M. A. King, Effects of motion, attenuation, and scatter corrections on gated cardiac SPECT reconstruction, Med. Phys., 38 (2011), 6571–6584. https://doi.org/10.1118/1.3660328 doi: 10.1118/1.3660328
![]() |
[4] |
R. G. Wells, Dose reduction is good but it is image quality that matters, J. Nucl. Cardiol., 1 (2018), 1–3. https://doi.org/10.1007/s12350-018-1378-5 doi: 10.1007/s12350-018-1378-5
![]() |
[5] |
J. Zhang, S. Li, A. Krol, C. R. Schmidtlein, E. Lipson, D. Feiglin, et al., Infimal convolution-based regularization for SPECT reconstruction, Med. Phys., 45 (2018), 5397–5410. https://doi.org/10.1002/mp.13226 doi: 10.1002/mp.13226
![]() |
[6] |
A. Krol, S. Li, L. Shen, Y. Xu, Preconditioned alternating projection algorithms for maximum a posteriori ECT reconstruction, Inverse Probl., 28 (2012), 115005. https://doi.org/10.1088/0266-5611/28/11/115005 doi: 10.1088/0266-5611/28/11/115005
![]() |
[7] |
Y. Luo, M. Wei, S. Li, J. Ling, G. Xie, S. Yao, An effective co-support guided analysis model for multi-contrast MRI reconstruction, IEEE J. Biomed. Health, (2023). https://doi.org/10.1109/JBHI.2023.3244669 doi: 10.1109/JBHI.2023.3244669
![]() |
[8] |
H. Zhang, B. Dong, A review on deep learning in medical image reconstruction, J. Oper. Res. Soc. China, 8 (2020), 311–340. https://doi.org/10.1007/s40305-019-00287-4 doi: 10.1007/s40305-019-00287-4
![]() |
[9] |
I. Häggström, C. R. Schmidtlein, G. Campanella, T. J. Fuchs, DeepPET: A deep encoder-decoder network for directly solving the PET reconstruction inverse problem, Med. Image Anal., 54 (2019), 253–262. https://doi.org/10.1016/j.media.2019.03.013 doi: 10.1016/j.media.2019.03.013
![]() |
[10] |
B. Zhu, J. Z. Liu, S. F. Cauley, B. R. Rosen, M. S. Rosen, Image reconstruction by domain-transform manifold learning, Nature, 555 (2018), 487–492. https://doi.org/10.1038/nature25988 doi: 10.1038/nature25988
![]() |
[11] | H. Zhang, B. Dong, B. Liu, JSR-Net: A deep network for Joint Spatial-Radon domain CT reconstruction from incomplete data, in IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), (2019), 3657–3661. https://doi.org/10.1109/ICASSP.2019.8682178 |
[12] | Y. Yang, H. Li, Z. Xu, J. Sun, Deep ADMM-Net for compressive sensing MRI, Adv. Neural Inf. Process. Syst., 29 (2016), 10–18. |
[13] |
Y. Yang, H. Li, Z. Xu, ADMM-CSNet: A deep learning approach for image compressive sensing, IEEE Trans. Pattern Anal. Mach. Intell., 42 (2018), 521–538. https://doi.org/10.1109/TPAMI.2018.2883941 doi: 10.1109/TPAMI.2018.2883941
![]() |
[14] |
J. Adler, O. Öktem, Learned primal-dual reconstruction, IEEE Trans. Med. Imaging, 37 (2018), 1322–1332. https://doi.org/10.1109/TMI.2018.2799231 doi: 10.1109/TMI.2018.2799231
![]() |
[15] |
B. Zhou, X. Chen, S. K. Zhou, J. S. Duncan, C. Liu. DuDoDR-Net: Dual-domain data consistent recurrent network for simultaneous sparse view and metal artifact reduction in computed tomography, Med. Image Anal., 75 (2022), 102289. https://doi.org/10.1016/j.media.2021.102289 doi: 10.1016/j.media.2021.102289
![]() |
[16] |
M. Li, W. Hsu, X. Xie, J. Cong, W. Gao. SACNN: Self-attention convolutional neural network for low-dose CT denoising with self-supervised perceptual loss network, IEEE Trans. Med. Imaging, 39 (2020), 2289–2301. https://doi.org/10.1109/TMI.2020.2968472 doi: 10.1109/TMI.2020.2968472
![]() |
[17] | W. Bae, J. J. Yoo, J. C. Ye, Beyond deep residual learning for image restoration: persistent homology-guided manifold simplification, in IEEE Conference on Computer Vision and Pattern Recognition (CVPR) Workshops, (2017), 145–153. https://doi.org/10.1109/CVPRW.2017.152 |
[18] |
H. Chen, Y. Zhang, M. K. Kalra, F. Lin, Y. Chen, P. Liao, et al., Low-dose CT with a residual encoder-decoder convolutional neural network, IEEE Trans. Med. Imaging, 36 (2017), 2524–2535. https://doi.org/10.1109/TMI.2017.2715284 doi: 10.1109/TMI.2017.2715284
![]() |
[19] |
Z. Zhang, X. Liang, X. Dong, Y. Xie, G. Cao, A sparse-view CT reconstruction method based on combination of DenseNet and deconvolution, IEEE Trans. Med. Imaging, 37 (2018), 1407–1417. https://doi.org/10.1109/TMI.2018.2823338 doi: 10.1109/TMI.2018.2823338
![]() |
[20] |
J. C. Ye, Y. Han, E. Cha, Deep convolutional framelets: A general deep learning framework for inverse problems, SIAM J. Imaging Sci., 11 (2018), 991–1048. https://doi.org/10.1137/17M1141771 doi: 10.1137/17M1141771
![]() |
[21] |
Y. Han, J. C. Ye, Framing U-Net via deep convolutional framelets: Application to sparse-view CT, IEEE Trans. Med. Imaging, 37 (2018), 1418–1429. https://doi.org/10.1109/TMI.2018.2823768 doi: 10.1109/TMI.2018.2823768
![]() |
[22] | C. Szegedy, W. Liu, Y. Jia, P. Sermanet, S. Reed, D. Reed, et al., Going deeper with convolutions, in IEEE Conference on Computer Vision and Pattern Recognition, (2015), 1–9. https://doi.org/10.1109/CVPR.2015.7298594 |
[23] |
S. Xie, X. Zheng, Y. Chen, L. Xie, J. liu, Y. Zhang, et al., Artifact removal using improved GoogLeNet for sparse-view CT reconstruction, Sci. Rep., 8 (2018), 1–9. https://doi.org/10.1038/s41598-018-25153-w doi: 10.1038/s41598-018-25153-w
![]() |
[24] |
H. Lee, J. Lee, H. Kim, B. Cho, S. Cho, Deep-neural-network-based sinogram synthesis for sparse-view CT image reconstruction, IEEE Trans. Radiat. Plasma Med. Sci., 3 (2018), 109–119. https://doi.org/10.1109/TRPMS.2018.2867611 doi: 10.1109/TRPMS.2018.2867611
![]() |
[25] |
B. Pan, N. Qi, Q. Meng, J. Wang, S. Peng, C. Qi, et al., Ultra high speed SPECT bone imaging enabled by a deep learning enhancement method: A proof of concept, EJNMMI Phys., 9 (2022), 1–15. https://doi.org/10.1186/s40658-022-00472-0 doi: 10.1186/s40658-022-00472-0
![]() |
[26] | H. Yuan, J. Jia, Z. Zhu, Sipid: A deep learning framework for sinogram interpolation and image denoising in low-dose CT reconstruction, in IEEE 15th International Symposium on Biomedical Imaging (ISBI), (2018), 1521–1524. https://doi.org/10.1109/ISBI.2018.8363862 |
[27] | X. Dong, S. Vekhande, G. Cao, Sinogram interpolation for sparse-view micro-CT with deep learning neural network, in SPIE Medical Imaging 2019: Physics of Medical Imaging, (2019), 109482O. https://doi.org/10.1117/12.2512979 |
[28] | C. Chrysostomou, L. Koutsantonis, C. Lemesios, C. N. Papanicolas, SPECT angle interpolation based on deep learning methodologies, in IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC), (2020), 1–4. https://doi.org/10.1109/NSS/MIC42677.2020.9507966 |
[29] | I. Shiri, P. Sheikhzadeh, M. R. Ay, Deep-fill: Deep learning based sinogram domain gap filling in positron emission tomography, preprint, arXiv: 1906.07168. https://doi.org/10.48550/arXiv.1906.07168 |
[30] |
S. Li, W. Ye, F. Li, LU-Net: Combining LSTM and U-Net for sinogram synthesis in sparse-view SPECT reconstruction, Math. Biosci. Eng., 19 (2022), 4320–4340. https://doi.org/10.3934/mbe.2022200 doi: 10.3934/mbe.2022200
![]() |
[31] |
E. Xie, P. Ni, R. Zhang, X. Li. Limited-angle CT Reconstruction with generative adversarial network sinogram inpainting and unsupervised artifact removal, Appl. Sci., 12 (2022), 6268. https://doi.org/10.3390/app12126268 doi: 10.3390/app12126268
![]() |
[32] |
I. Goodfellow, J. Pouget-Abadie, M. Mirza, B. Xu, D. Warde‐Farley, S. Ozair, et al., Generative adversarial nets, Commun. ACM, 63 (2020), 139–144. https://doi.org/10.1145/3422622 doi: 10.1145/3422622
![]() |
[33] |
C. Tang, W. Zhang, L. Wang, A. Cai, N. Liang, L. Li, et al., Generative adversarial network-based sinogram super-resolution for computed tomography imaging, Phys. Med. Biol., 65 (2020), 235006. https://doi.org/10.1088/1361-6560/abc12f doi: 10.1088/1361-6560/abc12f
![]() |
[34] | Z. Li, W. Zhang, L. Wang, A. Cai, N. Liang, B. Yan, et al., A sinogram inpainting method based on generative adversarial network for limited-angle computed tomography, in 15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine, (2019), 345–349. https://doi.org/10.1117/12.2533757 |
[35] |
Y. Wang, W. Zhang, A. Cai, L. Wang, C. Tang, Z. Feng, et al., An effective sinogram inpainting for complementary limited-angle dual-energy computed tomography imaging using generative adversarial networks, J. X-Ray Sci. Technol., 29 (2021), 37–61. https://doi.org/10.3233/XST-200736 doi: 10.3233/XST-200736
![]() |
[36] | M. Mirza, S. Osindero, Conditional generative adversarial nets, preprint, arXiv: 1411.1784. https://doi.org/10.48550/arXiv.1411.1784 |
[37] |
G. E. Hinton, R. R. Salakhutdinov, Reducing the dimensionality of data with neural networks, Science, 313 (2006), 504–507. https://doi.org/10.1126/science.1127647 doi: 10.1126/science.1127647
![]() |
[38] | P. Isola, J. Y. Zhu, T. Zhou, A. A. Efros, Image-to-image translation with conditional adversarial networks, in IEEE Conference on Computer Vision and Pattern Recognition (CVPR), (2017), 1125–1134. https://doi.org/10.1109/CVPR.2017.632 |
[39] | K. Simonyan, A. Zisserman, Very deep convolutional networks for large-scale image recognition, preprint, arXiv: 1409.1556. https://doi.org/10.48550/arXiv.1409.1556 |
[40] | A. Odena, C. Odena, J. Odena, Conditional image synthesis with auxiliary classifier GANs, in The 34th International Conference on Machine Learning, 70 (2017), 2642–2651. |
[41] |
C. A. Micchelli, L. Shen, Y. Xu, Proximity algorithms for image models: denoising, Inverse Probl., 27 (2011), 045009. https://doi.org/10.1088/0266-5611/27/4/045009 doi: 10.1088/0266-5611/27/4/045009
![]() |
[42] | M. Ljungberg, S. E. Strand, M. A. King, Monte Carlo Calculations in Nuclear Medicine: Applications In Diagnostic Imaging, CRC Press, 2012. https://doi.org/10.1201/b13073 |
[43] |
M. Morphis, J. Staden, H. D. Raan, M. Ljungberg, Modelling of energy-dependent spectral resolution for SPECT Monte Carlo simulations using SIMIND, Heliyon, 7 (2021), e06097. https://doi.org/10.1016/j.heliyon.2021.e06097 doi: 10.1016/j.heliyon.2021.e06097
![]() |
[44] | S. Peltonen, U. Tuna, E. Sanchez-Monge, U. Ruotsalainen, PET sinogram denoising by block-matching and 3D filtering, in 2011 IEEE Nuclear Science Symposium Conference Record, (2011), 3125–3129. https://doi.org/10.1109/NSSMIC.2011.6152568 |
Agent i | 1 | 2 | 3 | 4 | 5 | 6 |
[22] | 260 | 1897 | 1914 | 1704 | 1635 | 1861 |
[23] | 1632 | 1899 | 1917 | 1694 | 1633 | 1869 |
Our DETC | 50 | 65 | 70 | 106 | 42 | 6 |