JP4802333B2 - 画像合成における露光比決定方法 - Google Patents
画像合成における露光比決定方法 Download PDFInfo
- Publication number
- JP4802333B2 JP4802333B2 JP2007041053A JP2007041053A JP4802333B2 JP 4802333 B2 JP4802333 B2 JP 4802333B2 JP 2007041053 A JP2007041053 A JP 2007041053A JP 2007041053 A JP2007041053 A JP 2007041053A JP 4802333 B2 JP4802333 B2 JP 4802333B2
- Authority
- JP
- Japan
- Prior art keywords
- exposure ratio
- image
- signal
- exposure
- random noise
- 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 34
- 239000000203 mixture Substances 0.000 title claims description 11
- 238000011156 evaluation Methods 0.000 claims description 32
- 230000015572 biosynthetic process Effects 0.000 claims description 23
- 238000003786 synthesis reaction Methods 0.000 claims description 23
- 239000002131 composite material Substances 0.000 claims description 17
- 230000006870 function Effects 0.000 description 11
- 230000002194 synthesizing effect Effects 0.000 description 9
- 238000002474 experimental method Methods 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 229920006395 saturated elastomer Polymers 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241000519995 Stachys sylvatica Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
- Studio Devices (AREA)
Description
本発明は、ダイナミックレンジ拡大のための画像合成における露光比決定方法に関する。
明暗の差の大きい被写体をカメラで撮影し、撮影した画像を表示すると、表示された画像に明るい部分の白とびや暗い部分の黒つぶれが見られることがある。これはカメラが撮影可能な輝度の範囲、いわゆるダイナミックレンジが不足していることにより生じている。
従来、ダイナミックレンジを拡大するため、大きくわけてふたつの工夫が行われている。一つは、撮像素子そのものに手を加える手法である。例えば画素の半分を高輝度用、残りの半分を低輝度用とし、双方の画像出力を回路またはメモリ上で合成する手法である。対数増幅器を読出し用CCD(Charge Coupled Devices)の前につけている例もある。
ダイナミックレンジを拡大するためのもう一の工夫は、通常の撮像素子を用いて露光時間の異なる複数の画像を撮影し、これを回路またはメモリ上で合成する手法である。このような画像の合成によってダイナミックレンジを拡大させる場合、各画像の露光時間比は合成された画像の画質を左右する重要なパラメータである。
ここで、画像合成によるダイナミックレンジの拡大について、具体的に説明する。まず、低輝度領域に露光を合わせた画像(この画像を以下、画像Aとする)と、高輝度領域に露光を合わせた画像(この画像を以下、画像Bとする)を撮影する。画像Aの露光時間は、画像Bの露光時間のk倍である。
撮影した画像A、画像Bの両方の画像の入力輝度に対する信号値をグラフにすると図6に示すとおりである。図6に示すように、画像Bは撮影シーンの全輝度範囲に対して飽和のない信号が得られているが、画像Aはある輝度以上になると信号が飽和している。画像Bの入力輝度に対する信号値の傾きをmとすると、画像Aの傾きはk×mである。
そして、画像Aから低輝度部の情報を、画像Bから高輝度部の情報をそれぞれ取り出して合成画像(この画像を以下、画像Cとする)を作成する。この場合、画像Aのうちの飽和している信号部分を、画像Bの信号に画像A、Bの露光時間の比である露光比kをかけ、傾きをk×mとした画像で置き換えて合成することにより、画像Cが得られる。
以上のようにして、全入力輝度に対し、信号が飽和していないダイナミックレンジの広い画像が得られる。
しかしながら、合成された画像Cの元の2つの画像(画像A、B)の境目、すなわち、図6に示した入力輝度X1に相当する領域に、擬似的な輪郭が生じることが多い。これを回避するための従来技術として、例えば特許文献1,2,3に開示されている技術があった。
特開平7−131708号公報
特開平8−214211号公報
特開2004−48445号公報
上記特許文献1には、擬似輪郭発生の原因を画像合成時の2枚の画像の輝度レベル(上記図6での傾きに相当)のずれにあるとし、このずれをなくすために合成する画像信号を解析し、より正確な輝度レベルと合成するときの入力輝度を算出するようにした画像合成装置が開示されている。しかし、この画像合成装置によると、擬似輪郭は減少するものの、まったくなくなるというわけではない。
また、特許文献2には、撮影した画像の信号レベルのヒストグラムを解析し、その結果から2枚目の撮影のためのパラメータ、レンズの絞り値とシャッター時間を決定して撮影するようにした画像合成装置が開示されている。しかし、この画像合成装置では、レンズの絞りが変わるため、1枚目と2枚目の画像の被写界深度が異なってしまい、画像の合成境界付近に不自然さが見られるようになってしまう。また、発生するであろう擬似輪郭については、まったく考慮されていない。
さらに、特許文献3には、画像を滑らかにつなぐために合成境界領域についてローパスフィルター処理を行い、その後、輪郭強調処理を行うようにした画像合成装置および画像合成方法が開示されている。しかし、この画像合成装置および画像合成方法では、高い空間周波数を有する被写体の場合、合成された画像に不自然さが発生するおそれがあった。
以上のように、上述した従来のダイナミックレンジ拡大のための画像合成技術では、合成画像の境界(合成境界)付近における擬似輪郭または不自然さの発生が解消されていないという課題があった。
そこで、本発明は上記課題を解決するためになされたもので、ダイナミックレンジ拡大のための画像合成を行う場合の合成される2つの画像の露光比を決める露光比決定方法において、合成境界付近における擬似輪郭および不自然さの双方の発生を解消できるようにすることを目的とする。
上記課題を解決するため、本発明は、ダイナミックレンジ拡大のために画像合成を行う場合の合成される2つの画像の露光比を決める露光比決定方法であって、高輝度部が飽和しないように露光量を調整して撮影した画像について、各信号の頻度分布を作成し、作成した頻度分布と、ランダムノイズの差を視覚的に認識できる認識限界面積との関係から、2つの画像の合成境界における擬似輪郭の発生を抑制可能な信号範囲を求め、その求めた信号範囲に収まる2つの画像の輝度境界値を与える露光比のうち、各信号に含まれるランダムノイズの均等色空間上での変化量の総和として定義される評価値を最小とする露光比を最適露光比として決定する露光比決定方法を特徴とする。
上記の評価値が小さいほど、画像全体としてランダムノイズによる均等色空間上での変化が少ないことになるので、評価値を最小とする露光比を最適露光比として決定することによって、擬似輪郭が認識され難く、しかも色再現性への影響も小さい露光比が得られる。
また、上記評価値は、ランダムノイズの均等色空間上での変化量を考慮した重み付けを行える合成評価関数を用いて定義されているようにすることができる。
ランダムノイズは、絶対値ではなく基準となる信号値に対する相対的な量として知覚されるため、上記のような合成評価関数を用いることにより、評価値にランダムノイズ特性による均等色空間を考慮した重み付けを行うことができる。
さらに、認識限界面積よりも頻度の低い信号値の範囲を信号範囲とすることができる。
このような信号範囲に含まれる信号値を合成境界に用いると、合成境界部分の擬似輪郭を抑制することが可能となる。
以上詳述したように、本発明によれば、ダイナミックレンジ拡大のための画像合成を行う場合の合成される2つの画像の露光比を決める露光比決定方法において、合成境界付近における擬似輪郭および不自然さの双方の発生を解消できるようになる。
以下、本発明の実施の形態について説明する。なお、同一要素には同一符号を用い、重複する説明は省略する。
ダイナミックレンジ拡大のために画像合成を行う場合において、合成される2つの画像の露光比を大きくとった場合、合成境界を境にランダムノイズ特性が急激に変化する。そのため、連続的に信号が変化する領域に合成境界が存在すると、ランダムノイズ特性の急変により、合成境界に擬似輪郭が発生する可能性がある。このように、合成境界付近に擬似輪郭が発生する原因は、合成境界を境にランダムノイズ特性が急変することにある。
したがって、合成境界付近における擬似輪郭の発生を回避するには、合成境界の前後でランダムノイズが視覚的に認識できる限界(本実施の形態では、これを認識限界面積とする)を求める必要があり、本実施の形態では、この認識限界面積を後述するようにして定量的に把握するようにしている。
そして、本実施の形態では、合成境界付近における擬似輪郭および不自然さの双方の発生を解消可能な露光比で、撮影するシーン(撮影シーン)に適応した最適な露光比を最適露光比といい、最適露光比を決定する露光比決定方法は、以下の手順1〜手順6を有している。
手順1 まず、撮影シーンの高輝度部が飽和しないように露光時間を選んで露光量を調整し、撮像装置で画像を一枚撮影する。この撮像装置は、例えば、単板面順次方式のカラー撮像装置を用いることができる。また、この撮像装置では、3バンドの出力信号値を3×3の変換Matrixを用いてXYZ三刺激値を取得している。
手順2 次に手順1で撮影した画像(高輝度画像ともいい、図6の画像Bに相当する画像)についての撮像装置の出力信号(撮影した画像を構成するための信号)の頻度分布(ヒストグラムHist(Sig))を作成する。
手順3 擬似輪郭の認識限界面積を後述のようにして求める。
手順4 次に、手順2で作成したヒストグラムHist(Sig)に手順3で求めた認識限界面積を重ねた上で、ヒストグラムと認識限界面積との関係から、後述するようにして、合成してもランダムノイズを認識することがなく、合成境界において擬似輪郭の発生を抑制可能な信号の範囲(境界信号範囲という)を求める。
手順5 そして、手順4で求めた境界信号範囲内に収まる合成露光比(図6の入力輝度X1に相当する入力輝度の境界値を与える露光比)を用いて、ランダムノイズによる均等色空間上での変化量を考慮した重み付けを行える後述する合成評価関数ΔEi(Sig)を求める。ここで、iは擬似輪郭の発生しないi番目の露光比を示している。
式(1)に示すように、この評価値sumEiは、合成評価関数ΔEi(Sig)と、合成対象となる画像のヒストグラムHist(Sig)により定義されるものである。合成評価関数は画像合成時におけるランダムノイズの色視覚特性への影響の与え方を示すもので、詳しくは後述するが、ランダムノイズにCIELABの均等色空間(CIELAB色空間)を考慮した重み付けを行ったものである。
式(1)は、各画素の信号に含まれるノイズの、CIELAB色空間上での変化量について、全画素の総和を求めることを意味している。つまり、評価値sumEiが小さいほど、画像全体としてランダムノイズによる均等色空間上での変化が少ないということになる。したがって、この評価値sumEiが小さいほど、画像合成時の均等色空間上でのランダムノイズの影響が小さいため、擬似輪郭が認識され難くなり、しかも色再現性への影響も小さくなるということになる。
したがって、求めた評価値sumEiの最も値が小さくなる露光比を求めれば、その露光比が撮影シーンに適応した最適な露光比(最適露光比)ということになる。
ここで、上述の合成評価関数について説明する。撮影装置のランダムノイズの半値幅σが、信号値Sigに対し、図1のような特性を示したとする。なお、図1では、σは信号値Sigを変数とした直線(図1に示す直線P1)として近似できるとしている。
ランダムノイズは、絶対値ではなく基準となる信号値に対する相対的な量として知覚されるため、このノイズ特性(ランダムノイズ特性)による均等色空間を考慮した重み付けを行う。
ランダムノイズが基準となる信号値に対し十分に小さいとすれば、ランダムノイズが基準となる信号に加わった場合のCIELAB色空間における変化量ΔEはCIELAB色空間への変換式の微分で与えられる。そのため、変化量ΔEはXYZ三刺激値からの変換項である1/3乗項を用いて下記の式(2)で表すような傾向を示すといえる。
上記式(2)より、ランダムノイズの均等色空間上での変化量を考慮した重み付けが可能となる。
次に、式(2)より得られる評価関数を基本形とし、これを基に各露光比を用いた場合の合成画像におけるΔEを作成する。こうして、合成評価関数ΔEi(Sig)を求めることができる。
合成される画像は、高輝度部と低輝度部でそれぞれ別の画像信号を参照するため、ランダムノイズについても合成境界を挟んで、露光量の異なる画像の特性を引き継ぐことになる。そのため、合成により得られる画像のCIELAB色空間におけるランダムノイズ特性は図2に示すようになる。図2に示すように、大きい露光比とした場合ほどノイズによる影響を改善する信号範囲は狭くなるが、CIELAB色空間上での変化量をより抑制することが可能となる。
続いて、合成境界における擬似輪郭の認識限界について説明する。本実施の形態では、合成境界の前後でランダムノイズが視覚的に認識できる限界を調べるにあたり、ディスプレイにおける輝度ムラの認識限界に関する文献1(森由美、田村徹、吉武良治、森口喜代、棚橋高成、辻智、「認識限界コントラストに基づいた輝度ムラの定量化手法」、映像情報メディア学会誌、vol.56,no.11,pp.1837-1840,2002)を参考にした。
文献1より、均一背景中に呈示される輝度ムラの呈示面積の変化に対する認識限界コントラストは、下記の式(3)に示すように、輝度ムラの呈示面積Sの1/3乗に比例して変化する。ここで、Lは背景輝度を示し、ΔLはLに対する輝度ムラの差を示している。a,bは主観評価実験(具体的には後述する)より求められる定数である。
この式(3)に対し、本実施の形態では、輝度ムラの差ΔLをランダムノイズの半値幅σに置き換え、呈示面積Sを合成境界信号の頻度と置き換えても同様の関係式が得られると仮定している。
そして、背景輝度の変化に対し、認識限界コントラストがウェーバー則(基礎刺激Lによらず刺激の弁別閾はΔL/L=一定)に従い一定であるとすれば、式(3)から式(4)が得られる。式(4)より、各基準信号における擬似輪郭の認識限界面積Sが求められる。
式(4)によって認識限界面積を求めると、手順2で求めた撮影シーンのヒストグラムHist(Sig)を用い次のようにして境界信号範囲を求めることができる。例えば、ヒストグラムHist(Sig)に認識限界面積を重ねた結果、図3に示すようになったとする。
この場合、認識限界面積よりも頻度の低い信号値(Sig0.1以上の信号値)を合成境界に用いるとすると、合成境界部分の擬似輪郭を抑制することが可能となるから、この場合はSig0.1以上の信号値が境界信号範囲になる。
次に、認識限界面積を求めるための主観評価実験について説明する。主観評価実験は、例えば以下の実験条件によって行うことができる。
Display : 21インチ、広色域LCD(Liquid Crystal Display)
画素数 : 1600×1200
視距離 : 1m
水平面照度 : 124lx(室内照明下)
被験者 : 5名
Display : 21インチ、広色域LCD(Liquid Crystal Display)
画素数 : 1600×1200
視距離 : 1m
水平面照度 : 124lx(室内照明下)
被験者 : 5名
使用するノイズ画像は、信号値Sigに半値幅σのランダムノイズを加え、マトリクス変換とガンマ変換(γ=2.2)によりRGB階調値に変換して作成した。
評価には調整法1(大串健吾,他,「画質と音質の評価技術」、昭晃堂,1991)を用い、10種のノイズ面積についてノイズの分散量を徐々に上げる、または下げる方向に調整しながら、各被験者にノイズがちょうど認識できなくなるレベルを指示させた。また、実験は各条件の評価を3回ずつ試行し、その平均を測定値とした.
この実験結果は図4に示すとおりである。図4では、プロット点が測定点を表し、上述した文献1に従い、直線近似を行っている。図4に示す近似直線P2より、式(4)のa,bが決まる。そのため、決定されたa,bを用い、式(4)により、認識限界面積を求めることができる。
そして、本発明の有効性を確認するため、図示しない撮影シーンに対して、本実施の形態の露光比決定方法により最適露光比を決定し、その最適露光比にしたがって画像の合成を行った。各露光比の評価結果は図5に示すとおりである。図5に示す数値は各バンド(S1、S2、S3)におけるsumEiを表している。
図5において、各バンドにおける最小値を表示している欄に網掛けを施し、擬似輪郭が発生すると判断された数値には取り消し線を表示している。
sumEiの最小値はバンドS1、S2、S3毎にそれぞれ64倍、32倍、64倍の露光比を指している。ところが、これらは上述した主観評価実験の結果で擬似輪郭が発生すると判断されたために、取り消し線の表示されていない範囲での最小値である16倍、16倍、8倍が撮影シーンに対する最適露光比として決定されている。
こうして求められた最適露光比を用いて画像の合成を行うと、図示はしないが、原画像に対して、低輝度部でノイズが少なく、しかも合成境界付近の擬似輪郭が無い自然で鮮明な合成画像を得ることができる。
以上詳述したように、本実施の形態における露光比決定方法により決定した最適露光比を用いて画像の合成を行うと、原画像に対し低輝度部でノイズが少なく、しかも合成境界の擬似輪郭が無く自然で鮮明な合成画像を得ることができる。
また、上記露光比決定方法はランダムノイズという物理現象と人の平均的な視覚特性に基づいているため汎用性の高いものである。
また、上記露光比決定方法はランダムノイズという物理現象と人の平均的な視覚特性に基づいているため汎用性の高いものである。
以上の説明は、本発明の実施の形態についての説明であって、この発明を限定するものではなく、様々な変形例を容易に実施することができる。又、各実施形態における構成要素、機能、特徴あるいは方法ステップを適宜組み合わせて構成される装置又は方法も本発明に含まれるものである。
P1…直線、P2…近似直線
Claims (4)
- ダイナミックレンジ拡大のために画像合成を行う場合の合成される2つの画像の露光比を決める露光比決定方法であって、
高輝度部が飽和しないように露光量を調整して撮影した画像について、各信号の頻度分布を作成し、
前記作成した頻度分布と、ランダムノイズの差を視覚的に認識できる認識限界面積との関係から、前記2つの画像の合成境界における擬似輪郭の発生を抑制可能な信号範囲を求め、
その求めた前記信号範囲に収まる前記2つの画像の輝度境界値を与える露光比のうち、前記各信号に含まれるランダムノイズの均等色空間上での変化量の総和として定義される評価値を最小とする露光比を最適露光比として決定する露光比決定方法。 - 前記評価値は、ランダムノイズによる均等色空間上での変化量を考慮した重み付けを行える合成評価関数を用いて定義されている請求項1記載の露光比決定方法。
- 前記認識限界面積よりも頻度の低い信号値の範囲を前記信号範囲とする請求項1または2記載の露光比決定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007041053A JP4802333B2 (ja) | 2007-02-21 | 2007-02-21 | 画像合成における露光比決定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007041053A JP4802333B2 (ja) | 2007-02-21 | 2007-02-21 | 画像合成における露光比決定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2008205935A JP2008205935A (ja) | 2008-09-04 |
JP4802333B2 true JP4802333B2 (ja) | 2011-10-26 |
Family
ID=39782916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007041053A Active JP4802333B2 (ja) | 2007-02-21 | 2007-02-21 | 画像合成における露光比決定方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4802333B2 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6377258B2 (ja) | 2015-04-16 | 2018-08-22 | 三菱電機株式会社 | 画像処理装置及び方法、並びにプログラム及び記録媒体 |
GB2569593B (en) * | 2017-12-20 | 2021-05-26 | Apical Ltd | Exposure ratio control |
JP6810299B2 (ja) * | 2018-04-26 | 2021-01-06 | 富士フイルム株式会社 | 画像処理装置、方法、及びプログラム並びに撮像装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07131718A (ja) * | 1993-10-29 | 1995-05-19 | Canon Inc | 画像合成装置 |
JP2002165138A (ja) * | 2000-11-24 | 2002-06-07 | Olympus Optical Co Ltd | 撮像装置 |
JP4097980B2 (ja) * | 2002-04-23 | 2008-06-11 | オリンパス株式会社 | 画像合成装置 |
-
2007
- 2007-02-21 JP JP2007041053A patent/JP4802333B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2008205935A (ja) | 2008-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4218723B2 (ja) | 画像処理装置、撮像装置、画像処理方法およびプログラム | |
EP2624204B1 (en) | Image processing apparatus and method of controlling the same | |
KR101247646B1 (ko) | 화상 처리 장치, 화상 처리 방법, 및 기록 매체 | |
JP5083046B2 (ja) | 撮像装置及び撮像方法 | |
JP4208909B2 (ja) | 画像処理装置と撮影装置 | |
US7791656B2 (en) | Image sensing apparatus and image processing method | |
JP4645921B2 (ja) | 画像信号処理装置および方法、プログラム、並びに撮像装置 | |
JP4946795B2 (ja) | 画像処理装置及び画像処理方法 | |
JP2008104009A (ja) | 撮像装置および撮像方法 | |
JP6097588B2 (ja) | 画像処理装置及び画像処理方法 | |
JP5699482B2 (ja) | 画像処理装置、画像処理方法及び撮像装置 | |
US8085316B2 (en) | Image-processing unit, imaging apparatus, and computer program product | |
KR20090028040A (ko) | 노출 시간 설정 장치 및 방법 | |
JP2007082181A (ja) | 撮像装置及び画像処理方法 | |
US10218953B2 (en) | Image processing apparatus, image processing method, and storage medium | |
KR100933556B1 (ko) | 다이내믹 레인지를 확장하는 칼라 영상 처리장치 및 방법 | |
JP4802333B2 (ja) | 画像合成における露光比決定方法 | |
JP2011100204A (ja) | 画像処理装置、画像処理方法、画像処理プログラム、撮像装置及び電子機器 | |
JP5932392B2 (ja) | 画像処理装置および画像処理方法 | |
JP2016201733A (ja) | 撮像装置、撮像方法、並びにプログラム | |
JP2011199785A (ja) | 撮像装置及び方法 | |
KR101923162B1 (ko) | 액정 패널을 이용한 hdri 영상 획득 장치 및 방법 | |
EP1976273B1 (en) | Image processing apparatus, imaging apparatus and image processing program | |
JP2008219230A (ja) | 撮像装置及び画像処理方法 | |
JP7346077B2 (ja) | 画像処理装置、画像処理方法及びプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100209 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20110704 |
|
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: 20110712 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |