图像处理方法、图像处理装置、以及程序与流程

1.本公开的技术涉及图像处理方法、图像处理装置、以及程序。


背景技术：

2.期待利用眼底图像分析涡静脉(美国专利第8636364号说明书)。


技术实现要素：

3.本公开的技术的第一方面的图像处理方法为由处理器执行的图像处理方法，其包括：从被检查眼的眼底图像检测多个涡静脉的位置，算出检测出的所述多个涡静脉的位置的分布中心。
4.本公开的技术的第二方面的图像处理装置具备存储器、以及与所述存储器连接的处理器，所述处理器从被检查眼的眼底图像检测多个涡静脉的位置，算出检测出的所述多个涡静脉的位置的分布中心。
5.本公开的技术的第三方面的程序使计算机执行如下的步骤：从被检查眼的眼底图像检测多个涡静脉的位置，算出检测出的所述多个涡静脉的位置的分布中心。
附图说明
6.图1是眼科系统100的框图。
7.图2是示出眼科装置110的整体构成的概略构成图。
8.图3a是示出眼科装置110拍摄被检查眼12的眼底的拍摄范围的第一图。
9.图3b是示出眼科装置110拍摄被检查眼12的眼底的拍摄范围的第二图，是利用该拍摄得到的眼底的图像。
10.图3c是示出眼科装置110拍摄被检查眼12的眼底的拍摄范围的第三图。
11.图4是示出服务器140的构成的框图。
12.图5是服务器140的cpu162的功能框图。
13.图6是由服务器140执行的图像处理的流程图。
14.图7是图6的步骤210的投影处理的流程图。
15.图8是图6的步骤212的vv中心算出处理的流程图。
16.图9是图6的步骤214的特征量算出处理的流程图。
17.图10a是示出彩色眼底图像g1的图。
18.图10b是示出脉络膜血管图像g2的图。
19.图10c是示出检测出四个涡静脉vv1～vv4的状况的图。
20.图11是示出了向假想球面(假想眼球面)投影了视神经乳头onh
が
经度0度、视神经乳头onh、黄斑m、以及五个vv的位置pvv1～pvv5的眼球模型的图。
21.图12a示出从球体的中心c朝向各vv的位置pvv1～pvv5的矢量vc1～vc5，是倾斜观察眼球模型的图。
22.图12b示出从球体的中心c朝向各vv的位置pvv1～pvv5的矢量vc1～vc5，是从横向观察眼球模型的图。
23.图12c示出从球体的中心c朝向各vv的位置pvv1～pvv5的矢量vc1～vc5，是从上方观察眼球模型的图。
24.图13a示出通过将朝向各vv的矢量vc1～vc5合成而算出的合成矢量vct，是倾斜观察眼球模型的图。
25.图13b示出通过将朝向各vv的矢量vc1～vc5合成而算出的合成矢量vct，是从横向观察眼球模型的图。
26.图13c示出通过将朝向各vv的矢量vc1～vc5合成而算出的合成矢量vct，是从上方观察眼球模型的图。
27.图14a示出以使合成矢量vct的长度成为矢量vc1的长度(例如，1)的方式归一化而得到的归一化合成矢量vctn，是倾斜观察眼球模型的图。
28.图14b示出以使合成矢量vct的长度成为矢量vc1的长度(例如，1)的方式归一化而得到的归一化合成矢量vctn，是从横向观察眼球模型的图。
29.图14c示出以使合成矢量vct的长度成为矢量vc1的长度(例如，1)的方式归一化而得到的归一化合成矢量vctn，是从上方观察眼球模型的图。
30.图15是说明表示视神经乳头onh与vv中心的位置关系的特征量的概念的图。
31.图16是示出存储表示视神经乳头onh与vv中心的位置关系的特征量的特征量表的图。
32.图17是示出用于显示对象以及特征量的第1显示画面300a的图。
33.图18是示出用于显示对象以及特征量的第2显示画面300b的图。
具体实施方式
34.以下，参照附图，详细说明本公开的技术的实施方式。
35.参照图1说明眼科系统100的构成。如图1所示，眼科系统100具备眼科装置110、眼轴长度测定器120、管理服务器装置(以下，称为“服务器”)140、以及图像显示装置(以下，称为“观察器”)150。眼科装置110获取眼底图像。眼轴长度测定器120测定患者的眼轴长度。服务器140将利用眼科装置110拍摄患者的眼底而得到的眼底图像与患者的id对应地进行存储。观察器150显示从服务器140获取的眼底图像等的医疗信息。
36.眼科装置110、眼轴长度测定器120、服务器140、以及150经由网络130相互连接在一起。
37.接下来，参照图2说明眼科装置110的构成。
38.为了便于说明，将扫描型激光检眼镜(scanning laser ophthalmoscope)称为“slo”。另外，将光学相干断层扫描仪(optical coherence tomography)称为“oct”。
39.此外，将在眼科装置110在水平面设置的情况下的水平方向设为“x方向”，将相对于水平面的垂直方向设为“y方向”，将连结被检查眼12的前眼部的瞳孔的中心和眼球的中心的方向设为“z方向”。因此，x方向、y方向、以及z方向彼此垂直。
40.眼科装置110包括拍摄装置14以及控制装置16。拍摄装置14具备slo单元18、oct单元20、以及拍摄光学系统19，获取被检查眼12的眼底的眼底图像。以下，将利用slo单元18获
取的二维眼底图像称为slo图像。另外，将基于利用oct单元20获取的oct数据而创建的视网膜的断层图像或正面图像(en-face图像)等称为oct图像。
41.控制装置16具备具有cpu(centralprocessingunit：中央处理装置)16a、ram(random access memory：随机存储器)16b、rom(read-only memory：只读存储器)16c、以及输入输出(i/o)端口16d的计算机。
42.控制装置16具备经由i/o端口16d与cpu16a连接的输入/显示装置16e。输入/显示装置16e具有显示被检查眼12的图像、或者从用户受理各种指示的图形用户界面。作为图形用户界面而举出触控面板、显示器。
43.另外，控制装置16具备连接于i/o端口16d的图像处理器16g。图像处理器16g基于通过拍摄装置14得到的数据生成被检查眼12的图像。此外，也可以省略图像处理器16g，而由cpu16a基于通过拍摄装置14得到的数据生成被检查眼12的图像。控制装置16具备连接于i/o端口16d的通信界面(i/f)16f。眼科装置110经由通信界面(i/f)16f以及网络130与眼轴长度测定器120、服务器140、以及观察器150连接。
44.像上述那样，在图2中，眼科装置110的控制装置16具备输入/显示装置16e，但本公开的技术不限于此。例如，可以设为眼科装置110的控制装置16不具备输入/显示装置16e，而具备与眼科装置110物理上独立的其他的输入/显示装置。在该情况下，该显示装置具备在控制装置16的cpu16a的控制下动作的图像处理处理器单元。也可以设为图像处理处理器单元基于由cpu16a输出指示的图像信号来显示slo图像等。
45.拍摄装置14在控制装置16的cpu16a的控制下工作。拍摄装置14包括slo单元18、拍摄光学系统19、以及oct单元20。拍摄光学系统19包括第1光学扫描仪22、第2光学扫描仪24、以及广角光学系统30。
46.第1光学扫描仪22沿x方向、以及y方向二维扫描从slo单元18射出的光。第2光学扫描仪24沿x方向、以及y方向二维扫描从oct单元20射出的光。第1光学扫描仪22以及第2光学扫描仪24只要为能够使光束偏转的光学元件即可，例如，能够使用多面镜或电流镜等。另外，也可以为这些的组合。
47.广角光学系统30包括具有共用光学系统28的对物光学系统(在图2中未图示)、以及将来自slo单元18的光和来自oct单元20的光合成的合成部26。
48.此外，共用光学系统28的对物光学系统也可以为使用椭圆镜等凹面镜的反射光学系统、或使用广角透镜等的折射光学系统、或者将凹面镜或透镜组合的反射折射光学系统。通过使用利用了椭圆镜或广角透镜等的广角光学系统，不仅能够拍摄存在视神经乳头或黄斑的眼底中心部，还能够拍摄存在眼球的赤道部或涡静脉的眼底周边部的视网膜。
49.在使用包括椭圆镜在内的系统的情况下，也可以使用利用在国际公开wo2016/103484或者国际公开wo2016/103489中记载的椭圆镜的系统的构成。国际公开wo2016/103484的公开以及国际公开wo2016/103489的各公开整体通过参照而援引至本说明书。
50.利用广角光学系统30，在眼底实现基于广视角(fov：field of view)12a的观察。fov12a示出了能够利用拍摄装置14拍摄的范围。fov12a可表现为视角。视角在本实施方式中可利用内部照射角和外部照射角来规定。外部照射角是指，针对从眼科装置110向被检查眼12照射的光束的照射角而以瞳孔27为基准规定的照射角。另外，内部照射角是指，针对向眼底照射的光束的照射角而以眼球中心o为基准规定的照射角。外部照射角与内部照射角
处于对应关系。例如，在外部照射角为120度的情况下，内部照射角大约相当于160度。在本实施方式中，内部照射角设为200度。
51.内部照射角的200度为本公开的技术的“规定值”的一例。
52.在此，将以内部照射角且为160度以上的拍摄视角拍摄而得到的slo眼底图像称为uwf-slo眼底图像。此外，uwf是指超广角(ultra wide field)的简称。当然，通过以内部照射角且以小于160度的拍摄视角拍摄眼底，能够获取不为uwf的slo图像。
53.slo系统通过图2示出的控制装置16、slo单元18、以及摄影光学系统19来实现。slo系统具备广角光学系统30，因此，能够实现基于广fov12a的眼底摄影，具体来说，能够拍摄从被检查眼12的眼底的后极部超过赤道部的区域。
54.使用图3a说明赤道部174。眼球(被检查眼12)为设为直径大约为24mm的眼球中心170的球状的构造物。将该前极175和后极176相连的直线称为眼球轴172，将与眼球轴172正交的平面相交于眼球表面的线称为纬线，纬线长度最大的纬线为赤道174。将与赤道174的位置相当的视网膜或脉络膜的部分称为赤道部178。赤道部178为眼底周边部的一部分。
55.眼科装置110能够将被检查眼12的眼球中心170设为基准位置来拍摄内部照射角为200
°
的区域。此外，在200
°
的内部照射角的情况下，以被检查眼12的眼球的瞳孔为基准的外部照射角为167
°
。也就是说，广角光学系统80以外部照射角167
°
的视角从瞳孔照射激光，以内部照射角拍摄200
°
的眼底区域。
56.在图3b示出了利用能够以200
°
的内部照射角扫描的眼科装置110拍摄而得到的slo图像179。如图3b所示，赤道部174相当于内部照射角为180
°
，在slo图像179中，用虚线178a表示的部位相当于赤道部178。像这样，眼科装置110能够一并(一次拍摄、或者一次扫描)拍摄从包括后极176的后极部超过了赤道部178的眼底周边区域。也就是说，眼科装置110能够一并拍摄从眼底中心部到眼底周边部。
57.详细在后面进行说明，图3c是示出眼球中的脉络膜12m与涡静脉(vortex vein)12v1、v2的位置关系的图。
58.在图3c中，网格状的图案示出了脉络膜12m的脉络膜血管。脉络膜血管使血液在脉络膜整体中循环。然后，血液从在被检查眼12存在多个的涡静脉流向眼球之外。在图3c中，示出了在眼球的单侧存在的上侧涡静脉v1和下侧涡静脉v2。涡静脉大多存在于赤道部178附近。因此，为了拍摄在被检查眼12存在的涡静脉以及涡静脉周边的脉络膜血管，使用上述的能够以200
°
的内部照射角对眼底周边部进行宽范围的扫描的眼科装置110来进行。
59.作为具备广角光学系统30的眼科装置110的构成，可以使用在国际申请pct/ep2017/075852中记载的构成。通过参照将在2017年10月10日提交国际申请的国际申请pct/ep2017/075852(国际公开wo2018/069346)的公开整体援引至本说明书。
60.返回至图2，slo单元18具备多个光源、例如，b光(蓝色光)的光源40、g光(绿色光)的光源42、r光(红色光)的光源44、以及ir光(红外线(例如，近红外光))的光源46、以及对来自光源40、42、44、46的光进行反射或者透射而导入一个光路的光学系统48、50、52、54、56。光学系统48、50、56为反射/透射镜，光学系统52、54为分束器。b光被光学系统48反射，将光学系统50透射，被光学系统54反射，g光被光学系统50、54反射，r光将光学系统52、54透射，ir光被光学系统56、52反射，分别导向一个光路。
61.slo单元18构成为能够将发出g光、r光、以及b光的模式、和发出红外线的模式等发
出波长不同的激光的光源或者发光的光源的组合切换。在图2示出的例子中，具备b光(蓝色光)的光源40、g光的光源42、r光的光源44、以及ir光的光源46这四个光源，但本公开的技术不限于此。例如，slo单元18还可以构成为还具备白色光的光源，利用仅发出白色光的模式等各种模式发出光。
62.从slo单元18入射至拍摄光学系统19的光利用第1光学扫描仪22而沿x方向以及y方向扫描。扫描光经由广角光学系统30以及瞳孔27而照射至被检查眼12的后眼部。被眼底反射的反射光经由广角光学系统30以及第1光学扫描仪22而向slo单元18入射。
63.slo单元18具备来自被检查眼12的后眼部(例如，眼底)的光中的、反射b光且透射b光以外的分束器64、透射了分束器64的光中的、反射g光且透射g光以外的分束器58。slo单元18具备将透射了分束器58的光中的、r光反射且且将r光以外的光透射的分束器60。slo单元18具备将透射了分束器60的光中的、ir光反射的分束器62。
64.slo单元18与多个光源对应地具备多个光检测元件。slo单元18具备检测由分束器64反射的b光的b光检测元件70、以及检测由分束器58反射的g光的g光检测元件72。slo单元18具备检测由分束器60反射的r光的r光检测元件74、以及检测由分束器62反射的ir光的ir光检测元件76。
65.在经由广角光学系统30以及第1光学扫描仪22向slo单元18入射的光(也就是说，由眼底反射的反射光)为b光的情况下，由分束器64反射且被b光检测元件70接收，在为g光的情况下，透射分束器64、由分束器58反射且由g光检测元件72接收。在上述入射的光为r光的情况下，透射分束器64、58，由分束器60反射且由r光检测元件74接收。在上述入射的光为ir光的情况下，透射分束器64、58、60，由分束器62反射且被ir光检测元件76接收。在cpu16a的控制下动作的图像处理器16g使用利用b光检测元件70、g光检测元件72、r光检测元件74、以及ir光检测元件76检测出的信号，生成uwf-slo图像。
66.uwf-slo图像(如后述那样也称为uwf眼底图像、原始图像)存在用g色拍摄眼底得到的uwf-slo图像(g色眼底图像)、以及用r色拍摄眼底得到的uwf-slo图像(r色眼底图像)。uwf-slo图像存在用b色拍摄眼底得到的uwf-slo图像(b色眼底图像)、以及用ir色拍摄眼底得到的uwf-slo图像(ir眼底图像)。
67.另外，控制装置16以同时发光的方式控制光源40、42、44。通过以b光、g光以及r光同时拍摄被检查眼12的眼底，得到各位置彼此对应的g色眼底图像、r色眼底图像、以及b色眼底图像。从g色眼底图像、r色眼底图像、以及b色眼底图像得到rgb彩色眼底图像。控制装置16以同时发光的方式控制光源42、44，以g光以及r光同时拍摄被检查眼12的眼底，由此，得到各位置彼此对应的g色眼底图像以及r色眼底图像。从g色眼底图像以及r色眼底图像得到rg彩色眼底图像。
68.像这样，作为uwf-slo图像，具体来说，存在b色眼底图像、g色眼底图像、r色眼底图像、ir眼底图像、rgb彩色眼底图像、rg彩色眼底图像。uwf-slo图像的各图像数据与经由输入/显示装置16e而输入的患者的信息一并经由通信界面(i/f)16f从眼科装置110向服务器140发送。uwf-slo图像的各图像数据和患者的信息对应地存储在存储器164内。此外，患者的信息例如存在患者名id、姓名、年龄、视力、右眼/左眼的区别等。操作员经由输入/显示装置16e输入患者的信息。
69.oct系统利用图2示出的控制装置16、oct单元20、以及拍摄光学系统19来实现。oct
系统具备广角光学系统30，因此，与上述对slo眼底图像的拍摄同样地，能够进行以广fov12a的眼底摄影。oct单元20包括光源20a、传感器(检测元件)20b、第一光耦合器20c、参照光学系统20d、准直透镜20e、以及第二光耦合器20f。
70.从光源20a射出的光由第一光耦合器20c分支。分支出的一方的光被作为测定光用准直透镜20e准直成平行光之后，入射至拍摄光学系统19。利用第2光学扫描仪24沿x方向以及y方向扫描测定光。扫描光经由广角光学系统30以及瞳孔27而照射至眼底。由眼底反射的测定光经由广角光学系统30以及第2光学扫描仪24而向oct单元20入射，经由准直透镜20e以及第一光耦合器20c而入射至第二光耦合器20f。
71.从光源20a射出且由第一光耦合器20c分支的另一方的光作为参照光向参照光学系统20d入射，经由参照光学系统20d而入射至第二光耦合器20f。
72.入射至第二光耦合器20f的这些光、即、由眼底反射的测定光、以及参照光被第二光耦合器20f干涉而生成干涉光。干涉光被传感器20b接收。在cpu16a的控制下动作的图像处理器16g基于由传感器20b检测出的oct数据而生成断层图像或en-face图像等的oct图像。此外，也可以省略图像处理器16g，由cpu16a基于由传感器20b检测出的oct数据生成oct图像。
73.在此，将以内部照射角为160度以上的拍摄视角拍摄得到的oct眼底图像称为uwf-oct图像。当然，能够以内部照射角少于160度的拍摄视角获取oct眼底图像数据。
74.uwf-oct图像的图像数据与患者的信息一并经由通信界面(i/f)16f从眼科装置110向服务器140发送。uwf-oct图像的图像数据与患者的信息对应地存储在存储器164内。
75.此外，在本实施方式中，例示光源20a为波长扫描类型的ss-oct(swept-source oct：扫频光源oct)，但也可以为sd-oct(spectral-domain oct：断层扫描oct)、td-oct(time-domain oct：时域oct)等各种方式的oct系统。
76.接下来，说明眼轴长度测定器120。眼轴长度测定器120具备测定被检查眼12的眼轴方向上的长度即眼轴长度的第一模式和第二模式这两个模式。第一模式在将来自未图示的光源的光向被检查眼12导光之后，接受来自眼底的反射光和来自眼角膜的反射光的干涉光，基于表示接受的干涉光的干涉信号来测定眼轴长度。第二模式为使用未图示的超声波来测定眼轴长度的模式。
77.眼轴长度测定器120将利用第一模式或者第二模式测定出的眼轴长度发送至服务器140。可以利用第一模式以及第二模式测定眼轴长度，在该情况下，将利用两个模式测定出的眼轴长度的平均值作为眼轴长度而发送至服务器140。服务器140将患者的眼轴长度与患者名id对应地进行存储。
78.接下来，参照图4说明服务器140的构成。如图4所示，服务器140具备控制单元160、以及显示/操作单元170。控制单元160具备包括cpu162的计算机、存储装置即存储器164、以及通信界面(i/f)166等。此外，在存储器164存储有图像处理程序。显示/操作单元170具备显示器172以及触控面板等的输入/指示设备174，该输入/指示设备174显示图像、或者受理各种指示的图形用户界面。
79.cpu162为本公开的技术的“处理器”的一例。存储器164为本公开的技术的“计算机可读存储介质”的一例。控制单元160为本公开的技术的“计算机程序产品”的一例。服务器140为本公开的技术的“图像处理装置”的一例。
80.观察器150的构成与服务器140相同，因此，省略对其进行说明。
81.接下来，参照图5，说明通过服务器140的cpu162执行图像处理程序而实现的各种功能。图像处理程序具备图像处理功能、显示控制功能、以及处理功能。cpu162执行具备各功能的图像处理程序，由此，cpu162如图4所示，作为图像处理部182、显示控制部184、以及处理部186发挥作用。
82.图像处理部182为本公开的技术的“检测部”的一例。
83.图像处理部182为本公开的技术的“算出部”、“投影部”、以及“生成部”的一例。
84.接下来，使用图6详细说明服务器140的图像处理。通过服务器140的cpu162执行图像处理程序，实现图6的流程图示出的图像处理。
85.在服务器140接收用眼科装置110拍摄到的眼底图像(例如，uwf眼底图像)的图像数据，基于接收到的眼底图像中的第1眼底图像(r色眼底图像)和第2眼底图像(g色眼底图像)的图像数据而生成脉络膜血管图像时，执行图像处理程序。图10a示出了接收的眼底图像中的rgb彩色眼底图像g1。此外，眼科装置110向服务器140一并发送图像数据和患者的信息(例如，患者名、患者id、年龄、以及视力)。
86.脉络膜血管图像以下述方式来生成。
87.首先，说明在第1眼底图像(r色眼底图像)和第2眼底图像(g色眼底图像)包含的信息。
88.眼的构造成为由构造不同的多个层覆盖玻璃体。多个层从玻璃体侧的最内侧向外侧包括视网膜、脉络膜、巩膜。r光从视网膜通过而到达脉络膜。因此，第1眼底图像(r色眼底图像)包括存在于视网膜的血管(视网膜血管)的信息和存在于脉络膜的血管(脉络膜血管)的信息。与之相对地，g光仅到达视网膜。因此，第2眼底图像(g色眼底图像)包含存在于视网膜的血管(视网膜血管)的信息。因此，从第2眼底图像(g色眼底图像)提取视网膜血管，从第1眼底图像(r色眼底图像)除去视网膜血管，由此，能够得到脉络膜血管图像。
89.接下来，说明脉络膜血管图像的生成方法。服务器140的图像处理部182通过对第2眼底图像(g色眼底图像)实施黑帽滤波处理，从第2眼底图像(g色眼底图像)提取视网膜血管。接下来，图像处理部182使用从第2眼底图像(g色眼底图像)提取的视网膜血管来进行描绘处理，从第1眼底图像(r色眼底图像)除去视网膜血管。也就是说，使用从第2眼底图像(g色眼底图像)提取的视网膜血管的位置信息进行对第1眼底图像(r色眼底图像)的视网膜血管构造填涂成与周围的像素相同的值的处理。然后，图像处理部182相对于被除去了视网膜血管的第1眼底图像(r色眼底图像)的图像数据适当实施直方图均等化处理(clahe、contrast limited adaptive histogram equalization)，由此，在第1眼底图像(r色眼底图像)中增强脉络膜血管。由此，得到图10b示出的脉络膜血管图像g2。生成的脉络膜血管图像与患者的信息对应地存储在存储器164内。
90.另外，从第1眼底图像(r色眼底图像)和第2眼底图像(g色眼底图像)生成脉络膜血管图像，但接下来，图像处理部182也可以使用第1眼底图像(r色眼底图像)或者用ir光拍摄到的ir眼底图像来生成脉络膜血管图像。就生成脉络膜眼底图像的方法而言，通过参照在2018年3月20日提交的日本特愿2018-052246以及国际公开wo2019/181981a1的公开的全部内容而援引至本说明书。
91.当图像处理程序开始时，在图5的步骤202中，处理部186从存储器164读出脉络膜
血管图像g2(参照图10b)和g色眼底图像，来作为眼底图像。
92.在步骤204中，图像处理部182从g色眼底图像推定黄斑的位置。具体来说，由于黄斑在g色眼底图像中为暗区域，所以图像处理部182检测出在上述读出的g色眼底图像中像素值最小的规定数量的像素的区域来作为黄斑的位置。表示黄斑的位置的座标与患者的信息对应地存储在存储器164内。
93.在步骤206中，图像处理部182从g色眼底图像检测视神经乳头onh的位置。具体来说，图像处理部182通过相对上述读出的g色眼底图像进行事先规定的视神经乳头的图像的图案匹配，而在g色眼底图像中检测视神经乳头onh的位置。表示视神经乳头onh的位置的座标与患者的信息对应地存储在存储器164内。
94.此外，由于视神经乳头在g色眼底图像中为最亮的区域，所以也可以设为检测出在上述读出的g色眼底图像中像素值最大的规定数量的像素的区域来作为视神经乳头的位置。
95.另外，脉络膜血管图像通过以上述方式处理r色眼底图像以及g色眼底图像来创建。因此，当将g眼底图像的座标系与脉络膜血管图像的座标系重叠时，g眼底图像的座标系的各位置与脉络膜血管图像的座标系的各位置相同。因此，从g色眼底图像检测出的相当于黄斑以及视神经乳头各自的位置的脉络膜血管图像上的各位置为黄斑以及视神经乳头各自的位置。
96.因此，在步骤204的处理中，可以取代g色眼底图像而从脉络膜血管图像检测黄斑的位置。同样地，在步骤206的处理中，可以取代g色眼底图像而从脉络膜眼底图像检测视神经乳头的位置。
97.在步骤208中，图像处理部182在脉络膜血管图像中检测涡静脉(vortex vein，以下称为“vv”))位置。在此，涡静脉vv是指向脉络膜流入的血流的流出路径，在靠近球的赤道部的后极有多个。
98.图像处理部182求出脉络膜血管图像中的各像素的血管走向(blood vesselrunning direction)。具体来说，图像处理部182相对于所有像素反复进行下述处理。即，图像处理部182设定由以像素为中心的周围的多个像素构成的区域(单元)。然后，基于计算对象像素的周围的像素的亮度值，来计算单元的各像素中的亮度的斜率方向(例如，用0度以上～低于180度的角度来表示。此外，0度被定义为直线(也就是说，水平线)的方向)。相对于单元内的所有像素进行该斜率方向的计算。
99.接下来，创建有斜率方向为0度、20度、40度、60度、80度、100度、120度、140度、160度的9个竖条区域(bin)(例如，各竖条区域的宽度为20度)的直方图，因此，对与各竖条区域对应的斜率方向的单元内的像素数进行计数。在直方图的一个竖条区域的宽度相当于20度，0度的竖条区域设定有具有0度以上低于10度和170度以上低于180度的斜率方向的单元内的像素数(也就是说，计数值)。20度的竖条区域被设定有具有10度以上低于30度的斜率方向的单元内的像素数(也就是说，计数值)。同样地，也被设定40度、60度、80度、100度、120度、140度、160度的竖条区域的计数值。由于直方图的竖条区域的数量为9，所以像素的血管走向用9种方向的某一种来定义。此外，通过使竖条区域的宽度变窄并增多竖条区域的数量，能够提高血管走向的分辨率。各竖条区域中的计数值(也就是说，直方图的纵轴)被标准化，创建了相对于解析点的直方图。
100.接下来，图像处理部182从直方图确定解析点的血管走向。具体来说，确定计数值最小的角度(例如，设为60度)的竖条区域，并将确定出的竖条区域的斜率方向即60度确定为像素的血管走向。此外，计数最少的斜率方向为血管走向是出于以下的理由。在血管走向上亮度斜率小，另一方面，在除此以外的方向上亮度斜率大(例如，血管与血管以外的部分的亮度之差很大)。因此，若创建各像素的亮度斜率的直方图，则相对于血管走向的竖条区域的计数值变少。同样地，相对于脉络膜血管图像中的各像素创建直方图，算出各像素的血管走向。算出的各像素的血管走向与患者的信息对应地存储在存储器164内。
101.然后，图像处理部182在脉络膜血管图像上等间隔地设定在纵向上有m个、在横方向上有n个的合计l个假想的颗粒的初始位置。例如，设定m＝10、n＝50、合计l＝500个的初始位置。
102.图像处理部182还获取最初位置(l个中的某一个)的血管走向，使假想的颗粒沿着获取出的血管走向仅移动规定距离，在移动后的位置再次获取血管走向，并使假想的颗粒沿着获取的血管走向仅移动规定距离。像这样，以事先设定了移动次数重复执行沿着血管走向移动规定距离移动。在l个位置全部执行以上的处理。在以设定的移动次数使l个所有假想的颗粒进行了移动后的时间点，将聚集了一定个数以上假想的颗粒的点检测为vv位置。图10c示出了检测出四个涡静脉vv1～vv4的状况。
103.然后，将检测出的vv的个数、vv位置信息(例如，表示脉络膜血管图像中的vv位置的座标)与患者的信息对应地存储在存储器164内。不仅可以使用脉络膜血管图像，也可以使用各种眼底图像来进行涡静脉vv的检测。可以为用各种可视光的波长拍摄到的r色眼底图像或g色眼底图像等的彩色眼底图像、基于荧光拍摄的荧光眼底图像、从三维的oct体积数据生成的二维图像(en face oct图像)、进一步对这些图像进行图像处理而得到的2值化图像、对血管进行增强处理而得到的血管图像。还可以为通过从血管图像提取脉络膜血管的处理而生成的脉络膜血管图像。
104.在步骤210中，图像处理部182将视神经乳头onh、黄斑m、各vv投影至眼球模型的球面。
105.图7示出了图6的步骤210的投影处理的流程图。在图7的步骤221中，处理部186从存储器164获取与患者的信息对应的脉络膜眼底图像中的视神经乳头onh、黄斑m、vv位置的眼底图像(uwf眼底图像)上的各座标(x，y)。
106.在步骤223中，图像处理部182进行将视神经乳头onh、黄斑m、vv的位置的各座标投影至与图11示出的眼球模型的眼底对应的假想球面(眼球面或者眼球模型的表面)的处理。该眼球模型为将眼球的中心设为c、将半径设为r(将眼轴长度设为2r)的球体模型。将该眼球模型的球面定义为眼球面。眼球面上的点利用纬度和经度来确定位置。
107.以使视神经乳头onh作为基准点(经度0度、纬度0度)、且使黄斑m的纬度成为0度的方式分别进行投影。各个vv位置也以视神经乳头onh为基准点在眼球面上投影在pvv1、pvv2、pvv3、pvv4、pvv5(在被检查眼检测到五个涡静脉vv的情况)。
108.因此，视神经乳头onh和黄斑m分别投影至0度的纬线，在球面上将视神经乳头onh和黄斑m连结的线段位于0度的纬线上。
109.在此，眼球模型也由三维空间(x，y，z)规定，由图像处理的存储器空间定义。像上述那样，x方向为在眼科装置110设在水平面的情况下的水平方向，患者位于面向眼科装置
110的位置，因此，x方向为被检查眼的左右方向。y方向为相对于水平面的垂直方向，为被检查眼的上下方向。z方向为将被检查眼12的前眼部的瞳孔的中心和眼球的中心连结的方向。
110.在存储器164中存储了用于对二维的眼底图像(uwf眼底图像)上的位置的各座标(x，y)向三维空间(x，y，z)进行反向立体投影转换的转换式。在本实施方式中，利用该转换式，以使视神经乳头onh的位置在眼球面上的基准点(经度以及纬度各自为0度的位置)、且使黄斑m的位置成为纬度为0度的位置的方式进行转换(也就是说，投影)。
111.还使用该转换式将二维的眼底图像上的各vv的位置(x、y)转换(也就是说，投影)为三维的眼球面的位置pvv1～pvv5。
112.眼球面上的涡静脉vv1～vv5的位置pvv1～pvv5的座标(x，y，z)和视神经乳头onh、黄斑m的三维座标一并与患者的识别信息建立关联地存储在存储器164内。
113.当图7的步骤223的处理结束时，图像处理前进至图6的步骤212。
114.在步骤212中，图像处理部182算出多个涡静脉的分布中心点(以下，设为vv中心)。vv中心是指多个vv的各个位置的分布中心点。分布中心点能够通过首先求出从眼球模型的中心点c朝向涡静脉位置(pvv1等)的矢量、然后求出将朝向各vv的多个矢量合成的合成矢量来求得。
115.图8示出了图6的步骤212的vv中心的算出处理的流程图。在图8中，进行在被检查眼检测到五个涡静脉vv的说明。
116.在图8的步骤232中，处理部186与患者的信息对应地从存储器164获取各vv的位置pvv1～pvv5的座标(x1，y1，z1)～(x5，y5，z5)。在此，眼球模型使用将半径设为单位长度1的单位球。
117.在步骤234中，图像处理部182如图12a～图12c所示，如数学式1那样，生成从眼球模型的球体的中心c朝向各vv的位置pvv1～pvv5的矢量vc1～vc5。
118.【数学式1】
[0119][0120][0121][0122][0123][0124]
各矢量vc的大小为单位长度1。图12a～图12c分别是从斜方、横向、上方观察单位球的图。
[0125]
在步骤236中，图像处理部182如图13a～图13c所示，通过将朝向各vv的矢量vc1～vc5合成来算出合成矢量vct(参照数学式2的vvcomposition)。具体来说，如数学式2那样，为求各矢量vc的座标值之和的处理。
[0126]
【数学式2】
[0127][0128]
从中心c利用合成矢量vct确定的点pvvc的座标用位置pvv1～pvv5的各个x、y、z的各个座标之和来确定。点pvvc为vv1～vv5的各个位置的三维空间中的分布中心。点pvvc为vv1～vv5的各个位置之间的距离(例如，欧几里德距离)的合计成为最小的位置。点pvvc通常位于眼球模型的球体内部。图13a～图13c分别是从斜方、横向、上方观察作为单位球的眼球模型的图。
[0129]
在步骤238中，图像处理部182如图14a～图14c所示，以使合成矢量的大小成为单位大小1的方式将合成矢量vct归一化。由此，如数学式3那样得到归一化的合成矢量vctn(参照数学式3的vvcenter)。
[0130]
【数学式3】
[0131][0132]
然后，在步骤240中从中心c利用归一化合成矢量vctn确定的点pvvcn位于作为单位球的眼球模型的眼球面上。点pvvcn为归一化合成矢量vctn与单位球的眼球面的交点。点pvvcn为多个vv的各个位置的眼球模型的球面上的分布中心。
[0133]
另外，还能够定义多个vv各自的位置与分布中心之间的距离(例如，眼球模型的球面上的大圆距离)的合计成为最小的位置。图14a～图14c分别是斜方、横向、上方观察单位球的眼球模型的图。
[0134]
在步骤240中，求出多个涡静脉中的、在单位球的球面上的分布中心位置。因此，判断涡静脉的分布中心的纬度和经度。
[0135]
在步骤241中，图像处理部182以使单位球成为眼轴长度2r的眼球模型的方式放大。也就是说，若被检查眼的眼轴长度为24mm(即，2＊12mm、r＝12)，则以成为半径为12mm的球的方式进行转换。
[0136]
在步骤241中，图像处理部182还提取将反映了被检查眼的眼轴长度的眼球模型中的视神经乳头onh、黄斑m、涡静脉vv1～vv5、以及涡静脉的分布中心pvvcn各自在眼球面上的座标、或者以视神经乳头作为基准算出的黄斑m、涡静脉、涡静脉的分布中心各自的纬度以及经度。提取的座标或者纬度以及经度的信息与患者的信息对应地存储在存储器164内。也可以将座标和纬度以及经度双方保存在存储器164内。
[0137]
当图8的步骤241的处理结束时，图像处理前进至图6的步骤214。
[0138]
在步骤214中，图像处理部182算出表示眼球模型的特定部位的位置与包含vv中心在内的多个对象的位置之间的位置关系的至少一个特征量。作为特定部位，例如有视神经乳头、黄斑、瞳孔中心(也就是说，眼角膜的顶点)、中央凹等。以下，作为特定部位，对视神经乳头进行说明。
[0139]
图9示出了图6的步骤214的特征量算出处理的详细流程图。在图9的步骤242中，图像处理部182将识别特征量算出的对象的变量t设为0，在步骤244中，图像处理部182对变量t增加1。例如，利用变量t＝1、2、3
……
来识别vv中心、vv1、vv2
……
。
[0140]
在步骤246中，图像处理部182以视神经乳头onh为基准点(例如，经度：0度、纬度：0
度)，算出到对象t为止的大圆距离和角度、对象t所处的经度和纬度、以及从基准点到对象t为止的经度距离和纬度距离来作为对象t的特征量。具体来说，如下。
[0141]
首先，说明变量t＝1的情况，即、说明以vv中心为对象算出特征量的情况。
[0142]
图像处理部182如图15所示，根据球面三角法的公式算出视神经乳头onh与vv中心(点pvvcn)的大圆距离gcl。将在以从眼球模型的球体的中心c通过的方式切断球体时的切口定位为大圆，大圆距离是指，将球面上的距离计测对象的2地点(视神经乳头onh和vv中心(点pvvcn))连结的大圆的弧的长度。
[0143]
接下来，图像处理部182如图15所示，根据正形投影或球面三角法算出将黄斑m的位置和视神经乳头onh的位置连结的大圆上的第1线段gcl与将视神经乳头onh的位置和vv中心pvvcn的位置连结的第2线段gl所成角度θ。
[0144]
上述计算距离以及角度的方法与在2019年04月18日提交国际申请的国际申请第pct/jp2019/016653中记载的方法相同。通过参照在2019年04月18日提交国际申请的国际申请第pct/jp2019/016653(国际公开第wo2019203310号(在2019年10月24日提交的国际公开))中记载的上述计算距离以及角度的方法整体而援引至本说明书。
[0145]
接下来，图像处理部182如图15所示算出vv中心pvvcn的经度(lg＝lgn)与纬度(lt＝ltn)。像上述那样，眼球模型由三维空间来规定，但以纬度以及经度分别为0度的位置作为基准点设置视神经乳头ohn的位置，以使黄斑m的位置位于纬度0度的纬线llt上的方式来配置。
[0146]
图像处理部182将vv中心pvvcn的三维空间中的座标(x，y，z)设为视神经乳头onh的位置(x，y，z)＝(1，0，0)，并利用以下这样的规定的转换式进行转换，由此来算出vv中心pvvcn的纬度(lg＝lgn)和经度(lt＝ltn)。
[0147]
经度lt＝arctan(y/x)
[0148]
纬度lg＝arctan(z/√(x2 y2))
[0149]
由此，算出vv中心pvvcn的(lg(纬度)，经度(lt))为(lgn，ltn)。
[0150]
然后，图像处理部182如图15所示根据球面三角法的公式，算出vv中心pvvcn的纬度距离tl、经度距离gl。
[0151]
vv中心pvvcn的纬度距离tl为沿着从vv中心pvvcn通过的经度线llg的中心pvvcn与该经度线llg和纬度0的纬度线llt的交点ptg之间的大圆距离。
[0152]
vv中心pvvcn的经度距离gl为沿着纬度0的纬度线llt的视神经乳头ohn的位置与交点ptg之间的大圆距离。
[0153]
若针对像上述那样用变量t识别出的对象算出上述各特征量，则特征量算出处理前进至步骤248。
[0154]
在步骤248中，图像处理部182判断变量t是否等于算出特征量的对象的总数t。
[0155]
在变量t不等于总数t的情况下，存在没有被算出特征量的对象，因此，特征量算出处理返回至步骤244，执行以上的处理(也就是说，步骤244～248)。
[0156]
在变量t＝2以后的情况下，图像处理部182取代vv中心pvvcn，将vv1以后的位置设为对象，来算出上述距离、角度、纬度、经度、经度距离、以及纬度距离。
[0157]
在变量t等于总数t的情况下，针对应算出特征量的全部对象算出特征量，因此，在步骤250中，处理部186如图16所示与患者的信息对应地将各对象的各特征量存储至存储器
164的特征量表。在图16示出的例子中，特征量表设有患者id、左右眼的识别信息以及与左右眼对应的眼轴长度的项目。然后，针对各被检查眼(左眼或者右眼)设有特征量存储区域。在特征量存储区域中存储了表示上述对象的种类(视神经乳头、黄斑、涡静脉vv1～vv5、vv中心等)的类别信息、和距离、角度、纬度、经度、经度距离、以及纬度距离的特征量信息。此外，视神经乳头onh、黄斑m、涡静脉vv1～vv5、vv中心的眼球面上的座标也可以包含在特征量表内。
[0158]
另外，在上述的例子中，说明了涡静脉的数量为五个，但当然能够应用在被检查眼检测到的涡静脉的个数n(n为自然数)。
[0159]
而且，可以求出所求出的视神经乳头位置与vv中心位置的偏移量、具体来说是视神经乳头位置与vv中心位置之间的距离(大圆距离等)、和由从视神经乳头到vv中心位置的方向规定的特征量(即，将视神经乳头位置和vv中心位置连结的矢量)，来作为该特征量的一部分。而且，同样可以求出黄斑位置与vv中心位置的偏移量、眼底上的病变位置与vv中心位置的偏移量，来作为特征量。该偏移量通常能够设为病况的数值的指标。
[0160]
当步骤250的处理结束时，图6的步骤214的处理结束，因此，图像处理结束。
[0161]
接下来，说明在由服务器140进行多个涡静脉的分布中心的算出处理之后利用观察器150显示分布中心等的处理。
[0162]
在观察器150的显示器的显示画面显示了用于指示生成后述的图像(图17、图18)的图标或按钮。在作为用户的眼科医生对患者进行诊断时，考虑到想要知道vv中心的位置，而点击规定的图标等，从观察器150向服务器140发送与点击的图标等对应的指示信号。
[0163]
接收来自观察器150的指示信号的服务器140生成与指示信号对应的图像(图17、图18)，将所生成的图像的图像数据经由网络130向观察器150发送。从服务器140接收图像数据的观察器150基于接收的图像数据而在显示器显示图像。在服务器140的显示画面的生成处理根据由cpu162(显示控制部184)动作的显示画面生成程序来进行。
[0164]
图17示出了用于显示对象以及特征量的第1显示画面300a。如图17所示，第1显示画面300a具有显示患者的个人信息的个人信息显示栏302、以及图像显示栏320。
[0165]
个人信息显示栏302具有患者id显示栏304、患者姓名显示栏306、年龄显示栏308、眼轴长度显示栏310、视力显示栏312、以及患者选择图标314。在患者id显示栏304、患者姓名显示栏306、年龄显示栏308、眼轴长度显示栏310、以及视力显示栏312显示各信息。此外，当点击患者选择图标314时，在观察器150的显示器172显示患者一览，使用户(眼科医等)选择要解析的患者。
[0166]
图像显示栏320具有拍摄日期显示栏322n1、右眼信息显示栏324r、左眼信息显示栏324l、倾斜观察眼球模型的第1眼球模型图像显示栏326a、从横向观察眼球模型的第2眼球模型图像显示栏328、以及信息显示栏342。在信息显示栏342显示了用户(眼科医等)进行诊察时的评论或笔记来作为文本。
[0167]
在图17示出的例子中，能够显示在拍摄日为2018年3月10日、2017年12月10日、2017年9月10拍摄利用患者id：123456识别出的患者的右眼的眼底(324r点亮)的情况下的眼球模型图像。此外，点击拍摄日期显示栏322n1和右眼信息显示栏324r，显示了在2018年3月10日拍摄得到的右眼的眼球模型图像。
[0168]
如图17所示，在第1眼球模型图像显示栏326a显示倾斜观察到的眼球模型，在第2
眼球模型图像显示栏328显示从横向观察眼球模型的眼球模型。在眼球模型显示各vv1～vv5各自的位置pvv1～pvv5、朝向各vv1～vv5各自的矢量vc1～vc5、从中心c利用归一化合成矢量vctn规定的点pvvcn、归一化合成矢量vctn。
[0169]
当在第1眼球模型图像显示栏326a以及第2眼球模型图像显示栏328的某一栏中点击各vv1～vv5各自的位置pvv1～pvv5、点pvvcn的位置的某一处时，显示与所点击的位置对应的各特征量。例如，当点击点pvvcn的位置时，各特征量显示在点pvvcn的位置附近。具体来说，显示视神经乳头onh与vv中心的大圆距离gc、黄斑m的位置-视神经乳头onh的位置-vv中心pvvcn的位置所成的角度θ、vv中心pvvcn的纬度与经度、以及vv中心pvvcn的纬度距离tl以及经度距离gl。
[0170]
在图18示出了用于显示对象以及特征量的第2显示画面300b。由于第2显示画面300b与第1显示画面300a大致相同，所以对相同的显示部分标记相同的附图标记，省略对其进行说明，仅说明不同的显示部分。此外，第1显示画面300a与第2显示画面300b利用未图示的切换按钮来切换。
[0171]
在第2显示画面300b中，取代第1眼球模型图像显示栏326a以及第2眼球模型图像显示栏328地显示对象位置图像显示栏326b以及对象以及断层图像显示栏328b。
[0172]
在图18示出的例子中，示出了检测出4个vv的例子。
[0173]
像上述那样，眼球模型的眼球上的位置由纬度以及经度来规定，在图6的步骤214的特征量算出处理中，算出了各vv以及vv中心的各个纬度以及经度。在对象位置图像显示栏326b中，基于纬度以及经度的位置关系显示各vv以及vv中心的各个位置。
[0174]
在对象以及断层图像显示栏328b具有原始图像显示区域1032a和vv的断层图像显示区域890。
[0175]
在原始图像显示区域1032a中，与原始图像重叠地显示有oct图像获取区域872a、872b、872c、872d。如上所述，在图6的步骤208中检测出vv位置，在本实施方式中，针对以检测出的vv位置为中心的oct图像获取区域872a、872b、872c、872d获取oct图像。
[0176]
在oct图像显示区域890中显示有与与原始图像显示区域1032a的原始图像重叠显示的oct图像获取区域872a、872b、872c、872d分别对应的oct图像892a、892b、892c、892d。oct图像892a、892b、892c、892d例如显示出通过对oct图像获取区域872a、872b、872c、872d进行扫描而作为断层图像可视化的涡静脉的形状。在oct图像892a、892b、892c、892d的上部标注的no.1、no.2、no.3、no.4的标记与对oct图像获取区域872a、872b、872c、872d标注的no.1、no.2、no.3、no.4的标记分别对应。因此，对象以及断层图像显示栏328b能够将oct图像获取区域872a、872b、872c、872d和表示可视化的涡静脉(或涡静脉和与涡静脉相连的脉络膜血管)的形状的oct图像892a、892b、892c、892d对应关联地进行显示。
[0177]
在oct图像显示区域890显示了表示涡静脉的形状的图像，但本实施方式不限于此。在获取了黄斑的oct图像的情况下，在黄斑的断面图像获取了视神经乳头的oct图像的情况下，视神经乳头的断面图像分别显示在oct图像显示区域890中。
[0178]
在oct图像显示区域890设有未图示的显示切换按钮，当点击显示切换按钮时，显示用于切换oct图像892a、892b、892c、892d的显示的下拉菜单。通过从显示出的下拉菜单选择例如3d多边形、en-face图像、oct像素或血流可视化图像等，能够切换oct图像892a、892b、892c、892d的显示。
[0179]
如以上说明的那样，在本实施方式中，能够得到vv中心，能够支援眼科医生使用vv中心进行的诊断。
[0180]
另外，在上述实施方式中，算出表示眼底的特定部位的位置与分布中心的位置的位置关系的多个特征量并进行显示，因此，能够进一步支援眼科医生使用vv中心进行的诊断。例如，针对特征量算出视神经乳头与vv中心的距离并进行显示，因此，眼科医生能够理解与vv中心的通常的位置(例如，视神经乳头ohn的位置)产生的偏移量(距离)。
[0181]
在以上说明的实施方式中，能够得到纬度以及经度各自为0度的位置为眼底的特定部位的位置的眼球模型的球面中的涡静脉的位置的纬度以及经度、特定部位的位置与涡静脉的位置之间的纬度距离以及经度距离。能够支援眼科医生使用涡静脉的位置的纬度以及经度、特定部位的位置与涡静脉的位置之间的纬度距离以及经度距离进行的诊断。
[0182]
在以上说明的实施方式中，各距离以及各线段为大圆距离，但本公开的技术不限于此。也可以取代大圆距离而算出欧几里德距离(直线距离)，或一并算出大圆距离和欧几里德距离(直线距离)。
[0183]
在上述实施方式中，对通过将朝向各vv的位置的矢量合成而得到的合成矢量进行归一化，利用归一化合成矢量求出vv中心，但本公开的技术不限于此。
[0184]
即，例如，第一，可以在眼球模型中用纬度以及经度特定各vv的位置，求出各vv的位置的纬度以及经度的各个平均值，根据求出的纬度平均值和经度平均值，求出眼球模型的表面的位置(纬度、经度)的vv中心。
[0185]
第二，在眼球模型中，针对预想为vv中心的多个中心候补位置的各点，求出与各vv的位置之间的大圆距离的合计，使用蒙特卡罗方法等来特定大圆距离的合计成为最小的中心候补位置。然后，也可以将该大圆距离的合计成为最小的中心候补位置设为vv中心。
[0186]
第三，为了求出vv中心，也可以不利用检测出的所有vv的位置，而是仅利用在检测出的vv的位置中的位于可能性高的位置的vv的位置。例如，在检测出了10个vv的情况下，并不需要全部利用，例如仅利用9个。选出可能性高的结果为9个，也有选出可能性高的结果为8个的情况或7个的情况。可能性高是指，例如能够通过从对象去掉与相邻的vv之间的距离低于规定值的vv、或者从对象去掉与vv相连的脉络膜血管的根数为规定值以下的vv来选择。在该情况下，也可以算出到10个vv的位置为止的矢量，在合成矢量时，合成到位于可能性高的位置的9个vv的位置为止的矢量，或者仅算出到位于可能性高的位置的9个vv的位置为止的矢量，并将算出的矢量合成。
[0187]
在上述实施方式中，以vv中心作为特征量的算出以及显示的对象，但本公开的技术不限于此。例如，也可以将各vv的眼球模型(三维空间)中的分布中心作为对象。
[0188]
在第一实施方式中，算出了在以视神经乳头为基准点(纬度以及经度各自为0度)的情况下的vv中心的特征量(距离、角度、纬度、经度、经度距离、纬度距离)，但本公开的技术不限于此。例如，也可以算出在以黄斑m为基准点(纬度以及经度各自为0度)时的vv中心的特征量(例如，距离、角度、纬度、经度、经度距离、以及纬度距离)。
[0189]
使用纬度以及经度求出上述各点的位置，但也可以使用极坐标来求出上述各点的位置。
[0190]
在上述实施方式中，说明了利用眼科装置110获取内部光照射角为200度左右的眼底图像的例子。本公开的技术不限于此，也可以为利用内部照射角为100度以下的眼科装置
拍摄到的眼底图像，合成了多个眼底图像的蒙太奇图像也可以应用本公开的技术。
[0191]
在上述实施方式中，利用具备slo拍摄单元的眼科装置110拍摄眼底图像，但也可以为利用能够拍摄脉络膜血管的眼底摄像头拍摄到的眼底图像，利用oct血管造影得到的图像也可以应用本公开的技术。
[0192]
在上述实施方式中，管理服务器140执行图像处理程序。本公开的技术不限于此。例如，也可以由设于眼科装置110、图像观察器150、网络130的其他图像处理装置来执行图像处理程序。在眼科装置110执行图像处理程序的情况下，图像处理程序存储在rom16c内。在图像观察器150执行图像处理程序的情况下，图像处理程序存储在图像观察器150的存储器164内。在其他图像处理装置执行图像处理程序的情况下，图像处理程序存储在其他图像处理装置的存储器内。
[0193]
在上述实施方式中，以具有眼科装置110、眼轴长度测定装置120、管理服务器140、以及图像观察器150的眼科系统100为例子进行了说明，但本公开的技术不限于此。例如，作为第一例，省略眼轴长度测定装置120，但眼科装置110还可以具有眼轴长度测定装置120的功能。另外，作为第二例，眼科装置110还可以具有管理服务器140以及图像观察器150的至少一方的功能。例如，在眼科装置110具有管理服务器140的功能的情况下，能够省略管理服务器140。在该情况下，图像处理程序由眼科装置110或图像观察器150执行。另外，在眼科装置110具有图像观察器150的功能的情况下，能够省略图像观察器150。作为第三例，可以省略管理服务器140，并设为由图像观察器150执行管理服务器140的功能。
[0194]
在本公开中，各构成要素(装置等)只要不产生矛盾，可以仅存在一个也可以存在两个以上。
[0195]
另外，在上述实施方式以及各变形例中，例示了通过利用了计算机的软件构成实现数据处理的情况，但本公开的技术不限于此。例如，可以取代使用了计算机的软件构成，仅通过fpga(field-programmable gate array：现场可编辑逻辑门阵列)或asic(application specific integrated circuit：专用集成电路)等的硬件构成来执行数据处理。可以构成为利用软件构成来执行数据处理中的一部分的处理，并利用硬件执行剩余的处理。
[0196]
像这样，本公开的技术包括通过使用了计算机的软件构成来实现图像处理的情况和并非通过使用了计算机的软件构成来实现图像处理的情况，因此，包括以下的技术。
[0197]
(第一技术)
[0198]
一种图像处理装置，其具备：
[0199]
检测部，其根据被检查眼的眼底图像检测多个涡静脉的位置；以及
[0200]
算出部，其算出检测出的所述多个涡静脉的位置的分布中心。
[0201]
(第二技术)
[0202]
一种图像处理方法，其包括：
[0203]
检测部从被检查眼的眼底图像检测多个涡静脉的位置，
[0204]
算出部算出检测出的所述多个涡静脉的位置的分布中心。
[0205]
(第三技术)
[0206]
一种图像处理装置，其具备：
[0207]
检测部，其从被检查眼的眼底图像检测多个涡静脉的位置、以及视神经乳头的位
置；
[0208]
投影部，其将所述多个涡静脉和所述视神经乳头投影至眼球模型的眼球面；
[0209]
生成部，其生成从所述眼球模型的中心朝向所述眼球面上的各涡静脉的位置的多个矢量；以及
[0210]
算出部，其算出将所述多个矢量合成后的合成矢量。
[0211]
(第四技术)
[0212]
一种图像处理方法，其包括：
[0213]
检测部根据被检查眼的眼底图像检测多个涡静脉的位置、以及视神经乳头的位置，
[0214]
投影部将所述多个涡静脉和所述视神经乳头投影至眼球模型的眼球面，
[0215]
生成部生成从所述眼球模型的中心朝向所述眼球面上的各涡静脉的位置的多个矢量，
[0216]
算出部算出将所述多个矢量合成后的合成矢量。
[0217]
从以上的公开内容提出以下的技术。
[0218]
(第五技术)
[0219]
一种用于进行图像处理的计算机程序产品，
[0220]
所述计算机程序产品具备其自身并非临时性的信号的计算机可读存储介质，
[0221]
在所述计算机可读存储介质保存有程序，
[0222]
所述程序使计算机执行如下的步骤：
[0223]
从被检查眼的眼底图像检测多个涡静脉的位置，
[0224]
算出检测出的所述多个涡静脉的位置的分布中心。
[0225]
(第六技术)
[0226]
一种用于进行图像处理的计算机程序产品，
[0227]
所述计算机程序产品具备其自身并非临时性的信号的计算机可读存储介质，
[0228]
在所述计算机可读存储介质保存有程序，
[0229]
所述程序使计算机执行如下的步骤：
[0230]
从被检查眼的眼底图像检测多个涡静脉的位置、以及视神经乳头的位置，
[0231]
将所述多个涡静脉和所述视神经乳头投影至眼球模型的眼球面，
[0232]
生成从所述眼球模型的中心朝向所述眼球面上的各涡静脉的位置的多个矢量，
[0233]
算出将所述多个矢量合成后的合成矢量。
[0234]
以上说明的各图像处理只不过为一例。因此，当然可以在不脱离主旨的范围内删除不要的步骤、或者追加新的步骤、互换处理顺序。
[0235]
在本说明书中记载的所有文献、专利申请(包括日本技术特愿2019-220285)以及技术规格与具体且分别记载通过参照而援引各个文献、专利申请、以及技术规格的情况同样地，通过参照而援引至本说明书中。