Design of Multi-Mode Readout Integrated Circuit for Linear Avalanche Photodiode Focal Plane
-
摘要:
线性雪崩光电二极管(Avalanche photodiode,APD)焦平面红外探测器有着广泛应用场景,APD探测器耦合具有多种模式的读出电路可在有限像元面积内实现多模式探测,提升探测系统集成度。本文设计了一种具有红外热成像模式、门控3D成像模式、激光测距模式和异步激光脉冲探测模式的APD读出电路,四种模式复用输入级电路。通过Krummenacher结构抑制背景辐射影响,扩展了光子飞行时间探测范围;提出一种改进型时刻鉴别电路,通过减小时刻鉴别误差提升距离测量精度。读出电路采用0.18 μm 3.3 V CMOS工艺设计,阵列规模128×128、像元中心距30 μm,最大电荷容量3.74 Me-。仿真结果表明,激光测距模式,在积分电容13 fF、背景电流1~150 nA条件下,背景电流响应幅值≤1.35 mV,远小于激光响应电流500 nA时280 mV的响应幅值;异步激光脉冲探测模式的幅值灵敏度约110 nA、脉宽灵敏度约4 ns;改进型时刻鉴别电路对于150~500 nA的激光脉冲响应,时刻鉴别误差约4 ns。本文设计的多模式复用APD读出电路具有一定工程应用价值。
Abstract:Linear avalanche photodiode (APD) focal-plane infrared detectors have a wide range of applications. Coupling APD detectors with multiple-mode readout integrated circuits can achieve multimode detection within a limited pixel area, thereby enhancing the integration of the detection system. In this study, we designed an APD readout integrated circuit with four modes: infrared thermal imaging, gated 3D imaging, laser ranging, and asynchronous laser pulse detection. The input-stage circuits of the four modes were multiplexed. A Krummenacher structure was used to suppress the influence of background radiation, thereby expanding the detection range of photon flight time. An improved time-discrimination circuit was proposed to enhance distance measurement accuracy by reducing the time-discrimination error. The readout integrated circuit was designed using a 0.18 μm, 3.3 V complementary metal-oxide semiconductor (CMOS), with an array size of 128×128, a pixel center distance of 30 μm, and a maximum charge capacity of 3.74 Me-. Simulation results show that in the laser ranging mode, at an integration capacitance of 13 fF and a background current of 1-150 nA, the background current response amplitude was ≤1.35 mV, which is much smaller than the response amplitude of 280 mV when the laser response current is 500 nA. The amplitude sensitivity of the asynchronous laser pulse detection mode was approximately 110 nA, and the pulse width sensitivity was approximately 4 ns. For laser pulses with a response of 150-500 nA, the time discrimination error of the improved time discrimination circuit was approximately 4 ns. The multimode multiplexed APD readout integrated circuit designed in this study has certain engineering application value.
-
0. 引言
线性APD红外探测器可以实现主/被动多模的工作方式,在被动模式下进行红外热成像;主动模式利用光子飞行时间(Time of Flight, ToF)对目标进行二维测距或3D成像,获取目标强度信息和距离信息。2011年,法国CEA-Leti公司开发了阵列规模320×256、像元间距30 μm的线性碲镉汞(HgCdTe)APD探测器,可以实现主/被动双模式成像[1],30 m范围内距离分辨率30 cm。同年,以色列SCD公司报道了一款具有4种工作模式的读出电路[2],分别为标准红外热成像模式、低噪声成像模式、异步激光脉冲探测模式和激光测距模式,阵列规模640×512,像元间距15 μm。2016年,SCD公司在上一款电路基础上进行升级[3],面阵规模扩展至1280×1024,像元中心距缩小至10 μm,最大测距距离为18 km,对应距离分辨率为9 m。2018年,德国AIM公司分别基于HgCdTe和铟镓砷(InGaAs)开发了APD探测器组件[4-5],阵列规模640×512、像元间距为15 μm,系统可用于距离选通门控成像。2022年,东南大学报道了一款多功能红外焦平面读出电路[6],面阵规模64×64,像元中心距为30 μm,每个像元可实现标准热成像模式、微光成像模式,2×2像元共享实现异步激光脉冲探测模式和激光测距模式。
目前国内多模式红外焦平面探测器读出电路尚未见到成熟产品报道,多功能读出电路并非将不同电路简单累加,因此,在有限像元面积内实现多种功能具有重要意义。本文采用电容反馈跨阻放大器(Capacitor Feedback Transimpedance Amplifier,CTIA)结构作为像元输入级,在30 μm的像元内集成了红外热成像模式、门控3D成像模式、激光测距模式和异步激光脉冲探测模式4种功能,4种模式共同复用CTIA电路,结构紧凑。激光测距模式和异步激光脉冲探测模式通过Krummenacher结构消除直流分量[7-9],仅检测脉冲信号,扩展ToF探测范围。本文提出一种适用于红外焦平面阵列、小像元读出电路的高精度时刻鉴别电路,提升ToF探测精度。
1. 读出电路结构
读出电路系统架构如图 1所示。
读出电路包括128×128、30 μm中心距的像元阵列、128列处理电路(128列并行处理电路)、4路并行输出缓冲器、片上偏置电路以及行选、列选和数字控制模块等电路,通过片外全局斜坡信号实现像素级ToF探测,4路输出中有两路输出2D信号、两路输出3D信号,每路输出缓冲器对应64列像元与列处理电路。
1.1 四种工作模式
红外热成像模式:与普通热成像模式读出电路无异,像元输入级电路对光生电荷积分,按照一定顺序通过列Buffer、输出Buffer读出每个像元中的信号。
门控3D成像模式:激光器向目标发射激光脉冲,脉冲信号在传输过程中,部分激光能量被散射或中间障碍反射,影响成像质量,通过控制探测器曝光时间,实现有选择性接收信号,从而对特定景深的区域进行成像,抑制前景或背景中障碍物或干扰源对成像的影响。门控3D成像模式下,可以同时获得目标距离信息和强度信息,如图 2所示,其中S为成像景深。
激光测距模式:激光测距模式下不受门控时间限制,可实现大动态范围测距。
异步激光脉冲检测模式:该模式可用于预警、光斑搜索和跟踪等[2, 6]。
1.2 像元电路结构
多模式读出电路的像元结构如图 3所示,主要包括CTIA输入级、Krummenacher回路、时刻鉴别电路、采样保持电路和输出缓冲器,采用CTIA结构作为信号探测电路,具有注入效率高、偏置稳定性好等优点。CTIA电路有两个反馈电容,C3D为ToF探测的感应电容,为了实现快速响应,提升探测器灵敏度,C3D容值很小,CINT为红外热成像模式的积分电容。
红外热成像模式如图 4所示,开关SW导通、SW1断开,此时大电容CINT与小电容C3D并联,RST复位结束后光生电流在积分电容上积分,INT信号控制积分时间,积分结束后,行选开关导通,信号Vout_2D通过Buffer读出像元。
门控3D成像模式如图 5所示,开关SW断开、SW1断开,此时仅有小电容C3D接入反馈环路,由于采用小电容感应的方案,激光脉冲未返回时,探测器背景电流可能将小电容提前积分至比较器翻转阈值,造成误触,严重制约了ToF探测的动态范围,通过控制探测器曝光时间可实现门控3D成像,同时获取特定门控范围内目标的强度信息和距离信息。
激光测距模式如图 6所示,开关SW断开、SW1导通,Krummenacher电路接入反馈环路,Krummenacher结构可以滤除环境光直流分量,只检测激光脉冲交流信号,从而扩展ToF探测范围。
异步激光脉冲检测模式与激光测距模式复用一套电路,如图 7所示,探测器检测随机时刻到来的激光脉冲。
门控3D成像模式和激光测距模式采用直接飞行时间探测技术(Direct Time of Flight,DToF),同时采用时间幅度转换器(Time to Amplitude Converter,TAC)标定ToF。以门控3D成像模式时序为例进行说明,如图 8所示,激光发射器向目标发射脉冲,同时参考斜坡信号开始上升,当反射脉冲到达接收阵列时,光生电流在小电容C3D上快速积分至时刻鉴别电路的阈值电压Vth,此时的参考斜坡电压信号Vout_3D被锁存,之后继续对剩余光生电流进行积分得到Vout_2D,经后续处理即可得到目标距离信息和强度信息。
ToF可由公式(1)计算:
$$ {\text{ToF}} = \frac{{{V_{{\text{out_3D}}}} - {V_{{\text{ref}}}}}}{{2 \cdot \alpha }} $$ (1) 式中:α为参考斜坡信号斜率;Vout_3D为采样电路采到的参考斜坡电压值;Vref为时间参考基准斜坡电压值起始点。
2. 读出电路影响ToF探测的关键因素
门控3D成像模式和激光测距模式由于采用了电容感应的方案,光子返回后需等待输出电压积分到比较器翻转阈值后才能触发锁存器锁存此时参考斜坡信号的电压值,导致计算得到的目标距离大于实际距离,通常比较器和锁存器的延时较为固定,如图 8所示,Tdelay通常为固定值,所以对ToF的精度影响小。
2.1 ToF探测误差
CTIA作为前端信号探测电路测距时,由于小积分电容易受背景辐射和探测器暗电流影响,造成比较器误触发,如图 9中的1st所示,复位信号RST结束后探测器开始曝光,由于背景辐射过大,激光脉冲信号还未到来时,CTIA的输出电压信号已经到达比较器翻转阈值,采样到的VO1为无效信号,这一缺点严重制约了ToF探测范围,所以通常采用门控成像方式,控制探测器曝光时间,对特定景深的区域进行探测。
除此之外,对相同距离的目标测距时,由于激光信号在传输过程中受不同环境光和目标反射率影响,回波信号大小不同,导致前端探测电路输出电压达到比较器翻转阈值电压有所差异,出现非线性的时刻鉴别误差,也称行走误差(Walk Error),时刻鉴别误差会影响ToF量化的精度,时刻鉴别误差可用图 9中的2nd、3rd和4th三种情况描述,对相同距离的物体进行探测时,脉冲信号同时返回,但回波脉冲幅度和背景干扰不同,导致采样到的参考斜坡信号不同,分别为VO2、VO3和VO4,从而引起ToF的量化出现差异。
回波脉冲到来后触发时间鉴别电路采样得到的电压值Vout_3D是ToF与行走误差、电路固定延时的和,所以实际的ToF应该为(2):
$$ {\text{ToF}} = \frac{{{V_{{\text{out_3D}}}}}}{{2 \cdot \alpha }} - {T_{{\text{walkerror}}}} - {T_{{\text{delay}}}} $$ (2) 式中:Twalkerror为行走误差;Tdelay为比较器和触发器的固定延时,又因为ToF与探测距离的关系可知实际目标距离D应该为(3):
$$ D = \frac{{\left( {\frac{{{V_{{\text{out_3D}}}}}}{{2 \cdot \alpha }} - {T_{{\text{walkerror}}}} - {T_{{\text{delay}}}}} \right) \cdot c}}{2} $$ (3) 式中:c为光速。
2.2 ToF探测精度
通常对于门控3D成像模式和激光测距模式,最关键的指标是ToF探测精度,也称时间分辨率,即探测器能够分辨的最小距离差或时间差。探测精度可推导如下,设两探测目标距离分别为D1、D2,最小距离分辨率可表示为(4):
$$ {\left( {{D_1} - {D_2}} \right)_{\min }} = \frac{{\frac{{{{\left( {{V_{3D1}} - {V_{3D2}}} \right)}_{\min }}}}{{2 \cdot \alpha }} - \Delta {T_{{\text{walkerror}}}}}}{{\frac{2}{c}}} $$ (4) 式中:ΔTwalkerror=Twalkerror1-Twalkerror2,表示两次测距的时刻鉴别误差的差值;(V3D1-V3D2)min=VLSB,即为参考斜坡信号的最小分辨率。由此可知,最小时间分辨率可表示为(5):
$$ {\left( {{\text{To}}{{\text{F}}_1} - {\text{To}}{{\text{F}}_2}} \right)_{\min }} = \frac{{{V_{{\text{LSB}}}}}}{{2 \cdot \alpha }} - \Delta {T_{{\text{walkerror}}}} $$ (5) 所以,时刻鉴别误差会直接影响探测精度,高精度ToF探测必须降低时刻鉴别误差。
2.3 改进型时刻鉴别电路
时刻鉴别电路对ToF量化精度起到关键作用,普通五管开环比较器作为时刻鉴别电路有着结构简单、功耗低和面积小等优点,适用于焦平面电路,但会引入严重的时刻鉴别误差,影响ToF探测精度。
可以利用双阈值鉴别、恒比定时鉴别、峰值辅助校正等方法减小时刻鉴别误差[10-12],不过往往会导致电路复杂,版图面积和功耗过大。针对红外焦平面探测器读出电路,本文提出了一种高精度时刻鉴别电路,如图 10所示,与传统开环五管比较器相比,在输出端增加了反馈充电支路,充电支路可以有效提升比较器从低电平到高电平的翻转速度,增加反馈限流管M7减小了不同幅值输入电压导致比较器上翻速度不一致的问题,从而降低时刻鉴别误差。其工作原理如下,设有两个不同幅值激光脉冲,前端信号探测电路输出分别为,V1和V2,且V1>V2。V1引起的比较器输出Vout1抬升速度大于V2引起的比较器输出Vout2抬升速度,比较器输出电压的抬升使得M7关断,Vout1使M7关断的速度大于Vout2,所以反馈充电支路可以减小由于输入幅度不同导致比较器上翻速度不一致的问题。
须保证通过M6的电流小于通过M5的电流,否则比较器将一直被充电至高电平而失去功能。
2.4 Krummenacher结构抑制背景
由2.1节分析可知,背景干扰严重制约了ToF探测范围。图 3展示了Krummenacher结构的具体电路,其中流经M3的电流为流经M6的一半,M5为MOS电容,保证环路的稳定性。由文献[7-9]的分析可知,图 3的Krummenacher结构有两条信号反馈路径,第一条为M1和M2,可以等效电阻Rf;第二条为M1、M5和M4,可以等效为电感Lf。激光测距模式下,Krummenacher回路接入电路,等效电路如图 11所示,此时CTIA电路的反馈回路为积分电容C3D、等效电阻Rf和等效电感Lf并联,激光脉冲未到来时,背景直流分量通过等效电感Lf泄放,激光脉冲到来时,电路的输出电压Vctia上升。此时电路仅检测交流信号,并且电路可以完成脉冲检测后输出电压的自复位,但也由于电路的自复位特性,激光测距模式下很难获得目标的强度信息。
3. 电路仿真结果
3.1 改进型时刻鉴别电路仿真
对普通五管开环比较器和改进型比较器进行对比仿真,验证改进型比较器对行走误差的抑制作用。设电流脉冲幅度从150 nA以50 nA步进增加至500 nA,ToF设置为100 ns,此时,普通五管开环比较器的时刻鉴别误差约为10.38 ns,如图 12(a)所示;改进型比较器的时刻鉴别误差仅为4 ns,如图 12(b)所示,时刻鉴别误差降低了约61.3%。所以,改进型时刻鉴别电路可以有效降低时刻鉴别误差,提升3D门控成像模式和激光测距模式的ToF探测精度。
3.2 激光测距模式仿真
为了更好观察激光测距模式对背景电流的抑制能力,激光测距模式对单个像元进行仿真。设背景电流在1~150 nA之间以10 nA步进增加,回波脉冲幅度为500 nA,脉宽8 ns,延时1 μs,模拟激光回波到来时刻。仿真结果如图 13所示,RST信号结束后读出电路处于待测模式,从图中可以看到待测脉冲没有到来之前,背景直流分量并没有将CTIA积分至比较器翻转阈值引起误触发,而是等到激光脉冲信号到来后CTIA输出电压才上升至比较器翻转阈值,所以Krummenacher结构可以有效滤除背景直流分量,扩展ToF探测范围。
图 14为激光脉冲未到来前电路输出电压随着背景直流电流变化的变化,背景电流从1 nA增加到150 nA,输出电压的变化仅有1.35 mV,相对于激光脉冲到来时280 mV的输出电压,变化率约0.48%,直流背景消除电路可以抑制不同背景对测距的干扰,有效扩展ToF探测范围,同时降低背景干扰引起的行走误差。
3.3 门控3D成像模式仿真
对像元阵列第一行的1~9列像元输入脉冲电流,脉冲电流大小在200~500 nA之间以50 nA步进增加,设电流脉冲延时从180~20 ns之间以20 ns步进递减,以此来模拟不同的光子到达时刻,结果如图 15所示,Vout_2D_1、Vout_2D_2、Vout_3D_1和Vout_3D_2,四路信号并行输出,同时获取目标强度信息与距离信息,经计算,读出电路输出的2D信号和3D信号的线性度大于99%。
3.4 红外热成像模式仿真
对像元阵列的第一行的1~11列像元输入光生电流,大小在300~500 pA之间以20 pA步进增加,仿真结果如图 16所示,其中MC为主时钟信号,Vout_2D_1和Vout_2D_2分别为读出电路的两路并行输出,经计算,读出电路线性度大于99.6%,红外热成像模式下Vout_3D_1和Vout_3D_2不输出信号。
3.5 异步激光脉冲检测模式仿真
对电路输入脉冲信号,模拟随机时刻到来的激光脉冲,仿真结果如图 17所示,RST信号为低,代表读出电路复位结束,此时电路处于待测模式,当脉冲信号到来后,读出电路输出D_OUT电平为高,代表此时探测器探测到激光信号。
幅值灵敏度和脉宽灵敏度是激光检测电路的两个重要指标。设激光脉冲宽度为8 ns,仿真电路幅值灵敏度,激光脉冲信号幅值从100 nA以10 nA步进增加至200 nA,仿真结果如图 18所示,当脉冲信号最小幅度为110 nA时,输出信号D_OUT可以发生跳变,所以在8 ns脉宽下,电路的幅值灵敏度约为110 nA;设激光脉冲幅值为200 nA,仿真电路的脉宽灵敏度,激光脉冲信号脉宽从1 ns以1 ns步进增加至8 ns,仿真结果如图 19所示,当脉冲信号最小脉宽为4 ns时,输出信号D_OUT可以发生跳变,所以在200 nA幅值下电路的脉宽灵敏度约为4 ns。
time/μs
3.6 电路版图
电路版图如图 20所示,左侧为128×128版图示意,右侧为单个像元版图,由于像元面积仅有30 μm,为了尽可能提高像元版图利用率,积分电容使用高层金属堆叠在有源区和poly上方;数字电路布局在像元最下方,尽可能与放大器、时刻鉴别电路等模拟电路隔开,降低串扰;不同层金属之间尽量垂直走线,无法垂直走线的地方尽可能不交叠,降低金属线间的寄生电容。
表 1对本文电路与已经报道的国内外先进电路进行对比,本文提出的多模式读出电路在30 μm中心距像元中集成了4种功能,具有一定实用价值。
表 1 本文电路与国内外已经报道的先进电路对比Table 1. The circuits in this paper are compared with the advanced circuits reported at home and abroadLiterature [1] Literature [3] Literature [5] Literature [6] Ours Array size 320×256 1280×1024 640×512 64×64 128×128 Pixel pitch/μm 30 10 15 30 30 Supply Voltage/V - - - 3.3 3.3 Technology CMOS 0.18 CMOS 0.18 - CMOS 0.35 CMOS 0.18 Well capacity 3.6Me- 1Me- 120ke- 1.63Me- 3.74Me- Imaging mode 1)Infrared thermal imaging
2)3D imaging1)Daylight imaging
2)Low light level imaging
3)Active imaging
4)Asynchronous laser pulse detection1)Gated viewing
2)Gated viewing(CDS)
3)ITR snapshot1)Daylight standard imaging
2)Low-light-level imaging
3)Asynchronous laser pulse detection
4)Two-dimensional laser range finder1)Infrared thermal imaging mode
2)Gated 3D
3)Laser ranging
4)Asynchronous laser pulse detectionPost-Silicon or Simulation Post-Silicon Post-Silicon Post-Silicon Post-Silicon Simulation 4. 结论
本文提出一种线性APD红外焦平面阵列多模式读出电路,具备红外热成像模式、门控3D成像模式、激光测距模式和异步激光脉冲探测模式4种功能。激光测距模式引入Krummenacher结构抑制背景对测距范围和精度的影响,实现高精度、大动态ToF探测,进一步扩展了电路功能;提出一种改进型比较器作为时刻鉴别电路,降低ToF探测的时刻鉴别误差。
-
表 1 本文电路与国内外已经报道的先进电路对比
Table 1 The circuits in this paper are compared with the advanced circuits reported at home and abroad
Literature [1] Literature [3] Literature [5] Literature [6] Ours Array size 320×256 1280×1024 640×512 64×64 128×128 Pixel pitch/μm 30 10 15 30 30 Supply Voltage/V - - - 3.3 3.3 Technology CMOS 0.18 CMOS 0.18 - CMOS 0.35 CMOS 0.18 Well capacity 3.6Me- 1Me- 120ke- 1.63Me- 3.74Me- Imaging mode 1)Infrared thermal imaging
2)3D imaging1)Daylight imaging
2)Low light level imaging
3)Active imaging
4)Asynchronous laser pulse detection1)Gated viewing
2)Gated viewing(CDS)
3)ITR snapshot1)Daylight standard imaging
2)Low-light-level imaging
3)Asynchronous laser pulse detection
4)Two-dimensional laser range finder1)Infrared thermal imaging mode
2)Gated 3D
3)Laser ranging
4)Asynchronous laser pulse detectionPost-Silicon or Simulation Post-Silicon Post-Silicon Post-Silicon Post-Silicon Simulation -
[1] Eric de Borniol, Fabrice Guellec, Johan Rothman, et al. HgCdTe-based APD focal plane array for 2D and 3D active imaging: first results on a 320×256 with 30µm pitch demonstrator[C]//Proc. of SPIE, Infrared Technology and Applications XXXVI, 2010, 7660: 76603D.
[2] Lidia Langof, Dan Nussinson, Elad Ilan, et al. Advanced multi-function infrared detector with on-chip processing[C]//Proc. SPIE, Infrared Technology and Applications XXXVII, 2011, 8012: 8012F.
[3] Fraenkel R, Berkowicz E, Bykov L, et al. High definition 10μm pitch InGaAs detector with asynchronous laser pulse detection mode[C]//Proc. SPIE, Infrared Technology and Applications XLII, 2016: 981903.
[4] Rutz F, Aidam R, Bächle A, et al. InGaAs-based SWIR photodetectors for night vision and gated viewing[C]//Proc. SPIE, Electro-Optical and Infrared Systems: Technology and Applications XV, 2018: 1079503.
[5] Sieck A, Benecke M, Eich D, et al. Short-Wave infrared HgCdTe electron avalanche photodiodes for gated viewing[J]. Electronic Materials, 2018, 47(10): 5705-5714. DOI: 10.1007/s11664-018-6425-0
[6] 郑丽霞, 刘高龙, 吴金, 等. 多功能红外焦平面读出电路设计[J]. 红外与激光工程, 2022, 51(12): 83-89. ZHENG Lixia, LIU Gaolong, WU Jin, et al. Design of multifunctional infrared FPA ROIC[J]. Infrared and Laser Engineering, 2022, 51(12): 83-89.
[7] Krummenacher F. Pixel detectors with local intelligence: An IC designer point of view[J]. Nuclear Instruments & Methods in Physics Research Section A-accelerators Spectrometers Detectors and Associated Equipment, 1991, 305: 527-532. http://www.sciencedirect.com/science/article/pii/016890029190152G
[8] Szczygiel R. Krummenacher feedback analysis for high-count-rate semiconductor pixel detector readout[C]//Proceedings of the 17th International Conference Mixed Design of Integrated Circuits and Systems - MIXDES 2010, 2010: 412-415.
[9] 叶茂, 冯书涵, 赵毅强, 等. 激光红外主被动像素读出电路设计[J]. 激光与光电子学进展, 2024, 61(15): 325-332. YE Mao, FENG Shuhan, ZHAO Yiqiang, et al. Design of pixel readout circuit for active laser and passive infrared detection[J]. Laser & Optoelectronics Progress, 2024, 61(15): 325-332.
[10] Hong Soo Cho, Chung Hwan Kim, Sang Gug Lee, A high-sensitivity and low-walk error LADAR receiver for military application[J]. IEEE Transactions on Circuits and Systems I: Regular Papers, 2014, 61(10): 3007-3015. DOI: 10.1109/TCSI.2014.2327282
[11] Lim H. Constant fraction discriminator involving automatic gain control to reduce time walk[J]. IEEE Transactions on Nuclear Science, 2014, 61(4): 2351-2356. DOI: 10.1109/TNS.2014.2339362
[12] WANG Xiayu, MA Rui, LI Dong, et al. A low walk error analog front-end circuit with intensity compensation for direct ToF LiDAR[J]. IEEE Transactions on Circuits and Systems I: Regular Papers, 2020, 67(12): 4309-4321. DOI: 10.1109/TCSI.2020.3022714