CN117368122A - 一种基于比色卡的frd宫颈染色实时比对方法 - Google Patents

一种基于比色卡的frd宫颈染色实时比对方法 Download PDF

Info

Publication number
CN117368122A
CN117368122A CN202311666024.6A CN202311666024A CN117368122A CN 117368122 A CN117368122 A CN 117368122A CN 202311666024 A CN202311666024 A CN 202311666024A CN 117368122 A CN117368122 A CN 117368122A
Authority
CN
China
Prior art keywords
sliver
rgb image
image
cervical
absorbent
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
Application number
CN202311666024.6A
Other languages
English (en)
Other versions
CN117368122B (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.)
Jintai Tianjin Medical Instrument Co ltd
Original Assignee
Jintai Tianjin Medical Instrument Co ltd
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 Jintai Tianjin Medical Instrument Co ltd filed Critical Jintai Tianjin Medical Instrument Co ltd
Priority to CN202311666024.6A priority Critical patent/CN117368122B/zh
Publication of CN117368122A publication Critical patent/CN117368122A/zh
Application granted granted Critical
Publication of CN117368122B publication Critical patent/CN117368122B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B10/00—Other methods or instruments for diagnosis, e.g. instruments for taking a cell sample, for biopsy, for vaccination diagnosis; Sex determination; Ovulation-period determination; Throat striking implements
    • A61B10/02—Instruments for taking cell samples or for biopsy
    • A61B10/0291—Instruments for taking cell samples or for biopsy for uterus
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00—Sampling; Preparing specimens for investigation
    • G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • G01N1/30—Staining; Impregnating ; Fixation; Dehydration; Multistep processes for preparing samples of tissue, cell or nucleic acid material and the like for analysis
    • G01N1/31—Apparatus therefor
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/75—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated
    • G01N21/77—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator
    • G01N21/78—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator producing a change of colour
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/40—Extraction of image or video features
    • G06V10/52—Scale-space analysis, e.g. wavelet analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/40—Extraction of image or video features
    • G06V10/56—Extraction of image or video features relating to colour
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B10/00—Other methods or instruments for diagnosis, e.g. instruments for taking a cell sample, for biopsy, for vaccination diagnosis; Sex determination; Ovulation-period determination; Throat striking implements
    • A61B10/02—Instruments for taking cell samples or for biopsy
    • A61B2010/0216—Sampling brushes

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • Molecular Biology (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Biomedical Technology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Gynecology & Obstetrics (AREA)
  • Reproductive Health (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Plasma & Fusion (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Image Analysis (AREA)

Abstract

本发明涉及图像处理技术领域，具体涉及一种基于比色卡的FRD宫颈染色实时比对方法，包括：通过可视宫颈采样刷获取吸水棉条RGB图像；根据吸水棉条RGB图像获取吸水棉条光照图像和吸水棉条反射率图像；根据吸水棉条RGB图像中任意两个像素点的颜色相似度，获取吸水棉条RGB图像的每个窗口的高斯环绕函数的尺度参数，进而得到每个窗口的高斯环绕函数；将不受光照影响的吸水棉条图像与比色卡进行比对得到测试结果。本发明提高吸水棉条判断结果的准确性。

Description

一种基于比色卡的FRD宫颈染色实时比对方法
技术领域
本发明涉及图像处理技术领域，具体涉及一种基于比色卡的FRD宫颈染色实时比对方法。
背景技术
宫颈染色实时比对可以帮助早期检测宫颈疾病，如宫颈癌前病变。颜色变化可能是疾病发展的早期指标之一，通过比对颜色，可以及早发现潜在的问题。在实施实时比对之前，需要对比色卡进行标定。这意味着将比色卡放在相同的光线和环境条件下，使用专业的仪器或相机来拍摄比色卡，以建立比色卡颜色的准确参考。基于比对的结果，可以为宫颈样本提供评估或分类。这可能涉及到判断样本是否正常或存在异常，并采取进一步的医学行动。这个过程需要专业的设备和培训有关的医疗专业人员，以确保准确性和可靠性。同时，为了保持一致性，标定比色卡以及在分析宫颈样本图像时保持一致的光照和环境条件都是非常重要的，因此对于拍摄环境的要求比较高，不太好将此类产品普及开来。
通过Retinex算法可以解决上述此类产品不好普及的问题，Retinex算法主要步骤是需要对原图像进行模糊化处理，但是模糊化处理时使用的高斯环绕函数中的尺度参数的选择尤为重要，因此如何根据图像中的特征量化得到合适的尺度参数是需要解决的重要问题。
发明内容
为了解决上述问题，本发明提供一种基于比色卡的FRD宫颈染色实时比对方法，所述方法包括：
通过可视宫颈采样刷获取吸水棉条RGB图像；
根据吸水棉条RGB图像获取吸水棉条光照图像；
获取吸水棉条RGB图像中任意两个像素点的RGB颜色通道平均差异因子，根据RGB颜色通道平均差异因子获取吸水棉条RGB图像中任意两个像素点的颜色相似度；获取吸水棉条RGB图像的若干个窗口；根据吸水棉条RGB图像中每个窗口中任意两个像素点的颜色相似度，获取每个窗口的高斯环绕函数的尺度参数；根据每个窗口的高斯环绕函数的尺度参数，获取每个窗口的高斯环绕函数；
根据吸水棉条光照图像和每个窗口的高斯环绕函数，获取不受光照影响的吸水棉条图像。
优选的，所述根据吸水棉条RGB图像获取吸水棉条光照图像，包括的具体方法为：
根据RetinexL算法将吸水棉条RGB图像分解为吸水棉条光照图像和吸水棉条反射率图像；且三张图像存在的关系表达式为：
式中，表示吸水棉条RGB图像中坐标为/>的像素点灰度值；/>表示吸水棉条光照图像中坐标为/>的像素点灰度值；/>表示吸水棉条反射率图像中坐标为/>的像素点灰度值。
优选的，所述吸水棉条RGB图像中任意两个像素点的颜色相似度的具体公式为：
式中，表示吸水棉条RGB图像中第/>个像素点和第/>个像素点的颜色相似度；/>表示吸水棉条RGB图像中第/>个像素点和第/>个像素点的第一颜色通道差值；/>表示吸水棉条RGB图像中所有像素点中最大的第一颜色通道值；/>表示吸水棉条RGB图像中所有像素点中最小的第一颜色通道值；/>表示第/>个像素点和第/>个像素点的RGB颜色通道平均差异因子；/>表示以自然常数为底的指数函数。
优选的，所述吸水棉条RGB图像中第个像素点和第/>个像素点的第一颜色通道的具体获取方法为：
将第个像素点和第/>个像素点三个颜色通道差值最大的颜色通道，记为第/>个像素点和第/>个像素点的第一颜色通道。
优选的，所述获取吸水棉条RGB图像中任意两个像素点的RGB颜色通道平均差异因子的具体公式为：
式中，表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示取绝对值。
优选的，所述获取吸水棉条RGB图像的若干个窗口，包括的具体方法为：
对于吸水棉条RGB图像中任意一个像素点，将所述像素点作为窗口中心像素点，获取吸水棉条RGB图像的一个窗口，且所述窗口大小为预设参数；进而获得吸水棉条RGB图像的若干个窗口。
优选的，所述每个窗口的高斯环绕函数的尺度参数的具体公式为：
式中，表示吸水棉条RGB图像的第/>个窗口的高斯环绕函数的尺度参数；/>表示吸水棉条RGB图像的第/>个窗口中所有像素点的总数量；/>表示吸水棉条RGB图像的第/>个窗口的中心像素点灰度值；/>表示吸水棉条RGB图像的第/>个窗口中第/>个像素点灰度值；表示吸水棉条RGB图像中第/>个窗口的中心像素点和第/>个像素点的颜色相似度；/>表示以自然常数为底的指数函数。
优选的，所述每个窗口的高斯环绕函数的具体计算表达式为：
式中，表示吸水棉条RGB图像的第/>个窗口中坐标为/>像素点的高斯环绕函数值；/>表示坐标为/>像素点的横坐标；/>表示坐标为/>像素点的纵坐标；/>表示吸水棉条RGB图像的第/>个窗口的高斯环绕函数的尺度参数；/>表示以自然常数为底的指数函数。
优选的，所述根据吸水棉条光照图像和每个窗口的高斯环绕函数，获取不受光照影响的吸水棉条图像，包括的具体方法为：
将吸水棉条光照图像转换到对数域，则获取不受光照影响的吸水棉条图像的计算表达式为：
式中，表示不受光照影响的吸水棉条图像中坐标为/>的像素点灰度值；表示吸水棉条RGB图像的所有窗口总数量；/>表示吸水棉条RGB图像的第/>个窗口中坐标为/>像素点的高斯环绕函数值；/>表示吸水棉条光照图像中坐标为/>的像素点灰度值。
优选的，所述通过可视宫颈采样刷获取吸水棉条RGB图像，包括的具体方法为：
用户将可视宫颈采样刷的帽子打开，将刮头推出，然后将染色剂滴在吸水棉条上，再将刮头收回套筒；打开刮头的摄像头隔离套筒中的针孔摄像头，连接手机或电脑，将套筒插入阴道，推出刮头，根据图像显示，将宫颈管刮头插入宫颈口，待摄像头刚进入宫颈口后，再向前推两个刻度，使宫颈壁刮头充分与宫颈壁接触，旋转五到八圈，完成染色，通过可视宫颈采样刷中刮头的摄像头隔离套筒中的针孔摄像头拍摄获得染色之后的可视宫颈采样刷中宫颈管刮头上的吸水棉条图像，进而得到吸水棉条RGB图像；
所述可视宫颈采样刷由5个零件组成，分别为循环使用部分：可视手柄-主体，一次性使用部分：档位护管，套筒，刮头，帽子；刮头，其中，刮头有闭合状态和展开状态，刮头由宫颈管刮头和宫颈壁刮头、宫颈管比色卡和宫颈壁比色卡、吸水棉条、宫颈壁刮头支撑件、刮头支撑套管和摄像头隔离套筒。
本发明的技术方案的有益效果是：根据吸水棉条RGB图像中任意两个像素点的颜色相似度，获取吸水棉条RGB图像的每个窗口的高斯环绕函数的尺度参数，进而得到每个窗口的高斯环绕函数，以此得到不受光照影响的吸水棉条图像，解决了采集到的吸水棉条染色结果受到光照环境影响的问题，进而提高吸水棉条判断结果的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明一种基于比色卡的FRD宫颈染色实时比对方法的步骤流程图；
图2为本发明一种基于比色卡的FRD宫颈染色实时比对方法的可视宫颈采样刷示意图；
图3为本发明一种基于比色卡的FRD宫颈染色实时比对方法的刮头闭合状态示意图；
图4为本发明一种基于比色卡的FRD宫颈染色实时比对方法的刮头展开状态示意图。
图2、图3和图4中的标号为：1、可视采样刷-主体；2、套筒；3、档位护管；4、帽子；5、刷头；6、宫颈管刮头与比色卡；7、宫颈壁刮头支撑件；8、刮头支撑套管；9、摄像头隔离套筒；10、宫颈管刮头；11、宫颈管比色卡；12、宫颈壁比色卡；13、宫颈壁刮头。
具体实施方式
为了更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效，以下结合附图及较佳实施例，对依据本发明提出的一种基于比色卡的FRD宫颈染色实时比对方法，其具体实施方式、结构、特征及其功效，详细说明如下。在下述说明中，不同的“一个实施例”或“另一个实施例”指的不一定是同一实施例。此外，一或多个实施例中的特定特征、结构或特点可由任何合适形式组合。
除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。
下面结合附图具体的说明本发明所提供的一种基于比色卡的FRD宫颈染色实时比对方法的具体方案。
请参阅图1，其示出了本发明一个实施例提供的一种基于比色卡的FRD宫颈染色实时比对方法的步骤流程图，该方法包括以下步骤：
步骤S001：通过可视宫颈采样刷获取吸水棉条RGB图像。
需要说明的是，在用户采集吸水棉条图像时，由于用户所处的灯光环境种类各异，因此采集到的吸水棉条染色结果也会受到环境的影响从而导致对比结果不准确；利用对吸水棉条染色情况进行具体分析，量化出吸水棉条具体的着色情况，除去光照对于结果的干扰从而提高吸水棉条判断结果的准确性。
具体的，为了实现本实施例提出的一种基于比色卡的FRD宫颈染色实时比对方法，首先需要采集吸水棉条RGB图像，具体过程为：
用户使用可视宫颈采样刷进行宫颈癌FRD染色检测，在光亮的环境下，通过可视宫颈采样刷中刮头的摄像头隔离套筒中的针孔摄像头拍摄获得检测之后的可视宫颈采样刷中宫颈管刮头上的吸水棉条图像，使用机器识别将图像中的吸水棉条部分图像进行分割，得到吸水棉条图像；对吸水棉条图像进行中值滤波去噪得到吸水棉条RGB图像。其中，中值滤波为现有技术，此处不做过多赘述。
本实施例所使用的可视宫颈采样刷如图2所示，其中，可视宫颈采样刷由5个零件组成，分别为循环使用部分：可视手柄-主体，一次性使用部分：档位护管，套筒，刮头，帽子；刮头闭合状态如图3所示，刮头展开状态如图4所示，其中刮头由宫颈管刮头和宫颈壁刮头、宫颈管比色卡和宫颈壁比色卡、吸水棉条、宫颈壁刮头支撑件、刮头支撑套管和摄像头隔离套筒。
至此，通过上述方法得到吸水棉条RGB图像。
步骤S002：根据吸水棉条RGB图像获取吸水棉条光照图像和吸水棉条反射率图像。
需要说明的是，通过相机拍摄得到的吸水棉条图像会受到色温和周围环境的影响从而导致拍摄到的图像中的颜色并不会是吸水棉条真实的颜色从而导致于比色卡对比之后结果不会特别的准确，因此需要去除外在光照的影响来使拍摄到的吸水棉条颜色更加真实。通过Retinex理论可以知道，物体的颜色一般是具有一致性的，也就是说物体如果受到了不均匀的光照照射也不会影响物体本身的颜色，物体本身的颜色并不是由反射光的强度来决定，而是通过其自身对于各种波长光线的反射能力决定的。
具体的，根据RetinexL算法将吸水棉条RGB图像分解为吸水棉条光照图像和吸水棉条反射率图像；且三张图像存在的关系表达式为：
式中，表示吸水棉条RGB图像中坐标为/>的像素点灰度值；/>表示吸水棉条光照图像中坐标为/>的像素点灰度值；/>表示吸水棉条反射率图像中坐标为/>的像素点灰度值。
其中，根据RetinexL算法将吸水棉条RGB图像分解为吸水棉条光照图像和吸水棉条反射率图像为现有技术，本实施例不作过多赘述。
至此，通过上述方法得到吸水棉条光照图像和吸水棉条反射率图像。
步骤S003：根据吸水棉条RGB图像中任意两个像素点的颜色相似度，获取吸水棉条RGB图像的每个窗口的高斯环绕函数的尺度参数，进而得到每个窗口的高斯环绕函数。
需要说明的是，因为是将吸水棉条部分图像进行提取，而吸水棉条本身是纯白色的，可以根据反应之后图像中灰度值最高点来得到入射光的分量，在此处有现有技术可以用于提取入射光分量。因为入射光分量相对应的是图像中的低频分量，因此在提取出入射光分量之后对图像使用低通滤波器来得到图像中的低频分量，此手段可以用高斯环绕函数来达到此目的。因为吸水棉条在反应的过程中容易出现反应不均匀的情况，从而导致反应之后吸水棉条的颜色是不均匀的，存在一定差异的，因此需要在建立环绕函数时尽可能的降低这类问题。因为在模糊图像的过程中，如果某个区域的图像颜色分布不均匀，那么就需要模糊的更多一点，反之则可以少模糊一点从而保留更多的颜色特征，这样模糊之后的图像可以尽可能的保证颜色的一致性从而可以减少对于反应结果的误判。
1.获取吸水棉条RGB图像中任意两个像素点的颜色相似度。
需要说明的是，想要得到吸水棉条中的某个区域反应是否均匀，则需要更具各个像素点颜色通道之间的差值分布规则来得到像素点之间的颜色相似度，对于RGB图像来说，三个颜色通道的总差值一定，但RGB各通道的差值分布方式不同，会影响图像的视觉效果和颜色变化的感知。假设两种情况，第一种，三个颜色通道的总差值在60，RGB各通道差值都为20，在这种情况下，每个颜色通道的变化幅度相对较小，但是三个通道的总变化值相等，所以整体的颜色变化较为均匀。这可能导致图像颜色的整体平衡保持得相对较好，视觉效果可能会被认为是比较自然和一致的。第二种，假设R、G通道都一样，而B通道相差60，在这种情况下，蓝色通道的变化幅度明显大于红色和绿色通道。这可能导致图像的颜色偏向蓝色，而红色和绿色的颜色相对不变。结果可能是图像的整体色调偏向蓝色，可能会被视为冷色调。视觉效果可能会因为颜色偏移而产生较大的变化感知。总的来说，两种情况下的视觉效果会有所不同。第一种情况下，颜色变化相对均匀，可能会被认为更加自然和平衡。而第二种情况下，蓝色通道的较大变化可能会导致颜色偏移，可能会被认为是更明显的颜色变化。
具体的，使用imread函数来读取吸水棉条RGB图像，并获取吸水棉条RGB图像中每个像素点的RGB三个颜色通道值；对于吸水棉条RGB图像中第个像素点和第/>个像素点，将第/>个像素点和第/>个像素点三个颜色通道差值最大的颜色通道，记为第/>个像素点和第/>个像素点的第一颜色通道；则吸水棉条RGB图像中第/>个像素点和第/>个像素点的颜色相似度的计算表达式：
式中，表示吸水棉条RGB图像中第/>个像素点和第/>个像素点的颜色相似度；/>表示吸水棉条RGB图像中第/>个像素点和第/>个像素点的第一颜色通道差值；/>表示吸水棉条RGB图像中所有像素点中最大的第一颜色通道值；/>表示吸水棉条RGB图像中所有像素点中最小的第一颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示以自然常数为底的指数函数；/>表示取绝对值。
至此，获得吸水棉条RGB图像中任意两个像素点的颜色相似度。
2.获取吸水棉条RGB图像的每个窗口的高斯环绕函数。
需要说明的是，在得到像素点之间的颜色相似性之后，可以根据颜色相似性来得到相应的卷积核，卷积核的作用范围会受到尺度范围决定，卷积核的作用范围越大，那么估计的入射分量就比较平滑，能够更加有效的去除环境光的影响，则增强后的图像色彩就更加的真实。相反，如果尺度范围越小，则卷积核的作用范围越小，增强后图像中的细节就更加的突出。因此针对于吸水棉条中反应不均匀的区域，则需要增加相应的尺度参数来增大卷积核的范围从而模糊掉反应不均匀的区域使其色彩更加切合与充分反应区域。
预设一个参数，其中本实施例以/>为例进行叙述，本实施例不进行具体限定，其中/>根据具体实施情况而定。
具体的，对于吸水棉条RGB图像中任意一个像素点，将所述像素点作为窗口中心像素点，获取吸水棉条RGB图像的一个窗口，且所述窗口大小为预设参数；进而获得吸水棉条RGB图像的若干个窗口。其中，若所述像素点周围实际存在的像素点数量不满足预设参数/>时，那么以所述像素点周围实际存在的像素点数量为准，获取以所述像素点作为窗口中心像素点的窗口。
至此，获得吸水棉条RGB图像的若干个窗口。
进一步，吸水棉条RGB图像的第个窗口的高斯环绕函数的尺度参数的计算表达式为：
式中，表示吸水棉条RGB图像的第/>个窗口的高斯环绕函数的尺度参数；/>表示吸水棉条RGB图像的第/>个窗口中所有像素点的总数量；/>表示吸水棉条RGB图像的第/>个窗口的中心像素点灰度值；/>表示吸水棉条RGB图像的第/>个窗口中第/>个像素点灰度值；表示吸水棉条RGB图像中第/>个窗口的中心像素点和第/>个像素点的颜色相似度；/>表示以自然常数为底的指数函数。
根据吸水棉条RGB图像的第个窗口的高斯环绕函数的尺度参数，获取第/>个窗口的高斯环绕函数的计算表达式为：
式中，表示吸水棉条RGB图像的第/>个窗口中坐标为/>像素点的高斯环绕函数值；/>表示坐标为/>像素点的横坐标；/>表示坐标为/>像素点的纵坐标；/>表示吸水棉条RGB图像的第/>个窗口的高斯环绕函数的尺度参数；/>表示吸水棉条RGB图像的第/>个窗口的高斯环绕函数的尺度参数；/>表示以自然常数为底的指数函数。
至此，获得吸水棉条RGB图像的每个窗口的高斯环绕函数。
步骤S004：将不受光照影响的吸水棉条图像与比色卡进行比对得到测试结果。
需要说明的是，图像自身的光照图像决定了像素的动态范围，反射率图像反映了物体本身对各种波长的光的反射特性，也即为不受光照影响的图像；而物体本身对于各种波长的光的反射特性只会于物体表面的特性所决定，因此需要得到吸水棉条RGB图像中的反射率图像，则需要在计算过程中将吸水棉条RGB图像中的光照图像去掉即可。为了能够使得运算之后的结果更加符合人眼的亮度感知能力并且简化相关的计算量，将吸水棉条RGB图像换到对数域。
具体的，使用Python中的PIL库将吸水棉条光照图像转换到对数域，也即，其中，/>表示吸水棉条RGB图像中坐标为的像素点灰度值；/>表示吸水棉条光照图像中坐标为/>的像素点灰度值；表示吸水棉条反射率图像中坐标为/>的像素点灰度值；则获取不受光照影响的吸水棉条图像的计算表达式为：
式中，表示不受光照影响的吸水棉条图像中坐标为/>的像素点灰度值；表示吸水棉条RGB图像的所有窗口总数量；/>表示吸水棉条RGB图像的第/>个窗口中坐标为/>像素点的高斯环绕函数值；/>表示吸水棉条光照图像中坐标为/>的像素点灰度值。
至此，获得不受光照影响的吸水棉条图像。
进一步，在得到不受光照影响的吸水棉条图像之后，在利用机器视觉对不受光照影响的吸水棉条图像的颜色与比色卡之间进行配对，得到的测试结果。
至此，本实施例完成。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种基于比色卡的FRD宫颈染色实时比对方法，其特征在于，该方法包括以下步骤：
通过可视宫颈采样刷获取吸水棉条RGB图像；
根据吸水棉条RGB图像获取吸水棉条光照图像；
获取吸水棉条RGB图像中任意两个像素点的RGB颜色通道平均差异因子，根据RGB颜色通道平均差异因子获取吸水棉条RGB图像中任意两个像素点的颜色相似度；获取吸水棉条RGB图像的若干个窗口；根据吸水棉条RGB图像中每个窗口中任意两个像素点的颜色相似度，获取每个窗口的高斯环绕函数的尺度参数；根据每个窗口的高斯环绕函数的尺度参数，获取每个窗口的高斯环绕函数；
根据吸水棉条光照图像和每个窗口的高斯环绕函数，获取不受光照影响的吸水棉条图像。
2.根据权利要求1所述一种基于比色卡的FRD宫颈染色实时比对方法，其特征在于，所述根据吸水棉条RGB图像获取吸水棉条光照图像，包括的具体方法为：
根据RetinexL算法将吸水棉条RGB图像分解为吸水棉条光照图像和吸水棉条反射率图像；且三张图像存在的关系表达式为：
式中，表示吸水棉条RGB图像中坐标为/>的像素点灰度值；/>表示吸水棉条光照图像中坐标为/>的像素点灰度值；/>表示吸水棉条反射率图像中坐标为的像素点灰度值。
3.根据权利要求1所述一种基于比色卡的FRD宫颈染色实时比对方法，其特征在于，所述吸水棉条RGB图像中任意两个像素点的颜色相似度的具体公式为：
式中，表示吸水棉条RGB图像中第/>个像素点和第/>个像素点的颜色相似度；表示吸水棉条RGB图像中第/>个像素点和第/>个像素点的第一颜色通道差值；/>表示吸水棉条RGB图像中所有像素点中最大的第一颜色通道值；/>表示吸水棉条RGB图像中所有像素点中最小的第一颜色通道值；/>表示第/>个像素点和第/>个像素点的RGB颜色通道平均差异因子；/>表示以自然常数为底的指数函数。
4.根据权利要求3所述一种基于比色卡的FRD宫颈染色实时比对方法，其特征在于，所述吸水棉条RGB图像中第个像素点和第/>个像素点的第一颜色通道的具体获取方法为：
将第个像素点和第/>个像素点三个颜色通道差值最大的颜色通道，记为第/>个像素点和第/>个像素点的第一颜色通道。
5.根据权利要求1所述一种基于比色卡的FRD宫颈染色实时比对方法，其特征在于，所述获取吸水棉条RGB图像中任意两个像素点的RGB颜色通道平均差异因子的具体公式为：
式中，表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；/>表示吸水棉条RGB图像中第/>个像素点的/>颜色通道值；表示取绝对值。
6.根据权利要求1所述一种基于比色卡的FRD宫颈染色实时比对方法，其特征在于，所述获取吸水棉条RGB图像的若干个窗口，包括的具体方法为：
对于吸水棉条RGB图像中任意一个像素点，将所述像素点作为窗口中心像素点，获取吸水棉条RGB图像的一个窗口，且所述窗口大小为预设参数；进而获得吸水棉条RGB图像的若干个窗口。
7.根据权利要求1所述一种基于比色卡的FRD宫颈染色实时比对方法，其特征在于，所述每个窗口的高斯环绕函数的尺度参数的具体公式为：
式中，表示吸水棉条RGB图像的第/>个窗口的高斯环绕函数的尺度参数；/>表示吸水棉条RGB图像的第/>个窗口中所有像素点的总数量；/>表示吸水棉条RGB图像的第/>个窗口的中心像素点灰度值；/>表示吸水棉条RGB图像的第/>个窗口中第/>个像素点灰度值；/>表示吸水棉条RGB图像中第/>个窗口的中心像素点和第/>个像素点的颜色相似度；/>表示以自然常数为底的指数函数。
8.根据权利要求1所述一种基于比色卡的FRD宫颈染色实时比对方法，其特征在于，所述每个窗口的高斯环绕函数的具体计算表达式为：
式中，表示吸水棉条RGB图像的第/>个窗口中坐标为/>像素点的高斯环绕函数值；/>表示坐标为/>像素点的横坐标；/>表示坐标为/>像素点的纵坐标；/>表示吸水棉条RGB图像的第/>个窗口的高斯环绕函数的尺度参数；/>表示以自然常数为底的指数函数。
9.根据权利要求1所述一种基于比色卡的FRD宫颈染色实时比对方法，其特征在于，所述根据吸水棉条光照图像和每个窗口的高斯环绕函数，获取不受光照影响的吸水棉条图像，包括的具体方法为：
将吸水棉条光照图像转换到对数域，则获取不受光照影响的吸水棉条图像的计算表达式为：
式中，表示不受光照影响的吸水棉条图像中坐标为/>的像素点灰度值；/>表示吸水棉条RGB图像的所有窗口总数量；/>表示吸水棉条RGB图像的第/>个窗口中坐标为/>像素点的高斯环绕函数值；/>表示吸水棉条光照图像中坐标为/>的像素点灰度值。
10.根据权利要求1所述一种基于比色卡的FRD宫颈染色实时比对方法，其特征在于，所述通过可视宫颈采样刷获取吸水棉条RGB图像，包括的具体方法为：
用户将可视宫颈采样刷的帽子打开，将刮头推出，然后将染色剂滴在吸水棉条上，再将刮头收回套筒；打开刮头的摄像头隔离套筒中的针孔摄像头，连接手机或电脑，将套筒插入阴道，推出刮头，根据图像显示，将宫颈管刮头插入宫颈口，待摄像头刚进入宫颈口后，再向前推两个刻度，使宫颈壁刮头充分与宫颈壁接触，旋转五到八圈，完成染色，通过可视宫颈采样刷中刮头的摄像头隔离套筒中的针孔摄像头拍摄获得染色之后的可视宫颈采样刷中宫颈管刮头上的吸水棉条图像，进而得到吸水棉条RGB图像；
所述可视宫颈采样刷由5个零件组成，分别为循环使用部分：可视手柄-主体，一次性使用部分：档位护管，套筒，刮头，帽子；刮头，其中，刮头有闭合状态和展开状态，刮头由宫颈管刮头和宫颈壁刮头、宫颈管比色卡和宫颈壁比色卡、吸水棉条、宫颈壁刮头支撑件、刮头支撑套管和摄像头隔离套筒。
CN202311666024.6A 2023-12-07 2023-12-07 一种基于比色卡的frd宫颈染色实时比对方法 Active CN117368122B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311666024.6A CN117368122B (zh) 2023-12-07 2023-12-07 一种基于比色卡的frd宫颈染色实时比对方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311666024.6A CN117368122B (zh) 2023-12-07 2023-12-07 一种基于比色卡的frd宫颈染色实时比对方法

Publications (2)

Publication Number Publication Date
CN117368122A true CN117368122A (zh) 2024-01-09
CN117368122B CN117368122B (zh) 2024-02-13

Family

ID=89402609

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311666024.6A Active CN117368122B (zh) 2023-12-07 2023-12-07 一种基于比色卡的frd宫颈染色实时比对方法

Country Status (1)

Country Link
CN (1) CN117368122B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060072814A1 (en) * 2004-10-04 2006-04-06 Stmicroelectronics S.R.L. Color interpolation method of an image acquired by a digital sensor by directional filtering
US20110242126A1 (en) * 2010-04-05 2011-10-06 Microsoft Corporation Capturing image structure detail from a first image and color from a second image
US20110274338A1 (en) * 2010-05-03 2011-11-10 Sti Medical Systems, Llc Image analysis for cervical neoplasia detection and diagnosis
CN103578084A (zh) * 2013-12-09 2014-02-12 西安电子科技大学 基于亮通道滤波的彩色图像增强方法
CN103679157A (zh) * 2013-12-31 2014-03-26 电子科技大学 一种基于视网膜模型的人脸图像光照处理方法
US20140254931A1 (en) * 2013-03-06 2014-09-11 Yi-Shuan Lai Image Recovery Method
CN105072341A (zh) * 2015-08-12 2015-11-18 旗瀚科技股份有限公司 一种可用于机器视觉的高动态范围现实场景信息重构方法
CN108961222A (zh) * 2018-06-19 2018-12-07 江西大福医疗科技股份有限公司 一种基于阴道镜图像的宫颈癌早期筛查识别方法
CN111986120A (zh) * 2020-09-15 2020-11-24 天津师范大学 一种基于帧累加和多尺度Retinex的低光照图像增强优化方法
CN114403942A (zh) * 2022-02-16 2022-04-29 津泰(天津)医疗器械有限公司 一种可视的自采宫颈细胞取样器
CN115293989A (zh) * 2022-08-24 2022-11-04 东南大学 一种综合MsRcR和自动色阶的图像增强方法
WO2023070911A1 (zh) * 2021-10-27 2023-05-04 西安工程大学 一种基于自注意力的彩色纹理织物缺陷区域的检测方法

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060072814A1 (en) * 2004-10-04 2006-04-06 Stmicroelectronics S.R.L. Color interpolation method of an image acquired by a digital sensor by directional filtering
US20110242126A1 (en) * 2010-04-05 2011-10-06 Microsoft Corporation Capturing image structure detail from a first image and color from a second image
US20110274338A1 (en) * 2010-05-03 2011-11-10 Sti Medical Systems, Llc Image analysis for cervical neoplasia detection and diagnosis
US20140254931A1 (en) * 2013-03-06 2014-09-11 Yi-Shuan Lai Image Recovery Method
CN103578084A (zh) * 2013-12-09 2014-02-12 西安电子科技大学 基于亮通道滤波的彩色图像增强方法
CN103679157A (zh) * 2013-12-31 2014-03-26 电子科技大学 一种基于视网膜模型的人脸图像光照处理方法
CN105072341A (zh) * 2015-08-12 2015-11-18 旗瀚科技股份有限公司 一种可用于机器视觉的高动态范围现实场景信息重构方法
CN108961222A (zh) * 2018-06-19 2018-12-07 江西大福医疗科技股份有限公司 一种基于阴道镜图像的宫颈癌早期筛查识别方法
CN111986120A (zh) * 2020-09-15 2020-11-24 天津师范大学 一种基于帧累加和多尺度Retinex的低光照图像增强优化方法
WO2023070911A1 (zh) * 2021-10-27 2023-05-04 西安工程大学 一种基于自注意力的彩色纹理织物缺陷区域的检测方法
CN114403942A (zh) * 2022-02-16 2022-04-29 津泰(天津)医疗器械有限公司 一种可视的自采宫颈细胞取样器
WO2023155335A1 (zh) * 2022-02-16 2023-08-24 津泰(天津)医疗器械有限公司 一种可视的自采宫颈细胞取样器
CN115293989A (zh) * 2022-08-24 2022-11-04 东南大学 一种综合MsRcR和自动色阶的图像增强方法

Also Published As

Publication number Publication date
CN117368122B (zh) 2024-02-13

Similar Documents

Publication Publication Date Title
US9750326B2 (en) Transparency evaluation device, transparency evaluation method and transparency evaluation program
US20220237811A1 (en) Method for Testing Skin Texture, Method for Classifying Skin Texture and Device for Testing Skin Texture
US20230190180A1 (en) Enhancing pigmentation in dermoscopy images
CN110111245B (zh) 图像处理方法、装置、终端及计算机可读存储介质
WO2008042831A2 (en) Apparatus and method for analyzing skin using l*a*b* colorspace
KR20170083066A (ko) 기미 평가 장치, 기미 평가 방법 및 프로그램
US20170319147A1 (en) Cancerous lesion identifying method via hyper-spectral imaging technique
CN104000593B (zh) 一种测试皮肤的方法
CN112767392A (zh) 图像清晰度确定方法、装置、设备及存储介质
CN111062891A (zh) 图像处理方法、装置、终端及计算机可读存储介质
CN110495888B (zh) 一种基于中医舌、面图像的标准色卡及其应用
JPH0719839A (ja) 表面状態解析システム
CN110458787A (zh) 一种图像融合方法、装置及计算机存储介质
CN108601509A (zh) 图像处理装置、图像处理方法以及程序
CN114693676A (zh) 一种新材料纺织品漂白缺陷光学检测方法及装置
JP7141509B2 (ja) 脈波検出装置、脈波検出方法、および情報処理プログラム
CN109583330B (zh) 一种用于人脸照片的毛孔检测方法
CN117368122B (zh) 一种基于比色卡的frd宫颈染色实时比对方法
KR20120117604A (ko) 피부 상태 분석 방법 및 장치
CN116703798B (zh) 基于自适应干扰抑制的食管多模态内镜图像增强融合方法
CN111738984B (zh) 基于分水岭和种子填充的皮肤图像斑点评估方法及系统
CN113642358B (zh) 肤色检测方法、装置、终端和存储介质
CN116893134A (zh) 一种牛仔布耐色牢度测试方法
JPH0755447A (ja) 肌の表面状態の解析システム
JP2021058361A (ja) 生体情報取得装置およびプログラム

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