模拟面部皮肤老化和去老化的方法及装置

摘要

本发明公开了一种逼真地模拟促成皮肤的总体外观和状况的面部皮肤特征的进展或恶化的新方法及系统。所述方法利用面部的两种特写照片，一种是用数字照相机以标准白光采集，另一种是用同一个照相机以紫外光采集。接着，所述方法处理这些图像以模拟下列主要皮肤特征的进展或恶化：色素沉着斑点、皱纹和小纹理特征。这些特征的恶化将模拟由于长时间地暴露于阳光下而产生的面部皮肤老化、生物性老化或皮肤健康的恶化。这些特征的进展模拟在总体外观和健康方面面部皮肤的改善，如同患者已经经过治疗一样。因此，本发明公开了用于皮肤病学、化妆品和计算机动画的一系列方法。

说明书

发明领域

本发明涉及图像处理和模拟领域，具体地讲涉及用于描绘模拟皮肤老化或去老化的图像生成。

发明背景

皮肤老化对人脸外观的影响在皮肤病学中已得到很好的研究和证明。每个人的皮肤老化进展取决于内在和外在两种因素。内在因素例如性别、种族和肤色是由遗传所控制并且对于每个人是独一无二的，并可影响表皮变薄的速率、机械弹性的丧失和随着年龄而发生的其它具有鲜明特征的组织学变化和生物机械变化。内在因素影响阳光照不到的和暴露于阳光下的身体部位。外在因素包括个体的饮食、生活方式、皮肤护理习惯和日晒历史。已熟知长期日晒会加速皮肤老化的起始时间和严重程度。包括面部在内的所有暴露的身体部位均具有某种程度的皮肤光老化。(Gilchrest.，B.Photodamage，Blackwell Science，Inc.1995年)。

光老化的皮肤的视觉上最突出的特征之一是斑点纹状和不规则色素沉着，此在皮肤上表现为具有深棕色的斑点(Griffiths C.E.M.，“Theclinical identification and quantification of photodamage”，Brit.J.Derm.，第127卷(增刊41)，第37至42页，1992年；K.Miyamoto等人的“Utilization of a high-resolution digitalimaging system for the objective and quantitative assessment ofhyperpigmented spots on the face”，Skin Research andTechnology，第8卷，第2期，第73至78页，2002年5月，后文称为“Miyamoto参考文献”)。这些色素沉着损害称为老年斑、肝斑、老年斑痣或光化痣。采用会显示出在标准白光下不可见的亚表面色素沉着的方法可较好地观测光损坏的皮肤上的色素沉着过度。一种称为紫外激发荧光照相术的方法，其最初由Kollias提出(ollias等人的“Fluorescencephotography in the evaluation of hyperpigmentation inphotodamaged skin”，J Am Acad Dermatol.，第36卷，第226至230页，1997年)，涉及在中心为365nm的窄波段UVA下对皮肤进行成像。表皮黑素在这个UVA范围内强烈吸收，约为它在可见光谱中吸收性的3至5倍。未被表皮黑素所吸收的任何UVA进入真皮，在此处它被胶原和弹性蛋白纤维散射和吸收，所述胶原和弹性蛋白纤维将所吸收的一些能量转换成荧光。最大胶原放射的波长发生在中心为420nm的可见光谱中。黑素在420nm处的体内吸收性比在540nm处大一倍。因此，进入皮肤并达到真皮的UVA的总量被表皮黑素衰减5折并且可见荧光的量被同一表皮黑素衰减约2折。换句话讲，与可见光相比，利用紫外激发荧光的表皮黑素检测敏感约10倍。敏感性上的这种增强为检测在正常白光成像方法下不能看见的色素沉着点创造了条件。用可见光不能观察的色素沉着点在生命的后半段在正常可见光下将不受干预地变得更暗并更明显可见。

老化的皮肤的其它突出特征是部分地由真皮结缔组织例如胶原的逐渐蚀变和损失引起的粗糙纹理和皮肤皱纹(Leyden J.J.“Clinicalfeatures of ageing skin”，Br.J.Dermatol.第122卷，增刊35，第1至3页，1990年)，特别是在身体的日晒区域(Bailey，Molecularmechanisms of aging in connective tissues，Mech.Aging Dev.，第112卷，第7期，第735至755页，2001年)。色素沉着过度、皱纹和粗糙纹理是可见的皮肤特征，其对皮肤的总体外观和健康起着重要的作用。

能够准确地模拟老化过程具有实际价值。老化模拟具有几种有用的应用，例如计算机动画、面部识别、失踪人员鉴定、娱乐、医药和化妆品。已经利用各种模型来实现老化面部的逼真模拟，包括几何模型、基于物理的模型、基于图像的模型或生物-机械模型(Hussein，K.H，Towardrealistic facial modeling and re-rendering of human skin aginganimation，Proceedings of the Shape Modeling International2002，IEEE计算机协会，2002年)。已经进行了尝试来定制老化模拟，以便它更准确地描绘特定个人将来的老化的外观。例如，已经根据图像的人群队列结合与老化有关的有关面部变化的公布数据开发了老化算法以模拟个体的老化的外观(Hysert PE等人的“At Face Value”：ageprogression software provides personalized demonstration of theeffects of smoking on appearance”，Tobacco Control，第12卷，第238至240页，2003年)。这种方法的局限性在于，老化的图像是人群标准的反映，并且不一定反映个体独特的老化过程。

Boissiux等人开发了一种基于图像的模型以模拟皮肤老化，其中用预先计算好皱纹的一般蒙版作为在个人面部的3D模型上的纹理。利用了八个基本蒙版并且所用的具体蒙版与人的性别、面部形状和要模拟的表情种类相匹配(Boissiux等人的“Simulation of skin aging andwrinkle with cosmetic insight”，Computer Animation andSimulation，第15至27页，2000年)。因为它依赖于人群平均值，这种方法在其准确预报每个人独特的将随年龄显现的皮肤特征方面受到限制。

Zhang等人描述了一种用于将老年脸的几何细节转移到年轻脸的几何细节上以便使年轻脸看起来年老的方法(Zhang等人的“System andmethod for image-based surface detail transfer”，US7020347B2，2006年)。反过来，可将年轻脸的表面细节转移到年老的几何细节上使年老脸看起来年轻。这种方法受以下事实限制，即年老脸的老化特征将与年轻脸所将实际显示的特征不是完全一样。

发明概述

本发明涉及处理面部图像以检测和操纵皮肤特征例如色素沉着斑点、皱纹和细小纹理特征以便克服前述限制性的处理方法和装置。在本发明的一个方面，提供了可计算机执行的方法来检测和描绘在其中检测前述皮肤特征的数字面部图像的相关部分。此外，采用可计算机执行的方法来检测皮肤特征并用来操纵它们，例如通过强调或者和不再强调它们的外观，以模拟皮肤的老化和/或去老化。

在本发明的又一方面，处理在紫外照射下所采集的数字图像以检测在标准光照条件下不可见的斑点的存在以及预报它们的生长和潜在的可见性。

提供了用于辨别各种类型的面部特征(例如，斑点-皱纹-纹理-其它特征)并根据面部特征的类型正确地模拟面部特征的老化和去老化的方法。

根据下面的附图和详细说明，本发明的以上和其它方面及特征将显而易见。

附图概述

图1为显示根据本发明的用于面部皮肤的斑点、皱纹和纹理的老化、去老化模拟的示例性方法的高级流程图。

图2为显示根据本发明的示例性面部皮肤检测方法的流程图。

图3A显示基于整个面部斜视图图像生成的示例性面部皮肤蒙版；图3B显示示例性的斑点/皱纹老化蒙版(在黑线条内的区域)；以及图3C显示根据本发明的示例性实施方案所生成的示例性的纹理老化蒙版(在水平黑线下方以及到垂直黑线左侧的区域)。

图4为根据本发明的示例性实施方案的斑点老化模拟过程的流程图。

图5为根据本发明的检测紫外斑点并计算对比度的示例性过程的流程图。

图6为根据本发明的示例性斑点去老化过程的流程图。

图7A和7B显示根据本发明的示例性斑点检测算法的流程图。

图8为根据本发明的示例性皱纹老化和去老化模拟过程的流程图。

图9为根据本发明的示例性皱纹检测过程的流程图。

图10为根据本发明的示例性棱线检测过程的流程图。

图11为用于执行本发明的系统的示例性实施方案的方框图。

图12为根据本发明的示例性纹理老化过程的流程图。

图13为根据本发明的示例性纹理去老化过程的流程图。

图14为根据本发明的用于结合由斑点、皱纹和纹理所指示的面部皮肤老化模拟的示例性过程的流程图。

图15为根据本发明的用于结合由斑点、皱纹和纹理所指示的面部皮肤去老化模拟的示例性过程的流程图。

发明详述

示例性实施方案综述

图1为显示根据本发明的用于面部皮肤的斑点、皱纹和纹理的老化/去老化模拟的示例性方法的高级流程图。在101中，提供在标准光照下例如用常规的数字照相机采集的特写面部照片作为输入。在111中，在紫外光照模态下(在照相机前面的具有紫外滤光器的紫外光源)采集的同一对象的照片也被提供作为输入。为了提供标准化的和可重复产生的照明条件和图像基准，这两个图像优选地用自动控制的面部图像采集系统例如购自Canfield Scientific，Inc.的VISIA Complexion Analysis System(后文称为VISIA)进行采集。此外，这两个图像优选地应当由斜视图采集以更好地显示具有大皮肤斑点的脸颊区。

一般来讲，101中的标准光照图像输入将表示成RGB(红、绿、蓝)彩色图像。然而请注意，本发明并非仅限于任何特定格式。在步骤105中，将RGB图像转换成1976CIE L*a*b*颜色空间。这样一种颜色转换在本领域中通常用来分离图像的亮度和色度分量。后文中L*a*b*转换将称作LAB转换，并且转换的图像将称为LAB图像。LAB图像的L通道代表亮度，而A和B分量代表色度。本文所述的几种皮肤特征分析和再合成运算在LAB图像上执行。尽管所述的各种实施方案显示使用LAB颜色空间格式，但采用包括亮度和色度分量的其它颜色空间格式来执行本发明也是可能的。

在103中，执行面部皮肤检测，其要求确定来自整个面部图像的代表皮肤(而不是毛发、眼睛、嘴唇、鼻唇沟等)的那些像素。面部皮肤检测过程如下所述。

根据在103中确定的皮肤像素，操作然后前进到107，在此处描绘面部的具体区域或“蒙版”以执行斑点、皱纹和纹理的老化模拟。针对斑点和皱纹模拟生成覆盖面部的某些部分的第一蒙版，以及针对纹理模拟生成覆盖面部的某些部分的第二蒙版。下面详细描述了蒙版生成过程。

在113、115和117中分别进行斑点、皱纹和纹理的老化和去老化模拟。在113中的斑点老化/去老化模拟接收经LAB转换的标准图像(来自105)和在RGB域中的紫外图像(来自111)以及“斑点和皱纹老化蒙版”(来自107)，并在121中生成斑点老化的图像以及在122中生成斑点去老化的图像。

115中的皱纹老化/去老化模拟接收经LAB转换的图像(来自105)以及“斑点和皱纹老化蒙版”(来自107)，并分别在123和124中生成皱纹老化的图像和去老化的图像。

117中的纹理老化/去老化模拟接收经LAB转换的图像(来自105)以及纹理老化蒙版(来自107)，并分别在125和126中生成纹理老化的图像和去老化的图像。

下面非常详细地描述了斑点(113)、皱纹(115)和纹理(117)的老化和去老化模拟的具体实施以及其中将个体老化的图像和去老化的图像相结合的复合图像的生成。下面也描述了用于在计算机监视器上演示老化的图像和去老化的图像之间的过渡的一个交互式滑块应用软件。

图11是可用来执行本发明的一种系统1100的示例性实施方案的方框图。如图11所示，系统1100包括图像采集子系统1110例如前述的VISIA Complexion Analysis System等，其连接到通用计算机1120上，继而又连接到输出装置1130上。计算机1120可为被编程用于依照本发明进行操作的个人计算机等。输出装置1130可包括多种装置中的一种或多种，例如：常规计算机监视器等，计算机1120控制其显示图像例如依照本发明执行的各种模拟的结果；打印装置；存储装置；和通信装置等等。应当理解，可用各种各样的硬件配置实施本发明，并且不受图11系统的限制。

面部皮肤检测

基于皮肤特征的老化模拟应当在面部的皮肤区域上执行。在本发明的示例性实施方案中，面部的非皮肤区域例如嘴唇、毛发、眼睛、眉毛、鼻孔等被排除于模拟之外。由标准面部图像确定面部的皮肤区域。已经开发了几种皮肤检测算法用于各种用途，包括面部检测在内。(参见例如R.L.Hsu等人的“Face detection in color images”，IEEE Transactionson Pattern Analysis and Machine Intelligence，第24卷，第5期，第696至707页，2002年5月)。如果此类皮肤检测算法提供合适的颗粒性级别，则它们可用于根据本发明的面部皮肤老化模拟。

作为备选方案，可手工即用使用者输入执行皮肤检测(和随后的蒙版生成)。在给定一个面部图像的情况下，使用者可采用常规的基于计算机的绘图技术描绘面部的皮肤区域的轮廓。轮廓因此将限定老化/去老化模拟所要用的蒙版。尽管计算上简单，这种方法具有几个缺点。它具有在模拟过程中会包括面部的非皮肤部分的危险，并且也引入了人本身固有的主观性，可能导致结果上的巨大变化。

在一个优选的实施方案中，采用了新的皮肤检测算法，其根据斜视图或正视图仅截取面部皮肤的均匀亮度部分并排除非皮肤区域(眼睛、眉毛、毛发、髭和胡须)以及阴影皮肤区域(例如颈部区域)。皮肤检测根据从L、A和B测量值计算出的Individual Typology Angle(ITA)量度执行。(参见G.N.Stamatas等人的“Non-Invasive Measurements ofSkin Pigmentation In Situ”。Pigment Cell Research，第17卷，第618至626页，2004年)。ITA对于每个图像像素(i，j)被定义为arctan((L[i，j]-50)/B[i，j])并且其与皮肤上的黑素浓度有关。假设皮肤像素的ITA值将绕着某个值聚拢而非皮肤像素的ITA值明显远离皮肤像素的ITA值。

图2为显示说明根据本发明的示例性面部皮肤检测过程的流程图，其利用前述的ITA制。在皮肤检测之前，执行粗略的面部检测来从包含面部、毛发、颈部和背景的整个图像截取面部区域。所检测的面部区域应当包括面部的所有皮肤区域，但也可包括所有的面部特征(眼睛、眉毛、鼻孔、嘴唇、毛发)。为此，利用LUX颜色空间由特写图像分割面部区域。(参见M.Levin等人的“Nonlinear color space and spatiotemporalMRF for hierarchical segmentation of face features in video”，IEEE Transactions in Image Processing，第13卷，第1期，2004年1月)。

如图2所示，处理开始于标准的RGB全头部图像，例如以上在101中所提供的。在203中用以上在Levin参考文献中所描述的技术将所述图像由RGB转换到LUX空间。

在205中，分割面部区域。例如，这可通过在LUX图像的U通道上应用Otsu阈值法来进行。(参见N.Otsu的“A Threshold SelectionMethod from Gray-Level Histograms”。IEEE Transactions onSystems，Man，and Cybernetics，第9卷，第1期，第62至66页，1979年，后文“Otsu参考文献”)。在205中生成描绘面部区域的面部蒙版。面部皮肤检测过程的其余部分然后可仅在面部区域执行，从而减少搜索空间并降低计算成本。

在207中，根据在205中执行的分割所遮罩的原始RGB图像被转换到LAB空间。因此，在面部区域内执行随后的ITA制计算以进一步分割面部非皮肤部分。因为ITA制计算的除法和反切运算对噪声敏感，首先平滑L和B通道是优选的。如图所示，此类平滑可通过用二维Gaussian滤波器或其它类似技术过滤L和B图像而分别在209L和209B中进行。对于220PPI的工作分辨率，该滤波器的L通道方差被选为5，B通道方差被选为1.5。

在211中，根据以下公式在面部区域内计算每个像素的ITA：arctan((L[i，j]-50)/B[i，j])。ITA图像为处于[090]范围内的灰度图像，其中较小的ITA值与皮肤像素相对应以及较大的值与非皮肤像素相对应。将这个灰度图像在213中用Otsu阈值法分割成两个区域。为此，仅在面部区域中计算ITA图像的直方图。根据直方图，Otsu阈值算法返回一个阈值，其将把这个图像分割成具有最小类间差异的两类。此外，在这种阈值法中可加入与皮肤区域相对于整个面部图像的比率有关的先验信息。(参见Q.Hu等人的“Supervised range-constrainedthresholding”，IEEE Transactions in Image Processing，第15卷，第1期，第228至240页，2006年1月，后文称为“Hu参考文献”)。对于典型的斜视图图像，面部像素的至少25％应当属于皮肤像素。Hu参考文献描述了如何将该信息引入到基于Otsu的分割方法中。在由阈值法算法计算出最佳阈值之后，ITA值小于该阈值的像素被归类为皮肤像素。这以后，生成二元(黑和白)图像，其中用白色显示皮肤像素以及用黑色显示非皮肤像素。

在213中生成的分割皮肤区域可包括形成小岛的孤立的非皮肤像素。此类非皮肤岛屿可在215中通过采用盘结构元素的形态学闭运算或其它此类技术而消除。例如，对于220PPI的图像分辨率，选择这个盘的周长为10。作为另外一种选择，可具有在非皮肤面部特征(例如眉毛、毛发等)中所检测的皮肤碎片。这些小碎片也通过采用相同的盘结构元素用形态学开运算而消除。此外，由于独特的皮肤特征例如大色斑的缘故，一些个体可具有大的非皮肤碎片。这些也可通过应用形态学填充运算而消除。目标是要在一个包括脸颊、前额和鼻子在内但不包括鼻孔、有阴影的鼻唇沟、眼孔、眉毛和毛发(包括任何髭和胡须)的连续区域中检测面部皮肤。对于VISIA系统中采集的斜视图图像，图3A中显示了有效的面部皮肤蒙版的一个实例。这样一个面部皮肤蒙版对执行根据本发明的老化模拟是理想的。

老化模拟蒙版的设计

对每个皮肤特征(斑点、皱纹和纹理)的老化模拟可在与该特定老化模拟更为相关的面部皮肤区域的较小子集上执行。例如，对脸颊区域(在眼睛高度的下方以及嘴唇高度的上方)执行皱纹和斑点模拟比在其它面部皮肤区域这样做更加有效。为此，如图1所示，根据在103中生成的整个面部皮肤蒙版，在107中生成了两种不同的蒙版--斑点和皱纹老化蒙版和纹理老化蒙版。此类蒙版的实例示于图3B和3C中。这些蒙版根据眼睛、嘴唇和鼻子位置而设计。图3B中所示的皱纹和斑点包括从眼睛高度到嘴唇高度以及从鼻子高度到颊末的所有皮肤。图3C中所示的纹理蒙版可从眼睛高度向下延伸到颏的末端。眼睛和嘴唇区域被清晰地描绘在正面皮肤蒙版中。这些特征的位置可用该图像的垂直和水平投影进行计算。垂直投影的一个局部极小值提供眼睛的中心行，而第二局部极小值提供嘴唇的中心行。一旦确定了这些坐标，便据此修剪整个面部皮肤图像以生成两种前述的老化模拟蒙版。

斑点老化模拟

图4为根据本发明的示例性实施方案的斑点老化模拟过程的流程图。根据经L*a*b*转换的标准图像(来自105，图1)、紫外图像(来自111，图1)以及斑点和皱纹老化蒙版(来自107，图1)，该过程生成斑点老化的图像。如图4所示，提供了紫外图像以及如上所述所生成的斑点/皱纹蒙版作为斑点老化模拟过程的输入。采用荧光光谱技术采集的紫外图像对于色素沉着斑点表现出清晰可辨的标记(参见Miyamoto参考文献)。这种光照模态通常用于皮肤病学中以清楚地显示在标准图像中原本不可见的色素沉着损害。有力的证据表明，仅在紫外图像下可见的这些色素沉着斑点将由于光老化的缘故，随着色素沉着变得更加严重(即，随着黑素沉积增加)而变得可见。应当注意，尽管参照紫外图像和“紫外斑点”显示示例性实施方案，但紫外光并不是使亚皮肤表面特征能够可视化的唯一光谱。一般来讲，本发明的该方面适用于不能被肉眼容易看见的任何亚皮肤表面斑点，而与其中采集它们的光照模态的光谱无关。

图4中所示的本发明的示例性过程模拟上述过程。从紫外图像检测到的色素沉着斑点以及它们的对比度信息可用来调节标准图像中对应位置的强度和颜色对比，从而模拟“老化斑点”随时间的发展。

如上所述，标准和紫外图像理想地应当在模拟之前进行校准以便在显示器中最佳实现。例如用VISIA系统顺序地按最小延迟采集标准图像和紫外图像可减轻或消除对于校准的需要。然而，可采用几种熟知的校正技术中的任何一种校正未被正确校正的图像。(参见例如B.Srinivasa等人的“An FFT-Based Technique for Translation，Rotation and Scale-Invariant Image Registration”，IEEE Transactions on ImageProcessing，第5卷，第8期，1996年8月)。

假设图像被适当地校正，在403处执行基于紫外图像的紫外斑点检测。下面详细描述了一种根据本发明的示例性紫外斑点检测算法。紫外斑点检测算法返回紫外斑点的所有像素坐标以及它们的对比度信息。斑点被加上索引并且每个斑点附连一个具体标签(例如，号码)。加索引可通过逐行或逐列扫描代表紫外斑点的黑白图像并给每个斑点按次序赋予一个号码来进行。

在405处，紫外斑点抽取过程抽取一个斑点附近相邻的斑点以便并非紫外图像中的所有斑点在标准图像中均变成可见的。这种抽取过程可通过以下事实被证明是正确的，即仅所有紫外斑点的子集将变成在标准图像中可见。抽取过程可通过在带索引斑点列表中每隔一个或每隔两个选择一个斑点来进行。这将提供所有检测过的紫外斑点的稀疏子集。

在抽取之后，生成了剩余斑点的紫外对比度图像。紫外对比度图像是一种在紫外斑点的剩余子集中具有每个像素的紫外对比度强度的强度图像。在409处，剩余的紫外斑点的紫外对比度图像被扩大以放大紫外斑点。这对在标准和紫外图像中均可见的实际的色斑点将具有放大效应。紫外斑点的扩大可通过使紫外对比度图像变模糊来执行。这种操作可通过用一个二维Gaussian滤波器过滤紫外对比度图像来进行。Gaussian滤波器对于工作分辨率的方差设定为5并可被增大或减小以调节扩大效应。作为另外一种选择，可不进行扩大，因为模拟斑点老化而不扩大是可能的。

在411A、B和L处，采用扩大的紫外斑点对比度图像来修改原始标准图像亮度分量(L通道)和颜色分量(A和B通道)。紫外斑点对比度图像因此在被加到原始图像的L、A和B分量之前被加权。众所周知，色斑的效果随着强度变化而在L、A和B中是可见的。在一个示例性实施方案中，紫外斑点对比度在被加到L通道之前被乘以1.5(即，e_L＝1.5)，在被加到A通道之前乘以0.5以及在被加到B通道之前乘以0.5(即，e_A＝e_B＝-0.5)。这些数字的符号和绝对值根据研究发现和经验观察而确定。在将老化对比度相加之后，在413处通过执行LAB至RGB转换而合成斑点老化的图像。应当注意，如上所述，本发明并不仅限于任何特定的颜色或图像格式。例如，如果所得的图像要进行打印，则在413处的转换可为LAB至CMY转换(即，转换到熟知的典型用于打印的青-品红-黄颜色空间)。如所理解的那样，本发明所生成的图像可被打印、显示、存储、传输或经受任何进一步处理。此外，如所理解的那样，例如如果所得的图像要以LAB格式进行存储或传输，则在413处的转换可被免除或延缓。

为此，通过智能地将紫外对比度信息的因子加到强度(L)和颜色(A和B)分量中而在LAB色域中进行老化模拟。通常用色度计在LAB域中广泛研究和量化色素沉着过度(S.Alaluf等人的“The impact ofepidermal melanin on objective measurements of human skincolour”，Pigment Cell Research，第15卷，第119至126页，2002年，后文称为“Alaluf参考文献”)并在LAB域中分析图像(N.Kollias等人的“Optical Non-invasive Approaches to Diagnoses ofSkin Diseases”，Journal of Investigative DermatologyProceedings，第7卷，第1期，第64至75页，2002年)。一种研究涉及用LAB色度计对人皮肤的正常区域和色素沉着区域进行颜色测量，该研究显示所有的L、A和B值均随色素沉着程度(黑素含量)而变化。(参见Alaluf参考文献)。据报告，L值将随着黑素含量增加而变小，而A和B值将随黑素含量增加而变大。这解释了色素沉着斑点的暗棕色外观。

紫外斑点检测和对比度计算

图5为检测紫外斑点和计算对比度的示例性过程的流程图。该过程取前述紫外图像的蓝色通道(例如来自111，图1)，并返回紫外斑点以及紫外对比度图像。紫外图像的蓝色通道在各通道(R、G和B)中表现出最好的对比度，因为紫外荧光在蓝色光谱中较强。图5过程的目标是要从这个灰度图像中提取紫外斑点损害。

在503处，使蓝色通道紫外图像经受噪声过滤，其中图像上的小变化得到平滑。为此，已经发现[5x5]中值滤波器对于紫外图像(具有220PPI的工作分辨率)很有效果。因为光源的不均匀强度场和面部的三维形状，不是所有的图像像素均接收到等量的光，因此促使图像在面部的不同区域具有变化的强度。这种强度上的变化使得无法使用固定阈值来分离视觉上比背景更暗的紫外斑点损害。为补偿不均匀强度，在505中估算一种缓变的背景强度并从每个像素的已过滤的强度上移除。缓变的背景强度可通过利用具有用大滤波器支持的低通滤波器进行估算。这样一种滤波器可采用Wiener滤波器，即根据局部平均和局部方差估计低频二维强度表面的自适应低通滤波器，来执行。可用于该目的的一种示例性Wiener滤波器通过一组图像像素更新公式描述于附录A-1中。例如将Wiener滤波器的支持(尺寸)选为[41x41]，大得足以包封平均大尺寸的紫外斑点，假设为220PPI的工作分辨率。

当将背景强度水平从原始强度图像的噪声过滤形式移除时，便获得对比度图像。这个对比度图像包括正负两个分量。紫外斑点位于负对比度区域的子集中。因此，在507处，紫外斑点通过用固定阈值分割负对比度图像而获得。在一个示例性实施方案中，选择这个阈值处在约-3.5至-5.0的范围内。紫外斑点的判据是它的对比度值应当小于这个阈值。这种斑点分割与平均的人类感知吻合良好。

作为507中的分割运算的结果，获得了二元(黑色和白色)图像，其中白色损害代表紫外斑点以及黑像素代表背景。在509中平滑这个图像，例如用[5x5]中值滤波器。在511中，把紫外斑点加上索引并标号，并且计算与每个紫外斑点相关的面积(例如，像素的数目)。在513中，消除面积小于阈值(例如，150个像素)的小紫外斑点和面积大于阈值(例如，600个像素)的大紫外斑点。对于每个像素，即在过程506中生成的紫外斑点图像，返回剩下的紫外斑点以及对比度值。重要的是应记得这些对比度值是负值并代表深对比度。任选地，根据紫外对比度图像(ID)通过在515中进行对比度加权的评分而生成严重性评分。这个评分通过相加在有效紫外斑点内部的所有ID值进行计算。这个评分与色素沉着程度有关并可用来监测色素沉着的恶化和改善。此外，在517中计算所检测的紫外斑点周长以便可将它们重叠在紫外图像上来显示紫外斑点。

斑点去老化模拟

在一个示例性实施方案中，利用标准图像的L、A和B通道在LAB颜色空间中执行斑点去老化模拟。连同色素沉着的斑点一起，在这些通道中可辨别出红色斑点(因结疤和皮肤病例如痤疮而引起的小红肿区域)。为使模拟更加逼真，优选地根据本发明的示例性实施方案去除此类皮肤特征。

图6为根据本发明的一个示例性斑点去老化过程的流程图。采用LAB图像(例如来自105，图1)以及斑点和皱纹蒙版(例如来自107，图1)，在603中执行斑点检测和对比度计算。下面参照图7A和7B详细描述了一种示例性的斑点检测和对比度计算算法。对比度是指相对于由局部邻域计算出的低通背景强度的像素强度差异。

在609L中将斑点损害内的L上的对比度值乘以值e_L并在611中将它加到原始L通道上来使L的负对比度与背景水平齐平。类似地，在609A处将在斑点损害内A中的对比度值乘以值e_A并在611A处将它加到原始A通道上来使A的色差与背景颜色齐平。类似地，在609B中将斑点损害内B中的对比度值乘以值e_B并在611B中将它加到原始B通道上来使B的色差与背景颜色齐平。应当注意，L中的对比度用来修改斑点的暗度而A和B中的对比度用来修改斑点的颜色。在斑点损害内去除L、A和B通道中的对比度将会使斑点损害的强度和颜色与背景皮肤的强度和颜色齐平。这将具有移除斑点损害和面部皮肤外观更光滑的视觉效果。因此，斑点去老化的图像可用来预报有效治疗所期待的结果。

斑点检测算法

图7A和7B示出了一种示例性斑点检测算法。在703中，将标准RGB图像701转换到LAB颜色空间，并在705L、705A和705B中分别对L、A和B通道单独地应用噪声过滤。在所示的示例性实施方案中，用具有较小滤波器支持例如[5x5]的Wiener滤波器执行噪声过滤，如上所述。然后，在707L、707A和707B中，将具有例如[61x61]支持的Wiener滤波器分别应用到经噪声滤波的L、A和B图像上来估算在斑点和皱纹模拟蒙版中每个像素的背景强度和颜色。对于L、A和B通道的每个像素，通过分别在708L、708A和708B中用每个像素的经噪声滤波的L、A和B值减去每个像素的低通L、A和B值来计算对比度值。

对比度图像(此处“对比度图像”是指所有以对比度值作为强度的图像像素的集合)为良好的斑点指示物。在皮肤病学研究中已得到确认，斑点损害的强度小于背景皮肤的强度，并且它们在A和B中的颜色分量大于背景皮肤的颜色读数。(应当注意，背景皮肤在本文中被认为是健康的和光滑的，以及斑点损害被认为在背景里面是稀疏的)。根据这些判据，分别在709L、709A和709B中在通道L中的负对比度区域和通道A和B中的正对比度区域中选择斑点损害。此外，由708L、708A和708B获得的对比度图像通过709L、709A和709B运算被细化以产生更有意义的对比度图像。在这些运算后对比度图像被用于斑点的去老化模拟。

在711中，根据来自L、A和B通道的对比度值对于每个像素计算斑点色差制(DE)。CIE L*a*b*感知的色差制常常用于颜色学以量化人的视觉对两种色片间差异的敏感性。在示例性的实施方案中，采用这个尺度将斑点颜色与背景皮肤颜色区分开来以便斑点分割与人的感知相符。一般来讲，如果这个尺度大于3.5，则双眼可辨别颜色上的区别。

前进到图7B，在713中通过比较DE与阈值例如4.5以分割出斑点来分割斑点。根据所期望的敏感度，这个阈值可从3.5改变到5。在这个阈值运算之后，获得了斑点损害的二元(黑色和白色)图像，其中白色小岛代表斑点。在715中任选地平滑这个二元图像以具有平滑形状的斑点损害。

在717中，通过分配编号来对分割对象进行标记。

基于例如如上所述的阈值DE的斑点分割步骤一般将分割出皱纹部分和大毛孔的子群以及斑点。在719中，小对象例如毛孔(一般小于斑点)通过将最小面积约束应用到分割对象上而消除。例如，在给定分辨率(220PPI)下100个像素的面积阈值是令人满意的。

在721中，为了消除皱纹和皱纹状特征，然后计算剩余斑点损害的某些形状性质。示例性的性质可包括面积、纵横比、实心度、长轴长度、短轴长度、偏心度和范围。这些均是本领域所常用的二维形状性质并在附录A-3：Definitions of Shape Properties中加以定义。为消除皱纹和皱纹状特征，用纵横比(短轴长度/长轴长度)作为判据。例如，具有小于0.25的纵横比的对象可视为皱纹并作为斑点予以消除。例如同样可采用0.3的范围阈值来消除变形的和模糊的形状特征。(范围是密实度测量值，在[01]范围内变化，其中高值对应致密对象)。在721中应用这些形状和尺寸约束之后，记录剩下的对象和其像素位置以及前面在709L、709A和709B中对于这些斑点位置计算的对比度值。对比度值被用于斑点的去老化模拟。任选地，在723中根据在711中计算的整个对比度图像(DE)生成严重性评分。该评分通过相加有效斑点内的所有DE值进行计算。该评分与色素沉着过度和皮肤不均匀性的程度有关并且可用来监测皮肤状况的恶化和改善。此外，在725中计算所检测的斑点周长以便可将它们覆盖在图像上来显示斑点。

皱纹老化和去老化模拟

图8是根据本发明的一个示例性皱纹老化和去老化模拟过程的流程图。在801中，采用由以上生成的斑点和皱纹蒙版所遮罩的LAB图像的亮度(L)通道执行皱纹检测步骤。皱纹特征的颜色分析证明，细皱纹中的颜色并不明显不同于背景皮肤的颜色。然而皱纹的亮度明显不同于背景亮度。为此，采用L通道来检测和模拟皱纹老化/去老化。801中的皱纹检测步骤提供皱纹特征以及它们的“皱纹强度”值。“皱纹强度”是一个与对比度不同的量度并根据定向滤波器进行计算。(参见W.T.Freeman等人的“The design and use of steerable filters”，IEEE Trans.Pattern Analysis and Machine Intelligence，第13卷，第9期，第891至906页，1991年，后文称为“Freeman参考文献”)。下面详细描述了一种示例性的皱纹检测算法。

在皱纹检测之后，在803中进行伪皱纹消除步骤。801中的皱纹检测步骤所生成的候选皱纹被分割出，作为在暗皮肤背景上的白色对象。将这个黑白图像称作棱线对象图像。大多数棱线对象是由于皱纹和细纹的缘故，但有一些可来自其它面部特征例如大斑点的边界、对齐的毛孔、皮肤上的黑毛发和细血管丝等等。这些伪特征中的大部分可根据一组形状、尺寸和颜色判据进行消除。下面更详细描述了这种示例性过程。该过程返回有效皱纹以及它们的强度图像。皱纹强度图像在有效皱纹像素上取棱线图谱的值并在别处取零。该皱纹强度图像将被用于皱纹老化和去老化模拟。

对于皱纹的老化模拟，在807中扩大皱纹强度(后文称为皱纹对比度)以得到增浓效果，其将伴随老化一直发生。例如，可用例如滤波器方差为2的二维Gaussian滤波器执行扩大操作。上面相对于紫外斑点扩大描述了这个步骤。在809A中，将扩大的对比度接着乘以增强因子e_L(例如，2)并加到L通道上。这些运算的净效应是在原始图像中所见的皱纹显得更暗和更密，并且在原始图像中不是清楚可见的淡皱纹变得可见。最后，在811中通过LAB至RGB转换接着合成了皱纹老化的图像。重要的是要注意到，皱纹将随年龄生长，并且这个过程的模拟可通过延伸所检测的皱纹而进行。皱纹的延伸可除了扩大运算之外或作为其另外一种选择来进行。

对于皱纹的去老化，不执行对比度扩大并且在809D处从L通道移除皱纹对比度以便将皱纹的强度水平带到周围背景皮肤的强度水平上。最后，在813中通过LAB至RGB转换合成了皱纹去老化的图像。

皱纹检测算法

现在将参照图9描述皱纹检测过程。在901中，转换被斑点和皱纹蒙版所遮罩的标准RGB图像以获得LAB图像。在示例性的实施方案中，仅L通道用来检测皱纹。在903中，在皱纹老化模拟蒙版内将采用如上所述的Wiener滤波器的噪声过滤过程应用于L通道。滤波器具有例如[3x3]的支持。在905中，将具有例如[21x21]的支持的又一个Wiener滤波器应用到噪声过滤过的L通道上以估算背景照明强度。优选地，这个滤波器的规格应当足够大以覆盖工作分辨率中的皱纹。在907中，通过在皱纹蒙版内用噪声过滤过的L值减去低通L值计算每个像素的对比度值。

在909中，选择具有负对比度值的区域(例如，暗区域)用于皱纹检测。这是因为根据观察，细皱纹比背景要暗(在L中较低)。在911中，将棱线检测步骤应用到负对比度图像上来检测细长结构。在下面更详细描述的一个示例性实施方案中，棱线检测步骤采用定向滤波器(参见Freeman参考文献以及J.Staal等人的“Ridge-based segmentation incolor images of retina”，IEEE Transaction on Medical Imaging，第23卷，第4期，第501至509页，2004年4月，后文称为“Staal参考文献”)。棱线检测步骤接受对比度图像并返回“棱线强度”和“棱线取向”图像。这两个图像在913中进行进一步处理以获得经过修改的棱线强度图像或“棱线图谱”。下面详细描述了棱线图谱计算步骤。棱线图谱为代表曲线结构的灰色亮度图像并表现出对皱纹的强烈反应。

为从棱线图谱图像中确定皱纹结构，在915中应用滞后阈限(参见F.J.Canny的“A computational approach to edge detection”，IEEETransactions on Pattern Analysis and Machine Intelligence，第8卷，第6期，第679至698页，1986年)。滞后阈值法是将弱结构连接到强结构上的更软形式的阈值法并涉及低和高阈值。这些阈值的示例性值为4和8。

棱线检测

皱纹在标准图像中将本身显示成细长的结构。它们在强度上(L通道)相对于背景皮肤强度水平大部分是可见的并且与背景皮肤颜色相比在颜色方面(A和B通道)几乎没有差别。因此，可通过利用为细长结构所设计的检测器从L通道中将它们提取出来。

高斯核的二阶方向导数通常用来在图像处理中检测细长结构。(参见例如Staal参考文献。这些导数实际上为Freeman参考文献所描述的一类导向滤波器。)这些滤波器是对棱线特征敏感的基本滤波器并具有垂直、水平和对角取向。图10是利用此类导向滤波器的示例性棱线检测步骤的流程图。

如图10所示，分别在1001A、1001B和1001C中用例如以上所述的第一、第二和第三滤波器使DE图像经受二维卷积。为了分析一个图像中结构的取向和强度，在1003中对于每个像素形成了Hessian矩阵，其元素为基础滤波器响应。接着，在1005中，对每个像素的[2x2]Hessian矩阵执行Eigen解析。Eigen解析返回两个特征值(e1，e2)和彼此正交的两个特征向量(v1，v2)。棱线强度被定义为正特征值(例如，e1)-如果它的绝对值大于第二特征值(例如，e2)，并且棱线取向为第二特征向量(v2)。在这种方法的一个变型中，当e1＞0和|e1|＞|e2|时棱线强度被定义为(e1-e2)。已经观察到这种定义更好地强调皱纹结构。

棱线图谱计算

如上所述，示例性的棱线检测过程给每个像素返回两个有用的参数：棱线强度，为指示皱纹多深的标量值；和棱线取向，为指定皱纹在特定像素位置的方向的向量。根据这两个参数生成棱线图谱图像。在这种情况下，对于每个像素，定义新的棱线强度，其将原始棱线强度和取决于相邻像素取向的强度项考虑在内。这个强度项通过将当前像素的方向向量与当前像素的8个相连的邻域中的每个像素的方向向量的内积相加进行计算。这个过程用附录A-2中的一组公式进行描述。

伪皱纹消除

伪皱纹消除过程的目标是根据形状和尺寸性质消除伪正值(伪皱纹)。为此，把滞后阈限(915)之后的所有候选皱纹标号并且对于每个计算多个形状性质。这些形状性质可包括：短轴长度、长轴长度、面积、实心度和偏心度。这些二维形状性质的定义是标准的并在附录A-3中给出。

根据这些形状性质，将棱线对象分成四类：短皱纹、长皱纹、网络皱纹和非皱纹。要落在前三类之一中，棱线对象的性质必须满足对应的判据组。例如，对于是短皱纹的棱线形状，其长度必须介于最小(例如，30个像素)和最大长度(例如，50个像素)阈值之间；其纵横比(短轴长度/长轴长度)必须小于纵横阈值(例如，0.25)；其偏心度必须大于偏心度阈值(例如，0.97)，以及其实心度必须大于最小实心度阈值(例如，0.49)。类似地，对于长皱纹有一组判据以及对于网络皱纹有另一组判据。这些阈值根据检查一组训练图像上的皱纹而凭经验确定。未归类为这些皱纹类型中的一种的棱线对象被归类成非皱纹。剩余的棱线对象被称作有效皱纹并返回到皱纹检测算法。

纹理老化模拟

本文所用的术语“纹理”是指干扰皮肤总体光滑度的小皮肤特征。纹理老化和去老化模拟是根据纹理特征和对比度的检测。纹理特征包括毛孔、小白色斑和小的粗糙破坏。在纹理蒙版内执行纹理老化和去老化模拟。图3C示出了典型的纹理蒙版。

图12为根据本发明的一个示例性纹理老化过程的流程图。采用被纹理蒙版(例如，来自107的纹理蒙版，图1)所遮罩的标准面部图像(例如，来自105的标准面部图像，图1)的亮度(L)通道执行纹理老化模拟。在1201中，从L通道去除低通背景强度。为此，通过应用具有例如[21x21]的滤波器支持的例如如上所述的Wiener滤波器计算背景强度水平。从L通道减去该项以生成对比度图像。对比度图像具有正负两分量。具有负对比度值的区域被称作低纹理区域以及具有高对比度值的区域被称作高纹理区域。低纹理区域的实例为毛孔，而高纹理区域的实例为很小的白斑点。

在1203L中通过用负阈值(例如，-2.5)对对比度图像进行阈限，即通过选择对比度小于这个阈值的像素，进行低纹理区域的分割。此外，将分割纹理损害标号并且也记录下损害面积。应用小面积阈值(例如，10)来去除主要由于噪声产生的很小损害。应用大面积阈值(例如，120)来去除由于小斑点和皱纹的缘故产生的大损害。

在1205L处记录下每个像素的剩余纹理损害和它们的对比度值(低纹理对比度图像)。此外，在1207处通过采用具有方差1的二维Gaussian滤波器，扩大低纹理对比度图像。这种扩大运算的净效应是放大面部图像中的毛孔。毛孔放大随着年龄或者皮肤健康的恶化而自然地发生。可增加方差值来增加放大程度。

类似地，为分割高纹理区域，在1203H中将正阈值(一般为2.5)应用到正对比度图像上，即通过取大于这个阈值的像素。接下来，将分割纹理损害标号并记录下损害面积。应用小面积区域(一般10个像素)以去除主要因噪声的缘故产生的很小损害。应用大面积阈值(一般100个像素)来去除由于日照即过量的光照射在脸上产生的大损害。在1205H中记录下每个像素的剩余的纹理损害和它们的对比度值。

在1209L中，将扩大的低纹理对比度乘以增强因子e₁并在1211中加到L通道上。类似地，将高纹理对比度图像乘以增强因子e_h并在1211中加到L通道上。增强因子e₁和e_h的示例性值分别是1.0和0.5。在1213中，通过LAB至RGB转换合成纹理老化的图像。

纹理去老化模拟

示例性纹理去老化模拟旨在降低纹理特征例如毛孔和小白斑点的尺寸和强度。在斑点或皱纹去老化模拟中完全去除纹理特征将引起过于光滑的外观并且不会提供逼真的皮肤图像。

图13是根据本发明的示例性纹理去老化过程的流程图。也在亮度通道中执行纹理去老化。在1301中，从L通道移除低通背景强度。为此，通过应用具有例如[21x21]的滤波器支持的如上所述的Wiener滤波器计算背景强度水平。从L通道减去背景强度水平以生成对比度图像。对比度图像具有负和正两分量。具有负对比度值的区域被称作低纹理区域以及具有高对比度值的区域被称作高纹理区域。

在1303L处，为分割低纹理区域(即，大毛孔)，将负阈值(为一般-2.5)应用到负对比度图像。分割纹理损害被加上索引并标号并且也记录下损害的面积。此外，应用小面积阈值(一般为50)来去除小毛孔以及应用大面积阈值(一般为120)来去除由于斑点和皱纹的缘故产生的大损害。在1305L中计算并记录下每个像素的剩余纹理损害(其大多数为大毛孔)以及它们的对比度值(低纹理对比度图像)。在1307L中，通过对具有例如周长2的盘结构元件的低纹理对比度图像应用形态学扩大运算而使低纹理区域经受收缩。这种运算的净效应是毛孔的收缩以及在面部皮肤上毛孔的颜色变深减轻，与在有效治疗后改善皮肤状况有关。

类似地，在1303H处，为分割高纹理区域(很小的白斑点)，通过取大于某个阈值(例如，2.5)的像素将正阈值应用到正对比度图像。分割纹理损害被标号并也记录下损害的面积。应用小面积阈值(例如，30个像素)来去除很小的损害以及应用大面积阈值(例如，300个像素)来去除因收缩的缘故产生的大损害。在1305H处计算并记录下每个像素的剩余纹理损害以及它们的对比度值(高纹理对比度图像)。在1307H处，通过对具有例如周长2的盘结构元件的高纹理对比度图像应用形态侵蚀操作而使高纹理区域经受收缩。这种操作的净效应是明显的大白斑点毛孔的收缩以及在面部皮肤上这些特征的强度变小，同样与在有效治疗后改善皮肤状况有关。

在1309L处，将降低的低纹理对比度乘以增强因子e₁并在1211处加到L通道上。类似地，将高纹理对比度图像乘以增强因子e_h并在1311中加到L通道上。增强因子e₁和e_h的示例性值分别为1.0和1.0。在1313中，通过LAB至RGB转换合成纹理老化的图像。

总体皮肤老化和去老化模拟

可将上述对面部皮肤由于斑点、皱纹和纹理而出现的老化的模拟相组合来模拟面部皮肤的总体老化。图14为用于这样做的示例性过程的流程图。通过用斑点、皱纹和纹理的老化对比度修改L、A和B通道，在LAB色域中合成总体老化的图像。如图14所示，为了将斑点、皱纹和纹理合并进整体过程，在1401S中生成在L、A和B通道中的老化对比度图像；在1401W中生成在L通道中的皱纹老化对比度图像；以及在1401T中生成L通道中的纹理老化对比度图像。分别在1403SA、1403SB和1403SL中通过w_s因子加权每个A、B和通道斑点老化图像；在1403WL中用因子w_w加权L通道皱纹老化图像；以及在1403TL中用因子w_t加权L通道纹理老化图像。选择三个加权因子w_s、w_w和w_t来强调或不再强调总体老化的图像的各个分量的贡献。分别在1405SA和1405SB中将所加权过的A和B通道斑点老化的图像加到最终图像的A和B图像上。在1405L中将所加权过的L通道斑点、皱纹和纹理图像相组合并在1407L中加到最终图像的L通道上。在1409中使如此修改的L、A和B通道经受LAB至RGB转换以生成处在RGB域中的总体老化的图像。

以类似的方式，可将上述对由于斑点、皱纹和纹理而出现的面部皮肤去老化的模拟相组合以模拟面部皮肤的总体去老化。图15为用于这样做的示例性过程的流程图。分别在1501SL、1501SA、1501SB以及1501W和1501T处通过各自的斑点、皱纹和纹理去老化模拟生成在L、A和B中的去老化对比度。在1503SL、1503W和1503T中用各自的权重因子w_s、w_w和w_t加权L中的每个对比度图像来强调或不再强调各自分量对总体去老化的图像的贡献。这些权重因子的优选值全都是1。在1505中结合所加权过的L通道斑点、皱纹和纹理对比度图像并在1507中加到最终图像的L通道上。类似地，在1503SA处对A中的斑点对比度用w_s进行加权并在1507SA中加到A通道上。在1503SB处对B中的斑点对比度用w_s进行加权并在1507SB加到B通道上以得到最终的A和B通道。在1509处进行的LAB至RGB转换生成处于RGB域中的总体去老化的图像。在最终的图像中，消除了突出的皮肤特征并减少了小皮肤特征(毛孔)。这样一种图像对预报在对色素沉着过度、皱纹或皮肤纹理采用治疗后受试者皮肤面部可能看起来如何非常有用。

用于皮肤老化/去老化模拟的交互工具

可在计算机监视器上通过并列显示原始图像和模拟的图像并提供交互式滑块控制以使观察者能够调节老化程度来证明根据本发明的一个示例性实施方案的皮肤老化/去老化模拟。取决于所期望的模拟(斑点、皱纹、纹理或它们的任何组合)，将老化的或去老化的图像与原始图像相混合，其中混合程度取决于滑块位置。当滑块处在中间位置上时，原始图像被显示在左和右两个面板中。当使用者上移滑块时，通过将原始图像与去老化的图像进行α混合，去老化模拟图像被显示在右面板上。类似地，当使用者下移滑块时，通过将原始图像与老化的图像进行α混合，显示老化模拟图像。α混合是两个图像的线性加权并且是本领域常用来混合两个图像的标准运算。对于本专利申请而言，可用α混合离线生成斑点、皱纹和纹理的各种老化的和去老化的图像并且优选地实时执行图像渲染。

应当注意，在上述的每个老化和去老化模拟中，要进行模拟的老化或去老化的程度优选地在合适的期限内是使用者可选择的，例如，5至10年，以证明自然老化，或者几个月，以证明因治疗产生的去老化。

应当理解，上述实施方案示出了仅几个可代表本发明应用的具体实施方案。本领域的技术人员可作出很多不同的其它布置而不脱离本发明的实质和范围。

附录

A-1.WIENER滤波器

给定一个[MxN]灰度图像g，通过(i，j)给出在坐标g(i，j)处的值，下面的步骤采用中心在(i，j)处的局部[KxK]分析窗口执行Wiener滤波器，其中K是奇数。

1.在位于像素值所在之处的当前像素的[KxK]邻域中计算局部均差μ(i，j)和局部方差σ²(i，j)：(i，j)g(i，j)

L＝(K-1)/2

$\mu (i,j)={\Sigma }_{m=-L}^L{\Sigma }_{n=-L}^{L}g(i+m,j+n)$

${\sigma }^2(i,j)=-\mu {(i,j)}^2+{\Sigma }_{m=-L}^L{\Sigma }_{n=-L}^L{g}^2(i+m,j+n)$

2.通过平均横跨整个图像的局部方差计算噪声方差σ_w²

3.用下面的更新公式计算经过滤波的图像像素值f(i，j)：如果σ²(i，j)＞σ_w²

$f(i,j)=\mu (i,j)+\frac{({\sigma }^2(i,j)-{{\sigma }_w}^2)}{{\sigma }^2(i,j)}(g(i,j)-\mu (i,j))$

否则

f(i，j)＝μ(i，j)

4.对于图像中的所有像素重复步骤3。

A-2.棱线图谱生成

对于处在感兴趣区域(ROI)中的每个像素坐标(i，j)执行下面的计算步骤：

1.由棱线检测器获得下面的数量：

R(i，j)：棱线强度，正实数

V(i，j)：棱线取向向量，具有实数的2元向量。

2.根据这些数量，按处于八相连的邻域中的定向向量的内积之和计算方向强度：

$\mathrm{Ds}(i,j)={\Sigma }_{n=1}^8<V_c,V_n>,$其中<.>表示内积运算，以及V_c表示当前像素的棱线定向向量，以及V_n表示邻域中第n个像素的棱线定向向量。

3.将一部分方向强度加到棱线强度上来计算棱线图谱：

Rm(i，j)＝R(i，j)+αDs(i，j)其中α为处在[0.20.5]范围内的权重因子。

A-3.形状性质的定义

性质            定义

面积            对象中的像素数目

长轴长度        与对象具有同一正规化二阶中心矩的椭圆的长轴长度(单位                    为像素)。

短轴长度        与对象具有同一正规化二阶中心矩的椭圆的短轴长度(单位                    为像素)。

范围            也处在对象中的像素在边框中的比例。边框是包含对象的最                    小的矩形。

偏心度          与对象具有同一二阶矩的椭圆的偏心度。偏心度是在椭圆焦                    点之间的距离和它的长轴长度的比率。

实心度    也处在对象中的像素在凸壳中的比例。凸壳是包含对象的最              小凸多边形。