WO2022060099A1 - 저전압을 이용한 ct 영상에서 혈관 질환 위험도 계산 방법 및 장치 - Google Patents

저전압을 이용한 ct 영상에서 혈관 질환 위험도 계산 방법 및 장치 Download PDF

Info

Publication number
WO2022060099A1
WO2022060099A1 PCT/KR2021/012623 KR2021012623W WO2022060099A1 WO 2022060099 A1 WO2022060099 A1 WO 2022060099A1 KR 2021012623 W KR2021012623 W KR 2021012623W WO 2022060099 A1 WO2022060099 A1 WO 2022060099A1
Authority
WO
WIPO (PCT)
Prior art keywords
value
image
threshold
tube voltage
pixels
Prior art date
Application number
PCT/KR2021/012623
Other languages
English (en)
French (fr)
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 WO2022060099A1 publication Critical patent/WO2022060099A1/ko

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/504Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • 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
    • 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
    • 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
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • 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
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • 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 present invention relates to a method and apparatus for calculating the risk of vascular disease in a CT image using a low voltage. More specifically, it relates to a method for accurately calculating a Coronary Artery Calcium Score (CACS) at a voltage lower than a reference voltage, an apparatus for performing the same, and a CT scanner.
  • CACS Coronary Artery Calcium Score
  • CAC score which can be confirmed by computed tomography (CT) is the most reliable non-invasive predictor of cardiovascular disease and coronary artery disease.
  • the CAC score obtained through CT scan was defined as the Agatston Score calculated as ⁇ (area x concentration factor x section thickness/3 of each selected site). It is determined based on the CT image obtained when tube voltage is applied.
  • the method of calculating the CAC score also tried to calculate the CAC score at a voltage lower than the existing optimized standard (120 kV), but at a lower voltage, the actual CAC It is calculated as an exaggerated score rather than a score, so there is a problem in that people who do not need treatment are treated.
  • the inventors of the present invention found that even with a CT image taken at a tube voltage lower than 120 kV, a threshold value for calculating the CAC score can be obtained in the same way as the CAC score calculated under the reference tube voltage (120 kV tube voltage). An attempt was made to develop a risk calculation method and device.
  • the inventors of the present invention constructed a method and apparatus to obtain a threshold value for calculating a CAC score by comparing/matching a histogram of a reference CT image and a CT image taken at a lower tube voltage.
  • the inventors of the present invention took a CT image using a tube voltage of 100 kV in order to reduce the radiation dose proportional to the size of the tube voltage, and based on this, a method and apparatus for obtaining the patient's coronary calcium score. reached
  • the method includes: acquiring a CT image of a blood vessel using a tube voltage value set by a user; The HU values are counted in ascending order, determining pixels equal to the total number of pixels up to a preset threshold, and setting a HU value corresponding to the determined number of pixels as a new HU threshold. do.
  • the tube voltage value set by the user may have a lower value than the reference tube voltage value set to obtain the reference CT image.
  • the histogram of the reference CT image is divided into at least four sections based on four preset threshold values, and the setting of the new HU threshold value includes the tube voltage set by the user. It may be a step of setting the first to fourth HU threshold values according to
  • the setting of the new HU threshold value may include setting the lowest HU value in the extracted pixel as the first HU threshold value and based on the first HU threshold value.
  • the method may further include sequentially setting the second to fourth HU threshold values.
  • the tube voltage value set by the user may be a tube voltage value of less than 120 kV.
  • the setting of the new HU threshold may further include displaying pixels classified according to the new HU threshold in the histogram of the acquired CT image.
  • the method includes obtaining a CT image by performing computed tomography with a tube voltage value of 100 kV for the coronary artery, and the coronary calcium score based on the first to fourth HU threshold values in the histogram of the obtained CT image. is configured to include the step of calculating
  • the first HU threshold is 141HU or more and 145HU or less
  • the second HU threshold is 218HU or more and 222HU or less
  • the third HU threshold is 318HU or more and 322HU or less
  • the fourth HU The threshold value may have a value of 437HU or more and 441HU or less.
  • each of the first to fourth HU threshold values may correspond to 130HU, 200HU, 300HU, and 400HU threshold values set in a histogram of a CT image taken with a tube voltage value of 120 kV.
  • the apparatus includes a receiver configured to receive a CT image of a blood vessel, and a processor operatively connected to the receiver, wherein the processor acquires a CT image of the blood vessel using a tube voltage value set by a user, and obtains In the histogram of the CT image, pixels are arranged according to the HU value (Hounsfield Unit Value) for each pixel, the sorted pixels are counted in the order of the highest HU value, and the pixel equal to the total number of pixels up to a preset threshold is determined. and set a HU value corresponding to the determined number of pixels as a new HU threshold value.
  • HU value Hounsfield Unit Value
  • the apparatus includes a receiver configured to receive a CT image of a blood vessel, and a processor operatively connected to the receiver, wherein the processor performs computed tomography with a tube voltage value of 100 kV for a coronary artery to obtain a CT image. and calculate the coronary calcium score based on the first to fourth HU threshold values in the histogram of the acquired CT image.
  • the present invention can calculate an accurate CAC score using the HU threshold applicable under the tube voltage condition of 100 kV. In particular, since the tube voltage is lower than 120 kV, radiation exposure by CT scan can be minimized.
  • the appropriate HU threshold value according to the imaging conditions is simplified through comparison/matching with the CT image taken using the reference tube voltage.
  • the present invention does not use a complicated algorithm or a method of calculating the kV photon energy attenuation ratio, but intuitively calculates the threshold value for calculating the CAC score through the histogram pixel-HU value matching of the reference CT image and the photographed CT image. can be obtained. Since the CT scanner and the processor of the device do not need to perform complex calculations, overload of the device can be prevented and the CAC score can be calculated quickly.
  • the effect according to the present invention is not limited by the contents exemplified above, and more various effects are included in the present invention.
  • FIG. 1 is a schematic diagram for explaining the principle of a method for calculating the risk of vascular disease according to an embodiment of the present invention.
  • FIG. 2A is a schematic diagram exemplarily showing the configuration of a vascular disease risk calculation system according to an embodiment of the present invention.
  • 2B is a block diagram illustrating the configuration of an apparatus for calculating a risk of vascular disease according to an embodiment of the present invention.
  • FIG. 3 is an exemplary diagram of an interface screen output to the risk calculation device using a MATLAB program according to an embodiment of the present invention.
  • 4A to 4D are schematic diagrams exemplarily illustrating a histogram matching process according to an embodiment of the present invention.
  • FIG. 5 is a flowchart of a method for calculating a vascular disease risk according to an embodiment of the present invention.
  • FIG. 6 is a flowchart of a method of calculating a CAC score of a CT image taken at a tube voltage value of 100 kV according to an embodiment of the present invention.
  • FIG. 7 is a graph comparing a CAC score of a CT image taken under a low voltage condition with a reference CAC score according to an embodiment of the present invention.
  • expressions such as “have,” “may have,” “includes,” or “may include” refer to the presence of a corresponding characteristic (eg, a numerical value, function, operation, or component such as a part). and does not exclude the presence of additional features.
  • expressions such as “A or B,” “at least one of A or/and B,” or “one or more of A or/and B” may include all possible combinations of the items listed together.
  • “A or B,” “at least one of A and B,” or “at least one of A or B” means (1) includes at least one A, (2) includes at least one B; Or (3) it may refer to all cases including both at least one A and at least one B.
  • first may modify various elements, regardless of order and/or importance, and refer to one element. It is used only to distinguish it from other components, and does not limit the components.
  • first user equipment and the second user equipment may represent different user equipment regardless of order or importance.
  • the first component may be named as the second component, and similarly, the second component may also be renamed as the first component.
  • a component eg, a first component is "coupled with/to (operatively or communicatively)" to another component (eg, a second component);
  • another component eg, a second component
  • the certain element may be directly connected to the other element or may be connected through another element (eg, a third element).
  • a component eg, a first component
  • another component eg, a second component
  • a device configured to may mean that the device is “capable of” with other devices or parts.
  • a processor configured (or configured to perform) A, B, and C refers to a dedicated processor (eg, an embedded processor) for performing the corresponding operations, or by executing one or more software programs stored in a memory device.
  • a generic-purpose processor eg, a CPU or an application processor
  • FIG. 1 is a schematic diagram for explaining the principle of a method for calculating the risk of vascular disease according to an embodiment of the present invention.
  • a computed tomography (CT) image is acquired with a tube voltage value of 100 kV, and the CAC score is adjusted based on the existing optimal standard value. The amount of exposure can be reduced and an accurate diagnosis can be performed.
  • CT computed tomography
  • the reference value for calculating the CAC score at 100 kV is derived by comparing the result of CT scanning with the 120 kV tube voltage value, which is the existing optimal standard, by deriving the reference value calculated under the 120 kV condition. The same diagnostic result can be obtained.
  • FIG. 2A is a schematic diagram exemplarily showing the configuration of a vascular disease risk calculation system according to an embodiment of the present invention
  • FIG. 2B is a block diagram showing the configuration of an apparatus for calculating vascular disease risk according to an embodiment of the present invention.
  • the vascular disease risk calculation system 1000 includes an imaging apparatus 100 for capturing a medical image of an object 10 and an apparatus 200 for calculating a CAC score based on the medical image. may include
  • the imaging apparatus 100 includes a cylindrical bore 110 into which the object 10 is carried, and the object 10 in order to obtain a medical image of a target portion of the object 10 such as a person or an animal. It may include a transfer device 130 that is seated and carried into the bore 110 .
  • the target site includes various organs of the object 10 such as the brain, liver, and varicose veins, and in an embodiment of the present invention, the target site may be a blood vessel of the heart.
  • the imaging apparatus 100 may obtain a cross-sectional medical image in which the object 10 is horizontally cut by irradiating X-rays capable of projecting the object 10 toward the object.
  • the medical image of the object 10 obtainable by the imaging apparatus 100 may include a two-dimensional image, a three-dimensional volume image, a still image of one cut, a video composed of a plurality of cuts, and various cross-sectional images. It may include a plurality of images having
  • the vascular disease risk calculation device 200 (hereinafter referred to as “risk calculation device”) includes a receiving unit 210, an input unit 220, an output unit 230, a storage unit 240, and a processor ( 250) may be included.
  • the receiver 210 may be connected to the imaging apparatus 100 to receive a medical image of a target portion of the object 10 from the imaging apparatus 100 .
  • the receiver 210 may receive a CT image of a blood vessel, and the CT image may be an image taken using a tube voltage value of 120 kV or less.
  • the input unit 220 may receive various setting information from the user.
  • the input unit 220 may receive a photographing condition of the image photographing apparatus 100 from a user and may receive operation information of the image photographing apparatus 100 .
  • a medical practitioner may input a tube voltage value for taking a CT image using the input unit 220 , and may move the transfer device 130 to fit a target portion of the object 10 .
  • the medical practitioner may select a portion suspected of having a disease in the CT image output to the output unit 230 .
  • the output unit 230 may output an image photographed by the image photographing apparatus 100 .
  • the output unit 230 may output a CT image of a blood vessel of the object 10 .
  • the output unit 230 may display a portion in the CT image of the blood vessel in which the Hounsfield unit (HU) value is equal to or greater than a preset value.
  • HU Hounsfield unit
  • FIG. 3 is an exemplary diagram of an interface screen output to the risk calculation device using a MATLAB program according to an embodiment of the present invention.
  • the output unit 230 may output a CT image captured by the image capturing apparatus 100 .
  • the output unit 230 may display a portion having a HU value of 130HU or more in a CT image of a blood vessel in a chromatic color. Accordingly, the medical personnel can check the CT image according to the tube voltage values of 100kV and 120kV, see the area of 130HU or more displayed in chromatic colors, and diagnose whether a disease occurs.
  • the processor 250 may check and adjust the HU threshold value according to the imaging condition, where the HU threshold value is used for diagnosing a disease. It means the reference value for For example, if the HU value of the region/pixel selected by the medical practitioner is lower than the HU threshold value according to the imaging condition, the processor 250 may lower the HU threshold value to the HU value of the region/pixel selected by the medical practitioner.
  • the output unit 230 may receive a user's touch input, and in this case, the input unit 220 and the output unit 230 may be formed of one physical component. Accordingly, the output unit 230 may include various sensing means capable of recognizing various inputs generated in an image display area.
  • the storage unit 240 may store various data used in the image capturing apparatus 100 and the risk calculating apparatus 200 .
  • the storage unit 240 may store a CT image taken at a tube voltage value of 120 kV as a control, the average value of the CAC score, and the HU threshold value. , the number of pixels according to the HU value, etc. may be stored.
  • the storage unit 240 may match and store the HU threshold value for disease diagnosis with the low voltage value. That is, if the HU threshold value is initially set in the CT image taken with the low voltage value, the CAC score can be calculated based on the stored HU threshold value when shooting with the same low voltage value again.
  • the storage unit 240 may include a volatile or nonvolatile recording medium capable of storing various data, commands, and information.
  • the storage unit 240 may include a flash memory type, a hard disk type, a multimedia card micro type, a card type memory (eg, SD or XD memory, etc.), RAM, SRAM, ROM, EEPROM, PROM, network storage. It may include a storage medium of at least one type of storage, cloud, and block chain database.
  • the storage unit 240 may record instructions for operation with the risk calculating device 200 .
  • an application (not shown) for calculating the risk of vascular disease according to the tube voltage value may be recorded in the storage unit 240 .
  • the processor 250 is operatively connected to the receiving unit 210 , the input unit 220 , the output unit 230 , and the storage unit 240 to control the overall operation of the risk calculation device 200 , and the storage unit
  • various commands for setting a HU threshold value which is a criterion for diagnosing a disease, may be performed according to a tube voltage value input by a medical practitioner.
  • the processor 250 may correspond to a computing device such as a central processing unit (CPU) or an application processor (AP).
  • the processor 250 may be implemented in the form of an integrated chip (IC), such as a system on chip (SoC) in which various computing devices are integrated.
  • IC integrated chip
  • SoC system on chip
  • the processor 250 may set a new HU threshold value for diagnosing a disease based on a CT image taken at a tube voltage value of less than 120 kV. Specifically, in the case of a CT image taken at a low voltage value, since it has a larger number of pixels than necessary, the processor 250 may obtain a new HU threshold value through a process of reducing the number of pixels in the CT image, and the HU threshold value At least four values may be set to divide the HU value into four sections.
  • the processor 250 may acquire a CT image of a blood vessel using a tube voltage value set by a user, and may align pixels according to each HU value in a histogram of the acquired CT image.
  • the processor 250 may count the sorted pixels in the order of increasing the HU value, and determine a pixel equal to the total number of pixels up to a preset threshold. That is, the processor 250 may adjust the number of pixels for histogram matching with the optimal standard value, and accordingly, the preset threshold is 130HU, which is the minimum HU value that is the standard for diagnosing diseases in the CT image taken with the 120kV tube voltage value.
  • 4A to 4D are schematic diagrams exemplarily illustrating a histogram matching process according to an embodiment of the present invention.
  • each graph is an arrangement of pixels according to HU values in a histogram of a CT image taken with 120 kV and 80 kV tube voltage values.
  • the total number of pixels may be 260,104, and in the 80kV imaging condition, the total number of pixels may be 280,737.
  • values corresponding to the graph A may be previously stored as reference values in the storage unit 240 .
  • the processor 250 may sequentially count up to 260 and 104 from the highest HU value in the CT image taken at the tube voltage value of 80 kV. Referring to the graph of FIG. 4B , the processor 250 counts inversely from the number of pixels having a value of approximately 1000HU, but may end counting at 260 and 104 th. Accordingly, after counting 260,104 pixels, the processor 250 excludes the remaining 20,633 pixels from the data for calculating the CAC score, and sets the HU value of the 260,104th pixel as the minimum HU threshold that is the standard for diagnosing the disease. can For example, if the HU value of the 260,104th pixel is 151HU, the minimum HU threshold value for calculating the CAC score in the CT image taken at the 80kV tube voltage value may be 151HU.
  • the processor 250 may set a new first HU threshold value for the HU value of the last pixel as the pixels are counted in the reverse order. can be matched with
  • the processor 250 sets the remaining HU threshold values sequentially based on the newly set minimum HU value.
  • the reference HU threshold values (130HU, 200HU, 300HU, 400HU) shown in circles in FIG. 4A are changed and set as shown in red circles in FIG. 4D. can be
  • the processor 250 may calculate a CAC score for diagnosing a coronary artery disease (eg, coronary atherosclerosis) based on a CT image taken at a tube voltage value of 100 kV among tube voltage values of less than 120 kV. Specifically, the processor 250 acquires a CT image of the coronary artery photographed with a tube voltage value of 100 kV from the imaging apparatus 100, and based on the first to fourth HU threshold values in the histogram of the CT image, the coronary artery Calcium score can be calculated.
  • a coronary artery disease eg, coronary atherosclerosis
  • the first HU threshold is 141HU or more and 145HU or less
  • the second HU threshold is 218HU or more and 222HU or less
  • the third HU threshold is 318HU or more and 322HU or less
  • the fourth HU threshold is 437HU or more and 441HU or less , which may correspond to the 130HU, 200HU, 300HU, and 400HU threshold values set in the histogram of the CT image taken with the tube voltage value of 120 kV.
  • the apparatus 200 for calculating the risk of vascular disease according to an embodiment of the present invention has been described. According to the present invention, even if the risk calculation device 200 takes a CT image using a tube voltage lower than 120 kV, such as 70 kV, 80 kV, etc., through comparison/matching with the CT image taken using the reference tube voltage, the An appropriate HU threshold can be obtained simply.
  • FIG. 5 is a flowchart of a method for calculating a vascular disease risk according to an embodiment of the present invention.
  • the risk calculation apparatus 200 obtains a CT image of a blood vessel using a tube voltage value set by a user ( S110 ).
  • the risk calculation apparatus 200 may transmit the tube voltage value to the image photographing apparatus 100 and receive a CT image according to the specified tube voltage value from the image photographing apparatus 100 .
  • the risk calculation apparatus 200 aligns pixels according to each HU value in the histogram of the acquired CT image (S120). In other words, the risk calculating device 200 may count the number of pixels having the same HU value. For example, as the HU value increases as a result of the alignment by the risk calculation device 200, the number of pixels corresponding thereto may be small.
  • the risk calculating device 200 counts the sorted pixels in the order of the HU value, and determines the pixels equal to the total number of pixels up to a preset threshold (S130). Specifically, in the case of a CT image taken at a low voltage value, it is common to have a larger number of pixels than necessary. Accordingly, the risk calculation apparatus 200 may adjust the number of pixels to match the histogram with the optimal standard value.
  • the preset threshold may be 130HU, which is the minimum HU value serving as a criterion for diagnosing a disease in a CT image taken with a 120kV tube voltage value.
  • the risk calculation device 200 sets the HU value corresponding to the determined number of pixels as a new HU threshold value (S140).
  • the new HU threshold set through the number of pixels is the first reference value for calculating the CAC score.
  • the new HU threshold value may have a value higher than 130HU, the minimum HU threshold value at 120 kV, which is the optimal standard. .
  • the risk calculation device 200 may sequentially set the remaining HU threshold values based on the set minimum HU value, and finally obtain a new HU threshold value for dividing the histogram into four sections.
  • the threshold value for calculating the CAC score is intuitively determined through the histogram pixel-HU value matching of the reference CT image and the photographed CT image without using a complicated algorithm or a method of calculating the kV photon energy attenuation ratio.
  • FIG. 6 is a flowchart of a method of calculating a CAC score of a CT image taken at a tube voltage value of 100 kV according to an embodiment of the present invention.
  • the risk calculation apparatus 200 obtains a CT image of the coronary artery photographed with a tube voltage value of 100 kV from the imaging apparatus 100 ( S210 ). To this end, the risk calculation apparatus 200 may transmit information for controlling components in the image capturing apparatus 100 .
  • the risk calculating device 200 calculates a coronary calcium score based on the first to fourth HU threshold values in the histogram of the CT image (S220).
  • the first HU threshold is 141HU or more and 145HU or less
  • the second HU threshold is 218HU or more and 222HU or less
  • the third HU threshold is 318HU or more and 322HU or less
  • the fourth HU threshold is 437HU or more and 441HU or less can have a value of
  • the risk calculation device 200 may calculate a score according to the new HU threshold using the Agatston score calculation method used to calculate the coronary calcium score under the condition of 120 kV.
  • the risk calculation device 200 has a value of 1, 221HU or more and 320HU or less in the case of a value of 143HU or more and 220HU or less, 2 if it has a value of 321HU or more and 439HU or less If it has a value of 3, 440HU or more 4 , and the coronary calcium score can be calculated according to a formula such as “ ⁇ (area x concentration factor x section thickness/3) of each selected site”.
  • FIG. 7 is a graph comparing a CAC score of a CT image taken under a low voltage condition with a reference CAC score according to an embodiment of the present invention.
  • the CAC score was calculated using the vascular disease risk calculation method according to an embodiment of the present invention after taking pictures at a low voltage of less than 120 kV and a reference value of 120 kV for the same blood vessel. , it can be seen that the values are almost the same.
  • an accurate CAC score can be calculated using the HU threshold value applicable under the tube voltage condition of 100 kV.
  • the tube voltage condition is lower than 120 kV, radiation exposure by CT scan can be minimized.

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Physics & Mathematics (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biophysics (AREA)
  • Optics & Photonics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Dentistry (AREA)
  • Vascular Medicine (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Quality & Reliability (AREA)
  • Pulmonology (AREA)
  • Physiology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

본 발명은, 혈관 질환 위험도 계산방법에 관한 것으로서, 상기 계산방법은, 사용자에 의해 설정된 관전압 값을 이용하여 혈관에 대한 CT 영상을 획득하는 단계, 획득된 CT 영상의 히스토그램에서 픽셀 별 HU 값(Hounsfield Unit Value)에 따라 픽셀들을 정렬하는 단계, 정렬된 픽셀들을 상기 HU 값이 높은 순서대로 카운트하고, 기 설정된 임계치까지의 총 픽셀 개수와 동일해지는 픽셀을 결정하는 단계 및 결정된 픽셀의 개수와 대응되는 HU 값을 새로운 HU 임계 값으로 설정하는 단계를 포함하도록 구성된다.

Description

저전압을 이용한 CT 영상에서 혈관 질환 위험도 계산 방법 및 장치
본 발명은 저전압을 이용한 CT 영상에서 혈관 질환 위험도 계산 방법 및 장치에 관한 것이다. 보다 구체적으로는, 기준 전압보다 낮은 전압에서 관상동맥 칼슘 점수(Coronary Artery Calcium Score, CACS)를 정확하게 계산하는 방법과 이를 수행하는 장치, CT 스캐너에 관한 것이다.
컴퓨터 단층촬영(Computed Tomography, CT)을 통해 확인 가능한 CAC 점수는 심혈관 질환, 관상동맥 질환에 대한 가장 신뢰할 수 있는 비 침습적 예측 인자이다.
일반적으로, CT 촬영을 통해 얻는 CAC 점수는 ∑(각 선택 부위의 면적x농도인자 x절편두께/3)으로 계산되는 Agatston Score로 정의되었으며, 해당 계산식에 적용하기 위한 농도인자는 CT 스캐너에 120kV의 관전압을 인가하였을 때 얻어지는 CT 영상을 기반으로 결정된다.
최근, 영상의학검사에서 방사선 노출을 점차 낮추기 위한 다양한 연구가 진행됨에 따라, CAC 점수를 계산하는 방법에서도 기존의 최적화된 기준(120kV) 보다 낮은 전압에서 CAC 점수를 계산하고자 하였으나, 낮은 전압에서는 실제 CAC 점수 보다 과장된 점수로 계산되어, 치료가 불필요한 사람들이 치료 대상이 되는 문제점이 있다.
발명의 배경이 되는 기술은 본 발명에 대한 이해를 보다 용이하게 하기 위해 작성되었다. 발명의 배경이 되는 기술에 기재된 사항들이 선행기술로 존재한다고 인정하는 것으로 이해되어서는 안 된다.
이에, CAC 점수를 계산하는 가장 최적의 전압 값보다 낮은 전압 값에서 CAC 점수를 정확하게 계산할 수 있는 방법이 요구되고 있는 실정이다.
그 결과, 본 발명의 발명자들은, 120kV 보다 낮은 관전압에서 촬영된 CT 영상으로도, 기준 관전압(120kV 관전압) 조건에서 계산된 CAC 점수와 동일하게 CAC 점수 계산이 가능한 임계값을 얻을 수 있는 새로운 혈관 질환 위험도 계산 방법 및 장치를 개발하고자 하였다.
이때, 본 발명의 발명자들은, 기준 CT 영상과 보다 낮은 관전압에서 촬영된 CT 영상의 히스토그램을 비교/매칭하여, CAC 점수를 계산하기 위한 임계값을 구할 수 있도록 방법 및 장치를 구성하였다.
특히, 본 발명의 발명자들은, 관전압의 크기와 비례하는 방사선 피폭량을 줄이기 위해, 100kV의 관전압을 이용하여 CT 영상을 찍고, 이를 기초로 환자의 관상동맥 칼슘 점수를 구할 수 있는 방법 및 장치를 개발하기에 이르렀다.
본 발명의 과제들은 이상에서 언급한 과제들로 제한되지 않으며, 언급되지 않은 또 다른 과제들은 아래의 기재로부터 당업자에게 명확하게 이해될 수 있을 것이다.
전술한 바와 같은 과제를 해결하기 위하여 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산방법을 제공한다. 상기 방법은, 사용자에 의해 설정된 관전압 값을 이용하여 혈관에 대한 CT 영상을 획득하는 단계, 획득된 CT 영상의 히스토그램에서 픽셀 별 HU 값(Hounsfield Unit Value)에 따라 픽셀들을 정렬하는 단계, 정렬된 픽셀들을 상기 HU 값이 높은 순서대로 카운트하고, 기 설정된 임계치까지의 총 픽셀 개수와 동일해지는 픽셀을 결정하는 단계 및 결정된 픽셀의 개수와 대응되는 HU 값을 새로운 HU 임계 값으로 설정하는 단계를 포함하도록 구성된다.
본 발명의 특징에 따르면, 상기 사용자에 의해 설정된 관전압 값은, 기준 CT 영상을 획득하기 위해 설정된 기준 관전압 값보다 낮은 값을 가질 수 있다.
본 발명의 다른 특징에 따르면, 상기 기준 CT 영상의 히스토그램은, 기 설정된 4개의 임계 값을 기준으로 적어도 4개의 구간으로 나누어지며, 상기 새로운 HU 임계 값으로 설정하는 단계는, 상기 사용자에 의해 설정된 관전압에 따른 제1 내지 제4 HU 임계 값을 설정하는 단계일 수 있다.
본 발명의 또 다른 특징에 따르면, 상기 새로운 HU 임계 값으로 설정하는 단계는, 상기 추출된 픽셀에서 가장 낮은 HU 값을 상기 제1 HU 임계 값으로 설정하는 단계와 상기 제1 HU 임계 값을 기초로 상기 제2 내지 제4 HU 임계 값을 순차적으로 설정하는 단계를 더 포함할 수 있다.
본 발명의 또 다른 특징에 따르면, 상기 사용자에 의해 설정된 관전압 값은, 120kV 미만의 관전압 값일 수 있다.
본 발명의 또 다른 특징에 따르면, 상기 새로운 HU 임계 값으로 설정하는 단계는, 상기 획득된 CT 영상의 히스토그램에서 상기 새로운 HU 임계 값에 따라 구분된 픽셀들을 표시하는 단계를 더 포함할 수 있다.
전술한 바와 같은 과제를 해결하기 위하여 본 발명의 다른 실시 예에 따른 혈관 질환 위험도 계산방법을 제공한다. 상기 방법은, 관상동맥에 대하여 100kV의 관전압 값으로 컴퓨터 단층촬영을 수행하여, CT 영상을 획득하는 단계 및 획득된 CT 영상의 히스토그램에서 제1 내지 제4 HU 임계 값을 기초로 상기 관상동맥 칼슘 점수를 계산하는 단계를 포함하도록 구성된다.
본 발명의 특징에 따르면, 상기 제1 HU 임계 값은, 141HU 이상 145HU 이하, 상기 제2 HU 임계 값은, 218HU 이상 222HU 이하, 상기 제3 HU 임계 값은, 318HU 이상 322HU 이하, 상기 제4 HU 임계 값은, 437HU 이상 441HU 이하의 값을 가질 수 있다.
본 발명의 다른 특징에 따르면, 상기 제1 내지 제4 HU 임계 값 각각은, 120kV의 관전압 값으로 촬영된 CT 영상의 히스토그램에 설정된 130HU, 200HU, 300HU 및 400HU 임계 값과 대응될 수 있다.
전술한 바와 같은 과제를 해결하기 위하여 본 발명의 또 다른 실시 예에 따른 혈관 질환 위험도 계산 장치를 제공한다. 상기 장치는, 혈관에 대한 CT 영상을 수신하도록 구성된 수신부 및 상기 수신부와 동작 가능하게 연결된 프로세서를 포함하고, 상기 프로세서는, 사용자에 의해 설정된 관전압 값을 이용하여 혈관에 대한 CT 영상을 획득하고, 획득된 CT 영상의 히스토그램에서 픽셀 별 HU 값(Hounsfield Unit Value)에 따라 픽셀들을 정렬하며, 정렬된 픽셀들을 상기 HU 값이 높은 순서대로 카운트하고, 기 설정된 임계치까지의 총 픽셀 개수와 동일해지는 픽셀을 결정하고, 결정된 픽셀의 개수와 대응되는 HU 값을 새로운 HU 임계 값으로 설정하도록 구성된다.
전술한 바와 같은 과제를 해결하기 위하여 본 발명의 또 다른 실시 예에 따른 혈관 질환 위험도 계산 장치를 제공한다. 상기 장치는, 혈관에 대한 CT 영상을 수신하도록 구성된 수신부 및 상기 수신부와 동작 가능하게 연결된 프로세서를 포함하고, 상기 프로세서는, 관상동맥에 대하여 100kV의 관전압 값으로 컴퓨터 단층촬영을 수행하여, CT 영상을 획득하고, 획득된 CT 영상의 히스토그램에서 제1 내지 제4 HU 임계 값을 기초로 상기 관상동맥 칼슘 점수를 계산하도록 구성된다.
기타 실시 예의 구체적인 사항들은 상세한 설명 및 도면들에 포함되어 있다.
본 발명은 100kV의 관전압 조건에서 적용 가능한 HU 임계 값을 이용하여 정확한 CAC 점수를 계산할 수 있다. 특히, 120kV 보다 낮은 관전압 조건이므로, CT 촬영에 의한 방사선 노출을 최소화할 수 있다.
나아가, 100kV가 아닌, 70, 80kV 등 120kV 보다 낮은 관전압을 이용하여 CT 영상을 찍더라도, 기준 관전압을 이용하여 촬영된 CT 영상과의 비교/매칭을 통해, 촬영 조건에 따른 적절한 HU 임계 값을 간단하게 얻을 수 있다.
또한, 본 발명은, 복잡한 알고리즘이나, kV 광자 에너지 감쇠 비율을 계산 방법을 이용하지 않고, 기준 CT 영상과 촬영된 CT 영상의 히스토그램 픽셀-HU 값 매칭을 통해 CAC 점수 계산을 위한 임계값을 직관적으로 얻을 수 있다. 특히, CT 스캐너, 장치의 프로세서가 복잡한 연산을 수행하지 않아도 되므로, 장치의 과부하를 예방하고, CAC 점수를 빠르게 계산해낼 수 있다.
본 발명에 따른 효과는 이상에서 예시된 내용에 의해 제한되지 않으며, 더욱 다양한 효과들이 본 발명 내에 포함되어 있다.
도 1은 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산방법의 원리를 설명하기 위한 개략도이다.
도 2a는 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산 시스템의 구성을 예시적으로 나타낸 개략도이다.
도 2b는 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산 장치의 구성을 나타낸 블록도이다.
도 3은 본 발명의 일 실시 예에 따른 매트랩(MATLAB) 프로그램을 이용하여 위험도 계산 장치에 출력되는 인터페이스 화면의 예시도이다.
도 4a 내지 4d는 본 발명의 일 실시 예에 따른 히스토그램 매칭 과정을 예시적으로 나타낸 개략도이다.
도 5는 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산 방법의 순서도이다.
도 6은 본 발명의 일 실시 예에 따른 100kV의 관전압 값에서 촬영한 CT 영상의 CAC 점수를 계산하는 방법의 순서도이다.
도 7은 본 발명의 일 실시 예에 따른 저전압 조건에서 촬영된 CT 영상의 CAC 점수를 기준 CAC 점수와 비교한 그래프이다.
본 발명의 이점 및 특징, 그리고 그것들을 달성하는 방법은 첨부되는 도면과 함께 상세하게 후술되어 있는 실시 예들을 참조하면 명확해질 것이다. 그러나, 본 발명은 이하에서 개시되는 실시 예들에 한정되는 것이 아니라 서로 다른 다양한 형태로 구현될 것이며, 단지 본 실시 예들은 본 발명의 개시가 완전하도록 하며, 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자에게 발명의 범주를 완전하게 알려주기 위해 제공되는 것이며, 본 발명은 청구항의 범주에 의해 정의될 뿐이다. 도면의 설명과 관련하여, 유사한 구성요소에 대해서는 유사한 참조부호가 사용될 수 있다.
본 문서에서, "가진다," "가질 수 있다," "포함한다," 또는 "포함할 수 있다" 등의 표현은 해당 특징(예: 수치, 기능, 동작, 또는 부품 등의 구성요소)의 존재를 가리키며, 추가적인 특징의 존재를 배제하지 않는다.
본 문서에서, "A 또는 B," "A 또는/및 B 중 적어도 하나," 또는 "A 또는/및 B 중 하나 또는 그 이상" 등의 표현은 함께 나열된 항목들의 모든 가능한 조합을 포함할 수 있다. 예를 들면, "A 또는 B," "A 및 B 중 적어도 하나," 또는 "A 또는 B 중 적어도 하나"는, (1) 적어도 하나의 A를 포함, (2) 적어도 하나의 B를 포함, 또는(3) 적어도 하나의 A 및 적어도 하나의 B 모두를 포함하는 경우를 모두 지칭할 수 있다.
본 문서에서 사용된 "제1," "제2," "첫째," 또는 "둘째," 등의 표현들은 다양한 구성요소들을, 순서 및/또는 중요도에 상관없이 수식할 수 있고, 한 구성요소를 다른 구성요소와 구분하기 위해 사용될 뿐 해당 구성요소들을 한정하지 않는다. 예를 들면, 제1 사용자 기기와 제2 사용자 기기는, 순서 또는 중요도와 무관하게, 서로 다른 사용자 기기를 나타낼 수 있다. 예를 들면, 본 문서에 기재된 권리범위를 벗어나지 않으면서 제1 구성요소는 제2 구성요소로 명명될 수 있고, 유사하게 제2 구성요소도 제1 구성요소로 바꾸어 명명될 수 있다.
어떤 구성요소(예: 제1 구성요소)가 다른 구성요소(예: 제2 구성요소)에 "(기능적으로 또는 통신적으로) 연결되어((operatively or communicatively) coupled with/to)" 있다거나 "접속되어(connected to)" 있다고 언급된 때에는, 상기 어떤 구성요소가 상기 다른 구성요소에 직접적으로 연결되거나, 다른 구성요소(예: 제3 구성요소)를 통하여 연결될 수 있다고 이해되어야 할 것이다. 반면에, 어떤 구성요소(예: 제1 구성요소)가 다른 구성요소(예: 제2 구성요소)에 "직접 연결되어" 있다거나 "직접 접속되어" 있다고 언급된 때에는, 상기 어떤 구성요소와 상기 다른 구성요소 사이에 다른 구성요소(예: 제3 구성요소)가 존재하지 않는 것으로 이해될 수 있다.
본 문서에서 사용된 표현 "~하도록 구성된(또는 설정된)(configured to)"은 상황에 따라, 예를 들면, "~에 적합한(suitable for)," "~하는 능력을 가지는(having the capacity to)," "~하도록 설계된(designed to)," "~하도록 변경된(adapted to)," "~하도록 만들어진(made to)," 또는 "~ 를 할 수 있는(capable of)"과 바꾸어 사용될 수 있다. 용어 "~하도록 구성된(또는 설정된)"은 하드웨어적으로 "특별히 설계된(specifically designed to)" 것만을 반드시 의미하지 않을 수 있다. 대신, 어떤 상황에서는, "~하도록 구성된 장치"라는 표현은, 그 장치가 다른 장치 또는 부품들과 함께 "~할 수 있는" 것을 의미할 수 있다. 예를 들면, 문구 "A, B, 및 C를 수행하도록 구성된(또는 설정된)프로세서"는 해당 동작을 수행하기 위한 전용 프로세서(예: 임베디드 프로세서), 또는 메모리 장치에 저장된 하나 이상의 소프트웨어 프로그램들을 실행함으로써, 해당 동작들을 수행할 수 있는 범용 프로세서(generic-purpose processor)(예: CPU 또는 application processor)를 의미할 수 있다.
본 문서에서 사용된 용어들은 단지 특정한 실시 예를 설명하기 위해 사용된 것으로, 다른 실시 예의 범위를 한정하려는 의도가 아닐 수 있다. 단수의 표현은 문맥상 명백하게 다르게 뜻하지 않는 한, 복수의 표현을 포함할 수 있다. 기술적이거나 과학적인 용어를 포함해서 여기서 사용되는 용어들은 본 문서에 기재된 기술분야에서 통상의 지식을 가진 자에 의해 일반적으로 이해되는 것과 동일한 의미를 가질 수 있다. 본 문서에 사용된 용어들 중 일반적인 사전에 정의된 용어들은, 관련 기술의 문맥상 가지는 의미와 동일 또는 유사한 의미로 해석될 수 있으며, 본 문서에서 명백하게 정의되지 않는 한, 이상적이거나 과도하게 형식적인 의미로 해석되지 않는다. 경우에 따라서, 본 문서에서 정의된 용어일지라도 본 문서의 실시 예들을 배제하도록 해석될 수 없다.
본 발명의 여러 실시 예들의 각각 특징들이 부분적으로 또는 전체적으로 서로 결합 또는 조합 가능하며, 당업자가 충분히 이해할 수 있듯이 기술적으로 다양한 연동 및 구동이 가능하며, 각 실시 예들이 서로에 대하여 독립적으로 실시 가능할 수도 있고 연관 관계로 함께 실시 가능할 수도 있다.
본 명세서의 해석의 명확함을 위해, 이하에서는 본 명세서에서 사용되는 용어들을 정의하기로 한다.
도 1은 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산방법의 원리를 설명하기 위한 개략도이다.
도 1을 참조하면, 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산방법에서는 100kV의 관전압 값으로 CT(Computed Tomography) 영상을 획득하고, CAC 점수를 기존의 최적 표준 값을 기준으로 조정함으로써, 방사선 노출양은 줄이고 정확한 진단을 수행할 수 있다.
구체적으로, 본 발명의 일 실시 예에서는 기존의 최적 표준인 120kV 관전압 값으로 CT 스캐닝 수행한 결과와의 비교를 통해, 100kV에서 CAC 점수를 계산하기 위한 기준 값을 도출해냄으로써, 120kV 조건에서 계산한 것과 같은 진단 결과를 얻을 수 있다.
도 2a는 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산 시스템의 구성을 예시적으로 나타낸 개략도이고, 도 2b는 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산 장치의 구성을 나타낸 블록도이다.
도 2a를 참조하면, 혈관 질환 위험도 계산 시스템(1000)은 대상체(10)의 의료 영상을 촬영하는 영상 촬영 장치(100) 및 의료 영상을 기초로 CAC 점수를 계산하는 혈관 질환 위험도 계산 장치(200)를 포함할 수 있다.
영상 촬영 장치(100)는 사람, 동물 등의 대상체(10)의 목적 부위에 대한 의료 영상을 획득하기 위해, 대상체(10)가 반입되는 원통형의 보어(bore)(110), 대상체(10)가 안착되고 보어(110) 내부로 반입시키는 이송장치(130)를 포함할 수 있다. 여기서, 목적 부위는 뇌, 간, 정맥류 등 대상체(10)의 다양한 장기를 포함하며, 본 발명의 일 실시 예에서 목적 부위는 심장의 혈관일 수 있다.
또한, 영상 촬영 장치(100)는 대상체(10)를 투영할 수 있는 엑스레이(X-Ray)를 피사체를 향하여 조사함으로써, 대상체(10)를 가로로 자른 횡단면 의료 영상을 획득할 수 있다. 다만 실시 예에 따라, 영상 촬영 장치(100)가 획득할 수 있는 대상체(10)에 대한 의료 영상은 2차원 영상, 3차원 볼륨 영상, 한 컷의 스틸 영상, 복수 개의 컷으로 구성된 동영상, 다양한 단면상을 가지는 복수 개의 영상을 포함할 수 있다.
도 2b를 참조하면, 혈관 질환 위험도 계산 장치(200)(이하, “위험도 계산 장치”라 함)는 수신부(210), 입력부(220), 출력부(230), 저장부(240) 및 프로세서(250)를 포함할 수 있다.
수신부(210)는 영상 촬영 장치(100)와 연결되어, 영상 촬영 장치(100)로부터 대상체(10)의 목적 부위에 대한 의료 영상을 수신할 수 있다. 실시 예에 따라, 수신부(210)는 혈관에 대한 CT 영상을 수신할 수 있으며, 여기서의 CT 영상은 120kV 이하의 관전압 값을 이용하여 촬영된 영상일 수 있다.
입력부(220)는 사용자로부터 다양한 설정 정보를 입력 받을 수 있다. 입력부(220)는 사용자로부터 영상 촬영 장치(100)의 촬영 조건을 입력 받고, 영상 촬영 장치(100)의 동작 정보를 입력 받을 수 있다. 예를 들어, 의료인은 입력부(220)를 이용하여 CT 영상을 촬영하기 위한 관전압 값을 입력하고, 대상체(10)의 목적 부위에 맞게 이송장치(130)를 이동시킬 수 있다. 다른 예를 들어, 의료인은 출력부(230)에 출력되는 CT 영상에서 질환이 의심되는 부분을 선택할 수 있다.
출력부(230)는 영상 촬영 장치(100)가 촬영한 영상을 출력할 수 있다. 실시 예에 따라, 출력부(230)는 대상체(10)의 혈관에 대한 CT 영상을 출력할 수 있다. 출력부(230)는 혈관에 대한 CT 영상에서 하운스필드 유닛(Hounsfield unit, HU) 값이 기 설정된 값 이상인 부분을 표시할 수 있다.
도 3은 본 발명의 일 실시 예에 따른 매트랩(MATLAB) 프로그램을 이용하여 위험도 계산 장치에 출력되는 인터페이스 화면의 예시도이다.
도 3을 참조하면, 출력부(230)는 영상 촬영 장치(100)에서 촬영된 CT 영상을 출력할 수 있다. 이때, 출력부(230)에는 혈관에 대한 CT 영상에서 HU 값이 130HU 이상인 부분을 유채색으로 표시될 수 있다. 그에 따라, 의료인은 100kV, 120kV의 관전압 값에 따른 CT 영상을 확인하고, 이 중 유채색으로 표시된 130HU 이상의 영역을 보고, 질환의 발생 여부를 진단할 수 있다.
만약, 의료인이 입력부(220)를 이용하여 해당 영역에 질환이 발생한 것으로 체크한다면, 프로세서(250)는 촬영 조건에 따른 HU 임계 값을 확인 및 조정할 수 있으며, 여기서 HU 임계 값이란, 질환을 진단하기 위한 기준 값을 의미한다. 예를 들어, 의료인이 선택한 영역/픽셀의 HU 값이 촬영 조건에 따른 HU 임계 값 보다 낮다면, 프로세서(250)는 HU 임계 값을 의료인이 선택한 영역/픽셀의 HU 값으로 하향 조정할 수 있다.
다시 도 2b를 참조하여 설명하도록 한다.
출력부(230)는 사용자의 터치 입력을 수신할 수 있으며, 이러한 경우, 입력부(220)와 출력부(230)는 하나의 물리적 구성 요소로 이루어질 수 있다. 그에 따라, 출력부(230)는 영상을 표시하는 영역에서 발생하는 다양한 입력을 인식할 수 있는 다양한 센싱 수단들을 포함할 수 있다.
저장부(240)는 영상 촬영 장치(100) 및 위험도 계산 장치(200)에서 사용되는 다양한 데이터를 저장할 수 있다. 예를 들어, 저장부(240)는 120kV 보다 낮은 관전압 값에서 촬영된 CT 영상의 HU 임계 값을 구하기 위해, 대조군이 되는 120kV의 관전압 값에서 촬영된 CT 영상, CAC 점수의 평균 값, HU 임계 값, HU 값에 따른 픽셀 개수 등을 저장할 수 있다.
또한, 저장부(240)는 저전압 값과 함께 질환 진단을 위한 HU 임계 값을 매칭하여 저장할 수 있다. 즉, 저전압 값으로 촬영된 CT 영상에서 HU 임계 값을 최초 설정하였다면, 이후 또 다시 동일한 저전압 값으로의 촬영 시, 저장해둔 HU 임계 값을 기초로 CAC 점수를 계산할 수 있다.
다양한 실시 예에서 저장부(240)는 각종 데이터, 명령 및 정보를 저장할 수 있는 휘발성 또는 비휘발성 기록 매체를 포함할 수 있다. 예를 들어, 저장부(240)는 플래시 메모리 타입, 하드디스크 타입, 멀티미디어 카드 마이크로 타입, 카드 타입의 메모리(예를 들어 SD 또는 XD 메모리 등), 램, SRAM, 롬, EEPROM, PROM, 네트워크 저장 스토리지, 클라우드, 블록체인 데이터베이스 중 적어도 하나의 타입의 저장매체를 포함할 수 있다.
또한, 저장부(240)는 위험도 계산 장치(200)와의 동작을 위한 명령어들이 기록되어 있을 수 있다. 다양한 실시 예에서, 저장부(240)는 관전압 값에 따른 혈관 질환 위험도를 계산하기 위한 어플리케이션(미도시)이 기록되어 있을 수 있다.
프로세서(250)는 수신부(210), 입력부(220), 출력부(230) 및 저장부(240)와 동작 가능하게 연결되어, 위험도 계산 장치(200)의 전반적인 동작을 제어할 수 있으며, 저장부(240)에 저장된 어플리케이션 또는 프로그램을 구동하여 의료인이 입력한 관전압 값에 따라 질환을 진단하는 기준이 되는 HU 임계 값을 설정하기 위한 다양한 명령들을 수행할 수 있다.
한편, 프로세서(250)는 CPU(Central Processing Unit)나 AP(Application Processor)와 같은 연산 장치에 해당할 수 있다. 또한, 프로세서(250)는 다양한 연산 장치가 통합된 SoC(System on Chip)와 같은 통합 칩(Integrated Chip (IC))의 형태로 구현될 수 있다.
실시 예에 따라, 프로세서(250)는 120kV 미만의 관전압 값에서 촬영한 CT 영상을 기초로 질환을 진단하기 위한 새로운 HU 임계 값을 설정할 수 있다. 구체적으로, 저전압 값에서 촬영한 CT 영상의 경우, 필요 이상의 많은 픽셀 수를 가지고 있으므로, 프로세서(250)는 CT 영상 내 픽셀의 수를 줄이는 과정을 통해서 새로운 HU 임계 값을 얻을 수 있으며, HU 임계 값은 HU 값을 4구간으로 나누기 위해 적어도 4개의 값이 설정될 수 있다.
보다 구체적으로, 프로세서(250)는 사용자에 의해 설정된 관전압 값을 이용하여 혈관에 대한 CT 영상을 획득하고, 획득된 CT 영상의 히스토그램에서 각각의 HU 값에 따라 픽셀들을 정렬할 수 있다. 프로세서(250)는 정렬된 픽셀들을 HU 값이 높은 순서대로 카운트하고, 기 설정된 임계치까지의 총 픽셀 개수와 동일해지는 픽셀을 결정할 수 있다. 즉, 프로세서(250)는 최적 표준 값과의 히스토그램 매칭을 위해 픽셀의 개수를 조정할 수 있으며 그에 따라, 기 설정된 임계치는 120kV 관전압 값으로 촬영된 CT 영상에서 질환을 진단하는 기준이 되는 최소 HU 값인 130HU 일 수 있다.
도 4a 내지 도 4d는 본 발명의 일 실시 예에 따른 히스토그램 매칭 과정을 예시적으로 나타낸 개략도이다.
도 4a 및 4b를 참조하면, 각각의 그래프는 120kV, 80kV 관전압 값으로 촬영된 CT 영상의 히스토그램에서 HU 값에 따라 픽셀들을 정렬한 것이다. 120kV 촬영 조건에서는 픽셀 개수가 총 260,104개이고, 80kV 촬영 조건에서는 픽셀 개수가 총 280,737개일 수 있다.
프로세서(250)가 새로운 HU 임계 값을 설정하기 전에, 저장부(240)에는 미리 그래프 A에 해당하는 값들이 기준 값들로 저장되어 있을 수 있다.
프로세서(250)는 80kV의 관전압 값에서 촬영한 CT 영상에서 HU 값이 높은 순서부터 순차적으로 260,104까지 카운트할 수 있다. 도 4b의 그래프를 보면, 프로세서(250)는 대략 1000HU 값을 가지는 픽셀 개수부터 역으로 카운트하되, 260,104번째에서 카운트를 종료할 수 있다. 그에 따라, 프로세서(250)는 260,104개를 카운트 한 후 나머지 20,633개의 픽셀은 CAC 점수 계산을 위한 데이터에서 제외시키고, 260,104번째 픽셀의 HU 값을 질환을 진단하는 기준이 되는 최소 HU 임계 값으로 새롭게 설정할 수 있다. 예를 들어, 260,104번째 픽셀의 HU 값이 151HU 이라면, 80kV 관전압 값에서 촬영되는 CT 영상에서 CAC 점수를 계산하기 위한 최소 HU 임계 값은 151HU가 될 수 있다.
즉, 프로세서(250)는 픽셀을 역순으로 카운트함에 따라 가장 마지막에 해당하는 픽셀의 HU 값을 새로운 제1 HU 임계 값을 설정할 수 있으며 여기서의 제1 임계 값은 질환을 진단하는 기준이 되는 130HU 값과 대응될 수 있다.
도 4c의 그래프를 살펴보면, 히스토그램 매칭에 따라, 동일 HU 값을 가지는 픽셀의 개수가 상이한 것을 확인할 수 있으며, 프로세서(250)는 새롭게 설정한 최소 HU 값을 기준으로 나머지 HU 임계 값들을 순차적으로 새로이 설정할 수 있다.
도 4d의 그래프를 살펴보면, 프로세서(250)의 히스토그램 매칭에 따라 도 4a에 원으로 도시된 기준 HU 임계 값(130HU, 200HU, 300HU, 400HU)들이 도 4d에 빨간색 원으로 도시된 바와 같이 변경되어 설정될 수 있다.
다시 도 2b를 참조하여 설명하도록 한다.
실시 예에 따라, 프로세서(250)는 120kV 미만의 관전압 값 중에서도 100kV의 관전압 값에서 촬영한 CT 영상을 기초로 관상동맥 질환(예. 관상동맥경화증)을 진단하기 위한 CAC 점수를 계산할 수 있다. 구체적으로, 프로세서(250)는 영상 촬영 장치(100)로부터 100kV의 관전압 값으로 촬영된 관상동맥에 대한 CT 영상을 획득하고, CT 영상의 히스토그램에서 제1 내지 제4 HU 임계 값을 기초로 관상동맥 칼슘 점수를 계산할 수 있다.
구체적으로, 제1 HU 임계 값은, 141HU 이상 145HU 이하, 제2 HU 임계 값은, 218HU 이상 222HU 이하, 제3 HU 임계 값은, 318HU 이상 322HU 이하, 제4 HU 임계 값은, 437HU 이상 441HU 이하의 값을 가지며, 이는 120kV의 관전압 값으로 촬영된 CT 영상의 히스토그램에 설정된 130HU, 200HU, 300HU 및 400HU 임계 값과 대응될 수 있다.
지금까지 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산 장치(200)에 대하여 설명하였다. 본 발명에 따르면, 위험도 계산 장치(200)는 70kV, 80kV 등 120kV 보다 낮은 관전압을 이용하여 CT 영상을 찍더라도, 기준 관전압을 이용하여 촬영된 CT 영상과의 비교/매칭을 통해, 촬영 조건에 따른 적절한 HU 임계 값을 간단하게 얻을 수 있다.
이하, 상술한 위험도 계산 장치(200)가 저전압을 이용한 CT 영상에서 CAC 점수를 계산하기 위한 HU 임계 값을 설정하는 방법에 대하여 설명하도록 한다.
도 5는 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산 방법의 순서도이다.
도 5를 참조하면, 위험도 계산 장치(200)는 사용자에 의해 설정된 관전압 값을 이용하여 혈관에 대한 CT 영상을 획득한다(S110). 위험도 계산 장치(200)는 관전압 값을 영상 촬영 장치(100)로 송신하고, 영상 촬영 장치(100)로부터 지정된 관전압 값에 따른 CT 영상을 수신할 수 있다.
S110 단계 이후, 위험도 계산 장치(200)는 획득된 CT 영상의 히스토그램에서 각각의 HU 값에 따라 픽셀들을 정렬한다(S120). 다시 말해서, 위험도 계산 장치(200)는 동일한 HU 값을 가지는 픽셀들의 개수를 카운트할 수 있다. 예를 들어, 위험도 계산 장치(200)가 정렬한 결과 HU 값이 높아질수록 이에 해당하는 픽셀들의 개수는 적을 수 있다.
S120 단계 이후, 위험도 계산 장치(200)는 정렬된 픽셀들을 HU 값이 높은 순서대로 카운트하고, 기 설정된 임계치까지의 총 픽셀 개수와 동일해지는 픽셀을 결정한다(S130). 구체적으로, 저전압 값에서 촬영한 CT 영상의 경우, 필요 이상의 많은 픽셀 수를 가지고 있는 것이 일반적이다. 그에 따라, 위험도 계산 장치(200)는 최적 표준 값과의 히스토그램 매칭을 위해 픽셀의 개수를 조정할 수 있다. 여기서, 기 설정된 임계치는 120kV 관전압 값으로 촬영된 CT 영상에서 질환을 진단하는 기준이 되는 최소 HU 값인 130HU 일 수 있다.
S130 단계 이후, 위험도 계산 장치(200)는 결정된 픽셀의 개수와 대응되는 HU 값을 새로운 HU 임계 값으로 설정한다(S140). 픽셀의 개수를 통해 설정되는 새로운 HU 임계 값은 CAC 점수를 계산하기 위한 첫번째 기준 값이며, 예를 들어, 새로운 HU 임계 값은 최적 표준인 120kV에서의 최소 HU 임계 값인 130HU 보다 높은 값을 가질 수 있다.
이후, 위험도 계산 장치(200)는 설정한 최소 HU 값을 기준으로 나머지 HU 임계 값들을 순차적으로 새롭게 설정할 수 있으며, 최종적으로 히스토그램을 4개의 구간으로 나눌 수 있는 새로운 HU 임계 값을 얻을 수 있다.
지금까지 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산 방법에 대하여 설명하였다. 본 방법 발명에 따르면, 복잡한 알고리즘이나, kV 광자 에너지 감쇠 비율을 계산 방법을 이용하지 않고, 기준 CT 영상과 촬영된 CT 영상의 히스토그램 픽셀-HU 값 매칭을 통해 CAC 점수 계산을 위한 임계 값을 직관적으로 얻을 수 있다.
도 6은 본 발명의 일 실시 예에 따른 100kV의 관전압 값에서 촬영한 CT 영상의 CAC 점수를 계산하는 방법의 순서도이다.
도 6을 참조하면, 위험도 계산 장치(200)는 영상 촬영 장치(100)로부터 100kV의 관전압 값으로 촬영된 관상동맥에 대한 CT 영상을 획득한다(S210). 이를 위해, 위험도 계산 장치(200)는 영상 촬영 장치(100) 내 구성 요소를 제어하는 정보를 전달할 수 있다.
S210 단계 이후, 위험도 계산 장치(200)는 CT 영상의 히스토그램에서 제1 내지 제4 HU 임계 값을 기초로 관상동맥 칼슘 점수를 계산한다(S220). 구체적으로, 제1 HU 임계 값은, 141HU 이상 145HU 이하, 제2 HU 임계 값은, 218HU 이상 222HU 이하, 제3 HU 임계 값은, 318HU 이상 322HU 이하, 제4 HU 임계 값은, 437HU 이상 441HU 이하의 값을 가질 수 있다.
위험도 계산 장치(200)는 120kV 조건에서 관상동맥 칼슘 점수를 계산하기 위해 사용하는 Agatston 점수 계산 법을 이용하여 새로운 HU 임계 값에 따른 점수를 계산할 수 있다. 예를 들어, 위험도 계산 장치(200)는 143HU이상 220HU이하의 값의 경우 1, 221HU이상 320HU이하의 값을 가지는 경우 2, 321HU이상 439HU이하의 값을 가지는 경우 3, 440HU이상의 값을 가지는 경우 4로 설정하고, “∑(각 선택 부위의 면적x농도인자 x절편두께/3)”와 같은 계산식에 따라 관상동맥 칼슘 점수를 계산할 수 있다.
도 7은 본 발명의 일 실시 예에 따른 저전압 조건에서 촬영된 CT 영상의 CAC 점수를 기준 CAC 점수와 비교한 그래프이다.
한편, 도 7을 설명하기에 앞서, 전압 조건에 따라 서로 다른 CT 영상에서 CAC 점수를 비교하기 위한 환자 그룹은 연령 및 체질량 지수(BMI)의 평균값에 유의미한 차이가 존재하지 않는 환자들을 기준으로 모집되었다.
도 7의 (a)를 참조하면, 동일한 혈관에 대해 120kV 미만의 저전압에서, 120kV의 기준 값에서 촬영한 후, 본 발명의 일 실시 예에 따른 혈관 질환 위험도 계산방법을 이용하여 CAC 점수를 계산했을 때, 그 값의 거의 동일함을 확인할 수 있다.
또한, (b)를 살펴보면, 저전압 조건에서 촬영된 CT 영상에 대해 본 방법 발명을 적용하여 CAC 점수를 계산했을 때, 120kV 기준에서 계산된 CAC 값과의 차이가 신뢰 구간을 크게 벗어나지 않는 것을 확인할 수도 있다.
지금까지 본 발명의 일 실시 예에 따른 저전압을 이용한 CT 영상에서 혈관 질환 위험도를 계산하는 방법에 대하여 설명하였다. 본 발명에 따르면, 100kV의 관전압 조건에서 적용 가능한 HU 임계 값을 이용하여 정확한 CAC 점수를 계산할 수 있으며, 특히, 120kV 보다 낮은 관전압 조건이므로, CT 촬영에 의한 방사선 노출을 최소화할 수 있다.
이상 첨부된 도면을 참조하여 본 발명의 일 실시 예들을 더욱 상세하게 설명하였으나, 본 발명은 반드시 이러한 실시 예로 국한되는 것은 아니고, 본 발명의 기술사상을 벗어나지 않는 범위 내에서 다양하게 변형 실시될 수 있다. 따라서, 본 발명에 개시된 실시 예들은 본 발명의 기술 사상을 한정하기 위한 것이 아니라 설명하기 위한 것이고, 이러한 실시 예에 의하여 본 발명의 기술 사상의 범위가 한정되는 것은 아니다. 그러므로, 이상에서 기술한 실시 예들은 모든 면에서 예시적인 것이며 한정적이 아닌 것으로 이해해야만 한다. 본 발명의 보호 범위는 아래의 청구범위에 의하여 해석되어야 하며, 그와 동등한 범위 내에 있는 모든 기술 사상은 본 발명의 권리범위에 포함되는 것으로 해석되어야 할 것이다.
[이 발명을 지원한 국가연구개발사업]
[과제고유번호] 1345317188
[과제번호] 2018R1D1A1B07048179
[부처명] 교육부
[과제관리(전문)기관명] 한국연구재단
[연구사업명] 개인기초연구(교육부)(R&D)
[연구과제명] 직장암의 항암방사선치료 후 자기공명영상에서 다중 시퀀스 융합영상의 딥러닝 분석을 이용한 잔류 종양 위험도 지도 생성 및 예후예측 시스템 개발
[기여율] 1/2
[과제수행기관명] 연세대학교
[연구기간] 20200301 ~ 20210228
[이 발명을 지원한 국가연구개발사업]
[과제고유번호] 1711115069
[과제번호] 2020R1A2C1003762
[부처명] 과학기술정보통신부
[과제관리(전문)기관명] 한국연구재단
[연구사업명] 개인기초연구(과기정통부)(R&D)
[연구과제명] 폐동맥 고혈압 조기 진단을 위한 4차원 혈류 자기공명영상 생체표지자 개발
[기여율] 1/2
[과제수행기관명] 연세대학교
[연구기간] 20200301 ~ 20210228

Claims (11)

  1. 사용자에 의해 설정된 관전압 값을 이용하여 혈관에 대한 CT 영상을 획득하는 단계;
    획득된 CT 영상의 히스토그램에서 픽셀 별 HU 값(Hounsfield Unit Value)에 따라 픽셀들을 정렬하는 단계;
    정렬된 픽셀들을 상기 HU 값이 높은 순서대로 카운트하고, 기 설정된 임계치까지의 총 픽셀 개수와 동일해지는 픽셀을 결정하는 단계; 및
    결정된 픽셀의 개수와 대응되는 HU 값을 새로운 HU 임계 값으로 설정하는 단계; 를 포함하는 혈관 질환 위험도 계산방법.
  2. 제1항에 있어서,
    상기 사용자에 의해 설정된 관전압 값은,
    기준 CT 영상을 획득하기 위해 설정된 기준 관전압 값보다 낮은 값을 가지는, 혈관 질환 위험도 계산방법.
  3. 제2항에 있어서,
    상기 기준 CT 영상의 히스토그램은,
    기 설정된 4개의 임계 값을 기준으로 적어도 4개의 구간으로 나누어지며,
    상기 새로운 HU 임계 값으로 설정하는 단계는,
    상기 사용자에 의해 설정된 관전압에 따라서 제1 내지 제4 HU 임계 값을 설정하는 단계인, 혈관 질환 위험도 계산방법.
  4. 제3항에 있어서,
    상기 새로운 HU 임계 값으로 설정하는 단계는,
    상기 정렬된 픽셀에서 가장 낮은 HU 값을 상기 제1 HU 임계 값으로 설정하는 단계, 와
    상기 제1 HU 임계 값을 기초로 상기 제2 내지 제4 HU 임계 값을 순차적으로 설정하는 단계, 를 더 포함하는 혈관 질환 위험도 계산방법.
  5. 제2항에 있어서,
    상기 사용자에 의해 설정된 관전압 값은,
    120kV 미만의 관전압 값인, 혈관 질환 위험도 계산방법.
  6. 제1항에 있어서,
    상기 새로운 HU 임계 값으로 설정하는 단계는,
    상기 획득된 CT 영상의 히스토그램에서 상기 새로운 HU 임계 값에 따라 구분된 픽셀들을 표시하는 단계, 를 더 포함하는 혈관 질환 위험도 계산방법.
  7. 관상동맥에 대하여 100kV의 관전압 값으로 컴퓨터 단층촬영을 수행하여, CT 영상을 획득하는 단계; 및
    획득된 CT 영상의 히스토그램에서 제1 내지 제4 HU 임계 값을 기초로 상기 관상동맥 칼슘 점수를 계산하는 단계; 를 포함하는 혈관 질환 위험도 계산방법.
  8. 제7항에 있어서,
    상기 제1 HU 임계 값은, 141HU 이상 145HU 이하,
    상기 제2 HU 임계 값은, 218HU 이상 222HU 이하,
    상기 제3 HU 임계 값은, 318HU 이상 322HU 이하,
    상기 제4 HU 임계 값은, 437HU 이상 441HU 이하의 값을 가지는, 혈관 질환 위험도 계산방법.
  9. 제8항에 있어서,
    상기 제1 내지 제4 HU 임계 값 각각은,
    120kV의 관전압 값으로 촬영된 CT 영상의 히스토그램에 설정된 130HU, 200HU, 300HU 및 400HU 임계 값과 대응되는, 혈관 질환 위험도 계산방법.
  10. 혈관에 대한 CT 영상을 수신하도록 구성된 수신부; 및
    상기 수신부와 동작 가능하게 연결된 프로세서; 를 포함하고,
    상기 프로세서는,
    사용자에 의해 설정된 관전압 값을 이용하여 혈관에 대한 CT 영상을 획득하고, 획득된 CT 영상의 히스토그램에서 픽셀 별 HU 값(Hounsfield Unit Value)에 따라 픽셀들을 정렬하며, 정렬된 픽셀들을 상기 HU 값이 높은 순서대로 카운트하고, 기 설정된 임계치까지의 총 픽셀 개수와 동일해지는 픽셀을 결정하고, 결정된 픽셀의 개수와 대응되는 HU 값을 새로운 HU 임계 값으로 설정하는, 혈관 질환 위험도 계산 장치.
  11. 혈관에 대한 CT 영상을 수신하도록 구성된 수신부; 및
    상기 수신부와 동작 가능하게 연결된 프로세서; 를 포함하고,
    상기 프로세서는,
    관상동맥에 대하여 100kV의 관전압 값으로 컴퓨터 단층촬영을 수행하여, CT 영상을 획득하고, 획득된 CT 영상의 히스토그램에서 제1 내지 제4 HU 임계 값을 기초로 상기 관상동맥 칼슘 점수를 계산하는, 혈관 질환 위험도 계산 장치.
PCT/KR2021/012623 2020-09-16 2021-09-15 저전압을 이용한 ct 영상에서 혈관 질환 위험도 계산 방법 및 장치 WO2022060099A1 (ko)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR10-2020-0118773 2020-09-16
KR1020200118773A KR102418483B1 (ko) 2020-09-16 2020-09-16 저전압을 이용한 ct 영상에서 혈관 질환 위험도 계산 방법 및 장치

Publications (1)

Publication Number Publication Date
WO2022060099A1 true WO2022060099A1 (ko) 2022-03-24

Family

ID=80777280

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/KR2021/012623 WO2022060099A1 (ko) 2020-09-16 2021-09-15 저전압을 이용한 ct 영상에서 혈관 질환 위험도 계산 방법 및 장치

Country Status (2)

Country Link
KR (1) KR102418483B1 (ko)
WO (1) WO2022060099A1 (ko)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100156898A1 (en) * 2008-12-19 2010-06-24 Piedmont Healthcare, Inc. System and method for lesion-specific coronary artery calcium quantification
WO2013150884A1 (ja) * 2012-04-02 2013-10-10 株式会社 日立メディコ X線画像診断装置及びx線発生装置の制御方法
KR20160120707A (ko) * 2016-10-06 2016-10-18 사회복지법인 삼성생명공익재단 의료영상을 이용한 진단 보조 정보 제공 방법 및 그 시스템
KR20200056105A (ko) * 2018-11-14 2020-05-22 울산대학교 산학협력단 관상동맥석회화점수 산정 방법 및 장치

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100156898A1 (en) * 2008-12-19 2010-06-24 Piedmont Healthcare, Inc. System and method for lesion-specific coronary artery calcium quantification
WO2013150884A1 (ja) * 2012-04-02 2013-10-10 株式会社 日立メディコ X線画像診断装置及びx線発生装置の制御方法
KR20160120707A (ko) * 2016-10-06 2016-10-18 사회복지법인 삼성생명공익재단 의료영상을 이용한 진단 보조 정보 제공 방법 및 그 시스템
KR20200056105A (ko) * 2018-11-14 2020-05-22 울산대학교 산학협력단 관상동맥석회화점수 산정 방법 및 장치

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
HOU KUEI‐YUAN, TSUJIOKA KATSUMI, YANG CHING‐CHING: "Optimization of HU threshold for coronary artery calcium scans reconstructed at 0.5‐mm slice thickness using iterative reconstruction", JOURNAL OF APPLIED CLINICAL MEDICAL PHYSICS, vol. 21, no. 2, 1 February 2020 (2020-02-01), US , pages 111 - 120, XP055912262, ISSN: 1526-9914, DOI: 10.1002/acm2.12806 *

Also Published As

Publication number Publication date
KR20220036501A (ko) 2022-03-23
KR102418483B1 (ko) 2022-07-06

Similar Documents

Publication Publication Date Title
JP6436182B2 (ja) 動態画像解析装置
WO2017051944A1 (ko) 의료 영상 판독 과정에서 사용자의 시선 정보를 이용한 판독 효율 증대 방법 및 그 장치
JP6597548B2 (ja) 動態解析システム
JP6772873B2 (ja) 動態解析装置及び動態解析システム
WO2019103440A1 (ko) 피검체의 의료 영상의 판독을 지원하는 방법 및 이를 이용한 장치
JP6418091B2 (ja) 胸部画像表示システム及び画像処理装置
WO2019172621A1 (ko) 질환 예측 방법 및 이를 이용한 질환 예측 디바이스
WO2020071678A2 (ko) 실시간으로 획득되는 위 내시경 이미지를 기반으로 위 병변을 진단하는 내시경 장치 및 방법
WO2018155765A1 (ko) 컴퓨터 단층촬영 영상에서 플라크를 분석하기 위한 방법 및 장치
WO2021015507A2 (ko) 인공지능을 이용한 조영증강 영상 및 비조영증강 영상 변환 방법 및 장치
WO2016137157A1 (ko) 의료 영상 장치 및 의료 영상 처리 방법
WO2023063646A1 (ko) 기계학습 기반 고관절 골절 진단을 위한 골밀도 도출 방법 및 이를 이용한 골밀도 도출 프로그램
US20180018772A1 (en) Dynamic analysis apparatus
WO2020101264A1 (ko) 관상동맥석회화점수 산정 방법 및 장치
WO2010128818A2 (ko) 의료영상 처리 시스템 및 처리 방법
JP6825229B2 (ja) 動態解析システム
WO2019198981A1 (ko) 촬영 이미지 기반의 건강상태 분석 및 정보 제공 방법, 그의 장치 및 그의 기록 매체
WO2020179950A1 (ko) 딥러닝 기반 뇌 질환 진행 예측 방법 및 장치
WO2022060099A1 (ko) 저전압을 이용한 ct 영상에서 혈관 질환 위험도 계산 방법 및 장치
WO2022010075A1 (ko) 의료영상을 기초로 인체 조직을 분석하는 방법 및 그 장치
WO2024053925A1 (ko) 아밀로이드 혈관병증 의심 환자에 대한 아밀로이드 펫 양성 예측 방법 및 장치
WO2021054700A1 (ko) 치아 병변 정보 제공 방법 및 이를 이용한 장치
CN111753831A (zh) 一种图像分析方法、装置、图像采集设备及存储介质
JP5315686B2 (ja) 動態画像診断支援システム
WO2016080736A1 (ko) 의료영상을 이용한 진단 보조 정보 제공 방법 및 그 시스템

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: 21869721

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: 21869721

Country of ref document: EP

Kind code of ref document: A1