KR101484051B1 - Apparatus and method for preprocessing for CAD system using active contour method - Google Patents

Apparatus and method for preprocessing for CAD system using active contour method Download PDF

Info

Publication number
KR101484051B1
KR101484051B1 KR20130071785A KR20130071785A KR101484051B1 KR 101484051 B1 KR101484051 B1 KR 101484051B1 KR 20130071785 A KR20130071785 A KR 20130071785A KR 20130071785 A KR20130071785 A KR 20130071785A KR 101484051 B1 KR101484051 B1 KR 101484051B1
Authority
KR
South Korea
Prior art keywords
contour
image
computer
starting point
preprocessing
Prior art date
Application number
KR20130071785A
Other languages
Korean (ko)
Other versions
KR20150000067A (en
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 중앙대학교 산학협력단
Priority to KR20130071785A priority Critical patent/KR101484051B1/en
Publication of KR20150000067A publication Critical patent/KR20150000067A/en
Application granted granted Critical
Publication of KR101484051B1 publication Critical patent/KR101484051B1/en

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/40Picture signal circuits
    • H04N1/409Edge or detail enhancement; Noise or error suppression
    • H04N1/4092Edge or detail enhancement
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N2201/00Indexing scheme relating to scanning, transmission or reproduction of documents or the like, and to details thereof
    • H04N2201/0001Diagnosis, testing or measuring; Detecting, analysis or monitoring not otherwise provided for
    • H04N2201/0003Method used

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치가 개시된다. 인공물 제거부는 입력된 이미지에 대해 스레쉬홀드를 이용하여 2치 영상(binary image)을 획득하고, 2치 영상을 사전에 설정된 개수로 나누어 분할된 영역에 대해 픽셀 합을 계산하여 포인터의 시작점을 결정하고 시작점을 이용하여 입력 이미지에 대해 불필요한 인공물을 검색하여 제거한다. 윤곽선 추출부는 흉근 경계 정보를 포함하는 윤곽선을 획득하여 새로운 영역 기반의 에너지 함수를 정의하고, 입력된 이미지의 최외곽선을 기초로 점차적으로 내부로 이동하면서 사전에 설정된 기준에 해당하는 점들을 만나면 이동을 정지하고 정지된 지점들을 연결하여 하나 이상의 윤곽선을 생성한다. 흉근 제거부는 윤곽선의 마스크 텀(mask term)에 대한 변수를 이용하여 흉근 이치 영상(binary image of the pectoral muscle)을 추출하여 흉근을 제거한다. 본 발명에 따르면, 종양 조직(tumor tissue)의 정확한 검출을 위해 디지털 유방조영상 이미지(digital mammogram image) 내에 존재하는 서로 다른 종류의 인공물(artifact)을 제거할 수 있다. 이를 통해 의료 영상을 효과적으로 분석할 수 있는 툴을 제공할 수 있다. A preprocessing apparatus for a CAD system using an active contour method is disclosed. The artifact removal unit obtains a binary image using the threshold for the input image, divides the binary image into a predetermined number, calculates a pixel sum for the divided area, and determines the starting point of the pointer And the unnecessary artifacts are retrieved and removed from the input image using the starting point. The contour extractor obtains an outline including the ribbed boundary information, defines an energy function based on the new region, moves gradually inward based on the outline of the input image, and moves when it meets a predetermined reference point One or more outlines are created by connecting the stopped and stopped points. The pectoralis removal unit removes the pectoral muscle by extracting a binary image of the pectoral muscle using the parameters for the mask term of the contour. According to the present invention, different types of artifacts present in a digital mammogram image can be removed for accurate detection of tumor tissue. This can provide a tool for effectively analyzing medical images.

Description

능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치 및 방법{Apparatus and method for preprocessing for CAD system using active contour method}TECHNICAL FIELD The present invention relates to a preprocessing apparatus and method for a CAD system using an active contour method,

본 발명은 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치 및 방법에 관한 것으로, 더욱 상세하게는, 능동적 윤곽선 방법을 이용하여 라벨과 스캐닝된 인공물 및 흉근을 제거하는 유방 조영상 CAD 시스템을 위한 전처리 알고리즘에 관한 것이다. The present invention relates to a preprocessing apparatus and method for a CAD system using an active contour method, and more particularly, to a preprocessing algorithm for a mammographic CAD system that removes labels and scanned artifacts and stomach muscles using an active contour method .

엑스레이, 단층촬영(tomography), MRI에서의 관심 영역(ROI, region-of-interest) 분석은 의료 영상 CAD(computer-aided diagnosis)에서 필수적인 과정이다. 이를 위한 의료 영상의 등질성 지도는 관심 영역 분석에 대한 정보를 제공한다. 따라서 의료 영상에서의 관심 영역을 검출하기 위해 이미지를 분할(segmentation)하는 것에 대해 오랜 시간 동안 많은 연구가 진행되어 왔다. Analysis of x-rays, tomography, and region-of-interest (ROI) in MRI is an essential process in computer-aided diagnosis (CAD). For this purpose, homogeneity maps of medical images provide information on analysis of area of interest. Therefore, a lot of research has been conducted over a long period of time for image segmentation in order to detect a region of interest in a medical image.

상술한 분야의 연구의 하나로, 이미지 내의 관심 물체를 추출하기 위해 액티브 컨투어 모델(active contour model)로 알려진 편미분 방정식(PDE, partial differential equation)에 기반을 둔 방법이 있다. 즉, 액티브 컨투어 모델은 이미지 분할에 사용되는 기술 중 하나로서, 컨투어의 커브(curve)는 타겟 물체 경계를 향해 진행되며 이는 그 경계와 만나게 되면 멈추게 된다. 따라서 의료 영상 뿐만 아니라, 이미지 처리 분야에서 다양한 액티브 컨투어 모델과 업그레이드된 버전의 액티브 컨투어 모델이 적용되어 왔다. 그러나 종래 기법은 등질선 생성시 매번 기준 레벨을 설정하는 재초기화(re-initialization)과정을 거쳐야 하므로, 시간이 오래 걸리고 연산량이 많아지는 문제점이 있었다. As one of the researches in the field mentioned above, there is a method based on a partial differential equation (PDE) known as an active contour model for extracting objects of interest in an image. That is, the active contour model is one of the techniques used for image segmentation, in which the curve of the contour proceeds toward the target object boundary, which stops when it encounters the boundary. Therefore, various active contour models and upgraded versions of active contour models have been applied in image processing as well as medical images. However, the conventional technique requires a re-initialization process to set a reference level every time the generation of a homogeneous line, so that it takes a long time and a large amount of computation.

이와 관련된 선행기술을 구체적으로 살펴보면, 미국공개공보 제2010-0215238호(발명의 명칭 : Method for automatic segmentation of images)에는 이미지를 자동 분할하는 방법에 대해 개시하고 있다. 이는 심장(heart) 이미지에 대해 이미지 데이터를 재구성하여 내외부 표면을 가진 벽을 포함하는 심장 챔버(chamber)를 구별하여 챔버 벽의 내부 표면을 가리키는 내부 등고선을 생성하여 챔버의 의사극지도(pseudo-polar map)을 생성하고 이를 통해 챔버의 외부 표면을 가리키는 외부 등고선을 결정함으로서 내외부 등고선을 이용하여 재구성된 이미지를 구획한다.In detail, a method for automatic segmentation of images is disclosed in U.S. Publication No. 2010-0215238 (entitled Method for automatic segmentation of images). This reconstructs the image data for a heart image to distinguish a heart chamber that includes a wall with an inner and outer surface to create an inner contour pointing to the inner surface of the chamber wall to create a pseudo-polar map to determine an outer contour line that points to the outer surface of the chamber, thereby partitioning the reconstructed image using inner and outer contour lines.

또한, 미국공개공보 제2010-0226542호(발명의 명칭 : Method of bioimage data processing for revealing more meaningful anatomic features of diseased tissues)에는 이미지 데이터 세트에서 이미지 맵을 생성하는 방법을 개시하고 있다. 이는 3D 이미지 데이터 세트 내에서 경계 표면을 구별하고 경계 표면으로부터 참조 표면을 생성하여 3D 이미지 데이터 세트에서 두 번째 표면을 구별하여 참조 표면 위의 지점과 두 번째 표면에 해당하는 지점 간의 거리를 계산하여 이를 토대로 이미지 맵을 생성한다. In addition, US-A-2010-0226542 discloses a method for generating an image map from an image data set in a method of revealing more meaningful anatomical features of diseased tissues. This distinguishes the boundary surface within a 3D image data set and creates a reference surface from the boundary surface to distinguish the second surface from the 3D image data set to calculate the distance between the point on the reference surface and the point on the second surface, Create an image map based on this.

본 발명이 이루고자 하는 기술적 과제는, 종양 조직(tumor tissue)의 정확한 검출을 위해 디지털 유방조영상 이미지(digital mammogram image) 내에 존재하는 서로 다른 종류의 인공물(artifact)을 제거하는 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치 및 방법을 제공하는 데 있다.SUMMARY OF THE INVENTION The present invention is directed to a CAD system using an active contour method for removing different types of artifacts present in a digital mammogram image for accurate detection of tumor tissue. And to provide a preprocessing apparatus and method for the same.

본 발명이 이루고자 하는 다른 기술적 과제는, 종양 조직(tumor tissue)의 정확한 검출을 위해 디지털 유방조영상 이미지(digital mammogram image) 내에 존재하는 서로 다른 종류의 인공물(artifact)을 제거하는 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체를 제공하는 데 있다.Another object of the present invention is to provide a method and apparatus for detecting a tumor tissue by using an active contour method that removes different kinds of artifacts existing in a digital mammogram image, There is provided a computer-readable recording medium recording a program for causing a computer to execute a preprocessing method for a system.

상기의 기술적 과제를 달성하기 위한 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치는, 입력된 이미지에 대해 스레쉬홀드를 이용하여 2치 영상(binary image)을 획득하고, 상기 2치 영상을 사전에 설정된 개수로 나누어 분할된 영역에 대해 픽셀 합을 계산하여 포인터의 시작점을 결정하고 상기 시작점을 이용하여 상기 입력 이미지에 대해 불필요한 인공물을 검색하여 제거하는 인공물 제거부; 흉근 경계 정보를 포함하는 윤곽선을 획득하여 에너지 함수를 정의하고, 상기 입력된 이미지의 최외곽선을 기초로 점차적으로 내부로 이동하면서 사전에 설정된 기준에 해당하는 점들을 만나면 이동을 정지하고 상기 정지된 지점들을 연결하여 하나 이상의 윤곽선을 생성하는 윤곽선 추출부; 및 상기 윤곽선의 마스크 텀(mask term)에 대한 변수를 이용하여 흉근 이치 영상(binary image of the pectoral muscle)을 추출하여 흉근을 제거하는 흉근 제거부;를 구비한다. According to another aspect of the present invention, there is provided a preprocessing apparatus for a CAD system using an active contour method, the apparatus comprising: a binary image obtaining unit that obtains a binary image using a threshold for an input image, An artefact removal unit for dividing an image into a predetermined number and determining a starting point of a pointer by calculating a pixel sum with respect to the divided area, and searching for and removing unnecessary artifacts from the input image using the starting point; And an energy function is defined by acquiring an outline including the rib boundary information, and gradually moves to the inside based on the outermost line of the input image, stops moving when the point corresponding to the predetermined reference is encountered, A contour extractor for concatenating the contour lines to generate at least one contour; And a sternum removing section for extracting a binary image of the pectoral muscle using a parameter for a mask term of the outline to remove the sternum.

상기의 다른 기술적 과제를 달성하기 위한 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법은, 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치에 의해 수행되는 전처리 방법에 있어서, (a) 입력된 이미지에 대해 스레쉬홀드를 이용하여 2치 영상(binary image)을 획득하고, 상기 2치 영상을 사전에 설정된 개수로 나누어 분할된 영역에 대해 픽셀 합을 계산하여 포인터의 시작점을 결정하고 상기 시작점을 이용하여 상기 입력 이미지에 대해 불필요한 인공물을 검색하여 제거하는 인공물 제거 단계; (b) 흉근 경계 정보를 포함하는 윤곽선을 획득하여 에너지 함수를 정의하고, 상기 입력된 이미지의 최외곽선을 기초로 점차적으로 내부로 이동하면서 사전에 설정된 기준에 해당하는 점들을 만나면 이동을 정지하고 상기 정지된 지점들을 연결하여 하나 이상의 윤곽선을 생성하는 윤곽선 추출 단계; 및 (c) 상기 윤곽선의 마스크 텀(mask term)에 대한 변수를 이용하여 흉근 이치 영상(binary image of the pectoral muscle)을 추출하여 흉근을 제거하는 흉근 제거 단계;를 갖는다. According to another aspect of the present invention, there is provided a preprocessing method for a CAD system using an active contour method, the preprocessing method being performed by a preprocessing apparatus for a CAD system using an active contour method, the method comprising the steps of: A binary image is obtained using a threshold, a binary sum is divided into a predetermined number, and a pixel sum is calculated for a divided area to determine a starting point of the pointer, An artifact removal step of searching for and removing an unnecessary artifact to the input image; (b) defining an energy function by acquiring an outline including the ribbed boundary information, moving gradually toward the inside based on the outermost line of the input image, stopping the movement when a point corresponding to a predetermined reference is encountered, An outline extraction step of connecting the stopped points to generate at least one outline; And (c) extracting a binary image of the pectoral muscle using a parameter for the mask term of the contour line to remove the sternal muscle.

상기의 다른 기술적 과제를 달성하기 위한 본 발명에 따른 컴퓨터로 읽을 수 있는 매체는 상기한 방법 중 어느 하나를 컴퓨터에서 실행시키기 위한 프로그램을 기록한다. According to another aspect of the present invention, there is provided a computer readable medium storing a program for causing a computer to execute any one of the above methods.

본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치 및 방법에 의하면, 종양 조직(tumor tissue)의 정확한 검출을 위해 디지털 유방조영상 이미지(digital mammogram image) 내에 존재하는 서로 다른 종류의 인공물(artifact)을 제거할 수 있다. 이를 통해 의료 영상을 효과적으로 분석할 수 있는 툴을 제공할 수 있다. According to the preprocessing apparatus and method for a CAD system using the active contour method according to the present invention, in order to accurately detect a tumor tissue, artifacts of different kinds existing in a digital mammogram image Can be removed. This can provide a tool for effectively analyzing medical images.

도 1은 일반적인 유방조영상 이미지를 도시한 도면,
도 2는 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치의 구성을 도시한 블록도,
도 3은 본 발명의 일실시예에 따라 인공물을 제거하는 과정을 도시한 결과 영상,
도 4는 능동적 윤곽선 방법을 이용하여 다상 분할 과정을 도시한 도면,
도 5는 본 발명의 일실시예에 따라 능동적 윤곽선 분할 과정을 도시한 도면,
도 6은 본 발명의 일실시예에 따라 흉근을 제거하는 과정을 도시한 결과 영상,
도 7은 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법의 과정을 도시한 흐름도,
도 8은 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법과 다른 기존 2가지 방법을 적용한 결과 영상을 비교한 도면, 그리고,
도 9 내지 도 10은 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법을 이용한 분할 결과 영상을 도시한 도면이다.
1 shows a general mammogram image,
FIG. 2 is a block diagram showing the construction of a preprocessing apparatus for a CAD system using the active contour method according to the present invention.
FIG. 3 is a flowchart illustrating a process of removing artifacts according to an embodiment of the present invention.
FIG. 4 is a view showing a polyphase division process using an active contour method;
FIG. 5 illustrates an active contour splitting process according to an embodiment of the present invention.
FIG. 6 is a flowchart illustrating a process of removing a stomach muscle according to an embodiment of the present invention.
FIG. 7 is a flowchart illustrating a preprocessing method for a CAD system using the active contour method according to the present invention.
FIG. 8 is a diagram illustrating a comparison between a pre-processing method for a CAD system using the active contour method according to the present invention,
9 to 10 are diagrams showing divided result images using a preprocessing method for a CAD system using the active contour method according to the present invention.

이하에서 첨부의 도면들을 참조하여 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치 및 방법의 바람직한 실시예에 대해 상세하게 설명한다.DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Hereinafter, preferred embodiments of a pre-processing apparatus and method for a CAD system using an active contour method according to the present invention will be described in detail with reference to the accompanying drawings.

도 1은 인공물(artifact)과 흉근(pectoral muscle)을 포함하는 유방조영상 이미지(mammogram image)를 도시한 도면이다. 도 1을 참조하면, 디지털 유방 조영상 이미지는 서로 다른 종류의 인공물을 가질 수 있다. 이때, 인공물은 고강도 라벨(high intensity label), 저강도 라벨(low intensity label), 테이핑되거나 스캐닝된 인공물(taping or scanning artifacts), 그리고 흉근(pectoral muscle)을 포함할 수 있다. 유방조영상 이미지 내에 상술한 인공물을 포함되어 있으면 검출 과정에 문제가 발생할 수 있고 관심 영역(ROI : Region Of Interest) 검출의 정확성에 영향을 끼칠 수 있다. 따라서 이러한 문제가 되는 인공물을 제거할 수 있는 전처리 알고리즘이 중요하다. BRIEF DESCRIPTION OF THE DRAWINGS Figure 1 is a diagram showing a mammogram image including artifacts and pectoral muscles. Referring to FIG. 1, a digital mammogram image may have artifacts of different kinds. Artifacts may include high intensity labels, low intensity labels, taping or scanning artifacts, and pectoral muscles. If the artifacts are included in the mammogram image, a problem may occur in the detection process and the accuracy of ROI (Region Of Interest) detection may be affected. Therefore, a preprocessing algorithm that can remove these artifacts is important.

도 2는 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치(100)의 구성을 도시한 블록도이다. FIG. 2 is a block diagram showing a configuration of a preprocessing apparatus 100 for a CAD system using an active contour method according to the present invention.

도 2를 참조하면, 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치(100)는 인공물 제거부(200), 윤곽선 추출부(300) 및 흉근 제거부(400)를 포함한다. Referring to FIG. 2, a preprocessing apparatus 100 for a CAD system using an active contour method according to the present invention includes an artifact removing unit 200, an outline extracting unit 300, and a suture removing unit 400.

인공물 제거부(200)는 입력된 이미지에 대해 스레쉬홀드를 이용하여 2치 영상(binary image)을 획득하고, 2치 영상을 사전에 설정된 개수로 나누어 분할된 영역에 대해 픽셀 합을 계산하여 포인터의 시작점을 결정하고 시작점을 이용하여 입력 이미지에 대해 불필요한 인공물을 검색하여 제거한다. The artificial material removing unit 200 obtains a binary image using the threshold for the input image, divides the binary image into a preset number, calculates a pixel sum with respect to the divided area, And the unnecessary artifacts are retrieved and removed from the input image using the starting point.

즉, 첫 단계로, 라벨과 스캐닝된 인공물은 스레쉬홀드(thresholding)를 이용하여 제거된다. 먼저, 유방조영상 이미지를 가지고 와서, 도 3의 (b)에 도시된 바와 같이 2치영상(binary image)을 얻기 위해 스레쉬홀드 값을

Figure 112013055704111-pat00001
를 이용하여 스레쉬홀드를 적용한다. 그리고나서 도 3의 (c)에 도시된 바와 같이 이미지를 수직으로 절반으로 나누어 2개의 그룹을 얻는다. 그 후에 각 그룹에 대해 픽셀 합(pixel sum)을 계산하고, 포인터의 시작점(starting point)을 결정하기 위해 어떤 것이 더 큰 값을 가지는지 비교한다. That is, in the first step, the label and the scanned artifact are removed using thresholding. First, a mammogram image is fetched, and a threshold value is obtained to obtain a binary image as shown in FIG. 3 (b)
Figure 112013055704111-pat00001
To apply the threshold. Then, two groups are obtained by dividing the image vertically by half as shown in FIG. 3 (c). Then calculate the pixel sum for each group and compare what has a larger value to determine the starting point of the pointer.

다음 수학식 1을 통해, 포인터는 유방조영상 이미지 내의 왼쪽 코너에 있거나, 오른쪽 코너에 있을 수 있다. Through the following equation (1), the pointer can be located in the left corner of the mammogram, or in the right corner.

Figure 112013055704111-pat00002
Figure 112013055704111-pat00002

여기서,

Figure 112013055704111-pat00003
은 행(rows)의 최대치 수,
Figure 112013055704111-pat00004
는 열(columns)의 최대치 수, 그리고,
Figure 112013055704111-pat00005
는 이치영상을 나타낸다. here,
Figure 112013055704111-pat00003
Is the maximum number of rows,
Figure 112013055704111-pat00004
Is the maximum number of columns,
Figure 112013055704111-pat00005
Represents a binary image.

도 3의 (e)에 도시된 바와 같이, 시작점을 찾은 후엔, 이미지의 모든 행(row)에 대해 돌아다니면서 이미지 내의 불필요한 데이터(예를 들어, 라벨이나 스캐닝된 인공물 등)를 제거한다. 따라서 불필요한 데이터를 제거한 후, 도 3의 (f)에 도시된 바와 같이, 바람직한 결과를 얻기 위해 원본 유방조영상 이미지를 크게 증가시킨다. After finding the starting point, as shown in FIG. 3 (e), unnecessary data (for example, a label or a scanned artifact, etc.) in the image is removed by moving around all the rows of the image. Therefore, after removing unnecessary data, the original mammogram image is greatly increased to obtain a desired result as shown in FIG. 3 (f).

윤곽선 추출부(300)는 흉근 경계 정보를 포함하는 윤곽선을 획득하여 새로운 영역 기반의 에너지 함수를 정의하고, 상기 입력된 이미지의 최외곽선을 기초로 점차적으로 내부로 이동하면서 사전에 설정된 기준에 해당하는 점들을 만나면 이동을 정지하고 상기 정지된 지점들을 연결하여 하나 이상의 윤곽선을 생성한다. The contour extracting unit 300 acquires an outline including the ribbone boundary information to define an energy function based on the new region, and gradually moves to the inside based on the outermost line of the input image, Upon encountering the points, the movement is stopped and the stopped points are connected to create one or more contour lines.

즉, 두 번째 단계로, 능동적 윤곽선 방법을 이용하여 흉근 경계 정보를 포함하는 윤곽선(contour)을 획득한다. Chan&Vese(CV 모델)는 Mumford-Shah model을 기본으로 능동적 윤곽선(active contour) 방법을 제안하였다. 입력 이미지를

Figure 112013055704111-pat00006
로 하고, 폐곡선(closed curve)을
Figure 112013055704111-pat00007
로 하면, 다음 수학식 2와 같은 CV 모델을 이용한 새로운 영역 기반 에너지 함수를 정의할 수 있다. In other words, as a second step, the active contour method is used to obtain a contour that includes the information of the pectoral muscle boundary. Chan & Vese (CV model) proposed an active contour method based on the Mumford-Shah model. Input image
Figure 112013055704111-pat00006
, And the closed curve
Figure 112013055704111-pat00007
, It is possible to define a new region-based energy function using the CV model as shown in Equation (2).

Figure 112013055704111-pat00008
Figure 112013055704111-pat00008

Figure 112013055704111-pat00009
Figure 112013055704111-pat00009

여기서,

Figure 112013055704111-pat00010
,
Figure 112013055704111-pat00011
,
Figure 112013055704111-pat00012
, 그리고,
Figure 112013055704111-pat00013
은 고정된 변수들이다.
Figure 112013055704111-pat00014
은 Heaviside 함수이고,
Figure 112013055704111-pat00015
는 새로 소개된 마스크 텀(mask term)이며, 이는 후술할 수학식 3 및 수학식 4에서 윤곽선이 내부로 이동하는 데 도움을 준다.
Figure 112013055704111-pat00016
는 윤곽선의 정지점(stopping point)을 결정하기 위해 2개의 연속되는 윤곽선 사이의 차이를 계산하는 데 사용된다. 이는 CV 2상(two-phase) 능동적 윤곽선 모델을 다상(multiphase) 능동적 윤곽선 모델로 전환한다. here,
Figure 112013055704111-pat00010
,
Figure 112013055704111-pat00011
,
Figure 112013055704111-pat00012
, And,
Figure 112013055704111-pat00013
Are fixed variables.
Figure 112013055704111-pat00014
Is a Heaviside function,
Figure 112013055704111-pat00015
Is a newly introduced mask term, which helps the outline to move inward in equations (3) and (4) to be described later.
Figure 112013055704111-pat00016
Is used to calculate the difference between two consecutive contours to determine the stopping point of the contour. This translates the CV two-phase active contour model into a multiphase active contour model.

Figure 112013055704111-pat00017
Figure 112013055704111-pat00017

Figure 112013055704111-pat00018
Figure 112013055704111-pat00018

레벨 셋 함수(level set function)와 그것의 발전을 수반하여 더욱 정확한 계산을 위해 부호화 거리 함수(signed distance function)로부터 편차(deviation)를 페널라이징(penalizing)하는 레벨 셋 함수를 규칙화할 필요가 있다. 페널라이징 텀의 에너지 함수(

Figure 112013055704111-pat00019
)와 토탈 에너지 수식(
Figure 112013055704111-pat00020
)은 아래 수학식 5 및 수학식 6에 기재되어 있다. There is a need to regularize a level set function that penalizes deviations from a signed distance function for more accurate calculations accompanied by a level set function and its evolution. The energy function of the penalizing term (
Figure 112013055704111-pat00019
) And total energy formula (
Figure 112013055704111-pat00020
) Are described in the following equations (5) and (6).

Figure 112013055704111-pat00021
Figure 112013055704111-pat00021

Figure 112013055704111-pat00022
Figure 112013055704111-pat00022

다음 수학식 8 및 수학식 9에 기재된

Figure 112013055704111-pat00023
Figure 112013055704111-pat00024
는 윤곽선(
Figure 112013055704111-pat00025
)의 내부(inside)와 외부(outside) 각각의 이미지 인텐서티(intensity)를 계산해낸다.
Figure 112013055704111-pat00026
에 의해 초기화된 레벨 셋 함수
Figure 112013055704111-pat00027
는 다음 수학식 7을 참고한다. The following equations (8) and
Figure 112013055704111-pat00023
and
Figure 112013055704111-pat00024
Is the outline (
Figure 112013055704111-pat00025
And calculates the image intensities of the inside and outside of the image.
Figure 112013055704111-pat00026
Level set function initialized by
Figure 112013055704111-pat00027
(7).

Figure 112013055704111-pat00028
Figure 112013055704111-pat00028

최속 강하법(steepest descent method)을 이용하여 수학식 6에 기재된 토탈 에너지 함수를 최소화하는 것과 제로 레벨 셋을 가지고 윤곽선

Figure 112013055704111-pat00029
를 나타내면서, 다음 수학식 8 내지 수학식 10에 각각 기재된 변화 공식화(variational formulation)를 얻을 수 있다. The steepest descent method is used to minimize the total energy function described in equation (6)
Figure 112013055704111-pat00029
The variational formulation described in each of the following equations (8) to (10) can be obtained.

Figure 112013055704111-pat00030
Figure 112013055704111-pat00030

Figure 112013055704111-pat00031
Figure 112013055704111-pat00031

Figure 112013055704111-pat00032
Figure 112013055704111-pat00032

여기서,

Figure 112013055704111-pat00033
은 라플라시안 연산자(laplacian operator)이다. 따라서 함수
Figure 112013055704111-pat00034
는 오일러-라그랑지 방정식
Figure 112013055704111-pat00035
을 만족시키는 함수를 최소화시킨다. 수학식 10에 기재된 마지막 텀은 구배법(gradient descent method)를 이용하여 수학식 5를 해결하면서 획득된 규칙화 텀이다. 데이터 피팅 텀
Figure 112013055704111-pat00036
는 커브 전개(curve evolution)에 중요한 역할을 하며,
Figure 112013055704111-pat00037
Figure 112013055704111-pat00038
는 최초 텀와 두 번째 텀 사이에서 트레이드오프(tradeoff)를 지배한다. 대부분의 경우,
Figure 112013055704111-pat00039
,
Figure 112013055704111-pat00040
,
Figure 112013055704111-pat00041
는 스케일링 변수(scaling parameter)로 설정하고,
Figure 112013055704111-pat00042
는 univariate Dirac 함수를 의미한다.
Figure 112013055704111-pat00043
함수를 최소화하는 고전적인 반복 과정은 인공 시간(artificial time)
Figure 112013055704111-pat00044
를 이용한 다음 수학식 11의 그라디언트 흐름(gradient flow)일 수 있다. here,
Figure 112013055704111-pat00033
Is a laplacian operator. Therefore,
Figure 112013055704111-pat00034
Euler-Lagrange equation
Figure 112013055704111-pat00035
Is minimized. The end term defined in Equation (10) is a rule term obtained by solving Equation (5) using a gradient descent method. Data fittings Term
Figure 112013055704111-pat00036
Plays an important role in curve evolution,
Figure 112013055704111-pat00037
and
Figure 112013055704111-pat00038
Dominates the tradeoff between the first term and the second term. In most cases,
Figure 112013055704111-pat00039
,
Figure 112013055704111-pat00040
,
Figure 112013055704111-pat00041
Is set as a scaling parameter,
Figure 112013055704111-pat00042
Means the univariate Dirac function.
Figure 112013055704111-pat00043
The classic iterative process of minimizing the function is artificial time,
Figure 112013055704111-pat00044
And can be a gradient flow of Equation (11).

Figure 112013055704111-pat00045
Figure 112013055704111-pat00045

이러한 알고리즘은 객체를 분할하는 데 세컨드 영역 분할 접근법을 이용한다. 이는 이미지를 2개의 서브 영역으로 나누고, 내부 영역이 비거나(empty) 정지 조건(stopping condition)이 충족될 때까지 2개의 서브 영역 내부를 계속하여 2개의 서브 영역으로 나눈다. 최종 결과를 얻기 위한 반복 횟수를 n으로 설정하면, n+1은 분할된 영역의 전체 수를 나타내게 된다. 도 4는 능동적 윤곽선 방법(active contour method)을 이용하여 다상(multiphase) 분할되는 과정을 나타내고 있다. 도 4의 (a)에 도시된 윤곽선 지도(contour map)는 베이스 윤곽선으로부터 최종 윤곽선을 나타내고 있다. 분할 과정을 통해 생성된 윤곽선 지도로부터 획득된 내부 영역은 도 4의 (b)에 도시된 바와 같이 실선(solid line)에 의해 생성되는 유색 노드(colored node)로 나타나고 있다. 내부 영역과 외부 영역의 관계는 다음 수학식 12와 같이 나타낼 수 있다. These algorithms use a second domain segmentation approach to partition objects. This divides the image into two sub-regions and continues to divide the two sub-regions into two sub-regions until an empty stopping condition is met. When the number of repetitions for obtaining the final result is set to n, n + 1 represents the total number of divided areas. FIG. 4 shows a process of multiphase segmentation using an active contour method. The contour map shown in Fig. 4 (a) shows the final contour line from the base contour line. The inner region obtained from the contour map generated through the division process is represented as a colored node generated by a solid line as shown in FIG. 4 (b). The relationship between the inner and outer regions can be expressed as Equation (12).

Figure 112013055704111-pat00046
Figure 112013055704111-pat00046

본 발명에 따른 윤곽선 추출부(300)는 외부 영역을 버리고, 다음 단계에서 다시 한번 2개의 부분으로 나누기 위해 현 내부 윤곽선을 이용한다. 정지 알고리즘을 이용하여 임의의 지점에서 윤곽선을 정지시키기 위해 2개의 연속적인 윤곽선 내에서 각 픽셀의 유사성을 다음과 같이 체크한다. The contour extractor 300 according to the present invention discards the outer region and uses the inner contour to divide the contour into two portions in the next step. The similarity of each pixel within two consecutive contours is checked as follows to stop the contour at any point using the stop algorithm.

Figure 112013055704111-pat00047
Figure 112013055704111-pat00047

여기서,

Figure 112013055704111-pat00048
이고,
Figure 112013055704111-pat00049
는 최종 계산된 윤곽선의 마스크 텀(mask term)이고,
Figure 112013055704111-pat00050
는 현재 계산된 윤곽선의 마스크 텀이고,
Figure 112013055704111-pat00051
은 원본 유방조영상 이미지 내 행(row)의 최대 수를 나타내고, 그리고,
Figure 112013055704111-pat00052
는 열(column)의 최대 수를 나타낸다. 도 5는 계산된 윤곽선을 나타내고, 이는 흉근 경계를 포함하고 있다. here,
Figure 112013055704111-pat00048
ego,
Figure 112013055704111-pat00049
Is the mask term of the last calculated contour,
Figure 112013055704111-pat00050
Is the mask term of the currently calculated contour,
Figure 112013055704111-pat00051
Represents the maximum number of rows in the original mammogram image,
Figure 112013055704111-pat00052
Represents the maximum number of columns. Figure 5 shows the computed contours, which include the bust breech border.

흉근 제거부(400)는 윤곽선의 마스크 텀(mask term)에 대한 변수를 이용하여 흉근 이치 영상(binary image of the pectoral muscle)을 추출하여 흉근을 제거한다. The sternal descender 400 removes the pectoral muscle by extracting a binary image of the pectoral muscle using the parameters for the mask term of the contour.

즉, 세 번째 단계로, 본 발명은 상술한

Figure 112013055704111-pat00053
값을 이용하여 흉근을 제거한다.
Figure 112013055704111-pat00054
는 마지막에 계산된 윤곽선의 이치영상(binary image)을 제공한다. 상술한 이치 영상을 이용하여, 흉근 이치 영상(binary image of the pectoral muscle)을 추출한다. 먼저, 도 3의 (c)에 사용된 동일한 룰(rule)을 이용하여 포인터를 이동하여 시작점(starting point)을 찾는다. 시작점을 찾은 후, 흉근은 항상 시작점과 연결되어 있다는 정보를 고려하면서 객체가 시작점과 연결되도록 하여 남아있는 영역을 제거한다. 도 6은 유방조영상으로부터 흉근이 제거되는 과정을 도시하고 있다. That is, as a third step,
Figure 112013055704111-pat00053
Values are used to remove the pectoral muscle.
Figure 112013055704111-pat00054
Provides a binary image of the last calculated contour. A binary image of the pectoral muscle is extracted using the above-mentioned binary image. First, a pointer is moved by using the same rule used in (c) of FIG. 3 to find a starting point. After finding the starting point, we take the information that the pectoral muscle is always connected to the starting point, and connect the object to the starting point to remove the remaining area. FIG. 6 shows a process in which the pectoral muscle is removed from the mammogram.

도 7은 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법의 과정을 도시한 흐름도이다. 도 7을 참조하면, 먼저 유방조영상 이미지가 입력되면(S100) 사용자가 검출하고자 하는 관심 영역(ROI : Region Of Interest) 검출(S300)의 효율성을 증대시키기 위해 전처리 과정(S200)을 거치게 된다. 구체적으로, 먼저 입력된 유방조영상 이미지에 대해 스레쉬홀드를 이용하여 2치영상(binary image)를 획득하고, 2치 영상을 사전에 설정된 개수로 나누어 분할된 영역 각각에 대해 픽셀 값을 계산하여 포인터의 시작점을 결정한다. 결정된 포인터의 시작점을 이용하여 입력된 유방조영상 이미지에 포함된 인공물(artifact)를 제거한다(S210). 그리고 나서 능동적 윤곽선 방법을 이용하여 흉근 경계 정보를 포함하는 윤곽선을 추출하고(S220) 추출된 윤곽선의 마스크 텀 정보를 이용하여 최종적으로 흉근을 제거하게 된다(S230). 7 is a flowchart illustrating a preprocessing method for a CAD system using the active contour method according to the present invention. Referring to FIG. 7, when a mammogram image is input (S100), a preprocessing process (S200) is performed to increase the efficiency of detecting a region of interest (ROI) (S300) the user wants to detect. Specifically, a binary image is first obtained using a threshold for the input mammography image, a pixel value is calculated for each of the divided regions by dividing the binary image into a predetermined number, As shown in FIG. The artifact included in the inputted mammography image is removed using the determined starting point of the pointer (S210). Then, an active contour line method is used to extract contour lines including the contour line information (S220), and ultrasound contours are removed using the mask contour information of the extracted contour lines (S230).

이하에서는, 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법을 다른 기존 방법과 비교하여 더욱 정확한 결과 영상이 어떠한 것인지 살펴보도록 한다. 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법은 앞서 설명한 바와 같이 전처리 단계에서 원치 않은 인공물을 제거한다. Hereinafter, the preprocessing method for the CAD system using the active contour method according to the present invention will be compared with other existing methods, and a more accurate result image will be examined. The pre-processing method for the CAD system using the active contour method according to the present invention removes unwanted artifacts in the pre-processing step as described above.

도 8은 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법과 다른 기존 2가지 방법을 적용한 결과 영상을 비교한 것이다. 첫 번째 행(row 1)은 각각 mdb003, mdb005, mdb028, mdb044, mdb095, mdb123로 넘버링되어 mini-MIAS 데이터베이스로부터 얻은 원본 유방조영상 이미지를 나타낸다. 두 번째 행(row 2)은 흉근 영역의 지상검증자료(ground truth)를 나타낸다. 세 번째 행(row 3) 및 네 번째 행(row 4)은 각각 기존 2가지 방법을 각각 적용하여 흉근을 분할한 결과 영상을 나타낸다. 다섯 번째 행(row5)은 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법을 적용한 분할 결과를 나타낸다. 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법에 사용된 윤곽선 정지 값은 mdb003, mdb028, mdb123은 0.999, mdb005, mdb044, mdb095는 0.9995 이다. 다른 변수는 모두 동일하게 설정했다. FIG. 8 is a graph comparing a pre-processing method for a CAD system using the active contour method according to the present invention and a result image obtained by applying the other two conventional methods. The first row (row 1) represents the original mammograms obtained from the mini-MIAS database, numbered as mdb003, mdb005, mdb028, mdb044, mdb095, and mdb123, respectively. The second row (row 2) represents the ground truth of the pectoral region. The third row (row 3) and the fourth row (row 4) show images obtained by dividing the pectoral muscle by applying the two existing methods, respectively. The fifth row (row 5) shows the result of the division using the preprocessing method for the CAD system using the active contour method according to the present invention. The contour stop values used in the preprocessing method for the CAD system using the active contour method according to the present invention are mdb003, mdb028, mdb123 are 0.999, mdb005, mdb044, and mdb095 are 0.9995. All other variables were set the same.

유방조영상 mdb003, mdb005, mdb028, mdb044, mdb095의 흉근 분할 결과를 비교하면, 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법을 이용한 분할 결과가 다른 기존 2가지 방법을 적용할 때와 비교하여 더욱 흉근 지상검증자료와 가깝다는 것을 알 수 있다. 기존 방법은 흉근이 선형으로 거의 직선 형태로 분할되어 굴곡진(curved) 흉근 영역을 놓치게 되고, 또한, 흉근의 러프한 보더(rough border)를 찾기 위해 Otsu 스레쉬홀드 기법과 수학 형태학 처리를 이용하게 된다. 그러나 그 어떠한 기존 방법도 본 발명처럼 정확하지 않다. 도 8의 세 번째 행(row 3)을 보면, 분할 시 흉근의 곡선 경계는 없음을 알 수 있다. 따라서 결과 영상은 선형 직선 분할의 false positive, true negative 영역을 포함하기 때문에 정확성의 문제를 가지게 된다. mdb003의 분할 결과 영상에서, 다섯 번째 열(row 5)에 도시된 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법은 잘 분할함에 비해, 네 번째 열(row 4)을 살펴보면 흉근의 임의의 부분은 놓치고 있음을 알 수 있다. mdb005의 분할 결과 영상 중 네 번째 열(row 4)과 다섯 번째 열(row 5)에 도시된 결과 영상은 거의 비슷함을 알 수 있다. mdb 028의 분할 결과 영상을 살펴보면, 흉근 상의 조그만 혹(hump)은 다른 방법에서는 놓치고 있음에 반해,본 발명에서는 잘 분할되어 있음을 알 수 있다. mdb 044 분할 결과 영상 중 네 번째 열(row 4)과 다섯 번째 열(row 5)에 도시된 결과 영상은 거의 비슷함을 알 수 있다. mdb095 분할 결과 영상에서, 흉근 경계의 조그만 디테일들은 다른 방법이 적용된 경우 거의 놓치고 있음에 반해, 본 발명에서는 잘 분할되고 있음을 알 수 있다. Comparing the results of sternotomy segmentation of the mammograms mdb003, mdb005, mdb028, mdb044, mdb095, the results of the segmentation using the preprocessing method for the CAD system using the active contour method according to the present invention are compared with those of the conventional two methods And it is closer to the pesterial ground test data. In the conventional method, the pectoral muscle is linearly divided into a nearly linear shape to lose the curved pectoral region, and the Otsu threshold technique and the mathematical morphological processing are used to find a rough border of the pectoral muscle do. However, any existing method is not as accurate as the present invention. In the third row (row 3) of FIG. 8, it can be seen that there is no curve boundary of the sternum at the time of division. Therefore, the resultant image has a problem of accuracy because it includes false positive and true negative regions of linear straight line segmentation. In the segmentation result image of mdb003, the preprocessing method for the CAD system using the active contour method according to the present invention shown in the fifth column (row 5) is well divided, whereas in the fourth row (row 4) Is missing. It can be seen that the resultant images shown in the fourth column (row 4) and the fifth column (row 5) of the divided image of mdb005 are almost similar. Looking at the segmentation result image of mdb 028, it can be seen that a small hump on the sternum is missed in other methods, but is well divided in the present invention. The result images of the fourth row (row 4) and the fifth row (row 5) of the mdb 044 division result image are almost similar. In the mdb095 segmentation result image, the small details of the pectoral muscle boundaries are almost missed when other methods are applied, but are well divided in the present invention.

본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법은 mini-MIAS 데이터베이스로부터 획득된 유방조영상 이미지에 적용된다. 모든 이미지의 인텐서티 범위(range of intensity)는 0에서 255이며, 각 이미지는 1024x1024 픽셀을 가진다. 이때,

Figure 112013055704111-pat00055
,
Figure 112013055704111-pat00056
,
Figure 112013055704111-pat00057
,
Figure 112013055704111-pat00058
,
Figure 112013055704111-pat00059
,
Figure 112013055704111-pat00060
,
Figure 112013055704111-pat00061
, 그리고,
Figure 112013055704111-pat00062
을 이용한다. 도 9 내지 도 10은 본 발명에 따른 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 방법을 이용한 분할 결과 영상을 나타낸다. 상술한 바와 같이, mini-MIAS 데이터베이스로부터 획득된 이미지를 이용하여 본 발명을 적용하며, 정확성과 효율성에 영향을 끼치는 2가지 문제점이 있음을 전제한다. 본 발명의 정확성은 흉근의 인텐서티와 유방 영역의 나머지(rest of the breast region)는 대부분 유사하거나, 흉근의 인텐서티는 남아있는 유방 영역보다 더 낮은 경우 영향을 미친다. 도 9는 false positive와 false negative의 결과 영상을 보여주고 있다. 본 발명의 효율성은 윤곽선이 높은 인텐서티 영역 내에서 매우 천천히 움직임에 따라 높은 인텐서티 스무스(smooth) 영역내에서 유방 조영상에 영향을 미친다. 영상의 나머지 영역에서는, 본 발명은 도 10에 도시된 바와 같이 흉근의 정확한 분할 결과를 제공한다. A preprocessing method for a CAD system using the active contour method according to the present invention is applied to a mammogram image obtained from a mini-MIAS database. The intensities range of all images is 0 to 255, and each image has 1024 x 1024 pixels. At this time,
Figure 112013055704111-pat00055
,
Figure 112013055704111-pat00056
,
Figure 112013055704111-pat00057
,
Figure 112013055704111-pat00058
,
Figure 112013055704111-pat00059
,
Figure 112013055704111-pat00060
,
Figure 112013055704111-pat00061
, And,
Figure 112013055704111-pat00062
. 9 to 10 show divided result images using a preprocessing method for a CAD system using the active contour method according to the present invention. As described above, it is assumed that the present invention is applied to images obtained from the mini-MIAS database, and that there are two problems affecting accuracy and efficiency. The accuracy of the present invention is influenced by the intensities of the pectoral muscle and the rest of the breast region are mostly similar, or the intensities of the pectoral muscle are lower than the remaining breast regions. Figure 9 shows the result image of false positive and false negative. The efficiency of the present invention affects the mammogram in a high intensitized smooth region with very slow motion within the contour high intensities region. In the remaining region of the image, the present invention provides an accurate segmentation result of the sternum as shown in Fig.

표 1 및 표 2에 기재된 본 발명에 따른 방법의 정확도는 다음 수학식 14를 이용하여 계산된 것이다. The accuracy of the method according to the present invention shown in Tables 1 and 2 is calculated using the following Equation (14).

Figure 112013055704111-pat00063
Figure 112013055704111-pat00063

여기서,

Figure 112013055704111-pat00064
는 흉근 영역의 지상검증자료이고, 그리고,
Figure 112013055704111-pat00065
는 본 발명에 의해 계산된 흉근 영역을 나타낸다. here,
Figure 112013055704111-pat00064
Is the ground validation data of the pectoral region,
Figure 112013055704111-pat00065
Represents the pectoral region calculated by the present invention.

유방조영상 넘버Mammography 정확도 퍼센트(%)Accuracy percent (%) 018018 86.9986.99 030030 84.7984.79 134134 71.9471.94 253253 64.6864.68

유방조영상 넘버Mammography 정확도 퍼센트(%)Accuracy percent (%) 006006 97.7097.70 009009 98.2598.25 132132 97.3597.35 249249 98.0798.07

표 1 및 표 2로부터 계산된 정확도의 평균 확률은 각각 77.10%, 97.84%이다. 본 발명에 따른 방법의 정확도는 윤곽선의 정지 값의 선택에 영향을 받는다. 고 인텐서티 영역을 가진 이미지, 그리고 남아있는 유방 영역의 인텐서티에 가까운 흉근의 인텐서티를 가진 이미지에서는 정지 값이 1에 가깝게 된다. 그러나 흉근 영역과 남아있는 유방 영역 사이의 상당한 인텐서티 차이를 가진 이미지는 정지 값이 0.99에 가깝게 된다. The average probabilities of the accuracies calculated from Table 1 and Table 2 are 77.10% and 97.84%, respectively. The accuracy of the method according to the invention is influenced by the selection of the outline of the outline. In images with intensity intensities and images with intensities of the sternum close to the intensities of the remaining breast regions, the stop value is close to unity. However, images with significant intensities differences between the pectoral region and the remaining breast region are close to a stop value of 0.99.

본 발명은 또한 컴퓨터로 읽을 수 있는 기록매체에 컴퓨터가 읽을 수 있는 코드로서 구현하는 것이 가능하다. 컴퓨터가 읽을 수 있는 기록매체는 컴퓨터 시스템에 의하여 읽혀질 수 있는 데이터가 저장되는 모든 종류의 기록장치를 포함한다. 컴퓨터가 읽을 수 있는 기록매체의 예로는 ROM, RAM, CD-ROM, 자기 테이프, 플로피 디스크, 광 데이터 저장장치 등이 있으며, 또한 캐리어 웨이브(인터넷을 통한 전송)의 형태로 구현되는 것도 포함한다. 또한, 컴퓨터가 읽을 수 있는 기록매체는 유무선 통신망으로 연결된 컴퓨터 시스템에 분산되어 분산방식으로 컴퓨터가 읽을 수 있는 코드가 저장되고 실행될 수 있다. The present invention can also be embodied as computer-readable codes on a computer-readable recording medium. A computer-readable recording medium includes all kinds of recording apparatuses in which data that can be read by a computer system is stored. Examples of the computer-readable recording medium include a ROM, a RAM, a CD-ROM, a magnetic tape, a floppy disk, an optical data storage device, and the like, and a carrier wave (transmission via the Internet). In addition, the computer-readable recording medium may be distributed to a computer system connected to a wired / wireless communication network, and a computer-readable code may be stored and executed in a distributed manner.

이상에서 본 발명의 바람직한 실시예에 대해 도시하고 설명하였으나, 본 발명은 상술한 특정의 바람직한 실시예에 한정되지 아니하며, 청구범위에서 청구하는 본 발명의 요지를 벗어남이 없이 당해 발명이 속하는 기술분야에서 통상의 지식을 가진 자라면 누구든지 다양한 변형 실시가 가능한 것은 물론이고, 그와 같은 변경은 청구범위 기재의 범위 내에 있게 된다. While the present invention has been particularly shown and described with reference to exemplary embodiments thereof, it is clearly understood that the same is by way of illustration and example only and is not to be taken by way of limitation in the embodiment in which said invention is directed. It will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the scope of the appended claims.

100 : 능동적 윤곽선 방법을 이용한 CAD 시스템을 위한 전처리 장치
200 : 인공물 제거부
300 : 윤곽선 추출부
400 : 흉근 제거부
100: Preprocessing device for CAD system using active contour method
200: Artifact removal
300: contour extractor
400:

Claims (13)

입력된 이미지에 대해 스레쉬홀드를 이용하여 2치 영상(binary image)을 획득하고, 상기 2치 영상을 사전에 설정된 개수로 나누어 분할된 영역에 대해 픽셀 합을 계산하여 포인터의 시작점을 결정하고 상기 시작점을 이용하여 상기 입력 이미지에 대해 인공물을 검색하여 제거하는 인공물 제거부;
흉근 경계 정보를 포함하는 윤곽선을 획득하여 에너지 함수를 정의하고, 상기 입력된 이미지의 최외곽선을 기초로 점차적으로 내부로 이동하면서 사전에 설정된 기준에 해당하는 점들을 만나면 이동을 정지하고 상기 정지된 지점들을 연결하여 하나 이상의 윤곽선을 생성하는 윤곽선 추출부; 및
상기 윤곽선의 마스크 텀(mask term)에 대한 변수를 이용하여 흉근 2치 영상(binary image of the pectoral muscle)을 추출하여 흉근을 제거하는 흉근 제거부;를 포함하는 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 장치.
Determining a starting point of the pointer by calculating a pixel sum with respect to the divided area by dividing the binary image into a predetermined number, obtaining a binary image using a threshold for the input image, An artifact removal unit for searching for and removing artifacts from the input image using a starting point;
And an energy function is defined by acquiring an outline including the rib boundary information, and gradually moves to the inside based on the outermost line of the input image, stops moving when the point corresponding to the predetermined reference is encountered, A contour extractor for concatenating the contour lines to generate at least one contour; And
And removing the pectoral muscle by extracting a binary image of the pectoral muscle using a parameter for the mask term of the contour. A preprocessing device for medical imaging CAD (computer-aided diagnosis) system.
제 1항에 있어서,
상기 인공물 제거부는 상기 입력된 이미지를 분할한 영역별로 픽셀 합을 계산하여 상기 픽셀 합이 가장 큰 영역에서 포인터의 시작점을 결정하는 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 장치.
The method according to claim 1,
Wherein the artifact removal unit computes a pixel sum for each of the divided areas of the input image and determines a starting point of the pointer in the area where the pixel sum is largest. The computer-aided diagnosis (CAD) Preprocessing device for system.
제 2항에 있어서,
상기 포인터의 시작점은 상기 입력된 이미지의 모든 행(row)에 대해 순차적으로 상기 인공물을 검색하는 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 장치.
3. The method of claim 2,
Wherein the starting point of the pointer sequentially retrieves the artifacts for all the rows of the input image. A preprocessing apparatus for a computer-aided diagnosis (CAD) system using the active contour method.
제 1항에 있어서,
상기 윤곽선 추출부는 연속되는 윤곽선 간 픽셀 유사도를 통해 상기 윤곽선의 정지점을 결정하는 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 장치.
The method according to claim 1,
Wherein the contour extracting unit determines a stop point of the contour through a successive contour line pixel similarity. The preprocessing apparatus for a computer-aided diagnosis (CAD) system using the active contour method.
제 1항에 있어서,
상기 흉근 제거부는 상기 흉근 2치 영상을 사전에 설정된 개수로 나누어 분할된 영역에 대해 픽셀 합을 계산하여 포인터의 시작점을 결정하고 상기 시작점을 이용하여 상기 흉근을 제거하는 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 장치.
The method according to claim 1,
Wherein the stomach muscle removing unit calculates a pixel sum of the divided regions by dividing the stereoscopic two-dimensional image into a predetermined number, determines a starting point of the pointer, and removes the stomach muscles using the starting point. A preprocessing system for a computer-aided diagnosis (CAD) system.
제 1항에 있어서,
상기 입력된 이미지는 유방조영상 이미지인 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 장치.
The method according to claim 1,
Wherein the input image is a mammogram image. A preprocessing apparatus for a computer-aided diagnosis (CAD) system using an active contour method.
능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 장치에 의해 수행되는 전처리 방법에 있어서,
(a) 입력된 이미지에 대해 스레쉬홀드를 이용하여 2치 영상(binary image)을 획득하고, 상기 2치 영상을 사전에 설정된 개수로 나누어 분할된 영역에 대해 픽셀 합을 계산하여 포인터의 시작점을 결정하고 상기 시작점을 이용하여 상기 입력 이미지에 대해 인공물을 검색하여 제거하는 인공물 제거 단계;
(b) 흉근 경계 정보를 포함하는 윤곽선을 획득하여 에너지 함수를 정의하고, 상기 입력된 이미지의 최외곽선을 기초로 점차적으로 내부로 이동하면서 사전에 설정된 기준에 해당하는 점들을 만나면 이동을 정지하고 상기 정지된 지점들을 연결하여 하나 이상의 윤곽선을 생성하는 윤곽선 추출 단계; 및
(c) 상기 윤곽선의 마스크 텀(mask term)에 대한 변수를 이용하여 흉근 2치 영상(binary image of the pectoral muscle)을 추출하여 흉근을 제거하는 흉근 제거 단계;를 포함하는 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 방법.
A preprocessing method performed by a preprocessing apparatus for a medical-image CAD system using an active contour method,
(a) obtaining a binary image using a threshold for the input image, dividing the binary image into a predetermined number, calculating a pixel sum for the divided area, and calculating a starting point of the pointer An artifact removal step of searching for and removing artifacts from the input image using the starting point;
(b) defining an energy function by acquiring an outline including the ribbed boundary information, moving gradually toward the inside based on the outermost line of the input image, stopping the movement when a point corresponding to a predetermined reference is encountered, An outline extraction step of connecting the stopped points to generate at least one outline; And
(c) extracting a binary image of the pectoral muscle using a parameter for the mask term of the contour line to remove the pectoral muscle. A pre-processing method for a computer-aided diagnosis (CAD) system using a method.
제 7항에 있어서,
상기 (a) 단계는 상기 입력된 이미지를 분할한 영역별로 픽셀 합을 계산하여 상기 픽셀 합이 가장 큰 영역에서 포인터의 시작점을 결정하는 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 방법.
8. The method of claim 7,
Wherein the step (a) comprises: calculating a pixel sum for each of the divided areas of the input image, and determining a starting point of the pointer in the area having the largest pixel sum; (Preprocessing method for.
제 8항에 있어서,
상기 포인터의 시작점은 상기 입력된 이미지의 모든 행(row)에 대해 순차적으로 상기 인공물을 검색하는 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 방법.
9. The method of claim 8,
Wherein the starting point of the pointer sequentially retrieves the artifacts for all the rows of the input image. A preprocessing method for a computer-aided diagnosis (CAD) system using an active contour method.
제 7항에 있어서,
상기 (b) 단계는 연속되는 윤곽선 간 픽셀 유사도를 통해 상기 윤곽선의 정지점을 결정하는 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 방법.
8. The method of claim 7,
Wherein the step (b) comprises determining a stop point of the contour through successive contour line pixel similarities. A method of pre-processing for a computer-aided diagnosis (CAD) system using an active contour method.
제 7항에 있어서,
상기 (c) 단계는 상기 흉근 2치 영상을 사전에 설정된 개수로 나누어 분할된 영역에 대해 픽셀 합을 계산하여 포인터의 시작점을 결정하고 상기 시작점을 이용하여 상기 흉근을 제거하는 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 방법.
8. The method of claim 7,
Wherein the step (c) comprises the steps of: calculating a sum of pixels for the divided regions by dividing the stereoscopic binary image into a predetermined number, determining a starting point of the pointer, and removing the sick muscles using the starting point; A pre-processing method for a computer-aided diagnosis (CAD) system using a method.
제 7항에 있어서,
상기 입력된 이미지는 유방조영상 이미지인 것을 특징으로 하는 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 방법.
8. The method of claim 7,
Characterized in that the input image is a mammogram image. A preprocessing method for a computer-aided diagnosis (CAD) system using an active contour method.
제 7항 내지 제 12항 중 어느 한 항에 기재된 능동적 윤곽선 방법을 이용한 의료 영상 CAD(computer-aided diagnosis) 시스템을 위한 전처리 방법을 컴퓨터에서 실행시키기 위한 프로그램을 기록한 컴퓨터로 읽을 수 있는 기록매체. A computer-readable recording medium having recorded thereon a program for causing a computer to execute a preprocessing method for a medical image CAD (computer-aided diagnosis) system using the active contour method according to any one of claims 7 to 12.
KR20130071785A 2013-06-21 2013-06-21 Apparatus and method for preprocessing for CAD system using active contour method KR101484051B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR20130071785A KR101484051B1 (en) 2013-06-21 2013-06-21 Apparatus and method for preprocessing for CAD system using active contour method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR20130071785A KR101484051B1 (en) 2013-06-21 2013-06-21 Apparatus and method for preprocessing for CAD system using active contour method

Publications (2)

Publication Number Publication Date
KR20150000067A KR20150000067A (en) 2015-01-02
KR101484051B1 true KR101484051B1 (en) 2015-01-20

Family

ID=52474272

Family Applications (1)

Application Number Title Priority Date Filing Date
KR20130071785A KR101484051B1 (en) 2013-06-21 2013-06-21 Apparatus and method for preprocessing for CAD system using active contour method

Country Status (1)

Country Link
KR (1) KR101484051B1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596191A (en) * 2018-04-23 2018-09-28 福建师范大学 A kind of simple target extracting method for having weak edge

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101874950B1 (en) * 2016-11-23 2018-07-05 연세대학교 산학협력단 Recursive active contour segmentation for 3d cone beam CT and CT device using with the same
CN110956632B (en) * 2020-01-02 2023-04-18 广州柏视医疗科技有限公司 Method and device for automatically detecting pectoralis major region in molybdenum target image

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001522488A (en) 1997-03-27 2001-11-13 アール2 テクノロジー インコーポレイテッド Method and apparatus for automatic muscle segmentation of digital mammograms
JP2005137640A (en) 2003-11-07 2005-06-02 Fuji Photo Film Co Ltd Breast region detecting method, apparatus, and program
KR20120079931A (en) * 2011-01-06 2012-07-16 건국대학교 산학협력단 Method and apparatus for processing segmentation of imaging in medicine, and recording medium having computer program therefor
JP2013118914A (en) 2011-12-07 2013-06-17 Hitachi Medical Corp Medical image diagnosis support device and method of operating the same

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001522488A (en) 1997-03-27 2001-11-13 アール2 テクノロジー インコーポレイテッド Method and apparatus for automatic muscle segmentation of digital mammograms
JP2005137640A (en) 2003-11-07 2005-06-02 Fuji Photo Film Co Ltd Breast region detecting method, apparatus, and program
KR20120079931A (en) * 2011-01-06 2012-07-16 건국대학교 산학협력단 Method and apparatus for processing segmentation of imaging in medicine, and recording medium having computer program therefor
JP2013118914A (en) 2011-12-07 2013-06-17 Hitachi Medical Corp Medical image diagnosis support device and method of operating the same

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108596191A (en) * 2018-04-23 2018-09-28 福建师范大学 A kind of simple target extracting method for having weak edge
CN108596191B (en) * 2018-04-23 2021-06-29 福建师范大学 Method for extracting single target with weak edge

Also Published As

Publication number Publication date
KR20150000067A (en) 2015-01-02

Similar Documents

Publication Publication Date Title
US9996922B2 (en) Image processing of organs depending on organ intensity characteristics
US7315639B2 (en) Method of lung lobe segmentation and computer system
US8958618B2 (en) Method and system for identification of calcification in imaged blood vessels
CN109064476B (en) CT chest radiography lung tissue image segmentation method based on level set
US8090178B2 (en) System and method for automatic detection of internal structures in medical images
EP3352135B1 (en) Method and apparatus for segmentation of blood vessels
Pardo et al. A snake for CT image segmentation integrating region and edge information
CN107103605B (en) Breast tissue segmentation method
CN109478327B (en) Method for automatic detection of systemic arteries in Computed Tomography Angiography (CTA) of arbitrary field of view
US7787673B2 (en) Method and apparatus for airway detection and segmentation using 3D morphological operators
KR20150045885A (en) Systems and methods for registration of ultrasound and ct images
CN110458859B (en) Multi-sequence MRI-based multiple myeloma focus segmentation system
Lee et al. Segmentation of individual ribs from low-dose chest CT
CN105139377A (en) Rapid robustness auto-partitioning method for abdomen computed tomography (CT) sequence image of liver
CN110782489B (en) Image data matching method, device and equipment and computer readable storage medium
JP7278224B2 (en) Quantified aspects of lesions in medical images
WO2022164374A1 (en) Automated measurement of morphometric and geometric parameters of large vessels in computed tomography pulmonary angiography
Maklad et al. Blood vessel‐based liver segmentation using the portal phase of an abdominal CT dataset
KR101484051B1 (en) Apparatus and method for preprocessing for CAD system using active contour method
EP1447772B1 (en) A method of lung lobe segmentation and computer system
Freiman et al. Vessels-cut: a graph based approach to patient-specific carotid arteries modeling
Brahim et al. Semi-automated rib cage segmentation in CT images for mesothelioma detection
Leonardi et al. 3D reconstruction from CT-scan volume dataset application to kidney modeling
CN109564685B (en) Robust lobe segmentation
Wei et al. Segmentation of lung lobes in volumetric CT images for surgical planning of treating lung cancer

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20171227

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20181122

Year of fee payment: 5