JPH0516231B2 - - Google Patents
Info
- Publication number
- JPH0516231B2 JPH0516231B2 JP58157832A JP15783283A JPH0516231B2 JP H0516231 B2 JPH0516231 B2 JP H0516231B2 JP 58157832 A JP58157832 A JP 58157832A JP 15783283 A JP15783283 A JP 15783283A JP H0516231 B2 JPH0516231 B2 JP H0516231B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- linear array
- detector
- intensity
- statistics
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- 238000004364 calculation method Methods 0.000 claims description 16
- 230000035945 sensitivity Effects 0.000 claims description 11
- 238000000034 method Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 description 7
- 238000003384 imaging method Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000003702 image correction Methods 0.000 description 2
- 229920006395 saturated elastomer Polymers 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N25/00—Circuitry of solid-state image sensors [SSIS]; Control thereof
- H04N25/60—Noise processing, e.g. detecting, correcting, reducing or removing noise
- H04N25/67—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response
- H04N25/671—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response for non-uniformity detection or correction
- H04N25/677—Noise processing, e.g. detecting, correcting, reducing or removing noise applied to fixed-pattern noise, e.g. non-uniformity of response for non-uniformity detection or correction for reducing the column or line fixed pattern noise
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
Description
【発明の詳細な説明】
〔発明の利用分野〕
本発明は、リニアアレイセンサにより撮像され
た画像データの補正方式にかかわり、特にリニア
アレイセンサの画像検出器間相互の感度のばらつ
きを補正する好適な補正方式に関する。
た画像データの補正方式にかかわり、特にリニア
アレイセンサの画像検出器間相互の感度のばらつ
きを補正する好適な補正方式に関する。
第１図はリニアアレイセンサによる地表面撮像
の様子を示す。
の様子を示す。
第１図において、リニアアレイセンサ１０は、
多数の半導体検出器が直線状に並んだものであ
る。リニアアレイセンサと撮像対象１１を、相対
的に矢印１２の方向に移動させつつレンズ１３を
通して撮像し、画像検出器の出力を順次サンプル
すると、撮像対象の幅Ｗにわたる画像データを得
ることができる。
多数の半導体検出器が直線状に並んだものであ
る。リニアアレイセンサと撮像対象１１を、相対
的に矢印１２の方向に移動させつつレンズ１３を
通して撮像し、画像検出器の出力を順次サンプル
すると、撮像対象の幅Ｗにわたる画像データを得
ることができる。
一般に、リニアアレイセンサは、製造上の問題
から各検出器間には特性のばらつきがある。その
ため、人工衛星に搭載して地球を撮像するような
場合には、画像上に濃淡の縞模様が生じ画質を低
下させる。この例を模式的に第２図に示す。第２
図において、横軸ｉは画像検出器の番号を示し、
縦軸ｊは走査ライン番号を示す。
から各検出器間には特性のばらつきがある。その
ため、人工衛星に搭載して地球を撮像するような
場合には、画像上に濃淡の縞模様が生じ画質を低
下させる。この例を模式的に第２図に示す。第２
図において、横軸ｉは画像検出器の番号を示し、
縦軸ｊは走査ライン番号を示す。
複数検出器（これを＃１〜＃ｎとする）から成
るセンサの強度補正には絶対補正と相対補正があ
る。絶対補正とは、各検出器独立に行なうもの
で、較正用基準光源を用いて第３図のような入出
力関係を測定し、理想的特性に合わせて、それぞ
れの出力値を変換するものである。本補正は一般
的なものであり、リニアアレイセンサについても
適用される。
るセンサの強度補正には絶対補正と相対補正があ
る。絶対補正とは、各検出器独立に行なうもの
で、較正用基準光源を用いて第３図のような入出
力関係を測定し、理想的特性に合わせて、それぞ
れの出力値を変換するものである。本補正は一般
的なものであり、リニアアレイセンサについても
適用される。
相対補正とは、上記絶対補正によつても残る入
出力特性のばらつき、および特性の経時変化を補
正し、わずかな縞模様も除去しようとするもので
ある。
出力特性のばらつき、および特性の経時変化を補
正し、わずかな縞模様も除去しようとするもので
ある。
従来、LANDSAT衛星MSS（Multi Spectral
Scanner）センサの場合には、各検出器ごとの撮
像データの頻度分布を合わせこむ手法がとられて
きた。６個の検出器で６ライン平行に地表を走査
するため、どの検出器も一様に地上データをサン
プルしていることになり、入力強度の頻度分布は
統計的に等しいと言えるからである。ところが、
リニアアレイセンサでは、各検出器が走査する画
面上の位置は異なつているため、頻度分布を合わ
せ込んでも、対象とする地域依存した補正誤差を
生ずるという欠点があつた。
Scanner）センサの場合には、各検出器ごとの撮
像データの頻度分布を合わせこむ手法がとられて
きた。６個の検出器で６ライン平行に地表を走査
するため、どの検出器も一様に地上データをサン
プルしていることになり、入力強度の頻度分布は
統計的に等しいと言えるからである。ところが、
リニアアレイセンサでは、各検出器が走査する画
面上の位置は異なつているため、頻度分布を合わ
せ込んでも、対象とする地域依存した補正誤差を
生ずるという欠点があつた。
本発明は上記欠点を取り除くためになされたも
ので、その目的とするところは、リニアアレイセ
ンサ画像の検出器間の感度特性の微少なばらつき
に起因する画像上の縞模様ノイズを除去補正する
方式を提供することにある。
ので、その目的とするところは、リニアアレイセ
ンサ画像の検出器間の感度特性の微少なばらつき
に起因する画像上の縞模様ノイズを除去補正する
方式を提供することにある。
上記目的を達成するため本発明においては、リ
ニアアレイセンサの各画像検出器ごとに出力画像
強度の統計量を算出し、各画像検出器別統計量か
ら、それが隣接検出器間では急激に変化しないこ
とを前提に各画像検出器ごとの強度補正量を求
め、求められた強度補正量に応じてリニアアレイ
センサ撮像画像を補正する点に特徴がある。
ニアアレイセンサの各画像検出器ごとに出力画像
強度の統計量を算出し、各画像検出器別統計量か
ら、それが隣接検出器間では急激に変化しないこ
とを前提に各画像検出器ごとの強度補正量を求
め、求められた強度補正量に応じてリニアアレイ
センサ撮像画像を補正する点に特徴がある。
まず、本発明の原理を詳しく述べる。
本発明による補正方式では、まず検出器別に実
際の地表撮像データの統計量を充分な量のデータ
について求め、次に隣接検出器の統計量は本来急
激に変化しないという仮定から検出器の感度相対
誤差を求めようとするものである。統計量として
は、平均値と標準偏差を用いる。
際の地表撮像データの統計量を充分な量のデータ
について求め、次に隣接検出器の統計量は本来急
激に変化しないという仮定から検出器の感度相対
誤差を求めようとするものである。統計量として
は、平均値と標準偏差を用いる。
第２図に示すように、ｉ番目の検出器に関する
ｊライン目の真の画像データをｘ（ｉ，ｊ）とす
る。強度補正を行う前には、画像データx′（ｉ，
ｊ）は検出器ごとの感度，オフセツトの誤差を持
つている。
ｊライン目の真の画像データをｘ（ｉ，ｊ）とす
る。強度補正を行う前には、画像データx′（ｉ，
ｊ）は検出器ごとの感度，オフセツトの誤差を持
つている。
x′（ｉ，ｊ）＝α(i)・ｘ（ｉ，ｊ）＋β(i)
α(i)，β(i)０ (1)
従つて、充分なライン数Ｎに関する各検出器ご
との画像データの平均値m′(i)，標準偏差σ′(i)は、 となる。ここでｍ(i)，σ(i)はそれぞれ真の画像デ
ータの平均と分散である。LANDSATによる実
際の画像データの解析結果などからｍ(i)，α(i)は
検出器ｉに関して急激には変化しないと考えられ
る。そこで、画像データx′（ｉ，ｊ）の平均
m′(i)，標準偏差σ′(i)の急激な変化はすべて検出器
の感度相対誤差α(i)，β(i)によるものと考え、真
の値ｍ(i)，α(i)は下記低域フイルタで推定を行な
う。
との画像データの平均値m′(i)，標準偏差σ′(i)は、 となる。ここでｍ(i)，σ(i)はそれぞれ真の画像デ
ータの平均と分散である。LANDSATによる実
際の画像データの解析結果などからｍ(i)，α(i)は
検出器ｉに関して急激には変化しないと考えられ
る。そこで、画像データx′（ｉ，ｊ）の平均
m′(i)，標準偏差σ′(i)の急激な変化はすべて検出器
の感度相対誤差α(i)，β(i)によるものと考え、真
の値ｍ(i)，α(i)は下記低域フイルタで推定を行な
う。
ここで、vk，wkは荷重であり、p
〓k=-p
vk＝１，p
〓k=
−ｐ wk＝１である。
−ｐ wk＝１である。
式(3)を式(2)のｍ(i)，σ(i)に代入しα(i)，β(i)を
求めれば、式(1)より真の画像データｘ（ｉ，ｊ）
をつぎのように求めることができる。
求めれば、式(1)より真の画像データｘ（ｉ，ｊ）
をつぎのように求めることができる。
ｘ（ｉ，ｊ）＝１／α(i)｛x′（ｉ，ｊ）−β(i)
α(i)＝σ′(i)／σ／＾(i)，β(i)＝m′(i)−σ′(i
)／σ／＾(i)・ｍ(i) ……(4) 以上、本発明の基本的原理を示したが、画像検
出器のダイナミツクレンジを越えるデータが入力
した場合には、このような基本原理だけでは正し
く補正を行なえない場合が生ずる。すなわち、雲
や雪などの明るいデータ、水面などの暗いデータ
が検出器ダイナミツクレンジの不足から上下限値
に飽和した場合には、真値から誤差を持つ平均，
分散が求まることになる。しかもこの誤差は、検
出器の感度に依存するため、前記基本原理では画
像の正しい補正ができないことになる。
)／σ／＾(i)・ｍ(i) ……(4) 以上、本発明の基本的原理を示したが、画像検
出器のダイナミツクレンジを越えるデータが入力
した場合には、このような基本原理だけでは正し
く補正を行なえない場合が生ずる。すなわち、雲
や雪などの明るいデータ、水面などの暗いデータ
が検出器ダイナミツクレンジの不足から上下限値
に飽和した場合には、真値から誤差を持つ平均，
分散が求まることになる。しかもこの誤差は、検
出器の感度に依存するため、前記基本原理では画
像の正しい補正ができないことになる。
ここでは、上記問題を解決するため、飽和値を
平均，分散の算出から除く方式を示す。第４図
は、各々の検出器（＃１〜＃Ｎ）に関する撮像デ
ータのヒストグラムであり、それぞれ上下限値で
異なる飽和の様子を示している。今、式(3)のフイ
ルタ処理の区間〔ｉ−ｐ，ｉ＋ｐ〕を考え、フイ
ルタ処理で入力する平均m′（ｉ＋ｋ），分散σ′2（ｉ
＋ｋ）｛ｋ＝−ｐ〜ｐ｝を、飽和値のないヒスト
グラムから算出したものに修正することを考え
る。手順は、 (1) 区間〔ｉ−ｐ，ｉ＋ｐ〕で下限飽和値の個数
のもつとも多いものを見出し、その個数分のデ
ータを区間内の各ヒストグラムの下限値の方か
ら取除いていく。
平均，分散の算出から除く方式を示す。第４図
は、各々の検出器（＃１〜＃Ｎ）に関する撮像デ
ータのヒストグラムであり、それぞれ上下限値で
異なる飽和の様子を示している。今、式(3)のフイ
ルタ処理の区間〔ｉ−ｐ，ｉ＋ｐ〕を考え、フイ
ルタ処理で入力する平均m′（ｉ＋ｋ），分散σ′2（ｉ
＋ｋ）｛ｋ＝−ｐ〜ｐ｝を、飽和値のないヒスト
グラムから算出したものに修正することを考え
る。手順は、 (1) 区間〔ｉ−ｐ，ｉ＋ｐ〕で下限飽和値の個数
のもつとも多いものを見出し、その個数分のデ
ータを区間内の各ヒストグラムの下限値の方か
ら取除いていく。
(2) 同様の操作を上限飽和値について実行する。
修正後のヒストグラム（第４図ｂ）の平均，分
散により式(3)，(4)を計算すれば、検出器感度の
ばらつきと飽和値の存在による画像補正誤差を
除くことができる。
修正後のヒストグラム（第４図ｂ）の平均，分
散により式(3)，(4)を計算すれば、検出器感度の
ばらつきと飽和値の存在による画像補正誤差を
除くことができる。
以下、本発明の実施例を第５図により説明す
る。衛星１に搭載されたリニアアレイセンサによ
る地球観測データは、アンテナ２により受信され
衛星の上空通過時にリアルタイムで記録システム
３により高密度磁気テープ４に記録される。高密
度磁気テープ４は画像処理システムにおいて強度
的歪みと幾何学的歪みが除かれ補正画像データ５
として出力される。画像処理システムでは、まず
高密度磁気テープ４の未補正画像データは高密度
磁気テープ装置６により読込まれ、磁気デイスク
装置９に格納されると同時に演算装置７によりリ
ニアアレイセンサの各検出器ごとの統計量が計算
される。統計量としては最終的に平均値と分散値
が用いられるが演算に要するメモリ容量を少なく
するため、逐次型の演算とする。すなわち記憶装
置８には、第ｒライン目までの画像加算値r 〓j=1
x′（ｉ，ｊ）自乗加算値r 〓j=1 x′（ｉ，ｊ）2のみを記憶
し、高密度磁気テープ４から１ライン入力する度
に下記演算をおこなう。
る。衛星１に搭載されたリニアアレイセンサによ
る地球観測データは、アンテナ２により受信され
衛星の上空通過時にリアルタイムで記録システム
３により高密度磁気テープ４に記録される。高密
度磁気テープ４は画像処理システムにおいて強度
的歪みと幾何学的歪みが除かれ補正画像データ５
として出力される。画像処理システムでは、まず
高密度磁気テープ４の未補正画像データは高密度
磁気テープ装置６により読込まれ、磁気デイスク
装置９に格納されると同時に演算装置７によりリ
ニアアレイセンサの各検出器ごとの統計量が計算
される。統計量としては最終的に平均値と分散値
が用いられるが演算に要するメモリ容量を少なく
するため、逐次型の演算とする。すなわち記憶装
置８には、第ｒライン目までの画像加算値r 〓j=1
x′（ｉ，ｊ）自乗加算値r 〓j=1 x′（ｉ，ｊ）2のみを記憶
し、高密度磁気テープ４から１ライン入力する度
に下記演算をおこなう。
r+1
〓j=1
x′（ｉ，ｊ）＝r
〓j=1
x′（ｉ，ｊ）＋x′（ｉ，ｒ＋１）
r+1
〓j=1
x′（ｉ，ｊ）2＝r
〓j=1
x′（ｉ，ｊ）2＋x′（ｉ，ｒ＋１）2……(5)
したがつて統計量記憶装置８の記憶容量は、検
出器数を4000とすると、32KByteですむことに
なる。第６図に検出器別統計量を求める演算装置
７の構成を示す。今、第ｒ＋１番目の画像データ
列x′（ｉ，ｒ＋１）（ｉ＝１，……，Ｍ：Ｍ画素
数）が装置７に転送されると、乗算器７１と加算
器７３に逐次入力される。乗算器７１では画像値
の自乗演算がなされ、加算器７２によりデータメ
モリ８より読み出された第ｒ番目までの自乗加算
値列r 〓j=1 x′（ｉ，ｊ）2と順次加算され、結果は第ｒ
＋１番目までの自乗加算値r+1 〓j=1 x′（ｉ，ｊ）2として
再度データメモリ８に格納される。加算器７３で
は、データメモリ８より読み出された第ｒライン
目までの画像加算値列r 〓j=1 x′（ｉ，ｊ）と第ｒ＋１
番目の画像データ列x′（ｉ，ｒ＋１）が順次加算
され、結果は再度データメモリ８に格納される。
磁気デイスク装置９に全画像データが格納された
段階で、記憶装置８上の各検出器ごとの画像加算
値と自乗加算値は、強度補正量演算装置１０によ
り、下式で平均値m′(i)，標準偏差値σ′(i)に変換さ
れる。
出器数を4000とすると、32KByteですむことに
なる。第６図に検出器別統計量を求める演算装置
７の構成を示す。今、第ｒ＋１番目の画像データ
列x′（ｉ，ｒ＋１）（ｉ＝１，……，Ｍ：Ｍ画素
数）が装置７に転送されると、乗算器７１と加算
器７３に逐次入力される。乗算器７１では画像値
の自乗演算がなされ、加算器７２によりデータメ
モリ８より読み出された第ｒ番目までの自乗加算
値列r 〓j=1 x′（ｉ，ｊ）2と順次加算され、結果は第ｒ
＋１番目までの自乗加算値r+1 〓j=1 x′（ｉ，ｊ）2として
再度データメモリ８に格納される。加算器７３で
は、データメモリ８より読み出された第ｒライン
目までの画像加算値列r 〓j=1 x′（ｉ，ｊ）と第ｒ＋１
番目の画像データ列x′（ｉ，ｒ＋１）が順次加算
され、結果は再度データメモリ８に格納される。
磁気デイスク装置９に全画像データが格納された
段階で、記憶装置８上の各検出器ごとの画像加算
値と自乗加算値は、強度補正量演算装置１０によ
り、下式で平均値m′(i)，標準偏差値σ′(i)に変換さ
れる。
m′(i)＝１／ＮN
〓j=1
x′（ｉ，ｊ）σ′(i)＝｛１／ＮN
〓j=1
x′（ｉ，ｊ）2−m′(i)2｝ ……(6)
ｍ(i)とσ′(i)を用いて式(4)より計算された１／α(
i) とβ(i)が強度補正係数記憶装置１１に格納され
る。ここでＮは、入力した画像データのライン数
である。強度補正量演算装置１０は、装置７，１
２に比べ、画素ごとの処理でないため高速性は要
求されず、ミニコンピユータで実現できる。装置
１０の演算フローを第７図に示す。演算フロー
は、検出器ごとの画像平均値，標準偏差の算出ス
テツプ71、真の平均値，標準偏差の推定ステツプ
72、補正係数の算出ステツプ73からなる。画像強
度補正装置１２は、画像データを磁気デイスク装
置から１ラインずつ読出し、記憶装置１１上の強
度補正係数により、式(4)の強度補正演算を行う。
i) とβ(i)が強度補正係数記憶装置１１に格納され
る。ここでＮは、入力した画像データのライン数
である。強度補正量演算装置１０は、装置７，１
２に比べ、画素ごとの処理でないため高速性は要
求されず、ミニコンピユータで実現できる。装置
１０の演算フローを第７図に示す。演算フロー
は、検出器ごとの画像平均値，標準偏差の算出ス
テツプ71、真の平均値，標準偏差の推定ステツプ
72、補正係数の算出ステツプ73からなる。画像強
度補正装置１２は、画像データを磁気デイスク装
置から１ラインずつ読出し、記憶装置１１上の強
度補正係数により、式(4)の強度補正演算を行う。
第８図に示すように、第ｒ番目の画像データ列
x′（ｉ，ｒ）は加算器１２０により、検出器対応
のオフセツト補正係数β(i)を減ぜられ、続いて乗
算器１２１によりゲイン係数１／α(i)を乗ぜられ、 パイプライン的に出力される。強度補正済画像は
磁気テープ装置１３により磁気テープ５に出力さ
れる。
x′（ｉ，ｒ）は加算器１２０により、検出器対応
のオフセツト補正係数β(i)を減ぜられ、続いて乗
算器１２１によりゲイン係数１／α(i)を乗ぜられ、 パイプライン的に出力される。強度補正済画像は
磁気テープ装置１３により磁気テープ５に出力さ
れる。
以上の実施例において、入力画像データと検出
器のダイナミツクレンジの関係から飽和値を生ず
る場合には、検出器ごとの画像加算値と自乗加算
値ばかりでなく、上下限値付近の頻度値も求める
必要がある。従つて、検出器別統計量演算装置７
では上下限値付近の頻度の算出も行ない、統計量
記憶装置８は上下限値付近の頻度値も記憶する。
また、画像強度補正量演算装置１０では、式(6)の
算出を次のように修正して行なう。ただし、式(3)
のフイルタ処理区間〔ｉ−ｐ，ｉ＋ｐ〕中の最大
飽和値を上下限値xp，xnそれぞれに対してN0 i，
Ni nとし、各検出器の算出頻度をｈ(x)とする。
器のダイナミツクレンジの関係から飽和値を生ず
る場合には、検出器ごとの画像加算値と自乗加算
値ばかりでなく、上下限値付近の頻度値も求める
必要がある。従つて、検出器別統計量演算装置７
では上下限値付近の頻度の算出も行ない、統計量
記憶装置８は上下限値付近の頻度値も記憶する。
また、画像強度補正量演算装置１０では、式(6)の
算出を次のように修正して行なう。ただし、式(3)
のフイルタ処理区間〔ｉ−ｐ，ｉ＋ｐ〕中の最大
飽和値を上下限値xp，xnそれぞれに対してN0 i，
Ni nとし、各検出器の算出頻度をｈ(x)とする。
m′(i)＝１／Ｎ｛N
〓j=1
x′（ｉ，ｊ）−l0
〓l=0
（x0＋ｌ）・ｈ（x0＋ｌ）−ln
〓l=0
（xn−ｌ）・ｈ（xn−ｌ）｝
σ′2(i)＝１／Ｎ｛N
〓j=1
x′（ｉ，ｊ）2−l0
〓l=0
（x0＋ｌ）2・ｈ（x0＋ｌ）−ln
〓l=0
（xn−ｌ）2h（xn−ｌ）｝−m′(i)2
ただし、N0 i＝lp
〓l=0
ｈ（x0＋ｌ），Ni n＝ln
〓l=0
ｈ（xn−ｌ）
以下、補正係数α(i)，β(i)の算出手段、画像デ
ータの補正手段は、前記基本的実施例と同一であ
る。
ータの補正手段は、前記基本的実施例と同一であ
る。
本発明によれば、リニアアレイセンサの検出器
間の感度特性の微少なばらつき、経時変化に起因
する画像上の縞模様ノイズを除去補正できる効果
がある。
間の感度特性の微少なばらつき、経時変化に起因
する画像上の縞模様ノイズを除去補正できる効果
がある。
第１図はリニアアレイセンサによる地表面撮像
の様子を示す図、第２図は撮像された画像中の縞
パタンを示す図、第３図は各検出器の入出力特性
のばらつきの様子を示す図、第４図は各検出器に
よる撮像データのヒストグラム、第５図は本発明
の方式による衛星搭載リニアアレイセンサの地上
局画像補正システムの構成図、第６図は検出器別
の統計量を求める演算装置の構成図、第７図は強
度補正量演算装置における演算のフローチヤー
ト、第８図は画像強度補正装置の構成図である。 １……リニアアレイセンサ搭載衛星。
の様子を示す図、第２図は撮像された画像中の縞
パタンを示す図、第３図は各検出器の入出力特性
のばらつきの様子を示す図、第４図は各検出器に
よる撮像データのヒストグラム、第５図は本発明
の方式による衛星搭載リニアアレイセンサの地上
局画像補正システムの構成図、第６図は検出器別
の統計量を求める演算装置の構成図、第７図は強
度補正量演算装置における演算のフローチヤー
ト、第８図は画像強度補正装置の構成図である。 １……リニアアレイセンサ搭載衛星。
Claims (1)
- 【特許請求の範囲】 １ リニアアレイセンサの画像検出器ごとに所定
数の走査ラインにおける出力画像強度の統計量を
算出するステツプと、算出された各画像検出器別
統計量から該統計量が隣接検出器間では急激に変
化しないことを前提に画像検出器ごとの強度補正
量を算出するステツプと、算出された強度補正量
に応じてリニアアレイセンサの出力画像を補正す
るステツプとからなることを特徴とするリニアア
レイセンサ素子間感度補正方式。 ２ 上記統計量を算出するステツプは画像検出器
ごとの出力画像強度と該強度の自乗を上記走査ラ
インにわたつて累加する演算からなることを特徴
とする特許請求の範囲第１項のリニアアレイセン
サ素子間感度補正方式。 ３ 上記統計量を算出するステツプは画像検出器
ごとの上下限値画素数および上下限値付近の画像
値区間頻度を計算する演算を含むことを特徴とす
る特許請求の範囲第１項のリニアアレイセンサ素
子間感度補正方式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58157832A JPS6051377A (ja) | 1983-08-31 | 1983-08-31 | リニアアレイセンサ素子間感度補正方式 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58157832A JPS6051377A (ja) | 1983-08-31 | 1983-08-31 | リニアアレイセンサ素子間感度補正方式 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS6051377A JPS6051377A (ja) | 1985-03-22 |
| JPH0516231B2 true JPH0516231B2 (ja) | 1993-03-03 |
Family
ID=15658301
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP58157832A Granted JPS6051377A (ja) | 1983-08-31 | 1983-08-31 | リニアアレイセンサ素子間感度補正方式 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS6051377A (ja) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007028007A2 (en) | 2005-08-31 | 2007-03-08 | Venetec International, Inc. | Anchoring system for a catheter |
| US7879013B2 (en) | 2005-12-21 | 2011-02-01 | Venetec International, Inc. | Intravenous catheter anchoring device |
| US8218868B2 (en) * | 2007-06-13 | 2012-07-10 | Sensors Unlimited, Inc. | Method and apparatus for enhancing images |
| US8394067B2 (en) | 2009-05-21 | 2013-03-12 | C.R. Bard, Inc. | Medical device securement system |
| US10537714B2 (en) | 2009-11-11 | 2020-01-21 | Venetec International, Inc. | Stabilizing device for an extension set |
| CN110585567B (zh) | 2013-03-15 | 2022-04-12 | 维尼泰克国际股份有限公司 | 具有整体搭接件和敷料的固定设备 |
| WO2020160318A1 (en) | 2019-02-01 | 2020-08-06 | Becton, Dickinson And Company | Stabilization device, system, and methods thereof for integrated catheters |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS57152268A (en) * | 1981-03-14 | 1982-09-20 | Nec Corp | Reference information reading system |
-
1983
- 1983-08-31 JP JP58157832A patent/JPS6051377A/ja active Granted
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS57152268A (en) * | 1981-03-14 | 1982-09-20 | Nec Corp | Reference information reading system |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS6051377A (ja) | 1985-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Scribner et al. | Adaptive nonuniformity correction for IR focal-plane arrays using neural networks | |
| CN101873428B (zh) | 图像处理装置及图像处理方法 | |
| US6646246B1 (en) | Method and system of noise removal for a sparsely sampled extended dynamic range image sensing device | |
| US20150116539A1 (en) | Apparatus and method for high dynamic range imaging using spatially varying exposures | |
| US20080043115A1 (en) | Image Processing Device and Image Processing Program | |
| CN113763271B (zh) | 一种基于物理噪声模型的高质量光谱去噪方法 | |
| JPH04330872A (ja) | 白色シェーディング補正装置および補正方法 | |
| CN112529807B (zh) | 卫星影像的相对辐射校正方法和装置 | |
| JP2004005675A (ja) | 勾配解析を利用したデジタル画像のノイズ予測方法 | |
| JP2002525722A (ja) | 画像処理方法とシステム | |
| CN110782429B (zh) | 一种基于星载遥感相机的成像质量评估方法 | |
| Tanbakuchi et al. | Adaptive pixel defect correction | |
| JPH0516231B2 (ja) | ||
| CN112200848A (zh) | 低光照弱对比复杂环境下的深度相机视觉增强方法及系统 | |
| Cao et al. | Vignetting correction based on a two-dimensional Gaussian filter with harmony for area array sensors | |
| CN114519676A (zh) | 基于Bayer格式的raw图像的去噪装置及方法 | |
| CN116109491A (zh) | 利用明暗均匀区域进行成像光谱仪图像非均匀性校正方法 | |
| Jonckheere et al. | Derivative analysis for in situ high dynamic range hemispherical photography and its application in forest stands | |
| Koren | Measuring camera information capacity with slanted-edges | |
| CN113436096A (zh) | 一种基于像元标定的推扫式高光谱成像条带噪声消除方法 | |
| Balovsyak et al. | Hardware and Software Complex for Automatic Level Estimation and Removal of Gaussian Noise in Images | |
| JP2590095B2 (ja) | 衛星画像撮影装置および衛星画像補正方法 | |
| JP6552520B2 (ja) | デジタル撮像装置を較正するための較正方法および該デジタル撮像装置 | |
| JP3531256B2 (ja) | 赤外線検知素子のバラツキの自動補正方法 | |
| JP2874872B2 (ja) | 監視装置 |