WO2022109907A1 - Method and storage medium for precise acquisition of stenotic lesion range - Google Patents

Method and storage medium for precise acquisition of stenotic lesion range Download PDF

Info

Publication number
WO2022109907A1
WO2022109907A1 PCT/CN2020/131703 CN2020131703W WO2022109907A1 WO 2022109907 A1 WO2022109907 A1 WO 2022109907A1 CN 2020131703 W CN2020131703 W CN 2020131703W WO 2022109907 A1 WO2022109907 A1 WO 2022109907A1
Authority
WO
WIPO (PCT)
Prior art keywords
blood vessel
interval
pipe diameter
stenotic lesion
curve
Prior art date
Application number
PCT/CN2020/131703
Other languages
French (fr)
Chinese (zh)
Inventor
王鹏
刘广志
王之元
徐磊
Original Assignee
苏州润迈德医疗科技有限公司
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 苏州润迈德医疗科技有限公司 filed Critical 苏州润迈德医疗科技有限公司
Publication of WO2022109907A1 publication Critical patent/WO2022109907A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/40ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Definitions

  • the invention relates to the technical field of coronary medicine, in particular to a method and a storage medium for accurately acquiring a stenotic lesion interval.
  • FFR is one of the coronary vascular evaluation parameters, and the microcirculation resistance index IMR belongs to the coronary vascular evaluation parameters.
  • the existing technology can obtain the stenotic lesion interval through different methods, all of them are by defining the stenotic position, extending from the stenotic position to both ends, and extending to the preset parameters to obtain the starting position and ending position of the stenotic lesion, thereby obtaining the stenotic lesion. interval.
  • the fixed preset parameters are unable to cope with the diversity of blood vessels, resulting in insufficient intensive reading of the stenotic lesion interval.
  • the present invention provides a method and a storage medium for accurately acquiring stenotic lesion intervals, so as to solve the problem of insufficient precision of stenotic lesion intervals due to the inability of fixed preset parameters to deal with the diversity of blood vessels.
  • the present application provides a method for accurately obtaining a narrow lesion interval, including:
  • the misjudged stenotic region is removed from the first stenotic lesion interval to obtain a second stenotic lesion interval.
  • the preliminary judgment, the method for obtaining the first narrow lesion interval includes:
  • the first stenosis lesion interval is obtained.
  • the method for fitting a normal blood vessel diameter and obtaining a fitting diameter curve includes:
  • i represents the curve sampling point of the ith pipe diameter
  • n represents the sum of the number of pipe diameter curve sampling
  • xi represents the length of the curve sampling point of the ith pipe diameter
  • y i represents the pipe diameter at xi ;
  • the method for obtaining the first stenosis position according to the fitted diameter curve and the real diameter curve includes:
  • the intersection is the first entry point of the narrow area, otherwise, the intersection is the first exit of the narrow area point;
  • the curve between the first entry point and the first exit point is the preliminarily determined stenosis position, that is, the first stenotic lesion interval.
  • the method for removing a misjudged narrow lesion area from the first narrow lesion interval to obtain a second narrow lesion interval includes;
  • the misjudged stenosis region is removed from the first stenotic lesion interval according to the stenosis degree and the length of the blood vessel center line to obtain a second stenotic lesion interval.
  • the method for calculating a stenosis degree includes:
  • A represents the stenosis degree of the blood vessel
  • D min represents the minimum diameter of the blood vessel between the first entry point and the first exit point
  • D in and D out represent the vessel diameter of the first entry point and the first exit point, respectively vascular diameter.
  • the method for obtaining a second stenotic lesion interval by removing a misjudged stenotic area from the first stenotic lesion interval according to the stenosis degree and the length of the blood vessel centerline include:
  • the stenotic area obtained after removing the misjudged area is the second stenotic lesion interval.
  • the above-mentioned method for accurately obtaining a narrow lesion interval further includes:
  • the two points where the fitted pipe diameter interval curve and the real pipe diameter curve intersect are a second entry point and a second exit point, and the second entry point and the second The interval between the exit points is the third stenotic lesion interval.
  • the above-mentioned method for accurately obtaining the stenotic lesion interval includes:
  • 3D modeling is carried out according to the real-time diameter D t of the blood vessel, the length L of the blood vessel center line and the stenosis interval to form a 3D vessel model with stenotic lesion interval;
  • N-edge meshing is performed along the circumference of the three-dimensional blood vessel model with stenotic lesions to form a single-layer mesh model, where N ⁇ 6;
  • the surface layering process is performed on the single-layer grid model to form a double-layer grid model, that is, a blood vessel mathematical model.
  • N-edge meshing is performed along the circumference of the three-dimensional blood vessel model with stenotic lesion intervals to form a single-layer grid model, where N ⁇ 6 methods include:
  • every N triangle combination is converted into an N-sided shape to form an N-sided initial mesh
  • the method for meshing along the circumference of the three-dimensional blood vessel model with stenotic lesion intervals using triangles as the smallest unit includes:
  • the method for performing surface layering processing on the single-layer grid model to form a double-layer grid model includes:
  • Three-dimensional modeling is performed according to the blood vessel wall thickness h, the blood vessel starting diameter D, the blood vessel ending diameter D and the blood vessel centerline length L, and a three-dimensional truncated truncated model is formed on the inner surface or outer surface of the single-layer mesh model;
  • N-edge grid division is performed along the circumference of the three-dimensional truncated truncated model to form another single-layer grid model;
  • Two layers of the single-layer grid model and the blood vessel wall thickness h form the double-layer grid model, that is, the blood vessel mathematical model.
  • the present application provides a computer storage medium, and when the computer program is executed by a processor, the above-mentioned method for accurately acquiring a narrow lesion interval can be implemented.
  • the present application provides a method and a storage medium for accurately obtaining a stenosis lesion interval, and further corrects the vascular stenosis with coarse precision by refitting the three-dimensional modeling of the vessel wall to ensure the accuracy of the stenosis lesion interval.
  • FIG. 1 is a flowchart of an embodiment of a method for accurately obtaining a narrow lesion interval according to the present application
  • Fig. 2 is the flow chart of S100 of this application.
  • Fig. 5 is the flow chart of S200 of this application.
  • FIG. 7 is a flowchart of S230 of the application.
  • FIG. 10 is a flowchart of another embodiment of the method for accurately obtaining a stenotic lesion interval according to the present application.
  • the present application provides a method for accurately obtaining a narrow lesion interval, including:
  • N-edge meshing is performed along the circumference of the three-dimensional blood vessel model with stenotic lesions to form a single-layer mesh model, where N ⁇ 6, including:
  • S121 along the circumference of the three-dimensional blood vessel model with the stenotic lesion interval, perform grid division with triangles as the smallest unit, including: dividing the three-dimensional blood vessel model with the narrow lesion interval into K segments, and in each segment of the three-dimensional blood vessel model On the circumferential surface, the triangle is used as the smallest element for mesh division.
  • every N triangle combination is converted into an N-sided shape to form an N-sided initial mesh
  • the surface layering process is performed on the single-layer grid model to form a double-layer grid model, that is, a blood vessel mathematical model, including:
  • the two-layer single-layer grid model and the blood vessel wall thickness h form a double-layer grid model, that is, a blood vessel mathematical model.
  • i represents the curve sampling point of the ith pipe diameter
  • n represents the sum of the number of pipe diameter curve sampling
  • xi represents the length of the curve sampling point of the ith pipe diameter
  • y i represents the pipe diameter at xi ;
  • the curve between the first entry point and the first exit point is the initially determined stenosis position, that is, the first stenotic lesion interval.
  • the misjudged stenosis area is removed from the first stenotic lesion interval to obtain a second stenotic lesion interval, including:
  • A represents the stenosis degree of the blood vessel
  • D min represents the minimum diameter of the blood vessel between the first entry point and the first exit point
  • D in and D out represent the vessel diameter of the first entry point and the first exit point, respectively vascular diameter.
  • the misjudged stenosis area is removed from the first stenotic lesion interval, and the second stenotic lesion interval is obtained, including:
  • the stenotic area obtained again after removing the misjudged area is the second stenotic lesion interval.
  • a method for accurately obtaining a narrow lesion interval provided by the present application further includes:
  • the present application provides a computer storage medium, and when the computer program is executed by a processor, the above-mentioned method for accurately acquiring a narrow lesion interval can be implemented.
  • aspects of the present invention may be embodied as a system, method or computer program product. Accordingly, various aspects of the present invention may be embodied in the form of an entirely hardware implementation, an entirely software implementation (including firmware, resident software, microcode, etc.), or a combination of hardware and software aspects, It may be collectively referred to herein as a "circuit,” "module,” or “system.” Furthermore, in some embodiments, various aspects of the present invention may also be implemented in the form of a computer program product on one or more computer-readable media having computer-readable program code embodied thereon. Implementation of the method and/or system of embodiments of the invention may involve performing or completing selected tasks manually, automatically, or a combination thereof.
  • a data processor such as a computing platform for executing a plurality of instructions.
  • the data processor includes volatile storage for storing instructions and/or data and/or non-volatile storage for storing instructions and/or data, such as a magnetic hard disk and/or a Move media.
  • a network connection is also provided.
  • a display and/or user input device such as a keyboard or mouse, is optionally also provided.
  • the computer-readable medium may be a computer-readable signal medium or a computer-readable storage medium.
  • the computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus or device, or a combination of any of the above. More specific examples (non-exhaustive list) of computer-readable storage media would include the following:
  • a computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.
  • a computer-readable signal medium may include a propagated data signal in baseband or as part of a carrier wave, with computer-readable program code embodied thereon. Such propagated data signals may take a variety of forms including, but not limited to, electromagnetic signals, optical signals, or any suitable combination of the foregoing.
  • a computer-readable signal medium can also be any computer-readable medium other than a computer-readable storage medium that can transmit, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device .
  • Program code embodied on a computer-readable medium may be transmitted using any suitable medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
  • computer program code for performing operations for various aspects of the present invention may be written in any combination of one or more programming languages, including object-oriented programming languages such as Java, Smalltalk, C++, and conventional procedural programming languages, such as The "C" programming language or similar programming language.
  • the program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer, or entirely on the remote computer or server.
  • the remote computer may be connected to the user's computer through any kind of network - including a local area network (LAN) or wide area network (WAN) - or may be connected to an external computer (eg using an Internet service provider via Internet connection).
  • LAN local area network
  • WAN wide area network
  • These computer program instructions can also be stored on a computer readable medium, the instructions cause a computer, other programmable data processing apparatus, or other device to operate in a particular manner, whereby the instructions stored on the computer readable medium produce the An article of manufacture of instructions implementing the functions/acts specified in one or more blocks of the flowcharts and/or block diagrams.
  • Computer program instructions can also be loaded on a computer (eg, a coronary artery analysis system) or other programmable data processing device to cause a series of operational steps to be performed on the computer, other programmable data processing device or other device to produce a computer-implemented process , such that instructions executing on a computer, other programmable apparatus, or other device provide a process for implementing the functions/acts specified in the flowchart and/or one or more block diagram blocks.
  • a computer eg, a coronary artery analysis system
  • other programmable data processing device to produce a computer-implemented process , such that instructions executing on a computer, other programmable apparatus, or other device provide a process for implementing the functions/acts specified in the flowchart and/or one or more block diagram blocks.

Abstract

Provided are a method and storage medium for the precise acquisition of a stenotic lesion range, comprising: performing three-dimensional modeling according to a real-time vessel diameter Dt, a vessel centerline length L, and a stenosis range to form a mathematical model having a narrow lesion range (S100); preliminary determination to obtain a first stenotic lesion range (S200); removing erroneously determined stenotic regions from the first stenotic lesion range to obtain a second stenotic lesion range (S300). In the method and storage medium for the precise acquisition of a stenotic lesion range, by means of modeling a vascular wall in three dimensions after refitting, and further correcting for coarse-precision vascular stenosis, the accuracy of a narrow lesion range is ensured.

Description

精确获取狭窄病变区间的方法及存储介质Method and storage medium for accurately obtaining stenotic lesion interval 技术领域technical field
本发明涉及冠状动脉医学技术领域,特别是涉及一种精确获取狭窄病变区间的方法及存储介质。The invention relates to the technical field of coronary medicine, in particular to a method and a storage medium for accurately acquiring a stenotic lesion interval.
背景技术Background technique
人体血液中的脂类及糖类物质在血管壁上的沉积将在血管壁上形成斑块,继而导致血管狭窄;特别是发生在心脏冠脉附近的血管狭窄将导致心肌供血不足,诱发冠心病、心绞痛等病症,对人类的健康造成严重威胁。据统计,我国现有冠心病患者约1100万人,心血管介入手术治疗患者数量每年增长大于10%。The deposition of lipids and carbohydrates in human blood on the blood vessel wall will form plaques on the blood vessel wall, which will then lead to vascular stenosis; especially the vascular stenosis near the coronary arteries of the heart will lead to insufficient blood supply to the heart muscle and induce coronary heart disease. , angina pectoris and other diseases pose a serious threat to human health. According to statistics, there are about 11 million coronary heart disease patients in my country, and the number of patients treated with cardiovascular interventional surgery is increasing by more than 10% every year.
冠脉造影CAG、计算机断层扫描CT等常规医用检测手段虽然可以显示心脏冠脉血管狭窄的严重程度,但是并不能准确评价冠脉的缺血情况。为提高冠脉血管功能评价的准确性,1993年Pijls提出了通过压力测定推算冠脉血管功能的新指标——血流储备分数(Fractional Flow Reserve,FFR),经过长期的基础与临床研究,FFR已成为冠脉狭窄功能性评价的金标准。Although conventional medical detection methods such as coronary angiography (CAG) and computed tomography (CT) can display the severity of coronary artery stenosis, they cannot accurately evaluate coronary ischemia. In order to improve the accuracy of coronary vascular function evaluation, in 1993, Pijls proposed a new index for calculating coronary vascular function through pressure measurement - Fractional Flow Reserve (FFR). After long-term basic and clinical research, FFR It has become the gold standard for functional evaluation of coronary stenosis.
FFR是冠状动脉血管评定参数的一种,微循环阻力指数IMR等属于冠状动脉血管评定参数。FFR is one of the coronary vascular evaluation parameters, and the microcirculation resistance index IMR belongs to the coronary vascular evaluation parameters.
现有技术虽然能通过不同方法获取狭窄病变区间,但均是通过定义狭窄位置,从狭窄位置向两端扩展,扩展到预设参数为止,得到狭窄病变起始位置与结束位置,从而得到狭窄病变区间。而在实际情况中,固定的预设参数是无法应对血管的多样性,导致狭窄病变区间精读不足的问题。Although the existing technology can obtain the stenotic lesion interval through different methods, all of them are by defining the stenotic position, extending from the stenotic position to both ends, and extending to the preset parameters to obtain the starting position and ending position of the stenotic lesion, thereby obtaining the stenotic lesion. interval. In practice, the fixed preset parameters are unable to cope with the diversity of blood vessels, resulting in insufficient intensive reading of the stenotic lesion interval.
发明内容SUMMARY OF THE INVENTION
本发明提供了一种精确获取狭窄病变区间的方法及存储介质,以解决现有 由于固定的预设参数是无法应对血管的多样性,导致狭窄病变区间精度不足的问题。The present invention provides a method and a storage medium for accurately acquiring stenotic lesion intervals, so as to solve the problem of insufficient precision of stenotic lesion intervals due to the inability of fixed preset parameters to deal with the diversity of blood vessels.
为实现上述目的,第一方面,本申请提供了一种精确获取狭窄病变区间的方法,包括:In order to achieve the above object, in a first aspect, the present application provides a method for accurately obtaining a narrow lesion interval, including:
根据血管实时直径D t、血管中心线长度L和狭窄区间进行三维建模,形成具有狭窄病变区间的数学模型; According to the real-time diameter D t of the blood vessel, the length L of the center line of the blood vessel and the stenosis interval, three-dimensional modeling is performed to form a mathematical model with stenotic lesion interval;
初步判断,获得第一狭窄病变区间;Preliminary judgment to obtain the first stenotic lesion interval;
从所述第一狭窄病变区间去除误判的狭窄区域,获得第二狭窄病变区间。The misjudged stenotic region is removed from the first stenotic lesion interval to obtain a second stenotic lesion interval.
可选地是,上述的精确获取狭窄病变区间的方法,所述初步判断,获得第一狭窄病变区间的方法,包括:Optionally, in the above-mentioned method for accurately obtaining a narrow lesion interval, the preliminary judgment, the method for obtaining the first narrow lesion interval, includes:
拟合正常血管管径,获取拟合管径曲线;Fit the normal blood vessel diameter and obtain the fitted diameter curve;
根据所述数学模型获取真实管径曲线;Obtain the real pipe diameter curve according to the mathematical model;
根据所述拟合管径曲线和真实管径曲线,获取第一狭窄病变区间。According to the fitted diameter curve and the real diameter curve, the first stenosis lesion interval is obtained.
可选地是,上述的精确获取狭窄病变区间的方法,所述拟合正常血管管径,获取拟合管径曲线的方法包括:Optionally, in the above-mentioned method for accurately obtaining a stenotic lesion interval, the method for fitting a normal blood vessel diameter and obtaining a fitting diameter curve includes:
根据拟合代价函数,获取拟合管径,具体公式为
Figure PCTCN2020131703-appb-000001
According to the fitting cost function, the fitted pipe diameter is obtained, and the specific formula is:
Figure PCTCN2020131703-appb-000001
其中,i表示第i个管径的曲线采样点;n表示管径曲线采样数量的总和;x i表示第i个管径的曲线采样点的长度;y i表示在x i处的管径; Among them, i represents the curve sampling point of the ith pipe diameter; n represents the sum of the number of pipe diameter curve sampling; xi represents the length of the curve sampling point of the ith pipe diameter; y i represents the pipe diameter at xi ;
将每个所述拟合管径对应到坐标系中获得相应的管径点,将所述管径点依次平滑连接,获取所述拟合管径曲线。Corresponding each of the fitted pipe diameters to a coordinate system to obtain a corresponding pipe diameter point, and smoothly connecting the pipe diameter points in turn to obtain the fitted pipe diameter curve.
可选地是,上述的精确获取狭窄病变区间的方法,所述根据所述拟合管径曲线和真实管径曲线,获取第一狭窄位置的方法包括:Optionally, in the above-mentioned method for accurately obtaining a stenotic lesion interval, the method for obtaining the first stenosis position according to the fitted diameter curve and the real diameter curve includes:
获取所述血管的真实管径;Obtain the true diameter of the blood vessel;
将所述真实管径对应到拟合管径曲线的坐标系中;Corresponding the real pipe diameter to the coordinate system of the fitted pipe diameter curve;
获取所述真实管径曲线,以及所述真实管径曲线与所述拟合管径曲线的交叉点;Obtain the true pipe diameter curve, and the intersection of the true pipe diameter curve and the fitted pipe diameter curve;
如果交叉点前一点的所述真实管径大于所述拟合管径,则所述交叉点为所述狭窄区域的第一入口点,反之,所述交叉点为所述狭窄区域的第一出口点;If the real pipe diameter at a point before the intersection is larger than the fitted pipe diameter, the intersection is the first entry point of the narrow area, otherwise, the intersection is the first exit of the narrow area point;
所述第一入口点和所述第一出口点之间的曲线为初步判定的狭窄位置,即所述第一狭窄病变区间。The curve between the first entry point and the first exit point is the preliminarily determined stenosis position, that is, the first stenotic lesion interval.
可选地是,上述的精确获取狭窄病变区间的方法,所述从所述第一狭窄病变区间去除误判的狭窄区域,获得第二狭窄病变区间的方法包括;Optionally, in the above-mentioned method for accurately obtaining a narrow lesion interval, the method for removing a misjudged narrow lesion area from the first narrow lesion interval to obtain a second narrow lesion interval includes;
计算狭窄度;Calculate stenosis;
计算所述第一狭窄病变区间的血管中心线长度L;calculating the length L of the blood vessel centerline of the first stenotic lesion interval;
根据所述狭窄度和血管中心线长度从所述第一狭窄病变区间去除误判的狭窄区域,获得第二狭窄病变区间。The misjudged stenosis region is removed from the first stenotic lesion interval according to the stenosis degree and the length of the blood vessel center line to obtain a second stenotic lesion interval.
可选地是,上述的精确获取狭窄病变区间的方法,所述计算狭窄度的方法包括:Optionally, in the above-mentioned method for accurately obtaining a stenotic lesion interval, the method for calculating a stenosis degree includes:
Figure PCTCN2020131703-appb-000002
Figure PCTCN2020131703-appb-000002
其中,A表示血管的狭窄度,D min表示第一入口点与第一出口点之间的血管的最小管径,D 、D 分别表示第一入口点的血管管径和第一出口点的血管管径。 Among them, A represents the stenosis degree of the blood vessel, D min represents the minimum diameter of the blood vessel between the first entry point and the first exit point, D in and D out represent the vessel diameter of the first entry point and the first exit point, respectively vascular diameter.
可选地是,上述的精确获取狭窄病变区间的方法,所述根据所述狭窄度和血管中心线长度从所述第一狭窄病变区间去除误判的狭窄区域,获得第二狭窄病变区间的方法包括:Optionally, in the above-mentioned method for accurately obtaining a stenotic lesion interval, the method for obtaining a second stenotic lesion interval by removing a misjudged stenotic area from the first stenotic lesion interval according to the stenosis degree and the length of the blood vessel centerline include:
如果A<0.2,则判断为误判区域,采用本区域的拟合管径曲线替代所述误判区域内的真实管径曲线;If A<0.2, it is judged as a misjudged area, and the fitted pipe diameter curve in this area is used to replace the real pipe diameter curve in the misjudged area;
如果L<5mm,则判断为误判区域,采用本区域的拟合管径曲线替代所述误判区域内的真实管径曲线;If L<5mm, it is judged as a misjudged area, and the fitted pipe diameter curve in this area is used to replace the real pipe diameter curve in the misjudged area;
去除误判区域后重新获得的狭窄区域,即为所述第二狭窄病变区间。The stenotic area obtained after removing the misjudged area is the second stenotic lesion interval.
可选地是,上述的精确获取狭窄病变区间的方法,还包括:Optionally, the above-mentioned method for accurately obtaining a narrow lesion interval further includes:
取所述第一入口点前1~3cm,所述第一出口点后1~3cm之间的区域,根据所述拟合代价函数重新拟合管径区间曲线;Take the area between 1-3 cm before the first entry point and 1-3 cm after the first exit point, and re-fit the pipe diameter interval curve according to the fitting cost function;
获取所述重新拟合管径区间曲线中管径最小的点为狭窄点;Obtain the point with the smallest pipe diameter in the refitted pipe diameter interval curve as the narrow point;
在所述狭窄点两侧,所述拟合管径区间曲线与所述真实管径曲线交叉的两个点为第二入口点和第二出口点,所述第二入口点和所述第二出口点之间的区间为第三狭窄病变区间。On both sides of the narrow point, the two points where the fitted pipe diameter interval curve and the real pipe diameter curve intersect are a second entry point and a second exit point, and the second entry point and the second The interval between the exit points is the third stenotic lesion interval.
可选地是,上述的精确获取狭窄病变区间的方法,所述根据血管实时直径D t、血管中心线长度L和狭窄区间进行三维建模,形成具有狭窄病变区间的数学模型的方法,包括: Optionally, the above-mentioned method for accurately obtaining the stenotic lesion interval, the method for performing three-dimensional modeling according to the real-time diameter D t of the blood vessel, the blood vessel centerline length L and the stenotic interval to form a mathematical model with the stenotic lesion interval, includes:
根据血管实时直径D t、血管中心线长度L和狭窄区间进行三维建模,形成具有狭窄病变区间的三维血管模型; 3D modeling is carried out according to the real-time diameter D t of the blood vessel, the length L of the blood vessel center line and the stenosis interval to form a 3D vessel model with stenotic lesion interval;
沿着具有狭窄病变区间的所述三维血管模型的圆周面进行N边型网格划分,形成单层网格模型,其中N≥6;N-edge meshing is performed along the circumference of the three-dimensional blood vessel model with stenotic lesions to form a single-layer mesh model, where N≥6;
对所述单层网格模型进行表面分层化处理,形成双层网格模型,即血管数学模型。The surface layering process is performed on the single-layer grid model to form a double-layer grid model, that is, a blood vessel mathematical model.
可选地是,上述的精确获取狭窄病变区间的方法,所述沿着具有狭窄病变区间的所述三维血管模型的圆周面进行N边型网格划分,形成单层网格模型,其中N≥6的方法包括:Optionally, in the above-mentioned method for accurately obtaining stenotic lesion intervals, N-edge meshing is performed along the circumference of the three-dimensional blood vessel model with stenotic lesion intervals to form a single-layer grid model, where N≥ 6 methods include:
沿着所述具有狭窄病变区间的所述三维血管模型的圆周面,以三角形为最小单元进行网格划分;along the circumference of the three-dimensional blood vessel model with the stenotic lesion interval, meshing is performed with triangles as the smallest unit;
按照顺序,每N个三角形组合转换成1个N边形,形成N边形初始网格;According to the sequence, every N triangle combination is converted into an N-sided shape to form an N-sided initial mesh;
删除所述N边形初始网格中每个N边形内部的连接线,形成单层N边形网格模型,其中N≥6。The connecting lines inside each N-gon in the N-sided initial grid are deleted to form a single-layer N-sided grid model, where N≥6.
可选地是,上述的精确获取狭窄病变区间的方法,所述沿着所述具有狭窄病变区间的所述三维血管模型的圆周面,以三角形为最小单元进行网格划分的方法包括:Optionally, in the above-mentioned method for accurately acquiring stenotic lesion intervals, the method for meshing along the circumference of the three-dimensional blood vessel model with stenotic lesion intervals using triangles as the smallest unit includes:
将具有狭窄病变区间的所述三维血管模型分割成K段,segmenting the three-dimensional vessel model with stenotic lesion intervals into K segments,
在每分段所述三维血管模型的圆周面上,以三角形为最小单元进行网格划分。On the circumferential surface of the three-dimensional blood vessel model in each segment, a triangle is used as the smallest unit for mesh division.
可选地是,上述的精确获取狭窄病变区间的方法,所述对所述单层网格模型进行表面分层化处理,形成双层网格模型,即血管数学模型的方法包括:Optionally, in the above-mentioned method for accurately obtaining a stenotic lesion interval, the method for performing surface layering processing on the single-layer grid model to form a double-layer grid model, that is, a blood vessel mathematical model, includes:
获取血管壁厚h;Obtain the vessel wall thickness h;
根据所述血管壁厚h、血管起始直径D 、血管结束直径D 和血管中心线长度L进行三维建模,在所述单层网格模型内表面或者外表面形成圆台三维模型; Three-dimensional modeling is performed according to the blood vessel wall thickness h, the blood vessel starting diameter D, the blood vessel ending diameter D and the blood vessel centerline length L, and a three-dimensional truncated truncated model is formed on the inner surface or outer surface of the single-layer mesh model;
根据所述单层网格模型的获取方法,沿着所述圆台三维模型的圆周面进行N边型网格划分,形成另一单层网格模型;According to the method for obtaining the single-layer grid model, N-edge grid division is performed along the circumference of the three-dimensional truncated truncated model to form another single-layer grid model;
两层所述单层网格模型与所述血管壁厚h,形成所述双层网格模型,即所述血管数学模型。Two layers of the single-layer grid model and the blood vessel wall thickness h form the double-layer grid model, that is, the blood vessel mathematical model.
第二方面,本申请提供了一种计算机存储介质,计算机程序被处理器执行时实现上述的精确获取狭窄病变区间的方法。In a second aspect, the present application provides a computer storage medium, and when the computer program is executed by a processor, the above-mentioned method for accurately acquiring a narrow lesion interval can be implemented.
本申请实施例提供的方案带来的有益效果至少包括:The beneficial effects brought by the solutions provided in the embodiments of the present application include at least:
本申请提供了提供一种精确获取狭窄病变区间的方法及存储介质,通过对血管壁重新拟合后的三维建模,对粗精度的血管狭窄进一步修正,保证狭窄病 变区间精度。The present application provides a method and a storage medium for accurately obtaining a stenosis lesion interval, and further corrects the vascular stenosis with coarse precision by refitting the three-dimensional modeling of the vessel wall to ensure the accuracy of the stenosis lesion interval.
附图说明Description of drawings
此处所说明的附图用来提供对本发明的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:The accompanying drawings described herein are used to provide further understanding of the present invention and constitute a part of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the attached image:
图1为本申请的精确获取狭窄病变区间的方法的一个实施例的流程图;FIG. 1 is a flowchart of an embodiment of a method for accurately obtaining a narrow lesion interval according to the present application;
图2为本申请的S100的流程图;Fig. 2 is the flow chart of S100 of this application;
图3为本申请的S120的流程图;3 is a flowchart of S120 of the application;
图4为本申请的S130的流程图;4 is a flowchart of S130 of the application;
图5为本申请的S200的流程图;Fig. 5 is the flow chart of S200 of this application;
图6为本申请的S210的流程图;6 is a flowchart of S210 of the application;
图7为本申请的S230的流程图;FIG. 7 is a flowchart of S230 of the application;
图8为本申请的S300的流程图;8 is a flowchart of S300 of the application;
图9为本申请的S330的流程图;9 is a flowchart of S330 of the application;
图10为本申请的精确获取狭窄病变区间的方法的另一实施例的流程图。FIG. 10 is a flowchart of another embodiment of the method for accurately obtaining a stenotic lesion interval according to the present application.
具体实施方式Detailed ways
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to make the objectives, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the specific embodiments of the present invention and the corresponding drawings. Obviously, the described embodiments are only some, but not all, embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
以下将以图式揭露本发明的多个实施方式,为明确说明起见,许多实务上的细节将在以下叙述中一并说明。然而,应了解到,这些实务上的细节不应用以限制本发明。也就是说,在本发明的部分实施方式中,这些实务上的细节是非必要的。此外,为简化图式起见,一些习知惯用的结构与组件在图式中将以 简单的示意的方式绘示之。Various embodiments of the present invention will be disclosed in the drawings below, and for the sake of clarity, many practical details will be described together in the following description. It should be understood, however, that these practical details should not be used to limit the invention. That is, in some embodiments of the invention, these practical details are unnecessary. In addition, for the purpose of simplifying the drawings, some well-known structures and components will be shown in a simple schematic manner in the drawings.
实施例1:Example 1:
如图1所示,本申请提供了一种精确获取狭窄病变区间的方法,包括:As shown in FIG. 1 , the present application provides a method for accurately obtaining a narrow lesion interval, including:
S100,如图2所示,根据血管实时直径D t、血管中心线长度L和狭窄区间进行三维建模,形成具有狭窄病变区间的数学模型,包括: S100 , as shown in FIG. 2 , three-dimensional modeling is performed according to the real-time diameter D t of the blood vessel, the length L of the blood vessel centerline and the stenosis interval, to form a mathematical model with a stenotic lesion interval, including:
S110,根据血管实时直径D t、血管中心线长度L和狭窄区间进行三维建模,形成具有狭窄病变区间的三维血管模型; S110, perform three-dimensional modeling according to the real-time diameter D t of the blood vessel, the length L of the blood vessel center line and the stenosis interval, to form a three-dimensional blood vessel model with a stenotic lesion interval;
S120,如图3所示,沿着具有狭窄病变区间的三维血管模型的圆周面进行N边型网格划分,形成单层网格模型,其中N≥6,包括:S120 , as shown in FIG. 3 , N-edge meshing is performed along the circumference of the three-dimensional blood vessel model with stenotic lesions to form a single-layer mesh model, where N≥6, including:
S121,沿着具有狭窄病变区间的三维血管模型的圆周面,以三角形为最小单元进行网格划分,包括:将具有狭窄病变区间的三维血管模型分割成K段,在每分段三维血管模型的圆周面上,以三角形为最小单元进行网格划分。S121 , along the circumference of the three-dimensional blood vessel model with the stenotic lesion interval, perform grid division with triangles as the smallest unit, including: dividing the three-dimensional blood vessel model with the narrow lesion interval into K segments, and in each segment of the three-dimensional blood vessel model On the circumferential surface, the triangle is used as the smallest element for mesh division.
S122,按照顺序,每N个三角形组合转换成1个N边形,形成N边形初始网格;S122, according to the sequence, every N triangle combination is converted into an N-sided shape to form an N-sided initial mesh;
S123,删除N边形初始网格中每个N边形内部的连接线,形成单层N边形网格模型,其中N≥6。S123 , delete the connecting lines inside each N-gon in the N-sided initial grid to form a single-layer N-sided grid model, where N≥6.
S130,如图4所示,对单层网格模型进行表面分层化处理,形成双层网格模型,即血管数学模型,包括:S130, as shown in FIG. 4, the surface layering process is performed on the single-layer grid model to form a double-layer grid model, that is, a blood vessel mathematical model, including:
S131,获取血管壁厚h;S131, obtaining the blood vessel wall thickness h;
S132,根据血管壁厚h、血管起始直径D 、血管结束直径D 和血管中心线长度L进行三维建模,在单层网格模型内表面或者外表面形成圆台三维模型; S132, perform three-dimensional modeling according to the blood vessel wall thickness h, the blood vessel starting diameter D, the blood vessel ending diameter D and the blood vessel centerline length L, and form a three-dimensional truncated truncated model on the inner surface or outer surface of the single-layer mesh model;
S133,根据单层网格模型的获取方法,沿着圆台三维模型的圆周面进行N边型网格划分,形成另一单层网格模型;S133, according to the obtaining method of the single-layer mesh model, carry out N-edge mesh division along the circumferential surface of the three-dimensional truncated truncated model to form another single-layer mesh model;
S134,两层单层网格模型与血管壁厚h,形成双层网格模型,即血管数学 模型。S134, the two-layer single-layer grid model and the blood vessel wall thickness h form a double-layer grid model, that is, a blood vessel mathematical model.
S200,如图5所示,初步判断,获得第一狭窄病变区间,包括:S200, as shown in Fig. 5, preliminarily judges and obtains the first stenotic lesion interval, including:
S210,如图6所示,拟合正常血管管径,获取拟合管径曲线,包括:S210, as shown in FIG. 6, fitting a normal blood vessel diameter, and obtaining a fitting diameter curve, including:
S211,根据拟合代价函数,获取拟合管径,具体公式为:S211, according to the fitting cost function, obtain the fitting pipe diameter, and the specific formula is:
Figure PCTCN2020131703-appb-000003
Figure PCTCN2020131703-appb-000003
其中,i表示第i个管径的曲线采样点;n表示管径曲线采样数量的总和;x i表示第i个管径的曲线采样点的长度;y i表示在x i处的管径; Among them, i represents the curve sampling point of the ith pipe diameter; n represents the sum of the number of pipe diameter curve sampling; xi represents the length of the curve sampling point of the ith pipe diameter; y i represents the pipe diameter at xi ;
S212,将每个拟合管径对应到坐标系中获得相应的管径点,将管径点依次平滑连接,获取拟合管径曲线。S212: Corresponding each fitted pipe diameter to a coordinate system to obtain a corresponding pipe diameter point, and smoothly connecting the pipe diameter points in turn to obtain a fitted pipe diameter curve.
S220,根据数学模型获取真实管径曲线;S220, obtaining the real pipe diameter curve according to the mathematical model;
S230,如图7所示,根据拟合管径曲线和真实管径曲线,获取第一狭窄病变区间,包括:S230, as shown in FIG. 7, according to the fitted diameter curve and the real diameter curve, obtain a first stenosis lesion interval, including:
S231,获取血管的真实管径;S231, obtain the real diameter of the blood vessel;
S232,将真实管径对应到拟合管径曲线的坐标系中;S232, corresponding the real pipe diameter to the coordinate system of the fitted pipe diameter curve;
S233,获取真实管径曲线,以及真实管径曲线与拟合管径曲线的交叉点;S233, obtaining the real pipe diameter curve and the intersection point of the real pipe diameter curve and the fitted pipe diameter curve;
S234,如果交叉点前一点的真实管径大于拟合管径,则交叉点为狭窄区域的第一入口点,反之,交叉点为狭窄区域的第一出口点;S234, if the real pipe diameter at the point before the intersection is larger than the fitted pipe diameter, the intersection is the first entry point of the narrow area, otherwise, the intersection is the first exit point of the narrow area;
S235,第一入口点和第一出口点之间的曲线为初步判定的狭窄位置,即第一狭窄病变区间。S235, the curve between the first entry point and the first exit point is the initially determined stenosis position, that is, the first stenotic lesion interval.
S300,如图8所示,从第一狭窄病变区间去除误判的狭窄区域,获得第二狭窄病变区间,包括:S300 , as shown in FIG. 8 , the misjudged stenosis area is removed from the first stenotic lesion interval to obtain a second stenotic lesion interval, including:
S310,计算狭窄度,具体公式为:
Figure PCTCN2020131703-appb-000004
其中,A表示血管的狭窄度,D min表示第一入口点与第一出口点之间的血管的最小管径,D 、D 分别表示第一入口点的血管管径和第一出口点的血管管径。
S310, calculate the stenosis, and the specific formula is:
Figure PCTCN2020131703-appb-000004
Among them, A represents the stenosis degree of the blood vessel, D min represents the minimum diameter of the blood vessel between the first entry point and the first exit point, D in and D out represent the vessel diameter of the first entry point and the first exit point, respectively vascular diameter.
S320,计算第一狭窄病变区间的血管中心线长度L;S320, calculating the length L of the center line of the blood vessel in the first stenotic lesion interval;
S330,如图9所示,根据狭窄度和血管中心线长度从第一狭窄病变区间去除误判的狭窄区域,获得第二狭窄病变区间,包括:S330 , as shown in FIG. 9 , according to the stenosis degree and the length of the blood vessel centerline, the misjudged stenosis area is removed from the first stenotic lesion interval, and the second stenotic lesion interval is obtained, including:
S331,如果A<0.2,则判断为误判区域,采用本区域的拟合管径曲线替代误判区域内的真实管径曲线;S331, if A<0.2, it is judged as a misjudged area, and the fitted pipe diameter curve in this area is used to replace the real pipe diameter curve in the misjudged area;
S332,如果L<5mm,则判断为误判区域,采用本区域的拟合管径曲线替代误判区域内的真实管径曲线;S332, if L<5mm, it is judged as a misjudged area, and the fitted pipe diameter curve in this area is used to replace the real pipe diameter curve in the misjudged area;
S333,去除误判区域后重新获得的狭窄区域,即为第二狭窄病变区间。S333 , the stenotic area obtained again after removing the misjudged area is the second stenotic lesion interval.
实施例2:Example 2:
如图10所示,在实施例1的基础上,本申请提供的一种精确获取狭窄病变区间的方法,还包括:As shown in FIG. 10 , on the basis of Embodiment 1, a method for accurately obtaining a narrow lesion interval provided by the present application further includes:
S400,取第一入口点前1~3cm,第一出口点后1~3cm之间的区域,根据拟合代价函数重新拟合管径区间曲线;S400, take the area between 1-3 cm before the first entry point and 1-3 cm after the first exit point, and re-fit the pipe diameter interval curve according to the fitting cost function;
S500,获取重新拟合管径区间曲线中管径最小的点为狭窄点;S500, obtaining the point with the smallest pipe diameter in the re-fitted pipe diameter interval curve as the narrow point;
S600,在狭窄点两侧,拟合管径区间曲线与真实管径曲线交叉的两个点为第二入口点和第二出口点,第二入口点和第二出口点之间的区间为第三狭窄病变区间。S600, on both sides of the narrow point, the two points where the fitted pipe diameter interval curve and the real pipe diameter curve intersect are the second entry point and the second exit point, and the interval between the second entry point and the second exit point is the second entry point and the second exit point. Three stenotic lesions.
本申请提供了一种计算机存储介质,计算机程序被处理器执行时实现上述的精确获取狭窄病变区间的方法。The present application provides a computer storage medium, and when the computer program is executed by a processor, the above-mentioned method for accurately acquiring a narrow lesion interval can be implemented.
所属技术领域的技术人员知道,本发明的各个方面可以实现为系统、方法或计算机程序产品。因此,本发明的各个方面可以具体实现为以下形式,即:完全的硬件实施方式、完全的软件实施方式(包括固件、驻留软件、微代码等),或硬件和软件方面结合的实施方式,这里可以统称为“电路”、“模块”或“系统”。此外,在一些实施例中,本发明的各个方面还可以实现为在一个或多个 计算机可读介质中的计算机程序产品的形式,该计算机可读介质中包含计算机可读的程序代码。本发明的实施例的方法和/或系统的实施方式可以涉及到手动地、自动地或以其组合的方式执行或完成所选任务。As will be appreciated by one skilled in the art, various aspects of the present invention may be embodied as a system, method or computer program product. Accordingly, various aspects of the present invention may be embodied in the form of an entirely hardware implementation, an entirely software implementation (including firmware, resident software, microcode, etc.), or a combination of hardware and software aspects, It may be collectively referred to herein as a "circuit," "module," or "system." Furthermore, in some embodiments, various aspects of the present invention may also be implemented in the form of a computer program product on one or more computer-readable media having computer-readable program code embodied thereon. Implementation of the method and/or system of embodiments of the invention may involve performing or completing selected tasks manually, automatically, or a combination thereof.
例如,可以将用于执行根据本发明的实施例的所选任务的硬件实现为芯片或电路。作为软件,可以将根据本发明的实施例的所选任务实现为由计算机使用任何适当操作系统执行的多个软件指令。在本发明的示例性实施例中,由数据处理器来执行如本文的根据方法和/或系统的示例性实施例的一个或多个任务,诸如用于执行多个指令的计算平台。可选地,该数据处理器包括用于存储指令和/或数据的易失性储存器和/或用于存储指令和/或数据的非易失性储存器,例如,磁硬盘和/或可移动介质。可选地,也提供了一种网络连接。可选地也提供显示器和/或用户输入设备,诸如键盘或鼠标。For example, hardware for performing selected tasks according to embodiments of the invention may be implemented as a chip or a circuit. As software, selected tasks according to embodiments of the invention may be implemented as a plurality of software instructions executed by a computer using any suitable operating system. In an exemplary embodiment of the invention, one or more tasks in accordance with exemplary embodiments of methods and/or systems as herein are performed by a data processor, such as a computing platform for executing a plurality of instructions. Optionally, the data processor includes volatile storage for storing instructions and/or data and/or non-volatile storage for storing instructions and/or data, such as a magnetic hard disk and/or a Move media. Optionally, a network connection is also provided. A display and/or user input device, such as a keyboard or mouse, is optionally also provided.
可利用一个或多个计算机可读的任何组合。计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举列表)将包括以下各项:Any combination of one or more computer readable may be utilized. The computer-readable medium may be a computer-readable signal medium or a computer-readable storage medium. The computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus or device, or a combination of any of the above. More specific examples (non-exhaustive list) of computer-readable storage media would include the following:
具有一个或多个导线的电连接、便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。Electrical connection with one or more wires, portable computer disk, hard disk, random access memory (RAM), read only memory (ROM), erasable programmable read only memory (EPROM or flash memory), optical fiber, portable compact disk Read only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination of the above. In this document, a computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多 种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。A computer-readable signal medium may include a propagated data signal in baseband or as part of a carrier wave, with computer-readable program code embodied thereon. Such propagated data signals may take a variety of forms including, but not limited to, electromagnetic signals, optical signals, or any suitable combination of the foregoing. A computer-readable signal medium can also be any computer-readable medium other than a computer-readable storage medium that can transmit, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device .
计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括(但不限于)无线、有线、光缆、RF等等,或者上述的任意合适的组合。Program code embodied on a computer-readable medium may be transmitted using any suitable medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
例如,可用一个或多个编程语言的任何组合来编写用于执行用于本发明的各方面的操作的计算机程序代码,包括诸如Java、Smalltalk、C++等面向对象编程语言和常规过程编程语言,诸如"C"编程语言或类似编程语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络--包括局域网(LAN)或广域网(WAN)-连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。For example, computer program code for performing operations for various aspects of the present invention may be written in any combination of one or more programming languages, including object-oriented programming languages such as Java, Smalltalk, C++, and conventional procedural programming languages, such as The "C" programming language or similar programming language. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer, or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any kind of network - including a local area network (LAN) or wide area network (WAN) - or may be connected to an external computer (eg using an Internet service provider via Internet connection).
应当理解,流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合,都可以由计算机程序指令实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器,从而生产出一种机器,使得这些计算机程序指令在通过计算机或其它可编程数据处理装置的处理器执行时,产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to the processor of a general purpose computer, special purpose computer or other programmable data processing apparatus to produce a machine that causes the computer program instructions when executed by the processor of the computer or other programmable data processing apparatus , resulting in means for implementing the functions/acts specified in one or more blocks of the flowchart and/or block diagrams.
也可以把这些计算机程序指令存储在计算机可读介质中,这些指令使得计算机、其它可编程数据处理装置、或其它设备以特定方式工作,从而,存储在计算机可读介质中的指令就产生出包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的指令的制造品(article of manufacture)。These computer program instructions can also be stored on a computer readable medium, the instructions cause a computer, other programmable data processing apparatus, or other device to operate in a particular manner, whereby the instructions stored on the computer readable medium produce the An article of manufacture of instructions implementing the functions/acts specified in one or more blocks of the flowcharts and/or block diagrams.
还可将计算机程序指令加载到计算机(例如,冠状动脉分析系统)或其它可编程数据处理设备上以促使在计算机、其它可编程数据处理设备或其它设备上执行一系列操作步骤以产生计算机实现过程,使得在计算机、其它可编程装置或其它设备上执行的指令提供用于实现在流程图和/或一个或多个框图方框中指定的功能/动作的过程。Computer program instructions can also be loaded on a computer (eg, a coronary artery analysis system) or other programmable data processing device to cause a series of operational steps to be performed on the computer, other programmable data processing device or other device to produce a computer-implemented process , such that instructions executing on a computer, other programmable apparatus, or other device provide a process for implementing the functions/acts specified in the flowchart and/or one or more block diagram blocks.
本发明的以上的具体实例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above specific examples of the present invention further describe the purpose, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above are only specific embodiments of the present invention, and are not intended to limit the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims (13)

  1. 一种精确获取狭窄病变区间的方法,其特征在于,包括:A method for accurately obtaining a stenotic lesion interval, comprising:
    根据血管实时直径D t、血管中心线长度L和狭窄区间进行三维建模,形成具有狭窄病变区间的数学模型; According to the real-time diameter D t of the blood vessel, the length L of the center line of the blood vessel and the stenosis interval, three-dimensional modeling is performed to form a mathematical model with stenotic lesion interval;
    初步判断,获得第一狭窄病变区间;Preliminary judgment to obtain the first stenotic lesion interval;
    从所述第一狭窄病变区间去除误判的狭窄区域,获得第二狭窄病变区间。The misjudged stenotic region is removed from the first stenotic lesion section to obtain a second stenotic lesion section.
  2. 根据权利要求1所述的精确获取狭窄病变区间的方法,其特征在于,所述初步判断,获得第一狭窄病变区间的方法,包括:The method for accurately obtaining a narrow lesion interval according to claim 1, wherein the preliminary judgment, the method for obtaining the first narrow lesion interval, comprises:
    拟合正常血管管径,获取拟合管径曲线;Fit the normal blood vessel diameter and obtain the fitted diameter curve;
    根据所述数学模型获取真实管径曲线;Obtain the real pipe diameter curve according to the mathematical model;
    根据所述拟合管径曲线和真实管径曲线,获取第一狭窄病变区间。According to the fitted diameter curve and the real diameter curve, the first stenosis lesion interval is obtained.
  3. 根据权利要求2所述的精确获取狭窄病变区间的方法,其特征在于,所述拟合正常血管管径,获取拟合管径曲线的方法包括:The method for accurately obtaining a stenotic lesion interval according to claim 2, wherein the method for fitting a normal blood vessel diameter and obtaining a fitting diameter curve comprises:
    根据拟合代价函数,获取拟合管径,具体公式为
    Figure PCTCN2020131703-appb-100001
    According to the fitting cost function, the fitted pipe diameter is obtained, and the specific formula is:
    Figure PCTCN2020131703-appb-100001
    其中,i表示第i个管径的曲线采样点;n表示管径曲线采样数量的总和;x i表示第i个管径的曲线采样点的长度;y i表示在x i处的管径; Among them, i represents the curve sampling point of the ith pipe diameter; n represents the sum of the number of pipe diameter curve sampling; xi represents the length of the curve sampling point of the ith pipe diameter; y i represents the pipe diameter at xi ;
    将每个所述拟合管径对应到坐标系中获得相应的管径点,将所述管径点依次平滑连接,获取所述拟合管径曲线。Corresponding each of the fitted pipe diameters to a coordinate system to obtain a corresponding pipe diameter point, and smoothly connecting the pipe diameter points in turn to obtain the fitted pipe diameter curve.
  4. 根据权利要求3所述的精确获取狭窄病变区间的方法,其特征在于,所述根据所述拟合管径曲线和真实管径曲线,获取第一狭窄位置的方法包括:The method for accurately obtaining a stenotic lesion interval according to claim 3, wherein the method for obtaining the first stenosis position according to the fitted diameter curve and the real diameter curve comprises:
    获取所述血管的真实管径;Obtain the true diameter of the blood vessel;
    将所述真实管径对应到拟合管径曲线的坐标系中;Corresponding the real pipe diameter to the coordinate system of the fitted pipe diameter curve;
    获取所述真实管径曲线,以及所述真实管径曲线与所述拟合管径曲线的交 叉点;Obtain the true pipe diameter curve, and the intersection of the true pipe diameter curve and the fitted pipe diameter curve;
    如果交叉点前一点的所述真实管径大于所述拟合管径,则所述交叉点为所述狭窄区域的第一入口点,反之,所述交叉点为所述狭窄区域的第一出口点;If the real pipe diameter at a point before the intersection is larger than the fitted pipe diameter, the intersection is the first entry point of the narrow area, otherwise, the intersection is the first exit of the narrow area point;
    所述第一入口点和所述第一出口点之间的曲线为初步判定的狭窄位置,即所述第一狭窄病变区间。The curve between the first entry point and the first exit point is the preliminarily determined stenosis position, that is, the first stenotic lesion interval.
  5. 根据权利要求4所述的精确获取狭窄病变区间的方法,其特征在于,所述从所述第一狭窄病变区间去除误判的狭窄区域,获得第二狭窄病变区间的方法包括;The method for accurately obtaining a stenotic lesion interval according to claim 4, wherein the method for obtaining a second stenotic lesion interval by removing a misjudged stenotic area from the first stenotic lesion interval comprises;
    计算狭窄度;Calculate stenosis;
    计算所述第一狭窄病变区间的血管中心线长度L;calculating the length L of the blood vessel centerline of the first stenotic lesion interval;
    根据所述狭窄度和血管中心线长度从所述第一狭窄病变区间去除误判的狭窄区域,获得第二狭窄病变区间。The misjudged stenosis region is removed from the first stenotic lesion interval according to the stenosis degree and the length of the blood vessel center line to obtain a second stenotic lesion interval.
  6. 根据权利要求5所述的精确获取狭窄病变区间的方法,其特征在于,所述计算狭窄度的方法包括:The method for accurately obtaining a stenotic lesion interval according to claim 5, wherein the method for calculating the stenosis degree comprises:
    Figure PCTCN2020131703-appb-100002
    Figure PCTCN2020131703-appb-100002
    其中,A表示血管的狭窄度,D min表示第一入口点与第一出口点之间的血管的最小管径,D 、D 分别表示第一入口点的血管管径和第一出口点的血管管径。 Among them, A represents the stenosis degree of the blood vessel, D min represents the minimum diameter of the blood vessel between the first entry point and the first exit point, D in and D out represent the vessel diameter of the first entry point and the first exit point, respectively vascular diameter.
  7. 根据权利要求6所述的精确获取狭窄病变区间的方法,其特征在于,所述根据所述狭窄度和血管中心线长度从所述第一狭窄病变区间去除误判的狭窄区域,获得第二狭窄病变区间的方法包括:The method for accurately obtaining a stenosis lesion section according to claim 6, wherein the stenosis region that is misjudged from the first stenosis lesion section is removed according to the stenosis degree and the length of the blood vessel centerline to obtain the second stenosis Methods of lesion interval include:
    如果A<0.2,则判断为误判区域,采用本区域的拟合管径曲线替代所述误判区域内的真实管径曲线;If A<0.2, it is judged as a misjudged area, and the fitted pipe diameter curve in this area is used to replace the real pipe diameter curve in the misjudged area;
    如果L<5mm,则判断为误判区域,采用本区域的拟合管径曲线替代所述 误判区域内的真实管径曲线;If L<5mm, it is judged as a misjudged area, and the fitted pipe diameter curve in this area is used to replace the real pipe diameter curve in the misjudged area;
    去除误判区域后重新获得的狭窄区域,即为所述第二狭窄病变区间。The stenotic area obtained after removing the misjudged area is the second stenotic lesion interval.
  8. 根据权利要求4所述的精确获取狭窄病变区间的方法,其特征在于,还包括:The method for accurately obtaining a narrow lesion interval according to claim 4, further comprising:
    取所述第一入口点前1~3cm,所述第一出口点后1~3cm之间的区域,根据所述拟合代价函数重新拟合管径区间曲线;Take the area between 1-3 cm before the first entry point and 1-3 cm after the first exit point, and re-fit the pipe diameter interval curve according to the fitting cost function;
    获取所述重新拟合管径区间曲线中管径最小的点为狭窄点;Obtain the point with the smallest pipe diameter in the refitted pipe diameter interval curve as the narrow point;
    在所述狭窄点两侧,所述拟合管径区间曲线与所述真实管径曲线交叉的两个点为第二入口点和第二出口点,所述第二入口点和所述第二出口点之间的区间为第三狭窄病变区间。On both sides of the narrow point, the two points where the fitted pipe diameter interval curve and the real pipe diameter curve intersect are a second entry point and a second exit point, and the second entry point and the second The interval between the exit points is the third stenotic lesion interval.
  9. 根据权利要求1所述的精确获取狭窄病变区间的方法,其特征在于,所述根据血管实时直径D t、血管中心线长度L和狭窄区间进行三维建模,形成具有狭窄病变区间的数学模型的方法,包括: The method for accurately obtaining a stenotic lesion interval according to claim 1, wherein the three-dimensional modeling is performed according to the real-time blood vessel diameter D t , the blood vessel centerline length L and the stenotic interval to form a mathematical model with a stenotic lesion interval. methods, including:
    根据血管实时直径D t、血管中心线长度L和狭窄区间进行三维建模,形成具有狭窄病变区间的三维血管模型; 3D modeling is carried out according to the real-time diameter D t of the blood vessel, the length L of the blood vessel center line and the stenosis interval to form a 3D vessel model with stenotic lesion interval;
    沿着具有狭窄病变区间的所述三维血管模型的圆周面进行N边型网格划分,形成单层网格模型,其中N≥6;N-edge meshing is performed along the circumference of the three-dimensional blood vessel model with stenotic lesions to form a single-layer mesh model, where N≥6;
    对所述单层网格模型进行表面分层化处理,形成双层网格模型,即血管数学模型。The surface layering process is performed on the single-layer grid model to form a double-layer grid model, that is, a blood vessel mathematical model.
  10. 根据权利要求9所述的精确获取狭窄病变区间的方法,其特征在于,所述沿着具有狭窄病变区间的所述三维血管模型的圆周面进行N边型网格划分,形成单层网格模型,其中N≥6的方法包括:The method for accurately obtaining a stenotic lesion interval according to claim 9, wherein the N-edge meshing is performed along the circumference of the three-dimensional blood vessel model with the stenotic lesion interval to form a single-layer grid model , where N≥6 methods include:
    沿着所述具有狭窄病变区间的所述三维血管模型的圆周面,以三角形为最小单元进行网格划分;along the circumference of the three-dimensional blood vessel model with the stenotic lesion interval, meshing is performed with triangles as the smallest unit;
    按照顺序,每N个三角形组合转换成1个N边形,形成N边形初始网格;According to the sequence, every N triangle combination is converted into an N-sided shape to form an N-sided initial mesh;
    删除所述N边形初始网格中每个N边形内部的连接线,形成单层N边形网格模型,其中N≥6。The connecting lines inside each N-gon in the N-sided initial grid are deleted to form a single-layer N-sided grid model, where N≥6.
  11. 根据权利要求10所述的精确获取狭窄病变区间的方法,其特征在于,所述沿着所述具有狭窄病变区间的所述三维血管模型的圆周面,以三角形为最小单元进行网格划分的方法包括:The method for accurately obtaining a stenotic lesion section according to claim 10, wherein the method for meshing along the circumference of the three-dimensional blood vessel model with the stenotic lesion section takes triangles as the smallest unit include:
    将具有狭窄病变区间的所述三维血管模型分割成K段,segmenting the three-dimensional vessel model with stenotic lesion intervals into K segments,
    在每分段所述三维血管模型的圆周面上,以三角形为最小单元进行网格划分。On the circumferential surface of the three-dimensional blood vessel model in each segment, a triangle is used as the smallest unit for mesh division.
  12. 根据权利要求11所述的精确获取狭窄病变区间的方法,其特征在于,所述对所述单层网格模型进行表面分层化处理,形成双层网格模型,即血管数学模型的方法包括:The method for accurately obtaining a stenotic lesion interval according to claim 11, wherein the method for performing surface layering processing on the single-layer grid model to form a double-layer grid model, that is, a blood vessel mathematical model, comprises the following steps: :
    获取血管壁厚h;Obtain the vessel wall thickness h;
    根据所述血管壁厚h、血管起始直径D 、血管结束直径D 和血管中心线长度L进行三维建模,在所述单层网格模型内表面或者外表面形成圆台三维模型; Three-dimensional modeling is performed according to the blood vessel wall thickness h, the blood vessel starting diameter D, the blood vessel ending diameter D and the blood vessel centerline length L, and a three-dimensional truncated truncated model is formed on the inner surface or outer surface of the single-layer mesh model;
    根据所述单层网格模型的获取方法,沿着所述圆台三维模型的圆周面进行N边型网格划分,形成另一单层网格模型;According to the method for obtaining the single-layer grid model, N-edge grid division is performed along the circumference of the three-dimensional truncated truncated model to form another single-layer grid model;
    两层所述单层网格模型与所述血管壁厚h,形成所述双层网格模型,即所述血管数学模型。Two layers of the single-layer grid model and the blood vessel wall thickness h form the double-layer grid model, that is, the blood vessel mathematical model.
  13. 一种计算机存储介质,其特征在于,计算机程序被处理器执行时实现权利要求1~12任一项所述的精确获取狭窄病变区间的方法。A computer storage medium, characterized in that, when the computer program is executed by a processor, the method for accurately acquiring a narrow lesion interval according to any one of claims 1 to 12 is implemented.
PCT/CN2020/131703 2020-11-25 2020-11-26 Method and storage medium for precise acquisition of stenotic lesion range WO2022109907A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202011342207.9A CN112419280A (en) 2020-11-25 2020-11-25 Method for accurately acquiring narrow lesion interval and storage medium
CN202011342207.9 2020-11-25

Publications (1)

Publication Number Publication Date
WO2022109907A1 true WO2022109907A1 (en) 2022-06-02

Family

ID=74842430

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/131703 WO2022109907A1 (en) 2020-11-25 2020-11-26 Method and storage medium for precise acquisition of stenotic lesion range

Country Status (2)

Country Link
CN (1) CN112419280A (en)
WO (1) WO2022109907A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115619979A (en) * 2022-11-22 2023-01-17 广州中望龙腾软件股份有限公司 Semi-elongated surface cleaning method, terminal, and storage medium

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150087956A1 (en) * 2012-09-25 2015-03-26 Toshiba Medical Systems Corporation X-ray diagnostic apparatus and medical image processing apparatus
CN108038848A (en) * 2017-12-07 2018-05-15 上海交通大学 Quick calculation method and system based on medical image sequence plaque stability index
CN109872321A (en) * 2019-02-26 2019-06-11 数坤(北京)网络科技有限公司 A kind of hemadostewnosis detection method and equipment
CN110367965A (en) * 2018-09-19 2019-10-25 苏州润迈德医疗科技有限公司 The method, apparatus and system of convenient measurement coronary artery assessment parameters
CN110889896A (en) * 2019-11-11 2020-03-17 苏州润迈德医疗科技有限公司 Method, device and system for obtaining angiostenosis lesion interval and three-dimensional synthesis

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110584639A (en) * 2019-09-04 2019-12-20 北京工业大学 Method for predicting FFR (fringe field response) by processing data of CTA (computed tomography angiography) coronary artery image

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150087956A1 (en) * 2012-09-25 2015-03-26 Toshiba Medical Systems Corporation X-ray diagnostic apparatus and medical image processing apparatus
CN108038848A (en) * 2017-12-07 2018-05-15 上海交通大学 Quick calculation method and system based on medical image sequence plaque stability index
CN110367965A (en) * 2018-09-19 2019-10-25 苏州润迈德医疗科技有限公司 The method, apparatus and system of convenient measurement coronary artery assessment parameters
CN109872321A (en) * 2019-02-26 2019-06-11 数坤(北京)网络科技有限公司 A kind of hemadostewnosis detection method and equipment
CN110889896A (en) * 2019-11-11 2020-03-17 苏州润迈德医疗科技有限公司 Method, device and system for obtaining angiostenosis lesion interval and three-dimensional synthesis

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115619979A (en) * 2022-11-22 2023-01-17 广州中望龙腾软件股份有限公司 Semi-elongated surface cleaning method, terminal, and storage medium
CN115619979B (en) * 2022-11-22 2023-06-02 广州中望龙腾软件股份有限公司 Semi-long surface cleaning method, terminal and storage medium

Also Published As

Publication number Publication date
CN112419280A (en) 2021-02-26

Similar Documents

Publication Publication Date Title
Tu et al. Diagnostic accuracy of quantitative flow ratio for assessment of coronary stenosis significance from a single angiographic view: a novel method based on bifurcation fractal law
JP6283399B2 (en) Method of operating a user interface based medical system for planning stent placement and apparatus for planning stent placement
Pellicano et al. Validation study of image-based fractional flow reserve during coronary angiography
JP6522175B2 (en) System and method for estimating blood flow characteristics from vessel shape and physiology
CN112132949B (en) Three-dimensional modeling method, device and system for blood vessel with narrow lesion region
CN110889896A (en) Method, device and system for obtaining angiostenosis lesion interval and three-dimensional synthesis
CN112419484B (en) Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium
US20220277447A1 (en) Method and apparatus for acquiring contour line of blood vessel according to centerline of blood vessel
CN112164020A (en) Method, device, analysis system and storage medium for accurately extracting blood vessel center line
WO2022109905A1 (en) Method for adjusting blood vessel contours and centreline, and storage medium
WO2022109906A1 (en) Method and system for rendering and synthesizing three-dimensional blood vessel, and storage medium
CN112494020B (en) Method for acquiring aortic pressure curve of interest and storage medium
US20220319116A1 (en) Method and device for synthesizing mathematical model of blood vessel having stenotic lesion
CN113902690A (en) Method, device, computing equipment and storage medium for computing fractional flow reserve based on intravascular images
WO2021097821A1 (en) Method and device for extracting vascular centerline from two-dimensional coronary angiographic image
Van Hamersvelt et al. Diagnostic performance of on-site coronary CT angiography–derived fractional flow reserve based on patient-specific lumped parameter models
WO2022109907A1 (en) Method and storage medium for precise acquisition of stenotic lesion range
WO2022109908A1 (en) Method and system for adjusting vascular stenosis zone, and storage medium
US10354349B2 (en) Systems and methods for using geometry sensitivity information for guiding workflow
WO2021017327A1 (en) Method and apparatus for correcting blood flow velocity on the basis of interval time between angiographic images
WO2022109904A1 (en) Method for two-dimensional image selection and three-dimensional blood vessel synthesis and storage medium
CN112116711B (en) Synthesizing method and device for truncated cone blood vessel mathematical model for hydrodynamic analysis
WO2022109909A1 (en) Method and system for accurately acquiring stenotic lesion interval, and storage medium
Hwang et al. A simple method for automatic 3D reconstruction of coronary arteries from X-ray angiography
CN115222665B (en) Plaque detection method and device, electronic equipment and readable storage medium

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20962798

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20962798

Country of ref document: EP

Kind code of ref document: A1