CN101806587A - 一种采用绝对相位测量的光学三维测量方法 - Google Patents
一种采用绝对相位测量的光学三维测量方法 Download PDFInfo
- Publication number
- CN101806587A CN101806587A CN 201010161971 CN201010161971A CN101806587A CN 101806587 A CN101806587 A CN 101806587A CN 201010161971 CN201010161971 CN 201010161971 CN 201010161971 A CN201010161971 A CN 201010161971A CN 101806587 A CN101806587 A CN 101806587A
- Authority
- CN
- China
- Prior art keywords
- pattern
- coded
- grating
- optical
- gray
- 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.)
- Granted
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 16
- 238000005259 measurement Methods 0.000 title claims abstract description 15
- 238000000691 measurement method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 23
- 238000010276 construction Methods 0.000 claims description 3
- 238000012217 deletion Methods 0.000 claims description 3
- 230000037430 deletion Effects 0.000 claims description 3
- 230000000977 initiatory effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 9
- 238000010586 diagram Methods 0.000 description 7
- 238000013461 design Methods 0.000 description 4
- 230000000903 blocking effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000000059 patterning Methods 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000008921 facial expression Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000011505 plaster Substances 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
Images
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明属于光学三维传感技术领域,具体涉及一种采用绝对相位测量的光学三维测量方法。本发明提供了一种采用单幅辅助编码图案的方法,它结合相位编码方案中截断相位分布特点及编码图案变形特点,使用3值灰度编码,提高解码可靠性,并且简化了解码算法。即采用正弦光栅或罗奇光栅图案加一幅空间灰度编码图案的方法实现绝对相位测量,进而计算出物体的三维形貌和空间位置。
Description
技术领域
本发明属于光学三维传感技术领域,具体涉及一种采用绝对相位测量的光学三维测量方法。
背景技术
基于结构光照明的光学三维传感技术在工业检测、质量控制、机器视觉、影视特技和生物医学等领域有广泛的应用。该技术的关键之一是投影平面与摄像机像面对应点的匹配,通常利用颜色、灰度、光栅相位或利用上述参数编码的符号或符号序列等实现对应点匹配。
已有的结构光照明三维传感技术,按是否与颜色有关,可分为基于颜色和基于灰度的编码方案;按用于编码的参数不同,可分为正弦光栅相位编码方案(又称正弦条纹投影方案)和光强编码方案。
相位编码方案主要有相位测量轮廓术(DMP)和傅立叶变换轮廓术(FTP)等。相位编码方案中需要获取光栅的绝对相位才能实现对应点匹配,则物体的三维形貌和空间位置都可以被测量。
目前常见的绝对相位获取方法有:
(1)空间相位展开法。投影单一频率的光栅,由某一相位展开起始点开始,采用适当的算法进行相位展开;
(2)时间相位展开法。投影两个以上不同频率的光栅,利用这些光栅相位之间的关系进行相位展开;
(3)利用辅助的编码图案确定光栅条纹级次,展开相位。
第(3)种方法中又有单幅和多幅编码图案之分。采用多幅编码图案可以获得可靠的测量结果,解码方法简单,但测量时间较长;而采用单幅编码图案可提高测量速度,但通常解码算法复杂,可靠性相对低。
发明内容
本发明提供了一种采用单幅辅助编码图案的方法,它结合相位编码方案中截断相位分布特点及编码图案变形特点,使用3值灰度编码,提高解码可靠性,并且简化了解码算法。即采用正弦光栅或罗奇光栅图案加一幅空间灰度编码图案的方法实现绝对相位测量,进而计算出物体的三维形貌和空间位置。
本发明的具体技术路线是:
一种采用绝对相位测量的光学三维测量方法,包括以下步骤:
1)、一配备光栅的投影仪向物体投射结构光图案;
2)、一摄像机在另一位置拍摄被物体反射回来的变形图像;
3)、对物体反射回来的变形图像进行解析,获得与摄像机像面上的点对应的投影平面上点的坐标;
4)、根据几何关系计算出物体测量点坐标进行三维重构,获得物体三维图像;
其特征在于
步骤1)所述的结构光图案包含光栅图案和空间灰度编码图案,
该空间灰度编码图案由以下步骤构建:
A)、空间灰度编码图案由若干与光栅条纹对应的条带构成,条带的宽度与条纹周期相同;
B、利用条带中沿条带方向的灰度值变化规律确定一组编码符号;
C)、由这组编码符号构建的伪随机序列确定空间灰度编码图案。
所述的伪随机序列构建方法包含以下步骤:
a)、从指定的一组编码符号中随机选取个数等于窗口长度的编码符号,构成初始序列;
b)、从该组编码符号中随机选取1个符号,添加到已构建的序列后形成新序列;
c)、判断该新序列是否满足以下条件;
大于等于窗口长度的子序列是唯一的,
在任意等于窗口长度的子序列内无重复编码符号,
如符合,则检查整个序列长度是否不小于光栅条纹数,如果满足,构建结束,否则转b)步,
如不符合,则删除刚添加的编码符号,转b)步。
所述空间灰度编码图案的灰度值包含3级。
图1为光栅投影三维测量***示意图。它由投影仪和摄像机构成。Op、Oc分别为投影仪和摄像机光心,(xp0,yp0)、(xc0,yc0)分别为投影仪和摄像机主点的图像坐标。Pp、Pc为一对投影平面和摄像机像面的对应点。P为物体表面上的测量点,其坐标为(Xw,Yw,Zw)。由投影仪投出正弦光栅或罗奇光栅,摄像机在另一位置拍摄被物体表面反射回来的变形条纹。通过求解拍摄到的图像绝对相位,获得与摄像机像面上的点对应的投影平面上点的坐标。最后在由POcOp构成的三角形中,根据几何关系计算出P点坐标,进而恢复出物体三维形貌和空间位置。
灰度空间编码图案的基本设计思想是:
首先将灰度空间编码图案划分成若干与投影条纹对应的条带,条带的宽度与条纹周期相同,即一条条带对应条纹的一个周期。
利用条带中沿条带方向的灰度分布规律来编码一组符号。
然后用由这组符号构成的伪随机序列来确定图案的总体结构。
其中,伪随机序列应满足如下条件:
(1)在序列中,大于等于特定长度(称为窗口长度)的子序列是唯一的;
(2)在任意长度等于窗口长度的子序列内无重复符号。比如窗口长度为4的序列“ABDECFADBE FDBECDABFECBDEFBDCEABCDAECFBDE”。
符号的编码方法为:根据沿条带方向灰度是否有变化、灰度值及变化幅值等规律进行编码。例如,采用3级灰度进行编码,表示为0,1,2。沿条带方向灰度无变化,也可以认为变化规律为:0-0,1-1,2-2,每种灰度值编码一个符号,可以编码3个符号。沿条带方向灰度周期性二值变化,根据这两个变化的灰度值不同,可以编码3个符号,分别是:0-1,0-2,1-2二值变化。综合上述两类情况,一共可以编码6个符号,可以分别表示为:A(0-0),B(1-1),C(2-2),D(0-1),E(0-2),F(1-2)。
图2为编码图案举例,由相邻的相位跳变绝对值超过π的点连接而成的光栅相位跳变线位于条带中心。
本发明采用的解码算法方案如下:
实际测量获得的编码图案如图3所示,图中颜色较亮的线为跳变线。由于相位跳变线位于条带中心,其灰度受相邻条带影响较小,而且条带随跳变线变形而变形。正弦光栅条纹对应的编码可以通过分析跳变线上编码图案灰度变化情况确定。解码算法如下:
(1)去除物体表面反射率对编码图像的影响;
(2)根据阈值对编码图像进行三值量化;
(3)沿跳变线计算灰度直方图,某一灰度的点所占比例大于阈值,就确认该灰度存在;
(4)根据编码规则恢复代码。
附图说明
图1为本发明光栅投影三维成像***示意图
图2为本发明举例的编码图案示意图
图3为本发明举例的实际测量编码图案示意图
图4为本发明正弦及编码光栅示意图
图5为本发明人脸表情变化测量结果示意图
图6为本发明复杂场景测量结果示意图
具体实施方式
实施例1
一种采用绝对相位测量的光学三维测量方法,包括以下步骤:
1)、一配备光栅的投影仪向物体投射结构光图案;
2)、一摄像机在另一位置拍摄被物体反射回来的变形图像;
3)、对物体反射回来的变形图像进行解析,获得与摄像机像面上的点对应的投影平面上点的坐标;
4)、根据几何关系计算出物体测量点坐标进行三维重构,获得物体三维图像;
其特征在于
步骤1)所述的结构光图案包含光栅图案和空间灰度编码图案,
该空间灰度编码图案由以下步骤构建:
A)、空间灰度编码图案由若干与投影条纹对应的条带构成,条带的宽度与条纹周期相同;
B、利用条带中沿条带方向的灰度值变化规律确定一组编码符号;
C)、由这组编码符号构建的伪随机序列确定空间灰度编码图案。
所述的伪随机序列构建方法包含以下步骤:
a)、从指定的一组编码符号中随机选取个数等于窗口长度的编码符号,构成初始序列;
b)、从该组编码符号中随机选取1个符号,添加到已构建的序列后形成新序列;
c)、判断该新序列是否满足以下条件;
大于等于窗口长度的子序列是唯一的,
在任意等于窗口长度的子序列内无重复编码符号,
如符合,则检查整个序列长度是否不小于光栅条纹数,如果满足,构建结束,否则转b)步,
如不符合,则删除刚添加的编码符号,转b)步。
所述空间灰度编码图案的灰度值包含3级。
具体实验步骤如下:
根据投影仪的像素列数、光栅周期计算伪随机序列长度,采用前述方法构造伪随机序列;
根据前述的空间灰度编码图案编码方法,设计与每个符号对应的编码条带;
按照构造的代码序列,将编码单元按顺序排列,最后得到与伪随机序列对应的编码图案;
将设计得到的光栅和编码图案投影到被测物体上,并同时用摄像机拍摄被物体形貌调制的变形光栅和编码图案;
在投影完所有图案后,求解光栅的截断相位,并结合得到的截断相位将拍摄到的编码图案解码,计算绝对相位,最后得到物体的三维形貌和空间位置。
采用窗口长度为3的序列:
“ABDECFACDCEFBEADFEBCAEDAFDBFCBADEFCABEDC F”,符号用3值灰度编码。首先用π相移傅立叶变换轮廓术结合编码图案进行动态过程测量实验,光栅周期为18像素。记录速度为每秒60帧三维数据。图4(a)、图4(b)、图4(c)为设计的一套光栅,图4(a)为正弦光栅,图4(b)为π相移后的正弦光栅,图4(c)为编码光栅。
图5为笑过程中人脸变化测量结果。然后采用五步相移法结合编码图案测量了由若干孤立物体组成的场景。图6(a)为被测场景照片,其中墙和石膏像距离约690mm;图6(b)为测量结果;图6(c)为去除墙后的结果。可见采用本方法能够可靠地测量复杂场景的三维形貌和空间位置,并有较大的测量景深。
Claims (3)
1.一种采用绝对相位测量的光学三维测量方法,包括以下步骤:
1)、一配备光栅的投影仪向物体投射结构光图案;
2)、一摄像机在另一位置拍摄被物体反射回来的变形图像;
3)、对物体反射回来的变形图像进行解析,获得与摄像机像面上的点对应的投影平面上点的坐标;
4)、根据几何关系计算出物体测量点坐标进行三维重构,获得物体三维图像;
其特征在于
步骤1)所述的结构光图案包含光栅图案和空间灰度编码图案,该空间灰度编码图案由以下步骤构建:
A)、空间灰度编码图案由若干与光栅条纹对应的条带构成,条带的宽度与条纹周期相同;
B)、利用条带中沿条带方向的灰度值变化规律确定一组编码符号;
C)、由这组编码符号构建的伪随机序列确定空间灰度编码图案。
2.根据权利要求1所述的光学三维测量方法,其特征在于所述的伪随机序列构建方法包含以下步骤:
a)、从指定的一组编码符号中随机选取个数等于窗口长度的编码符号,构成初始序列;
b)、从该组编码符号中随机选取1个符号,添加到已构建的序列后形成新序列;
c)、判断该新序列是否满足以下条件;
大于等于窗口长度的子序列是唯一的,
在任意等于窗口长度的子序列内无重复编码符号,
如符合,则检查整个序列长度是否不小于光栅条纹数,
如果满足,构建结束,否则转b)步,
如不符合,则删除刚添加的编码符号,转b)步。
3.根据权利要求1所述的光学三维测量方法,其特征在于所述空间灰度编码图案包含三级灰度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010161971 CN101806587B (zh) | 2010-04-29 | 2010-04-29 | 一种采用绝对相位测量的光学三维测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010161971 CN101806587B (zh) | 2010-04-29 | 2010-04-29 | 一种采用绝对相位测量的光学三维测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101806587A true CN101806587A (zh) | 2010-08-18 |
CN101806587B CN101806587B (zh) | 2012-12-12 |
Family
ID=42608478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010161971 Expired - Fee Related CN101806587B (zh) | 2010-04-29 | 2010-04-29 | 一种采用绝对相位测量的光学三维测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101806587B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102012892A (zh) * | 2010-12-13 | 2011-04-13 | 上海大学 | 可重构的三值光学处理器 |
CN102243138A (zh) * | 2011-08-05 | 2011-11-16 | 中国科学院光电技术研究所 | 一种用于投影光刻中的焦面检测装置 |
CN102322823A (zh) * | 2011-09-13 | 2012-01-18 | 四川大学 | 基于相位级次自编码的光学三维测量方法 |
WO2012071817A1 (zh) * | 2010-11-30 | 2012-06-07 | 东南大学 | 基于小波变换的三维测量方法 |
CN103398675A (zh) * | 2013-08-10 | 2013-11-20 | 南昌航空大学 | 一种基于条纹周期校正的复杂大物体三维测量方法 |
CN104132628A (zh) * | 2014-07-27 | 2014-11-05 | 四川大学 | 用相位计算实现线结构光三维测量的方法 |
CN104180773A (zh) * | 2013-05-20 | 2014-12-03 | 沈阳同联集团高新技术有限公司 | 一种三维形貌测量装置 |
CN104697468A (zh) * | 2015-03-30 | 2015-06-10 | 南昌航空大学 | 基于倾斜相位编码条纹的光学三维测量方法 |
CN106840039A (zh) * | 2017-03-10 | 2017-06-13 | 合肥工业大学 | 一种基于复合光栅投影的快速三维测量方法 |
CN106910213A (zh) * | 2017-01-23 | 2017-06-30 | 四川大学 | 一种基于计算成像的光学不可视场景三维信息获取方法 |
CN109945802A (zh) * | 2018-10-11 | 2019-06-28 | 宁波深浅优视智能科技有限公司 | 一种结构光三维测量方法 |
CN109974626A (zh) * | 2019-04-08 | 2019-07-05 | 四川大学 | 一种基于相移量编码条纹级次的结构光三维测量方法 |
CN110986828A (zh) * | 2019-12-13 | 2020-04-10 | 杨鑫鑫 | 一种新型真实场景三维彩色数据的采集显示方法 |
TWI699525B (zh) * | 2019-05-21 | 2020-07-21 | 中華學校財團法人中華科技大學 | 三維相位移瑕疵檢測方法及系統 |
CN111947601A (zh) * | 2020-08-12 | 2020-11-17 | 上海科技大学 | 一种灰度伪随机编码结构光条纹的投影解算方法 |
EP3577416A4 (en) * | 2017-02-06 | 2020-12-16 | Modit3d, Inc. | SYSTEM AND PROCEDURE FOR 3D SCANNING |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101504277A (zh) * | 2009-02-26 | 2009-08-12 | 浙江师范大学 | 一种光学三维传感获取物体三维图像的方法 |
-
2010
- 2010-04-29 CN CN 201010161971 patent/CN101806587B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101504277A (zh) * | 2009-02-26 | 2009-08-12 | 浙江师范大学 | 一种光学三维传感获取物体三维图像的方法 |
Non-Patent Citations (1)
Title |
---|
《中国优秀硕士学位论文全文数据库基础科学辑》 20080215 纪文丽 基于光栅投影方法的物体三维轮廓重构 第3.4节 1-3 , 第2期 2 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9163937B2 (en) | 2010-11-30 | 2015-10-20 | Southeast University | Three-dimensional measurement method based on wavelet transform |
WO2012071817A1 (zh) * | 2010-11-30 | 2012-06-07 | 东南大学 | 基于小波变换的三维测量方法 |
CN102012892A (zh) * | 2010-12-13 | 2011-04-13 | 上海大学 | 可重构的三值光学处理器 |
CN102243138A (zh) * | 2011-08-05 | 2011-11-16 | 中国科学院光电技术研究所 | 一种用于投影光刻中的焦面检测装置 |
CN102322823A (zh) * | 2011-09-13 | 2012-01-18 | 四川大学 | 基于相位级次自编码的光学三维测量方法 |
CN104180773A (zh) * | 2013-05-20 | 2014-12-03 | 沈阳同联集团高新技术有限公司 | 一种三维形貌测量装置 |
CN103398675B (zh) * | 2013-08-10 | 2016-06-15 | 南昌航空大学 | 一种基于条纹周期校正的复杂大物体三维测量方法 |
CN103398675A (zh) * | 2013-08-10 | 2013-11-20 | 南昌航空大学 | 一种基于条纹周期校正的复杂大物体三维测量方法 |
CN104132628A (zh) * | 2014-07-27 | 2014-11-05 | 四川大学 | 用相位计算实现线结构光三维测量的方法 |
CN104697468A (zh) * | 2015-03-30 | 2015-06-10 | 南昌航空大学 | 基于倾斜相位编码条纹的光学三维测量方法 |
CN106910213B (zh) * | 2017-01-23 | 2019-09-03 | 四川大学 | 一种基于计算成像的光学不可视场景三维信息获取方法 |
CN106910213A (zh) * | 2017-01-23 | 2017-06-30 | 四川大学 | 一种基于计算成像的光学不可视场景三维信息获取方法 |
US11330243B2 (en) | 2017-02-06 | 2022-05-10 | Riven, Inc. | System and method for 3D scanning |
EP3577416A4 (en) * | 2017-02-06 | 2020-12-16 | Modit3d, Inc. | SYSTEM AND PROCEDURE FOR 3D SCANNING |
CN106840039A (zh) * | 2017-03-10 | 2017-06-13 | 合肥工业大学 | 一种基于复合光栅投影的快速三维测量方法 |
CN109945802B (zh) * | 2018-10-11 | 2021-03-09 | 苏州深浅优视智能科技有限公司 | 一种结构光三维测量方法 |
CN109945802A (zh) * | 2018-10-11 | 2019-06-28 | 宁波深浅优视智能科技有限公司 | 一种结构光三维测量方法 |
CN109974626A (zh) * | 2019-04-08 | 2019-07-05 | 四川大学 | 一种基于相移量编码条纹级次的结构光三维测量方法 |
CN109974626B (zh) * | 2019-04-08 | 2021-02-09 | 四川大学 | 一种基于相移量编码条纹级次的结构光三维测量方法 |
TWI699525B (zh) * | 2019-05-21 | 2020-07-21 | 中華學校財團法人中華科技大學 | 三維相位移瑕疵檢測方法及系統 |
CN110986828A (zh) * | 2019-12-13 | 2020-04-10 | 杨鑫鑫 | 一种新型真实场景三维彩色数据的采集显示方法 |
CN110986828B (zh) * | 2019-12-13 | 2023-09-01 | 杨鑫鑫 | 一种新型真实场景三维彩色数据的采集显示方法 |
CN111947601A (zh) * | 2020-08-12 | 2020-11-17 | 上海科技大学 | 一种灰度伪随机编码结构光条纹的投影解算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101806587B (zh) | 2012-12-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101806587B (zh) | 一种采用绝对相位测量的光学三维测量方法 | |
CN101813461B (zh) | 基于复合彩色条纹投影的绝对相位测量方法 | |
CN101694375B (zh) | 一种用于强反射表面三维形貌测量的立体视觉检测方法 | |
CN101504277B (zh) | 一种光学三维传感获取物体三维图像的方法 | |
Pages et al. | Optimised De Bruijn patterns for one-shot shape acquisition | |
CN102322823B (zh) | 基于相位级次自编码的光学三维测量方法 | |
CN110645919B (zh) | 一种基于空域二值编码的结构光三维测量方法 | |
CN109341589B (zh) | 一种光栅图像投影方法、三维重建方法及三维重建*** | |
Song et al. | An accurate and robust strip-edge-based structured light means for shiny surface micromeasurement in 3-D | |
CN106643562B (zh) | 基于时域空域混合编码的结构光条纹投射方法 | |
CN110174079B (zh) | 一种基于四步相移编码型面结构光的三维重建方法 | |
EP2881702B1 (en) | Image processing system, and image processing method | |
Song et al. | A new phase unwrapping algorithm based on three wavelength phase shift profilometry method | |
CN104567730B (zh) | 一种时空二元编码产生正弦结构光场的方法 | |
CN108955571A (zh) | 双频外差与相移编码相结合的三维测量方法 | |
CN104390608A (zh) | 基于投影栅相位法的结构光三维形状构建方法 | |
CN105303572B (zh) | 基于主被动结合的深度信息获取方法 | |
CN110849290B (zh) | 基于形态学操作的分段量化编码强度的三维测量方法 | |
CN111207692A (zh) | 一种改进的分段阶梯相位编码三维测量方法 | |
CN103453852A (zh) | 三维扫描***中的快速相位展开方法 | |
CN113345039B (zh) | 三维重建量化结构光相位图像编码方法 | |
CN107421467B (zh) | 结构光三维测量方法及*** | |
US10801834B2 (en) | Fringe projection for determining topography of a body | |
Boukamcha et al. | Robust technique for 3D shape reconstruction | |
Zhang et al. | Determination of edge correspondence using color codes for one-shot shape acquisition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121212 Termination date: 20150429 |
|
EXPY | Termination of patent right or utility model |