JP5406499B2 - Chest X-ray image conversion method, apparatus, system, computer program, and medium - Google Patents
Chest X-ray image conversion method, apparatus, system, computer program, and medium Download PDFInfo
- Publication number
- JP5406499B2 JP5406499B2 JP2008255980A JP2008255980A JP5406499B2 JP 5406499 B2 JP5406499 B2 JP 5406499B2 JP 2008255980 A JP2008255980 A JP 2008255980A JP 2008255980 A JP2008255980 A JP 2008255980A JP 5406499 B2 JP5406499 B2 JP 5406499B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- chest
- ray
- ray image
- template
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 85
- 238000006243 chemical reaction Methods 0.000 title claims description 25
- 238000004590 computer program Methods 0.000 title claims description 11
- 238000012545 processing Methods 0.000 claims description 37
- 238000011976 chest X-ray Methods 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 19
- 238000000605 extraction Methods 0.000 claims description 16
- 238000009966 trimming Methods 0.000 claims description 11
- 238000003384 imaging method Methods 0.000 claims description 9
- 230000002708 enhancing effect Effects 0.000 claims description 5
- 238000003702 image correction Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 description 5
- 230000000877 morphologic effect Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000013500 data storage Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 230000004931 aggregating effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Description
本願発明は、フィルムレス化により映像がディジタル化された医療用デジタルX線画像の処理方法及び同画像の処理システム等に関する。特にデジタル胸部X線画像の縦位置、横位置あるいは左右を判別して適切な変換を行う胸部正面X線画像の変換方法、装置等に関する。 The present invention relates to a medical digital X-ray image processing method in which an image is digitized by filmless, a processing system for the image, and the like. In particular, the present invention relates to a chest front X-ray image conversion method, apparatus, and the like that perform appropriate conversion by determining the vertical position, horizontal position, or left and right of a digital chest X-ray image.
医療用X線画像撮影のフィルムレス化、デジタル化に伴い、画像の確認作業はモニター上での検像システムを用いた作業へ移行しつつある。画像の確認作業とは、医師がX線写真を利用して診断を行うに先だって、X線写真の映像が適正に撮影されているかを確認する作業である。この作業が必要なのは、人体の影響を考慮してX線線量をできるだけ少なくして撮影するため、ノイズが多いなど診断に不適切な映像を得る場合があるからである。また、何らかの理由によって画像の上下あるいは縦横、さらには前後が逆転してしまうこともある。フィルムを利用していたときはもちろん、モニター上でデジタル画像の確認を行うようになっても、検像作業は放射線技師自身がX線写真映像を自分の目で確認しているのが現状である。能率が良くないことはもちろん、フィルムレス化で映像がディジタル化されたメリットを十分享受していないと言える。
胸部X線画像については、画像の上下、縦横が誤って逆になったり、左右(あるいは表裏)が逆になったりすることも起こりえるが、このような画像の逆転を生じた場合、画像解析を行うX線技師や放射線科医師は、当然正常な向きに画像表示を変換した上で解析を行わなければならず、その作業には一定の時間と労力とを要する。 For chest X-ray images, the image may be reversed upside down, vertically and horizontally, or the left and right (or front and back) may be reversed. The X-ray engineer or radiologist performing the analysis must naturally perform the analysis after converting the image display in the normal direction, and the work requires a certain amount of time and labor.
このような問題に対して、本願発明は、X線撮影装置がディジタル化された利点をさらに活かすため、ディジタル化された医療用X線画像、とりわけ胸部X線画像の縦横、左右および前後逆転の有無を判別し、これらの逆転が起きていると判別された場合には正しい画像になるようにデータの変換を行なうものである。このような判別と変換を自動的に行うことによって、胸部X線画像を評価するX線技師および放射線科医師等は、常に縦横および左右の統一された画像に接することができ、作業の能率化と信頼性の向上を図ることができる。 In order to further take advantage of the digitalization of the X-ray imaging apparatus, the present invention addresses such problems, so that the digitized medical X-ray image, in particular, the chest X-ray image can be vertically / horizontally, left / right, and front / rear reversed. The presence / absence is determined, and when it is determined that the reverse is occurring, the data is converted so that a correct image is obtained. By automatically performing such discrimination and conversion, X-ray technicians and radiologists who evaluate chest X-ray images can always come into contact with unified images in the vertical and horizontal directions and the left and right, thereby improving work efficiency. And reliability can be improved.
上記の問題を解決するために、本発明は、行列配置された画素情報からなるデジタル胸部正面X線画像から対象領域を抽出する抽出ステップと、
前記画像の方向を判別する判別ステップと、
判別された前記画像の方向が所定の方向と異なる場合に、前記画像の方向を変換して所定方向の画像を得る画像修正ステップとを有する胸部正面X線画像の変換方法を提案する。
In order to solve the above problem, the present invention includes an extraction step of extracting a target region from a digital chest front X-ray image composed of pixel information arranged in a matrix,
A determination step of determining a direction of the image;
A method of converting a chest front X-ray image comprising an image correcting step of converting the direction of the image to obtain an image in a predetermined direction when the determined direction of the image is different from a predetermined direction is proposed.
本明細書において行列配置された画素情報とは、横方向(X方向)の座標値と縦方向(Y方向)の座標値に対して画像を構成する画素(ピクセル)の値が定義される形式の画像データの意味で用いる。画素情報とは、一般には白から黒に至る濃淡を示すデジタルデータであるが、デジタルデータは必ずしも濃淡に対応するものである必要はなく、透過放射線強度あるいはそれ以外のデータに対応するものであっても良い。対象領域とは、以後の処理の対象となる画像領域を指す。 The pixel information arranged in a matrix in this specification is a format in which values of pixels (pixels) constituting an image are defined with respect to coordinate values in the horizontal direction (X direction) and coordinate values in the vertical direction (Y direction). This is used to mean the image data. The pixel information is generally digital data indicating the shade from white to black, but the digital data does not necessarily correspond to the shade, and corresponds to the transmitted radiation intensity or other data. May be. The target area refers to an image area to be processed later.
画像の方向判別は、上下の逆転、縦向き横向きの違い、左右(裏表)の逆転が生じているか否かを判断することをいう。判別の結果、前記いずれかの逆転が生じている場合は画像の方向を変換して正しい向きの画像を得るが、何れの逆転も生じていない場合は、元の画素情報が正しい向きの画像であることはいうまでもない。上記の方法によって画像の変換方向を決定した後は、当該変換方向に従って、もとの画像の方向変換を行う。 Image orientation determination refers to determining whether vertical reversal, vertical / horizontal difference, and left / right (back and front) reversal have occurred. As a result of the determination, if any of the above inversions occurs, the direction of the image is converted to obtain an image with the correct orientation, but if no inversion occurs, the original pixel information is the image with the correct orientation. Needless to say. After the image conversion direction is determined by the above method, the original image is converted in accordance with the conversion direction.
上記の画像の変換方法において、抽出ステップは、
画像の濃淡を強調するコントラスト調整過程と、
画像を白または黒で表示する2値化過程と、
ノイズ除去過程と、
背景除去過程と、
トリミング過程との内の少なくとも1つを含むものであることが望ましい。
In the image conversion method, the extraction step includes:
A contrast adjustment process that enhances the contrast of the image,
A binarization process for displaying an image in white or black;
Noise removal process,
Background removal process,
It is desirable to include at least one of the trimming processes.
コントラスト調整には、ガンマ補正、階調変更、モルフォロジカルフィルタ、ボケマスク処理などを適宜用いることができる。ここで、ガンマ補正とは、白又は黒に近い画素値を一層白又は黒に近づける作業、あるいは、白又は黒の値で置き換える作業をいう。階調変更とは、例えば、256段階の階調で表現されているデータを32段階の階調に置き換える(近似する)ことをいう。モルフォロジカルフィルタとはモルフォロジー演算に基づくフィルター処理、ボケマスク処理とは、平滑化フィルター処理をさす。 For contrast adjustment, gamma correction, gradation change, morphological filter, blur mask processing, and the like can be used as appropriate. Here, the gamma correction refers to an operation of making a pixel value close to white or black closer to white or black, or an operation of replacing with a white or black value. Tone change means, for example, replacing (approximate) data expressed in 256 levels of gradation with 32 levels of gradation. The morphological filter means filter processing based on morphological operation, and the blur mask processing means smoothing filter processing.
2値化処理は、4段階(2ビット)以上の階調で表現されている画素データを、2段階(1ビット)データに圧縮する作業であり、基本的には白と黒の中間に一定の閾値を定めて当該閾値より黒よりの値を黒に、閾値よりも白よりの値を白に置き換えるものである。ノイズ除去処理は、特に2値化によって生じた固有のノイズを除去するために実施するもので、種々の技法が存在する。収縮処理(値0を有する画素に隣接する8つの画素の値を0とする処理)と膨張処理(値1を有する画素に隣接する8つの画素の値を1とする処理)の連続によってノイズ除去を行うことができるがこれらに限定されるわけではない。背景除去処理は、例えば、抽出対象の領域内に領域拡張回始点を設定し、その画素と連結したが疎である条件を満たすが疎を1つの領域として抽出する領域拡張法などがあるがこれに限定されない。トリミングは、不要な周辺部位を削除して解析対象を絞り込む操作である。 The binarization process is an operation to compress pixel data expressed in gradations of 4 steps (2 bits) or more into 2 steps (1 bit) data. Basically, it is constant between white and black. And a value from black is replaced with black, and a value from white is replaced with white. The noise removal processing is performed particularly for removing inherent noise caused by binarization, and various techniques exist. Noise removal by continuous shrinkage processing (processing that sets the value of 8 pixels adjacent to a pixel having a value of 0 to 0) and expansion processing (processing that sets the value of 8 pixels that are adjacent to a value of 1 to 1) However, the present invention is not limited to these. Background removal processing includes, for example, a region expansion method in which a region expansion start point is set in a region to be extracted and connected to the pixel but satisfies a sparse condition but sparse is extracted as one region. It is not limited to. Trimming is an operation of narrowing down the analysis target by deleting unnecessary peripheral parts.
ここで、前記判別ステップは、
画像を構成する行列の列ごとに黒画素の数を集計する過程と、
1以上の黒画素が存在する2つの列の間に黒がその存在しない列があれば前記画像は縦画像であると判別する縦横判別過程を含むのが好ましい。
Here, the determination step includes:
A process of counting the number of black pixels for each column of the matrix constituting the image;
It is preferable to include a vertical / horizontal discrimination process for discriminating that the image is a vertical image if there is a column where black does not exist between two columns where one or more black pixels exist.
当該判別ステップは、さらに、
胸部正面X線画像テンプレートおよび胸部背面X線画像テンプレートと前記画像の相関係数を算出し、
1)胸部正面X線画像テンプレートとの相関が第1の所定値以上であり、胸部背面X線画像テンプレートとの相関が第2の所定値以下であれば前記画像は正面画像であると判別し、
2)胸部正面X線画像テンプレートとの相関が第2の所定値以下であり、胸部背面X線画像テンプレートとの相関が第1の所定値以上であれば前記画像は背面画像であると判別し、
3)上記1)と2)に何れにも当てはまらなければ前記画像が正面画像か背面画像かは不明であると判別する前後判別過程を含んでもよい。
The determination step further includes
Calculating a correlation coefficient between the chest front X-ray image template and the chest back X-ray image template and the image;
1) If the correlation with the chest front X-ray image template is not less than the first predetermined value and the correlation with the chest back X-ray image template is not more than the second predetermined value, it is determined that the image is a front image. ,
2) If the correlation with the chest front X-ray image template is less than or equal to the second predetermined value and the correlation with the chest rear X-ray image template is greater than or equal to the first predetermined value, it is determined that the image is a rear image. ,
3) If neither of the above 1) and 2) applies, it may include a before-and-after discriminating process for discriminating whether the image is a front image or a back image.
前記胸部正面X線画像テンプレートとは、X線画像の正面及び背面の判別に用いるための代表的な正面画像と背面画像である。最大値が1に、最小値が0に基準化された状態で、例えば、第1の所定値は0.7、第2の所定値とは0.3という値である。第1及び第2の所定値は、テンプレートとの関係等を考慮して適宜選択することができる。前記1)と2)の何れにも当てはまらない場合は、その旨を出力して作業者の注意を喚起するのが好ましい。さらに、当該方法は、画像の逆転等を発見して変換が行われたか否かを表記して作業者の注意を喚起するのが一層好ましい。 The chest front X-ray image template is a typical front image and back image to be used for discrimination of the front and back of an X-ray image. In a state where the maximum value is normalized to 1 and the minimum value is normalized to 0, for example, the first predetermined value is 0.7 and the second predetermined value is 0.3. The first and second predetermined values can be appropriately selected in consideration of the relationship with the template and the like. If none of the above 1) and 2) apply, it is preferable to output that fact and alert the operator. Furthermore, it is more preferable that the method calls the operator's attention by indicating whether or not the conversion has been performed by finding the reverse of the image.
本発明はさらに、行列配置された画素情報からなるデジタル胸部正面X線画像から対象領域を抽出する抽出ステップと、
前記画像の方向を判別する判別ステップと、
判別された前記画像の方向が所定の方向と異なる場合に、前記画像の方向を変換して所定方向の画像を得る画像修正ステップとをコンピュータに実行させる胸部正面X線画像の変換用コンピュータプログラムを提供する。
The present invention further includes an extraction step of extracting a target region from a digital chest front X-ray image composed of pixel information arranged in a matrix;
A determination step of determining a direction of the image;
A computer program for converting a chest front X-ray image that causes a computer to execute an image correction step of converting the direction of the image to obtain an image in a predetermined direction when the determined direction of the image is different from a predetermined direction provide.
ここで、コンピュータは、データ演算処理能力を有する装置全体を指しており、いわゆるパソコンあるいはサーバと称する独立したコンピュータ及び、画像表示装置あるいはX線画像撮像装置と一体となった機器を含む用語として用いる。 Here, the computer refers to the entire apparatus having the data calculation processing capability, and is used as a term including an independent computer called a personal computer or a server and an apparatus integrated with an image display apparatus or an X-ray imaging apparatus. .
前記コンピュータプログラムが実行する抽出ステップは、
前記画像の濃淡を強調するコントラスト調整過程と、
前記画像を白または黒で表示する2値化過程と、
ノイズ除去過程と、
背景除去過程と、
トリミング過程との内の少なくとも1つを含み、
前記判別ステップは、
前記画像を構成する行列の列ごとに黒画素の数を集計する過程と、
1以上の黒画素が存在する2つの列の間に黒がその存在しない列があれば前記画像は縦画像であると判別する縦横判別過程と、
胸部正面X線画像テンプレートおよび胸部背面X線画像テンプレートと前記画像の相関係数を算出し、
1)胸部正面X線画像テンプレートとの相関が第1の所定値以上であり、胸部背面X線画像テンプレートとの相関が第2の所定値以下であれば前記画像は正面画像であると判別し、
2)胸部正面X線画像テンプレートとの相関が第2の所定値以下であり、胸部背面X線画像テンプレートとの相関が第1の所定値以上であれば前記画像は背面画像であると判別し、
3)上記1)と2)に何れにも当てはまらなければ前記画像が正面画像か背面画像かは不明であると判別する前後判別過程を含む、請求項5に記載の胸部正面X線画像の変換用コンピュータプログラム。
The extraction step executed by the computer program includes:
A contrast adjustment process for enhancing the shade of the image;
A binarization process for displaying the image in white or black;
Noise removal process,
Background removal process,
Including at least one of trimming processes;
The determination step includes
Totalizing the number of black pixels for each column of the matrix constituting the image;
A vertical / horizontal discrimination process for discriminating that the image is a vertical image if there is a column in which black does not exist between two columns in which one or more black pixels exist;
Calculating a correlation coefficient between the chest front X-ray image template and the chest back X-ray image template and the image;
1) If the correlation with the chest front X-ray image template is not less than the first predetermined value and the correlation with the chest back X-ray image template is not more than the second predetermined value, it is determined that the image is a front image. ,
2) If the correlation with the chest front X-ray image template is less than or equal to the second predetermined value and the correlation with the chest rear X-ray image template is greater than or equal to the first predetermined value, it is determined that the image is a rear image. ,
6) The conversion of the chest front X-ray image according to claim 5, further comprising a before-and-after determination process for determining whether the image is a front image or a back image if neither of the above 1) and 2) applies. Computer program.
上記プログラムは、変換が行われた場合は変換の結果得られた画像を、変換が行われなかった場合は元の画像を記憶し、また必要に応じて表示する機能を有するのが好ましい。 The program preferably has a function of storing an image obtained as a result of conversion when conversion is performed, or storing an original image when conversion is not performed, and displaying the image as necessary.
本発明はまた、前記胸部正面X線画像の変換用コンピュータプログラムを格納したコンピュータによって読み取り可能な記憶媒体を提供する。 The present invention also provides a computer-readable storage medium storing a computer program for converting the chest front X-ray image.
本発明はさらに、行列配置された画素情報からなるデジタル胸部正面X線画像を記憶する記憶媒体と、
前記画像から対象領域を抽出する抽出ステップと、前記画像の方向を判別する判別ステップと、判別された前記画像の方向が所定の方向と異なる場合に、前記画像の方向を変換して所定方向の画像を得る画像修正ステップとを実行する演算処理装置と、
前記演算処理装置が得た画像を記憶する記憶媒体とを有する胸部正面X線画像の変換装置を提供する。
The present invention further includes a storage medium for storing a digital chest front X-ray image comprising pixel information arranged in a matrix,
An extraction step for extracting a target region from the image; a determination step for determining the direction of the image; and if the determined direction of the image is different from a predetermined direction, the direction of the image is converted to a predetermined direction An arithmetic processing unit that executes an image correction step for obtaining an image;
A chest front X-ray image conversion device having a storage medium for storing an image obtained by the arithmetic processing device is provided.
前記変換装置は、全体として1つの機械的集合として捉えられる形態を有している装置と、分散しているために複数の部分の単なる集合の形態を有している装置とを含む。 The conversion device includes a device having a form that can be regarded as one mechanical set as a whole, and a device having a form of a simple set of a plurality of parts because they are distributed.
前記装置における抽出ステップは、前記方法及びコンピュータプログラムにおけるのと同様、
画像の濃淡を強調するコントラスト調整過程と、
前記画像を白または黒で表示する2値化過程と、
ノイズ除去過程と、
背景除去過程と、
トリミング過程との内の少なくとも1つを含み、
前記判別ステップは、
前記画像を構成する行列の列ごとに黒画素の数を集計する過程と、
1以上の黒画素が存在する2つの列の間に黒がその存在しない列があれば前記画像は縦画像であると判別する縦横判別過程と、
胸部正面X線画像テンプレートおよび胸部背面X線画像テンプレートと前記画像の相関係数を算出し、
1)胸部正面X線画像テンプレートとの相関が第1の所定値以上であり、胸部背面X線画像テンプレートとの相関が第2の所定値以下であれば前記画像は正面画像であると判別し、
2)胸部正面X線画像テンプレートとの相関が第2の所定値以下であり、胸部背面X線画像テンプレートとの相関が第1の所定値以上であれば前記画像は背面画像であると判別し、
3)上記1)と2)に何れにも当てはまらなければ前記画像が正面画像か背面画像かは不明であると判別する前後判別過程を含むのが好ましい。
The extraction step in the apparatus is the same as in the method and computer program.
A contrast adjustment process that enhances the contrast of the image,
A binarization process for displaying the image in white or black;
Noise removal process,
Background removal process,
Including at least one of trimming processes;
The determination step includes
Totalizing the number of black pixels for each column of the matrix constituting the image;
A vertical / horizontal discrimination process for discriminating that the image is a vertical image if there is a column in which black does not exist between two columns in which one or more black pixels exist;
Calculating a correlation coefficient between the chest front X-ray image template and the chest back X-ray image template and the image;
1) If the correlation with the chest front X-ray image template is not less than the first predetermined value and the correlation with the chest back X-ray image template is not more than the second predetermined value, it is determined that the image is a front image. ,
2) If the correlation with the chest front X-ray image template is less than or equal to the second predetermined value and the correlation with the chest rear X-ray image template is greater than or equal to the first predetermined value, it is determined that the image is a rear image. ,
3) It is preferable to include a before-and-after discriminating process for discriminating whether the image is a front image or a back image if neither of the above 1) and 2) applies.
本発明はさらに、行列配置された画素情報からなるデジタル胸部正面X線画像撮影装置と、
前記胸部正面X線画像の変換装置と、
前記演算処理装置が得た画像を表示する表示装置とを備えた胸部正面X線画像処理システムを提供する。
The present invention further includes a digital chest front X-ray imaging apparatus comprising pixel information arranged in a matrix,
The chest front X-ray image converter;
A chest front X-ray image processing system including a display device for displaying an image obtained by the arithmetic processing device is provided.
以下に実施例に基づいて本発明を説明するが、実施例は本発明の理解の助けのために記載するものであって、本発明が以下に記載する実施例に限定されるものではないことは自明である。 The present invention will be described below based on examples, but the examples are described for the purpose of helping understanding of the present invention, and the present invention is not limited to the examples described below. Is self-explanatory.
抽出ステップ
まず、実施例に基づいて、行列配置された画素情報からなるデジタル胸部正面X線画像から対象領域を抽出する抽出ステップについて以下に説明する。実施例によれば、抽出ステップは、画像の濃淡を強調するコントラスト調整過程と、画像を白または黒で表示する2値化過程と、ノイズ除去過程と、背景除去過程と、トリミング過程とを含む。ただし、これら全ての過程を含むことは必ずしも必要ではなく、以下に記載する過程の少なくとも1つを含むことで足りる場合もある。
Extraction Step First, an extraction step for extracting a target region from a digital chest front X-ray image made up of pixel information arranged in a matrix will be described based on an embodiment. According to the embodiment, the extraction step includes a contrast adjustment process for enhancing the shading of the image, a binarization process for displaying the image in white or black, a noise removal process, a background removal process, and a trimming process. . However, it is not always necessary to include all these processes, and it may be sufficient to include at least one of the processes described below.
図1(A)は、上下が逆転した胸部X線画像の例、図1(B)は縦横を取り違えた胸部X線画像の例である。これらに例示される、上下の転倒、縦横の取り違えおよび左右(表裏)の間違いを自動的に発見し、これらを修正して正しい胸部X線画像を提供することが本発明の目的である。 FIG. 1A is an example of a chest X-ray image with the top and bottom reversed, and FIG. 1B is an example of a chest X-ray image with the vertical and horizontal reversed. It is an object of the present invention to automatically detect upside down, upside down, left and right (front and back) mistakes, and correct these to provide a correct chest X-ray image.
作業の流れとしては、図2に示すように、胸部X線画像撮像装置または胸部X線画像データ読み出し装置10から画像データの送信を受け、演算装置100によって本発明によるコンピュータプログラムを実行し、画像表示装置または画像データ記憶装置20に変換画像または変換が不要な場合は元の画像を記憶させることになる。
As shown in FIG. 2, the work flow includes transmission of image data from the chest X-ray image capturing device or the chest X-ray image
図3を参照して、受信した画像について行うコントラスト調整等の各処理について以下に説明する。受信画像については、まずコントラスト調整110を行う。コントラスト調整は、前記のように、ガンマ補正、階調変更、モルフォロジカルフィルタ、ボケマスク処理の何れかあるいは複数を併用することができる。図3の110pがコントラスト調整を行った胸部X線画像の例である。コントラスト調整の結果、白と黒の間の所定の範囲に値が分散したデジタル画像データを得ることができる。
With reference to FIG. 3, each process such as contrast adjustment performed on the received image will be described below.
次に2値化処理120を行う。2値化処理の1例は判断分析法である。判断分析法とは、画像の濃度に対応する値のヒストグラムを求め、濃度の山が2つになった場合、片方の山を背景、他方の山を対象物と考え、2つの山の間を閾値とする方法である。例えば、このようにして決定した閾値より大きい値を有する画素の値を「1」(白色)に置換し、閾値よりも小さい画素の値を「0」(黒色)で置換することで、全てのがその値を1または0に置き換える作業である。2値化処理を行った結果の画像を120pに示す。
Next,
次に、2値化処理によって生じたノイズを除去するノイズ除去処理またはオープニング処理130を行う。ノイズ除去処理は、例えば、値0を有する画素に隣接する8つの画素の値を全て0とする収縮処理によって孤立点や小さな飛び出しを除去し、さらに、値1を有する画素に隣接する8つのがその値を全て1とする膨張処理によって小さな穴やくぼみを埋める処理である。130pはノイズ除去処理後の画像の例である。
Next, a noise removal process or an
次に、背景除去処理140を行う。背景除去処理140は、抽出対象の領域内に領域拡張回始点を設定し、その画素と連結したが疎である条件を満たすが疎を1つの領域として抽出する領域拡張方によることができる。その結果、背景を取り除いた対象領域を取り出すことが可能になる。140pは背景除去処理後の画像を例示するものである。
Next,
さらに、トリミング処理150を行う。トリミングは、任意の領域から肺領域の上下左右の最高点、最下点、最右点、最左点を決定し、当該点よりも外側を除去することによって行うことができる。150pはトリミング後の画像である。
Further, trimming
判別ステップ
次に、画像の方向を判別する判別ステップについて、実施例に基づいて以下に記載する。実施例の判別ステップは、画像を構成する行列の列ごとに黒画素の数を集計する過程と、1以上の黒画素が存在する2つの列の間に黒がその存在しない列があれば前記画像は縦画像であると判別する縦横判別過程、および、胸部正面X線画像テンプレートおよび胸部背面X線画像テンプレートと前記画像の相関係数を算出し、相関係数に基づいて画像の前後を判別する前後判別過程を含むが、縦横判別過程と前後判別過程の両方を含むことは必ずしも必須ではなく、これら一方を含むものであっても良い。
Determination Step Next, the determination step for determining the direction of the image will be described below based on the embodiment. The determination step of the embodiment includes the process of counting the number of black pixels for each column of the matrix constituting the image, and if there is a column where black does not exist between two columns where one or more black pixels exist Vertical / horizontal discrimination process for discriminating that the image is a vertical image, and the correlation coefficient between the chest front X-ray image template and the chest back X-ray image template and the image is calculated, and the front and back of the image are determined based on the correlation coefficient However, it is not always essential to include both the longitudinal and lateral discrimination process and the front and back discrimination process, and one of these may be included.
実施例に基づけば、図5に示す判別ステップ200は、
画像を構成する行列の列ごとに黒画素の数を集計する過程と、
1以上の黒画素が存在する2つの列の間に黒がその存在しない列があれば前記画像は縦画像であると判別する。
Based on the embodiment, the determination step 200 shown in FIG.
A process of counting the number of black pixels for each column of the matrix constituting the image;
If there is a column in which black does not exist between two columns in which one or more black pixels exist, it is determined that the image is a vertical image.
図4および図5の200p2は、トリミング後の画像200p1について、画像を構成する行列の列ごとに黒画素の数を集計した結果をグラフ表示したものである。図4の横軸は画像150pの横軸の対応し、縦軸は、その横軸の位置における白画素の数を表す。グラフは、左端から右に向かって、大きな値(白画素が多い)から徐々に下降し極小地をとった後、中央付近で最大値となり(全部白画素である)、再度極小値をとった後に、右端で大きな値にいたる曲線である。
200p2 in FIG. 4 and FIG. 5 is a graph display of the result of aggregating the number of black pixels for each column of the matrix constituting the image for the trimmed image 200p1. The horizontal axis in FIG. 4 corresponds to the horizontal axis of the
このグラフの特徴は、横軸の中央付近で最大値を生じている点、換言すれば、画像の左右中央位置を上下に切った位置では全ての画素が白であることを示している点である。画像の縦横が逆になっている(つまり画像が寝ている)状態では、どの位置においても上下方向の全ての画素が白になることはない。したがって、本実施例によれば、1以上の黒画素が存在する2つの列の間に黒がその存在しない列があれば前記画像は縦画像であると判別する。 The characteristic of this graph is that the maximum value is generated near the center of the horizontal axis, in other words, that all pixels are white at the position where the left and right center position of the image is cut up and down. is there. In a state where the vertical and horizontal directions of the image are reversed (that is, the image is lying down), all the pixels in the vertical direction are not white at any position. Therefore, according to the present embodiment, if there is a column in which black does not exist between two columns in which one or more black pixels exist, it is determined that the image is a vertical image.
また、テンプレートとの相関を用いて画像の上下および縦横を同時に判別することも可能である。
すなわち、図5に300p1〜300p4に例示するような、上下および縦横の向きがそれぞれことなる4種のテンプレートを準備し、評価対象となる胸部X線画像がテンプレートのうちのどれと最も相関が高いかに基づいて評価対象画像の向きを判断するものである。相関は、以下の式1によって算出することができる。
That is, four types of templates having different vertical and horizontal orientations as exemplified by 300p1 to 300p4 in FIG. 5 are prepared, and the chest X-ray image to be evaluated has the highest correlation with any of the templates. The direction of the evaluation target image is determined based on how. The correlation can be calculated by the following equation 1.
rが1に近い場合(数学的には−1に近い場合も)、評価対象画像はテンプレートと強い相関を有志、rが0に近い場合はテンプレートとの相関は弱い。通常は、前記4つのテンプレートのうちの1つのみと強い層間を示し、他の3つのテンプレートとの相関は弱い。いずれのテンプレートとも層間が十分強くない場合、あるいは複数のテンプレートと強い層間を示す場合は、注意を喚起するメッセージと共に評価者の判断を求めるのが好ましい。 When r is close to 1 (also mathematically close to -1), the evaluation target image has a strong correlation with the template, and when r is close to 0, the correlation with the template is weak. Usually, only one of the four templates and a strong interlayer are shown, and the correlation with the other three templates is weak. When the layers are not sufficiently strong in any template, or when showing a strong layer between a plurality of templates, it is preferable to request the evaluator's judgment together with a message for calling attention.
上記の手順によって、画像の上下あるいは縦横を判断を行った結果、対象画像が標準的な向き(正立で人体を正面から見た位置関係)と異なっている場合は、標準的な向きになるように画像の座標変換を行う。座標変換の方法自体は、純粋に数学的な処理であり、格別の困難性もないのでここでは割愛する。 If the target image is different from the standard orientation (position that is upright and the human body is viewed from the front) as a result of determining the top / bottom or length / width of the image by the above procedure, the orientation is standard. In this way, image coordinate conversion is performed. The coordinate transformation method itself is purely mathematical processing and has no particular difficulty, so it is omitted here.
次に、画像の前後判別を行い、必要な場合には画像の座標変換を行う。そのためには、 前記と同様の手法に従い、胸部正面X線画像テンプレートおよび胸部背面X線画像テンプレートと前記画像の相関係数を算出し、
1)胸部正面X線画像テンプレートとの相関が第1の所定値以上であり、胸部背面X線画像テンプレートとの相関が第2の所定値以下であれば前記画像は正面画像であると判別し、
2)胸部正面X線画像テンプレートとの相関が第2の所定値以下であり、胸部背面X線画像テンプレートとの相関が第1の所定値以上であれば前記画像は背面画像であると判別し、
3)上記1)と2)に何れにも当てはまらなければ前記画像が正面画像か背面画像かは不明であると判別する。
Next, an image front-rear discrimination is performed, and if necessary, image coordinate conversion is performed. For this purpose, according to the same method as described above, a correlation coefficient between the chest front X-ray image template and the chest back X-ray image template and the image is calculated,
1) If the correlation with the chest front X-ray image template is not less than the first predetermined value and the correlation with the chest back X-ray image template is not more than the second predetermined value, it is determined that the image is a front image. ,
2) If the correlation with the chest front X-ray image template is less than or equal to the second predetermined value and the correlation with the chest rear X-ray image template is greater than or equal to the first predetermined value, it is determined that the image is a rear image. ,
3) If neither of the above 1) and 2) applies, it is determined that it is unknown whether the image is a front image or a back image.
画像の前後判断を行った結果、対象画像が標準的な向き(人体を正面から見た位置関係)と異なっている場合は、標準的な向きになるように画像の座標変換を行う点は前出の上下および縦横判断と同様である。 If the target image is different from the standard orientation (positional relationship when the human body is viewed from the front) as a result of judging the front and back of the image, the point that the image is transformed to the standard orientation is This is the same as determining the up / down and vertical / horizontal directions.
10 胸部X線画像データ読み出し装置
20 画像データ記憶装置
100 演算装置
110p コントラスト調整を行った胸部X線画像の例
120 2値化処理
120p 2値化処理を行った結果の画像
130 オープニング処理
130p ノイズ除去処理後の画像の例
140 背景除去処理
140p 背景除去処理後の画像の例
DESCRIPTION OF
Claims (9)
前記画像の方向を判別する判別ステップと、
判別された前記画像の方向が所定の方向と異なる場合に、前記画像の方向を変換して所定方向の画像を得る画像修正ステップとを有する胸部正面X線画像の変換方法であって、
前記判別ステップは、
胸部正面X線画像テンプレートおよび胸部背面X線画像テンプレートと前記画像の相関係数を算出し、
1)胸部正面X線画像テンプレートとの相関が第1の所定値以上であり、胸部背面X線画像テンプレートとの相関が第2の所定値以下であれば前記画像は正面画像であると判別し、
2)胸部正面X線画像テンプレートとの相関が第2の所定値以下であり、胸部背面X線画像テンプレートとの相関が第1の所定値以上であれば前記画像は背面画像であると判別し、
3)上記1)と2)に何れにも当てはまらなければ前記画像が正面画像か背面画像かは不明であると判別する前後判別過程を含む胸部正面X線画像の変換方法。 An extraction step of extracting a target region from a digital chest X-ray image comprising pixel information arranged in a matrix;
A determination step of determining a direction of the image;
An image correction step of converting a direction of the image to obtain an image in a predetermined direction when the determined direction of the image is different from a predetermined direction ,
The determination step includes
Calculating a correlation coefficient between the chest front X-ray image template and the chest back X-ray image template and the image;
1) If the correlation with the chest front X-ray image template is not less than the first predetermined value and the correlation with the chest back X-ray image template is not more than the second predetermined value, it is determined that the image is a front image. ,
2) If the correlation with the chest front X-ray image template is less than or equal to the second predetermined value and the correlation with the chest rear X-ray image template is greater than or equal to the first predetermined value, it is determined that the image is a rear image. ,
3) A method for converting a chest front X-ray image including a before-and-after determination process for determining whether the image is a front image or a back image if neither of the above 1) and 2) applies.
前記画像の濃淡を強調するコントラスト調整過程と、
前記画像を白または黒で表示する2値化過程と、
ノイズ除去過程と、
背景除去過程と、
トリミング過程とを含む請求項1に記載の胸部正面X線画像の変換方法。 The extraction step includes
A contrast adjustment process for enhancing the shade of the image;
A binarization process for displaying the image in white or black;
Noise removal process,
Background removal process,
The method of converting a chest front X-ray image according to claim 1 including a trimming process.
前記画像を構成する行列の列ごとに黒画素の数を集計する過程と、
1以上の黒画素が存在する2つの列の間に黒がその存在しない列があれば前記画像は縦画像であると判別する縦横判別過程を含む請求項1又は2に記載の胸部正面X線画像の変換方法。 The determination step includes
Totalizing the number of black pixels for each column of the matrix constituting the image;
The chest front X-ray according to claim 1 or 2, further comprising a vertical / horizontal discrimination process for discriminating that the image is a vertical image if there is a column in which black does not exist between two columns in which one or more black pixels exist. Image conversion method.
前記画像の方向を判別する判別ステップと、
判別された前記画像の方向が所定の方向と異なる場合に、前記画像の方向を変換して所定方向の画像を得る画像修正ステップとをコンピュータに実行させる胸部正面X線画像の変換用コンピュータプログラムであって、
前記判別ステップは、
胸部正面X線画像テンプレートおよび胸部背面X線画像テンプレートと前記画像の相関係数を算出し、
1)胸部正面X線画像テンプレートとの相関が第1の所定値以上であり、胸部背面X線画像テンプレートとの相関が第2の所定値以下であれば前記画像は正面画像であると判別し、
2)胸部正面X線画像テンプレートとの相関が第2の所定値以下であり、胸部背面X線画像テンプレートとの相関が第1の所定値以上であれば前記画像は背面画像であると判別し、
3)上記1)と2)に何れにも当てはまらなければ前記画像が正面画像か背面画像かは不明であると判別する前後判別過程を含む、胸部正面X線画像の変換用コンピュータプログラム。 An extraction step of extracting a target region from a digital chest X-ray image comprising pixel information arranged in a matrix;
A determination step of determining a direction of the image;
If the direction of the determined the image is different from the predetermined direction, the conversion computer program frontal chest X-ray image to execute an image correction step of obtaining an image of a predetermined direction by changing the direction of the image to the computer There,
The determination step includes
Calculating a correlation coefficient between the chest front X-ray image template and the chest back X-ray image template and the image;
1) If the correlation with the chest front X-ray image template is not less than the first predetermined value and the correlation with the chest back X-ray image template is not more than the second predetermined value, it is determined that the image is a front image. ,
2) If the correlation with the chest front X-ray image template is less than or equal to the second predetermined value and the correlation with the chest rear X-ray image template is greater than or equal to the first predetermined value, it is determined that the image is a rear image. ,
3) A computer program for conversion of a chest front X-ray image including a before-and-after determination process for determining whether the image is a front image or a back image if neither of the above 1) and 2) applies.
前記画像の濃淡を強調するコントラスト調整過程と、
前記画像を白または黒で表示する2値化過程と、
ノイズ除去過程と、
背景除去過程と、
トリミング過程との内の少なくとも1つを含み、
前記判別ステップは、
前記画像を構成する行列の列ごとに黒画素の数を集計する過程と、
1以上の黒画素が存在する2つの列の間に黒がその存在しない列があれば前記画像は縦画像であると判別する縦横判別過程と、を含む、請求項4に記載の胸部正面X線画像の変換用コンピュータプログラム。 The extraction step includes
A contrast adjustment process for enhancing the shade of the image;
A binarization process for displaying the image in white or black;
Noise removal process,
Background removal process,
Including at least one of trimming processes;
The determination step includes
Totalizing the number of black pixels for each column of the matrix constituting the image;
The chest front X according to claim 4 , further comprising: a vertical / horizontal discrimination process for discriminating that the image is a vertical image if there is a black non-existent column between two columns having one or more black pixels. Computer program for line image conversion.
前記画像から対象領域を抽出する抽出ステップと、前記画像の方向を判別する判別ステップと、判別された前記画像の方向が所定の方向と異なる場合に、前記画像の方向を変換して所定方向の画像を得る画像修正ステップとを実行する演算処理装置と、
前記演算処理装置が得た画像を記憶する記憶媒体とを有する胸部正面X線画像の変換装置であって、
前記判別ステップは、
胸部正面X線画像テンプレートおよび胸部背面X線画像テンプレートと前記画像の相関係数を算出し、
1)胸部正面X線画像テンプレートとの相関が第1の所定値以上であり、胸部背面X線画像テンプレートとの相関が第2の所定値以下であれば前記画像は正面画像であると判別し、
2)胸部正面X線画像テンプレートとの相関が第2の所定値以下であり、胸部背面X線画像テンプレートとの相関が第1の所定値以上であれば前記画像は背面画像であると判別し、
3)上記1)と2)に何れにも当てはまらなければ前記画像が正面画像か背面画像かは不明であると判別する前後判別過程を含む、胸部正面X線画像の変換装置。 A storage medium for storing a digital chest front X-ray image composed of pixel information arranged in a matrix;
An extraction step for extracting a target region from the image; a determination step for determining the direction of the image; and if the determined direction of the image is different from a predetermined direction, the direction of the image is converted to a predetermined direction An arithmetic processing unit that executes an image correction step for obtaining an image;
A chest front X-ray image conversion device having a storage medium for storing an image obtained by the arithmetic processing device ,
The determination step includes
Calculating a correlation coefficient between the chest front X-ray image template and the chest back X-ray image template and the image;
1) If the correlation with the chest front X-ray image template is not less than the first predetermined value and the correlation with the chest back X-ray image template is not more than the second predetermined value, it is determined that the image is a front image. ,
2) If the correlation with the chest front X-ray image template is less than or equal to the second predetermined value and the correlation with the chest rear X-ray image template is greater than or equal to the first predetermined value, it is determined that the image is a rear image. ,
3) An apparatus for converting a chest front X-ray image, which includes a before-and-after determination process for determining whether the image is a front image or a back image if neither of the above 1) and 2) applies.
前記画像の濃淡を強調するコントラスト調整過程と、
前記画像を白または黒で表示する2値化過程と、
ノイズ除去過程と、
背景除去過程と、
トリミング過程との内の少なくとも1つを含み、
前記判別ステップは、
前記画像を構成する行列の列ごとに黒画素の数を集計する過程と、
1以上の黒画素が存在する2つの列の間に黒がその存在しない列があれば前記画像は縦画像であると判別する縦横判別過程と、
を含む、請求項7に記載の胸部正面X線画像の変換装置。 The extraction step includes
A contrast adjustment process for enhancing the shade of the image;
A binarization process for displaying the image in white or black;
Noise removal process,
Background removal process,
Including at least one of trimming processes;
The determination step includes
Totalizing the number of black pixels for each column of the matrix constituting the image;
A vertical / horizontal discrimination process for discriminating that the image is a vertical image if there is a column in which black does not exist between two columns in which one or more black pixels exist ;
The chest front X-ray image conversion device according to claim 7 , comprising:
前記請求項7又は8に記載の胸部正面X線画像の変換装置と、
前記演算処理装置が得た画像を表示する表示装置とを備えた胸部正面X線画像処理システム。 A digital chest front X-ray imaging apparatus comprising pixel information arranged in a matrix;
The chest front X-ray image conversion device according to claim 7 or 8 ,
A chest front X-ray image processing system comprising: a display device that displays an image obtained by the arithmetic processing device.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008255980A JP5406499B2 (en) | 2008-10-01 | 2008-10-01 | Chest X-ray image conversion method, apparatus, system, computer program, and medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008255980A JP5406499B2 (en) | 2008-10-01 | 2008-10-01 | Chest X-ray image conversion method, apparatus, system, computer program, and medium |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010082296A JP2010082296A (en) | 2010-04-15 |
JP5406499B2 true JP5406499B2 (en) | 2014-02-05 |
Family
ID=42246830
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008255980A Active JP5406499B2 (en) | 2008-10-01 | 2008-10-01 | Chest X-ray image conversion method, apparatus, system, computer program, and medium |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5406499B2 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012086134A1 (en) * | 2010-12-21 | 2012-06-28 | パナソニック株式会社 | Optical tomography image acquisition device |
JP5286395B2 (en) * | 2011-09-12 | 2013-09-11 | 富士フイルムRiファーマ株式会社 | Receptor binding ability imaging program, recording medium and method |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06123B2 (en) * | 1985-07-26 | 1994-01-05 | コニカ株式会社 | Chest X-ray examination support device |
JPH07101446B2 (en) * | 1987-06-04 | 1995-11-01 | 富士写真フイルム株式会社 | Image rollover determination method |
JP2808773B2 (en) * | 1990-01-09 | 1998-10-08 | 株式会社日立製作所 | Automatic gradation conversion device |
JP2929030B2 (en) * | 1990-06-28 | 1999-08-03 | コニカ株式会社 | Digital radiographic image position determination device |
JP4143149B2 (en) * | 1997-11-20 | 2008-09-03 | キヤノン株式会社 | Region extraction apparatus, region extraction method, and computer-readable recording medium |
JP2000157519A (en) * | 1998-11-25 | 2000-06-13 | Konica Corp | Image processor |
JP2003220055A (en) * | 2001-11-20 | 2003-08-05 | Konica Corp | Characteristic quantity extracting method, object recognizing method, image processing apparatus and image processing program |
JP4708740B2 (en) * | 2004-06-08 | 2011-06-22 | キヤノン株式会社 | Image processing apparatus and image processing method |
-
2008
- 2008-10-01 JP JP2008255980A patent/JP5406499B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2010082296A (en) | 2010-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5039310B2 (en) | Cerebral hemorrhage segmentation device | |
US7231522B2 (en) | Information processing apparatus, method for controlling the same, and storage medium | |
EP1501048B1 (en) | Method of segmenting a radiographic image into diagnostically relevant and diagnostically irrelevant regions | |
WO2015076406A1 (en) | Device for assisting in diagnosis of osteoporosis | |
JPH03206572A (en) | Automatizing system for gradation conversion | |
CN101076811A (en) | Acceptability of automatic determining digital picture | |
CN108510489B (en) | Pneumoconiosis detection method and system based on deep learning | |
KR20200087427A (en) | The diagnostic method of lymph node metastasis in thyroid cancer using deep learning | |
CN111369557B (en) | Image processing method, device, computing equipment and storage medium | |
JP5406499B2 (en) | Chest X-ray image conversion method, apparatus, system, computer program, and medium | |
JP2011076224A (en) | Device, system, method, and program for information processing | |
WO2013136093A2 (en) | Image data storage and sharing | |
JP6935663B1 (en) | Oral mucosal disease diagnosis support system, method and program | |
KR20190067551A (en) | Method for contrast enhancement of medical diagnostic image using statistical estimation and mapping function modeling technique | |
US8131101B2 (en) | Image processing method and image processing apparatus | |
US10552970B2 (en) | Efficient guide filter for depth refinement | |
CN112565623A (en) | Dynamic image display system | |
JP2000067223A (en) | Method and device for preparing file, method and device for controlling file transfer and recording medium | |
CN111429364A (en) | Image restoration method and device | |
JP2006263224A (en) | Radiation image processing apparatus, the radiation image processing method, program, and computer readable media | |
KR102633823B1 (en) | Apparatus for discriminating medical image and method thereof | |
CN117974475B (en) | Focus image fusion method and system under four-dimensional ultrasonic endoscope observation | |
KR102026085B1 (en) | Method and apparatus for surface completion in ct data | |
JP4303061B2 (en) | Image processing apparatus and method, and program | |
JP4741289B2 (en) | Image processing apparatus and image processing method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20110921 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130709 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130828 |
|
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: 20131008 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20131101 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5406499 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |