CN113034532B - 一种基于无网格模型的整形手术术后软组织形变预测方法 - Google Patents
一种基于无网格模型的整形手术术后软组织形变预测方法 Download PDFInfo
- Publication number
- CN113034532B CN113034532B CN202110232153.9A CN202110232153A CN113034532B CN 113034532 B CN113034532 B CN 113034532B CN 202110232153 A CN202110232153 A CN 202110232153A CN 113034532 B CN113034532 B CN 113034532B
- Authority
- CN
- China
- Prior art keywords
- soft tissue
- mesh
- plastic
- extracting
- deformation
- 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
Links
- 210000004872 soft tissue Anatomy 0.000 title claims abstract description 109
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000002316 cosmetic surgery Methods 0.000 title claims abstract description 28
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 23
- 238000006073 displacement reaction Methods 0.000 claims abstract description 18
- 238000004088 simulation Methods 0.000 claims abstract description 13
- 238000009499 grossing Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 4
- 230000037182 bone density Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 9
- 238000001356 surgical procedure Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000002980 postoperative effect Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007493 shaping process Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 210000002615 epidermis Anatomy 0.000 description 1
- 230000003340 mental effect Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000003491 skin Anatomy 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4023—Scaling of whole images or parts thereof, e.g. expanding or contracting based on decimating pixels or lines of pixels; based on inserting pixels or lines of pixels
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及医学整形预测领域,公开了一种基于无网格模型的整形手术术后软组织形变预测方法,用以提升整形手术效果预测的精准程度。本发明包括;S1.从患者术前整形部CT数据中提取受手术影响区域的软组织的内、外表面;S2.在软组织内、外表面之间提取手术影响范围内的无网格软组织节点;S3.采用无网格建模方法,构建软组织节点的无网格形函数及其偏导数;S4.根据骨骼位移量对软组织节点进行无网格形变仿真;S5.通过对软组织外表面的节点位移进行插值和平滑,得到整个软组织外表面的位置仿真预测值。本发明适用于颌面整形手术、正畸手术等各种通过矫正骨骼来引起软组织外形改变的整形手术术后软组织形变预测。
Description
技术领域
本发明涉及医学整形预测领域,特别涉及一种基于无网格模型的整形手术术后软组织形变预测方法。
背景技术
随着经济发展和人们精神文化水平的提高,要求进行体态整形手术的患者数量逐年增多,颌面骨骼整形手术就是其中比较有代表性的一种。在颌面骨骼整形手术中,患者最关心的问题之一就是术后自己的面貌将得到多大的改善,而目前商用的整形手术仿真软件主要致力于截骨、正畸等硬组织的手术仿真及方案设计,而对于术后的软组织和面貌改变仅基于计算机图形学进行面部3D图像的轮廓形变预测,与真实的术后效果存在较大差异。经发明人研究,发现造成前述问题的主要原因是,现有技术缺乏对患者术后效果的客观量化的定性分析预测,而且容易出现术后真实效果与患者自己期望的改善效果相距甚远、引发患者对治疗结果不满等情况。
目前,公开号为CN107134010的专利涉及到使用有限元方法对***弹性软组织进行形貌预测,但是有限元方法基于网格模型建立,网格模型在在处理网格的初始化、大变形(如手术方案对骨骼进行较大位移)、网格分布不连续(如颌面部的口鼻处会破坏软组织网格的连续性)等问题时存在很大的困难,很容易出现病态网格,影响计算的精度和稳定性。
发明内容
本发明要解决的技术问题是:提供一种基于无网格模型的整形手术术后软组织形变预测方法,用以提升整形手术效果预测的精准程度。
为实现上述目的,本发明采用的技术方案是:一种基于无网格模型的整形手术术后软组织形变预测方法,包括以下步骤:
S1.从患者术前整形部CT数据中提取受手术影响区域的软组织的内、外表面;
S2.在软组织内、外表面之间提取手术影响范围内的无网格软组织节点;
S3.采用无网格建模方法,构建软组织节点的无网格形函数及其偏导数;
S4.根据骨骼位移量对软组织节点进行无网格形变仿真;
S5.通过对软组织外表面的节点位移进行插值和平滑,得到整个软组织外表面的位置仿真预测值。
具体的,所述整形部可以为颌面部。
进一步的,步骤S1具体包括:
S1-1.从患者术前整形部CT影像数据文件中提取三维体素数据,并将三维体素的各个维度调整为统一的分辨率;
S1-2.根据软组织和骨骼在整形部CT中的密度差异,设定骨骼和软组织的密度分割阈值,利用三维等值面提取算法,提取三维体素在骨骼密度阈值下的等值面作为软组织内表面surface1,提取三维体素在软组织密度阈值下的等值面作为软组织外表面surface2。
进一步的,步骤S2具体包括:
S2-1.根据手术方案,在三维体素数据中确定接受整形部手术及受该手术影响的骨骼区域area1的三维坐标范围;
S2-2.设定受手术影响的软组织距离范围阈值threshold1,将步骤S2-1中得到的骨骼区域三维坐标范围向外扩充距离threshold1,得到受整形部手术影响的软组织区域坐标范围area2;
S2-3.在上述受整形部手术影响的软组织区域坐标范围area2内,按照设定的节点间距阈值threshold2提取位于步骤S1-2中所述软组织内表面上、软组织外表面上以及软组织内、外表面之间的软组织节点。
进一步的,步骤S3中,形函数及其偏导数的计算采用无网格移动最小二乘法或无网格径向基插值法。
进一步的,步骤S4具体包括:
S4-1.利用无网格法构建全部软组织节点的弹性力学总体刚度矩阵K;
S4-2.以软组织内表面上的节点位移量作为本质边界条件构建无网格法的弹性力学惩罚刚度矩阵Ka和惩罚外力矩阵Fa;
S4-3.构建弹性力学离散控制方程(K+Ka)U=Fa,利用最小二乘法解算出式中的全部节点位移向量U。
本发明的有益效果是:本发明基于无网格建模方法,通过对颌面整形后面部软组织的位移变化进行精确的生物力学建模仿真,具有方法新颖、结果可靠、预测精准等优点,可提升整形手术效果预测的精准程度,为医患提供重要的术前参考。本发明可适用于颌面整形手术、正畸手术等各种通过矫正骨骼来引起软组织外形改变的整形手术术后软组织形变预测。
附图说明
图1为实施例提供的一种基于无网格模型的颌面整形手术术后软组织形变预测流程图。
图2为实施例从dicom中提取得到的软组织上表面surface1和软组织下表面surface2示意图。
图3为实施例提取出的软组织无网格节点示意图。
图4为实施例根据仿真结果预测的软组织上表面示意图。
具体实施方式
为了提升整形手术效果预测的精准程度,本发明提供了一种基于无网格模型的整形手术术后软组织形变预测方法,包括以下步骤:
S1.从患者术前整形部CT数据中提取受手术影响区域的软组织的内、外表面;
S2.在软组织内、外表面之间提取手术影响范围内的无网格软组织节点;
S3.采用无网格建模方法,构建软组织节点的无网格形函数及其偏导数;
S4.根据骨骼位移量对软组织节点进行无网格形变仿真;
S5.通过对软组织外表面的节点位移进行插值和平滑,得到整个软组织外表面的位置仿真预测值。
该方法通过无网格方法对颌面部软组织在术后的形变进行仿真,提升手术效果预测的精准程度,为医患提供重要的术前参考。
下面以颌面整形手术术后软组织形变预测为例对本发明技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一种优选实施例,并不用于限制本发明,对于本领域的技术人员来说,对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
如图1所示,实施例提供的一种基于无网格模型的颌面整形手术术后软组织形变预测方法,其具体实施步骤如下:
S1,从患者术前CT数据中提取受手术影响区域的软组织的内、外表面:
S1-1,从患者术前头部CT影像数据dicom文件中提取三维体素数据,并根据dicom数据中的像素分辨率信息和层厚信息,将三维体素的各个维度调整为统一的分辨率,如1mm*1mm*1mm;
S1-2,根据软组织和骨骼在CT中的密度差异,设定骨骼和软组织的密度分割阈值,典型的骨骼密度阈值如2000,典型的软组织密度分割阈值如800,利用三维等值面提取算法如Dual Contouring算法,提取三维体素在骨骼密度阈值下的等值面作为软组织内表面surface1,提取三维体素在软组织密度阈值下的等值面作为软组织外表面surface2,surface1对应附着在骨骼表面的最靠内的软组织层,surface2对应皮肤的表皮层,如图2所示。
S2,在软组织内、外表面之间提取手术影响范围内的无网格软组织节点:
S2-1,根据手术方案,在三维体素数据中确定接受颌面部手术及受该手术影响的骨骼区域area1的三维坐标范围,例如,在颏成型手术中,该坐标范围即包含手术方案所切割和移动的颏顶部骨骼,而在下颌支倒L型骨切开手术中,该坐标范围应不仅包含切开和移动的倒L型骨部分,还应该包含随着倒L型骨同时产生位移的整个下颌骨部分;
S2-2,设定受手术影响的软组织距离范围阈值threshold1,典型的阈值为30mm~50mm,将步骤S2-1中得到的骨骼区域三维坐标范围向外扩充距离threshold1,得到受颌面部手术影响的软组织区域坐标范围area2;
S2-3,在上述受颌面部手术影响的软组织区域坐标范围area2内,按照设定的节点间距阈值threshold2提取位于步骤S1-2中所述软组织内表面上、软组织外表面上以及软组织内、外表面之间的软组织节点P,如图3所示。
S3,采用无网格建模方法,构建软组织节点的无网格形函数及其偏导数,形函数及其偏导数的计算通常可采用无网格移动最小二乘法(MLS)或无网格径向基插值法(RPIM);
S4,根据骨骼位移量对软组织节点进行无网格形变仿真:
S4-1,利用无网格法构建全部软组织节点的弹性力学总体刚度矩阵K;
S4-2,以软组织内表面上的节点位移量作为本质边界条件构建无网格法的弹性力学惩罚刚度矩阵Ka和惩罚外力矩阵Fa;
S4-3,构建弹性力学离散控制方程(K+Ka)U=Fa,利用最小二乘法解算出式中的全部节点位移向量U。
S5,通过对软组织外表面的节点位移进行插值和平滑,得到整个软组织外表面的位置仿真预测值,得到的预测软组织外表面如图4所示。
Claims (6)
1.一种基于无网格模型的整形手术术后软组织形变预测方法,其特征在于,包括以下步骤:
S1.从患者术前整形部CT数据中提取受手术影响区域的软组织的内、外表面;
S2.在软组织内、外表面之间提取手术影响范围内的无网格软组织节点;
S3.采用无网格建模方法,构建软组织节点的无网格形函数及其偏导数;
S4.根据骨骼位移量对软组织节点进行无网格形变仿真;
S5.通过对软组织外表面的节点位移进行插值和平滑,得到整个软组织外表面的位置仿真预测值。
2.如权利要求1所述的一种基于无网格模型的整形手术术后软组织形变预测方法,其特征在于,所述整形部为颌面部。
3.如权利要求1或2所述的一种基于无网格模型的整形手术术后软组织形变预测方法,其特征在于,步骤S1具体包括:
S1-1.从患者术前整形部CT影像数据文件中提取三维体素数据,并将三维体素的各个维度调整为统一的分辨率;
S1-2.根据软组织和骨骼在整形部CT中的密度差异,设定骨骼和软组织的密度分割阈值,利用三维等值面提取算法,提取三维体素在骨骼密度阈值下的等值面作为软组织内表面surface1,提取三维体素在软组织密度阈值下的等值面作为软组织外表面surface2。
4.如权利要求3所述的一种基于无网格模型的整形手术术后软组织形变预测方法,其特征在于,步骤S2具体包括:
S2-1.根据手术方案,在三维体素数据中确定接受整形部手术及受该手术影响的骨骼区域area1的三维坐标范围;
S2-2.设定受手术影响的软组织距离范围阈值threshold1,将步骤S2-1中得到的骨骼区域三维坐标范围向外扩充距离threshold1,得到受整形部手术影响的软组织区域坐标范围area2;
S2-3.在上述受整形部手术影响的软组织区域坐标范围area2内,按照设定的节点间距阈值threshold2提取位于步骤S1-2中所述软组织内表面上、软组织外表面上以及软组织内、外表面之间的软组织节点。
5.如权利要求1或2所述的一种基于无网格模型的整形手术术后软组织形变预测方法,其特征在于,步骤S3中,形函数及其偏导数的计算采用无网格移动最小二乘法或无网格径向基插值法。
6.如权利要求1或2所述的一种基于无网格模型的整形手术术后软组织形变预测方法,其特征在于,步骤S4具体包括:
S4-1.利用无网格法构建全部软组织节点的弹性力学总体刚度矩阵K;
S4-2.以软组织内表面上的节点位移量作为本质边界条件构建无网格法的弹性力学惩罚刚度矩阵Ka和惩罚外力矩阵Fa;
S4-3.构建弹性力学离散控制方程(K+Ka)U=Fa,利用最小二乘法解算出式中的全部节点位移向量U。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110232153.9A CN113034532B (zh) | 2021-03-02 | 2021-03-02 | 一种基于无网格模型的整形手术术后软组织形变预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110232153.9A CN113034532B (zh) | 2021-03-02 | 2021-03-02 | 一种基于无网格模型的整形手术术后软组织形变预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113034532A CN113034532A (zh) | 2021-06-25 |
CN113034532B true CN113034532B (zh) | 2023-02-03 |
Family
ID=76465974
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110232153.9A Active CN113034532B (zh) | 2021-03-02 | 2021-03-02 | 一种基于无网格模型的整形手术术后软组织形变预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113034532B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102044086A (zh) * | 2010-11-30 | 2011-05-04 | 华北水利水电学院 | 一种软组织形变仿真方法 |
EP2789308A1 (en) * | 2013-04-12 | 2014-10-15 | Stryker Leibinger GmbH & Co. KG | Computer-implemented technique for generating a data set that geometrically defines a bone cut configuration |
CN106504329A (zh) * | 2016-09-27 | 2017-03-15 | 西安科技大学 | 一种基于牙体长轴的质点弹簧模型的牙龈变形仿真方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6898263B2 (en) * | 2002-11-27 | 2005-05-24 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for soft-tissue volume visualization |
CN105303605B (zh) * | 2015-10-26 | 2017-10-13 | 哈尔滨理工大学 | 一种基于力反馈的骨外科手术仿真*** |
CN105608741A (zh) * | 2015-12-17 | 2016-05-25 | 四川大学 | 预测颌面骨整形术后软组织容貌变化的计算机仿真方法 |
CN106570341B (zh) * | 2016-11-14 | 2019-06-18 | 南昌大学 | 一种基于无网格径向基数据拟合的软组织受力形变模型建模方法 |
CN106875462B (zh) * | 2017-01-13 | 2020-07-07 | 北京航空航天大学 | 一种基于元球模型和混合驱动方法的实时数字器官切割方法 |
CN107134010A (zh) * | 2017-04-27 | 2017-09-05 | 杭州电子科技大学 | 一种基于有限元的弹性软组织的形貌效果预测方法 |
CN108836522B (zh) * | 2018-03-20 | 2020-07-07 | 北京正齐口腔医疗技术有限公司 | 牙龈变形的仿真方法及装置 |
CN108536936A (zh) * | 2018-03-27 | 2018-09-14 | 南京信息工程大学 | 一种多重优化的无网格软组织形变模拟方法 |
CN108710735A (zh) * | 2018-05-08 | 2018-10-26 | 南京信息工程大学 | 一种实时交互的无网格软组织形变模拟方法 |
CN111127477B (zh) * | 2019-12-12 | 2023-10-20 | 明峰医疗***股份有限公司 | 一种ct多色谱体素仿真方法 |
-
2021
- 2021-03-02 CN CN202110232153.9A patent/CN113034532B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102044086A (zh) * | 2010-11-30 | 2011-05-04 | 华北水利水电学院 | 一种软组织形变仿真方法 |
EP2789308A1 (en) * | 2013-04-12 | 2014-10-15 | Stryker Leibinger GmbH & Co. KG | Computer-implemented technique for generating a data set that geometrically defines a bone cut configuration |
CN106504329A (zh) * | 2016-09-27 | 2017-03-15 | 西安科技大学 | 一种基于牙体长轴的质点弹簧模型的牙龈变形仿真方法 |
Non-Patent Citations (1)
Title |
---|
虚拟手术中软组织形变仿真研究综述;郭煜等;《集成技术》;20130315(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113034532A (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Marchetti et al. | Validation of new soft tissue software in orthognathic surgery planning | |
Chabanas et al. | Patient specific finite element model of the face soft tissues for computer-assisted maxillofacial surgery | |
JP5976730B2 (ja) | 形状情報を引き出すための方法 | |
Schendel et al. | 3-dimensional facial simulation in orthognathic surgery: is it accurate? | |
CN105608741A (zh) | 预测颌面骨整形术后软组织容貌变化的计算机仿真方法 | |
Makram et al. | Reeb graph for automatic 3D cephalometry | |
CN102254317A (zh) | 牙种植导航中牙弓曲面的自动提取方法 | |
CN110236673B (zh) | 一种基于数据库的双侧颌骨缺损重建术前设计方法及装置 | |
Moiduddin et al. | Computer assisted design and analysis of customized porous plate for mandibular reconstruction | |
Zachow | Computational planning in facial surgery | |
CN107874831A (zh) | 一种基于隐函数的颅颌面部导板设计方法 | |
Zachow et al. | Computer assisted planning in cranio-maxillofacial surgery | |
CN113034532B (zh) | 一种基于无网格模型的整形手术术后软组织形变预测方法 | |
CN114863056B (zh) | 基于面型参数生成颞下颌关节髁突运动包络面方法和装置 | |
CN112017275B (zh) | 一种面部正中矢状参考平面的辅助定位方法 | |
CN113052864B (zh) | 一种基于机器学习的整形手术术后体貌预测方法 | |
Li et al. | Accuracy of three-dimensional facial soft tissue simulation in post-traumatic zygoma reconstruction | |
KR101863956B1 (ko) | 두개악안면 결손 보철 모델 생성 방법, 장치, 컴퓨터 프로그램 및 두개악안면 결손 보철 제조 방법 | |
Chabanas et al. | Evaluating soft tissue simulation in maxillofacial surgery using preoperative and postoperative CT scans | |
Ulusoy et al. | Use of the dynamic volume spline method to predict facial soft tissue changes associated with orthognathic surgery | |
Liao et al. | Physical modeling with orthotropic material based on harmonic fields | |
Chabanas et al. | Models for planning and simulation in computer assisted orthognatic surgery | |
Fang et al. | Tracked motion-based dental occlusion surface estimation for crown restoration | |
Chabanas et al. | Computer aided planning for orthognatic surgery | |
Gladilin et al. | A framework for biomechanical simulation of cranio-maxillofacial surgery interventions |
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 |