Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP4634872B2 - Ultrasonic diagnostic equipment - Google Patents
[go: Go Back, main page]

JP4634872B2 - Ultrasonic diagnostic equipment - Google Patents

Ultrasonic diagnostic equipment Download PDF

Info

Publication number
JP4634872B2
JP4634872B2 JP2005172030A JP2005172030A JP4634872B2 JP 4634872 B2 JP4634872 B2 JP 4634872B2 JP 2005172030 A JP2005172030 A JP 2005172030A JP 2005172030 A JP2005172030 A JP 2005172030A JP 4634872 B2 JP4634872 B2 JP 4634872B2
Authority
JP
Japan
Prior art keywords
initial
endpoint
end point
pixel data
ultrasonic diagnostic
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.)
Expired - Fee Related
Application number
JP2005172030A
Other languages
Japanese (ja)
Other versions
JP2006345893A (en
Inventor
賢 村下
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Aloka Co Ltd
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 Aloka Co Ltd filed Critical Aloka Co Ltd
Priority to JP2005172030A priority Critical patent/JP4634872B2/en
Publication of JP2006345893A publication Critical patent/JP2006345893A/en
Application granted granted Critical
Publication of JP4634872B2 publication Critical patent/JP4634872B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

本発明は、超音波診断装置に関し、特に超音波画像内において対象組織の境界を求める超音波診断装置に関する。   The present invention relates to an ultrasonic diagnostic apparatus, and more particularly to an ultrasonic diagnostic apparatus that obtains a boundary of a target tissue in an ultrasonic image.

超音波画像に映し出された臓器や腫瘍などの対象組織の大きさを計測することは診断上において重要である。計測項目としては、例えば、対象組織の二点間の距離、周囲長、面積、体積などが挙げられる。また、股関節角度、産科計測、ヒストグラム、LV計測なども知られている。   It is important for diagnosis to measure the size of a target tissue such as an organ or a tumor projected on an ultrasonic image. Examples of the measurement item include a distance between two points of the target tissue, a perimeter, an area, and a volume. Also known are hip joint angles, obstetric measurements, histograms, LV measurements, and the like.

超音波診断装置によって臓器や腫瘍などの内径を計測する場合、モニタ上に映し出された画像を医師が目視しながら、トラックボールなどのポインティングデバイスを利用して内径端と思われる場所、つまり、臓器境界にポイントを設定する手法が知られている。この場合、目視によって正確な境界を検出するのは困難であり、また、正確に境界を設定しようとすると、医師などのユーザによる位置指定のための操作に時間がかかるなどの問題がある。   When measuring the inner diameter of an organ, tumor, etc. with an ultrasound diagnostic device, the doctor looks at the image displayed on the monitor and uses a pointing device such as a trackball to see the inner diameter end, that is, the organ A method for setting points at boundaries is known. In this case, it is difficult to detect an accurate boundary by visual observation, and there is a problem that it takes time for an operation for specifying a position by a user such as a doctor if an attempt is made to set the boundary accurately.

このため、従来から、対象組織の境界を検出するための補助的な機能を備えた超音波診断装置が知られている。   For this reason, conventionally, an ultrasonic diagnostic apparatus having an auxiliary function for detecting a boundary of a target tissue is known.

例えば、特許文献1には、臓器を横切る線分上で画像情報の変化度合いを検出し、検出された変化度合いから臓器の境界を判定する技術が示されている。つまり、臓器の境界部分において画像情報が大きく変化することを利用して境界を判定しようとする技術が示されている。   For example, Patent Document 1 discloses a technique for detecting the degree of change in image information on a line segment crossing an organ and determining the boundary of the organ from the detected degree of change. That is, a technique for determining a boundary by utilizing a large change in image information in a boundary portion of an organ is shown.

また、特許文献2には、境界抽出のための開始点を設定してスキャンラインを放射状に延出し、そのライン上で手動による境界との交点を中心とした検出範囲内において、画像データを閾値で2値化処理して補正すべき境界位置を検出する技術が示されている。   Patent Document 2 discloses that a start point for boundary extraction is set, a scan line is radially extended, and image data is set as a threshold value within a detection range centered on an intersection with the boundary manually on the line. Shows a technique for detecting a boundary position to be corrected by binarization processing.

特開平8−627号公報JP-A-8-627 特開平11−164834号公報JP-A-11-164834

臓器によっては、画像情報の変化度合いを利用する手法や画像データを閾値で2値化処理する手法では、正確に所望の境界を抽出することができない場合がある。例えば、所望の境界の近傍に他の境界が存在する場合などである。この場合、特許文献1に記載の技術によって画像情報の変化度合いを利用して境界を検出しても、複数の境界が検出されてしまい、それら複数の境界のうちから検査者が求める境界を判断するのは容易ではない。この問題は、特許文献2に記載の技術においても、検出範囲内で境界位置が複数検出された場合に起こり得ることである。このため、臓器や腫瘍などの対象組織の境界を正確に設定することができる技術が望まれていた。   Depending on the organ, a method that uses the degree of change in image information or a method that binarizes image data with a threshold may not be able to accurately extract a desired boundary. For example, there is a case where another boundary exists in the vicinity of a desired boundary. In this case, even if the boundary is detected using the change degree of the image information by the technique described in Patent Document 1, a plurality of boundaries are detected, and the boundary required by the inspector is determined from the plurality of boundaries. It is not easy to do. This problem is also possible in the technique described in Patent Document 2 when a plurality of boundary positions are detected within the detection range. For this reason, the technique which can set the boundary of target tissues, such as an organ and a tumor correctly, was desired.

本発明は、このような背景において成されたものであり、その目的は、超音波画像内に対象組織の境界を正確に設定する技術を提供することにある。   The present invention has been made in such a background, and an object thereof is to provide a technique for accurately setting a boundary of a target tissue in an ultrasonic image.

上記目的を達成するために、本発明の好適な態様である超音波診断装置は、対象組織を含む空間内に超音波を送受波してエコーデータを取得する送受波手段と、前記エコーデータに基づいて対象組織を含む超音波画像の画像データを形成する画像形成手段と、対象組織の境界に対応した端点として前記画像データ内に初期端点を設定する初期端点設定手段と、前記画像データ内の複数の画素データについて、各画素データごとにその近傍の画素データとの間における相違の度合いを示した指標値を算出する指標値算出手段と、前記複数の画素データについての複数の指標値に対して、前記初期端点からの距離に応じた重み付け処理を施す重み付け処理手段と、前記重み付け処理された複数の指標値に基づいて、前記初期端点の近傍領域内において指標値が最大となる画素データを検出し、検出された画素データの位置に修正された端点として最終端点を設定する最終端点設定手段と、を有する、ことを特徴とする。   In order to achieve the above object, an ultrasonic diagnostic apparatus according to a preferred aspect of the present invention includes a transmission / reception unit that transmits and receives an ultrasonic wave in a space including a target tissue to acquire echo data, and the echo data An image forming means for forming image data of an ultrasound image including the target tissue based on the image, an initial end point setting means for setting an initial end point in the image data as an end point corresponding to the boundary of the target tissue, For a plurality of pixel data, an index value calculation means for calculating an index value indicating the degree of difference between the pixel data for each pixel data and its neighboring pixel data, and a plurality of index values for the plurality of pixel data A weighting processing means for performing weighting processing in accordance with a distance from the initial endpoint, and a plurality of the weighted index values in a region near the initial endpoint. Index value detects the pixel data having the largest has a final end point setting means for setting a final end point as were fixed to the detected pixel data position endpoints, and it is characterized by Te.

上記構成では、重み付け処理手段によって、複数の指標値に対して初期端点からの距離に応じた重み付け処理が施されている。例えば、初期端点の位置において最大の重み付け係数とし、初期端点から離れるに従って重み付け係数を小さくすることによって、初期端点に近い位置に存在する画素データの指標値ほど大きく重み付けされる。このため、最終端点設定手段において、初期端点の近傍領域内で指標値が最大となる画素データを検出することによって、初期端点に最も近い位置に存在する指標値のピークを検出することができる。つまり、初期端点の近傍に複数の指標値のピーク(極値)が存在する場合でも、初期端点に最も近いピークに大きな重み付けがなされているため、初期端点に最も近い指標値のピークを的確に検出することができる。その結果、対象組織の境界を正確に設定することができる。   In the above-described configuration, the weighting processing unit performs weighting processing corresponding to the distance from the initial endpoint with respect to the plurality of index values. For example, by setting the maximum weighting coefficient at the position of the initial end point and decreasing the weighting coefficient as the distance from the initial end point increases, the index value of the pixel data existing at a position close to the initial end point is weighted more. For this reason, the final endpoint setting means can detect the peak of the index value present at the position closest to the initial endpoint by detecting the pixel data having the maximum index value in the region near the initial endpoint. In other words, even when there are multiple index value peaks (extreme values) near the initial endpoint, the peak closest to the initial endpoint is heavily weighted, so the peak of the index value closest to the initial endpoint is accurately determined. Can be detected. As a result, the boundary of the target tissue can be set accurately.

望ましくは、前記指標値算出手段は、前記指標値として、各画素データについての空間的な微分値を算出する、ことを特徴とする。   Preferably, the index value calculation means calculates a spatial differential value for each pixel data as the index value.

望ましくは、前記初期端点設定手段は、ユーザからの指示に基づいて前記端点として初期端点Aおよび初期端点Bを設定し、前記指標値算出手段は、前記初期端点Aおよび初期端点Bを通る直線上の各画素データについての空間的な微分値を算出し、前記重み付け処理手段は、前記直線上の複数の画素データに対応した複数の微分値に対して、前記初期端点Aからの距離に応じた重み付け処理および前記初期端点Bからの距離に応じた重み付け処理を施し、前記最終端点設定手段は、前記重み付け処理された複数の微分値に基づいて、前記初期端点Aの近傍領域内において微分値が最大となる画素データを検出してその位置に最終端点Cを設定し、さらに、前記初期端点Bの近傍領域内において微分値が最大となる画素データを検出してその位置に最終端点Dを設定し、前記最終端点Cと最終端点Dとの間の距離を計測する距離計測手段をさらに有する、ことを特徴とする。   Preferably, the initial endpoint setting means sets an initial endpoint A and an initial endpoint B as the endpoint based on an instruction from a user, and the index value calculation means is a straight line passing through the initial endpoint A and the initial endpoint B. The spatial differential value for each of the pixel data is calculated, and the weighting processing unit responds to the distance from the initial end point A for the plurality of differential values corresponding to the plurality of pixel data on the straight line. A weighting process and a weighting process corresponding to the distance from the initial end point B are performed, and the final end point setting means has a differential value within a region near the initial end point A based on the plurality of differential values subjected to the weighting process. The maximum pixel data is detected and the final end point C is set at the position, and further, the pixel data having the maximum differential value in the region near the initial end point B is detected and Set the final end point D to the location, the further comprising a distance measuring means for measuring a distance between the final end point C and the final end point D, characterized in that.

望ましくは、前記指標値算出手段は、前記指標値として、各画素データとその近傍の画素データからなる複数の画素データについての分散値を算出する、ことを特徴とする。   Preferably, the index value calculation means calculates a variance value for a plurality of pixel data including each pixel data and neighboring pixel data as the index value.

望ましくは、前記初期端点設定手段は、ユーザからの指示に基づいて前記端点として初期端点Aおよび初期端点Bを設定し、前記重み付け処理手段は、前記複数の画素データについての複数の分散値に対して、前記初期端点Aからの距離に応じた重み付け処理および前記初期端点Bからの距離に応じた重み付け処理を施し、前記最終端点設定手段は、前記重み付け処理された複数の分散値に基づいて、前記初期端点Aの近傍領域内において分散値が最大となる画素データを検出してその位置に最終端点Cを設定し、さらに、前記初期端点Bの近傍領域内において分散値が最大となる画素データを検出してその位置に最終端点Dを設定し、前記最終端点Cと最終端点Dとの間の距離を計測する距離計測手段をさらに有する、ことを特徴とする。   Preferably, the initial end point setting means sets an initial end point A and an initial end point B as the end points based on an instruction from a user, and the weighting processing means applies to a plurality of variance values for the plurality of pixel data. The weighting process according to the distance from the initial endpoint A and the weighting process according to the distance from the initial endpoint B are performed, and the final endpoint setting means is based on the plurality of weighted variance values, Pixel data having a maximum variance value in the region near the initial endpoint A is detected, the final endpoint C is set at that position, and pixel data having a maximum variance value in the region near the initial endpoint B And a distance measuring means for setting the final end point D at the position and measuring the distance between the final end point C and the final end point D.

本発明により、超音波画像内に対象組織の境界を正確に設定することができる。   According to the present invention, the boundary of the target tissue can be accurately set in the ultrasonic image.

以下、本発明の好適な実施形態を図面に基づいて説明する。   DESCRIPTION OF EXEMPLARY EMBODIMENTS Hereinafter, preferred embodiments of the invention will be described with reference to the drawings.

図1には、本発明に係る超音波診断装置の好適な実施形態が示されており、図1はその全体構成を示すブロック図である。   FIG. 1 shows a preferred embodiment of an ultrasonic diagnostic apparatus according to the present invention, and FIG. 1 is a block diagram showing the overall configuration thereof.

プローブ10は、臓器や腫瘍などの対象組織を含む空間内に超音波を送受波する超音波探触子である。プローブ10の内部には、図示しない複数の振動素子が含まれており、それら複数の振動素子によって超音波ビームが電子的に走査される。プローブ10は、三次元空間内のエコーデータ(ボリュームデータ)を取得するタイプ、二次元平面内のエコーデータを取得するタイプのいずれでもよい。また、プローブ10は、生体の体表に接して用いられるタイプ、生体の体腔内に挿入して用いられるタイプのいずれでもよい。   The probe 10 is an ultrasonic probe that transmits and receives ultrasonic waves in a space including target tissues such as organs and tumors. The probe 10 includes a plurality of vibration elements (not shown), and the ultrasonic beam is electronically scanned by the plurality of vibration elements. The probe 10 may be either a type that acquires echo data (volume data) in a three-dimensional space or a type that acquires echo data in a two-dimensional plane. The probe 10 may be either a type used in contact with a body surface of a living body or a type used by being inserted into a body cavity of a living body.

送受信部20は、プローブ10を制御して超音波ビームを走査して受信信号(エコーデータ)を取得する。送受信部20は、送信ビームフォーマおよび受信ビームフォーマとして機能し、複数の振動素子の各々に対して送信信号を供給し、また、複数の振動素子の各々からの受信信号に対して整相加算処理を実行する。   The transmission / reception unit 20 controls the probe 10 to scan the ultrasonic beam and obtain a reception signal (echo data). The transmission / reception unit 20 functions as a transmission beamformer and a reception beamformer, supplies a transmission signal to each of the plurality of vibration elements, and performs phasing addition processing on the reception signal from each of the plurality of vibration elements Execute.

画像形成部30は、送受信部20から供給される受信信号(エコーデータ)に基づいて対象組織を含む超音波画像の画像データを形成する。送受信部20から供給される受信信号は、各超音波ビームごとの受信信号として供給され、超音波の送受波の座標系(例えば、rθφ極座標系)上のデータとして供給される。このため、画像形成部30は、受信信号の座標変換処理を行って、超音波画像の表示座標系(例えば、xyz直交座標系)に変換する処理を行う。この際、必要に応じて補間処理なども実行する。   The image forming unit 30 forms image data of an ultrasonic image including the target tissue based on the reception signal (echo data) supplied from the transmission / reception unit 20. The reception signal supplied from the transmission / reception unit 20 is supplied as a reception signal for each ultrasonic beam, and is supplied as data on a coordinate system (for example, rθφ polar coordinate system) of ultrasonic transmission / reception. For this reason, the image forming unit 30 performs a process of converting the received signal into a display coordinate system (for example, an xyz orthogonal coordinate system) by performing a coordinate conversion process on the received signal. At this time, interpolation processing or the like is also executed as necessary.

このように、画像形成部30は、超音波の送受波の座標系から超音波画像の表示座標系への変換処理を実行するデジタルスキャンコンバータとしての機能を備えている。なお、プローブ10から三次元空間内のエコーデータ(ボリュームデータ)が取得されている場合、画像形成部30は、三次元空間内における所定の断面の超音波画像の画像データを形成する。例えば、直交三断面の断層画像データを形成する。また、プローブ10から二次元平面内のエコーデータが取得されている場合、画像形成部30は、その二次元平面に対応する断層画像データを形成する。   In this manner, the image forming unit 30 has a function as a digital scan converter that executes a conversion process from the coordinate system of ultrasonic transmission / reception to the display coordinate system of ultrasonic images. When echo data (volume data) in a three-dimensional space is acquired from the probe 10, the image forming unit 30 forms image data of an ultrasonic image of a predetermined cross section in the three-dimensional space. For example, tomographic image data having three orthogonal cross sections is formed. When echo data in a two-dimensional plane is acquired from the probe 10, the image forming unit 30 forms tomographic image data corresponding to the two-dimensional plane.

画像形成部30で形成された画像データ(断層画像データ)は、モニタ40に供給され、モニタ40に当該画像データに対応した超音波画像が表示される。また、画像形成部30で形成された画像データは、境界検出ブロック50にも供給される。   Image data (tomographic image data) formed by the image forming unit 30 is supplied to the monitor 40, and an ultrasonic image corresponding to the image data is displayed on the monitor 40. The image data formed by the image forming unit 30 is also supplied to the boundary detection block 50.

境界検出ブロック50は、超音波画像の画像データに含まれる対象組織の境界に対応した端点を検出して設定する。前述のとおり、プローブ10により対象組織を含む空間内に超音波が送受波されるため、画像形成部30において、臓器や腫瘍などの対象組織が含まれた断層画像データが形成される。境界検出ブロック50は、断層画像データ上において対象組織の境界を正確に抽出して端点を設定する。そこで、以下、図2から図12を利用して、境界検出ブロック50の動作を説明する。なお、以下において図1に示した部分には図1の符号を付して説明する。   The boundary detection block 50 detects and sets an end point corresponding to the boundary of the target tissue included in the image data of the ultrasonic image. As described above, since ultrasonic waves are transmitted and received by the probe 10 in the space including the target tissue, the image forming unit 30 forms tomographic image data including the target tissue such as an organ or a tumor. The boundary detection block 50 accurately extracts the boundary of the target tissue on the tomographic image data and sets the end points. Therefore, the operation of the boundary detection block 50 will be described below with reference to FIGS. In the following description, the parts shown in FIG.

図2は、境界検出ブロック50の動作を説明するためのフローチャートである。図2のフローチャートの各ステップにおける処理は次のとおりである。   FIG. 2 is a flowchart for explaining the operation of the boundary detection block 50. The processing in each step of the flowchart of FIG. 2 is as follows.

まず、初期端点設定部51によって初期端点が設定される(S201)。医師などのユーザは、モニタ40に表示された対象組織の断層画像を見ながら、トラックボールやマウスなどのデバイスを利用して、断層画像内に対象組織の境界に対応した端点の位置を指定する。   First, an initial endpoint is set by the initial endpoint setting unit 51 (S201). A user such as a doctor designates the position of the end point corresponding to the boundary of the target tissue in the tomographic image using a device such as a trackball or a mouse while viewing the tomographic image of the target tissue displayed on the monitor 40. .

図3は、モニタ40に表示された断層画像の表示例を示しており、臓器や腫瘍などの対象組織60,60´が映し出されている。図3の断層画像において、左側の対象組織60の長軸長を計測したい場合、ユーザは、トラックボールやマウスなどのデバイスを利用して、断層画像内に対象組織60の境界に対応した端点の位置にA点およびB点を指定する。図3においては、図示の都合上、対象組織60とその周囲の組織との境界がはっきり示されているが、実際にモニタ40に表示される断層画像においては、境界部分がぼやけていることもあり、また表示画像が小さいなどの要因から、モニタ40を目視することのみで境界を正確に特定することは困難である。このため、図3においては、本来の境界から少しずれた位置にA点およびB点が指定されたと仮定する。ちなみに本実施形態においては、これらA点およびB点を基準に正確な境界の位置が探索される。初期端点設定部51は、ユーザからの指示に基づいて、端点として初期端点Aおよび初期端点Bを設定する。   FIG. 3 shows a display example of a tomographic image displayed on the monitor 40, and target tissues 60 and 60 'such as an organ and a tumor are displayed. In the tomographic image of FIG. 3, when the user wants to measure the long axis length of the left target tissue 60, the user uses a device such as a trackball or a mouse to find the end point corresponding to the boundary of the target tissue 60 in the tomographic image. Designate points A and B as the positions. In FIG. 3, for the convenience of illustration, the boundary between the target tissue 60 and the surrounding tissue is clearly shown. However, in the tomographic image actually displayed on the monitor 40, the boundary portion may be blurred. In addition, due to factors such as a small display image, it is difficult to accurately identify the boundary only by visually observing the monitor 40. For this reason, in FIG. 3, it is assumed that the points A and B are designated at positions slightly deviated from the original boundary. Incidentally, in the present embodiment, an accurate boundary position is searched based on these points A and B. The initial endpoint setting unit 51 sets an initial endpoint A and an initial endpoint B as endpoints based on an instruction from the user.

図2に戻り、初期端点A,Bが設定されると、画像データ取得部52は、A点とB点を結ぶ直線上の画素データ(輝度値)を画像形成部30内のメモリ32から読み込む(S202)。つまり、画像データ取得部52は、図3に示す直線70上に存在する画素の画素データ(輝度値)をメモリ32から読み込む。この際、メモリ32から読み込んだデータに補間処理などを施して、直線70上に相当する画素の輝度値を求めてもよい。   Returning to FIG. 2, when the initial end points A and B are set, the image data acquisition unit 52 reads pixel data (luminance value) on a straight line connecting the points A and B from the memory 32 in the image forming unit 30. (S202). That is, the image data acquisition unit 52 reads pixel data (luminance values) of pixels existing on the straight line 70 shown in FIG. At this time, the luminance value of the pixel corresponding to the straight line 70 may be obtained by performing interpolation processing or the like on the data read from the memory 32.

図4は、画像データ取得部52によって読み込まれた輝度値を説明するための図であり、図4において、横軸は直線(図3の符号70)上の位置を示しており、縦軸は直線(図3の符号70)上の各位置の画素の輝度値を示している。図4において横軸の位置A,Bは、各々、図3における点A,Bの位置を示している。   FIG. 4 is a diagram for explaining the luminance value read by the image data acquisition unit 52. In FIG. 4, the horizontal axis indicates a position on a straight line (reference numeral 70 in FIG. 3), and the vertical axis indicates The luminance value of the pixel of each position on the straight line (code | symbol 70 of FIG. 3) is shown. In FIG. 4, positions A and B on the horizontal axis indicate the positions of points A and B in FIG. 3, respectively.

図2に戻り、A点とB点を結ぶ直線上の輝度値が読み込まれると、平滑化処理部53は、読み込まれた直線上の輝度値(直線データ)に対して平滑化処理を施す(S203)。画像データ取得部52によって読み込まれた図4の直線データは、送受波状態の微小変動やノイズなどの影響により、局所的な変動を伴っている。そこで、図4の直線データに対して平滑化処理を施して、その局所的な変動をキャンセルする。   Returning to FIG. 2, when the luminance value on the straight line connecting the points A and B is read, the smoothing processing unit 53 performs a smoothing process on the luminance value (straight line data) on the read straight line ( S203). The straight line data of FIG. 4 read by the image data acquisition unit 52 is accompanied by local fluctuations due to the influence of minute fluctuations in the transmission / reception state, noise, and the like. Therefore, a smoothing process is performed on the straight line data in FIG. 4 to cancel the local variation.

図5は、平滑化処理が施された直線上の輝度値(直線データ)を説明するための図であり、図4と同様に、横軸は直線(図3の符号70)上の位置を示しており、縦軸は直線(図3の符号70)上の各位置の画素の平滑処理後の輝度値を示している。平滑化は、例えば、平均値処理やメディアン処理によって実行される。   FIG. 5 is a diagram for explaining a luminance value (straight line data) on a straight line on which smoothing processing has been performed. As in FIG. 4, the horizontal axis indicates a position on a straight line (reference numeral 70 in FIG. 3). The vertical axis indicates the luminance value after smoothing processing of pixels at each position on a straight line (reference numeral 70 in FIG. 3). Smoothing is performed by, for example, average value processing or median processing.

図2に戻り、平滑化処理が施されると、指標値算出部54は、平滑化後の直線データに対して線形一次微分フィルタ処理を施す(S204)。図9は、指標値算出部54によって実行される微分処理を説明するための図である。図9において、Yi−3,Yi−2,Yi−1,Y,・・・は、各々、平滑化処理が施された直線上の輝度値を示している。微分値は、対象画素が輝度値Yの場合、輝度値Yの画素に隣接する画素の輝度値から|Yi+1−Yi−1|で与えられる。指標値算出部54は、平滑化処理が施された直線上の全データに関する微分値を演算する。 Returning to FIG. 2, when the smoothing process is performed, the index value calculation unit 54 performs a linear first-order differential filtering process on the smoothed straight line data (S204). FIG. 9 is a diagram for explaining differentiation processing executed by the index value calculation unit 54. In FIG. 9, Y i−3 , Y i−2 , Y i−1 , Y i ,... Each indicate a luminance value on a straight line subjected to smoothing processing. Differential value, the target pixel is the case of luminance values Y i, from the luminance values of pixels adjacent to the pixel of the luminance value Y i | is given by | Y i + 1 -Y i- 1. The index value calculation unit 54 calculates differential values for all data on the straight line that has been subjected to the smoothing process.

図6は、微分処理が施された直線上の輝度値(直線データ)を説明するための図であり、横軸は直線(図3の符号70)上の位置を示しており、縦軸は直線(図3の符号70)上の各位置の微分値を示している。微分処理の性質上、輝度値の濃淡差が大きい位置において微分値が大きくなっている。つまり、図3に示す対象組織60および対象組織60´の各々とその周囲の組織との境界部分で微分値が大きい値を示している。ただし、ユーザによって指定されたA点,B点は、本来の境界部分からずれた位置に設定されている。   FIG. 6 is a diagram for explaining a luminance value (straight line data) on a straight line on which differentiation processing has been performed. The horizontal axis indicates a position on the straight line (reference numeral 70 in FIG. 3), and the vertical axis indicates The differential value of each position on the straight line (reference numeral 70 in FIG. 3) is shown. Due to the nature of the differential processing, the differential value is large at the position where the brightness difference of the luminance value is large. That is, the differential value is large at the boundary portion between each of the target tissue 60 and the target tissue 60 ′ shown in FIG. 3 and the surrounding tissue. However, the points A and B designated by the user are set at positions shifted from the original boundary portion.

なお、指標値算出部54において次のような微分処理を実行してもよい。図10は、指標値算出部54によって実行される別の微分処理を説明するための図である。図10において、Yi−3,Yi−2,Yi−1,Y,・・・は、各々、平滑化処理が施された直線上の輝度値を示している。図10に示す微分処理では、直線上の輝度値に加えて、その直線に隣接する直線上の輝度値も利用される。つまり、対象となる直線の上側に位置する直線上の輝度値Zi−3,Zi−2,Zi−1,Z,・・・、下側に位置する直線上の輝度値Xi−3,Xi−2,Xi−1,X,・・・を利用する。そして、対象画素が輝度値Yの場合に、微分値は、|(Zi+1+Yi+1+Xi+1)−(Zi−1+Yi−1+Xi−1)|で与えられる。なお、輝度値Yが属する直線上の値に重み付けを施して、微分値を|(Zi+1+2Yi+1+Xi+1)−(Zi−1+2Yi−1+Xi−1)|で求めてもよい。 The index value calculation unit 54 may execute the following differentiation process. FIG. 10 is a diagram for explaining another differentiation process executed by the index value calculation unit 54. 10, Y i−3 , Y i−2 , Y i−1 , Y i ,... Each indicate a luminance value on a straight line subjected to smoothing processing. In the differentiation process shown in FIG. 10, in addition to the luminance value on the straight line, the luminance value on the straight line adjacent to the straight line is also used. That is, the luminance values Z i-3 , Z i-2 , Z i-1 , Z i ,... On the straight line located above the target straight line, the luminance values X i on the straight line located below. -3 , X i-2 , X i-1 , X i ,. When the target pixel has a luminance value Y i , the differential value is given by | (Z i + 1 + Y i + 1 + X i + 1 ) − (Z i−1 + Y i−1 + X i−1 ) |. Note that the differential value may be obtained by | (Z i + 1 + 2Y i + 1 + X i + 1 ) − (Z i−1 + 2Y i−1 + X i−1 ) | by weighting the value on the straight line to which the luminance value Y i belongs. Good.

図2に戻り、微分処理が施されると、重み付け処理部55は、微分処理後の直線データに対して端点からの距離に応じた重み付け処理を施す(S205)。重み付け処理部55は、図6に示す微分処理後の直線データに対して、ユーザによって指定されたB点を中心として、B点において最大の重み付け係数を与え、B点から離れるに従って重み付け係数を徐々に小さくする。例えば、B点において重み付け係数を「5」として、B点から1画素分だけ離れた位置において重み付け係数「4」を与える。さらに、1画素分だけ離れるに従って重み付け係数を「3」,「2」,「1」,「0」と段階的に小さくする。   Returning to FIG. 2, when the differentiation process is performed, the weighting processing unit 55 performs the weighting process according to the distance from the end point to the linear data after the differentiation process (S205). The weighting processing unit 55 gives the maximum weighting coefficient at the point B centered on the point B designated by the user to the linear data after the differentiation process shown in FIG. 6 and gradually increases the weighting coefficient as the distance from the point B increases. Make it smaller. For example, the weighting coefficient is “5” at point B, and the weighting coefficient “4” is given at a position one pixel away from point B. Further, the weighting coefficient is gradually reduced to “3”, “2”, “1”, “0” as the distance increases by one pixel.

重み付け処理部55は、ユーザによって指定されたA点を中心とした重み付け処理も行う。つまり、例えば、A点において重み付け係数を「5」として、A点から1画素分だけ離れた位置において重み付け係数「4」を与える。さらに、1画素分だけ離れるに従って重み付け係数を「3」,「2」,「1」,「0」と段階的に小さくする。   The weighting processing unit 55 also performs weighting processing around the point A designated by the user. In other words, for example, the weighting coefficient is “5” at point A, and the weighting coefficient “4” is given at a position one pixel away from point A. Further, the weighting coefficient is gradually reduced to “3”, “2”, “1”, “0” as the distance increases by one pixel.

図7は、重み付け処理が施された直線上の微分値を説明するための図であり、図6と同様に、横軸は直線(図3の符号70)上の位置を示しており、縦軸は直線(図3の符号70)上の各位置の微分値を示している。重み付け処理の結果、ユーザによって指定されたA点,B点の近くに位置する微分値のピークに比べて、A点,B点から遠くに位置する微分値のピーク(図7においてB点よりも右側に位置するピーク)が小さく抑えられている。   FIG. 7 is a diagram for explaining a differential value on a straight line on which weighting processing has been performed. Like FIG. 6, the horizontal axis indicates a position on the straight line (reference numeral 70 in FIG. 3), and the vertical axis The axis indicates the differential value at each position on the straight line (reference numeral 70 in FIG. 3). As a result of the weighting process, the peak of the differential value located farther from the points A and B than the peak of the differential value located near the points A and B specified by the user (in FIG. The peak located on the right side) is kept small.

なお、重み付け処理部55における重み付け処理は、重み付け係数「5」〜「0」の6段階の重み付けに限定されるものではなく、6段階よりも多く、あるいは、6段階よりも少ない段階設定でもよい。また、A点やB点を中心として1画素分だけ離れるに従って重み付け係数を段階的に小さくしていく手法に換えて、複数画素分を単位として重み係数を段階的に小さくするなどでもよい。   Note that the weighting processing in the weighting processing unit 55 is not limited to the six-stage weighting of the weighting coefficients “5” to “0”, and may be set with more stages than six stages or less than six stages. . Further, instead of the method of decreasing the weighting coefficient in steps as the distance from the A point and the B point becomes one pixel, the weighting coefficient may be decreased in steps for a plurality of pixels.

図2に戻り、重み付け処理が施されると、最終端点設定部56は、A点およびB点の各々の近傍領域内において微分値が最大となる画素データ(輝度値)を検出し、検出された画素データの位置に修正された端点として最終端点を設定する(S206)。   Returning to FIG. 2, when the weighting process is performed, the final endpoint setting unit 56 detects and detects pixel data (luminance value) having a maximum differential value in the vicinity of each of the points A and B. The final endpoint is set as the corrected endpoint at the position of the pixel data (S206).

図8は、最終端点を説明するための図であり、図7と同様に、横軸は直線(図3の符号70)上の位置を示しており、縦軸は直線(図3の符号70)上の各位置の重み付け処理後の微分値を示している。最終端点設定部56は、A点(図7参照)の近傍領域内において微分値が最大となる輝度値の位置に最終端点Cを設定する。A点の近傍領域としては、例えば、A点から所定画素数以内の領域とする。同様に、最終端点設定部56は、B点(図7参照)の近傍領域内(所定画素数以内)において微分値が最大となる輝度値の位置に最終端点Dを設定する。   FIG. 8 is a diagram for explaining the final end point. Like FIG. 7, the horizontal axis indicates the position on the straight line (reference numeral 70 in FIG. 3), and the vertical axis indicates the straight line (reference numeral 70 in FIG. 3). ) The differential value after the weighting process of each position is shown. The final end point setting unit 56 sets the final end point C at the position of the luminance value at which the differential value is maximized in the vicinity region of the point A (see FIG. 7). For example, the vicinity of the point A is an area within a predetermined number of pixels from the point A. Similarly, the final end point setting unit 56 sets the final end point D at the position of the luminance value at which the differential value is maximum within the vicinity of the point B (see FIG. 7) (within a predetermined number of pixels).

図7を利用して説明したように、重み付け処理の結果、ユーザによって指定されたA点,B点の近くに位置する微分値のピークに比べて、A点,B点から遠くに位置する微分値のピーク(図7においてB点よりも右側に位置するピーク)が小さく抑えられている。このため、A点およびB点の各々の近傍領域内において微分値が最大となる輝度値を検出することによって、対象組織(図3の符号60)本来の境界部分に相当する最終端点Cおよび最終端点Dが適切に抽出される。   As described with reference to FIG. 7, as a result of the weighting process, the derivative located farther from the points A and B than the peak of the derivative located near the points A and B specified by the user. The value peak (the peak located on the right side of the point B in FIG. 7) is kept small. For this reason, by detecting the luminance value having the maximum differential value in the vicinity region of each of the points A and B, the final end point C and the final value corresponding to the original boundary portion of the target tissue (reference numeral 60 in FIG. 3) are detected. The end point D is appropriately extracted.

なお、図8において、最終端点Cおよび最終端点Dは、各々、微分値の山のピーク、つまり画像の濃淡差が大きい部分に設定されているが、例えば、濃淡差が大きくなり始める微分値の山の立ち上がり部分などに最終端点を設定してもよい。   In FIG. 8, the final end point C and the final end point D are each set to a peak of the differential value peak, that is, a portion having a large shade difference of the image. The final end point may be set at the rising part of the mountain.

図2に戻り、最終端点が設定されると、距離計測部57は、最終端点Cと最終端点Dとの間の距離を計測する(S207)。上述したように、対象組織(図3の符号60)本来の境界部分に相当する最終端点Cおよび最終端点Dが適切に抽出されているため、対象組織の境界間の距離(長軸長)が正確に計測される。つまり、ユーザが当初設定した端点A,Bの各々の位置が修正され、修正後の最終端点C,Dから正確な長軸長が計測される。そして、対象組織の画像や距離の計測結果がモニタ40に表示される(S208)。   Returning to FIG. 2, when the final endpoint is set, the distance measuring unit 57 measures the distance between the final endpoint C and the final endpoint D (S207). As described above, since the final end point C and the final end point D corresponding to the original boundary portion of the target tissue (reference numeral 60 in FIG. 3) are appropriately extracted, the distance (major axis length) between the target tissue boundaries is long. Accurately measured. That is, the positions of the end points A and B initially set by the user are corrected, and an accurate major axis length is measured from the corrected final end points C and D. Then, the image of the target tissue and the measurement result of the distance are displayed on the monitor 40 (S208).

図11は、最終端点C,Dが設定された後のモニタ40の表示例を示している。図11に示すように、対象組織60の断層画像上に対象組織60の長軸長72が表示されている。なお、図11に示す表示例に加えて、さらに長軸長72の計測値が表示されてもよい。   FIG. 11 shows a display example of the monitor 40 after the final end points C and D are set. As shown in FIG. 11, the long axis length 72 of the target tissue 60 is displayed on the tomographic image of the target tissue 60. In addition to the display example shown in FIG. 11, the measurement value of the long axis length 72 may be further displayed.

以上のようにして、境界検出ブロック50において対象組織の境界が正確に抽出され、その結果、対象組織の長軸長などを正確に計測することができる。なお、図2を利用して説明した境界検出動作においては、ユーザによって指定された点A,Bを通る直線上において最終端点C,Dを探索しているが、点A,Bを通る直線上に限定されずに、次のように、点A,Bに最も近い位置に最終端点を設定するようにしてもよい。   As described above, the boundary of the target tissue is accurately extracted in the boundary detection block 50. As a result, the major axis length of the target tissue can be accurately measured. In the boundary detection operation described with reference to FIG. 2, the final end points C and D are searched for on the straight line passing through the points A and B designated by the user. Without being limited thereto, the final end point may be set at a position closest to the points A and B as follows.

図12は、境界検出ブロック50の別の動作を説明するためのフローチャートである。 まず、初期端点設定部51によって初期端点が設定される(S1201)。つまり、図2のS201における処理と同じ動作によって、初期端点設定部51は、ユーザからの指示に基づいて、端点として初期端点Aおよび初期端点Bを設定する。   FIG. 12 is a flowchart for explaining another operation of the boundary detection block 50. First, an initial endpoint is set by the initial endpoint setting unit 51 (S1201). That is, the initial endpoint setting unit 51 sets the initial endpoint A and the initial endpoint B as endpoints based on an instruction from the user by the same operation as the processing in S201 of FIG.

初期端点A,Bが設定されると、指標値算出部54は、初期端点A,Bそれぞれの周辺各点における分散値を計算する(S1202)。指標値算出部54は、例えば、初期端点A,Bの各々を中心とする二次元的なサーチエリアを設定する。そして、初期端点Aを中心とするサーチエリア内の各点(各画素の位置)において分散値を計算し、また、初期端点Bを中心とするサーチエリア内の各点(各画素の位置)において分散値を計算する。サーチエリアは、本来の境界が含まれる程度の大きさであればよいが、超音波画像全体をサーチエリアとしてもよい。各点における分散値は、例えば、その点を中心として縦3画素、横3画素の9画素からなる参照エリアを設定し、その参照エリア内の9つの画素の輝度値の分散値とする。参照エリアは、縦5画素、横5画素の25画素からなるエリアなどでもよい。   When the initial end points A and B are set, the index value calculation unit 54 calculates a variance value at each of the peripheral points around the initial end points A and B (S1202). For example, the index value calculation unit 54 sets a two-dimensional search area centered on each of the initial end points A and B. Then, a variance value is calculated at each point (position of each pixel) in the search area centered on the initial end point A, and at each point (position of each pixel) in the search area centered on the initial end point B. Calculate the variance value. The search area may be large enough to include the original boundary, but the entire ultrasonic image may be used as the search area. For the variance value at each point, for example, a reference area consisting of 9 pixels of 3 pixels in the vertical direction and 3 pixels in the horizontal direction is set as the center, and the variance value of the luminance values of 9 pixels in the reference area is set. The reference area may be an area composed of 25 pixels of 5 pixels vertically and 5 pixels horizontally.

初期端点A,Bそれぞれの周辺各点における分散値が計算されると、重み付け処理部55は、各点における分散値に対して端点からの距離に応じた重み付け処理を施す(S1203)。つまり、初期端点Aを中心とするサーチエリア内において、A点に対して最大の重み付け係数を与え、A点から離れるに従って重み付け係数を徐々に小さくする。例えば、A点において重み付け係数を「5」として、A点を取り囲む8画素の位置において重み付け係数「4」を与える。さらに、A点から1画素分だけ離れるに従って重み付け係数を「3」,「2」,「1」,「0」と段階的に小さくして、二次元的に重み付け処理を施す。なお、初期端点Bを中心とするサーチエリア内においても、初期端点Aの場合と同様に二次元的に重み付け処理を施す。   When the variance values at the respective peripheral points of the initial endpoints A and B are calculated, the weighting processing unit 55 performs a weighting process according to the distance from the endpoint on the variance value at each point (S1203). That is, in the search area centered on the initial end point A, the maximum weighting coefficient is given to the point A, and the weighting coefficient is gradually reduced as the distance from the point A increases. For example, the weighting coefficient is “5” at point A, and the weighting coefficient “4” is given at the positions of 8 pixels surrounding point A. Further, the weighting coefficient is reduced stepwise as “3”, “2”, “1”, “0” as the distance from the point A by one pixel, and weighting processing is performed two-dimensionally. Even in the search area centered on the initial end point B, the weighting process is two-dimensionally performed as in the case of the initial end point A.

重み付け処理が施されると、最終端点設定部56は、A点およびB点の各々のサーチエリア内において分散値が最大となる点(画素データの位置)を検出し、検出された画素データの位置に修正された端点として最終端点を設定する(S1204)。重み付け処理部55による重み付け処理の結果、ユーザによって指定されたA点,B点の近くの点における分散値に比べて、A点,B点から遠い点における分散値が小さく抑えられている。つまり、A点,B点から遠い位置に存在する対象組織(図3の符号60´)の境界における分散値が小さく抑えられる。このため、A点およびB点の各々のサーチエリア内において分散値が最大となる点を二次元的に検索することによって、対象組織(図3の符号60)本来の境界部分が適切に抽出される。最終端点設定部56は、A点に対応する修正された端点として最終端点Cを設定し、B点に対応する修正された端点として最終端点Dを設定する。   When the weighting process is performed, the final endpoint setting unit 56 detects the point (pixel data position) at which the variance value is maximum in each of the search areas of the points A and B, and the detected pixel data The final end point is set as the end point corrected to the position (S1204). As a result of the weighting processing by the weighting processing unit 55, the dispersion values at points far from the points A and B are suppressed to be smaller than the dispersion values at points near the points A and B specified by the user. That is, the dispersion value at the boundary of the target tissue (reference numeral 60 ′ in FIG. 3) existing at a position far from the points A and B is suppressed. For this reason, the original boundary portion of the target tissue (reference numeral 60 in FIG. 3) is appropriately extracted by two-dimensionally searching for the point having the maximum variance value in the search area of each of the points A and B. The The final endpoint setting unit 56 sets the final endpoint C as the corrected endpoint corresponding to the point A, and sets the final endpoint D as the corrected endpoint corresponding to the point B.

最終端点が設定されると、距離計測部57は、最終端点Cと最終端点Dとの間の距離を計測する(S1205)。対象組織(図3の符号60)本来の境界部分に相当する最終端点Cおよび最終端点Dが適切に抽出されているため、対象組織の境界間の距離(長軸長)が正確に計測される。つまり、ユーザが当初設定した端点A,Bの各々の位置が二次元的に修正され、修正後の最終端点C,Dから正確な長軸長が計測される。そして、対象組織の画像や距離の計測結果がモニタ40に表示される(S1206)。   When the final endpoint is set, the distance measuring unit 57 measures the distance between the final endpoint C and the final endpoint D (S1205). Since the final endpoint C and the final endpoint D corresponding to the original boundary portion of the target tissue (reference numeral 60 in FIG. 3) are appropriately extracted, the distance (major axis length) between the boundaries of the target tissue is accurately measured. . That is, the positions of the end points A and B initially set by the user are two-dimensionally corrected, and an accurate major axis length is measured from the corrected final end points C and D. Then, the image of the target tissue and the measurement result of the distance are displayed on the monitor 40 (S1206).

図12を利用して説明した動作により、直線上に限定されずに、端点A,Bの各々に最も近い位置の濃淡差の大きい箇所を最終端点C,Dとすることができる。なお、図12のS1202において、サーチエリア内の各点の分散値に換えて、ラプラシアンフィルタなどを利用して各点における二次元的な勾配を求め、それを分散値に代わる指標値としてもよい。   By the operation described with reference to FIG. 12, a point having a large shade difference at a position closest to each of the end points A and B can be set as the final end points C and D without being limited to a straight line. In S1202 of FIG. 12, instead of the variance value of each point in the search area, a two-dimensional gradient at each point may be obtained using a Laplacian filter or the like, and this may be used as an index value instead of the variance value. .

以上、本発明の好適な実施形態を説明したが、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。   As mentioned above, although preferred embodiment of this invention was described, embodiment mentioned above is only a mere illustration in all the points, and does not limit the scope of the present invention.

本発明に係る超音波診断装置の全体構成を示すブロック図である。1 is a block diagram showing an overall configuration of an ultrasonic diagnostic apparatus according to the present invention. 境界検出ブロックの動作を説明するためのフローチャートである。It is a flowchart for demonstrating operation | movement of a boundary detection block. モニタに表示された断層画像の表示例を示す図である。It is a figure which shows the example of a display of the tomographic image displayed on the monitor. 画像データ取得部によって読み込まれた輝度値を説明するための図である。It is a figure for demonstrating the luminance value read by the image data acquisition part. 平滑化処理が施された直線上の輝度値を説明するための図である。It is a figure for demonstrating the luminance value on the straight line in which the smoothing process was performed. 微分処理が施された直線上の輝度値を説明するための図である。It is a figure for demonstrating the luminance value on the straight line in which the differentiation process was performed. 重み付け処理が施された直線上の微分値を説明するための図である。It is a figure for demonstrating the differential value on the straight line to which the weighting process was performed. 最終端点を説明するための図である。It is a figure for demonstrating the last endpoint. 指標値算出部における微分処理を説明するための図である。It is a figure for demonstrating the differentiation process in an index value calculation part. 指標値算出部における別の微分処理を説明するための図である。It is a figure for demonstrating another differentiation process in an index value calculation part. 最終端点C,Dが設定された後のモニタの表示例を示す図である。It is a figure which shows the example of a display of the monitor after the last end points C and D are set. 境界検出ブロックの別の動作を説明するためのフローチャートである。It is a flowchart for demonstrating another operation | movement of a boundary detection block.

符号の説明Explanation of symbols

30 画像形成部、51 初期端点設定部、52 画像データ取得部、53 平滑化処理部、54 指標値算出部、55 重み付け処理部、56 最終端点設定部、57 距離計測部。   30 image forming unit, 51 initial endpoint setting unit, 52 image data acquisition unit, 53 smoothing processing unit, 54 index value calculation unit, 55 weighting processing unit, 56 final endpoint setting unit, 57 distance measurement unit.

Claims (5)

対象組織を含む空間内に超音波を送受波してエコーデータを取得する送受波手段と、
前記エコーデータに基づいて対象組織を含む超音波画像の画像データを形成する画像形成手段と、
対象組織の境界に対応した端点として前記画像データ内に初期端点Aおよび初期端点Bを設定する初期端点設定手段と、
前記初期端点Aおよび初期端点Bを通る直線上の複数の画素データについて、各画素データごとにその近傍の画素データとの間における相違の度合いを示した指標値を算出する指標値算出手段と、
前記複数の画素データについての複数の指標値に対して、前記初期端点からの距離に応じた重み付け処理および前記初期端点Bからの距離に応じた重み付け処理を施す重み付け処理手段と、
前記重み付け処理された複数の指標値に基づいて、前記初期端点の近傍領域内において指標値が最大となる前記直線上の画素データを検出してその位置に最終端点Cを設定し、さらに、前記初期端点Bの近傍領域内において指標値が最大となる前記直線上の画素データを検出してその位置に最終端点Dを設定する最終端点設定手段と、
前記最終端点Cと最終端点Dとの間の距離を計測する距離計測手段と、
を有する、
ことを特徴とする超音波診断装置。
Wave transmitting / receiving means for transmitting and receiving ultrasonic waves in a space including the target tissue to acquire echo data;
Image forming means for forming image data of an ultrasonic image including the target tissue based on the echo data;
An initial endpoint setting means for setting an initial endpoint A and an initial endpoint B in the image data as endpoints corresponding to the boundary of the target tissue;
An index value calculation means for calculating an index value indicating a degree of difference between each of the pixel data on the straight line passing through the initial endpoint A and the initial endpoint B and the neighboring pixel data for each pixel data;
Weighting processing means for applying a weighting process according to a distance from the initial endpoint A and a weighting process according to a distance from the initial endpoint B to a plurality of index values for the plurality of pixel data;
Based on the plurality of weighted index values, the pixel data on the straight line having the maximum index value in the vicinity of the initial endpoint A is detected, and the final endpoint C is set at that position. A final endpoint setting means for detecting pixel data on the straight line having the maximum index value in a region near the initial endpoint B and setting the final endpoint D at the position ;
Distance measuring means for measuring a distance between the final end point C and the final end point D;
Having
An ultrasonic diagnostic apparatus.
請求項1に記載の超音波診断装置において、
前記指標値算出手段は、前記指標値として、各画素データについての空間的な微分値を算出する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1,
The index value calculation means calculates a spatial differential value for each pixel data as the index value.
An ultrasonic diagnostic apparatus.
請求項1または2に記載の超音波診断装置において、
前記初期端点設定手段は、ユーザからの指示に基づいて前記初期端点Aおよび初期端点Bを設定する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1 or 2,
It said initial end point setting means sets the initial end point A and an initial end point B on the basis of an instruction from the user,
An ultrasonic diagnostic apparatus.
請求項1に記載の超音波診断装置において、
前記指標値算出手段は、前記指標値として、各画素データとその近傍の画素データからなる複数の画素データについての分散値を算出する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1,
The index value calculation means calculates a variance value for a plurality of pixel data including each pixel data and neighboring pixel data as the index value.
An ultrasonic diagnostic apparatus.
請求項4に記載の超音波診断装置において、
前記初期端点設定手段は、ユーザからの指示に基づいて前記初期端点Aおよび初期端点Bを設定する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 4,
It said initial end point setting means sets the initial end point A and an initial end point B on the basis of an instruction from the user,
An ultrasonic diagnostic apparatus.
JP2005172030A 2005-06-13 2005-06-13 Ultrasonic diagnostic equipment Expired - Fee Related JP4634872B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005172030A JP4634872B2 (en) 2005-06-13 2005-06-13 Ultrasonic diagnostic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005172030A JP4634872B2 (en) 2005-06-13 2005-06-13 Ultrasonic diagnostic equipment

Publications (2)

Publication Number Publication Date
JP2006345893A JP2006345893A (en) 2006-12-28
JP4634872B2 true JP4634872B2 (en) 2011-02-16

Family

ID=37642371

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005172030A Expired - Fee Related JP4634872B2 (en) 2005-06-13 2005-06-13 Ultrasonic diagnostic equipment

Country Status (1)

Country Link
JP (1) JP4634872B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2012088226A (en) * 2010-10-21 2012-05-10 Jtekt Corp Nondestructive inspection method and nondestructive inspection system
JP7157710B2 (en) * 2019-07-10 2022-10-20 富士フイルム株式会社 Measuring device, ultrasonic diagnostic device, measuring method, measuring program

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5620443A (en) * 1979-07-31 1981-02-26 Tokyo Shibaura Electric Co Ultrasonic diagnosing device
JP3325983B2 (en) * 1993-12-27 2002-09-17 アロカ株式会社 Ultrasonic image processing apparatus and method
JP3563904B2 (en) * 1996-12-16 2004-09-08 ジーイー横河メディカルシステム株式会社 Ultrasound diagnostic apparatus and wall thickness measuring method
JP3355298B2 (en) * 1998-02-05 2002-12-09 松下電器産業株式会社 Ultrasound diagnostic equipment
JP3330090B2 (en) * 1998-09-30 2002-09-30 松下電器産業株式会社 Organ boundary extraction method and apparatus
JP4219456B2 (en) * 1998-12-07 2009-02-04 株式会社日立製作所 Region extraction method and region extraction device
JP3689627B2 (en) * 2000-09-19 2005-08-31 アロカ株式会社 Ultrasonic image processing device

Also Published As

Publication number Publication date
JP2006345893A (en) 2006-12-28

Similar Documents

Publication Publication Date Title
JP5606203B2 (en) Ultrasonic system and sensor coordinate calibration method
JP5579527B2 (en) System and method for providing 2D CT images corresponding to 2D ultrasound images
JP5498299B2 (en) System and method for providing 2D CT images corresponding to 2D ultrasound images
JP5160227B2 (en) Ultrasonic diagnostic apparatus and ultrasonic image display method
JP5587332B2 (en) Ultrasonic imaging apparatus and program for ultrasonic imaging
JP4751282B2 (en) Ultrasonic diagnostic equipment
EP1913874B1 (en) Ultrasound diagnostic apparatus and method for measuring a size of a target object
JP2014525328A (en) How to detect and track a needle
KR101175426B1 (en) Ultrasound system and method for providing three-dimensional ultrasound image
JP4879263B2 (en) Ultrasonic diagnostic apparatus and ultrasonic diagnostic method
JP2001286474A (en) Dynamic measurement of subject's parameter
CN107690312B (en) Ultrasonic imaging apparatus
JP4963278B2 (en) Boundary extraction by weighted least squares method using separation measure
JP5642997B2 (en) Ultrasound system and method for providing multiple slice images
JP5523019B2 (en) Ultrasonic system for forming elastic image and elastic image forming method
JP2011031023A (en) Ultrasonic system and method for setting sagittal view
CN109452954B (en) Ultrasonic imaging method and device
JP2010119847A (en) Ultrasound system and method of providing standard view information
JP4634872B2 (en) Ultrasonic diagnostic equipment
JP5826527B2 (en) Ultrasound system and method for providing elastic images
KR20120046539A (en) Ultrasound system and method for providing body mark
JP6761365B2 (en) Ultrasound diagnostic equipment and programs
KR20090013199A (en) Geometry Analysis of 3D Echocardiography
JP6731369B2 (en) Ultrasonic diagnostic device and program
JP2007518508A (en) Ultrasound imaging system for calculating compressibility

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20080318

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20100824

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20101019

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20101116

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20101119

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131126

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees