JP2003310589A - Dynamic range compressing method for radiation picture - Google Patents
Dynamic range compressing method for radiation pictureInfo
- Publication number
- JP2003310589A JP2003310589A JP2003074985A JP2003074985A JP2003310589A JP 2003310589 A JP2003310589 A JP 2003310589A JP 2003074985 A JP2003074985 A JP 2003074985A JP 2003074985 A JP2003074985 A JP 2003074985A JP 2003310589 A JP2003310589 A JP 2003310589A
- Authority
- JP
- Japan
- Prior art keywords
- dynamic range
- compression
- sus
- image
- subject
- 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.)
- Granted
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000007906 compression Methods 0.000 claims abstract description 116
- 230000006835 compression Effects 0.000 claims abstract description 114
- 238000012545 processing Methods 0.000 claims abstract description 21
- 238000012937 correction Methods 0.000 claims description 29
- 238000010586 diagram Methods 0.000 description 15
- 238000006243 chemical reaction Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 11
- 238000003384 imaging method Methods 0.000 description 10
- 101000920618 Homo sapiens Transcription and mRNA export factor ENY2 Proteins 0.000 description 9
- 102100031954 Transcription and mRNA export factor ENY2 Human genes 0.000 description 9
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 8
- 230000007423 decrease Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 230000005284 excitation Effects 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 3
- 210000000038 chest Anatomy 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 3
- 238000012886 linear function Methods 0.000 description 3
- 238000004020 luminiscence type Methods 0.000 description 3
- 210000000988 bone and bone Anatomy 0.000 description 2
- 239000011248 coating agent Substances 0.000 description 2
- 238000000576 coating method Methods 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 210000001370 mediastinum Anatomy 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 101150018759 CG10 gene Proteins 0.000 description 1
- 101000795655 Canis lupus familiaris Thymic stromal cotransporter homolog Proteins 0.000 description 1
- 102100033040 Carbonic anhydrase 12 Human genes 0.000 description 1
- 101000867855 Homo sapiens Carbonic anhydrase 12 Proteins 0.000 description 1
- 101150048609 RR21 gene Proteins 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 210000000709 aorta Anatomy 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000007740 vapor deposition Methods 0.000 description 1
Landscapes
- Radiography Using Non-Light Waves (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
Description
【０００１】[0001]
【発明の属する技術分野】本発明は放射線画像のダイナ
ミックレンジ圧縮方法に関し、詳しくは、オリジナル画
像信号を処理して該オリジナル画像よりもダイナミック
レンジの狭い画像を担持する画像信号を得る放射線画像
のダイナミックレンジ圧縮方法の改善技術に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method for compressing a dynamic range of a radiation image, and more particularly, to a dynamic image of a radiation image for processing an original image signal to obtain an image signal carrying an image having a narrower dynamic range than the original image. The present invention relates to a technique for improving a range compression method.
【０００２】[0002]
【従来の技術】従来、放射線画像において、画像領域内
の微細構造の観察適正を確保しつつ、濃度域を圧縮する
方法として、例えば特許文献１に開示されるようなもの
がある。2. Description of the Related Art Conventionally, as a method for compressing a density region of a radiation image while ensuring proper observation of a fine structure in the image region, there is a method disclosed in, for example, Japanese Patent Application Laid-Open No. 2004-242242.
【０００３】前記特許文献１に開示される圧縮方法は、
各画素点に対応して該各画素点を含む所定マスク領域内
のオリジナル画像信号Ｓorgを平均化することで非鮮鋭
マスク信号（ボケマスク信号）Ｓusを求め、この非鮮鋭
マスク信号Ｓusが増大するにつれて単調減少する関数を
ｆ（Ｓus）としたときに、ダイナミックレンジが圧縮さ
れた処理済み画像信号Ｓprocを、Ｓproc＝Ｓorg＋ｆ
（Ｓus）として得るものである。The compression method disclosed in Patent Document 1 is
An unsharp mask signal (blurring mask signal) Sus is obtained by averaging the original image signal Sorg in a predetermined mask area including each pixel point corresponding to each pixel point, and as the unsharp mask signal Sus increases. When the function that monotonically decreases is f (Sus), the processed image signal Sproc whose dynamic range is compressed is Sproc = Sorg + f
(Sus).
【０００４】[0004]
【特許文献１】特開平３−２２２５７７号公報[Patent Document 1] Japanese Patent Application Laid-Open No. 3-222577
【０００５】[0005]
【発明が解決しようとする課題】しかし、従来のダイナ
ミックレンジ圧縮方法によると、被写体のダイナミック
レンジに関わらず、一定の圧縮度合いでダイナミックレ
ンジを圧縮していたため、例えば、標準よりも痩せた人
の胸部画像における低濃度部分の圧縮処理では過剰圧縮
となって、縦隔部分の濃度が高くなり、診断性能を低下
させることになり、逆に、太った人の場合には、圧縮が
不足して縦隔部分が白く潰れるという問題が生じる可能
性があった。However, according to the conventional dynamic range compression method, the dynamic range is compressed with a certain degree of compression regardless of the dynamic range of the subject. The compression process of the low density part in the chest image results in excessive compression, which increases the density of the mediastinum part and reduces the diagnostic performance. Conversely, in the case of a fat person, the compression is insufficient and There was a possibility that the space would be crushed white.
【０００６】本発明は上記問題点に鑑みなされたもので
あり、放射線画像のダイナミックレンジの圧縮方法にお
いて、被写体の違いや撮影部位の違いによって圧縮度合
いに過不足が生じることを回避できるようにすることを
目的とする。The present invention has been made in view of the above problems, and in a method of compressing a dynamic range of a radiographic image, it is possible to avoid an excess or deficiency in the degree of compression due to a difference in subject and a difference in imaged region. The purpose is to
【０００７】[0007]
【課題を解決するための手段】そのため請求項１記載の
発明では、オリジナル画像信号Ｓorgを処理して前記オ
リジナル画像よりもダイナミックレンジの狭い画像を担
持する処理済み画像信号Ｓprocを得る放射線画像のダイ
ナミックレンジ圧縮方法において、前記被写体の情報に
応じてダイナミックレンジ圧縮方法と圧縮度合いとの少
なくと一方を変化させるようにした。Therefore, according to the first aspect of the invention, the dynamic image of the radiation image is obtained by processing the original image signal Sorg to obtain a processed image signal Sproc carrying an image having a narrower dynamic range than the original image. In the range compression method, at least one of the dynamic range compression method and the compression degree is changed according to the information of the subject.
【０００８】かかる構成によると、被写体の違いによっ
て圧縮度合いが過不足とならないように、被写体の情報
に応じてダイナミックレンジの圧縮方法や圧縮度合いを
変化させる。According to this structure, the compression method and the compression degree of the dynamic range are changed according to the information of the subject so that the degree of compression is not excessive or insufficient depending on the difference of the subject.
【０００９】請求項２記載の発明では、前記被写体の情
報として、被写体のダイナミックレンジの情報と、被写
体の撮影部位の情報との少なくとも一方を含む構成とし
た。請求項３記載の発明では、各画素点に対応して非鮮
鋭マスク信号Ｓusを求め、該非鮮鋭マスク信号Ｓusの関
数である補正値ｆ１（Ｓus）によってオリジナル画像信
号Ｓorgを補正することで前記処理済み画像信号Ｓproc
を得る構成とした。According to the second aspect of the present invention, the information of the subject includes at least one of the information of the dynamic range of the subject and the information of the photographing region of the subject. In the invention according to claim 3, the non-sharp mask signal Sus is obtained corresponding to each pixel point, and the original image signal Sorg is corrected by the correction value f1 (Sus) which is a function of the non-sharp mask signal Sus. Image signal Sproc
Was obtained.
【００１０】請求項４記載の発明では、前記補正値ｆ１
（Ｓus）が、前記非鮮鋭マスク信号Ｓusの増大に伴って
傾きが変化する特性を有する構成とした。請求項５記載
の発明では、前記ダイナミックレンジの圧縮度合いを被
写体のダイナミックレンジの増大に応じて増大させる構
成とした。According to a fourth aspect of the invention, the correction value f1
(Sus) has a characteristic in which the inclination changes as the non-sharp mask signal Sus increases. According to a fifth aspect of the invention, the degree of compression of the dynamic range is increased in accordance with the increase of the dynamic range of the subject.
【００１１】請求項６記載の発明では、処理済み画像信
号Ｓprocのダイナミックレンジが基準のダイナミックレ
ンジに略一致するように前記ダイナミックレンジの圧縮
度合いを変化させる構成とした。According to the sixth aspect of the invention, the degree of compression of the dynamic range is changed so that the dynamic range of the processed image signal Sproc substantially matches the reference dynamic range.
【００１２】請求項７記載の発明では、前記ダイナミッ
クレンジ圧縮度合いの最大値を設け、該最大値を越える
圧縮度合いの設定を禁止するよう構成した。According to a seventh aspect of the present invention, the maximum value of the dynamic range compression degree is provided, and the setting of the compression degree exceeding the maximum value is prohibited.
【００１３】[0013]
【発明の実施の形態】以下に本発明の実施形態を説明す
る。一実施形態を示す図１は、本発明にかかる放射線画
像のダイナミックレンジ圧縮方法が適用された画像処理
装置を含む放射線画像読取装置を示すものであり、医療
用としての人体の撮影を行う例を示す。BEST MODE FOR CARRYING OUT THE INVENTION Embodiments of the present invention will be described below. FIG. 1 showing an embodiment shows a radiographic image reading apparatus including an image processing apparatus to which a method for compressing a dynamic range of a radiographic image according to the present invention is applied, and an example of imaging a human body for medical use is shown. Show.
【００１４】ここで、放射線発生源１は、放射線制御装
置２によって制御されて、被写体（人体胸部等）Ｍに向
けて放射線（一般的にはＸ線）を照射する。記録読取装
置３は、被写体Ｍを挟んで放射線源１と対向する面に変
換パネル４を備え、該変換パネル４は放射線源１からの
照射放射線量に対する被写体Ｍの放射線透過率分布に従
ったエネルギーを輝尽層に蓄積し、そこに被写体Ｍの潜
像を形成する。The radiation source 1 is controlled by the radiation control device 2 to irradiate the subject (human chest or the like) M with radiation (generally X-rays). The recording / reading device 3 is provided with a conversion panel 4 on a surface facing the radiation source 1 with the subject M interposed therebetween, and the conversion panel 4 has energy according to a radiation transmittance distribution of the subject M with respect to an irradiation radiation amount from the radiation source 1. Are accumulated in the photostimulation layer, and a latent image of the subject M is formed there.
【００１５】前記変換パネル４は、支持体上に輝尽層
を、輝尽性蛍光体の気相堆積、或いは輝尽性蛍光体塗料
塗布によって設けてあり、該輝尽層は環境による悪影響
及び損傷を遮断するために保護部材によって遮蔽若しく
は被覆される。The conversion panel 4 is provided with a stimulable layer on a support by vapor deposition of a stimulable phosphor or coating of a stimulable phosphor coating, and the stimulable layer has an adverse effect on the environment and It is shielded or covered by a protective member to block damage.
【００１６】該輝尽性蛍光体材料としては、例えば、特
開昭６１−７２０９１号公報、或いは、特開昭５９−７
５２００号公報に開示されるような材料が使われる。光
ビーム発生部（ガスレーザ，固体レーザ，半導体レーザ
等）５は、出射強度が制御された光ビームを発生し、そ
の光ビームは種々の光学系を経由して走査器６に到達
し、そこで偏向を受け、更に、反射鏡７で光路を偏向さ
せて、変換パネル４に輝尽励起走査光として導かれる。Examples of the photostimulable phosphor material include, for example, JP-A-61-272091 and JP-A-59-7.
The material as disclosed in Japanese Patent No. 5200 is used. A light beam generator (gas laser, solid-state laser, semiconductor laser, etc.) 5 generates a light beam whose emission intensity is controlled, and the light beam reaches a scanner 6 via various optical systems and is deflected there. Further, the light path is deflected by the reflecting mirror 7 and is guided to the conversion panel 4 as stimulated excitation scanning light.
【００１７】集光体８は、輝尽励起光が走査される変換
パネル４に近接して光ファイバである集光端が位置さ
れ、上記光ビームで走査された変換パネル４からの潜像
エネルギーに比例した発光強度の輝尽発光を受光する。The condensing body 8 has a condensing end, which is an optical fiber, located in the vicinity of the conversion panel 4 which is scanned by the stimulated excitation light, and the latent image energy from the conversion panel 4 which is scanned by the above light beam. The stimulated luminescence having a luminescence intensity proportional to is received.
【００１８】９は、集光体８から導入された光から輝尽
発光波長領域の光のみを通過させるフィルタであり、該
フィルタ９を通過した光は、フォトマル10に入射して、
その入射光に対応した電流信号に光電変換される。Reference numeral 9 is a filter that allows only light in the stimulated emission wavelength range from the light introduced from the light collector 8 to pass through. The light that has passed through the filter 9 is incident on the photomultiplier 10.
It is photoelectrically converted into a current signal corresponding to the incident light.
【００１９】フォトマル10からの出力電流は、電流／電
圧変換器11で電圧信号に変換され、増幅器12で増幅され
た後、Ａ／Ｄ変換器13で画素毎のディジタルデータから
なる放射線画像信号に変換される。An output current from the photomultiplier 10 is converted into a voltage signal by a current / voltage converter 11, amplified by an amplifier 12, and then a radiation image signal composed of digital data for each pixel by an A / D converter 13. Is converted to.
【００２０】そして、このディジタル放射線画像信号
（オリジナル画像信号Ｓorg）は、マイクロコンピュー
タを内蔵した画像処理装置14に順次出力される。15は画
像信号を記憶させておくための画像メモリ（磁気ディス
ク装置）である。Then, the digital radiation image signal (original image signal Sorg) is sequentially output to the image processing device 14 having a microcomputer built therein. Reference numeral 15 is an image memory (magnetic disk device) for storing image signals.
【００２１】また、16は画像処理装置14から直接又は前
記画像メモリ15から読み出された放射線画像信号をプリ
ンタ17に伝送するためのインターフェイスである。18は
読取ゲイン調整回路であり、この読取ゲイン調整回路18
により光ビーム発生部５の光ビーム強度調整、フォトマ
ル用高圧電源19の電源電圧調整によるフォトマル10のゲ
イン調整、電流／電圧変換器11と増幅器12のゲイン調
整、及びＡ／Ｄ変換器13の入力ダイナミックレンジの調
整が行われ、放射線画像信号の読取ゲインが総合的に調
整される。Reference numeral 16 is an interface for transmitting a radiation image signal read from the image processing device 14 directly or from the image memory 15 to the printer 17. Reference numeral 18 is a read gain adjustment circuit.
Adjusts the light beam intensity of the light beam generator 5, adjusts the gain of the photomultiplier 10 by adjusting the power supply voltage of the high voltage power supply 19 for photomultiplier, adjusts the gain of the current / voltage converter 11 and the amplifier 12, and A / D converter 13. The input dynamic range is adjusted, and the reading gain of the radiation image signal is adjusted comprehensively.
【００２２】尚、前記画像処理装置14に出力させるオリ
ジナル放射線画像信号Ｓorgの取得方法を、輝尽性蛍光
体を励起光で走査して発光せしめた輝尽発光を光電変換
して得る方法に限定するものではなく、例えば放射線フ
ィルムの画像を、光電変換によって読み取る方法や、被
写体を透過した放射線を蛍光体に照射されて蛍光に変換
し、該蛍光を光電変換して読み取る方法などであっても
良い。The method of acquiring the original radiation image signal Sorg output to the image processing device 14 is limited to the method of photoelectrically converting the stimulated emission obtained by scanning the stimulable phosphor with excitation light to emit light. However, for example, a method of reading an image on a radiation film by photoelectric conversion, a method of irradiating a phosphor with radiation passing through an object to convert it into fluorescence, and photoelectrically converting the fluorescence to read good.
【００２３】オリジナル放射線画像信号Ｓorgは、検出
された放射線の強度に比例する形でも良いし、検出され
た放射線の強度の対数に比例する形でも良いが、後者の
方が好ましい。The original radiation image signal Sorg may be proportional to the intensity of the detected radiation or may be proportional to the logarithm of the intensity of the detected radiation, but the latter is preferable.
【００２４】ここで、前記画像処理装置14には、入力さ
れるオリジナル画像信号Ｓorgのダイナミックレンジを
圧縮して、オリジナル画像よりもダイナミックレンジの
狭い画像を担持する処理済みの画像信号Ｓprocを得る画
像処理機能が備えられており、かかるダイナミックレン
ジ圧縮のための画像処理は、以下の式に従って行われ
る。
Ｓproc＝Ｓorg＋ｆ１（Ｓus）
上式で、Ｓusは、各画素点に対応して該各画素点を含む
所定マスク領域内のオリジナル画像信号Ｓorg を平均化
することにより求められる非鮮鋭マスク信号である。Here, the image processing device 14 compresses the dynamic range of the input original image signal Sorg to obtain a processed image signal Sproc carrying an image having a narrower dynamic range than the original image. A processing function is provided, and image processing for such dynamic range compression is performed according to the following equation. Sproc = Sorg + f1 (Sus) In the above formula, Sus is a non-sharp mask signal obtained by averaging the original image signal Sorg within a predetermined mask area including each pixel point corresponding to each pixel point.
【００２５】また、オリジナル画像信号Ｓorgに加算さ
れるｆ１（Ｓus）は、非鮮鋭マスク信号Ｓusの関数とし
て求められる補正値であり、非鮮鋭マスク信号Ｓusの増
大に応じて単調減少する。Further, f1 (Sus) added to the original image signal Sorg is a correction value obtained as a function of the unsharp mask signal Sus, and monotonically decreases as the unsharp mask signal Sus increases.
【００２６】尚、「単調減少」とは、非鮮鋭マスク信号
Ｓusの増大に伴ってｆ１（Ｓus）が必ず減少するもので
ある必要はなく、部分的にはＳusが変化してもｆ１（Ｓ
us）が変化しない領域が存在しても良いこととする。Note that "monotonically decreasing" does not necessarily mean that f1 (Sus) decreases with an increase in the non-sharp mask signal Sus, and f1 (Sus) may partially change even if Sus changes.
There may be areas where us) does not change.
【００２７】図２は、前記補正値ｆ１（Ｓus）の一例を
示す図であり、この図２において補正値ｆ１（Ｓus）
は、Ｓus１点よりも大きい領域では零で、前記Ｓus１点
から非鮮鋭マスク信号Ｓusの減少に応じてプラス側に一
定割合で増加する特性となっている。FIG. 2 is a diagram showing an example of the correction value f1 (Sus). In FIG. 2, the correction value f1 (Sus) is shown.
Is zero in a region larger than the Sus1 point and increases at a constant rate from the Sus1 point to the plus side in accordance with the decrease of the unsharp mask signal Sus.
【００２８】かかる特性を式で示すと、以下のようにな
る。
ｆ１（Ｓus）＝β（Ｓus１−Ｓus）……（Ｓus＜Ｓus
１）
ｆ１（Ｓus）＝０……（Ｓus≧Ｓus１）
本実施例では、画像信号Ｓorg の増大は濃度の増大を示
すから、前記補正値ｆ１（Ｓus）をオリジナル画像信号
Ｓorg に加算することで、低濃度側が持ち上げられて、
低濃度側の圧縮が行われる。The characteristic can be expressed by the following equation. f1 (Sus) = β (Sus1−Sus) ... (Sus <Sus
1) f1 (Sus) = 0 ... (Sus ≧ Sus1) In this embodiment, since the increase of the image signal Sorg indicates the increase of the density, the correction value f1 (Sus) is added to the original image signal Sorg. , The low concentration side is lifted,
Compression on the low density side is performed.
【００２９】かかる構成によると、非鮮鋭マスク信号Ｓ
usによって平均濃度として低い領域のダイナミックレン
ジを圧縮するから、微細構造のコントラストをそのまま
にして低濃度側が圧縮された画像を得ることができる。According to this structure, the non-sharp mask signal S
Since us compresses the dynamic range of a region having a low average density, it is possible to obtain an image in which the low density side is compressed while maintaining the contrast of the fine structure.
【００３０】更に、低濃度側における圧縮と、高濃度側
における圧縮とを同時に実行させることも可能である
が、この場合、低濃度側で要求される圧縮度合いと、高
濃度側で要求される圧縮度合いとが異なることが多いの
で、それぞれに異なる圧縮度合いを設定して低・高濃度
側でそれぞれに最適なダイナミックレンジの圧縮を行わ
せることが好ましい。Further, it is possible to simultaneously perform the compression on the low density side and the compression on the high density side. In this case, the degree of compression required on the low density side and the compression degree on the high density side are required. Since the degree of compression is often different, it is preferable to set different degree of compression for each so as to perform compression in the optimum dynamic range on the low and high density sides respectively.
【００３１】図３は、頭部正面画像に対して適合させた
補正値ｆ１（Ｓus）の特性例を示すものであり、かかる
図３に示す補正値ｆ１（Ｓus）の特性を式で示すと以下
のようになる。FIG. 3 shows an example of the characteristic of the correction value f1 (Sus) adapted to the head front image. The characteristic of the correction value f1 (Sus) shown in FIG. 3 is expressed by an equation. It looks like this:
【００３２】ｆ１（Ｓus）＝β1（Ｓus１−Ｓus）……
（Ｓus≦Ｓus１）
ｆ１（Ｓus）＝０……（Ｓus１＜Ｓus≦Ｓus２）
ｆ１（Ｓus）＝β2（Ｓus２−Ｓus）……（Ｓus＞Ｓus
２）
ここで、圧縮度合いを決定する係数βを、β1＝0.6、β
2＝0.4とした。F1 (Sus) = β 1 (Sus1−Sus) ...
(Sus ≦ Sus1) f1 (Sus) = 0 ... (Sus1 <Sus ≦ Sus2) f1 (Sus) = β 2 (Sus2-Sus) .... (Sus> Sus)
2) Here, the coefficient β that determines the compression degree is β 1 = 0.6, β
2 = 0.4.
【００３３】更に、圧縮濃度領域を決定するＳus１とＳ
us２とは、画像内の関心領域の最大信号値をＳmax、最
小信号値をＳminとしたときに、最大信号値Ｓmax，最小
信号値Ｓminの加重平均の重みｋを変えて以下のように
して演算した。Further, Sus1 and S for determining the compression density area
Us2 is calculated as follows by changing the weight k of the weighted average of the maximum signal value Smax and the minimum signal value Smin, where Smax is the maximum signal value of the region of interest in the image and Smin is the minimum signal value. did.
【００３４】Ｓus１＝ｋ・Ｓmax＋（１−ｋ）・Ｓmin…
…（ｋ＝0.5）
Ｓus２＝ｋ・Ｓmax＋（１−ｋ）・Ｓmin……（ｋ＝0.
9）
上記のように、高濃度側と低濃度側とでそれぞれに異な
る圧縮度合いに基づいて圧縮を行わせる構成であれば、
一定の傾き（係数β）を与えた場合のように、一方の濃
度域では最適な圧縮がなされるが、他方の濃度域は所期
の圧縮効果が得られなくなるといった問題がなく、高濃
度側と低濃度側とでそれぞれに最適な圧縮を同時に行わ
せることができる。Sus1 = k.Smax + (1-k) .Smin ...
… (K = 0.5) Sus2 = k · Smax + (1-k) · Smin …… (k = 0.
9) As described above, if the compression is performed on the high density side and the low density side based on different compression degrees,
As in the case where a constant slope (coefficient β) is given, optimum compression is performed in one concentration range, but there is no problem that the desired compression effect cannot be obtained in the other concentration range, Optimum compression can be performed simultaneously on the low density side and the low density side.
【００３５】特に医療用放射線画像においては、高濃度
側の圧縮度合いを低濃度側に比して小さくすることがよ
り好ましい。尚、上記実施例では、補正値ｆ１（Ｓus）
はいずれも非鮮鋭マスク信号Ｓusの増大に伴って単調減
少する１次関数として与えられる構成であったが、例え
ば図４に示すように傾きが変化する特性であっても良
く、この場合も、低濃度側と高濃度側とでそれぞれの要
求に合った圧縮度合い（傾き）を設定させるようにすれ
ば良い。Particularly in a medical radiation image, it is more preferable that the degree of compression on the high density side is smaller than that on the low density side. In the above embodiment, the correction value f1 (Sus)
In each of the configurations, the linear function that monotonously decreases with an increase in the unsharp mask signal Sus is given. However, for example, the characteristic may be such that the slope changes as shown in FIG. It suffices to set the degree of compression (slope) that meets the respective requirements on the low density side and the high density side.
【００３６】また、圧縮度合いは、前記係数βによって
変化させる構成に限定されるものではなく、他のパラメ
ータによって圧縮度合いを異ならせる構成としても良
い。ところで、オリジナル放射線画像Ｓorgのダイナミ
ックレンジは、被写体に応じて変化するため、一定の圧
縮度合いでダイナミックレンジ圧縮を行わせると、圧縮
度合いに過不足を生じ、最適な圧縮処理とならない場合
がある。Further, the degree of compression is not limited to the configuration that is changed by the coefficient β, and the configuration may be such that the degree of compression is changed by other parameters. By the way, the dynamic range of the original radiation image Sorg changes depending on the subject. Therefore, if the dynamic range compression is performed at a constant compression degree, the compression degree may become excessive or insufficient, and the optimum compression process may not be performed.
【００３７】そこで、本実施形態は、被写体の情報に応
じて、ダイナミックレンジ圧縮の圧縮度合いや圧縮方法
を変化させる構成となっている。図５は、圧縮度合いを
被写体のダイナミックレンジ（被写体情報）に応じて変
化させる実施形態の基本構成を示すブロック図ある。Therefore, the present embodiment is configured to change the compression degree and the compression method of the dynamic range compression according to the information of the subject. FIG. 5 is a block diagram showing a basic configuration of an embodiment in which the degree of compression is changed according to the dynamic range (subject information) of the subject.
【００３８】ここで、オリジナル画像信号Ｓorg取得手
段Ａは、前述のように、輝尽性蛍光体を励起光で走査し
て発光せしめた輝尽発光を光電変換して得る方法や、放
射線フィルムの画像を光電変換によって読み取る方法
や、被写体を透過した放射線を蛍光体に照射されて蛍光
に変換し、該蛍光を光電変換して読み取る方法などであ
る。Here, as described above, the original image signal Sorg acquisition means A is obtained by photoelectrically converting the stimulated luminescence obtained by scanning the stimulable phosphor with the excitation light and causing the stimulable phosphor to emit light. There are a method of reading an image by photoelectric conversion, a method of irradiating a phosphor with radiation passing through an object to convert it into fluorescence, and photoelectrically converting the fluorescence to read it.
【００３９】また、ダイナミックレンジ計測手段Ｂは、
被写体のダイナミックレンジを計測する手段であり、例
えば、オリジナル画像信号Ｓorgのヒストグラム解析か
ら、被写体画像領域における最大信号値，最小信号値を
求め、これらから被写体のダイナミックレンジを決定す
る。Further, the dynamic range measuring means B is
It is a means for measuring the dynamic range of a subject, for example, the maximum signal value and the minimum signal value in the subject image area are obtained from the histogram analysis of the original image signal Sorg, and the dynamic range of the subject is determined from these.
【００４０】また、オリジナル画像のプロファイル情報
から最大信号値，最小信号値を求める構成であっても良
い。更に、画像信号を用いずに被写体のダイナミックレ
ンジを検知するものであっても良く、例えば、被写体を
透過した放射線の強度分布を直接測定したり、また、前
記放射線を一旦光に変換し、該光の強度分布を測定した
り、更に、被写体の厚さを測定するものであっても良
い。The maximum signal value and the minimum signal value may be obtained from the profile information of the original image. Further, the dynamic range of the subject may be detected without using an image signal, for example, the intensity distribution of the radiation transmitted through the subject is directly measured, or the radiation is once converted into light, The intensity distribution of light may be measured, or the thickness of the subject may be measured.
【００４１】ダイナミックレンジ設定手段Ｃは、処理済
み画像信号Ｓprocのダイナミックレンジをどのようにす
るかを設定するための手段であり、例えば、オペレータ
がキーボード操作によって入力・設定する構成であった
り、また、装置内に既定値として記憶されているもので
あったり、更には、複数種の既定値の中からオペレータ
が選択したり、既定値をオペレータが変更して設定する
構成とする。The dynamic range setting means C is means for setting how to set the dynamic range of the processed image signal Sproc. For example, the operator inputs or sets by keyboard operation, or The configuration is such that it is stored as a default value in the device, further, the operator selects from a plurality of types of default values, or the operator changes and sets the default value.
【００４２】圧縮率決定手段Ｄは、後述する圧縮手段Ｅ
によってダイナミックレンジの圧縮を行うときの圧縮率
を可変設定する手段である。ここで、圧縮率は、圧縮率
＝（処理済み画像信号Ｓprocのダイナミックレンジ）／
（オリジナル画像信号Ｓorgのダイナミックレンジ）と
して定義し、圧縮率の低下は圧縮度合いの増大を示すも
のとする。The compression rate determining means D is a compression means E described later.
Is a means for variably setting the compression ratio when performing compression of the dynamic range. Here, the compression rate is: compression rate = (dynamic range of processed image signal Sproc) /
(Dynamic range of the original image signal Sorg), and a decrease in compression rate indicates an increase in compression degree.
【００４３】前記圧縮率決定手段Ｄでは、被写体のダイ
ナミックレンジを圧縮率に変換する図６〜図８に示すよ
うな変換テーブルを備え、該変換テーブルを参照して圧
縮率を決定する構成とすることができる。The compression rate determination means D is provided with a conversion table as shown in FIGS. 6 to 8 for converting the dynamic range of the subject into a compression rate, and the compression rate is determined by referring to the conversion table. be able to.
【００４４】また、前記ダイナミックレンジ設定手段Ｃ
で設定される処理済み画像信号Ｓprocのダイナミックレ
ンジ（基準ダイナミックレンジ）と、前記ダイナミック
レンジ計測手段Ｂで計測される被写体のダイナミックレ
ンジとの比を演算して、既定されている処理済み画像信
号Ｓprocのダイナミックレンジに合わせるような圧縮率
の設定を行わせる構成であっても良い。The dynamic range setting means C
The calculated dynamic range (reference dynamic range) of the processed image signal Sproc and the dynamic range of the object measured by the dynamic range measuring means B are calculated to calculate a predetermined processed image signal Sproc. The compression rate may be set so as to match the dynamic range of the above.
【００４５】上記のように、被写体のダイナミックレン
ジに応じて圧縮率を変化させる構成とすれば、被写体
（人体）の体格の差によって圧縮率に過不足が生じた
り、処理済み画像信号Ｓprocのダイナミックレンジが大
きくばらつくことを回避でき、安定した処理済み画像を
提供できる。As described above, if the compression rate is changed according to the dynamic range of the subject, the compression rate may become excessive or insufficient due to the difference in the physique of the subject (human body), or the dynamics of the processed image signal Sproc may be generated. It is possible to avoid a large range variation and provide a stable processed image.
【００４６】圧縮手段Ｅは、非鮮鋭マスク信号Ｓusの関
数である補正値ｆ１（Ｓus）でオリジナル画像信号Ｓor
gを補正することで、ダイナミックレンジが圧縮された
処理済み画像信号Ｓprocを得る構成とすることができ、
Ｓproc＝Ｓorg−ｆ１（Ｓus）とした場合には、前記補
正値ｆ１（Ｓus）は非鮮鋭マスク信号Ｓusの増加に応じ
て単調増加する関数であり、また、Ｓproc＝Ｓorg＋ｆ
１（Ｓus）とした場合には、前記補正値ｆ１（Ｓus）は
非鮮鋭マスク信号Ｓusの増加に応じて単調減少する関数
である。The compression means E uses the correction value f1 (Sus) which is a function of the unsharp mask signal Sus to generate the original image signal Sor.
By correcting g, it is possible to obtain a processed image signal Sproc whose dynamic range is compressed,
When Sproc = Sorg-f1 (Sus), the correction value f1 (Sus) is a function that monotonically increases as the non-sharp mask signal Sus increases, and Sproc = Sorg + f.
When it is set to 1 (Sus), the correction value f1 (Sus) is a function that monotonically decreases according to the increase of the non-sharp mask signal Sus.
【００４７】前記補正値ｆ１（Ｓus）は、具体的には、
例えばｆ１（Ｓus）＝β（Ｓus−Ｓus１）として与えら
れる。ここで、前記係数β（０≦β≦1.0）を変化させ
ることで圧縮率が変化することになり、上記式では圧縮
率＝（１−β）／１となるから、設定された圧縮率から
前記係数βを求めることで、設定された圧縮率に応じた
処理を行わせることができる。The correction value f1 (Sus) is, concretely,
For example, it is given as f1 (Sus) = β (Sus-Sus1). Here, the compression rate is changed by changing the coefficient β (0 ≦ β ≦ 1.0), and in the above equation, the compression rate = (1−β) / 1. Therefore, from the set compression rate By obtaining the coefficient β, it is possible to perform processing according to the set compression rate.
【００４８】尚、前記圧縮率決定手段Ｄで圧縮率を決定
させる代わりに、前記係数βを直接決定させる構成とし
ても良い。また、前記圧縮手段Ｅは、非鮮鋭マスク信号
Ｓusを用いる構成の他、畳み込み演算によって圧縮を行
うものであっても良い。Instead of the compression rate determining means D determining the compression rate, the coefficient β may be directly determined. Further, the compression means E may use a non-sharp mask signal Sus or may perform compression by a convolution operation.
【００４９】図９は、被写体のダイナミックレンジに応
じて圧縮率を変更すると共に、圧縮条件を撮影部位（被
写体情報）の解析結果に基づいて決定する実施形態の基
本構成を示すブロック図である。FIG. 9 is a block diagram showing the basic construction of an embodiment in which the compression ratio is changed according to the dynamic range of the object and the compression condition is determined based on the analysis result of the imaged region (object information).
【００５０】尚、前記図５のブロック図と同一要素には
同一符合を付して詳細な説明を省略する。この図９にお
いて、撮影部位解析手段Ｆは、オリジナル画像信号Ｓor
gに基づいて被写体の撮影部位を解析する手段であり、
例えばヒストグラム解析，輪郭抽出，プロファイル解析
などの公知の撮影部位解析方法を用いることができる。The same elements as those in the block diagram of FIG. 5 are designated by the same reference numerals and detailed description thereof will be omitted. In FIG. 9, the imaged region analysis means F displays the original image signal Sor.
It is a means to analyze the imaged part of the subject based on g,
For example, a known imaging region analysis method such as histogram analysis, contour extraction, profile analysis or the like can be used.
【００５１】また、ニューラルネットワークを用いた認
識などを適用しても良く、更には、これらの解析方法の
組み合わせによって撮影部位を解析させる構成であって
も良い。Further, recognition using a neural network may be applied, and further, a configuration in which an imaged region is analyzed by a combination of these analysis methods may be used.
【００５２】圧縮条件決定手段Ｇは、前記撮影部位解析
手段Ｆで解析された撮影部位に応じてダイナミックレン
ジを圧縮する画像処理における各種の条件を決定する手
段であり、具体的には、以下のような条件を決定する。The compression condition determining means G is means for determining various conditions in the image processing for compressing the dynamic range according to the imaged region analyzed by the imaged region analyzing means F. Determine such conditions.
【００５３】ａ．圧縮手段Ｅで用いる圧縮方法……非鮮
鋭マスク方法又は畳み込み演算法（コンボリュージョン
法）のいずれを選択するかを、撮影部位に応じて選択す
る。ｂ．補正関数形……例えば非鮮鋭マスク法を用いる
場合で、信号補正関数ｆ１（Ｓus）がｆ１（Ｓus）＝β
（Ｓus１−Ｓus）として与えられる場合には、前記係数
β又は基準値Ｓus１の値を、撮影部位に応じて変化させ
る。尚、係数βやＳus１を変えることでダイナミックレ
ンジ圧縮の度合いやダイナミックレンジ圧縮の適用範囲
が変化する。A. The compression method used by the compression unit E ... Which of the non-sharp mask method and the convolution calculation method (convolution method) is selected is selected according to the imaging region. b. Correction function form ... For example, when the non-sharp mask method is used, the signal correction function f1 (Sus) is f1 (Sus) = β
When given as (Sus1-Sus), the value of the coefficient β or the reference value Sus1 is changed according to the imaging region. The degree of dynamic range compression and the applicable range of dynamic range compression are changed by changing the coefficient β and Sus1.
【００５４】ｃ．処理後の画像信号のダイナミックレン
ジ……撮影部位毎に、処理後のダイナミックレンジの既
定値を複数記憶しておき、解析された撮影部位に応じて
前記記憶値の中から選択する。又は、撮影部位毎に、被
写体のダイナミックレンジを圧縮率に変換するテーブル
を複数備え、解析された撮影部位に応じて前記複数の変
換テーブルの中から選択する。C. Dynamic range of image signal after processing ... A plurality of default values of the dynamic range after processing are stored for each imaging region, and selected from the stored values according to the analyzed imaging region. Alternatively, a plurality of tables for converting the dynamic range of the subject into the compression rate are provided for each imaging region, and the table is selected from the plurality of conversion tables according to the analyzed imaging region.
【００５５】上記実施形態では、撮影部位をオリジナル
画像信号Ｓorgの信号解析によって求めるようにした
が、オペレータの操作等によって撮影部位の情報を入力
する構成としても良く、かかる実施形態に対応する構成
を図１０のブロック図に示してある。In the above embodiment, the imaged region is obtained by signal analysis of the original image signal Sorg. However, the imaged region information may be input by an operator's operation or the like. This is shown in the block diagram of FIG.
【００５６】この図１０において、撮影部位入力手段Ｈ
は、圧縮処理の対象となっている放射線画像における撮
影部位がどこであるかを示すデータを入力するための手
段であり、オペレータがキーボードを操作して撮影部位
情報を直接入力したり、予め設定されている処理メニュ
ーをオペレータが選択する構成とする。In FIG. 10, the imaging region input means H
Is a means for inputting data indicating where an imaged part in the radiographic image to be compressed is located, and the operator operates the keyboard to directly input the imaged part information or is set in advance. The operator selects the existing processing menu.
【００５７】尚、前記撮影部位の情報と共に、撮影部位
に応じて画像処理条件等を同時に入力させる構成とする
こともできる。圧縮条件決定手段Ｇでは、前記実施形態
と同様に、撮影部位の情報に従って圧縮処理方法，関数
形，処理後のダイナミックレンジなどを決定する。It should be noted that it is also possible to have a structure in which the image processing conditions and the like are simultaneously input in accordance with the imaged region together with the information of the imaged region. The compression condition determining means G determines a compression processing method, a functional form, a processed dynamic range, and the like according to the information of the imaged region, as in the above embodiment.
【００５８】上記のように、被写体の撮影部位の情報に
基づいて圧縮条件を可変設定させれば、撮影部位に応じ
た要求圧縮率の変化に対応することができ、撮影部位の
変化によって圧縮率が過不足となることを回避でき、安
定して最適な圧縮処理を行わせることができる。As described above, if the compression condition is variably set on the basis of the information of the photographic part of the subject, it is possible to cope with the change in the required compression rate depending on the photographic part, and the compression rate changes depending on the change of the photographic part. It is possible to avoid the excess and deficiency of the value of "," and to perform stable and optimal compression processing.
【００５９】尚、上記のように、被写体のダイナミック
レンジ情報などに基づいて圧縮率（圧縮度合い）を変化
させる構成の場合には、圧縮率の最大値を予め設定して
おき、被写体のダイナミックレンジ情報等に基づいて前
記最大値を越える圧縮率が設定された場合には、圧縮率
を前記最大値に制限し、最大値を越える圧縮率の設定を
禁止する構成とすることが好ましい。When the compression rate (compression degree) is changed based on the dynamic range information of the subject as described above, the maximum value of the compression rate is set in advance and the dynamic range of the subject is set. When a compression rate exceeding the maximum value is set based on information or the like, it is preferable to limit the compression rate to the maximum value and prohibit the setting of the compression rate exceeding the maximum value.
【００６０】前記実施形態において、非鮮鋭マスクのマ
スクサイズ或いは周波数特性は画像の診断性を左右する
重要なパラメータである。ダイナミックレンジ圧縮処理
では、被写体の大まかな構造の変化（肺野部，縦隔部な
どの滑らかな信号差）に対応する超低周波数成分のみを
非鮮鋭マスク信号Ｓusとして抽出し、Ｓusに基づいて補
正値ｆ１（Ｓus）を設定することによって、細かい構造
物の変化（骨，血管など）を維持しつつ、全体の濃度範
囲を圧縮することができるのである。In the above embodiment, the mask size or frequency characteristic of the non-sharp mask is an important parameter that influences the diagnostic ability of the image. In the dynamic range compression processing, only the ultra-low frequency component corresponding to the rough structural change of the subject (smooth signal difference such as lung field and mediastinum) is extracted as the unsharp mask signal Sus, and based on Sus By setting the correction value f1 (Sus), it is possible to compress the entire density range while maintaining the small changes in the structure (bones, blood vessels, etc.).
【００６１】マスクサイズが小さいと、非鮮鋭マスク信
号Ｓusが被写体の大まかな変化に対応する超低周波数成
分だけでなく、細かい構造物の変化に相当する周波数成
分をも含んでしまい、非鮮鋭マスク信号Ｓusに基づく補
正値を加算することによって細かい構造物の変化が打ち
消されて骨や血管などのコントラストが低下してしま
う。When the mask size is small, the non-sharp mask signal Sus includes not only the ultra-low frequency component corresponding to the rough change of the object but also the frequency component corresponding to the change of the fine structure, so that the non-sharp mask signal is generated. By adding the correction value based on the signal Sus, the change of the fine structure is canceled and the contrast of the bone, the blood vessel, etc. is lowered.
【００６２】一方、マスクサイズが大きいと前述のよう
に信号値の変化が急激な部分での非鮮鋭画像のエッジ切
れが悪くなり、圧縮を行いたい領域と行いたくない領域
との境界付近で望まない圧縮が施されてしまう。On the other hand, when the mask size is large, the edge cut of the non-sharp image in the portion where the signal value changes abruptly becomes worse as described above, and it is desired near the boundary between the region to be compressed and the region not to be compressed. There is no compression.
【００６３】また、マスクサイズを更に大きくし過ぎる
と、被写体の大まかな変化に相当する周波数成分までも
失ってしまう（極端な場合全く平坦な画像になってしま
う）ので、非鮮鋭マスク信号Ｓusに基づく補正値を加算
してもダイナミックレンジ圧縮効果が得られなくなる。Further, if the mask size is made too large, the frequency component corresponding to the rough change of the object is also lost (in the extreme case, the image becomes a flat image). The dynamic range compression effect cannot be obtained even if the correction value based on the above is added.
【００６４】以上のような観点で発明者が検討した結
果、マスクサイズの大きさは等身大画像上の長さで10mm
から60mmが好ましく、より好ましくは15mmから30mmであ
り、最も好ましくは20mmから30mmであることを見出し
た。As a result of the inventor's examination from the above viewpoints, the size of the mask size is 10 mm in the length on the life-size image.
To 60 mm are preferred, more preferably from 15 mm to 30 mm, most preferably from 20 mm to 30 mm.
【００６５】マスクサイズが10mmよりも小さいと、細か
い構造物の変化に対応する周波数成分が急激に増加する
ので、このようなマスクサイズで求めた非鮮鋭マスク信
号Ｓusに基づいて補正値を設定すると、著しく診断性能
が低下してしまう。When the mask size is smaller than 10 mm, the frequency component corresponding to the change of the fine structure sharply increases. Therefore, when the correction value is set based on the non-sharp mask signal Sus obtained by such a mask size. , The diagnostic performance is significantly reduced.
【００６６】また、特に、胸部画像や腹部画像において
は、マスクサイズを15mm以上にすれば、Ｓusは大動脈な
どの太い血管に対応する周波数成分を持たなくなり、マ
スクサイズを20mm以上にすればＳusは肋骨などに対応す
る比較的低周波数であるが、コントラストを下げたくな
い周波数成分を含まなくなるので、診断性能の高い画像
が得られる。In particular, in the chest image and the abdominal image, if the mask size is 15 mm or more, Sus does not have frequency components corresponding to large blood vessels such as aorta, and if the mask size is 20 mm or more, Sus is Although it has a relatively low frequency corresponding to ribs and the like, it does not include a frequency component whose contrast is not desired to be lowered, so that an image with high diagnostic performance can be obtained.
【００６７】ここでマスクサイズは、矩形ならば短辺の
長さと長辺の長さの平均値、正方形なら一辺の長さ、円
ならば直径、楕円ならば長径と短径の平均値を指す。ま
た、マスクサイズのかわりに非鮮鋭マスクのもつ周波数
特性で記述すると、非鮮鋭マスクの変調伝達関数が0.01
サイクル／mmのとき0.5以上かつ0.06サイクル／mmのと
き0.5以下であることが好ましく、より好ましくは0.02
サイクル／mmのとき0.5以上かつ0.04サイクル／mmのと
き0.5以下、更に好ましくは0.02サイクル／mmのとき0.5
以上かつ0.03サイクル／mmのとき0.5以下である。Here, the mask size refers to the average value of the length of the short side and the length of the long side for a rectangle, the length of one side for a square, the diameter for a circle, and the average value of the major axis and the minor axis for an ellipse. . If the frequency characteristics of the unsharp mask are used instead of the mask size, the modulation transfer function of the unsharp mask is 0.01
Cycle / mm is preferably 0.5 or more and 0.06 cycle / mm is preferably 0.5 or less, more preferably 0.02
Cycle / mm 0.5 or more and 0.04 cycles / mm 0.5 or less, more preferably 0.02 cycles / mm 0.5
Above 0.5 and below at 0.03 cycles / mm.
【００６８】尚、本発明のように非鮮鋭マスク信号Ｓus
を求めるための平均化処理において、マスク領域内にお
ける中心画素と周辺画素との信号差の絶対値に応じた重
み付けを行う、及び／又は、マスク領域内における中心
画素と周辺画素との位置関係に応じた重み付けを行うこ
とにより、信号値の変化が急激な部分での非鮮鋭画像の
エッジ切れの劣化を防ぐことができ、好ましいマスクサ
イズの範囲が10mm〜60mmであったのを10mm〜80mmに広げ
ることができる。The non-sharp mask signal Sus as in the present invention is used.
In the averaging process for obtaining the value, weighting is performed according to the absolute value of the signal difference between the central pixel and the peripheral pixels in the mask area, and / or the positional relationship between the central pixel and the peripheral pixels in the mask area is determined. By performing the weighting according to, it is possible to prevent the deterioration of the edge cut of the non-sharp image in the portion where the change of the signal value is abrupt, and the preferable mask size range was 10 mm to 60 mm from 10 mm to 80 mm. Can be expanded.
【００６９】更に本発明において、非鮮鋭マスク信号Ｓ
usの関数である補正値ｆ１（Ｓus）の絶対値の最大値は
被写体の関心領域のダイナミックレンジの１／８から１
／２であることが好ましい。Further, in the present invention, the non-sharp mask signal S
The maximum absolute value of the correction value f1 (Sus), which is a function of us, is 1/8 to 1 of the dynamic range of the region of interest of the subject.
It is preferably / 2.
【００７０】例えば、被写体の関心領域のダイナミック
レンジが２桁であるとき、圧縮補正量の絶対値の最大値
は１／４桁から１桁であることが好ましい。また、補正
値ｆ１（Ｓus）が、β（Ｓus１−Ｓus）の様に非鮮鋭マ
スク信号Ｓusの一次関数で表されるとき、その傾きであ
って圧縮度合いを決めるβの好ましい範囲は0.2〜1.0で
あり、より好ましくは0.4〜0.8である。For example, when the dynamic range of the region of interest of the subject is 2 digits, the maximum absolute value of the compression correction amount is preferably 1/4 digit to 1 digit. Further, when the correction value f1 (Sus) is represented by a linear function of the unsharp mask signal Sus like β (Sus1-Sus), the preferable range of β which is the slope and determines the compression degree is 0.2 to 1.0. And more preferably 0.4 to 0.8.
【００７１】補正量が小さすぎるとダイナミックレンジ
圧縮効果が現れず、一方、補正量が大きすぎるとオリジ
ナル画像における領域毎の濃度の大小関係が逆転して
（例えば肺野の平均濃度よりも縦隔部の平均濃度の方が
高くなって）しまい、診断に耐えない画像になってしま
う。If the correction amount is too small, the dynamic range compression effect does not appear. On the other hand, if the correction amount is too large, the size relationship of the densities of the regions in the original image is reversed (for example, the median density is higher than the average density of the lung field). The average density of the area becomes higher), resulting in an image that cannot be diagnosed.
【００７２】例えば、上記一次関数の傾きβを１より大
きくした場合このような問題が起こる。For example, such a problem occurs when the slope β of the linear function is set to be larger than 1.
【００７３】[0073]
【発明の効果】以上説明したように本発明にかかる放射
線画像のダイナミックレンジ圧縮方法によると、被写体
のダイナミックレンジや撮影部位などの情報に基づい
て、圧縮特性を変化させるようにしたので、被写体の違
いや撮影部位の違いによって圧縮度合いの過不足を生じ
たりすることを回避できるという効果がある。As described above, according to the method of compressing a dynamic range of a radiation image according to the present invention, the compression characteristic is changed based on the information such as the dynamic range of the subject and the region to be imaged. There is an effect that it is possible to avoid the occurrence of excess or deficiency of the compression degree due to the difference or the difference in the imaging region.
【図１】本発明が適用される放射線画像読取処理装置を
示すブロック図。FIG. 1 is a block diagram showing a radiation image reading processing apparatus to which the present invention is applied.
【図２】非鮮鋭マスク信号に応じた補正値の関数形の一
例を示す線図。FIG. 2 is a diagram showing an example of a functional form of a correction value according to an unsharp mask signal.
【図３】高・低濃度側で圧縮を行う補正値の特性例を示
す線図。FIG. 3 is a diagram showing a characteristic example of a correction value for performing compression on the high and low density sides.
【図４】高・低濃度側で圧縮を行う補正値の特性例を示
す線図。FIG. 4 is a diagram showing a characteristic example of a correction value for performing compression on the high and low density sides.
【図５】被写体情報に応じて圧縮特性を変化させる実施
形態のブロック図。FIG. 5 is a block diagram of an embodiment in which a compression characteristic is changed according to subject information.
【図６】被写体のダイナミックレンジに応じた圧縮率特
性例を示す線図。FIG. 6 is a diagram showing a compression rate characteristic example according to a dynamic range of a subject.
【図７】被写体のダイナミックレンジに応じた圧縮率特
性例を示す線図。FIG. 7 is a diagram showing a compression rate characteristic example according to a dynamic range of a subject.
【図８】被写体のダイナミックレンジに応じた圧縮率特
性例を示す線図。FIG. 8 is a diagram showing a compression rate characteristic example according to a dynamic range of a subject.
【図９】被写体情報に応じて圧縮特性を変化させる実施
形態のブロック図。FIG. 9 is a block diagram of an embodiment in which a compression characteristic is changed according to subject information.
【図１０】被写体情報に応じて圧縮特性を変化させる実
施形態のブロック図。FIG. 10 is a block diagram of an embodiment in which a compression characteristic is changed according to subject information.
１…放射線発生源 ３…記録読取装置 14…画像処理装置 15…画像メモリ 16…インタフェイス 17…プリンタ 1 ... Radiation source 3 ... Recording / reading device 14 ... Image processing device 15 ... Image memory 16 ... Interface 17 ... Printer
───────────────────────────────────────────────────── フロントページの続き Ｆターム(参考） 2H013 CZ01 4C093 AA28 CA01 EB05 FF29 5B057 AA08 BA03 CA02 CA08 CA12 CA16 CB02 CB08 CB12 CB16 CC01 CE11 CG10 CH08 CH11 5C077 LL01 LL19 NN02 PP10 PP15 PP46 PP68 PQ08 PQ18 PQ22 RR21 SS03 ─────────────────────────────────────────────────── ─── Continued front page F-term (reference) 2H013 CZ01 4C093 AA28 CA01 EB05 FF29 5B057 AA08 BA03 CA02 CA08 CA12 CA16 CB02 CB08 CB12 CB16 CC01 CE11 CG10 CH08 CH11 5C077 LL01 LL19 NN02 PP10 PP15 PP46 PP68 PQ08 PQ18 PQ22 RR21 SS03
Claims (7)
オリジナル画像を表すオリジナル画像信号Ｓorgを処理
して前記オリジナル画像よりもダイナミックレンジの狭
い画像を担持する処理済み画像信号Ｓprocを得る放射線
画像のダイナミックレンジ圧縮方法において、 前記被写体の情報に応じてダイナミックレンジ圧縮方法
と圧縮度合いとの少なくと一方を変化させることを特徴
とする放射線画像のダイナミックレンジ圧縮方法。1. A dynamic image of a radiation image for processing an original image signal Sorg representing an original image based on radiation image information transmitted through an object to obtain a processed image signal Sproc carrying an image having a narrower dynamic range than the original image. In the range compression method, a dynamic range compression method for a radiation image, wherein at least one of a dynamic range compression method and a compression degree is changed according to the information of the subject.
ミックレンジの情報と、被写体の撮影部位の情報との少
なくとも一方を含むことを特徴とする請求項１記載の放
射線画像のダイナミックレンジ圧縮方法。2. The method of compressing a dynamic range of a radiation image according to claim 1, wherein the information on the subject includes at least one of information on a subject's dynamic range and information on an imaged region of the subject.
を求め、該非鮮鋭マスク信号Ｓusの関数である補正値ｆ
１（Ｓus）によってオリジナル画像信号Ｓorgを補正す
ることで前記処理済み画像信号Ｓprocを得る構成である
ことを特徴とする請求項１又は２記載の放射線画像のダ
イナミックレンジ圧縮方法。3. An unsharp mask signal Sus corresponding to each pixel point.
Correction value f which is a function of the unsharp mask signal Sus
3. The method for compressing a dynamic range of a radiation image according to claim 1, wherein the processed image signal Sproc is obtained by correcting the original image signal Sorg with 1 (Sus).
スク信号Ｓusの増大に伴って傾きが変化する特性を有す
ることを特徴とする請求項３記載の放射線画像のダイナ
ミックレンジ圧縮方法。4. The method of compressing a dynamic range of a radiation image according to claim 3, wherein the correction value f1 (Sus) has a characteristic that the inclination changes with an increase in the non-sharp mask signal Sus.
被写体のダイナミックレンジの増大に応じて増大させる
ことを特徴とする請求項１〜４のいずれか１つに記載の
放射線画像のダイナミックレンジ圧縮方法。5. The degree of compression of the dynamic range
The method for compressing a dynamic range of a radiation image according to claim 1, wherein the method is increased in accordance with an increase in a dynamic range of a subject.
クレンジが基準のダイナミックレンジに略一致するよう
に前記ダイナミックレンジの圧縮度合いを変化させるこ
とを特徴とする請求項１〜４のいずれか１つに記載の放
射線画像のダイナミックレンジ圧縮方法。6. The compression degree of the dynamic range is changed so that the dynamic range of the processed image signal Sproc substantially matches a reference dynamic range. A method for compressing a dynamic range of a radiographic image described.
値を設け、該最大値を越える圧縮度合いの設定を禁止す
ることを特徴とする請求項１〜６のいずれか１つに記載
の放射線画像のダイナミックレンジ圧縮方法。7. The dynamic of a radiation image according to claim 1, wherein a maximum value of the compression degree of the dynamic range is provided, and setting of the compression degree exceeding the maximum value is prohibited. Range compression method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003074985A JP3637049B2 (en) | 1993-03-30 | 2003-03-19 | Dynamic range compression method for radiographic images |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7203193 | 1993-03-30 | ||
| JP5-72031 | 1993-03-30 | ||
| JP2003074985A JP3637049B2 (en) | 1993-03-30 | 2003-03-19 | Dynamic range compression method for radiographic images |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP07649593A Division JP3425780B2 (en) | 1993-03-30 | 1993-04-02 | Dynamic range compression method for radiographic images |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003310589A true JP2003310589A (en) | 2003-11-05 |
| JP3637049B2 JP3637049B2 (en) | 2005-04-06 |
Family
ID=29551512
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003074985A Expired - Lifetime JP3637049B2 (en) | 1993-03-30 | 2003-03-19 | Dynamic range compression method for radiographic images |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3637049B2 (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011239804A (en) * | 2010-05-14 | 2011-12-01 | Shimadzu Corp | X-ray examination apparatus |
| JP2018033745A (en) * | 2016-08-31 | 2018-03-08 | 富士フイルム株式会社 | Image processing apparatus, method, and program |
| US9959614B2 (en) | 2013-09-30 | 2018-05-01 | Fujifilm Corporation | Image processing device and method |
-
2003
- 2003-03-19 JP JP2003074985A patent/JP3637049B2/en not_active Expired - Lifetime
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011239804A (en) * | 2010-05-14 | 2011-12-01 | Shimadzu Corp | X-ray examination apparatus |
| US9959614B2 (en) | 2013-09-30 | 2018-05-01 | Fujifilm Corporation | Image processing device and method |
| JP2018033745A (en) * | 2016-08-31 | 2018-03-08 | 富士フイルム株式会社 | Image processing apparatus, method, and program |
| US10430930B2 (en) | 2016-08-31 | 2019-10-01 | Fujifilm Corporation | Image processing apparatus, image processing method, and image processing program for performing dynamic range compression process |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3637049B2 (en) | 2005-04-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3467285B2 (en) | Radiation image processing method | |
| US5471987A (en) | Method of compressing a dynamic range for a radiation image | |
| EP0440166B1 (en) | Method for compressing dynamic ranges of images | |
| JP3785817B2 (en) | Radiographic image processing device | |
| US7483556B2 (en) | Energy subtraction processing method and apparatus | |
| JP4129598B2 (en) | Image processing apparatus and image processing method | |
| JP3814864B2 (en) | Radiographic image processing condition determining apparatus, image processing apparatus, radiographic image processing condition determining method, and image processing method | |
| JPH05174141A (en) | Method and device for image processing | |
| JPH06292008A (en) | Dynamic range compression processing unit for radiation picture | |
| JP2509503B2 (en) | Image processing method and apparatus | |
| US6118892A (en) | Method for automatic detection of region of interest for digital x-ray detectors using a filtered histogram | |
| JP2000079110A (en) | Picture processor | |
| JPH08294006A (en) | Method and device for applying dynamic range compression processing to image | |
| JPH06339025A (en) | Dynamic range compressing method for radiation picture | |
| JP3455566B2 (en) | Dynamic range compression method for radiographic images | |
| JP3637049B2 (en) | Dynamic range compression method for radiographic images | |
| JP3434532B2 (en) | Dynamic range compression method for radiographic images | |
| JP2000023950A (en) | Image processing device for radiation image | |
| JP2002008009A (en) | Image processing condition determination method and device | |
| JP2005021456A (en) | Device, method and program of image processing for radiographic image | |
| JP3765320B2 (en) | Image processing condition determination apparatus and image processing apparatus for radiation image | |
| JP2001086409A (en) | Device and method for image processing | |
| JP2852794B2 (en) | Digital radiation image signal processing device | |
| JP3586274B2 (en) | Apparatus and method for dynamic range compression processing of radiographic image | |
| JP4258092B2 (en) | Image processing apparatus and image processing method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040810 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20041008 |
|
| 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: 20050104 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050106 |
|
| 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: 20090114 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100114 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110114 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120114 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130114 Year of fee payment: 8 |