CN101877030B - 一种耳穴伪彩色分布图的计算机实现方法 - Google Patents
一种耳穴伪彩色分布图的计算机实现方法 Download PDFInfo
- Publication number
- CN101877030B CN101877030B CN 200910235475 CN200910235475A CN101877030B CN 101877030 B CN101877030 B CN 101877030B CN 200910235475 CN200910235475 CN 200910235475 CN 200910235475 A CN200910235475 A CN 200910235475A CN 101877030 B CN101877030 B CN 101877030B
- Authority
- CN
- China
- Prior art keywords
- point
- color
- acupuncture point
- ear
- resistance
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005259 measurement Methods 0.000 claims abstract description 36
- 238000012360 testing method Methods 0.000 claims abstract description 21
- 239000003086 colorant Substances 0.000 claims abstract description 17
- 238000001514 detection method Methods 0.000 claims abstract description 9
- 238000005192 partition Methods 0.000 claims abstract description 7
- 238000001467 acupuncture Methods 0.000 claims description 178
- 230000001537 neural effect Effects 0.000 claims description 18
- 206010048865 Hypoacusis Diseases 0.000 claims description 10
- 239000000523 sample Substances 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 7
- 101000579646 Penaeus vannamei Penaeidin-1 Proteins 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 2
- 238000012876 topography Methods 0.000 claims description 2
- 238000013507 mapping Methods 0.000 abstract description 8
- 238000012545 processing Methods 0.000 abstract description 6
- 238000004458 analytical method Methods 0.000 abstract description 4
- 230000006978 adaptation Effects 0.000 abstract 1
- 230000003902 lesion Effects 0.000 abstract 1
- 210000005036 nerve Anatomy 0.000 description 11
- 210000000624 ear auricle Anatomy 0.000 description 5
- 206010028980 Neoplasm Diseases 0.000 description 4
- 210000000056 organ Anatomy 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 210000001835 viscera Anatomy 0.000 description 4
- 241000746998 Tragus Species 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 210000003414 extremity Anatomy 0.000 description 3
- 210000000256 facial nerve Anatomy 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 210000001932 glossopharyngeal nerve Anatomy 0.000 description 3
- 210000002837 heart atrium Anatomy 0.000 description 3
- 210000000988 bone and bone Anatomy 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 210000004373 mandible Anatomy 0.000 description 2
- 230000036285 pathological change Effects 0.000 description 2
- 231100000915 pathological change Toxicity 0.000 description 2
- 230000001515 vagal effect Effects 0.000 description 2
- 240000006409 Acacia auriculiformis Species 0.000 description 1
- 101100433727 Caenorhabditis elegans got-1.2 gene Proteins 0.000 description 1
- 235000007926 Craterellus fallax Nutrition 0.000 description 1
- 240000007175 Datura inoxia Species 0.000 description 1
- 206010033733 Papule Diseases 0.000 description 1
- 230000004308 accommodation Effects 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 238000007630 basic procedure Methods 0.000 description 1
- 210000001217 buttock Anatomy 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 210000003792 cranial nerve Anatomy 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 210000005069 ears Anatomy 0.000 description 1
- 210000003754 fetus Anatomy 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000001802 infusion Methods 0.000 description 1
- 210000000653 nervous system Anatomy 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000009418 renovation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 210000002027 skeletal muscle Anatomy 0.000 description 1
- 210000001032 spinal nerve Anatomy 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
- 210000004291 uterus Anatomy 0.000 description 1
- 210000001186 vagus nerve Anatomy 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Images
Landscapes
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
一种耳穴伪彩色分布图的计算机实现方法,其特征在于,所述方法是在计算机中首先初始化,输入系统耳穴图;设立由测试界面子模块、测量控制子模块、结果分析子模块组成的信号检测模块;根据耳穴阻值进行伪彩色映射,并用插值的方式给出穴位点以外区域的颜色,经过图像处理,得到带有耳廓轮廓线、分区信息的耳穴伪彩色分布图。然后,按照规范的步骤进行穴位测量。本发明将耳穴电阻值通过伪彩色形式直观的显示出来,绘制耳廓伪彩色电阻分布图,便于进行病灶区域定位,同时还能根据实际需要,采取全部测量、部分测量和随机测量三种测量模式,从而提高了耳穴检测的适应范围。
Description
技术领域
本发明属于电子学、计算机和信号处理技术领域,特别涉及一种用耳穴伪彩色分布图的计算机实现方法,由此进一步诊断与耳穴相关的人体部位的功能状况。
背景技术
耳穴是耳廓皮肤表面与人体脏腑、经络、组织器官、四肢百骸相互沟通的部位,也是脉气输注的所在。在耳廓上反应人机体的生理功能和病理变化的部位统称为耳穴。耳穴的分布是有一定规律的,在耳廓前外侧面的排列像一个在子宫内倒置的胎儿,头部朝下,臀部及四肢朝上,胸部及躯干在中间。内脏器官在耳廓代表区的形态与器官自身的形态相似,往往呈“投影”的对应关系。耳前控制人体的前面、五脏六腑、组织器官和五官七窍,耳背控制人体的背面,神经系统,骨骼肌肉等运动系统。左耳控制人体左半身组织器官,右耳控制人体右半身组织器官。当人的脏腑、四肢、百骸、五官等器官患病时,即可导致与这些脏腑或器官相应的特定耳穴呈现病理阳性反应,如出现电阻改变、变色、变形、血管变化、丘疹、脱屑等。这些病理变化既包括颜色和形态的改变,也包括耳穴电学特性的变化,如:与疾病相关的耳穴电阻降低,因此可以通过测量耳穴电阻来进行疾病的定位。
1993年,《中华人民共和国国家标准GB/T13734-92耳穴名称与部位》的颁布,为耳穴的名称和定位给出了一个标准,2000年程红峰、程凯耳穴图在此基础上发展为耳廓上定位110个耳穴点,除了添加了少数耳穴外,基本同标准耳穴图。对耳廓结构的研究,将耳廓分为耳轮、耳舟、对耳轮、三角窝、耳甲艇、耳甲腔、耳屏、对耳屏、耳垂、耳背十个区域;而耳廓分布着枕小神经,三叉神经的耳颞神经下颌支,面神经、舌咽神经和迷走神经的混合支,耳大神经。因此,可以将上述110个耳穴划分到不同的解剖分区和神经分区中。将耳穴按照神经和解剖区域进行了划分是本发明的独特之处之一。
耳穴和疾病相关的中医理论早在中国古代四大经典著作之一的《黄帝内经》中就有提及,近年来,国内也有相关的一些发明专利。常义松、鲁官明、雷勇、张卫东、秦明的发明专利肿瘤耳穴探测电脑诊断仪(专利申请号CN91103392.0),采用具有汉字系统的中华学习机,直流放大的探测电路和利用中华学习机外设接口的现有引线扩展EPROM,使该诊断仪能够运行64K的肿瘤诊断专家系统,操作十分方便,运行可靠,速度快。它不仅可以定性诊断,而且可以确定肿瘤的部位,定性准确率90%,定位准确率80%,该仪器具有良好的功能价格,适于广大的医院、保健单位推广使用。但在计算机的使用上,此发明采用早已被淘汰的中华学习机且需要扩展EPROM,专家诊断系统的软件只能给出和肿瘤诊断相关的结果,而无法直观的看到整个耳廓电阻分布的信息。
和之前的专利相比,本发明克服了已有技术的不足之处,在对上述系统的处理、显示和分析方法上进行了技术革新改造,将耳穴电阻值通过伪彩色形式直观的显示出来,也就是将离散的耳穴点转化为了面,既能够更好的表现不同穴位点之间的联系,又能清楚的看到各个穴位点之间电阻变化的趋势,便于进行病灶区域定位。强大的分析功能还使对耳穴按照神经和解剖区域分区,观察各个区域耳穴阻值情况成为可能。同时,还能根据实际需要,采取全部测量、部分测量和随机测量三种测量模式,从而提高了耳穴检测的适应范围。
发明内容
一种耳穴伪彩色分布图的计算机实现方法,是在计算机中依次按以下步骤实现的:
步骤(1):计算机初始化:
输入已经规范过和分区后的系统耳穴图,分为正面、侧面和背面三个耳穴图。以像素为单位,正面耳穴图大小为110×200,背面耳穴图大小为60×155,侧面耳穴图大小为60×35;并以所述系统耳穴图的左上角为原点,给出所述系统耳穴图所有穴位的编号和坐标,编号为1~110号;各个解剖分区和神经分区中所述所有各穴位点的名称和标号,解剖分区编号为1~10号,依次为耳轮、耳舟、对耳轮、三角窝、耳甲艇、耳甲腔、耳屏、对耳屏、耳垂、耳背;神经分区编号为1~4号,依次为枕小神经,三叉神经的耳颞神经下颌支,面神经、舌咽神经和迷走神经的混合支,耳大神经分区;区域之间的轮廓线用黑色表示;
设立包含有测试界面子模块、测量控制子模块和结果分析子模块在内的信号检测模块:
测试界面子模块,设有基准点测量的界面、全部测量的界面、部分测量的界面和随机测量的界面,所述基准点是指耳背上耳根穴位,在用户选择所需测量方式后实时显示用耳穴电阻表示的测量数据;
测量控制子模块,设有:输入用户需要的测试数据的USB口,测量中断指令输入端口,测量中断指令输出端口,重测指令输入端口以及列表显示并更新指令的输入端口;
结果分析子模块,对测试结果进行统计分析,并建立以伪彩色形式绘制的耳穴电阻地形分布图,称为耳穴伪彩色分布图,在所属结果分析子模块中设定;
耳穴图内,穴位的最小和最大电阻值,在0kΩ和100kΩ之间,而耳穴的电阻用Rx表示,其中:
基准点的电阻Rr为30kΩ,用黄色表示;
当Rx<Rr时,电阻为Rx的穴位点对应的颜色R、G、B为:
当Rr≤Rx<50时,电阻为Rx的穴位点对应的颜色R、G、B为:
当50≤Rx<100-Rr时,电阻为Rx的穴位点对应的颜色R、G、B为:
当100-Rr≤Rx<100时,电阻为Rx的穴位点对应的颜色R、G、B为:
对需要加强的解剖分区和神经分区内,按照上述伪彩色方式进行加亮,而对不需要加强的区域内,各点的R、G、B为:
对于穴位点以外的区域,采用三角形面片插值方法进行插值,将相邻穴位点分组,每三点分为一组,穴位点作为三角形的三个顶点,构成三角形面片,插值计算其间每个像素的颜色,以便用伪彩色显示耳廓区域的颜色R、G、B的分布。
三角形面片插值方法,穴位点作为三角形的三个顶点,如果其中有一个穴位点未测量时,将其颜色赋值为三角形中与其编号最接近的穴位点的颜色。如果其中有两个穴位点未测量时,将这两个点的颜色赋值为三角形中已测的穴位点的颜色。
步骤(2),按以下步骤进行穴位检测:
步骤(2.1),在耳背的上耳根耳穴上放置探笔,测量其电阻,按下探笔上的开关,把结果送入计算机;
步骤(2.2),若选择全部测量则把探笔依次置于1~110号耳穴位置处,用步骤(2.1)所述的方法测量各穴位点的电阻值;
步骤(2.3),若选择部分测量则在所述的耳穴的1~10个解剖分区中任意选择一个或多个解剖分区,用步骤(2.1)所述的方法测量各穴位点的电阻值;
步骤(2.4),若选择随机测量则在所述的1~110号耳穴中任意选择一个或多个穴位点,用步骤(2.1)所述的方法测量各穴位点的电阻值;
步骤(2.5),用上述三种方式之一测量完毕后,在系统软件界面上点击“显示结果”按键,进行统计分析,显示测量结果;
步骤(2.6),存储测量数据,测试结束。
一种耳穴伪彩色分布图的计算机实现方法,经过三角形面片插值后的伪彩色分布图,需要对耳廓外的部分进行裁剪,显示耳廓轮廓线,进而得到带有耳廓轮廓线、分区信息的耳穴伪彩色分布图,图像处理方法是依次按以下步骤实现的:
步骤(1):区域编码:
步骤(1.1):对耳廓中的每个解剖分区给定一个分区代码,分别为0x1,0x2,0x4,0x8,0x10,0x20,0x40,0x80,0x100,0x200;
步骤(1.2):对于所述区域编号对应的至少两个分区代码进行按位“或”运算,得到区域标识码;若用户只选择了一个分区时,将这个分区的分区代码作为区域标识码;
步骤(2):区域解码:
对所述区域标识码与已知的至少两个解剖分区的分区代码进行按位“与”运算,结果每个解剖分区中的非零的区域为测量区域,得到测量的解剖分区的编号;
步骤(3):裁剪;
步骤(3.1):读入一张与待显示的伪彩色图一一对应的系统耳穴分区图,轮廓线用黑色表示,耳廓外的部分用白色表示,不同分区使用不同颜色标记;
步骤(3.2):对系统耳穴分区图,按x方向和y方向遍历,获取它在位置(x,y)处的颜色和在待显示的伪彩色图中位置(x,y)处的点的颜色;
步骤(3.3):若系统耳穴分区图中(x,y)处颜色为白色,将伪彩色图的对应点涂为白色;
步骤(3.4):若系统耳穴分区图中(x,y)处颜色为黑色,这一点是耳廓分区的轮廓线,将伪彩色图的对应点涂为黑色;
步骤(3.5):若系统耳穴分区图中(x,y)处的颜色不是白色或黑色,根据步骤(2)所述的解剖分区的编号,进一步得到需要显示分区对应的颜色;
步骤(3.6):根据步骤(3.5)所述的系统耳穴分区图中(x,y)处的颜色,按以下步骤执行;
步骤(3.6.1):若系统耳穴分区图中的颜色为该解剖分区的颜色时,点(x,y)在需要显示的解剖分区中,将伪彩色图对应位置保留;
步骤(3.6.2):若系统耳穴分区图中的颜色不是这个解剖分区的颜色时,点(x,y)不在需要显示的解剖分区中,将待显示的伪彩色图中的对应位置涂为白色;
步骤(3.7),对待显示的伪彩色图中每一点的颜色进行逐点处理,就得到了经过裁剪的在计算机上显示的含有测量电阻阻值信息的耳穴图。
附图内容
图1为本发明的系统组成框图。
图2为本发明的系统耳穴图。
图3为本实施例的耳穴解剖分区图:R、G、B分别为红、绿、蓝三种颜色成分。
图4为本实施例的耳穴神经分区图。
图5为本实施例的伪彩色映射曲线图。
图6为本实施例的灰度映射曲线图。
图7为本实施例的系统控制构成框图。
图8为本实施例的插值算法示意图。
图9为本实施例的耳穴检测基本流程图。
具体实施方式
本发明所述的规范化的系统耳穴图,分别为正面耳穴图、背面耳穴图和侧面耳穴图,三张耳穴图的大小,以像素为单位,分别为110×200,60×155和60×35,如图2所示。其中,正面耳穴图共有穴位点94个,背面耳穴图共有穴位点12个,侧面耳穴图共有穴位点4个。正面耳穴图、背面耳穴图和侧面耳穴图中,耳穴名称、编号和坐标对应关系如表1、2和表3所示。
表1 正面耳穴图穴位点编号与坐标对应表
表2 背面耳穴图穴位点坐标对应表
表3 侧面耳穴图穴位点坐标对应表
对耳廓解剖区域划分,则是在中医标准耳穴图将耳廓正面分为八个分区的基础上再进一步细分,将较大的耳甲区分为耳甲艇和耳甲腔,再加上耳背分为一个区,整个耳廓共分为十个区,即耳轮,耳舟,对耳轮,三角窝,耳甲艇,耳甲腔,耳屏,对耳屏,耳垂和耳背。具体分区如图3所示,各个解剖区包含的穴位点和标号如表4所示。
表4 各个解剖分区中穴位点列表
耳穴图神经分区是根据耳廓实际神经分布划分的。耳垂、耳轮、耳舟及对耳轮区,主要是脊神经,也就是耳大神经和枕小神经支配。耳甲区为脑神经,也就是耳颞神经和迷走神经耳支、舌咽神经与面神经的混合分布支支配。三角窝内神经分布极为丰富,几乎所有支配耳廓的神经都有分支至三角窝内。根据上述理论,对耳廓按照神经进行了划分,如图4所示。对全部110个穴位点进行分区,结果如表5所示。
表5 各个神经分区中穴位点列表
本实施例的软件部分实现对测试结果的显示、分析和管理的信号检测模块。
其组成如图7所示。
信号检测模块包括基准点、全部耳穴、部分耳穴和随机耳穴的检测、处理和结果显示,其具体功能为:
1.测试界面:包括基准点测量界面、全部测量界面、部分测量界面和随机测量界面,以及各种测量方式的测试数据的实时显示,如当前测试穴位点和测量阻值。在测试前,可以根据需求选择合适的测量方式,如选择部分测量或随机测量时,用户可以选择特定的耳穴分区或特定的穴位点进行测量。
2.测试控制:通过USB连接系统硬件,实现耳穴电阻的测量。开始测量后,当测量一个穴位点完毕后,用户需按下探笔上的按键,测量结果即通过USB传入计算机。此外,还能实现对测量的暂停和重测某些特定穴位点。测量数据实时在软件中进行列表显示并更新,用户能够随时看到最新一次的测量数据。
3.结果分析:包括对测试数据的统计分析和以伪彩色形式绘制的耳穴电阻分布图。数据统计分析能够给出所有测量点和按照解剖分区和神经分区的均值和方差,给出所有测量点中正常点、低阻点和超低阻点的点数,各阻值段电阻占总电阻数的百分比。耳穴电阻分布地形图使用不同的颜色标示耳穴电阻,按照伪彩色映射的方式,将各个穴位点的阻值对应为颜色。测试时,穴位的最小和最大电阻值在硬件电路中被限制在0和100kΩ之间,基准点上耳根处的电阻Rr约为30kΩ左右。电阻为30kΩ处映射为黄色,即黄色=红色r+绿色g,另外一侧的转折位置是与30kΩ这一点关于50kΩ对称的地方,即70kΩ处,映射为青色=蓝色b+绿色g,映射曲线如图5所示。其它穴位点阻值所对应的颜色值就可得出。
根据耳穴电阻Rx的取值范围不同,可分为以下四种情况:
1)当Rx<Rr时,电阻为Rx的穴位点对应的颜色R、G、B为:
2)当Rr≤Rx<50时,电阻为Rx的穴位点对应的颜色R、G、B为:
3)当50≤Rx<100-Rr时,电阻为Rx的穴位点对应的颜色R、G、B为:
4)当100-Rr≤Rx<100时,电阻为Rx的穴位点对应的颜色R、G、B为:
按照电阻的取值范围,就能把电阻阻值转化为颜色,并在软件中显示出来,使阻值信息的表示更为直观和明晰。
为了更好的表示各个神经分区和解剖分区的电阻分布情况,还能够根据需要,对某个解剖分区或神经分区进行加亮,对其他不需要加亮的分区用灰色表示,阻值越大则对应的颜色越深。映射曲线图如图6所示,映射关系如下:
公式(1)、(2)、(3)、(4)和(5)能够给出耳廓中所有穴位点的颜色,本发明在单个穴位点用伪彩色显示的基础上,还对穴位点以外的区域进行了插值,从而能够用伪彩色显示整个耳廓区域的情况。插值时,首先将相邻耳穴点分组,每三个点分为一组,先根据三个穴位点的阻值计算出三角形三个顶点的颜色,构成三角形面片,然后插值计算其间每个像素的颜色,得到彩色的三角形面片,进而显示整个耳廓上耳穴阻值的分布状况,使用户能够直观的得到耳廓阻值区间的信息。该部分针对不同情况,采用的算法包括线段插值算法和三角形面片插值算法。
当已知一条线段两个端点的颜色值时,使用线段插值算法计算线段上每一点的颜色值。图8中的直线AB所示,如果已知图中A、B两点的颜色值分别为(R1,G1,B1)和(R2,G2,B2),那么E点的颜色值(R3,G3,B3)可以使用公式(6)-(8)来进行计算。
当已知三角形三个顶点的颜色值时,则使用三角形面片插值算法来计算三角形内每一点的颜色值。如图8所示,首先在三角形三个顶点中找到y坐标位于中间的点(如图8中A点),沿此点作x方向的水平线,与另外一条边相交于点D,点D的颜色值可以利用线段插值算法,根据B、C点的颜色值计算出来。这样就得到了上下两个三角形,这两个三角形都有一条边和x方向平行,插值就变成分别对两个三角形进行插值。以下方的三角形ABD为例,三角形内任意一点颜色值的计算方法是,先找到经过该点的x方向的水平线,如图中线段EF;找到水平线和三角形非水平边的交点,如图中E点和F点,然后使用线段插值算法计算G点颜色值。根据A、B两点的颜色值计算E点的颜色值,根据B、D两点的颜色值插值计算F点的颜色值;然后根据E、F两点的颜色值插值计算出G点的颜色值,这样,整个三角形中每一点的颜色值都用这种算法计算得到了。
当三角形的三个顶点中有未测量点时,这就意味着该三角形是边界区域,将未测量点设置为和已测点相同的颜色,再进行面片插值。
由于三角形面片插值是对全图进行插值的,即对耳廓轮廓外和待测区域外都进行了插值操作,能够得到一个全平面的伪彩色图。为了只对图中需要的部分进行显示,需要对耳廓外的部分进行裁剪,显示耳廓轮廓线,进而得到带有耳廓轮廓线、分区信息的耳穴伪彩色分布图,图像处理过程是通过首先对耳廓解剖分区、神经分区进行编码,然后对图像进行裁剪实现的。
区域编码是对耳廓中的每个区都制定一个区域代码,1~10区的代码,用十六进制表示,分别为0x1,0x2,0x4,0x8,0x10,0x20,0x40,0x80,0x100,0x200。这些编码在转化为二进制时,都有一个共同的特点,即有一位为1,其余位都为0。不同的分区取1的位不同,最终测量区域的标识码为该区域中包含的若干分区代码进行按位“或”运算的结果。例如,当测量区域为1,2,5区时,分区代码分别为0x1,0x2,0x10。区域标识码为0x1|0x2|0x10=0x13。这个区域标识码将作为参数进行传递,用于后续绘图操作。若只选择了一个区,则区域标识码为该区的区域代码。
图像裁剪之前,需要先对区域标识码进行解码,即通过对区域标识码进行某种运算,获得测量区域的编号,这个过程与区域编码是一个相反的过程,然后才能进行裁剪。首先,将区域标识码依次与已知的各个区域的代码进位进行“与”运算,结果非零的区域为测量区域。以区域标识码0x13为例,计算的过程是分别对1~10区的分区代码进行“与”运算,以计算1~5区的结果为例:
1区:0x1 & 0x13=0x1
2区:0x2 & 0x13=0x2
3区:0x4 & 0x13=0
4区:0x8 & 0x13=0
5区:0x10 & 0x13=0x10
同理,6~10区的计算结果都为0。
这样,就获得了测量区域的编号为1,2和5了。这种先编码后解码的方式,能够大大减少测量区域编号作为参数传递的数据量,也便于软件进行统一的分析和处理。然后,根据前面解码出的区域标识码对全平面伪彩色图进行裁剪,只留下测量区域的伪彩色分布图信息。最后,开始裁剪,软件先读入一张系统耳穴分区图,这张分区图的特点是它和待显示的伪彩色图是一一对应的,轮廓线用黑色表示,耳廓外的部分用白色表示,不同的分区使用特定的某种颜色进行标记,分别对两张图的坐标进行遍历,获取该点在系统耳穴分区图中的颜色就能确定该点在待绘制的全平面伪彩色图中的分区。此外,还能利用系统耳穴分区图获取耳廓分区轮廓线。整个裁减的算法流程如下:
1)读入一张与待显示的伪彩色图一一对应的系统耳穴分区图,轮廓线用黑色表示,耳廓外的部分用白色表示,不同分区使用不同颜色标记;
2)对系统耳穴分区图,按x方向和y方向遍历,获取它在位置(x,y)处的颜色和在待显示的伪彩色图中位置(x,y)处的点的颜色;
3)若系统耳穴分区图中(x,y)处颜色为白色,将伪彩色图的对应点涂为白色;
4)若系统耳穴分区图中(x,y)处颜色为黑色,这一点是耳廓分区的轮廓线,将伪彩色图的对应点涂为黑色;
5)若系统耳穴分区图中(x,y)处的颜色不是白色或黑色,根据步骤(2)所述的解剖分区的编号,进一步得到需要显示分区对应的颜色;
6)根据步骤5)得到的系统耳穴分区图中(x,y)处的颜色,按以下步骤执行
a)若系统耳穴分区图中的颜色为该解剖分区的颜色时,点(x,y)在需要显示的解剖分区中,将伪彩色图对应位置保留;
b)若系统耳穴分区图中的颜色不是这个解剖分区的颜色时,点(x,y)不在需要显示的解剖分区中,将待显示的伪彩色图中的对应位置涂为白色;
举例进行说明如下:第1区标记的颜色是红色,(R1,G1,B1)=(255,0,0),也就是当系统耳穴分区图中坐标为(x,y)的点的颜色为红色,即(R0,G0,B0)=(255,0,0)时,这一点位于第1区。由于伪彩色图和原始分区图的一一对应关系,伪彩色图在(x,y)点处也是属于第1区。当第1区是测量区域时,在伪彩色图中保留这一点,即位置(x,y)处的点颜色(R,G,B)不变;当第1区不是测量区域时,在伪彩色图中将其设为白色。
这样,就获得了完整的带分区轮廓线的耳廓图和耳穴电阻地形图的伪彩色表征。
此外,在用户点击测量结果或将鼠标停留在耳穴地形图中耳穴的位置时,该点的测量结果将会在图中的对应位置显示出来。
本实施例的耳穴检测基本流程如图9所示,适用于各种方式的耳穴电阻测量,其具体步骤包括:
步骤(1):计算机初始化:输入已经规范过和分区后的系统耳穴图,设立包含有测试界面子模块、测量控制子模块和结果分析子模块在内的信号检测模块。
步骤(2):按以下步骤进行穴位检测
步骤(2.1),在耳背的上耳根耳穴上放置探笔,测量其电阻,按下探笔上的开关,把结果送入计算机;
步骤(2.2),若选择全部测量则把探笔依次置于1~110号耳穴位置处,用步骤(2.1)所述的方法测量各穴位点的电阻值;
步骤(2.3),若选择部分测量则在所述的耳穴的1~10个解剖分区中任意选择一个或多个解剖分区,用步骤(2.1)所述的方法测量各穴位点的电阻值;
步骤(2.4),若选择随机测量则在所述的1~110号耳穴中任意选择一个或多个穴位点,用步骤(2.1)所述的方法测量各穴位点的电阻值;
步骤(2.5),用上述三种方式之一测量完毕后,在系统软件界面上点击“显示结果”按键,进行统计分析,显示测量结果;
步骤(2.6),存储测量数据,测试结束。
Claims (4)
1.一种耳穴伪彩色分布图的计算机实现方法,其特征在于,所述方法是在计算机中依次按以下步骤实现的:
步骤(1):计算机初始化:
输入已经规范过和分区后的系统耳穴图,并以所述系统耳穴图的左上角为原点,给出所述系统耳穴图中所有穴位的编号和坐标,所述编号为1~110号;各个解剖分区和神经分区中所述所有各穴位点的名称和标号,解剖分区编号为1~10号,神经分区编号为1~4号;区域之间的轮廓线用黑色线表示;
设立:包含有测试界面子模块、测量控制子模块和结果分析子模块在内的信号检测模块,其中:
测试界面子模块,设有:基准点测量的界面、全部测量的界面、部分测量的界面和随机测量的界面,所述基准点是指耳背上耳根穴位,在用户选择所需测量方式后实时显示用耳穴电阻表示的测量数据;
测量控制子模块,设有:输入用户需要的测试数据的USB口、测量中断指令输入端口、测量中断指令输出端口、重测指令输入端口以及列表显示并更新指令的输入端口;
结果分析子模块,对测试结果进行统计分析,并建立以伪彩色形式绘制的耳穴电阻地形分布图,称为耳穴伪彩色分布图,在所述结果分析子模块中设定;耳穴图内,穴位的最小和最大电阻值,在0kΩ和100kΩ之间,而耳穴的电阻用Rx表示,其中:
基准点的电阻Rr为30kΩ,用黄色表示;
当Rx<Rr时,电阻为Rx的穴位点对应的颜色R、G、B为:
当Rr≤Rx<50时,电阻为Rx的穴位点对应的颜色R、G、B为:
当50≤Rx<100-Rr时,电阻为Rx的穴位点对应的颜色R、G、B为:
当100-Rr≤Rx<100时,电阻为Rx的穴位点对应的颜色R、G、B为:
对需要加强的解剖分区和神经分区,按照上述伪彩色方式进行加亮,而对不需要加强的区域内,各点的颜色R、G、B为:
对于穴位点以外的区域,采用三角形面片插值方法进行插值,将相邻穴位点分组,每三点分为一组,穴位点作为三角形的三个顶点,构成三角形面片,插值计算其间每个像素的颜色,以便用伪彩色显示耳廓区域的颜色R、G、B的分布;
步骤(2):经过三角形面片插值后的伪彩色分布图,需要对耳廓外的部分进行裁剪,显示耳廓轮廓线,进而得到带有耳廓轮廓线、分区信息的耳穴伪彩色分布图,图像处理方法是依次按以下步骤实现的:
步骤(2.1):区域编码:
步骤(2.1.1):对耳廓中的每个解剖分区给定一个分区代码,分别为0x1,0x2,0x4,0x8,0x10,0x20,0x40,0x80,0x100,0x200;
步骤(2.1.2):对于所述区域编号对应的至少两个分区代码进行按位“或”运算,得到区域标识码;
步骤(2.2):区域解码:
对所述区域标识码与已知的至少两个解剖分区的分区代码进行按位“与”运算,结果每个解剖分区中的非零的区域为测量区域,得到测量的解剖分区的编号;
步骤(2.3):裁剪:
步骤(2.3.1):读入一张与待显示的伪彩色图一一对应的系统耳穴分区图,轮廓线用黑色表示,耳廓外的部分用白色表示,不同分区使用不同颜色标记;
步骤(2.3.2):对系统耳穴分区图,按x方向和y方向遍历,获取它在位置(x,y)处的颜色和在待显示的伪彩色图中位置(x,y)处的点的颜色;
步骤(2.3.3):若系统耳穴分区图中(x,y)处颜色为白色,将伪彩色图的对应点涂为白色;
步骤(2.3.4):若系统耳穴分区图中(x,y)处颜色为黑色,这一点是耳廓分区的轮廓线,将伪彩色图的对应点涂为黑色;
步骤(2.3.5):若系统耳穴分区图中(x,y)处的颜色不是白色或黑色,根据步骤(2.2)所述的解剖分区的编号,进一步得到需要显示分区对应的颜色;
步骤(2.3.6):根据步骤(2.3.5)所述的系统耳穴分区图中(x,y)处的颜色,按以下步骤执行;
步骤(2.3.6.1):若系统耳穴分区图中的颜色为该解剖分区的颜色时,点(x,y)在需要显示的解剖分区中,将伪彩色图对应位置保留;
步骤(2.3.6.2):若系统耳穴分区图中的颜色不是这个解剖分区的颜色时,点(x,y)不在需要显示的解剖分区中,将待显示的伪彩色图中的对应位置涂为白色;
步骤(2.3.7),对待显示的伪彩色图中每一点的颜色进行逐点处理,就得到了进行了裁剪的在计算机上显示的含有测量电阻阻值信息的耳穴图;
步骤(3),按以下步骤进行穴位检测:
步骤(3.1),在耳背的上耳根耳穴上放置探笔,测量其电阻,按下探笔上的开关,把结果送入计算机;
步骤(3.2),若选择全部测量,则把探笔依次置于1~110号耳穴位置处,用步骤(3.1)所述的方法测量各穴位点的电阻值;
步骤(3.3),若选择部分测量,则在所述的耳穴的1~10个解剖分区中任意选择一个或多个解剖分区,用步骤(3.1)所述的方法测量各穴位点的电阻值;
步骤(3.4),若选择随机测量则在所述的1~110号耳穴中任意选择一个或多个穴位点,用步骤(3.1)所述的方法测量各穴位点的电阻值;
步骤(3.5),用上述三种方式之一测量完毕后,在系统软件界面上点击“显示结果”按键,进行统计分析,显示测量结果;
步骤(3.6),存储测量数据,测试结束。
2.根据权利要求1所述的耳穴伪彩色分布图的计算机实现方法,其特征在于,所述的三角形面片插值方法,穴位点作为三角形的三个顶点,若其中有一个穴位点未测量时,将其颜色赋值为三角形中与其编号最接近的穴位点的颜色。
3.根据权利要求1所述的耳穴伪彩色分布图的计算机实现方法,其特征在于,所述的三角形面片插值方法,穴位点作为三角形的三个顶点,若其中有两个穴位点未测量时,将这两个点的颜色赋值为三角形中已测的穴位点的颜色。
4.根据权利要求1所述耳穴伪彩色分布图的计算机实现方法,其特征在于,在步骤(2.1.2)中,当用户只选择了一个分区时,将这个分区的分区代码作为区域标识码。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910235475 CN101877030B (zh) | 2009-10-26 | 2009-10-26 | 一种耳穴伪彩色分布图的计算机实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910235475 CN101877030B (zh) | 2009-10-26 | 2009-10-26 | 一种耳穴伪彩色分布图的计算机实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101877030A CN101877030A (zh) | 2010-11-03 |
CN101877030B true CN101877030B (zh) | 2013-04-03 |
Family
ID=43019585
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910235475 Expired - Fee Related CN101877030B (zh) | 2009-10-26 | 2009-10-26 | 一种耳穴伪彩色分布图的计算机实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101877030B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102908140A (zh) * | 2012-10-31 | 2013-02-06 | 高少斌 | 可视化经络扫描装置及其视频采集与叠加方法 |
CN104083167B (zh) * | 2014-07-02 | 2016-05-04 | 高少斌 | 用于显示人体经络图像的对比方法 |
CN104083166A (zh) * | 2014-07-02 | 2014-10-08 | 高少斌 | 用于显示人体经络图像的采集及拼接方法 |
CN104799828B (zh) * | 2015-05-15 | 2018-11-02 | 春泉健康管理(上海)有限公司 | 耳纹学识别远程图像身体诊断系统 |
CN107145741B (zh) * | 2017-05-05 | 2020-06-05 | 必应(上海)医疗科技有限公司 | 基于图像分析的耳诊数据采集方法及装置 |
CN108078765A (zh) * | 2017-12-15 | 2018-05-29 | 泗洪县正心医疗技术有限公司 | 一种基于标记点投影实现腧穴定位的方法 |
CN109730923B (zh) * | 2019-03-04 | 2021-02-19 | 黑龙江中医药大学 | 辅助耳穴压籽的耳部穴位自动定位装置、定位系统及定位方法 |
CN112330802B (zh) * | 2020-11-05 | 2023-08-01 | 中国中医科学院中医药信息研究所 | 基于穴位的皮区模型的建立方法及装置 |
CN114099322B (zh) * | 2021-12-06 | 2023-05-26 | 贵州中医药大学第一附属医院 | 一种便于定位耳穴的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1072031A (zh) * | 1992-10-24 | 1993-05-12 | 中医研究院针灸研究所 | 热像图穴位诊断的计算机分区方法 |
CN1081780A (zh) * | 1993-06-05 | 1994-02-09 | 山东大学 | 生物全息电图诊断计算机处理系统 |
CN101292864A (zh) * | 2008-05-28 | 2008-10-29 | 青岛大学医学院附属医院 | 人体疾病信息手区成像的影像表达测定方法及其装置 |
-
2009
- 2009-10-26 CN CN 200910235475 patent/CN101877030B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1072031A (zh) * | 1992-10-24 | 1993-05-12 | 中医研究院针灸研究所 | 热像图穴位诊断的计算机分区方法 |
CN1081780A (zh) * | 1993-06-05 | 1994-02-09 | 山东大学 | 生物全息电图诊断计算机处理系统 |
CN101292864A (zh) * | 2008-05-28 | 2008-10-29 | 青岛大学医学院附属医院 | 人体疾病信息手区成像的影像表达测定方法及其装置 |
Also Published As
Publication number | Publication date |
---|---|
CN101877030A (zh) | 2010-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101877030B (zh) | 一种耳穴伪彩色分布图的计算机实现方法 | |
CN105824418B (zh) | 一种基于非对称视觉诱发电位的脑-机接口通讯系统 | |
CN204654887U (zh) | 基于云端的中医脉象辅助诊断系统 | |
CN100547626C (zh) | 人体穴位识别与考试系统 | |
CN102222153A (zh) | 中医机器问诊数量化辩证诊断方法 | |
CN109416939B (zh) | 面向群体应用的经颅脑图谱生成方法、预测方法及其装置 | |
CN101596128A (zh) | 具有健康状态指示符的患者监视系统 | |
CN103150023A (zh) | 一种基于脑机接口的光标控制系统及方法 | |
CN105868560B (zh) | 一种社会心理学调查系统 | |
CN101972149B (zh) | 视触觉测试仪及视触觉敏感性测试方法 | |
CN102609090A (zh) | 采用脑电时频成分双重定位范式的快速字符输入方法 | |
WO2018214520A1 (zh) | 视频数据和运动效果动画的同步方法及装置 | |
CN105718730B (zh) | 对受试者的疼痛量化评估方法以及实施该方法的系统 | |
CN106175747A (zh) | 利用导程的差动电压产生导程心电图信号的方法与系统 | |
EP2918226A1 (en) | Method and device for simplifying information obtained from electrical impedance tomography | |
CN102727996A (zh) | 智能高血压针灸诊疗方法及系统 | |
KR20140075258A (ko) | 스마트통신기기와 피부관리기기를 연동한 피부관리 시스템 및 그 제어방법 | |
CN111694425A (zh) | 一种基于ar-ssvep的目标识别方法及系统 | |
CN105943019A (zh) | 一种人体经络穴位生物电信号测量数据的三维分析方法 | |
CN201188625Y (zh) | 一种对体育锻炼的运动量进行网络化积分奖励的系统 | |
CN117711611B (zh) | 一种基于场景融合和mr的MDT远程会诊系统及方法 | |
Spence et al. | Using color to code quantity in spatial displays. | |
CN201111409Y (zh) | 量测显示装置 | |
CN103705216B (zh) | 等距与类别量表结合的中医舌色感知量化分级方法 | |
CN112656379B (zh) | 一种经穴脏腑平衡检测系统、方法及存储介质 |
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: 20130403 Termination date: 20141026 |
|
EXPY | Termination of patent right or utility model |