JP2006230910A5 - - Google Patents

Download PDF

Info

Publication number
JP2006230910A5
JP2006230910A5 JP2005053493A JP2005053493A JP2006230910A5 JP 2006230910 A5 JP2006230910 A5 JP 2006230910A5 JP 2005053493 A JP2005053493 A JP 2005053493A JP 2005053493 A JP2005053493 A JP 2005053493A JP 2006230910 A5 JP2006230910 A5 JP 2006230910A5
Authority
JP
Japan
Prior art keywords
feature amount
image
image processing
calculated
false positive
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2005053493A
Other languages
Japanese (ja)
Other versions
JP2006230910A (en
Filing date
Publication date
Application filed filed Critical
Priority to JP2005053493A priority Critical patent/JP2006230910A/en
Priority claimed from JP2005053493A external-priority patent/JP2006230910A/en
Publication of JP2006230910A publication Critical patent/JP2006230910A/en
Publication of JP2006230910A5 publication Critical patent/JP2006230910A5/ja
Pending legal-status Critical Current

Links

Claims (35)

医用画像の濃度分布を示す曲面の形状を表す特徴量を算出する第1の特徴量算出手段と、
前記第1の特徴量算出手段により算出された特徴量を用いて、前記医用画像から異常陰影候補領域を検出する検出手段と、
前記検出手段により検出された異常陰影候補領域の特徴量を計算する第2の特徴量算出手段と、
前記第2の特徴量算出手段により算出された特徴量に基づいて前記異常陰影候補領域の偽陽性を削除する偽陽性削除手段と、
を備えることを特徴とする画像処理装置。
A first feature amount calculating means for calculating a feature amount representing a shape of a curved surface indicating a density distribution of a medical image;
Detecting means for detecting an abnormal shadow candidate region from the medical image using the feature amount calculated by the first feature amount calculating unit;
Second feature amount calculating means for calculating a feature amount of the abnormal shadow candidate region detected by the detecting means;
False positive deletion means for deleting false positives of the abnormal shadow candidate region based on the feature quantity calculated by the second feature quantity calculation means;
An image processing apparatus comprising:
前記第1の特徴量算出手段は、前記曲面の形状を表す特徴量としてシェイプインデックスを算出することを特徴とする請求項1に記載の画像処理装置。   The image processing apparatus according to claim 1, wherein the first feature amount calculating unit calculates a shape index as a feature amount representing the shape of the curved surface. 前記第2の特徴量算出手段により算出される特徴量には、シェイプインデックスが含まれることを特徴とする請求項1又は2に記載の画像処理装置。   The image processing apparatus according to claim 1, wherein the feature quantity calculated by the second feature quantity calculation unit includes a shape index. 前記検出手段は、前記第1の特徴量算出手段により算出された特徴量に基づいて、予め設定された当該特徴量の閾値を用いて前記医用画像を二値化することにより、異常陰影候補領域を検出し、
前記偽陽性削除手段は、前記二値化により得られた異常陰影候補領域から偽陽性の画像領域を削除することを特徴とする請求項1〜3の何れか一項に記載の画像処理装置。
The detection unit binarizes the medical image using a preset threshold value of the feature amount based on the feature amount calculated by the first feature amount calculation unit, so that an abnormal shadow candidate region is obtained. Detect
The image processing apparatus according to claim 1, wherein the false positive deletion unit deletes a false positive image region from the abnormal shadow candidate region obtained by the binarization.
前記偽陽性削除手段による偽陽性の削除により得られた画像を表示する表示手段を備えることを特徴とする請求項1〜4の何れか一項に記載の画像処理装置。   The image processing apparatus according to claim 1, further comprising a display unit that displays an image obtained by deleting false positives by the false positive deleting unit. 医用画像における任意の注目画素を設定するとともに、その注目画素から所定範囲内の画像領域を設定する設定手段を備え、
前記第1の特徴量算出手段は、前記設定された画像領域における濃度分布の曲面形状を表す特徴量を算出することを特徴とする請求項1〜5の何れか一項に記載の画像処理装置。
A setting unit for setting an arbitrary pixel of interest in the medical image and setting an image region within a predetermined range from the pixel of interest,
The image processing apparatus according to claim 1, wherein the first feature amount calculation unit calculates a feature amount representing a curved surface shape of a density distribution in the set image region. .
前記医用画像の濃度分布を示す曲面上の注目画素における法線を回転軸として法平面を所定角度ずつ回転させ、その回転角度毎に、注目画素から所定範囲内の画像領域の曲面を法平面で切り出すことによって得られた曲線を近似する関数を算出する関数算出手段を備え、
前記第1の特徴量算出手段は、前記関数算出手段により算出された近似関数に基づいて特徴量を算出することを特徴とする請求項1〜6の何れか一項に記載の画像処理装置。
The normal plane is rotated by a predetermined angle around the normal line of the pixel of interest on the curved surface showing the density distribution of the medical image as a rotation axis, and the curved surface of the image area within the predetermined range from the pixel of interest is plotted on the normal plane for each rotation angle. A function calculating means for calculating a function approximating a curve obtained by cutting out,
The image processing apparatus according to claim 1, wherein the first feature amount calculating unit calculates a feature amount based on the approximate function calculated by the function calculating unit.
前記関数算出手段は、最小二乗法を用いて前記曲線を近似する関数を算出することを特徴とする請求項7に記載の画像処理装置。   The image processing apparatus according to claim 7, wherein the function calculating unit calculates a function that approximates the curve using a least square method. 前記関数算出手段は、前記近似関数として近似円を算出し、
前記第1の特徴量算出手段は、回転角度毎に、前記算出された近似円の半径から前記曲線の注目画素における曲率を算出し、回転角度毎に算出された曲率に基づいて特徴量を算出することを特徴とする請求項7又は8に記載の画像処理装置。
The function calculating means calculates an approximate circle as the approximate function,
The first feature amount calculation means calculates a curvature at the target pixel of the curve from the calculated approximate circle radius for each rotation angle, and calculates a feature amount based on the curvature calculated for each rotation angle. The image processing apparatus according to claim 7, wherein the image processing apparatus is an image processing apparatus.
前記第1の特徴量算出手段は、前記医用画像の全ての画素について特徴量を算出することを特徴とする請求項1〜9の何れか一項に記載の画像処理装置。   The image processing apparatus according to claim 1, wherein the first feature amount calculating unit calculates feature amounts for all pixels of the medical image. 前記偽陽性削除手段は、前記異常陰影候補領域のうち、大きさが所定範囲外の画像領域を偽陽性の画像領域として削除することを特徴とする請求項1〜10の何れか一項に記載の画像処理装置。   The said false positive deletion means deletes the image area | region outside a predetermined range as a false positive image area among the said abnormal shadow candidate area | regions, The Claim 1 characterized by the above-mentioned. Image processing apparatus. 前記偽陽性削除手段は、前記異常陰影候補領域のうち、形状が所定の形状以外の画像領域を偽陽性の画像領域として削除することを特徴とする請求項1〜10の何れか一項に記載の画像処理装置。   The false positive deletion means deletes an image area whose shape is other than a predetermined shape from the abnormal shadow candidate areas as a false positive image area. Image processing apparatus. 前記偽陽性削除手段は、前記第1及び/又は第2の特徴量算出手段で算出された特徴量を含む前記医用画像の各種特徴量を入力データとし、真陽性の異常陰影の度合いを示す指標値を出力データとする多変量解析を用いて、偽陽性の画像領域を削除することを特徴とする請求項1〜12の何れか一項に記載の画像処理装置。   The false positive deletion means uses various feature quantities of the medical image including the feature quantities calculated by the first and / or second feature quantity calculation means as input data, and indicates an index indicating the degree of true positive abnormal shadow The image processing apparatus according to claim 1, wherein false-positive image regions are deleted using multivariate analysis using values as output data. 前記第1の特徴量算出手段は、特徴量を算出する曲面の領域範囲を変更し、その変更された領域範囲で特徴量を算出することを特徴とする請求項1〜13の何れか一項に記載の画像処理装置。   The first feature amount calculation unit changes a region range of a curved surface from which a feature amount is calculated, and calculates a feature amount in the changed region range. An image processing apparatus according to 1. 前記第1の特徴量算出手段は、前記医用画像の解像度を変更し、その変更された解像度で特徴量を算出することを特徴とする請求項1〜13の何れか一項に記載の画像処理装置。   The image processing according to any one of claims 1 to 13, wherein the first feature amount calculating unit changes a resolution of the medical image and calculates a feature amount with the changed resolution. apparatus. 前記偽陽性削除手段による偽陽性の削除より、異常陰影としてガウス分布状の曲面を構成する画像領域が検出されることを特徴とする請求項1〜15の何れか一項に記載の画像処理装置。   The image processing apparatus according to claim 1, wherein an image area constituting a curved surface with a Gaussian distribution is detected as an abnormal shadow from the deletion of false positives by the false positive deletion unit. . 前記偽陽性削除手段による偽陽性の削除より、異常陰影として腫瘤陰影の画像領域が検出されることを特徴とする請求項16に記載の画像処理装置。   The image processing apparatus according to claim 16, wherein an image area of a tumor shadow is detected as an abnormal shadow from the false positive deletion by the false positive deletion means. 医用画像の濃度分布を示す曲面の形状を表す特徴量を算出する第1の特徴量算出工程と、
前記第1の特徴量算出工程において算出された特徴量を用いて、前記医用画像から異常陰影候補領域を検出する検出工程と、
前記検出工程において検出された異常陰影候補領域の特徴量を計算する第2の特徴量算出工程と、
前記第2の特徴量算出工程において算出された特徴量に基づいて前記異常陰影候補領域の偽陽性を削除する偽陽性削除工程と、
を含むことを特徴とする画像処理方法。
A first feature amount calculating step of calculating a feature amount representing a shape of a curved surface indicating a density distribution of a medical image;
A detection step of detecting an abnormal shadow candidate region from the medical image using the feature amount calculated in the first feature amount calculation step;
A second feature amount calculating step of calculating a feature amount of the abnormal shadow candidate region detected in the detection step;
A false positive deletion step of deleting false positives of the abnormal shadow candidate region based on the feature amount calculated in the second feature amount calculation step;
An image processing method comprising:
前記第1の特徴量算出工程では、前記曲面の形状を表す特徴量としてシェイプインデックスが算出されることを特徴とする請求項18に記載の画像処理方法。   The image processing method according to claim 18, wherein in the first feature amount calculating step, a shape index is calculated as a feature amount representing the shape of the curved surface. 前記第2の特徴量算出工程において算出される特徴量には、シェイプインデックスが含まれることを特徴とする請求項18又は19に記載の画像処理方法。   20. The image processing method according to claim 18, wherein the feature quantity calculated in the second feature quantity calculation step includes a shape index. 前記検出工程では、前記第1の特徴量算出工程において算出された特徴量に基づいて、予め設定された当該特徴量の閾値を用いて前記医用画像を二値化することにより、異常陰影候補領域が検出され、
前記偽陽性削除工程では、前記二値化により得られた異常陰影候補領域から偽陽性の画像領域が削除されることを特徴とする請求項18〜20の何れか一項に記載の画像処理方法。
In the detection step, based on the feature amount calculated in the first feature amount calculation step, the medical image is binarized using a preset threshold value of the feature amount, whereby an abnormal shadow candidate region is obtained. Is detected,
The image processing method according to any one of claims 18 to 20, wherein in the false positive deletion step, a false positive image region is deleted from the abnormal shadow candidate region obtained by the binarization. .
前記偽陽性削除工程による偽陽性の削除により得られた画像を表示する表示工程を含むことを特徴とする請求項18〜21の何れか一項に記載の画像処理方法。   The image processing method according to any one of claims 18 to 21, further comprising a display step of displaying an image obtained by deleting false positives in the false positive deletion step. 医用画像における任意の注目画素を設定するとともに、その注目画素から所定範囲内の画像領域を設定する設定工程を含み、
前記第1の特徴量算出工程では、前記設定された画像領域における濃度分布の曲面形状を表す特徴量が算出されることを特徴とする請求項18〜22の何れか一項に記載の画像処理方法。
A setting step of setting an arbitrary target pixel in the medical image and setting an image region within a predetermined range from the target pixel;
The image processing according to any one of claims 18 to 22, wherein, in the first feature amount calculation step, a feature amount representing a curved surface shape of a density distribution in the set image region is calculated. Method.
前記医用画像の濃度分布を示す曲面上の注目画素における法線を回転軸として法平面を所定角度ずつ回転させ、その回転角度毎に、注目画素から所定範囲内の画像領域の曲面を法平面で切り出すことによって得られた曲線を近似する関数を算出する関数算出工程を含み、
前記第1の特徴量算出工程では、前記関数算出工程において算出された近似関数に基づいて特徴量が算出されることを特徴とする請求項18〜23の何れか一項に記載の画像処理方法。
The normal plane is rotated by a predetermined angle around the normal line of the pixel of interest on the curved surface showing the density distribution of the medical image as a rotation axis, and the curved surface of the image area within the predetermined range from the pixel of interest is plotted on the normal plane for each rotation angle. Including a function calculation step of calculating a function approximating a curve obtained by cutting,
The image processing method according to any one of claims 18 to 23, wherein in the first feature amount calculation step, a feature amount is calculated based on the approximate function calculated in the function calculation step. .
前記関数算出工程では、最小二乗法を用いて前記曲線を近似する関数が算出されることを特徴とする請求項24に記載の画像処理方法。   The image processing method according to claim 24, wherein in the function calculating step, a function that approximates the curve is calculated using a least square method. 前記関数算出工程では、前記近似関数として近似円が算出され、
前記第1の特徴量算出工程では、回転角度毎に、前記算出された近似円の半径から前記曲線の注目画素における曲率が算出され、回転角度毎に算出された曲率に基づいて特徴量が算出されることを特徴とする請求項24又は25に記載の画像処理方法。
In the function calculating step, an approximate circle is calculated as the approximate function,
In the first feature amount calculation step, the curvature of the target pixel of the curve is calculated from the calculated radius of the approximate circle for each rotation angle, and the feature amount is calculated based on the curvature calculated for each rotation angle. The image processing method according to claim 24 or 25, wherein:
前記第1の特徴量算出工程では、前記医用画像の全ての画素について特徴量が算出されることを特徴とする請求項18〜26の何れか一項に記載の画像処理方法。   The image processing method according to any one of claims 18 to 26, wherein in the first feature amount calculation step, a feature amount is calculated for all pixels of the medical image. 前記偽陽性削除工程では、前記異常陰影候補領域のうち、大きさが所定範囲外の画像領域が偽陽性の画像領域として削除されることを特徴とする請求項18〜27の何れか一項に記載の画像処理方法。   28. The false positive deletion step, wherein, among the abnormal shadow candidate regions, an image region whose size is outside a predetermined range is deleted as a false positive image region. The image processing method as described. 前記偽陽性削除工程では、前記異常陰影候補領域のうち、形状が所定の形状以外の画像領域が偽陽性の画像領域として削除されることを特徴とする請求項18〜27の何れか一項に記載の画像処理方法。   28. The false positive deletion step, wherein an image region whose shape is other than a predetermined shape is deleted as a false positive image region among the abnormal shadow candidate regions. The image processing method as described. 前記偽陽性削除工程では、前記第1及び/又は第2の特徴量算出工程で算出された特徴量を含む前記医用画像の各種特徴量を入力データとし、真陽性の異常陰影の度合いを示す指標値を出力データとする多変量解析を用いて、偽陽性の画像領域が削除されることを特徴とする請求項18〜29の何れか一項に記載の画像処理方法。   In the false positive deletion step, various feature amounts of the medical image including the feature amount calculated in the first and / or second feature amount calculation step are used as input data, and an index indicating the degree of true positive abnormal shadow The image processing method according to any one of claims 18 to 29, wherein a false positive image region is deleted using multivariate analysis using a value as output data. 前記第1の特徴量算出工程では、特徴量を算出する曲面の領域範囲が変更され、その変更された領域範囲で特徴量が算出されることを特徴とする請求項18〜30の何れか一項に記載の画像処理方法。   31. The method according to claim 18, wherein in the first feature amount calculation step, a region range of the curved surface for calculating the feature amount is changed, and the feature amount is calculated in the changed region range. The image processing method according to item. 前記第1の特徴量算出工程では、前記医用画像の解像度が変更され、その変更された解像度で特徴量が算出されることを特徴とする請求項18〜30の何れか一項に記載の画像処理方法。   The image according to any one of claims 18 to 30, wherein, in the first feature amount calculation step, the resolution of the medical image is changed, and the feature amount is calculated with the changed resolution. Processing method. 前記偽陽性削除工程における偽陽性の削除より、異常陰影としてガウス分布状の曲面を構成する画像領域が検出されることを特徴とする請求項18〜32の何れか一項に記載の画像処理方法。   The image processing method according to any one of claims 18 to 32, wherein an image area constituting a curved surface with a Gaussian distribution is detected as an abnormal shadow from the false positive deletion in the false positive deletion step. . 前記偽陽性削除工程における偽陽性の削除より、異常陰影として腫瘤陰影の画像領域が検出されることを特徴とする請求項33に記載の画像処理方法。   34. The image processing method according to claim 33, wherein an image region of a tumor shadow is detected as an abnormal shadow from the false positive deletion in the false positive deletion step. コンピュータを、Computer
医用画像の濃度分布を示す曲面の形状を表す特徴量を算出する第1の特徴量算出手段、A first feature amount calculating means for calculating a feature amount representing a shape of a curved surface indicating a density distribution of a medical image;
前記第1の特徴量算出手段により算出された特徴量を用いて、前記医用画像から異常陰影候補領域を検出する検出手段、Detecting means for detecting an abnormal shadow candidate region from the medical image using the feature quantity calculated by the first feature quantity calculating means;
前記検出手段により検出された異常陰影候補領域の特徴量を計算する第2の特徴量算出手段、A second feature amount calculating means for calculating a feature amount of the abnormal shadow candidate region detected by the detecting means;
前記第2の特徴量算出手段により算出された特徴量に基づいて前記異常陰影候補領域の偽陽性を削除する偽陽性削除手段、False positive deletion means for deleting false positives of the abnormal shadow candidate region based on the feature quantity calculated by the second feature quantity calculation means;
として機能させるためのプログラム。Program to function as.
JP2005053493A 2005-02-28 2005-02-28 Image processor and image processing method Pending JP2006230910A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005053493A JP2006230910A (en) 2005-02-28 2005-02-28 Image processor and image processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005053493A JP2006230910A (en) 2005-02-28 2005-02-28 Image processor and image processing method

Publications (2)

Publication Number Publication Date
JP2006230910A JP2006230910A (en) 2006-09-07
JP2006230910A5 true JP2006230910A5 (en) 2008-04-10

Family

ID=37039175

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005053493A Pending JP2006230910A (en) 2005-02-28 2005-02-28 Image processor and image processing method

Country Status (1)

Country Link
JP (1) JP2006230910A (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5078486B2 (en) * 2007-07-26 2012-11-21 オリンパスメディカルシステムズ株式会社 Medical image processing apparatus and method of operating medical image processing apparatus
WO2010024331A1 (en) * 2008-09-01 2010-03-04 株式会社 日立メディコ Image processing device and method for processing image
JP5566299B2 (en) * 2008-10-20 2014-08-06 株式会社日立メディコ Medical image processing apparatus and medical image processing method
CN102298780B (en) * 2011-08-15 2012-12-12 天津大学 Method for detecting shadow of color image
JP6356528B2 (en) * 2014-08-06 2018-07-11 株式会社日立製作所 Ultrasonic diagnostic equipment
CN110070545B (en) * 2019-03-20 2023-05-26 重庆邮电大学 Method for automatically extracting urban built-up area by urban texture feature density

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001195566A (en) * 1999-10-26 2001-07-19 Hitachi Medical Corp Image sorting device
JP4124406B2 (en) * 2001-06-13 2008-07-23 富士フイルム株式会社 Abnormal shadow detection device
JP2005506140A (en) * 2001-10-16 2005-03-03 ザ・ユニバーシティー・オブ・シカゴ Computer-aided 3D lesion detection method
JP4604451B2 (en) * 2003-02-24 2011-01-05 コニカミノルタホールディングス株式会社 Medical image processing apparatus and malignancy determination method
US7333644B2 (en) * 2003-03-11 2008-02-19 Siemens Medical Solutions Usa, Inc. Systems and methods for providing automatic 3D lesion segmentation and measurements
JP2004344232A (en) * 2003-05-20 2004-12-09 Konica Minolta Medical & Graphic Inc Medical image processor and method of detecting abnormal shadow candidate
JP2004351056A (en) * 2003-05-30 2004-12-16 Konica Minolta Medical & Graphic Inc Medical image processing apparatus
JP2005080757A (en) * 2003-09-05 2005-03-31 Konica Minolta Medical & Graphic Inc Signal processor
JP2005080758A (en) * 2003-09-05 2005-03-31 Konica Minolta Medical & Graphic Inc Image processing apparatus

Similar Documents

Publication Publication Date Title
Wang et al. Automated estimation of reinforced precast concrete rebar positions using colored laser scan data
CN110378900B (en) Method, device and system for detecting product defects
Karasik et al. 3D scanning technology as a standard archaeological tool for pottery analysis: practice and theory
Arena et al. A new computational approach to cracks quantification from 2D image analysis: Application to micro-cracks description in rocks
EP3105664B1 (en) Determining pitch and yaw for touchscreen interactions
US9117281B2 (en) Surface segmentation from RGB and depth images
JP7319432B2 (en) LEARNING DATA COLLECTION DEVICE, LEARNING DATA COLLECTION METHOD, AND PROGRAM
JP2018112550A (en) Inspection apparatus and method for detecting firearms
Kim et al. Automated rebar diameter classification using point cloud data based machine learning
JP2006230910A5 (en)
Zhang Digital image processing method for estimating leaf length and width tested using kiwifruit leaves (Actinidia chinensis Planch)
JP2018181333A5 (en)
CN115937167A (en) Method and device for detecting defects of battery pole piece ceramic and electronic equipment
Lin et al. Surface defect detection of machined parts based on machining texture direction
Thurley et al. Identifying, visualizing, and comparing regions in irregularly spaced 3D surface data
JP6629685B2 (en) Crack analysis method and crack analysis system
KR20230042706A (en) Neural network analysis of LFA test strips
JP4505509B2 (en) Weather information detection method, weather information detection apparatus, and weather information detection program
US9305233B2 (en) Isotropic feature matching
JP2019190998A (en) Cavity identification system and cavity identification method
Liu et al. Automated building change detection using UltraCamD images and existing CAD data
CN101718674A (en) Method for measuring shape parameters of particles of particulate materials
JP2018523967A (en) Method and apparatus for determining the temporal behavior of an object
Li et al. Automatic crack detection method for loaded coal in vibration failure process
Noskov Computer vision approaches for big geo-spatial data: quality assessment of raster tiled web maps for smart city solutions