JP2022128882A - Image processing device, image processing program, and image processing method - Google Patents

Image processing device, image processing program, and image processing method Download PDF

Info

Publication number
JP2022128882A
JP2022128882A JP2021027339A JP2021027339A JP2022128882A JP 2022128882 A JP2022128882 A JP 2022128882A JP 2021027339 A JP2021027339 A JP 2021027339A JP 2021027339 A JP2021027339 A JP 2021027339A JP 2022128882 A JP2022128882 A JP 2022128882A
Authority
JP
Japan
Prior art keywords
pixels
pixel
correction
contour
image processing
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.)
Pending
Application number
JP2021027339A
Other languages
Japanese (ja)
Inventor
英治 内野
Eiji Uchino
陸 福田
Riku Fukuda
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.)
Yamaguchi University NUC
Original Assignee
Yamaguchi University NUC
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 Yamaguchi University NUC filed Critical Yamaguchi University NUC
Priority to JP2021027339A priority Critical patent/JP2022128882A/en
Publication of JP2022128882A publication Critical patent/JP2022128882A/en
Pending legal-status Critical Current

Links

Images

Abstract

To correct contour lines of blood vessels extracted from captured images.SOLUTION: An image processing device 1 corrects contour lines of blood vessels generated on the basis of captured images in which the blood vessels are captured. The contour lines include a first contour line and a second contour line. The image processing device comprises: a storage unit 12 which stores the captured image in association with position information of contour pixels corresponding to the contour lines of the blood vessels; a contour line selection unit 131 which selects one of the first contour line and the second contour line as a target contour line of correction; an extraction unit 132 which extracts the plurality of contour pixels as extraction pixels on the basis of the position information of the contour pixels; a pixel selection unit 133 which selects two extraction pixels having a prescribed positional relation as selection pixels; a correction pixel determination unit 134 which determines correction pixels corresponding to the target contour line after the correction on the basis of pixel information on a plurality of reference pixels arranged at prescribed positions from the two selection pixels; and a correction contour line identifying unit 136 which identifies the target contour line after the correction on the basis of the position information of the correction pixels.SELECTED DRAWING: Figure 2

Description

本発明は、画像処理装置と画像処理プログラムと画像処理方法とに関する。 The present invention relates to an image processing device, an image processing program, and an image processing method.

脳内の血管障害の早期発見には、CT(Computed Tomography)やMRI(Magnetic Resonance Imaging)のような人体の内側を可視化する医療機器が使用されている。しかし、これらの医療機器は、高額である。そのため、設備投資の負担が大きく、これらの医療機器を有する病院は限られる。また、医療費も高額になるため、患者の医療費負担も大きい。 Medical equipment for visualizing the inside of the human body, such as CT (Computed Tomography) and MRI (Magnetic Resonance Imaging), is used for early detection of intracerebral vascular disorders. However, these medical devices are expensive. Therefore, the burden of equipment investment is large, and the number of hospitals that have such medical equipment is limited. In addition, since the medical expenses are high, the burden of medical expenses on the patient is also large.

脳は、発生学的に眼と同一の器官である。そのため、脳内の血管の状態は、眼底の網膜血管の状態から推測できる。そこで、脳内の血管障害のリスク診断には、前述の医療機器を必要としない眼底検査が用いられてきている。近年、第3期特定健康診断の改訂により、眼底検査が内科の検査項目に取り入れられている。しかし、眼底検査には判定医の経験と技量とが大きく寄与する。そのため、専門医ではない内科医にとって、精度の良い眼底検査は困難である。したがって、眼底検査の自動化が求められている。眼底検査の自動化には、眼底を撮像した撮像画像から血管を正確に抽出する画像処理技術が必要である。 The brain is an embryologically identical organ to the eye. Therefore, the state of blood vessels in the brain can be estimated from the state of retinal blood vessels in the fundus. Therefore, fundus examination, which does not require the above-mentioned medical equipment, has been used for diagnosing the risk of intracerebral vascular disorders. In recent years, due to the revision of the 3rd specific health checkup, fundus examination has been included in internal medicine examination items. However, the experience and skill of the judgment doctor greatly contributes to fundus examination. Therefore, it is difficult for physicians who are not specialists to perform an accurate fundus examination. Therefore, there is a demand for automation of fundus examination. Automated fundus examination requires an image processing technique for accurately extracting blood vessels from captured images of the fundus.

これまでにも、撮像画像から血管を抽出する画像処理技術が提案されている(例えば、非特許文献1参照)。 Image processing techniques for extracting blood vessels from captured images have also been proposed (see, for example, Non-Patent Document 1).

非特許文献1に開示されている技術は、エネルギー関数を用いた探索により、血管壁の解候補を選出し、血管壁との類似度が最大の解候補を血管輪郭線として抽出する。しかしながら、眼底画像の視認性は悪く、非特許文献1に開示されている技術による血管輪郭線の抽出精度は必ずしも良好とは言えない。そのため、撮像画像から抽出された血管の輪郭線を補正する技術が求められている。 The technique disclosed in Non-Patent Document 1 selects solution candidates for the blood vessel wall by searching using an energy function, and extracts the solution candidate with the highest degree of similarity to the blood vessel wall as the blood vessel contour line. However, the visibility of the fundus image is poor, and the extraction accuracy of the blood vessel contour line by the technique disclosed in Non-Patent Document 1 is not necessarily good. Therefore, there is a demand for a technique for correcting the contour lines of blood vessels extracted from captured images.

岡見雄貴,他3名,“白内障患者の眼底血管の口径算出”,第50回日本人間工学会中国・四国支部大会講演論文集,2017/12/9,p.49-50Yuki Okami, 3 others, “Calculation of the caliber of fundus blood vessels in cataract patients”, Proceedings of the 50th Annual Meeting of the Chugoku-Shikoku Branch of the Japan Ergonomics Society, 2017/12/9, p. 49-50

本発明は、撮像画像から抽出された血管の輪郭線を補正する技術を提供することを目的とする。 SUMMARY OF THE INVENTION An object of the present invention is to provide a technique for correcting the contour line of a blood vessel extracted from a captured image.

本発明に係る画像処理装置は、血管が撮像された撮像画像に基づいて生成された血管の輪郭線を補正する画像処理装置であって、輪郭線は、血管の幅方向における一方側の輪郭を示す第1輪郭線と、幅方向における他方の輪郭を示す第2輪郭線と、を含み、撮像画像と、撮像画像に撮像された血管の輪郭線に対応する輪郭画素の位置情報と、を関連付けて記憶する記憶部と、第1輪郭線と第2輪郭線とのいずれか一方を補正の対象輪郭線として選択する輪郭線選択部と、輪郭画素の位置情報に基づいて、対象輪郭線に対応する輪郭画素のうち、複数の輪郭画素を抽出画素として抽出する抽出部と、複数の抽出画素のうち、所定の位置関係を有する2つの抽出画素を選択画素として選択する画素選択部と、選択された2つの選択画素ごとに、2つの選択画素から所定位置に配置される複数の参照画素それぞれの画素情報に基づいて、撮像画像を構成する画素の中から補正後の対象輪郭線に対応する補正画素を決定する補正画素決定部と、補正画素の位置情報に基づいて、補正後の対象輪郭線を特定する補正輪郭線特定部と、を有してなる、ことを特徴とする。 An image processing apparatus according to the present invention is an image processing apparatus that corrects a contour line of a blood vessel generated based on a captured image of the blood vessel. and a second contour line indicating the other contour in the width direction. a storage unit for storing the following: a contour selection unit for selecting either one of the first contour and the second contour as a correction target contour; an extraction unit for extracting a plurality of contour pixels as extracted pixels from among the contour pixels, and a pixel selection unit for selecting two extracted pixels having a predetermined positional relationship from among the plurality of extracted pixels as selected pixels; For each of the two selected pixels, based on the pixel information of each of a plurality of reference pixels arranged at predetermined positions from the two selected pixels, correction corresponding to the target contour line after correction from among the pixels constituting the captured image It is characterized by comprising a correction pixel determination unit that determines pixels, and a correction contour line identification unit that identifies a target contour line after correction based on position information of the correction pixels.

本発明によれば、撮像画像から抽出された血管の輪郭線を補正することができる。 According to the present invention, it is possible to correct the contour line of a blood vessel extracted from a captured image.

本発明に係る画像処理装置の実施の形態を示すネットワーク接続図である。1 is a network connection diagram showing an embodiment of an image processing apparatus according to the present invention; FIG. 本発明に係る画像処理装置の実施の形態を示す機能ブロック図である。1 is a functional block diagram showing an embodiment of an image processing device according to the present invention; FIG. 図2の画像処理装置が備える記憶部に記憶される撮像画像の例を示す模式図である。FIG. 3 is a schematic diagram showing an example of a captured image stored in a storage unit included in the image processing apparatus of FIG. 2; 本発明に係る画像処理方法の実施の形態を示すフローチャートである。1 is a flow chart showing an embodiment of an image processing method according to the present invention; 図4の画像処理方法に含まれる1次補正処理のフローチャートである。5 is a flowchart of primary correction processing included in the image processing method of FIG. 4; 図5の1次補正処理の各段階の撮像画像を示す模式図である。6A and 6B are schematic diagrams showing captured images at each stage of the primary correction processing in FIG. 5 ; FIG. 図5の1次補正処理の途中段階の撮像画像を示す模式図である。6 is a schematic diagram showing a captured image in the middle of the primary correction process of FIG. 5; FIG. 図5の1次補正処理の参照画素の輝度値とその変化量との例を示す模式図である。FIG. 6 is a schematic diagram showing an example of the luminance value of a reference pixel and the amount of change thereof in the primary correction process of FIG. 5; 図4の画像処理方法に含まれる2次補正処理のフローチャートである。5 is a flowchart of secondary correction processing included in the image processing method of FIG. 4; 図8の2次補正処理の各段階の撮像画像を示す模式図である。9A and 9B are schematic diagrams showing captured images at each stage of the secondary correction process of FIG. 8; 図8の2次補正処理で生成される仮想円を示す模式図である。FIG. 9 is a schematic diagram showing a virtual circle generated by the secondary correction processing of FIG. 8; 本発明に係る画像処理装置の別の実施の形態を示す機能ブロック図である。FIG. 3 is a functional block diagram showing another embodiment of an image processing device according to the present invention; 本発明に係る画像処理方法の別の実施の形態を示すフローチャートである。Fig. 4 is a flow chart showing another embodiment of the image processing method according to the present invention; 図12の画像処理方法に含まれる1次補正処理のフローチャートである。13 is a flowchart of primary correction processing included in the image processing method of FIG. 12; 図13の1次補正処理の各段階の撮像画像を示す模式図である。14A and 14B are schematic diagrams showing captured images at each stage of the primary correction processing of FIG. 13; 図13の1次補正処理の途中段階の撮像画像を示す模式図である。FIG. 14 is a schematic diagram showing a captured image in the middle of the primary correction process of FIG. 13;

以下、図面を参照しながら、本発明に係る画像処理装置(以下「本装置」という。)と画像処理プログラム(以下「本プログラム」という。)と画像処理方法(以下「本方法」という。)との実施の形態について説明する。 Hereinafter, with reference to the drawings, an image processing apparatus (hereinafter referred to as "this apparatus"), an image processing program (hereinafter referred to as "this program"), and an image processing method (hereinafter referred to as "this method") according to the present invention. and an embodiment will be described.

本発明は、血管が撮像された撮像画像に基づいて生成された血管の輪郭線を補正するものである。具体的には、本発明は、撮像画像を構成する画素の画素情報に基づいて、撮像画像を構成する画素の中から補正後の輪郭線に対応する画素(以下「補正画素」という。)を決定し、補正画素の位置情報に基づいて補正後の輪郭線を特定する。また、本発明は、補正画素ごとの曲率に基づいて、補正画素の位置情報を補正する。 The present invention corrects the contour line of a blood vessel generated based on a captured image of the blood vessel. Specifically, the present invention selects pixels (hereinafter referred to as “corrected pixels”) corresponding to the contour line after correction from among the pixels forming the captured image based on the pixel information of the pixels forming the captured image. Then, the post-correction outline is specified based on the positional information of the corrected pixels. Also, the present invention corrects the position information of the correction pixel based on the curvature of each correction pixel.

「撮像画像」は、例えば、眼底カメラで撮像された眼底画像、X線撮像装置で撮像されたX線画像、CT装置で撮像されたCT画像、の他、通常の光学カメラで撮像された画像も含む。本発明において、撮像画像には、少なくとも血管が被写体として撮像される。前述のとおり、本実施の形態において、撮像画像は、眼底画像である。 The “captured image” is, for example, a fundus image captured by a fundus camera, an X-ray image captured by an X-ray imaging device, a CT image captured by a CT device, or an image captured by a normal optical camera. Also includes In the present invention, at least a blood vessel is captured as a subject in the captured image. As described above, in the present embodiment, the captured image is a fundus image.

「血管」は、人の血管だけでなく、人を除く動物の血管も含む。本実施の形態において、血管は、人の眼底の血管である。 "Vessels" includes not only human blood vessels, but also non-human animal blood vessels. In this embodiment, the blood vessel is a human fundus blood vessel.

「輪郭線」は、撮像画像に撮像されている血管の輪郭をなぞる線である。本実施の形態において、輪郭線は、例えば、撮像画像に公知の画像処理が実行されることにより生成される。ここで、血管は、幅を有する管体であるため、撮像画像に撮像された1本の血管の輪郭線は、血管の幅方向の一方側の輪郭を示す輪郭線(以下「第1輪郭線」という。)と、同幅方向の他方側の輪郭を示す輪郭線(以下「第2輪郭線」という。)と、を含む。 A “contour line” is a line that traces the contour of the blood vessel captured in the captured image. In the present embodiment, the contour line is generated, for example, by performing known image processing on the captured image. Here, since the blood vessel is a tubular body having a width, the contour line of one blood vessel captured in the captured image is a contour line indicating the contour on one side of the blood vessel in the width direction (hereinafter referred to as "first contour line ”) and a contour line indicating the contour on the other side in the same width direction (hereinafter referred to as “second contour line”).

「画素情報」は、画素ごとの情報である。本実施の形態において、画素情報は、画素の輝度値と、隣接する画素との間の輝度値の変化量(輝度勾配)と、を含む。 "Pixel information" is information for each pixel. In this embodiment, the pixel information includes the brightness value of the pixel and the amount of change in brightness value (brightness gradient) between adjacent pixels.

「位置情報」は、画素の位置に関する情報、すなわち、画素の座標を示す情報である。 "Positional information" is information about the position of a pixel, that is, information indicating the coordinates of a pixel.

●画像処理装置(1)●
先ず、本装置の実施の形態(第1実施形態)について説明する。
●Image processing device (1)●
First, an embodiment (first embodiment) of this device will be described.

図1は、本装置の実施の形態を示すネットワーク接続図である。
同図は、本装置1が、有線通信方式または無線通信方式を利用するネットワーク(通信回線)Nを介して、情報記憶装置2に接続されていることを示す。
FIG. 1 is a network connection diagram showing an embodiment of this device.
The figure shows that the device 1 is connected to the information storage device 2 via a network (communication line) N using a wired communication system or a wireless communication system.

本装置1は、撮像画像i1(図3参照。以下同じ。)を構成する画素の画素情報に基づいて、撮像画像i1に撮像されている血管の輪郭線を補正する。本装置1の具体的な構成と動作とは、後述する。 The apparatus 1 corrects the contour line of the blood vessel captured in the captured image i1 based on the pixel information of the pixels forming the captured image i1 (see FIG. 3; the same shall apply hereinafter). The specific configuration and operation of this device 1 will be described later.

情報記憶装置2は、撮像画像i1と、撮像画像i1に撮像されている血管の輪郭線に対応する画素(以下「輪郭画素」という。)の位置情報と、を関連付けて記憶する。情報記憶装置2は、例えば、情報記憶サーバ、NAS、パーソナルコンピュータである。 The information storage device 2 associates and stores the captured image i1 and the position information of the pixels (hereinafter referred to as "contour pixels") corresponding to the contour lines of the blood vessels captured in the captured image i1. The information storage device 2 is, for example, an information storage server, NAS, or personal computer.

ネットワークNは、例えば、インターネット、移動体通信網、LAN(Local Area Network)、WAN(Wide Area Network)、Wi-Fi(登録商標)、Bluetooth(登録商標)のような通信網である。 The network N is, for example, a communication network such as the Internet, a mobile communication network, a LAN (Local Area Network), a WAN (Wide Area Network), Wi-Fi (registered trademark), or Bluetooth (registered trademark).

●画像処理装置(1)の構成
図2は、本装置1の実施の形態を示す機能ブロック図である。
●Configuration of Image Processing Apparatus (1) FIG. 2 is a functional block diagram showing an embodiment of the apparatus 1. As shown in FIG.

本装置1は、例えば、パーソナルコンピュータで実現される。本装置1では、本プログラムが動作して、本プログラムが本装置1のハードウェア資源と協働して、本方法を実現する。 The device 1 is realized by, for example, a personal computer. In this device 1, this program operates, and this program cooperates with the hardware resources of this device 1 to realize this method.

ここで、コンピュータ(不図示)に本プログラムを実行させることで、本プログラムは、同コンピュータを本装置1と同様に機能させて、同コンピュータに本方法を実行させることができる。 Here, by causing a computer (not shown) to execute this program, the program can cause the computer to function in the same manner as the apparatus 1 and cause the computer to execute the method.

本装置1は、通信部11と、記憶部12と、制御部13と、を備える。 The device 1 includes a communication section 11 , a storage section 12 and a control section 13 .

通信部11は、ネットワークNを介して、情報記憶装置2から撮像画像i1を取得する。通信部11は、例えば、通信モジュールとアンテナとにより構成される。撮像画像i1は、記憶部12に記憶される。 The communication unit 11 acquires the captured image i1 from the information storage device 2 via the network N. FIG. The communication unit 11 is composed of, for example, a communication module and an antenna. The captured image i<b>1 is stored in the storage unit 12 .

記憶部12は、本装置1が後述する本方法を実行するために必要な情報を記憶する。記憶部12は、例えば、本装置1が備えるHDD(Hard Disk Drive)またはSSD(Solid State Drive)のような記録装置、RAM(Random Access Memory)のような揮発性メモリ、および/または、フラッシュメモリのような可搬性記憶媒体、により構成される。本実施の形態において、記憶部12は、予め情報記憶装置2から受信した撮像画像i1と、撮像画像i1に撮像された血管の輪郭画素の位置情報と、を関連付けて記憶している。 The storage unit 12 stores information necessary for the apparatus 1 to execute the method described later. The storage unit 12 is, for example, a recording device such as a HDD (Hard Disk Drive) or an SSD (Solid State Drive) included in the apparatus 1, a volatile memory such as a RAM (Random Access Memory), and/or a flash memory. A portable storage medium such as In the present embodiment, the storage unit 12 associates and stores the captured image i1 received in advance from the information storage device 2 and the position information of the contour pixels of the blood vessels captured in the captured image i1.

制御部13は、本装置1全体の動作を制御すると共に、後述する本方法を実行する。制御部13は、例えば、本装置1が備えるCPU(Central Processing Unit)と、CPUの作業領域として機能するRAMのような揮発性メモリと、本プログラムなどの各種情報を記憶するROM(Read Only Memory)のような不揮発性メモリと、により構成される。制御部13は、輪郭線選択部131と、抽出部132と、画素選択部133と、補正画素決定部134と、2次補正部135と、補正輪郭線特定部136と、を備える。 The control unit 13 controls the operation of the entire device 1 and executes the method described later. The control unit 13 includes, for example, a CPU (Central Processing Unit) provided in the device 1, a volatile memory such as a RAM that functions as a work area for the CPU, and a ROM (Read Only Memory) that stores various information such as the program. ) and a non-volatile memory such as The control unit 13 includes a contour line selection unit 131 , an extraction unit 132 , a pixel selection unit 133 , a correction pixel determination unit 134 , a secondary correction unit 135 , and a correction contour line identification unit 136 .

輪郭線選択部131は、輪郭線のうち、第1輪郭線と第2輪郭線とのいずれか一方を補正の対象の輪郭線(以下「対象輪郭線」という。)として選択する。輪郭線選択部131の具体的な動作は、後述する。 The contour selection unit 131 selects either one of the first contour and the second contour as a contour to be corrected (hereinafter referred to as a “target contour”). A specific operation of the outline selection unit 131 will be described later.

抽出部132は、輪郭画素の位置情報に基づいて、対象輪郭線を構成する(対象輪郭線に対応する)輪郭画素のうち、複数の輪郭画素を抽出画素として抽出する。抽出部132の具体的な動作は、後述する。 The extraction unit 132 extracts a plurality of contour pixels as extracted pixels from contour pixels forming the target contour (corresponding to the target contour) based on the position information of the contour pixels. A specific operation of the extraction unit 132 will be described later.

画素選択部133は、複数の抽出画素のうち、所定の位置関係を有する2つの抽出画素を選択画素として選択する。画素選択部133の具体的な動作と所定の位置関係とは、後述する。 The pixel selection unit 133 selects two extracted pixels having a predetermined positional relationship among the plurality of extracted pixels as selected pixels. A specific operation of the pixel selection unit 133 and a predetermined positional relationship will be described later.

補正画素決定部134は、選択された2つの選択画素ごとに、2つの選択画素から同じ距離の位置に配置される複数の画素(以下「参照画素」という。)それぞれの画素情報に基づいて、撮像画像i1を構成する画素の中から補正後の対象輪郭線に対応する補正画素を決定する。補正画素決定部134の具体的な動作は、後述する。 Based on the pixel information of each of the pixels arranged at the same distance from the two selected pixels (hereinafter referred to as "reference pixels"), the correction pixel determination unit 134 determines: A corrected pixel corresponding to the target contour line after correction is determined from the pixels forming the captured image i1. A specific operation of the correction pixel determination unit 134 will be described later.

「同じ距離の位置」は、本発明における所定位置の例であり、撮像画像i1において、2つの選択画素からの直線距離が同じとなる位置であることを意味する。ここで、撮像画像i1は複数の画素から構成されるため、2つの選択画素からの直線距離が完全に同じとなる位置に単一の画素が存在しないことがあり得る。そのため、同じ距離の位置に配置される画素は、2つの選択画素からの直線距離が同じとなる位置に配置される画素だけでなく、同位置に隣接する画素(2つの選択画素からの直線距離が「略同じ」となる画素)を含み得る。すなわち、例えば、2つの選択画素の間に偶数(例えば、10)の画素が有る場合、1の選択画素から6画素分の位置にある画素(他の1の選択画素から4画素分の位置にある画素)は、2つの選択画素からの直線距離が同じとなる画素となり得る。 "A position at the same distance" is an example of a predetermined position in the present invention, and means a position at which the linear distance from two selected pixels is the same in the captured image i1. Here, since the captured image i1 is composed of a plurality of pixels, it is possible that a single pixel does not exist at a position where the linear distances from the two selected pixels are exactly the same. Therefore, pixels arranged at the same distance from two selected pixels are not only pixels arranged at positions with the same linear distance from two selected pixels, but also pixels adjacent to the same position (linear distance from two selected pixels). are “substantially the same”). That is, for example, when there are an even number (for example, 10) of pixels between two selected pixels, a pixel located at a position of 6 pixels from one selected pixel (a pixel located at a position of 4 pixels from another selected pixel). A pixel) can be the pixel that has the same linear distance from the two selected pixels.

2次補正部135は、補正画素の位置情報を補正する。2次補正部135の具体的な動作は、後述する。 A secondary correction unit 135 corrects the position information of the correction pixel. A specific operation of the secondary correction unit 135 will be described later.

補正輪郭線特定部136は、補正画素の位置情報に基づいて、補正後の対象輪郭線を特定する。補正輪郭線特定部136の具体的な動作は、後述する。 The corrected contour specifying unit 136 specifies the target contour after correction based on the position information of the corrected pixels. A specific operation of the correction contour specifying unit 136 will be described later.

●画像処理装置(1)の動作
次に、本装置1の動作、すなわち、本装置1が実行する本方法の実施の形態について、図2も参照しつつ説明する。
●Operation of Image Processing Apparatus (1) Next, the operation of the apparatus 1, that is, the embodiment of the method executed by the apparatus 1 will be described with reference to FIG.

以下の説明において、本装置1はネットワークNを介して情報記憶装置2から撮像画像i1と、撮像画像i1に撮像された血管の輪郭線に対応する輪郭画素の位置情報と、を取得し、両者を関連付けて記憶部12に記憶しているものとする。また、以下の説明で用いられる撮像画像i1には、輪郭線が重畳されているものとする。さらに、以下の説明において、輪郭線には符号「L」が付され、第1輪郭線には符号「L1」が付され、第2輪郭線には符号「L2」が付され、対象輪郭線には符号「Lt」が付されるものとする。さらに、以下の説明において、各画素を区別するため、各画素には、「P」を頭文字とする符号が付されるものとする。 In the following description, the device 1 acquires a captured image i1 from the information storage device 2 via the network N and position information of contour pixels corresponding to the contour lines of the blood vessels captured in the captured image i1. are stored in the storage unit 12 in association with each other. Also, it is assumed that a contour line is superimposed on the captured image i1 used in the following description. Further, in the following description, the contour is denoted by "L", the first contour is denoted by "L1", the second contour is denoted by "L2", and the target contour shall be given the code "Lt". Furthermore, in the following description, each pixel is assigned a reference numeral beginning with "P" in order to distinguish each pixel.

図3は、撮像画像i1の例を示す模式図である。
同図は、説明の便宜上、1つの血管の一部を拡大して示す。同図は、撮像画像i1において、撮像されている血管の一部に輪郭線L(第1輪郭線L1、第2輪郭線L2)が重畳されていることを示す。
FIG. 3 is a schematic diagram showing an example of the captured image i1.
For convenience of explanation, this figure shows an enlarged part of one blood vessel. The figure shows that contour lines L (first contour line L1, second contour line L2) are superimposed on a part of the imaged blood vessel in the captured image i1.

図4は、本方法の実施の形態を示すフローチャートである。 FIG. 4 is a flow chart illustrating an embodiment of the method.

本方法は、1次補正処理(S1)と2次補正処理(S2)とを実行する。 The method performs a primary correction process (S1) and a secondary correction process (S2).

●1次補正処理
先ず、本装置1は、1次補正処理(S1)を実行する。
● Primary Correction Processing First, the apparatus 1 executes the primary correction processing (S1).

図5は、1次補正処理(S1)のフローチャートである。
図6は、1次補正処理(S1)の各段階の撮像画像i1を示す模式図である。
FIG. 5 is a flow chart of the primary correction process (S1).
FIG. 6 is a schematic diagram showing a captured image i1 at each stage of the primary correction process (S1).

「1次補正処理(S1)」は、撮像画像i1を構成する画素の画素情報に基づいて、撮像画像i1に関連付けられている輪郭線Lに対応する輪郭画素の位置情報を補正する処理である。 The "primary correction process (S1)" is a process of correcting the position information of the contour pixels corresponding to the contour line L associated with the captured image i1 based on the pixel information of the pixels forming the captured image i1. .

先ず、輪郭線選択部131は、第1輪郭線L1と第2輪郭線L2とのいずれか一方を対象輪郭線Ltとして選択する(S101)。すなわち、例えば、輪郭線選択部131は、撮像画像i1における、第1輪郭線L1および第2輪郭線L2の位置関係(例えば、撮像画像i1を構成する画素の座標の原点(0,0)との距離)に基づいて、対象輪郭線Ltを選択する。 First, the contour line selection unit 131 selects one of the first contour line L1 and the second contour line L2 as the target contour line Lt (S101). That is, for example, the contour line selection unit 131 determines the positional relationship between the first contour line L1 and the second contour line L2 in the captured image i1 (for example, the origin (0, 0) of the coordinates of the pixels forming the captured image i1). distance), the target contour line Lt is selected.

なお、本発明における輪郭線選択部は、本装置の使用者が選択手段(不図示)を用いて選択した輪郭線を対象輪郭線として選択してもよく、あるいは、血管に対する相対的な位置関係(血管に対する上下左右方向など)に基づいて対象輪郭線を選択してもよい。 Note that the contour line selection unit in the present invention may select a contour line selected by the user of the apparatus using a selection means (not shown) as the target contour line, or a relative positional relationship to the blood vessel. The target contour line may be selected based on (eg, up, down, left, or right with respect to the blood vessel).

次いで、補正画素決定部134は、対象輪郭線Ltとして選択されていない輪郭線L(本実施の形態では、第2輪郭線L2)を含む所定領域に含まれる画素にマスク処理を実行して(S102)、同所定領域を後述する補正画素の対象から除外する。ここで、補正画素決定部134は、例えば、第2輪郭線L2に公知の膨張処理を実行して、膨張化された第2輪郭線L2に覆われる領域(所定領域)を画定する。 Next, the correction pixel determining unit 134 performs mask processing on pixels included in a predetermined region including the contour line L (the second contour line L2 in the present embodiment) that is not selected as the target contour line Lt ( S102), the same predetermined area is excluded from correction pixel targets, which will be described later. Here, the correction pixel determination unit 134 performs, for example, a known expansion process on the second contour line L2, and defines an area (predetermined area) covered by the expanded second contour line L2.

なお、マスク処理において補正画素の対象から除外される所定領域は、対象輪郭線として選択されていない輪郭線を基準とした一定の領域でもよい。また、マスク処理は、後述する参照画素の輝度値の変化量の算出までに実行されればよい。 Note that the predetermined area excluded from the correction pixel target in the mask processing may be a certain area based on the contour line not selected as the target contour line. Also, the masking process may be executed before calculating the amount of change in luminance value of a reference pixel, which will be described later.

次いで、抽出部132は、撮像画像i1を構成する画素のうち、対象輪郭線Ltに対応する輪郭画素の位置情報に基づいて、対象輪郭線Ltに対応する輪郭画素のうち、複数の輪郭画素を抽出画素P1(n)として抽出する(S103)。すなわち、例えば、抽出部132は、複数の輪郭画素から、略等間隔に一定数の輪郭画素を抽出画素P1(n)として抽出する。以下の説明において、抽出部132は、6つの抽出画素P1(1)~P1(6)を抽出したものとする。 Next, the extraction unit 132 extracts a plurality of contour pixels out of the contour pixels corresponding to the target contour line Lt based on the position information of the contour pixels corresponding to the target contour line Lt among the pixels forming the captured image i1. It is extracted as an extraction pixel P1(n) (S103). That is, for example, the extracting unit 132 extracts a certain number of contour pixels at approximately equal intervals as extracted pixels P1(n) from a plurality of contour pixels. In the following description, it is assumed that the extraction unit 132 has extracted six extraction pixels P1(1) to P1(6).

なお、本発明における抽出部は、対象輪郭線に対応する輪郭画素の全てを抽出画素として抽出してもよい。 Note that the extraction unit in the present invention may extract all contour pixels corresponding to the target contour line as extracted pixels.

次いで、画素選択部133は、複数の抽出画素P1(1)~P1(6)のうち、隣り合う2つの抽出画素P1(n),P1(n+1)を選択画素P2(1),P2(2)として選択する(S104)。すなわち、例えば、6つの抽出画素P1(1)~P1(6)が番号順に配列されているとき、画素選択部133は、抽出画素P1(1),P1(2)を選択画素P2(1),P2(2)として選択する。 Next, the pixel selection unit 133 selects adjacent two extracted pixels P1(n) and P1(n+1) among the plurality of extracted pixels P1(1) to P1(6) as selected pixels P2(1) and P2(2). ) (S104). That is, for example, when six extracted pixels P1(1) to P1(6) are arranged in numerical order, the pixel selection unit 133 selects the extracted pixels P1(1) and P1(2) as the selected pixel P2(1). , P2(2).

ここで、本実施の形態において、「隣り合う」位置関係は、本発明における所定の位置関係の例である。 Here, in the present embodiment, the "adjacent" positional relationship is an example of a predetermined positional relationship in the present invention.

図7は、1次補正処理(S1)の途中段階の撮像画像i1を示す模式図である。
同図は、説明の便宜上、2つの選択画素P2(1),P2(2)の近傍を拡大して示す。
FIG. 7 is a schematic diagram showing a captured image i1 in the middle of the primary correction process (S1).
For convenience of explanation, the drawing shows an enlarged view of the vicinity of two selected pixels P2(1) and P2(2).

図5~図7に戻る。
次いで、補正画素決定部134は、2つの選択画素P2(1),P2(2)同士を結ぶ仮想線分Lmを生成し、仮想線分Lmの中点において、仮想線分Lmに直交する仮想垂線Lnを生成する(S105)。ここで、仮想垂線Lnの長さは、例えば、所定の長さ(例えば、中点を中心画素とする9画素)に予め設定されている。
Returning to FIGS. 5-7.
Next, the correction pixel determining unit 134 generates a virtual line segment Lm connecting the two selected pixels P2(1) and P2(2), and creates a virtual line segment Lm that is orthogonal to the virtual line segment Lm at the midpoint of the virtual line segment Lm. A perpendicular line Ln is generated (S105). Here, the length of the virtual perpendicular line Ln is set in advance to, for example, a predetermined length (for example, 9 pixels with the midpoint as the central pixel).

なお、仮想垂線の長さは、例えば、第1輪郭線と第2輪郭線との間の距離(すなわち、血管の幅、口径)に基づいて、設定されてもよい。すなわち、例えば、仮想垂線の長さは、同距離の2倍未満に設定されてもよい。 Note that the length of the virtual perpendicular line may be set, for example, based on the distance between the first contour line and the second contour line (that is, the width and diameter of the blood vessel). That is, for example, the length of the virtual perpendicular may be set to less than twice the same distance.

次いで、補正画素決定部134は、撮像画像i1を構成する画素のうち、仮想垂線Lnに対応する複数の画素を参照画素P3(n)として特定する(S106)。すなわち、例えば、仮想垂線Lnの長さが9画素の場合、9つの画素を参照画素P3(1)~P3(9)として特定する。ここで、仮想垂線Lnに対応する画素は、2つの選択画素P2(1),P2(2)から略同じ距離の位置に配置される画素である。すなわち、仮想垂線Lnに対応する画素は、本発明における所定位置に配置される画素である。 Next, the correction pixel determination unit 134 identifies, as reference pixels P3(n), a plurality of pixels corresponding to the virtual perpendicular line Ln among the pixels forming the captured image i1 (S106). That is, for example, when the length of the virtual perpendicular line Ln is 9 pixels, the 9 pixels are specified as the reference pixels P3(1) to P3(9). Here, the pixels corresponding to the imaginary perpendicular line Ln are pixels arranged at substantially the same distance from the two selected pixels P2(1) and P2(2). That is, the pixels corresponding to the imaginary perpendicular line Ln are pixels arranged at predetermined positions in the present invention.

次いで、補正画素決定部134は、記憶部12から参照画素P3(1)~P3(9)それぞれの輝度値を取得して、複数の参照画素P3(1)~P3(9)それぞれの輝度値に基づいて、隣接する2つの参照画素P3(n),P3(n+1)ごとの輝度値の変化量を取得する(S107)。すなわち、例えば、補正画素決定部134は、参照画素P3(1)~P3(9)それぞれの1次輝度勾配と2次輝度勾配とを算出する。 Next, the correction pixel determining unit 134 acquires the luminance values of the reference pixels P3(1) to P3(9) from the storage unit 12, and calculates the luminance values of the plurality of reference pixels P3(1) to P3(9). Based on (S107), the amount of change in luminance value for each of two adjacent reference pixels P3(n) and P3(n+1) is obtained. That is, for example, the correction pixel determination unit 134 calculates the primary luminance gradient and the secondary luminance gradient of each of the reference pixels P3(1) to P3(9).

なお、本発明における補正画素決定部は、1次輝度勾配のみを算出してもよく、あるいは、他の計算により輝度値の変化量を算出してもよい。 Note that the correction pixel determination unit in the present invention may calculate only the primary luminance gradient, or may calculate the amount of change in the luminance value by other calculations.

図8は、各参照画素P3(n)の画素値と、その変化量(1次勾配)の例を示す模式図である。同図は、参照画素P3(6)が輝度値の変化量が最大であることを示す。 FIG. 8 is a schematic diagram showing an example of the pixel value of each reference pixel P3(n) and its variation (primary gradient). The figure shows that the reference pixel P3(6) has the largest amount of change in luminance value.

図5~図8に戻る。
次いで、補正画素決定部134は、複数の参照画素P3(1)~P3(9)のうち、輝度値の変化量が最大の参照画素P3(n)を補正画素P4(n)として決定する(S108)。通常、撮像画像i1において、被写体である血管は、背景よりも暗く撮像される。そのため、輝度値の変化量の大きい画素は、血管のエッジに相当する画素とみなすことができる。
Returning to FIGS. 5-8.
Next, the correction pixel determination unit 134 determines the reference pixel P3(n) having the largest luminance value variation among the plurality of reference pixels P3(1) to P3(9) as the correction pixel P4(n) ( S108). Usually, in the captured image i1, the blood vessel, which is the subject, is imaged darker than the background. Therefore, a pixel with a large amount of change in luminance value can be regarded as a pixel corresponding to the edge of a blood vessel.

ここで、仮想垂線Lnは、対象輪郭線Ltとして選択されていない輪郭線L(第2輪郭線L2)上を通過し得る。そのため、画素値の変化量が最大の参照画素P3(n)は、第2輪郭線L2近辺にも存在し得る。その結果、補正画素決定部134は、第2輪郭線L2近辺の参照画素P3(n)を補正画素P4(n)として決定し得る。しかしながら、本実施の形態では、処理(S102)のマスク処理において、第2輪郭線L2近辺の参照画素P3(n)は、補正画素P4(n)の対象から除外されている。そのため、処理(S108)において、輝度値の変化量が最大の画素は、単一の参照画素P3(n)となる。 Here, the virtual perpendicular line Ln can pass over the contour line L (the second contour line L2) that is not selected as the target contour line Lt. Therefore, the reference pixel P3(n) with the largest change in pixel value can also exist near the second contour line L2. As a result, the correction pixel determination unit 134 can determine the reference pixel P3(n) near the second contour line L2 as the correction pixel P4(n). However, in the present embodiment, the reference pixels P3(n) near the second contour line L2 are excluded from the correction pixels P4(n) in the masking process (S102). Therefore, in the process (S108), the single reference pixel P3(n) is the pixel with the largest change in luminance value.

なお、本発明における補正画素決定部は、前述のマスク処理を実行しなくてもよい。この場合、本発明における補正画素決定部は、例えば、画素値の変化量が大きい2つの参照画素のうち、対象輪郭線に近い画素を補正画素として決定する。 It should be noted that the correction pixel determination unit in the present invention does not need to perform the mask processing described above. In this case, the correction pixel determining unit of the present invention, for example, determines a pixel closer to the target contour line as a correction pixel, out of two reference pixels having a large amount of change in pixel value.

次いで、制御部13(たとえば、補正画素決定部134)は、選択画素P2(1),P2(2)として未選択の抽出画素P1(n)の有無を判定する(S109)。選択画素P2(1),P2(2)として未選択の抽出画素P1(n)が有るとき(S109の「Y」)、1次補正処理(S1)は、処理(S103)に戻る。すなわち、例えば、画素選択部133は、未選択の抽出画素P1(3)~P1(6)のうち、既に選択済みの抽出画素P1(2)に隣り合う抽出画素P1(3)と、既に選択済みの抽出画素P1(1),P1(2)のうち、未選択の抽出画素P1(3)に隣り合う抽出画素P1(2)と、を次の選択画素P2(1),P2(2)として選択する。 Next, the control unit 13 (for example, the correction pixel determination unit 134) determines whether there is an unselected extracted pixel P1(n) as the selected pixels P2(1) and P2(2) (S109). When there is an unselected extracted pixel P1(n) as the selected pixels P2(1) and P2(2) ("Y" in S109), the primary correction process (S1) returns to the process (S103). That is, for example, the pixel selection unit 133 selects the extracted pixel P1(3) adjacent to the already selected extracted pixel P1(2) among the unselected extracted pixels P1(3) to P1(6), and the already selected extracted pixel P1(3). Among the already extracted pixels P1(1) and P1(2), the extracted pixel P1(2) adjacent to the unselected extracted pixel P1(3) is used as the next selected pixels P2(1) and P2(2). Select as

一方、選択画素P2(1),P2(2)として未選択の抽出画素P1(n)が無いとき(S109の「N」)、補正画素決定部134は、全ての補正画素P4(1)~P4(5)の位置情報を、補正画素群の位置情報として撮像画像i1に関連付けて記憶部12に記憶する(S110)。 On the other hand, when there is no unselected extracted pixel P1(n) as the selected pixels P2(1) and P2(2) (“N” in S109), the correction pixel determination unit 134 determines all the correction pixels P4(1) to The position information of P4(5) is stored in the storage unit 12 as the position information of the corrected pixel group in association with the captured image i1 (S110).

次いで、制御部13(例えば、補正画素決定部134)は、対象輪郭線Ltとして未選択の輪郭線Lの有無を判定する(S111)。対象輪郭線Ltとして未選択の輪郭線Lが有るとき(S111の「Y」)、1次補正処理(S1)は、処理(S101)に戻る。一方、対象輪郭線Ltとして未選択の輪郭線Lが無いとき(S111の「N」)、1次補正処理(S1)は、終了する。 Next, the control unit 13 (for example, the correction pixel determination unit 134) determines whether or not there is an unselected contour line L as the target contour line Lt (S111). When there is a contour line L that has not been selected as the target contour line Lt ("Y" in S111), the primary correction process (S1) returns to the process (S101). On the other hand, when there is no contour line L that has not been selected as the target contour line Lt ("N" in S111), the primary correction process (S1) ends.

●2次補正処理
次いで、本装置1は、2次補正処理(S2)を実行する。
●Secondary Correction Processing Next, the apparatus 1 executes the secondary correction processing (S2).

図9は、2次補正処理(S2)のフローチャートである。
図10は、2次補正処理(S2)の各段階の撮像画像i1を示す模式図である。
図10は、静脈が動脈と交叉して、静脈の一部が動脈の一部に覆い隠されている状態を示す。
FIG. 9 is a flowchart of secondary correction processing (S2).
FIG. 10 is a schematic diagram showing captured images i1 at each stage of the secondary correction process (S2).
FIG. 10 shows a state where a vein crosses an artery and a portion of the vein is obscured by a portion of the artery.

「2次補正処理(S2)」は、1次補正処理(S1)で決定された補正画素P4(n)の位置情報を補正する処理である。前述のとおり、撮像画像i1において血管は背景よりも暗く撮像されるため、輝度値の変化量で補正する1次補正処理(S1)により決定される補正画素P4(n)は、略正確に血管の輪郭に対応する。しかしながら、例えば、動脈と静脈とが交叉する血管の交叉部では、動脈の輝度値と静脈との輝度値との間の差異が少ない。そのため、例えば、交叉部における静脈の補正画素P4(n)は、動脈側に突出しまたは静脈内に凹み得る。2次補正処理(S2)は、このように、輝度値の差異が少ない部位において決定された補正画素P4(n)の位置情報を、輝度値とは異なる要素(曲率)に基づいて補正する処理である。 The "secondary correction process (S2)" is a process of correcting the position information of the correction pixel P4(n) determined in the first correction process (S1). As described above, the blood vessels are imaged darker than the background in the captured image i1. corresponds to the contour of However, for example, at a blood vessel intersection where an artery and a vein intersect, there is little difference between the brightness value of the artery and the brightness value of the vein. Thus, for example, the venous correction pixel P4(n) at the intersection may protrude toward the artery or recede into the vein. The secondary correction process (S2) is a process of correcting the positional information of the correction pixel P4(n) determined in the region where the luminance value difference is small based on the element (curvature) different from the luminance value. is.

先ず、2次補正部135は、2次補正の対象となる補正画素群を記憶部12から読み出す(S201)。以下の説明において、補正画素群は、5つの補正画素P4(11)~P4(15)で構成されているものとする。 First, the secondary correction unit 135 reads a correction pixel group to be subjected to secondary correction from the storage unit 12 (S201). In the following description, it is assumed that the correction pixel group is composed of five correction pixels P4(11) to P4(15).

次いで、2次補正部135は、補正画素群を構成する補正画素P4(11)~P4(15)のうち、始点となる補正画素P4(11)と終点となる補正画素P4(15)とを除く3つの補正画素P4(12)~P4(14)それぞれに対応する曲率を算出する。ここで、始点と終点それぞれは、例えば、補正画素P4(n)の決定順などにより、予め決定されている。 Next, the secondary correction unit 135 selects the correction pixel P4(11) as the start point and the correction pixel P4(15) as the end point among the correction pixels P4(11) to P4(15) forming the correction pixel group. Curvatures corresponding to each of the three correction pixels P4(12) to P4(14) are calculated. Here, each of the start point and the end point is determined in advance by, for example, the order of determination of the correction pixel P4(n).

具体的には、2次補正部135は、補正画素P4(12)~P4(14)のうち、始点に最も近く、かつ、曲率が未算出の補正画素P4(12)を算出画素P5(1)として選択する(S202)。次いで、2次補正部135は、算出画素P5(1)の両隣に配置される2つの補正画素P4(11),P4(13)を算出補助画素P5(2),P5(3)として選択する(S203)。 Specifically, of the corrected pixels P4(12) to P4(14), the secondary correction unit 135 converts the corrected pixel P4(12), which is closest to the starting point and whose curvature has not been calculated, to the calculated pixel P5(1). ) (S202). Next, the secondary correction unit 135 selects two correction pixels P4(11) and P4(13) arranged on both sides of the calculation pixel P5(1) as calculation auxiliary pixels P5(2) and P5(3). (S203).

次いで、2次補正部135は、算出画素P5(1)と2つの算出補助画素P5(2),P5(3)とを結ぶ仮想円Cを生成する(S204)。 Next, the secondary correction unit 135 generates a virtual circle C that connects the calculated pixel P5(1) and the two calculated auxiliary pixels P5(2) and P5(3) (S204).

図11は、仮想円Cを示す模式図である。
同図は、説明の便宜上、仮想円Cの半径「R」の算出に用いるパラメータも示す。
FIG. 11 is a schematic diagram showing the virtual circle C. As shown in FIG.
For convenience of explanation, the figure also shows parameters used to calculate the radius "R" of the virtual circle C. As shown in FIG.

図9~図11に戻る。
次いで、2次補正部135は、次式(1)を用いて、仮想円Cの半径「R」を算出する(S205)。
Returning to FIGS. 9-11.
Next, the secondary correction unit 135 calculates the radius "R" of the virtual circle C using the following formula (1) (S205).

Figure 2022128882000002
ここで、
「||・||」はベクトルのノルムを意味する(以下同じ)。
「pi-1」は、算出補助画素P5(2)の座標である。
「pi+1」は、算出補助画素P5(3)の座標である。
「q」は、算出補助画素P5(2)と仮想円Cの中心点「O」とを結ぶ仮想線Loと、仮想円Cとの交点「Q」の座標である。
「θ」は、算出補助画素P5(3)と交点「Q」とを結ぶ仮想線Lpと、仮想線Loと、の間の角度である。
Figure 2022128882000002
here,
"||·||" means the norm of the vector (same below).
“p i−1 ” is the coordinates of the calculated auxiliary pixel P5(2).
“p i+1 ” is the coordinates of the calculated auxiliary pixel P5(3).
"q" is the coordinate of the intersection point "Q" between the virtual line Lo connecting the calculated auxiliary pixel P5(2) and the center point "O" of the virtual circle C and the virtual circle C;
“θ” is the angle between the virtual line Lp connecting the calculation auxiliary pixel P5(3) and the intersection point “Q” and the virtual line Lo.

次いで、2次補正部135は、次式(2)を用いて、仮想円Cの曲率「k」を算出する(S206)。 Next, the secondary correction unit 135 calculates the curvature "k i " of the virtual circle C using the following equation (2) (S206).

Figure 2022128882000003
Figure 2022128882000003

算出された曲率「k」は、対応する算出画素P5(1)(補正画素P4(12))に関連付けられて記憶部12に記憶される。 The calculated curvature “k i ” is stored in the storage unit 12 in association with the corresponding calculated pixel P5(1) (corrected pixel P4(12)).

次いで、2次補正部135は、曲率「k」が未算出の補正画素P4(n)の有無を判定する(S207)。曲率「k」が未算出の補正画素P4(n)が有るとき(S207の「Y」)、2次補正処理(S2)は、処理(S202)に戻る。 Next, the secondary correction unit 135 determines whether or not there is a corrected pixel P4(n) for which the curvature "k i " has not been calculated (S207). When there is a correction pixel P4(n) whose curvature "k i " has not been calculated ("Y" in S207), the secondary correction process (S2) returns to the process (S202).

一方、曲率「k」が未算出の補正画素P4(n)が無いとき(S207の「N」)、2次補正部135は、全ての曲率「k」の平均値「Vave.」と標準偏差「σ」とを算出する(S208)。 On the other hand, when there is no correction pixel P4(n) for which the curvature "k i " has not yet been calculated ("N" in S207), the secondary correction unit 135 calculates the average value "V ave. " of all the curvatures "k i ". and the standard deviation "σ" are calculated (S208).

次いで、2次補正部135は、算出された平均値「Vave.」と標準偏差「σ」とに基づいて、曲率「k」が相対的に大きい補正画素P4(n)の有無を判定する(S209)。すなわち、例えば、2次補正部135は、平均値「Vave.」を中心とする正規分布の「±2σ」から外れる曲率「k」を有する補正画素P4(n)の有無を判定する。すなわち、本実施の形態において、「曲率「k」が相対的に大きい」は、曲率「k」が平均値「Vave.」を中心とする正規分布の「±2σ」から外れることを意味する。 Next, the secondary correction unit 135 determines whether there is a corrected pixel P4(n) having a relatively large curvature “k i ” based on the calculated average value “V ave. ” and standard deviation “σ”. (S209). That is, for example, the secondary correction unit 135 determines whether or not there is a correction pixel P4(n) having a curvature "k i " that deviates from "±2σ" of the normal distribution centered on the average value "V ave. ". That is, in the present embodiment, "the curvature 'k i ' is relatively large" means that the curvature 'k i ' deviates from '±2σ' of the normal distribution centered on the average value 'V ave. means.

曲率「k」が相対的に大きい補正画素P4(n)が有るとき(S209の「Y」)、2次補正部135は、同補正画素P4(n)のうち、始点に最も近く、かつ、位置補正画素P6(n)として選択されていない補正画素P4(n)を位置補正画素P6(1)として選択する(S210)。以下の説明において、補正画素P4(13)が位置補正画素P6(1)として選択されたものとする。 When there is a corrected pixel P4(n) with a relatively large curvature “k i ” (“Y” in S209), the secondary correction unit 135 selects the corrected pixel P4(n) closest to the starting point and , the correction pixel P4(n) that has not been selected as the position correction pixel P6(n) is selected as the position correction pixel P6(1) (S210). In the following description, it is assumed that the correction pixel P4(13) is selected as the position correction pixel P6(1).

次いで、2次補正部135は、位置補正画素P6(1)の両隣に配置される2つの補正画素P4(12),P4(14)を基準画素P7(1),P7(2)として選択する(S211)。 Next, the secondary correction unit 135 selects two correction pixels P4(12) and P4(14) arranged on both sides of the position correction pixel P6(1) as reference pixels P7(1) and P7(2). (S211).

次いで、2次補正部135は、基準画素P7(1),P7(2)の中点に対応する画素(以下「中点画素」という。)Poの位置情報を算出する(S212)。 Next, the secondary correction unit 135 calculates the position information of the pixel (hereinafter referred to as "midpoint pixel") Po corresponding to the midpoint of the reference pixels P7(1) and P7(2) (S212).

次いで、2次補正部135は、位置補正画素P6(1)の位置情報を中点画素Poの位置情報に補正する(S213)。すなわち、2次補正部135は、位置補正画素P6(1)に対応する補正画素P4(13)の位置を中点画素Poの位置に移動させる。つまり、2次補正部135は、基準画素P7(1),P7(2)の位置情報に基づいて、位置補正画素P6(1)の位置情報を補正する。 Next, the secondary correction unit 135 corrects the position information of the position correction pixel P6(1) to the position information of the midpoint pixel Po (S213). That is, the secondary correction unit 135 moves the position of the correction pixel P4(13) corresponding to the position correction pixel P6(1) to the position of the midpoint pixel Po. That is, the secondary correction unit 135 corrects the position information of the position correction pixel P6(1) based on the position information of the reference pixels P7(1) and P7(2).

次いで、2次補正部135は、位置補正画素P6(n)として未選択の補正画素P4(n)の有無を判定する(S214)。位置補正画素P6(n)として未選択の補正画素P4(n)が有るとき(S214の「Y」)、2次補正処理(S2)は、処理(S210)に戻る。 Next, the secondary correction unit 135 determines whether or not there is an unselected correction pixel P4(n) as the position correction pixel P6(n) (S214). When there is an unselected correction pixel P4(n) as the position correction pixel P6(n) ("Y" in S214), the secondary correction process (S2) returns to the process (S210).

一方、曲率「k」が相対的に大きい補正画素P4(n)が無いとき(S209の「N」)、または、位置補正画素P6(n)として未選択の補正画素P4(n)が無いとき(S214の「N」)、2次補正部135は、位置情報が補正されていない補正画素P4(n)の位置情報と、位置情報が補正された補正画素P4(n)の位置情報と、を補正後の輪郭線Lxに対応する画素の位置情報として記憶部12に記憶する(S215)。 On the other hand, when there is no correction pixel P4(n) with a relatively large curvature "k i "("N" in S209), or there is no correction pixel P4(n) unselected as the position correction pixel P6(n) When (“N” in S214), the secondary correction unit 135 combines the position information of the corrected pixel P4(n) whose position information is not corrected and the position information of the corrected pixel P4(n) whose position information is corrected. , are stored in the storage unit 12 as the position information of the pixels corresponding to the contour line Lx after correction (S215).

次いで、補正輪郭線特定部136は、補正後の輪郭線Lxに対応する画素の位置情報に基づいて、補正後の輪郭線Lxを特定する(S216)。 Next, the corrected contour line identifying unit 136 identifies the corrected contour line Lx based on the position information of the pixels corresponding to the corrected contour line Lx (S216).

このように、本装置1は、1次補正処理(S1)において、輪郭線Lから抽出した抽出画素P1(n)の画素情報に基づいて、補正画素P4(n)を決定して輪郭線Lの位置情報を補正すると共に、2次補正処理(S2)において、さらに補正画素P4(n)の曲率「k」に基づいて、補正画素P4(n)の位置情報を補正する。すなわち、本装置1は、画素情報と曲率「k」という異なる2つの情報に基づいて、輪郭線Lの位置情報を2段階に補正する。その結果、本装置1は、撮像画像i1から抽出された血管の輪郭線Lを高精度に補正することができる。 In this way, in the primary correction process (S1), the apparatus 1 determines the correction pixel P4(n) based on the pixel information of the extraction pixel P1(n) extracted from the contour line L, and corrects the contour line L. is corrected, and in the secondary correction process (S2), the position information of the corrected pixel P4(n) is further corrected based on the curvature "k i " of the corrected pixel P4(n). That is, the apparatus 1 corrects the position information of the contour line L in two stages based on two different pieces of information, ie, the pixel information and the curvature "k i ". As a result, the apparatus 1 can highly accurately correct the contour line L of the blood vessel extracted from the captured image i1.

●まとめ(1)
以上説明した実施の形態によれば、本装置1は、第1輪郭線L1と第2輪郭線L2とのいずれか一方を補正の対象輪郭線Ltとして選択する輪郭線選択部131と、輪郭画素の位置情報に基づいて対象輪郭線Ltに対応する輪郭画素のうち、複数の輪郭画素を抽出画素P1(n)として抽出する抽出部132と、複数の抽出画素P1(n)のうち、所定の位置関係を有する2つの抽出画素P1(n)を選択画素P2(1),P2(2)として選択する画素選択部133と、選択された2つの選択画素P2(1),P2(2)ごとに、2つの選択画素P2(1),P2(2)から所定位置に配置される複数の参照画素P3(n)それぞれの画素情報に基づいて、撮像画像i1を構成する画素の中から補正後の対象輪郭線Ltに対応する補正画素P4(n)を決定する補正画素決定部134と、補正画素P4(n)の位置情報に基づいて、補正後の対象輪郭線Ltを特定する補正輪郭線特定部136と、を有してなる。この構成によれば、本装置1は、2つの抽出画素P1(n)ごとに参照画素P3(n)の画素情報(輝度値の変化量)に基づいて1つの補正画素P4(n)を決定することができる。そして、本装置1は、決定された補正画素P4(n)の位置情報に基づいて、補正後の対象輪郭線Lt(輪郭線L)を特定することができる。すなわち、本装置1は、撮像画像i1から抽出された血管の輪郭線Lを補正することができる。
●Summary (1)
According to the embodiment described above, the apparatus 1 includes the contour line selection unit 131 that selects one of the first contour line L1 and the second contour line L2 as the correction target contour line Lt, and the contour pixel An extraction unit 132 for extracting a plurality of contour pixels as extraction pixels P1(n) from among the contour pixels corresponding to the target contour line Lt based on the position information of the extraction unit 132, and a predetermined A pixel selection unit 133 that selects two extracted pixels P1(n) having a positional relationship as selected pixels P2(1) and P2(2), and for each of the two selected selected pixels P2(1) and P2(2) Then, based on the pixel information of each of the plurality of reference pixels P3(n) arranged at predetermined positions from the two selected pixels P2(1) and P2(2), the pixels forming the captured image i1 after correction are selected. A correction pixel determination unit 134 that determines the correction pixel P4(n) corresponding to the target contour line Lt of the correction contour line that specifies the target contour line Lt after correction based on the position information of the correction pixel P4(n) and a specifying unit 136 . According to this configuration, the apparatus 1 determines one correction pixel P4(n) based on pixel information (amount of change in luminance value) of the reference pixel P3(n) for each two extracted pixels P1(n). can do. Then, the apparatus 1 can specify the corrected target contour line Lt (contour line L) based on the determined positional information of the correction pixel P4(n). That is, the apparatus 1 can correct the contour line L of the blood vessel extracted from the captured image i1.

また、以上説明した実施の形態によれば、画素情報は、参照画素P3(n)の輝度値を含む。また、補正画素決定部134は、複数の参照画素P3(n)それぞれの輝度値に基づいて、隣接する2つの参照画素P3(n)ごとの輝度値の変化量を取得し、複数の参照画素P3(n)のうち、変化量が最大の参照画素P3(n)を補正画素P4(n)として決定する。通常、撮像画像i1において、被写体である血管は、背景よりも暗く撮像される。そのため、画素値の変化量の大きい画素は、血管のエッジに相当する画素とみなすことができる。この構成によれば、本装置1は、容易に血管のエッジに相当する画素(関連すれば、補正後の輪郭線Lxの画素)を特定することができる。 Moreover, according to the embodiment described above, the pixel information includes the luminance value of the reference pixel P3(n). Further, the correction pixel determination unit 134 acquires the amount of change in luminance value for each of two adjacent reference pixels P3(n) based on the luminance values of each of the plurality of reference pixels P3(n), Among P3(n), the reference pixel P3(n) having the largest amount of change is determined as the corrected pixel P4(n). Usually, in the captured image i1, the blood vessel, which is the subject, is imaged darker than the background. Therefore, a pixel with a large amount of change in pixel value can be regarded as a pixel corresponding to the edge of a blood vessel. According to this configuration, the device 1 can easily identify the pixels corresponding to the edge of the blood vessel (the pixels of the contour line Lx after correction, if related).

さらに、以上説明した実施の形態によれば、補正画素決定部134は、第1輪郭線L1が選択されているとき、第2輪郭線L2を含む所定領域に含まれる画素を補正画素P4(n)の対象から除外し、第2輪郭線L2が選択されているとき、第1輪郭線L1を含む所定領域に含まれる画素を補正画素P4(n)の対象から除外する。この構成によれば、仮想垂線Lnにおいて輝度値の変化量を算出するとき、対象輪郭線Ltに選択されていない輪郭線L付近の輝度値の変化の影響は、受けない。その結果、輝度値の変化量が最大の画素は、単一の参照画素P3(n)となり、容易に特定することができる。 Furthermore, according to the embodiment described above, when the first contour line L1 is selected, the correction pixel determining unit 134 selects the pixels included in the predetermined area including the second contour line L2 as the correction pixels P4(n ), and when the second contour line L2 is selected, the pixels included in the predetermined region including the first contour line L1 are excluded from the correction pixel P4(n). According to this configuration, when calculating the amount of change in luminance value on the virtual perpendicular line Ln, the change in luminance value near the contour line L not selected as the target contour line Lt is not affected. As a result, the pixel with the largest change in luminance value is the single reference pixel P3(n), which can be easily identified.

さらにまた、以上説明した実施の形態によれば、補正画素決定部134は、第1輪郭線L1または第2輪郭線L2に膨張処理を実行することにより所定領域を画定する。この構成によれば、本装置1は、マスクする所定領域を容易に確定することができる。 Furthermore, according to the embodiment described above, the correction pixel determination unit 134 demarcates a predetermined area by performing expansion processing on the first contour line L1 or the second contour line L2. According to this configuration, the device 1 can easily determine the predetermined region to be masked.

さらにまた、以上説明した実施の形態によれば、本装置1は、補正画素P4(n)の位置情報を補正する2次補正部135を有する。2次補正部135は、複数の補正画素P4(n)それぞれに対応する曲率を算出し、補正画素P4(n)それぞれの曲率に基づいて、補正画素P4(n)の位置情報を補正する。この構成によれば、本装置1は、輝度値の変化量に大きな差異がなく、補正画素P4(n)の位置情報に誤りが生じ易い部分の補正画素P4(n)の位置情報を、さらに補正することができる。すなわち、本装置1は、高精度に輪郭線Lを補正することができる。また、本装置1は、曲率、すなわち、他の補正画素P4(n)に対する相対的な位置情報の差異に基づいて、位置情報を補正する補正画素P4(n)を特定することができる。 Furthermore, according to the embodiment described above, the device 1 has the secondary correction section 135 that corrects the position information of the correction pixel P4(n). The secondary correction unit 135 calculates the curvature corresponding to each of the plurality of corrected pixels P4(n), and corrects the position information of the corrected pixel P4(n) based on the curvature of each corrected pixel P4(n). According to this configuration, the apparatus 1 can obtain the position information of the correction pixel P4(n) in the portion where there is no large difference in the amount of change in the luminance value and the position information of the correction pixel P4(n) is likely to be erroneous. can be corrected. That is, the device 1 can correct the contour line L with high accuracy. Further, the apparatus 1 can specify the correction pixel P4(n) whose position information is corrected based on the curvature, that is, the difference in position information relative to other correction pixels P4(n).

さらにまた、以上説明した実施の形態によれば、2次補正部135は、算出された曲率に基づいて、補正の対象となる補正画素P4(n)を位置補正画素P6(n)として選択する。この構成によれば、本装置1は、曲率に基づいて、他の補正画素P4(n)に対する相対的な位置情報が大きくずれている補正画素P4(n)を特定し、同補正画素P4(n)の位置情報のみを補正することができる。また、2次補正部135は、位置補正画素P6(n)の両隣に配置される2つの補正画素P4(n)を基準画素P7(1),P7(2)として選択し、2つの基準画素P7(1),P7(2)の位置情報に基づいて、位置補正画素P6(n)の位置情報を補正する。この構成によれば、本装置1は、位置補正画素P6(n)よりは位置情報が血管の輪郭からずれておらず、かつ、位置補正画素P6(n)の正しい位置に近い画素(基準画素P7(1),P7(2))の位置情報に基づいて、位置補正画素P6(n)の位置情報を補正することができる。すなわち、本装置1は、撮像画像i1から抽出された血管の輪郭線Lを高精度に補正することができる。 Furthermore, according to the embodiment described above, the secondary correction unit 135 selects the correction pixel P4(n) to be corrected as the position correction pixel P6(n) based on the calculated curvature. . According to this configuration, the apparatus 1 identifies the corrected pixel P4(n) whose position information relative to the other corrected pixel P4(n) is greatly deviated based on the curvature, and determines the corrected pixel P4(n). Only the position information of n) can be corrected. Further, the secondary correction unit 135 selects two correction pixels P4(n) arranged on both sides of the position correction pixel P6(n) as reference pixels P7(1) and P7(2), and selects two reference pixels P7(1) and P7(2). The position information of the position correction pixel P6(n) is corrected based on the position information of P7(1) and P7(2). According to this configuration, the apparatus 1 can detect pixels (reference pixels) whose position information is not deviated from the contour of the blood vessel and which are closer to the correct positions of the position-corrected pixels P6(n) than the position-corrected pixels P6(n). The position information of the position correction pixel P6(n) can be corrected based on the position information of P7(1), P7(2)). That is, the apparatus 1 can highly accurately correct the contour line L of the blood vessel extracted from the captured image i1.

さらにまた、以上説明した実施の形態によれば、2次補正部135は、曲率が相対的に大きい補正画素P4(n)を位置補正画素P6(n)として選択する。この構成によれば、本装置1は、曲率に基づいて、他の補正画素P4(n)に対する相対的な位置情報が大きくずれている補正画素P4(n)のみを特定し、同補正画素P4(n)の位置情報のみを補正することができる。 Furthermore, according to the embodiment described above, the secondary correction unit 135 selects the correction pixel P4(n) having a relatively large curvature as the position correction pixel P6(n). According to this configuration, the apparatus 1 identifies only the corrected pixel P4(n) whose relative positional information is greatly deviated with respect to the other corrected pixel P4(n) based on the curvature. Only the position information of (n) can be corrected.

さらにまた、以上説明した実施の形態によれば、2次補正部135は、複数の補正画素P4(n)のうち、曲率を算出する補正画素P4(n)を算出画素P5(1)として選択し、算出画素P5(1)の両隣に配置される2つの補正画素P4(n)を算出補助画素P5(2),P5(3)として選択し、算出画素P5(1)と2つの算出補助画素P5(2),P5(3)とを仮想的に結ぶ仮想円Cの半径に基づいて、算出画素P5(1)の曲率を算出する。この構成によれば、本装置1は、補正画素P4(n)の位置情報のみに基づいて、容易に曲率を算出することができる。 Furthermore, according to the embodiment described above, the secondary correction unit 135 selects the correction pixel P4(n) whose curvature is to be calculated from among the plurality of correction pixels P4(n) as the calculation pixel P5(1). Then, the two correction pixels P4(n) arranged on both sides of the calculation pixel P5(1) are selected as the calculation auxiliary pixels P5(2) and P5(3), and the calculation pixel P5(1) and the two calculation auxiliary pixels are selected. The curvature of the calculated pixel P5(1) is calculated based on the radius of the virtual circle C that virtually connects the pixels P5(2) and P5(3). With this configuration, the device 1 can easily calculate the curvature based only on the position information of the corrected pixel P4(n).

●画像処理装置(2)●
次に、本装置の別の実施の形態(第2実施形態)について、先に説明した第1実施形態とは異なる部分を中心に説明する。第2実施形態の本装置は、制御部の画素選択部と補正画素決定部とが第1実施形態と異なる。
●Image processing device (2)●
Next, another embodiment (second embodiment) of the present apparatus will be described, focusing on portions different from the first embodiment described above. The apparatus of the second embodiment differs from the first embodiment in the pixel selection section and correction pixel determination section of the control section.

●画像処理装置(2)の構成
図12は、本装置の別の実施の形態を示す機能ブロック図である。
●Configuration of Image Processing Apparatus (2) FIG. 12 is a functional block diagram showing another embodiment of this apparatus.

本装置1Aは、通信部11と、記憶部12と、制御部13Aと、を備える。 The device 1A includes a communication section 11, a storage section 12, and a control section 13A.

制御部13Aは、本装置1A全体の動作を制御すると共に、後述する本方法を実行する。制御部13Aの物理的な構成は、第1実施形態の制御部13の物理的な構成と同じである。制御部13Aは、輪郭線選択部131と、抽出部132と、画素選択部133Aと、補正画素決定部134Aと、2次補正部135と、補正輪郭線特定部136と、を備える。 The control unit 13A controls the operation of the entire device 1A and executes the method described later. The physical configuration of the controller 13A is the same as the physical configuration of the controller 13 of the first embodiment. The control unit 13A includes a contour line selection unit 131, an extraction unit 132, a pixel selection unit 133A, a correction pixel determination unit 134A, a secondary correction unit 135, and a correction contour line identification unit 136.

画素選択部133Aは、複数の抽出画素のうち、所定の位置関係を有する2つの抽出画素を選択画素として選択する。画素選択部133Aの具体的な動作と所定の位置関係とは、後述する。 133 A of pixel selection parts select two extraction pixels which have a predetermined positional relationship as a selection pixel among several extraction pixels. A specific operation of the pixel selection unit 133A and the predetermined positional relationship will be described later.

補正画素決定部134Aは、選択された2つの選択画素ごとに、2つの選択画素から略同じ距離の位置に配置される複数の画素(参照画素)それぞれの画素情報に基づいて、撮像画像i1を構成する画素の中から補正後の対象輪郭線に対応する補正画素を決定する。補正画素決定部134Aの具体的な動作は、後述する。 The correction pixel determination unit 134A determines the captured image i1 based on the pixel information of each of a plurality of pixels (reference pixels) arranged at positions approximately the same distance from the two selected pixels for each of the two selected pixels. A correction pixel corresponding to the target contour line after correction is determined from among the constituent pixels. A specific operation of the correction pixel determination unit 134A will be described later.

●画像処理装置(2)の動作
次に、本装置1Aの動作、すなわち、本装置1Aが実行する本方法の別の実施の形態について、図12も参照しつつ説明する。
●Operation of Image Processing Apparatus (2) Next, the operation of the apparatus 1A, that is, another embodiment of the method executed by the apparatus 1A will be described with reference to FIG.

図13は、本方法の別の実施の形態を示すフローチャートである。 FIG. 13 is a flow chart illustrating another embodiment of the method.

本方法は、1次補正処理(S1A)と2次補正処理(S2)とを実行する。本実施の形態において、2次補正処理(S2)は、第1実施形態の2次補正処理(S2)と共通する。 The method performs a primary correction process (S1A) and a secondary correction process (S2). In this embodiment, the secondary correction process (S2) is common to the secondary correction process (S2) of the first embodiment.

●1次補正処理
先ず、本装置1Aは、1次補正処理(S1A)を実行する。
● Primary Correction Processing First, the device 1A executes the primary correction processing (S1A).

図14は、1次補正処理(S1A)のフローチャートである。
図15は、1次補正処理(S1A)の各段階の撮像画像i1を示す模式図である。
FIG. 14 is a flow chart of the primary correction process (S1A).
FIG. 15 is a schematic diagram showing captured images i1 at each stage of the primary correction processing (S1A).

「1次補正処理(S1A)」は、撮像画像i1を構成する画素の画素情報に基づいて、撮像画像i1に関連付けられている輪郭線Lに対応する輪郭画素の位置情報を補正する処理である。 The "primary correction process (S1A)" is a process of correcting the position information of the contour pixels corresponding to the contour line L associated with the captured image i1, based on the pixel information of the pixels forming the captured image i1. .

先ず、輪郭線選択部131は、第1輪郭線L1と第2輪郭線L2とのいずれか一方を対象輪郭線Ltとして選択する(S101)。 First, the contour line selection unit 131 selects one of the first contour line L1 and the second contour line L2 as the target contour line Lt (S101).

次いで、抽出部132は、撮像画像i1を構成する画素のうち、対象輪郭線Ltに対応する輪郭画素の位置情報に基づいて、対象輪郭線Ltに対応する輪郭画素のうち、複数の輪郭画素を抽出画素P1(n)として抽出する(S102A)。すなわち、例えば、抽出部132は、複数の輪郭画素から、略等間隔に一定数の輪郭画素を抽出画素P1(n)として抽出する。ここで、処理(S102A)は、第1実施形態の処理(S103)と共通する。以下の説明において、抽出部132は、6つの抽出画素P1(1)~P1(6)を抽出したものとする。 Next, the extraction unit 132 extracts a plurality of contour pixels out of the contour pixels corresponding to the target contour line Lt based on the position information of the contour pixels corresponding to the target contour line Lt among the pixels forming the captured image i1. It is extracted as an extraction pixel P1(n) (S102A). That is, for example, the extracting unit 132 extracts a certain number of contour pixels at approximately equal intervals as extracted pixels P1(n) from a plurality of contour pixels. Here, the process (S102A) is common to the process (S103) of the first embodiment. In the following description, it is assumed that the extraction unit 132 has extracted six extraction pixels P1(1) to P1(6).

次いで、画素選択部133Aは、複数の抽出画素P1(1)~P1(6)のうち、補正の候補となる抽出画素P1(n)を候補画素P8(n)として選択する(S103A)。本方法では、候補画素P8(n)は、始点となる抽出画素P1(1)と終点となる抽出画素P1(6)とを除く4つの抽出画素P1(2)~P1(5)の中から選択される。具体的には、画素選択部133Aは、抽出画素P1(2)~P1(5)のうち、始点に最も近く、かつ、候補画素P8(n)として選択されていない抽出画素P1(2)を候補画素P8(1)として選択する。 Next, the pixel selection unit 133A selects the extraction pixel P1(n), which is a candidate for correction, from among the plurality of extraction pixels P1(1) to P1(6) as the candidate pixel P8(n) (S103A). In this method, the candidate pixel P8(n) is selected from among the four extracted pixels P1(2) to P1(5) excluding the starting pixel P1(1) and the ending pixel P1(6). selected. Specifically, the pixel selection unit 133A selects the extracted pixel P1(2), which is closest to the starting point and is not selected as the candidate pixel P8(n), from among the extracted pixels P1(2) to P1(5). Select as candidate pixel P8(1).

次いで、画素選択部133Aは、候補画素P8(1)の両隣に配置される2つの抽出画素P1(1),P1(3)を選択画素P2A(1),P2A(2)として選択する(S104A)。すなわち、画素選択部133Aは、複数の抽出画素P1(1)~P1(6)のうち、候補画素P8(n)の両隣(1つ飛ばし)の位置関係となる2つの抽出画素P1(n),P1(n+2)を選択画素P2A(1),P2A(2)として選択する(S104A)。 Next, the pixel selection unit 133A selects two extracted pixels P1(1) and P1(3) arranged on both sides of the candidate pixel P8(1) as selected pixels P2A(1) and P2A(2) (S104A ). That is, the pixel selection unit 133A selects two extracted pixels P1(n) that have a positional relationship on both sides (skipping one) of the candidate pixel P8(n) among the plurality of extracted pixels P1(1) to P1(6). , P1(n+2) are selected as selected pixels P2A(1) and P2A(2) (S104A).

ここで、本実施の形態において、「候補画素P8(n)の両隣(1つ飛ばし)」の位置関係は、本発明における所定の位置関係の例である。 Here, in the present embodiment, the positional relationship between "two neighboring pixels (skipping one) of the candidate pixel P8(n)" is an example of a predetermined positional relationship in the present invention.

図16は、1次補正処理(S1A)の途中段階の撮像画像i1を示す模式図である。
同図は、説明の便宜上、2つの選択画素P2A(1),PA2(2)の近傍を拡大して示す。
FIG. 16 is a schematic diagram showing a captured image i1 in the middle of the primary correction process (S1A).
For convenience of explanation, this figure shows the vicinity of two selected pixels P2A(1) and PA2(2) in an enlarged manner.

図14~図16に戻る。
次いで、補正画素決定部134Aは、2つの選択画素P2A(1),P2A(2)同士を結ぶ仮想線分Lmを生成し、仮想線分Lmの中点において、仮想線分Lmに直交する仮想垂線Lnを生成する(S105)。ここで、仮想垂線Lnの長さは、例えば、所定の長さ(例えば、中点を中心画素とする9画素)に予め設定されている。
Returning to FIGS. 14-16.
Next, the correction pixel determination unit 134A generates a virtual line segment Lm that connects the two selected pixels P2A(1) and P2A(2), and creates a virtual line segment Lm that is orthogonal to the virtual line segment Lm at the midpoint of the virtual line segment Lm. A perpendicular line Ln is generated (S105). Here, the length of the virtual perpendicular line Ln is set in advance to, for example, a predetermined length (for example, 9 pixels with the midpoint as the central pixel).

次いで、補正画素決定部134Aは、撮像画像i1を構成する画素のうち、仮想垂線Lnに対応する複数の画素を参照画素P3(n)として特定する(S106)。すなわち、例えば、仮想垂線Lnの長さが9画素の場合、9つの画素を参照画素P3(1)~P3(9)として決定する。ここで、仮想垂線Lnに対応する画素は、2つの選択画素P2A(1),P2A(2)から略同じ距離の位置に配置される画素である。すなわち、仮想垂線Lnに対応する画素は、本発明における所定位置に配置される画素である。 Next, the correction pixel determination unit 134A identifies, among the pixels forming the captured image i1, a plurality of pixels corresponding to the virtual perpendicular line Ln as reference pixels P3(n) (S106). That is, for example, when the length of the virtual perpendicular line Ln is 9 pixels, the 9 pixels are determined as the reference pixels P3(1) to P3(9). Here, the pixels corresponding to the imaginary perpendicular line Ln are the pixels arranged at substantially the same distance from the two selected pixels P2A(1) and P2A(2). That is, the pixels corresponding to the imaginary perpendicular line Ln are pixels arranged at predetermined positions in the present invention.

次いで、補正画素決定部134Aは、記憶部12から参照画素P3(1)~P3(9)それぞれの輝度値を取得して、複数の参照画素P3(1)~P3(9)それぞれの輝度値に基づいて、隣接する2つの参照画素P3(n),P3(n+1)ごとの輝度値の変化量を取得する(S107)。すなわち、例えば、補正画素決定部134は、参照画素P3(1)~P3(9)それぞれの1次輝度勾配と2次輝度勾配とを算出する。 Next, the correction pixel determination unit 134A acquires the luminance values of the reference pixels P3(1) to P3(9) from the storage unit 12, and calculates the luminance values of the plurality of reference pixels P3(1) to P3(9). Based on (S107), the amount of change in luminance value for each of two adjacent reference pixels P3(n) and P3(n+1) is acquired. That is, for example, the correction pixel determination unit 134 calculates the primary luminance gradient and the secondary luminance gradient of each of the reference pixels P3(1) to P3(9).

次いで、補正画素決定部134Aは、複数の参照画素P3(1)~P3(9)のうち、対象輪郭線に近く、かつ、変化量が最大の参照画素P3(n)を補正画素P4(n)として決定する(S108)。 Next, the correction pixel determination unit 134A selects the reference pixel P3(n), which is close to the target contour line and has the largest change amount among the plurality of reference pixels P3(1) to P3(9), as the correction pixel P4(n). ) (S108).

なお、本発明における補正画素決定部は、第1実施形態の本方法と同様のマスク処理を実行してもよい。 Note that the correction pixel determination unit in the present invention may perform mask processing similar to the present method of the first embodiment.

次いで、制御部13A(たとえば、補正画素決定部134A)は、候補画素P8(n)として未選択の抽出画素P1(n)の有無を判定する(S109A)。未選択の抽出画素P1(n)が有るとき(S109Aの「Y」)、1次補正処理(S1A)は、処理(S103A)に戻る。すなわち、例えば、画素選択部133Aは、未選択の抽出画素P1(3)~P1(5)のうち、既に選択済みの抽出画素P1(2)に隣り合う抽出画素P1(3)を次の候補画素P8(2)として選択する。 Next, the control unit 13A (for example, the correction pixel determination unit 134A) determines whether or not there is an unselected extracted pixel P1(n) as the candidate pixel P8(n) (S109A). When there is an unselected extraction pixel P1(n) ("Y" in S109A), the primary correction process (S1A) returns to the process (S103A). That is, for example, the pixel selection unit 133A selects the extracted pixel P1(3) adjacent to the already selected extracted pixel P1(2) among the unselected extracted pixels P1(3) to P1(5) as the next candidate. Select as pixel P8(2).

一方、未選択の抽出画素P1(n)が無いとき(S109Aの「N」)、補正画素決定部134Aは、全ての候補画素P8(1)~P8(4)(すなわち、抽出画素P1(2)~P1(5))の位置情報を対応する補正画素P4(1)~P4(4)の位置情報に補正する(S110A)。すなわち、補正画素決定部134Aは、候補画素P8(n)に対応する抽出画素P1(n)の位置を補正画素P4(n)の位置に移動させる。 On the other hand, when there is no unselected extraction pixel P1(n) (“N” in S109A), the correction pixel determination unit 134A selects all the candidate pixels P8(1) to P8(4) (that is, the extraction pixel P1(2) ) to P1(5)) are corrected to corresponding position information of the corrected pixels P4(1) to P4(4) (S110A). That is, the correction pixel determination unit 134A moves the position of the extraction pixel P1(n) corresponding to the candidate pixel P8(n) to the position of the correction pixel P4(n).

次いで、補正画素決定部134Aは、全ての補正画素P4(2)~P4(5)の位置情報を、補正画素群の位置情報として撮像画像i1に関連付けて記憶部12に記憶する(S111A)。 Next, the correction pixel determining unit 134A stores the position information of all the correction pixels P4(2) to P4(5) in the storage unit 12 in association with the captured image i1 as the position information of the correction pixel group (S111A).

次いで、制御部13A(例えば、補正画素決定部134A)は、対象輪郭線Ltとして未選択の輪郭線Lの有無を判定する(S112A)。未選択の輪郭線Lが有るとき(S112Aの「Y」)、1次補正処理(S1A)は、処理(S101)に戻る。一方、未選択の輪郭線Lが無いとき(S112Aの「N」)、1次補正処理(S1A)は、終了する。ここで、処理(S112A)は、第1実施形態の処理(S111)と共通する。 Next, the control unit 13A (for example, the correction pixel determination unit 134A) determines whether or not there is an unselected contour line L as the target contour line Lt (S112A). When there is an unselected contour line L ("Y" in S112A), the primary correction process (S1A) returns to the process (S101). On the other hand, when there is no unselected contour line L ("N" in S112A), the primary correction process (S1A) ends. Here, the process (S112A) is common to the process (S111) of the first embodiment.

●まとめ(2)
以上説明した実施の形態によれば、画素選択部133は、複数の抽出画素P1(n)のうち、補正の候補となる抽出画素P1(n)を候補画素P8(n)として選択し、候補画素P8(n)の両隣に配置される2つの抽出画素P1(n)を選択画素P2A(1),P2A(2)として選択する。補正画素決定部134は、候補画素P8(n)の位置情報を補正画素P4(n)の位置情報に補正する。この構成によれば、本装置1は、抽出画素P1(n)の位置情報そのものを補正することができる。
●Summary (2)
According to the embodiment described above, the pixel selection unit 133 selects the extracted pixel P1(n), which is a candidate for correction, from among the plurality of extracted pixels P1(n) as the candidate pixel P8(n). Two extracted pixels P1(n) arranged on both sides of the pixel P8(n) are selected as selected pixels P2A(1) and P2A(2). The correction pixel determination unit 134 corrects the position information of the candidate pixel P8(n) to the position information of the correction pixel P4(n). According to this configuration, the apparatus 1 can correct the position information itself of the extracted pixel P1(n).

●まとめ(その他)
なお、以上説明した各実施の形態では、本装置1,1Aは、2次補正処理(S2)を実行していた。これに代えて、本装置は、1次補正処理において輪郭線を所望の精度で補正可能であれば、2次補正処理を実行しなくてもよい。
●Summary (Others)
In each of the embodiments described above, the apparatus 1, 1A executes secondary correction processing (S2). Alternatively, if the contour line can be corrected with desired accuracy in the primary correction process, the apparatus does not need to perform the secondary correction process.

また、本装置は、1次補正処理において、第1実施形態の選択画素の選択と、第2実施形態の選択画素の選択と、を組み合わせてもよい。すなわち、例えば、抽出画素P1(1)~P1(6)に対して、本装置は、始点に対応する抽出画素P1(1)と、同抽出画素P1(1)に隣り合う抽出画素P1(2)と、を選択画素として選択し、次いで、抽出画素P1(2)と1つ飛ばしの抽出画素P1(4)とを選択画素として選択してもよい。 Further, the present apparatus may combine the selection of the selected pixels in the first embodiment and the selection of the selected pixels in the second embodiment in the primary correction processing. That is, for example, for the extracted pixels P1(1) to P1(6), the present apparatus extracts the extracted pixel P1(1) corresponding to the starting point and the extracted pixel P1(2) adjacent to the extracted pixel P1(1). ) and are selected as the selected pixels, and then the extracted pixel P1(2) and the extracted pixel P1(4) skipping by one may be selected as the selected pixels.

1 画像処理装置
12 記憶部
131 輪郭線選択部
132 抽出部
133 画素選択部
134 補正画素決定部
135 2次補正部
136 補正輪郭線特定部
1A 画像処理装置
133A 画素選択部
134A 補正画素決定部
1 image processing device 12 storage unit 131 contour line selection unit 132 extraction unit 133 pixel selection unit 134 correction pixel determination unit 135 secondary correction unit 136 correction contour line identification unit 1A image processing device 133A pixel selection unit 134A correction pixel determination unit

Claims (15)

血管が撮像された撮像画像に基づいて生成された前記血管の輪郭線を補正する画像処理装置であって、
前記輪郭線は、
前記血管の幅方向における一方側の輪郭を示す第1輪郭線と、
前記幅方向における他方の輪郭を示す第2輪郭線と、
を含み、
前記撮像画像と、前記撮像画像に撮像された前記血管の前記輪郭線に対応する輪郭画素の位置情報と、を関連付けて記憶する記憶部と、
前記第1輪郭線と前記第2輪郭線とのいずれか一方を補正の対象輪郭線として選択する輪郭線選択部と、
前記輪郭画素の位置情報に基づいて、前記対象輪郭線に対応する前記輪郭画素のうち、複数の前記輪郭画素を抽出画素として抽出する抽出部と、
複数の前記抽出画素のうち、所定の位置関係を有する2つの前記抽出画素を選択画素として選択する画素選択部と、
選択された2つの前記選択画素ごとに、2つの前記選択画素から所定位置に配置される複数の参照画素それぞれの画素情報に基づいて、前記撮像画像を構成する画素の中から補正後の前記対象輪郭線に対応する補正画素を決定する補正画素決定部と、
前記補正画素の位置情報に基づいて、補正後の前記対象輪郭線を特定する補正輪郭線特定部と、
を有してなる、
ことを特徴とする画像処理装置。
An image processing device for correcting a contour line of a blood vessel generated based on a captured image in which the blood vessel is captured,
The outline is
a first contour line indicating the contour of one side in the width direction of the blood vessel;
a second contour line indicating the other contour in the width direction;
including
a storage unit that associates and stores the captured image and position information of contour pixels corresponding to the contour line of the blood vessel captured in the captured image;
a contour selection unit that selects one of the first contour and the second contour as a contour to be corrected;
an extraction unit that extracts a plurality of the contour pixels as extracted pixels from the contour pixels corresponding to the target contour based on the position information of the contour pixels;
a pixel selection unit that selects, as selected pixels, two of the extracted pixels having a predetermined positional relationship from among the plurality of extracted pixels;
For each of the selected two selected pixels, the target after correction from among pixels constituting the captured image based on pixel information of each of a plurality of reference pixels arranged at predetermined positions from the two selected pixels. a correction pixel determination unit that determines correction pixels corresponding to the contour;
a corrected contour identifying unit that identifies the target contour after correction based on the position information of the corrected pixel;
comprising
An image processing apparatus characterized by:
前記画素情報は、
前記参照画素の輝度値、
を含み、
前記補正画素決定部は、
複数の前記参照画素それぞれの前記輝度値に基づいて、隣接する2つの前記参照画素ごとの前記輝度値の変化量を取得し、
複数の前記参照画素のうち、前記変化量が最大の前記参照画素を前記補正画素として決定する、
請求項1記載の画像処理装置。
The pixel information is
a luminance value of the reference pixel;
including
The correction pixel determination unit
Obtaining an amount of change in the luminance value for each of two adjacent reference pixels based on the luminance value of each of the plurality of reference pixels;
determining the reference pixel having the largest amount of change among the plurality of reference pixels as the correction pixel;
The image processing apparatus according to claim 1.
前記画素選択部は、
複数の前記抽出画素のうち、隣り合う2つの前記抽出画素を前記選択画素として選択する、
請求項2記載の画像処理装置。
The pixel selection unit
selecting two adjacent extracted pixels from among the plurality of extracted pixels as the selected pixels;
3. The image processing apparatus according to claim 2.
前記画素選択部は、
複数の前記抽出画素のうち、補正の候補となる前記抽出画素を候補画素として選択し、前記候補画素の両隣に配置される2つの前記抽出画素を前記選択画素として選択する、
請求項2記載の画像処理装置。
The pixel selection unit
selecting the extracted pixel that is a candidate for correction from among the plurality of extracted pixels as a candidate pixel, and selecting two of the extracted pixels arranged on both sides of the candidate pixel as the selected pixels;
3. The image processing apparatus according to claim 2.
前記補正画素決定部は、前記候補画素の位置情報を前記補正画素の位置情報に補正する、
請求項4記載の画像処理装置。
The correction pixel determination unit corrects the position information of the candidate pixel to the position information of the correction pixel.
5. The image processing apparatus according to claim 4.
前記補正画素決定部は、
前記第1輪郭線が選択されているとき、前記第2輪郭線を含む所定領域に含まれる画素を前記補正画素の対象から除外し、
前記第2輪郭線が選択されているとき、前記第1輪郭線を含む所定領域に含まれる画素を前記補正画素の対象から除外する、
請求項1記載の画像処理装置。
The correction pixel determination unit
when the first contour line is selected, excluding pixels included in a predetermined region including the second contour line from the target of the correction pixels;
when the second contour is selected, excluding pixels included in a predetermined region including the first contour from the target of the correction pixels;
The image processing apparatus according to claim 1.
前記補正画素の位置情報を補正する2次補正部、
を有してなり、
前記2次補正部は、
複数の前記補正画素それぞれに対応する曲率を算出し、
前記補正画素それぞれの前記曲率に基づいて、前記補正画素の位置情報を補正する、
請求項1乃至6のいずれか記載の画像処理装置。
a secondary correction unit that corrects the position information of the correction pixel;
and
The secondary correction unit is
calculating a curvature corresponding to each of the plurality of correction pixels;
correcting the position information of the correction pixel based on the curvature of each of the correction pixels;
7. The image processing apparatus according to claim 1.
前記2次補正部は、複数の前記補正画素のうち、補正後の前記対象輪郭線の始点および終点に対応する前記補正画素を除く、前記補正画素それぞれに対応する前記曲率を算出する、
請求項7記載の画像処理装置。
The secondary correction unit calculates the curvature corresponding to each of the corrected pixels excluding the corrected pixels corresponding to the start point and the end point of the target contour line after correction among the plurality of corrected pixels,
The image processing apparatus according to claim 7.
前記2次補正部は、
算出された前記曲率に基づいて、補正の対象となる前記補正画素を位置補正画素として選択し、
前記位置補正画素の両隣に配置される2つの前記補正画素を基準画素として選択し、
2つの前記基準画素の位置情報に基づいて、前記位置補正画素の位置情報を補正する、
請求項8記載の画像処理装置。
The secondary correction unit is
selecting the correction pixel to be corrected as a position correction pixel based on the calculated curvature;
selecting two of the correction pixels arranged on both sides of the position correction pixel as reference pixels;
correcting the position information of the position correction pixel based on the position information of the two reference pixels;
9. The image processing apparatus according to claim 8.
前記2次補正部は、前記曲率が相対的に大きい前記補正画素を前記位置補正画素として選択する、
請求項9記載の画像処理装置。
The secondary correction unit selects the correction pixel having the relatively large curvature as the position correction pixel.
10. The image processing apparatus according to claim 9.
前記2次補正部は、2つの前記基準画素の中点に対応する画素の位置情報に前記位置補正画素の位置情報を補正する、
請求項9記載の画像処理装置。
The secondary correction unit corrects the position information of the position correction pixel to the position information of the pixel corresponding to the midpoint of the two reference pixels.
10. The image processing apparatus according to claim 9.
前記2次補正部は、
複数の前記補正画素のうち、前記曲率を算出する前記補正画素を算出画素として選択し、
前記算出画素の両隣に配置される2つの前記補正画素を算出補助画素として選択し、
前記算出画素と2つの前記算出補助画素とを仮想的に結ぶ仮想円の半径に基づいて、前記算出画素の前記曲率を算出する、
請求項7記載の画像処理装置。
The secondary correction unit is
Selecting the correction pixel for calculating the curvature as a calculation pixel from among the plurality of correction pixels;
selecting the two correction pixels arranged on both sides of the calculation pixel as calculation auxiliary pixels;
calculating the curvature of the calculated pixel based on the radius of a virtual circle that virtually connects the calculated pixel and two of the calculated auxiliary pixels;
The image processing apparatus according to claim 7.
コンピュータを請求項1記載の画像処理装置として機能させる、
ことを特徴とする画像処理プログラム。
causing a computer to function as the image processing device according to claim 1;
An image processing program characterized by:
血管が撮像された撮像画像に基づいて生成された前記血管の輪郭線を補正する画像処理装置により実行される画像処理方法であって、
前記輪郭線は、
前記血管の幅方向における一方側の輪郭を示す第1輪郭線と、
前記幅方向における他方の輪郭を示す第2輪郭線と、
を含み、
前記画像処理装置は、
前記撮像画像と、前記撮像画像に撮像された前記血管の前記輪郭線に対応する輪郭画素の位置情報と、を関連付けて記憶する記憶部、
を備え、
前記画像処理装置が、
前記第1輪郭線と前記第2輪郭線とのいずれか一方を補正の対象輪郭線として選択するステップと、
前記輪郭画素の位置情報に基づいて、前記対象輪郭線に対応する前記輪郭画素のうち、複数の前記輪郭画素を抽出画素として抽出するステップと、
複数の前記抽出画素のうち、所定の位置関係を有する2つの前記抽出画素を選択画素として選択するステップと、
選択された2つの前記選択画素ごとに、2つの前記選択画素から所定位置に配置される複数の参照画素それぞれの画素情報に基づいて、前記撮像画像を構成する画素の中から補正後の前記対象輪郭線に対応する補正画素を決定するステップと、
前記補正画素の位置情報に基づいて、補正後の前記対象輪郭線を特定するステップと、
を含む、
ことを特徴とする画像処理方法。
An image processing method executed by an image processing device for correcting a contour line of a blood vessel generated based on a captured image in which the blood vessel is captured,
The outline is
a first contour line indicating a contour on one side in the width direction of the blood vessel;
a second contour line indicating the other contour in the width direction;
including
The image processing device is
a storage unit that associates and stores the captured image and position information of contour pixels corresponding to the contour line of the blood vessel captured in the captured image;
with
The image processing device is
selecting either one of the first contour and the second contour as a contour to be corrected;
a step of extracting a plurality of the contour pixels as extracted pixels from the contour pixels corresponding to the target contour line based on the positional information of the contour pixels;
a step of selecting two extracted pixels having a predetermined positional relationship from among the plurality of extracted pixels as selected pixels;
For each of the selected two selected pixels, the target after correction from among pixels constituting the captured image based on pixel information of each of a plurality of reference pixels arranged at predetermined positions from the two selected pixels. determining a correction pixel corresponding to the contour;
a step of identifying the corrected target contour line based on the position information of the corrected pixel;
including,
An image processing method characterized by:
前記補正画素の位置情報を補正するステップ、
を含む、
請求項14記載の画像処理方法。
correcting the position information of the correction pixel;
including,
15. The image processing method according to claim 14.
JP2021027339A 2021-02-24 2021-02-24 Image processing device, image processing program, and image processing method Pending JP2022128882A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021027339A JP2022128882A (en) 2021-02-24 2021-02-24 Image processing device, image processing program, and image processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021027339A JP2022128882A (en) 2021-02-24 2021-02-24 Image processing device, image processing program, and image processing method

Publications (1)

Publication Number Publication Date
JP2022128882A true JP2022128882A (en) 2022-09-05

Family

ID=83150621

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021027339A Pending JP2022128882A (en) 2021-02-24 2021-02-24 Image processing device, image processing program, and image processing method

Country Status (1)

Country Link
JP (1) JP2022128882A (en)

Similar Documents

Publication Publication Date Title
US6674883B1 (en) System and method for the detection of anatomic landmarks for total hip replacement
US10950204B2 (en) Diagnosis support apparatus and diagnosis support method
CN110766730B (en) Image registration and follow-up evaluation method, storage medium and computer equipment
JP2008259622A (en) Report writing supporting apparatus and its program
JP5812841B2 (en) Image diagnosis support apparatus, image diagnosis support method, and image diagnosis support program
US9198603B2 (en) Device, method and program for searching for the shortest path in a tubular structure
CN111311655B (en) Multi-mode image registration method, device, electronic equipment and storage medium
US20220207771A1 (en) Heart Position Estimation
CN111932492B (en) Medical image processing method and device and computer readable storage medium
JP2008079770A (en) Diagnostic reading support apparatus, diagnostic reading support method and program
KR102228087B1 (en) Method and apparatus for segmentation of specific cartilage in medical image
US9846939B2 (en) Image display apparatus and image display method
KR20190105784A (en) Method for diagnosing scoliosis using spatial coordinates of body shape and computer program therefor
KR20190080705A (en) Program and method for providing feedback about result of surgery
US20210068695A1 (en) Method Providing ECG Analysis Interface and System
CN114842003A (en) Medical image follow-up target pairing method, device and application
CN109350059A (en) For ancon self-aligning combined steering engine and boundary mark engine
Reddy et al. Anatomical Landmark Detection using Deep Appearance-Context Network
KR102216455B1 (en) Apparatus and method for measuring left atrium thickness of heart
JP2022128882A (en) Image processing device, image processing program, and image processing method
US20220249174A1 (en) Surgical navigation system, information processing device and information processing method
CN115700740A (en) Medical image processing method, apparatus, computer device and storage medium
US10299864B1 (en) Co-localization of multiple internal organs based on images obtained during surgery
KR102655362B1 (en) Method and apparatus for stitching medical images
JP2008061847A (en) Image pickup device and image analysis program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240118