CN110514305B - 火焰温度场测量系统光场相机数量和机位布置的优化方法 - Google Patents

火焰温度场测量系统光场相机数量和机位布置的优化方法 Download PDF

Info

Publication number
CN110514305B
CN110514305B CN201910776336.XA CN201910776336A CN110514305B CN 110514305 B CN110514305 B CN 110514305B CN 201910776336 A CN201910776336 A CN 201910776336A CN 110514305 B CN110514305 B CN 110514305B
Authority
CN
China
Prior art keywords
flame
light
cameras
arrangement
micro
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910776336.XA
Other languages
English (en)
Other versions
CN110514305A (zh
Inventor
张彪
王志博
魏泽坤
许传龙
王晨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201910776336.XA priority Critical patent/CN110514305B/zh
Publication of CN110514305A publication Critical patent/CN110514305A/zh
Application granted granted Critical
Publication of CN110514305B publication Critical patent/CN110514305B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/0014Radiation pyrometry, e.g. infrared or optical thermometry for sensing the radiation from gases, flames
    • G01J5/0018Flames, plasma or welding
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J5/02Constructional details
    • G01J5/03Arrangements for indicating or recording specially adapted for radiation pyrometers
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J2005/0077Imaging

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Strategic Management (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Biophysics (AREA)
  • Artificial Intelligence (AREA)
  • Development Economics (AREA)
  • Computational Linguistics (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • Plasma & Fusion (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)
  • Radiation Pyrometers (AREA)

Abstract

本发明公开了一种火焰温度场测量系统光场相机数量和机位布置的优化方法,包括如下步骤:对实际成像相机进行标定,获得成像相机参数;设定测量区域形状尺寸及火焰的衰减系数分布,根据分辨率要求进行网格划分;进行采样光线逆向追迹,获得火焰温度场测量的重建矩阵A;计算权重矩阵的Sen及其标准差σ,将其作为采样特性寻优算法的目标值;利用量子微粒群算法对测量系统的相机数量和机位布置角度进行优化,获得最佳的采样布置。本发明研究多相机光场成像系统对火焰辐射信息的获取能力,并优化成像系统的参数,提高光场采样质量,给多相机的布置方案提供依据。

Description

火焰温度场测量系统光场相机数量和机位布置的优化方法
技术领域
本发明涉及一种多光场相机火焰温度场测量系统相机数量和机位布置的优化方法,属于火焰温度测量技术领域。
背景技术
作为燃烧的重要参数之一,火焰温度影响火焰内部光与热的产生及传递,同时也影响着燃烧反应的进程。因此,火焰温度既是燃烧机理的研究基础,与火焰燃烧状态密切相关,也是火焰浓度、燃尽度和燃烧放热量分布的诊断依据,还是燃烧装置(如锅炉、发动机、流化床等)性能分析、故障诊断与优化设计的重要依据,是实现能源智慧管理不可或缺的数据。
目前,火焰温度测量主要分为接触式测量和非接触式测量两种。接触式测量方式简单、易于实现在工程上的应用,但是无法完全复现火焰温度的全貌,而且往往局限于单点测量,响应时间长。非接触式测量可以实现无扰、快速的火焰测量。其中,非接触式辐射成像技术利用火焰辐射图像处理温度测量,具有精度高、非接触、实时测量等优点,且能够获得火焰三维温度分布信息。
对光场相机而言,远离相机处的物点释放热辐射,辐射光线穿透靠近相机处的介质而被相机接收。由于成像探测面上每个像素的视角范围很小,相应的落在该范围内的火焰微元体数量很少,空间内许多微元体对单个像素辐射能的传输贡献为零,各微元体对某一像素的能量贡献份额差异很大,造成重建三维温度场的数学问题具有较强奇异性。另一方面,火焰介质的散射与吸收特性使得热辐射能沿着传播路径被逐步衰减,对于光学厚度较大的火焰,远离相机处物点的辐射信息甚至无法达到相机,远离光场相机的火焰微元体所发射的辐射能对光场成像的影响急剧减小,造成火焰辐射信息的损失,损害三维温度场的重建精度。因此,提高多光场相机系统的火焰辐射采样均匀性,对于降低火焰三维温度场可视化问题的病态特性具有重要意义。
发明内容
本发明所要解决的技术问题是针对上述现有技术存在的不足,而提供一种提高光场采样质量的火焰温度场测量系统光场相机数量和机位布置的优化方法。
为解决上述技术问题,本发明采用的技术方案是:
火焰温度场测量系统光场相机数量和机位布置的优化方法,其特征在于,包括如下步骤:
步骤一、对实际成像相机进行标定,获得成像相机参数;
步骤二、设定测量区域形状尺寸及火焰的衰减系数分布,根据分辨率要求进行网格划分;
步骤三、根据几何光学理论和光场成像模型进行采样光线逆向追迹,获得火焰温度场测量的重建矩阵A;
步骤四、计算权重矩阵的Sen及其标准差σ,将其作为采样特性寻优算法的目标值;
步骤五、利用量子微粒群算法对测量系统的相机数量和机位布置角度进行优化,获得最佳的采样布置。
步骤三中,假设火焰满足各向同性散射,可定义广义源项S(r)如下:
Figure BDA0002174754030000021
式中,ω是散射反射率,ω=σs/β。将火焰划分成V个微元体并编号(1,2,v,...,V),光场相机探测面共有M个像素并编号(1,2,m,...,M),将rm处像素的辐射光线穿过的火焰微元体个数记为n,此时该路径上各微元体的编号依次记为v1,v2,vi,(vj),...,vn,火焰热辐射传输方程的离散方程如下:
Figure BDA0002174754030000022
式中I(r,Ω)是火焰在r位置处Ω方向上的辐射强度,W/(m2·sr)],Ib(r)是火焰在r位置处的自身辐射强度,W/(m2·sr),可认为是黑体辐射强度,β(r)是衰减系数,m-1,κ(r)是吸收系数,m-1,σs(r)是散射系数,m-1,Φ(Ω′,Ω)是由Ω′方向入射向Ω方向散射的散射相函数。
根据离散方程得到某成像光线的出射辐射强度值如下:
Figure BDA0002174754030000023
式中Δs是该成像光线穿过当前火焰微元体所经过的路径长度。可以看出,成像光线的出射辐射强度与每个火焰微元体的广义源项线性相关,联立M个像素成像光线的火焰辐射强度计算式,得到线性方程组式
I=A·S (4)
式中,I是光场相机探测面上各像素接收到的辐射强度组成的矩阵[I1,I2,...,Im,...,IM],A是系数矩阵,[A1 1,A1 2,...,A1 V;A2 1,A2 2,...,A2 V;...;AM 1,AM 2,...,AM V],S是火焰微元体的广义源项组成的矩阵[S1,S2,...,Sv,...SV]。
步骤四中,统计每个微元体内全部采样光线的耗散系数Am,v之和,此和值表征了多相机系统对第v个火焰微元体的探测灵敏程度所占权重,定义和值为权重系数Sen(v),其表达式如式(6),所有微元体的权重系数组成权重矩阵Sen。
Figure BDA0002174754030000031
因此,对系统采样特性的均匀化可等价于对探测权重场的均匀化,用权重矩阵的系数的标准差进行表征,用σ表示,是评价光线性能特性的一项重要指标。
步骤五中,量子微粒群算法对多光场相机的光场采样特性进行均匀化分析,需要设置系统控制参数包括粒子总数、问题维数、最大迭代步数、目标函数最大容忍度。量子微粒群算法计算流程图如图1所示。
与现有技术相比,本发明的有益效果是:
本发明针对弥散介质的半透明特性,提出探测权重矩阵均匀度这一性能指标,来评价多光场相机火焰温度场测量系统相机数量和机位布置的优化方法,据此分析了火焰微元体划分数量、火焰的物理性质和多相机系统配置参数对光场采样特性的影响,并通过量子微粒群算法实现多光场相机系统的配置方案寻优,从而实现火焰辐射采样的均匀性。本发明可为研究多光场相机火焰温度场测量系统相机数量和机位布置的优化作了配置方案寻优,为多光场相机测温系统的研制提供设计依据。
附图说明
图1为改进的量子微粒群算法计算流程图。
图2为火焰辐射光场成像的采样模型图。
图3为2种不同的衰减系数分布图,其中a为第一种设置,b为第二种设置。
图4不同采样角度下火焰辐射采样的权重矩阵标准差。
图5两个机位位置寻优后相机布置以及正交布置图。
图6相同衰减系数分布下,火焰辐射采样权重分布及光场相机布置,其中a为机位位置寻优后最佳角度的权重矩阵标准差,b为两相机正交下的权重矩阵标准差。
具体实施方式
下面结合附图和具体施例,进一步阐述本发明。应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
火焰温度场测量系统光场相机数量和机位布置的优化方法,包括如下步骤:
步骤1、设置光场相机参数;
主透镜和微透镜的焦距分别为50mm和567μm,主透镜-微透镜距离为57mm,微透镜-CCD距离为794μm,物距为548.3mm,微透镜阵列中的微透镜个数为60个,微透镜直径为54.2μm,取CCD探测面的像素个数为500个,各像素的尺寸为5.5μm×5.5μm。
步骤二、设置火焰参数;
为了在在贴近燃烧诊断实际情况的同时考虑到研究效率,采用火焰的二维横截面作为模型进行研究,此时微透镜阵列与CCD探测面均为线阵。首先设置火焰横截面为边长22mm的正方形,以火焰中心为原点,探测面法向量方向为x轴方向,沿线阵向上为y轴方向建立全局坐标系,如图2所示。将火焰沿长度与宽度方向划分为18×18个微元体,设置2种不同的衰减系数β的分布如图3所示,其中两种设置的背景区域的衰减系数均为大小为6m-1,第一种设置中:矩形、圆形、三角形区域的衰减系数依次为9m-1、12m-1、15m-1,如图3(a)所示,第二种设置中:矩形、圆形、三角形区域的衰减系数依次为9m-1、15m-1、12m-1,如图3(b)所示。
步骤三、逆向追迹成像光线;
假设火焰满足各向同性散射,可定义广义源项S(r)如下:
Figure BDA0002174754030000041
式中,ω是散射反射率,ω=σs/β。将火焰划分成V个微元体并编号(1,2,v,...,V),光场相机探测面共有M个像素并编号(1,2,m,...,M),将rm处像素的辐射光线穿过的火焰微元体个数记为n,此时该路径上各微元体的编号依次记为v1,v2,vi,(vj),...,vn,火焰热辐射传输方程的离散方程如下:
Figure BDA0002174754030000042
式中I(r,Ω)是火焰在r位置处Ω方向上的辐射强度,W/(m2·sr)],Ib(r)是火焰在r位置处的自身辐射强度,W/(m2·sr),可认为是黑体辐射强度,β(r)是衰减系数,m-1,κ(r)是吸收系数,m-1,σs(r)是散射系数,m-1,Φ(Ω′,Ω)是由Ω′方向入射向Ω方向散射的散射相函数。
根据离散方程得到某成像光线的出射辐射强度值如下:
Figure BDA0002174754030000043
式中,Δs是该成像光线穿过当前火焰微元体所经过的路径长度。可以看出,成像光线的出射辐射强度与每个火焰微元体的广义源项线性相关,联立M个像素成像光线的火焰辐射强度计算式,得到线性方程组式
I=A·S (4)
式中,I是光场相机探测面上各像素接收到的辐射强度组成的矩阵[I1,I2,...,Im,...,IM],A是系数矩阵,[A1 1,A1 2,...,A1 V;A2 1,A2 2,...,A2 V;...;AM 1,AM 2,...,AM V],S是火焰微元体的广义源项组成的矩阵[S1,S2,...,Sv,...SV]。
步骤四、计算权重矩阵Sen和标准差σ,作为采样特性优化的目标值;
步骤五、利用量子微粒群算法对测量系统的相机数量和机位布置角度进行优化,获得最佳的采样布置。需要设置的系统控制参数包括粒子总数、问题维数、最大迭代步数、目标函数最大容忍度,量子微粒群算法计算流程图如图1所示。
在本实施例中,设置量子微粒群算法的初始参数具体为:粒子总数为20、问题维数为5、最大迭代步数为200、目标函数最大容忍度为1。用量子微粒群算法对图3(a)处理,得到不同采样角度下火焰辐射采样的权重矩阵标准差如图4所示。图4中一个相机布置角度为107.9°是为最优解。相同衰减系数分布下,对图(3)b进行两个机位位置的寻优,得到两个机位最佳角度为40.1°与98.7°,如图5所示。图6为相同衰减系数分布下,火焰辐射采样权重分布及光场相机布置,其中a为机位位置寻优后最佳角度的权重矩阵标准差,b为两相机正交下(图5中的0°和90°)的权重矩阵标准差。

Claims (3)

1.一种火焰温度场测量系统光场相机数量和机位布置的优化方法,其特征在于,包括如下步骤:
步骤一、对实际成像相机进行标定,获得成像相机参数;
步骤二、设定测量区域形状尺寸及火焰的衰减系数分布,根据分辨率要求进行网格划分;
步骤三、进行采样光线逆向追迹,获得火焰温度场测量的重建系数矩阵A;
步骤四、计算权重矩阵Sen及其标准差σ,将其作为采样特性寻优算法的目标值;
步骤五、利用量子微粒群算法对测量系统的相机数量和机位布置角度进行优化,获得最佳的采样布置;
步骤三中重建系数矩阵A的各项元素值为:
Figure FDA0002730749080000011
式中,m是有效光线的序号,v代表微元体的编号,j代表光线沿传输方向穿过的第j个火焰微元体,N代表每根光线穿过的火焰微元体的个数,βvj是光线沿传输方向穿过的第j个火焰微元体所对应的第v个微元体的衰减系数,Δsvj是成像光线穿过沿传输方向第j个火焰微元体所对应的第v个微元体所经过的路径长度;
步骤四中的权重矩阵Sen定义为:
Figure FDA0002730749080000012
式中,M表示光线总数。
2.根据权利要求1所述的优化方法,其特征在于,所述相机参数包括主透镜与微透镜焦距、主透镜-微透镜距离、微透镜-CCD距离、物距、微透镜个数与直径。
3.根据权利要求1所述的优化方法,其特征在于,步骤五中相机的数量需要遍历有限多个,根据量子微粒群算法的对每个数量下相机的布置角度进行优化,然后对比获得最优相机个数及其对应的机位布置。
CN201910776336.XA 2019-08-21 2019-08-21 火焰温度场测量系统光场相机数量和机位布置的优化方法 Active CN110514305B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910776336.XA CN110514305B (zh) 2019-08-21 2019-08-21 火焰温度场测量系统光场相机数量和机位布置的优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910776336.XA CN110514305B (zh) 2019-08-21 2019-08-21 火焰温度场测量系统光场相机数量和机位布置的优化方法

Publications (2)

Publication Number Publication Date
CN110514305A CN110514305A (zh) 2019-11-29
CN110514305B true CN110514305B (zh) 2020-12-11

Family

ID=68626168

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910776336.XA Active CN110514305B (zh) 2019-08-21 2019-08-21 火焰温度场测量系统光场相机数量和机位布置的优化方法

Country Status (1)

Country Link
CN (1) CN110514305B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112284546B (zh) * 2020-10-16 2022-11-01 中国航发四川燃气涡轮研究院 基于双目视觉的尾喷口温度场可视化装置及其识别方法
CN112418079B (zh) * 2020-11-20 2021-08-31 浙江科技学院 基于安全等级的地图分析与寻优路径算法结合的工业建筑视频监控安装方法
CN114723621B (zh) * 2022-03-07 2022-09-27 北京拙河科技有限公司 基于感知模型的光场成像方法及系统

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2390675A (en) * 2002-07-10 2004-01-14 Univ Greenwich Flame characteristic monitor using digitising image camera
CN101493980B (zh) * 2009-03-05 2011-07-20 中国科学技术大学 一种基于多特征融合的快速视频火焰探测方法
CN104240429B (zh) * 2014-08-27 2016-09-07 中国资源卫星应用中心 一种红外与多光谱波段组合监测森林火灾的方法
CN105606222B (zh) * 2015-09-06 2018-07-24 东南大学 一种火焰三维温度场的测量装置及测量方法
CN106600687B (zh) * 2016-12-08 2019-12-27 南京理工大学 一种多方向火焰发射层析系统
KR20180118956A (ko) * 2017-04-24 2018-11-01 티텍 주식회사 온도, 열유속 및 부식 측정 시스템
CN108225577B (zh) * 2017-12-19 2020-02-18 东南大学 基于双光谱辐射信息的火焰三维温度分布重建方法
CN109974858A (zh) * 2019-01-10 2019-07-05 山东理工大学 一种基于热化学发光层析成像技术的燃烧火焰三维轮廓快速重建方法
CN109916531B (zh) * 2019-03-04 2020-09-11 东南大学 一种基于光场重聚焦的半透明火焰三维温度场测量方法
CN109932420A (zh) * 2019-03-08 2019-06-25 李言钦 小波基函数用于炉膛断面流场声波法测量重建的方法

Also Published As

Publication number Publication date
CN110514305A (zh) 2019-11-29

Similar Documents

Publication Publication Date Title
CN110514305B (zh) 火焰温度场测量系统光场相机数量和机位布置的优化方法
Zhao et al. Optical sectioning tomographic reconstruction of three-dimensional flame temperature distribution using single light field camera
Sun et al. Three-dimensional temperature field measurement of flame using a single light field camera
Sun et al. A novel calibration method of focused light field camera for 3-D reconstruction of flame temperature
WO2018210072A1 (zh) 一种基于多目图像识别的定日镜面形检测系统及方法
Xiao et al. A model-based approach for optical performance assessment and optimization of a solar dish
CN105606222A (zh) 一种火焰三维温度场测量的成像装置、测量装置及测量方法
Ma et al. Combined soot optical characterization using 2-D multi-angle light scattering and spectrally resolved line-of-sight attenuation and its implication on soot color-ratio pyrometry
CN108225577B (zh) 基于双光谱辐射信息的火焰三维温度分布重建方法
CN109916531B (zh) 一种基于光场重聚焦的半透明火焰三维温度场测量方法
Abuseada et al. Characterization of a new 10 kWe high flux solar simulator via indirect radiation mapping technique
Yuan et al. Light-field-camera imaging simulation of participatory media using Monte Carlo method
Arqueros et al. A novel procedure for the optical characterization of solar concentrators
CN111795746B (zh) 基于主被动光学层析融合探测的火焰多参数场协同测量方法
Niu et al. Three-dimensional rapid visualization of flame temperature field via compression and noise reduction of light field imaging
Liu et al. Spatial resolution of light field sectioning pyrometry for flame temperature measurement
Zhang et al. Summary on calibration method of line-structured light sensor
Yuan et al. Influence analysis of radiative properties and flame temperature reconstruction based on optical tomography
Li et al. Joint method for reconstructing three-dimensional temperature of flame using Lucy-Richardson and nearest neighbor filtering using light-field imaging
Shi et al. Reconstruction of 3D flame temperature and absorption coefficient field by the hybrid light-field imaging and laser extinction technique
Qi et al. Simultaneous reconstruction of flame temperature and soot volume fraction through weighted non-negative least squares and light field imaging techniques
Wang et al. A method for in situ measurement of directional and spatial radiosity distributions from complex-shaped solar thermal receivers
CN102564614A (zh) 激光光斑动态测量方法及测量仪
Wen et al. A hybrid LSQP algorithm for simultaneous reconstruction of the temperature and absorption coefficient field from the light-field image
CN110470618A (zh) 基于大气选择透过特性的单色仪光波长偏移量的检测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant