Development of a handheld fluorescence imaging system for detecting food residues
-
摘要: 为了解决食品生产加工场所的卫生检查地点不固定、检测效果受环境因素影响大等问题, 根据不同种类食品残留物在紫外光照射下产生不同波长荧光的特性, 采用荧光成像技术和分块大津算法, 开发了一种能够协助现场检测员进行视觉卫生检查的手持式荧光成像设备。利用3种常见材料表面(木制案板、不锈钢板、聚乙烯塑料板)对奶粉(体积分数为50%、33%、20%、10%、2%)和菠菜残渣(体积分数为50%、33%、25%、20%、10%)进行了实验验证。结果表明, 该系统可以协助检测员进行卫生安全检查, 不同体积分数菠菜残渣在3种材料表面检出率为100%, 不同体积分数奶粉残留物在木制案板和不锈钢板表面检出率为100%, 在聚乙烯塑料板表面, 除体积分数为2%的奶粉残留物部分检出外, 其余体积分数奶粉残留物均能有效检出。该研究为食品加工场所手持式卫生检测设备开发提供了参考。Abstract: In order to solve the problems of food production and processing sites, such as the unfixed location of health inspection and the great influence of environmental factors on the inspection effect, according to the characteristics of different kinds of food residues under ultraviolet irradiation, a handheld fluorescence imaging device was developed by using fluorescence imaging technology and block Otsu algorithm, which can assist field inspectors in carrying out visual health inspection. The equipment detected milk powder (the volume fractions are 50%, 33%, 20%, 10%, and 2%) and spinach residue (the volume fractions are 50%, 33%, 25%, 20%, and 10%) on the surface of three common materials (wooden board, stainless steel board, polyethylene board). The results show that the system can assist inspectors in carrying out health and safety inspections. The detection rate of spinach residue with different volume fractions on the surface of the three materials is 100%, the detection rate of milk powder residue with different volume fractions on the surface of the wooden board and the stainless steel board is 100%, and on the surface of polyethylene plastic plate, except for the residual part of milk powder with the volume fractions of 2%, the remaining volume fraction of milk powder residue can be effectively detected. This study provides a reference for the development of handheld sanitary inspection equipment in food processing places.
-
引言
激光雷达(light detection and ranging, LiDAR)作为一种常用光学遥感设备在风场测量应用中扮演着愈发重要的作用[1-9]。相干测风技术通过本振光与大气后向散射信号光的拍频获取大气气溶胶的多普勒频移信息,从而实现风速测量。相较直接探测方式,其具有高增益、小体积、激光器指标易实现等优点,是远程测风的通用手段。相干效率表征了接收的信号光光场与本振光光场之间的匹配程度,是相干探测模式下的关键指标。系统中激光器、探测器、光学天线等关键光学单元的性能参数,以及大气扰动程度直接影响相干效率的大小。目前,相干测风激光雷达正朝着远程、小型化、低成本的趋势发展[10-12],对关键器件的优化设计提出了更高的要求,因此, 研究测风激光雷达光学参量对相干效率的影响具有重要意义。
相干测风激光雷达的相关技术在国内外得到了广泛研究。AUGERE等人将全光纤相干多普勒测风雷达搭载于机载平台进行了挂飞测量,水平探测距离为3 km,垂直距离分辨率100 m[13]。WITSCHAS等人基于重力波研究开展了机载激光测风雷达的飞行试验,测量数据经比对,风速均方根误差小于0.2 m/s[14]。WANG等人提出了相干多普勒测风雷达联合时频算法,与传统快速傅里叶变化法相比,提升了系统时空分辨率、信噪比和最大测程[15]。MARTIN等人针对脉冲相干激光雷达系统,研制了小型化集成电路,配合信息处理算法实现了小能量远距离的硬目标探测[16]。CLEMENT等人利用频移环自相关处理法实现了1550 nm激光相干探测系统发射光束的脉冲压缩[17]。HAYLOCK等人研制了微型亚微秒电光开关,并应用于近程激光相干测速系统[18]。以上研究主要集中于对相干测风激光雷达相关试验结果进行分析,或是对风场反演和信息处理算法进行优化,但是对于光学系统的设计仿真鲜有报道,无法从理论上分析不同湍流条件下重要光学参量对系统性能指标的影响。
本文中针对相干测风激光雷达,首先概述了光学系统的组成和工作原理,然后基于光场叠加法讨论分析了关键光学参量对系统相干效率的影响,最后针对不同参量开展了整机性能验证实验,可为激光器、光学天线、接收器件的设计优化提供理论支撑。
1. 相干测风光学系统
相干测风光学系统如图 1所示。1550 nm半导体激光种子器发射几百赫兹线宽的激光,经1×2光纤分束器分束后,一路作为本振光(local oscillator, LO),频率为f0,输入光纤耦合器;另一路经声光调制器(acousto-optic modulator, AOM)调制成脉冲激光,频率为f0+f1,再经脉冲掺铒光纤放大器(erbium-doped fiber amplifier, EDFA)放大输出,然后经光学环行器和天线射向目标空域。激光与大气气溶胶作用后产生多普勒频移,后向散射信号光频率为f0+f1+f2,依次经天线和环行器进入光纤耦合器,与本振光相干拍频。最后,拍频信号由光电平衡探测器(balance detector, BD)转换为光电流,实现后续的信号处理和风场反演。定义信号光与本振光在拍频过程中振幅和相位的匹配程度为相干效率。该指标直接影响相干测风激光雷达系统的信噪比,是重要性能指标。由于光学系统涉及的参量较多,包括本振光束腰、光纤耦合器芯径、光学系统像差等,所以对相干效率的影响需要逐一分析。
2. 光学参量分析
2.1 本振光束腰
本振光通过窄线宽半导体种子激光器分束获得,光场ELO满足高斯分布:
ELO=ALOexp(−r2w0 2)exp(iφLO) (1) 式中,ALO为本振光振幅强度,r为径向尺寸,w0为束腰半径,φLO为本振光初始相位。
回波信号光经光学系统耦合进入光纤耦合器,光场Es满足艾里分布:
\begin{gathered} E_{\mathrm{s}}=A_{\mathrm{s}} \frac{{\rm{ \mathsf{ π} }}\left(\frac{d}{2}\right)^2}{\mathrm{i} \lambda_{\mathrm{s}} f} \frac{2 \mathrm{~J}_1\left(\frac{{\rm{ \mathsf{ π} }} r}{F \lambda_{\mathrm{s}}}\right)}{\frac{{\rm{ \mathsf{ π} }} r}{F \lambda_{\mathrm{s}}}} \times \\ \exp \left\{\mathrm{i}\left[k_{\mathrm{s}}\left(f+\frac{r^2}{2 f}\right)+\varphi_{\mathrm{s}}\right]\right\} \end{gathered} (2) 式中,As为信号光振幅强度,J1(·)为第1类1阶贝塞尔函数,d为光学系统入瞳直径,λs为信号光波长,f为光学系统焦距,ks为信号光波数,φs为信号光初始相位,F为光学系统F数(F=f/d)。
通过计算本振光和信号光在光纤耦合器端面的叠加场,得到相干效率η:
\begin{gathered} \eta= \\ \frac{\left|\int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{r_0} E_{\mathrm{LO} }(r, \theta) E_{\mathrm{s}}(r, \theta) \cos \left(\varphi_{\mathrm{LO} }-\varphi_{\mathrm{s}}\right) r \mathrm{~d} r \mathrm{~d} \theta\right|^2}{\int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{\infty}\left|E_{\mathrm{LO}}(r, \theta)\right|^2 r \mathrm{~d} r \mathrm{~d} \theta \cdot \int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{\infty}\left|E_{\mathrm{s}}(r, \theta)\right|^2 r \mathrm{~d} r \mathrm{~d} \theta} \end{gathered} (3) 式中,r0为纤芯半径,r、θ为极坐标参数。
系统工作波长为1550 nm,光学系统焦距为300 mm,接收光纤芯径为10 μm,假设初始状态下信号光与本振光同相位,不同光瞳口径下相干效率与本振光束腰半径的变化曲线如图 2所示。对于相同的光瞳口径,相干效率随光纤端面本振光的增大而先增大再减小,具有最大值,最大相干效率为81.4%。光瞳口径从80 mm增大到100 mm时,最优束腰半径从4.25 μm减小到3.3 μm,这是因为在系统焦距一定的情况下,光瞳口径越大,信号光在焦平面处的艾里斑半径越小,所以最佳匹配的本振光束腰半径越小。束腰半径大于最优值时,相干效率急剧下降,当w0=20 μm时,相干效率只有6.8%,原因在于在束散角相同的条件下,束腰越大,光束质量越差,本振光光场振幅分布严重变形,有效相干面积大幅减小。
2.2 光纤耦合器芯径
本振光和信号光叠加后的光场与光纤耦合器的模场匹配时,才能近似空间信号光被全部无损接收,因此光纤耦合器端面有效叠加场的相干效率η′:
\begin{gathered} \eta^{\prime}= \\ \frac{\left|\int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{r_0} E_{\mathrm{LO}}(r, \theta) E_{\mathrm{s}}(r, \theta) M(r, \theta) r \mathrm{~d} r \mathrm{~d} \theta\right|^2}{\int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{\infty}\left|E_{\mathrm{LO}}(r, \theta) E_{\mathrm{s}}(r, \theta)\right|^2 r \mathrm{~d} r \mathrm{~d} \theta \cdot \int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{\infty}|M(r, \theta)|^2 r \mathrm{~d} r \mathrm{~d} \theta} \end{gathered} (4) 式中,M(r, θ)为光纤模场分布。光纤耦合器通常采用单模光纤,模场分布为:
M(r, \theta)=\sqrt{\frac{2}{{\rm{ \mathsf{ π} }} w_{\mathrm{m}}\text{ }^2}} \exp \left(-\frac{r^2}{w_{\mathrm{m}}\text{ }^2}\right) (5) 式中,wm为模场半径。根据Marcuse经验公式,wm=r0/(lnν)0.5(2 < ν < 2.4),ν为归一化频率。
目前市场上传输光纤的芯径一般为8.4 μm~10.2 μm,也有特殊的4 μm、6 μm、12 μm和25 μm。不同光学F数情况下,相干效率与光纤耦合器芯径的变化曲线如图 3所示。对于相同的光学系统,相干效率随光纤芯径的增大而迅速增大,然后趋于平稳。这是因为当耦合器芯径小于有效叠加场尺寸时,随着耦合器芯径的增大,接收视场将会增大,有效相干面积呈芯径的平方增大;当耦合器芯径大于有效叠加场尺寸时,即使耦合器芯径继续增大,但有效相干面积已不变,所以相干效率也保持稳定。不同光学系统的最大相干效率一致,为81.4%,但达到最大相干效率需要的最小芯径r0, min有差异。当F=2时,r0, min=4.5 μm,该类光纤耦合器属于常规货架产品;当F=4时,r0, min=6 μm,该类光纤耦合器需要定制。因此,光学系统的设计在满足最大相干效率前提下,需要权衡系统尺寸重量和光纤耦合器的选型及价格等因素。
2.3 光学系统像差
相位匹配是实现相干探测系统高效拍频的关键因素,光学系统像差会造成信号光波前畸变,产生相位误差,导致耦合光斑的斯特列尔比下降。利用Zernike多项式表征信号光波前畸变,光学系统产生的5种初级像差(倾斜、离焦、球差、彗差、像散)分别表示为[20]:
\left\{\begin{array}{l} W_{\mathrm{t}}(r, \theta)=W_{111} r \cos \theta \\ W_{\mathrm{d}}(r, \theta)=W_{020} r^2 \\ W_{\mathrm{SA}}(r, \theta)=W_{040} r^4 \\ W_{\mathrm{c}}(r, \theta)=W_{131} r^3 \cos \theta \\ W_{\mathrm{a}}(r, \theta)=W_{222} r^2 \cos \theta \end{array}\right. (6) 式中,W111、W020、W040、W131、W222分别为5种像差系数。
含有光学系统像差影响的相干效率η′表示为:
\begin{gathered} \eta^{\prime}= \\ \frac{\left|\int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{r_0} E_{\mathrm{LO}}(r, \theta) E_{\mathrm{s}}(r, \theta) \cos \left[k_{\mathrm{s}} W(r, \theta)\right] r \mathrm{~d} r \mathrm{~d} \theta\right|^2}{\int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{\infty}\left|E_{\mathrm{LO}}(r, \theta)\right|^2 r \mathrm{~d} r \mathrm{~d} \theta \cdot \int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{\infty}\left|E_{\mathrm{s}}(r, \theta)\right|^2 r \mathrm{~d} r \mathrm{~d} \theta} \end{gathered} (7) 式中,信号光波数ks=2π/λs。
光纤耦合器纤芯直径为10 μm,天线口径为100 mm,本振光束腰半径为3.3 μm,相干效率与光学系统像差的变化曲线如图 4所示。倾斜对相干效率的影响最大,当波像差为0.2λ(λ为工作波长)时,相干效率为71.8%,下降了9.6%,原因在于波前倾斜造成焦平面处信号光光斑偏移,与本振光同轴度降低,两光场重叠面积快速减小,损耗了有用的光能量。同等波像差,像散和离焦导致相干效率下降4.4%和2.1%,而彗差和球差造成的影响近似相等且最小,约0.3%。实际应用中,环境温度变化和机械振动必然存在,所以应采用主动或被动消热差设计,使整个光学系统波像差均方根值(root mean square,RMS)不大于0.06λ,同时对光学系统进行减震设计。
3. 大气湍流的影响
相干测风激光雷达以大气为传输通道,由于大气湍流的存在,回波信号光会引入波前相位差[19]。根据冯卡曼模型,平均相位延迟〈φ2〉表示为[20]:
\left\langle\varphi^2\right\rangle=0.782 L k_{\mathrm{s}}{ }^2 C_n{ }^2 L_0{ }^{5 / 3} (8) 式中,L为探测距离,Cn2为大气折射率结构常数,L05/3为湍流外尺度。
考虑大气湍流影响后,相干效率η′表示为:
\begin{gathered} \eta^{\prime}= \\ \frac{\left|\int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{r_0} E_{\mathrm{LO}}(r, \theta) E_{\mathrm{s}}(r, \theta) \cos \varphi r \mathrm{~d} r \mathrm{~d} \theta\right|^2}{\int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{\infty}\left|E_{\mathrm{LO}}(r, \theta)\right|^2 r \mathrm{~d} r \mathrm{~d} \theta \cdot \int_0^{2 {\rm{ \mathsf{ π} }}} \int_0^{\infty}\left|E_{\mathrm{s}}(r, \theta)\right|^2 r \mathrm{~d} r \mathrm{~d} \theta} \end{gathered} (9) 通常认为Cn2 < 10-17 m-2/3为弱湍流,10-17 m-2/3≤Cn2≤10-15 m-2/3为中湍流,Cn2>10-15 m-2/3为强湍流。不同湍流强度情况下,相干效率与探测距离的变化曲线如图 5所示。在相同天气条件下,相干效率随探测距离的增大而减小。对于强湍流情形,当L=8 km时,η′降至28.2%,再往远距离探测较困难;对于弱湍流情形,当L=30 km时,η′仍有41.6%,与目前报道的可实现探测最远距离相吻合。在同一探测距离处,随着湍流强度增大,相干效率减小,但减小量在不同距离区间有差异。当L < 3 km时,|Δη′| < 12.2%;当3 km≤L≤30 km时,η′急剧下降,特别地,在L=17 km处,减小量为最大值,|Δη′|max=48.1%。这是因为近距离探测时,大气散射回波信号强度较强,湍流影响较小;中远距离探测时,激光受大气扰动及退偏影响严重,散斑效应显著。
4. 实验与讨论
分束器本振光支路分别采用不同芯径的单模光纤跳线模拟本振光束腰w0,利用不同光瞳口径d的光学天线在同一测试地点,对相同环境的径向风速进行测量,如图 6所示。由于测风激光雷达相干效率表征信号光场与本振光场的匹配度,所以最终由系统获得的信号频谱强度I或信噪比α量化表示。
实验中,激光器工作波长为1550 nm,输出单脉冲能量为300 μJ,AOM中频为80 MHz,距离门长度设置为80 m。测量并记录不同天线口径下,不同本振光束腰拍频计算得到的第40个距离门的频谱强度,如表 1所示。当天线口径为80 mm时,随着本振光的增大,频谱强度先增大后减小,因此本振光束腰半径存在最优解;随着天线口径增大,最优解会减小,与仿真结果相符。
表 1 不同天线口径、不同本振光束腰下的频谱强度Table 1. Frequency spectrum intensity of different antenna aperture and LO beam waistd/mm 80 100 w0/μm 4 4.2 5 4 4.2 5 I/a.u. 1.7×1011 2×1011 1.8×1011 1.5×1011 1.3×1011 1.2×1011 为了验证接收光纤芯径对相干效率的影响,订制了纤芯直径D为12 μm和25 μm的光纤耦合器,与常用的9 μm芯径光纤耦合器进行比对实验。光学系统焦距为300 mm,口径80 mm和100 mm天线对应的F数分别为3.75和3, 测量结果如表 2所示。当光学F数不变时,频谱强度随光纤芯径的增大而增大,然后趋于平稳。在平稳状态前,对于相同的光纤芯径,F数越小,频谱强度越大。例如,对于9 μm光纤芯径,F数从3增加到3.75时,频谱强度下降25%。这是因为F数与光学口径成反比,大口径光学天线接收的回波信号更多。
表 2 不同光学F数、不同接收光纤芯径下的频谱强度Table 2. Frequency spectrum intensity of different F number and receive core diameterF 3 3.75 D/μm 9 12 25 9 12 25 I/a.u. 1.2×1011 1.4×1011 1.4×1011 0.9×1011 1.4×1011 1.4×1011 为了验证光学系统像差对相干效率的影响,加工了两套相同的光学系统,一套经过严格装调测试,满足波像差RMS为0.036λ,如图 7所示。另一套未做装调,并处于离焦状态,RMS远大于0.06λ。将两套光学系统装入测风激光雷达系统进行比对测量,记录在第40个距离门的频谱强度,前者为3×1011a.u.,信号细锐清晰,后者为2.5×107a.u.,尖峰噪声激增,信号完全淹没在噪声中。实验结果表明, 满足光学像差要求的系统在第40个距离门处的中频信号频谱可有效提取,相干效率高;反之,信噪比极低,相干效率很低,与仿真结果一致。
开展大气湍流影响实验时,将雷达安置于同一测试地点,针对典型天气条件(条件一:大雨或雾霾,能见度V≤5 km;条件二:晴朗,大气气溶胶含量高,能见度V≥20 km)进行测量,并利用系统在不同探测距离处的信噪比α表征相干效率,测量结果如图 8所示。光学天线口径为100 mm,束散角为40 μrad,波像差为0.036λ,本振光束腰为3.3 μm,光纤耦合器芯径为12 μm。晴空条件下,系统最远测程可达15 km,大雨或雾霾天气条件下,系统最远测程6 km,原因在于雨水和大气微粒对信号光的吸收和散射导致光强起伏和相位波动,造成相干效率下降。对于近程,即L < 3 km时,湍流对系统信噪比影响较小;当L>3 km时,在强湍流情况下,系统信噪比随探测距离的增加急剧下降;当L=6 km时,系统信噪比小于5,提取中频信号频谱变得困难。实验结果与仿真结果基本吻合, 实测的最远测程比理论值小,且信噪比曲线存在局部波动,这是因为建模仿真使用的冯卡曼模型将相位噪声理想化为平均相位延迟,忽略了实际情况中大气湍流的随机性。
5. 结论
针对相干测风激光雷达,阐述了光学系统组成和工作原理,定义了系统相干效率; 基于本振光和信号光的光场叠加理论,讨论分析了本振光束腰、光纤耦合器芯径、光学系统像差对相干效率的影响; 引入大气湍流平均相位延迟参量,建模分析了大气扰动对相干效率的影响,并开展了验证实验,测量了不同光学参量的系统在典型天气条件下的信噪比。结果表明,光瞳口径从80 mm增大到100 mm时,最优束腰半径从4.25 μm减小到3.3 μm。当光学F数不变时,相干效率随接收光纤芯径的增大而增大,最后趋于平稳。对于9 μm光纤芯径,F数从3增大到3.75时,频谱强度下降25%。为保证系统相干效率,光学系统波像差RMS不大于0.06λ。在相同天气条件下,相干效率随探测距离的增大而减小,晴朗天气条件下,能见度不小于20 km时,系统最远测程可达15 km。对于近程,湍流对相干效率的影响较小,当L>3 km时,在强湍流情况下,系统信噪比随探测距离的增加急剧下降;当L=6 km时,系统信噪比小于5。因为大气干扰和退偏效应客观存在,所以在优化光学参量的同时,应考虑信号光与本振光的偏振匹配。
本文中对相干测风激光雷达激光器和光学天线的设计,以及接收器件的选型和参数优化具有参考意义。
-
表 1 相机主要配置参数
Table 1 Main camera configuration parameters
register name address value TIMING DVPHO H 0x3808 0x03 TIMING DVPHO L 0x3809 0x20 TIMING DVPVO L 0x380a 0x01 TIMING DVPVO L 0x380b 0xe0 AEC PK MANUAL 0x3503 0x03 AEC PK EXPOSURE 0x3500 0x00 AEC PK EXPOSURE 0x3501 0x6e AEC PK EXPOSURE 0x3502 0x30 FORMAT CONTROL 0x4300 0x61 ISP FORMAT 0x501f 0x01 -
[1] 李红秋, 郭云昌, 宋壮志, 等. 2019年中国大陆食源性疾病暴发监测资料分析[J]. 中国食品卫生杂志, 2021, 33(6): 650-656. https://www.cnki.com.cn/Article/CJFDTOTAL-ZSPZ202106003.htm LI H Q, GUO Y Ch, SONG Zh Zh, et al. Analysis of foodborne di-sease outbreaks in China in 2019[J]. Chinese Journal of Food Hygiene, 2021, 33(6): 650-656(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZSPZ202106003.htm
[2] 李晓凤, 秦宇龙, 刘姜汝, 等. 即食果蔬中食源性致病菌风险评估的研究进展[J]. 食品与发酵工业, 2022, 48(18): 322-328. https://www.cnki.com.cn/Article/CJFDTOTAL-SPFX202218046.htm LI X F, QIN Y L, LIU J R, et al. Research progress in risk assessment of foodborne pathogens in ready-to-eat fruits and vegetables[J]. Food and Fermentation Industries, 2022, 48(18): 322-328 (in Ch-inese). https://www.cnki.com.cn/Article/CJFDTOTAL-SPFX202218046.htm
[3] 黎远鹏, 黄富荣, 董佳, 等. 荧光光谱成像技术结合主成分分析与Fisher判别快速鉴别肉苁蓉[J]. 光谱学与光谱分析, 2015, 35(3): 689-694. https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201503028.htm LI Y P, HUANG F R, DONG J, et al. Rapid identification of cistanche via fluorescence spectrum imaging technology combined with principal components analysis and Fisher distinction[J]. Spectroscopy and Spectral Analysis, 2015, 35(3): 689-694(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201503028.htm
[4] 宋鑫澍, 陈国庆, 朱焯炜, 等. 3维荧光光谱测定白酒年份酒中乙酸的体积分数[J]. 激光技术, 2018, 42(4): 531-535. DOI: 10.7510/jgjs.issn.1001-3806.2018.04.018 SONG X S, CHEN G Q, ZHU Zh W, et al. Determination of volume fraction of acetic acid in Chinese aged liquor by 3-D fluorescence spectrometry[J]. Laser Technology, 2018, 42(4): 531-535(in Ch-inese). DOI: 10.7510/jgjs.issn.1001-3806.2018.04.018
[5] SHIBATA M, CHEN J, OKADA K, et al. Detection of food residues on stainless steel surfaces using fluorescence fingerprint[J]. Food Science and Technology Research, 2020, 26(3): 389-397. DOI: 10.3136/fstr.26.389
[6] ABBAN S, JAKOBSEN M, JESPERSEN L. Assessment of interplay between UV wavelengths, material surfaces and food residues in open surface hygiene validation[J]. Journal of Food Science and Technology, 2014, 51(12): 3977-3983. DOI: 10.1007/s13197-013-0927-9
[7] 康震, 丁利, 程云辉, 等. 基于上转换发光技术的便携式食品安全检测仪设计[J]. 传感器与微系统, 2019, 38(11): 114-116. https://www.cnki.com.cn/Article/CJFDTOTAL-CGQJ201911032.htm KANG Zh, DING L, CHENG Y H, et al. Design of portable food safety detector based on up-conversion phosphor technology[J]. Transducer and Microsystem Technologies, 2019, 38(11): 114-116 (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-CGQJ201911032.htm
[8] 张利永. 基于智能手机黄曲霉菌毒素的快速检测[D]. 西安: 电子科技大学, 2019: 36-45. ZHANG L Y. Rapid detection of aflatoxin based on smartphone[D]. Xi'an: University of Electronic Science and Technology of China, 2019: 36-45(in Chinese).
[9] SEO Y, LEE H, MO C, et al. Multispectral fluorescence imaging technique for on-line inspection of fecal residues on poultry carcasses[J]. Sensors, 2019, 19(16): 1-16.
[10] BURFOOT D, TINKER D, THORN R, et al. Use of fluorescence imaging as a hygiene indicator for beef and lamb carcasses in UK slaughterhouses[J]. Biosystems Engineering, 2011, 109(3): 175-185.
[11] GORJI H T, SHAHABI S M, SHARMA A, et al. Combining deep learning and fluorescence imaging to automatically identify fecal contamination on meat carcasses[J]. Scientific Reports, 2022, 12(1): 1-11.
[12] SUEKER M, STROMSODT K, GORJI H T, et al. Handheld multispectral fluorescence imaging system to detect and disinfect surface contamination[J]. Sensors, 2021, 21(21): 1-15.
[13] BECK E A, LEFCOURT A M, LO Y M, et al. Use of a portable fluorescence imaging device to facilitate cleaning of deli slicers[J]. Food Control, 2015, 51(5): 256-262.
[14] WIEDERODER M S, LIU N T, LEFCOURT A M, et al. Use of a portable hyperspectral imaging system for monitoring the efficacy of sanitation procedures in produce processing plants[J]. Journal of Food Engineering, 2013, 117(2): 217-226.
[15] LEE H, KIM M S, CHAO K, et al. Development of fluorescence based handheld imaging devices for food safety inspection[C]// International Society for Optical Engineering. Sensing for Agriculture and Food Quality and Safety Ⅴ. Bellingham, USA: International Society for Optical Engineering, 2013: 45-52.
[16] OH M, LEE H, CHO H, et al. Detection of fecal contamination on beef meat surfaces using handheld fluorescence imaging device (HFID)[C]// International Society for Optical Engineering. Sensing for Agriculture and Food Quality and Safety Ⅷ. Bellingham, USA: International Society for Optical Engineering, 2016: 161-166.
[17] EVERARD C D, KIM M S, LEE H. Assessment of a handheld fluorescence imaging device as an aid for detection of food residues on processing surfaces[J]. Food Control, 2016, 59(1): 243-249.
[18] 孙颖馨. 基于FPGA红外成像光谱数据处理系统研究[J]. 激光技术, 2019, 43(6): 763-767. DOI: 10.7510/jgjs.issn.1001-3806.2019.06.006 SUN Y X. Research of data processing systems for infrared imaging spectrometer based on FPGA[J]. Laser Technology, 2019, 43(6): 763-767(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2019.06.006
[19] 周小萌, 吴静静, 安伟. 一种鲁棒的非均匀光场中SD卡形态识别算法[J]. 传感技术学报, 2019, 32(4): 549-554. https://www.cnki.com.cn/Article/CJFDTOTAL-CGJS201904013.htm ZHOU X M, WU J J, AN W. A robust SD card shape recognition algorithm in inhomogeneous optical field[J]. Chinese Journal of Sensors and Actuators, 2019, 32(4): 549-554(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-CGJS201904013.htm
[20] BARROS W K P, DIAS L A, FERNANDES M A C. Fully parallel implementation of otsu automatic image thresholding algorithm on FPGA[J]. Sensors, 2021, 21(12): 1-17.
[21] EVERARD C D, KIM M S, LEE H. A comparison of hyperspectral reflectance and fluorescence imaging techniques for detection of contaminants on spinach leaves[J]. Journal of Food Engineering, 2014, 143(24): 139-145.