Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JPH0691609B2 - Linear array sensor abnormality detection method - Google Patents
[go: Go Back, main page]

JPH0691609B2 - Linear array sensor abnormality detection method - Google Patents

Linear array sensor abnormality detection method

Info

Publication number
JPH0691609B2
JPH0691609B2 JP60239328A JP23932885A JPH0691609B2 JP H0691609 B2 JPH0691609 B2 JP H0691609B2 JP 60239328 A JP60239328 A JP 60239328A JP 23932885 A JP23932885 A JP 23932885A JP H0691609 B2 JPH0691609 B2 JP H0691609B2
Authority
JP
Japan
Prior art keywords
image
detector
linear array
array sensor
average value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP60239328A
Other languages
Japanese (ja)
Other versions
JPS62100070A (en
Inventor
弘一 本間
文伸 古村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP60239328A priority Critical patent/JPH0691609B2/en
Publication of JPS62100070A publication Critical patent/JPS62100070A/en
Publication of JPH0691609B2 publication Critical patent/JPH0691609B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Picture Signal Circuits (AREA)
  • Transforming Light Signals Into Electric Signals (AREA)
  • Facsimile Heads (AREA)
  • Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
  • Facsimiles In General (AREA)
  • Facsimile Scanning Arrangements (AREA)
  • Facsimile Image Signal Circuits (AREA)

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は、リニアアレイセンサの異常検出器の自動検出
にかかわり、特に検出器別の画像平均値を求めることに
より精度の高い異常検出を可能とするに好適なリニアア
レイセンサ異常検出方式に関する。
DETAILED DESCRIPTION OF THE INVENTION [Industrial field of application] The present invention relates to automatic detection of anomaly detectors of a linear array sensor, and in particular enables highly accurate anomaly detection by obtaining an image average value for each detector. A linear array sensor abnormality detection method suitable for

〔従来の技術〕[Conventional technology]

従来、検出器の異常によるスパイク状のノイズ検出し、
その検出器に関する画像データを周囲の画像データから
補間する方式については、安岡らの国立公害研究所報告
第39号(1983年)第119頁から第128頁において論じられ
ている。
Conventionally, spike-like noise is detected due to abnormality of the detector,
The method of interpolating the image data of the detector from the surrounding image data is discussed in Yasuoka et al., National Pollution Research Institute Report No. 39 (1983), pages 119 to 128.

すなわち、画像走査線垂直方向の座標をr、近傍区間を
N(r)とし、走査線垂直断面の強度分布として得られ
る原画像を、f(r)とする。近傍区間内で画像を一次
多項式、 f(r)=α+βr+ε(r) (1) で表わし、その係数α,βは最小二乗法により求める。
すなわち、推定式は となる。r=0の画像データが、スパイクノイズかどう
かの判定は、 の95%信頼区間からはずれるか否かの仮説検定すなわち
七分布表との比較により実現される。いったん、r=0
の画像データすなわち走査ラインデータが異常であるこ
とが検出できれば、該走査ラインデータは捨て周囲の画
像データで内挿補間する。
That is, let r be the coordinate in the direction perpendicular to the image scanning line, N (r) be the neighborhood section, and f (r) be the original image obtained as the intensity distribution of the cross section perpendicular to the scanning line. An image is represented by a first-order polynomial, f (r) = α + βr + ε (r) (1) in the neighborhood section, and the coefficients α and β are obtained by the least square method.
That is, the estimation formula is Becomes The determination as to whether the image data of r = 0 is spike noise is It is realized by a hypothesis test of whether or not it deviates from the 95% confidence interval of, that is, by comparison with the seven distribution table. Once, r = 0
If it is possible to detect that the image data, that is, the scanning line data, is abnormal, the scanning line data is discarded and interpolated with the surrounding image data.

〔発明が解決しようとする問題点〕[Problems to be solved by the invention]

上記従来技術は、走査データの異常が画像の途中から生
ずることも考慮しているため、異常の検出処理を画像上
の各点で行うことから、同一の異常走査ラインも背景
(平均値あるいは分散値)の違いにより、ある場所では
異常、別の場所では正常という判定を受けることがあ
り、従って異常走査データの検出精度が低下する恐れが
あるという問題点があった。
Since the above-mentioned conventional technique also considers that an abnormality in the scan data occurs in the middle of the image, the abnormality detection processing is performed at each point on the image. There is a problem in that, depending on the difference in the values), it may be determined that an abnormality is found at one place and normal at another place, and therefore the accuracy of detecting abnormal scan data may be reduced.

本発明の目的は、上記問題点を解決し、精度の高いリニ
アアレイセンサ異常検出方式を提供することにある。
An object of the present invention is to solve the above problems and provide a highly accurate linear array sensor abnormality detection method.

〔問題点を解決するための手段〕[Means for solving problems]

上記目的は、リニアアレイセンサ画像の検出器別の画像
平均値を算出し、該検出器別平均値上でスパイク状ノイ
ズを検出することにより、達成される。
The above object is achieved by calculating an image average value for each detector of the linear array sensor image and detecting spike noise on the average value for each detector.

〔作用〕[Action]

検出器別画像平均値を、検出器座標すなわち走査ライン
垂直方向座標rにたいし、F(r)と表わす。r=0で
検出器の異常が生じており、撮像対象にかかわらず固定
値bが出力されるとすれば、検出器別画像平均値F
(0)もやはりbとなる。ところで、異常検出器座標r
=0の近傍N(r)における画像強度の分散は、原画像
1ライン強度f(r)についての分散より、検出器別平
均値F(r)についての分散の方が小さくなる。原画像
の隣接画素間に相関のない特殊な場合では、平均値を計
算するライン数をMとすると、近傍N(r)における画
像強度分散は1/Mに比例して減少する。原画像の隣接画
素間に相関がある場合にも、ライン数Mが充分大きい場
合にはやはり1/Mで減少する。言いかえれば、原画像1
ライン強度f(r)に比べ、検出器別平均値F(r)は
滑らかに変化する。従って、異常検出器に対する画像平
均値bの異常判定の精度が大幅に向上する。(3)式で
みれば、分母が で増大し、その結果、tが小さくなるから周囲データか
らの推定値と前記bとの差が小さくとも、精度良く異常
判定を行えることになる。
The image average value for each detector is represented by F (r) with respect to the detector coordinate, that is, the coordinate r in the vertical direction of the scanning line. If a detector abnormality occurs at r = 0 and a fixed value b is output regardless of the imaging target, the detector average image value F
(0) is also b. By the way, the abnormality detector coordinate r
The variance of the image intensity in the vicinity N (r) of = 0 is smaller than the variance of the original image one-line intensity f (r) with respect to the average value F (r) for each detector. In a special case where there is no correlation between adjacent pixels of the original image, if the number of lines for calculating the average value is M, the image intensity variance in the neighborhood N (r) decreases in proportion to 1 / M. Even when there is a correlation between adjacent pixels of the original image, when the number of lines M is sufficiently large, it also decreases by 1 / M. In other words, original image 1
The average value F (r) for each detector changes more smoothly than the line intensity f (r). Therefore, the accuracy of the abnormality determination of the image average value b for the abnormality detector is significantly improved. According to equation (3), the denominator is Therefore, even if the difference between the estimated value from the surrounding data and b is small, the abnormality determination can be accurately performed.

〔実施例〕〔Example〕

以下、本発明の一実施例を第1図〜第4図により説明す
る。地球観測衛星1に搭載された、多数の検出器からな
るリニアアレイセンサ15は、衛星の移動17とともに、地
表16をWの幅で1ラインずつレンズ系18を通して観測
し、画像データを地上に送信する。アンテナ2は、受信
した画像データを高速データ記録システム3に送り、記
録システム3は該画像データを高密度磁気テープ4に記
録する。いったん高密度磁気テープ4に記録された画像
データは、画像処理プロセッサ6の能力に合わせ、低速
で磁気テープ装置5から1走査ラインずつ読み込まれ
る。jライン目に読み込まれたr画素目の画像データfj
(r)は、処理プロセッサ6により、累積画像データバ
ッファ7から読み出されたj−1ライン目までの累積画
像Sj-1(r)と次式により加算され、 Sj(r)=Sj-1(r)+fj(r) (4) 再度累積画像データバッファ7に書き込まれる。磁気テ
ープ装置5から読み込まれた画像データは、累積処理と
は別に、ディスク装置8に書き込まれる。所定の画像ラ
インMが読み込まれた後、スパイク状ノイズ検出手段9
は、累積画像データバッファ7から検出器別累積画像デ
ータSM(r)を読み出し、Mで徐し平均値f(r)とし
た後、スパイク状ノイズ検出処理を加える。第3図に本
処理手順のフローチャートを示す。第3図の処理は、各
検出器は番号rに対し実行される。ボックス19では、第
4図に示すような番号rの近傍N(r)の平均画像強度
f(i),i∈N(r){・印で示す}に対する最小自乗
直線近似式の係数 を求める。ここでnは近傍N(r)中のデータ点数であ
る。ボックス20では、該直線近似式からの近次誤差の自
乗和Se(r)を近傍N(r)について求める。ボックス
21では、問題にしている検出器番号rの検出器に対する
画像平均値f(r)とその推定値 との差 を、近傍の近似誤差分散 と比較した値を求める。ボックス22では95%信頼区間で
七分布と比較し、大きければ、ボックス23でrを異常検
出器として、異常検出器番号ファイル10に出力する。異
常データ補間処理プロセッサ11は、異常検出器ファイル
10より異常検出器番号rを読み出し、つぎの補間処理を
行う。
An embodiment of the present invention will be described below with reference to FIGS. The linear array sensor 15 consisting of a large number of detectors mounted on the earth observation satellite 1 observes the ground surface 16 line by line with a width of W through the lens system 18 as the satellite moves 17, and transmits image data to the ground. To do. The antenna 2 sends the received image data to the high-speed data recording system 3, and the recording system 3 records the image data on the high-density magnetic tape 4. The image data once recorded on the high-density magnetic tape 4 is read from the magnetic tape device 5 one scan line at a low speed in accordance with the capability of the image processor 6. Image data fj of the r pixel read in the jth line
(R) is added by the processor 6 with the cumulative image Sj -1 (r) up to the j-1th line read from the cumulative image data buffer 7 according to the following equation, and Sj (r) = Sj -1 (R) + fj (r) (4) It is written in the accumulated image data buffer 7 again. The image data read from the magnetic tape device 5 is written to the disk device 8 separately from the accumulating process. After the predetermined image line M is read, spike noise detecting means 9
Reads the accumulated image data S M (r) for each detector from the accumulated image data buffer 7, divides by M to obtain the average value f (r), and then performs spike noise detection processing. FIG. 3 shows a flowchart of this processing procedure. The process of FIG. 3 is performed for each detector with the number r. In box 19, the coefficient of the least-squares linear approximation formula for the average image intensity f (i), iεN (r) {indicated by a mark} in the neighborhood N (r) with the number r as shown in FIG. Ask for. Here, n is the number of data points in the neighborhood N (r). In box 20, the square sum Se (r) of near-order errors from the linear approximation formula is obtained for the neighborhood N (r). box
In 21, the image average value f (r) and its estimated value for the detector with the detector number r in question Difference from Is the approximation error variance of the neighborhood Calculate the value compared with. In box 22, the 95% confidence interval is compared with seven distributions, and if larger, r is output as an anomaly detector in the anomaly detector number file 10 in a box 23. Anomaly data interpolation processor 11 is an anomaly detector file
The abnormality detector number r is read from 10 and the following interpolation processing is performed.

fj(r)={fj(r−1)+fj(r+1)}/2 (5) 補間処理結果fj(r)は磁気ディスク8の異常データの
上書き込まれる。磁気ディスク8の画像データは、その
後画像補性処理システム12において、検出器の微小感度
むら補正、幾何学的歪補正等、従来の補正処理を受け、
磁気テープ装置13から磁気テープ14に出力される。
fj (r) = {fj (r-1) + fj (r + 1)} / 2 (5) The interpolation process result fj (r) is overwritten with abnormal data on the magnetic disk 8. The image data of the magnetic disk 8 is then subjected to conventional correction processing such as correction of minute sensitivity unevenness of the detector and correction of geometric distortion in the image complementation processing system 12,
The data is output from the magnetic tape device 13 to the magnetic tape 14.

また、本発明は上記説明のリニアアレイセンサ搭載衛星
用地上局の画像処理方式として実施できるばかりでな
く、X線CTやFax装置など、他のリニアアレイ型撮像セ
ンサの検出器異常検出に適用できる。X線CT装置では、
少数の検出器の異常は再構成画像上では画像全体の画質
劣化となって徐々に現われるため、本発明により原セン
サ画像データを処理し、異常検出器データを自動検出
し、補間処理を行なうことによる画質向上の効果は絶大
である。また、FAX装置では、上記実施例のごとく高精
度な異常検出アルゴリズムを用いなくても、検出器間で
平均値の差の絶対値を求め、両隣りとも設定閾値を越え
る差があるかどうかを判定すれば簡単に本発明を実施で
きる。
Further, the present invention can be implemented not only as an image processing method of the satellite ground station equipped with the linear array sensor described above, but also applicable to detector abnormality detection of other linear array type image sensor such as an X-ray CT or fax machine. . With an X-ray CT system,
Since the abnormality of a small number of detectors gradually appears on the reconstructed image as the deterioration of the image quality of the entire image, it is necessary to process the original sensor image data, automatically detect the abnormality detector data, and perform the interpolation process according to the present invention. The effect of improving the image quality is extremely large. Further, in the fax machine, even without using a highly accurate abnormality detection algorithm as in the above embodiment, the absolute value of the difference in average value between the detectors is determined, and whether or not there is a difference exceeding the set threshold value on both sides is checked. If determined, the present invention can be easily implemented.

〔発明の効果〕〔The invention's effect〕

本発明によれば、リニアアレイセンサの異常検出器の自
動検出を高い精度で行うことができ、また1次元データ
の処理による異常判定処理となるため2次元画像全面の
判定処理に比べ高速となる効果がある。
According to the present invention, the automatic detection of the abnormality detector of the linear array sensor can be performed with high accuracy, and since the abnormality determination processing is performed by processing the one-dimensional data, the speed is higher than the determination processing of the entire surface of the two-dimensional image. effective.

【図面の簡単な説明】[Brief description of drawings]

第1図は本発明の一実施例の全体構成図、第2図は衛星
搭載リニアアレイセンサによる地表走査の様子を示す
図、第3図は検出器別画像平均値からの異常検出器番号
の検出処理手順のフローチャート、第4図は近傍データ
の線型近似式に基づく異常データの仮説検出の概念図で
ある。
FIG. 1 is an overall configuration diagram of an embodiment of the present invention, FIG. 2 is a diagram showing a state of ground surface scanning by a satellite-mounted linear array sensor, and FIG. 3 is an abnormality detector number from an image average value for each detector. FIG. 4 is a flowchart of a detection processing procedure, and FIG. 4 is a conceptual diagram of hypothesis detection of abnormal data based on a linear approximation formula of neighborhood data.

Claims (1)

【特許請求の範囲】[Claims] 【請求項1】複数の検出器を備えたリニアアレイセンサ
と、該センサにより撮像された画像の入力装置と、画像
処理装置と、画像出力装置よりなる画像処理システムに
おいて、上記検出器別の画像平均値を算出し、算出され
た平均値にもとづき、該平均値中のスパイク状ノイズを
検出し、検出されたスパイク状ノイズに対応する検出器
出力を補正することを特徴とするリニアアレイセンサ異
常検出方式。
1. An image processing system comprising a linear array sensor having a plurality of detectors, an input device for an image picked up by the sensors, an image processing device, and an image output device. Abnormal linear array sensor characterized by calculating an average value, detecting spike-like noise in the average value based on the calculated average value, and correcting the detector output corresponding to the detected spike-like noise Detection method.
JP60239328A 1985-10-28 1985-10-28 Linear array sensor abnormality detection method Expired - Fee Related JPH0691609B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP60239328A JPH0691609B2 (en) 1985-10-28 1985-10-28 Linear array sensor abnormality detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP60239328A JPH0691609B2 (en) 1985-10-28 1985-10-28 Linear array sensor abnormality detection method

Publications (2)

Publication Number Publication Date
JPS62100070A JPS62100070A (en) 1987-05-09
JPH0691609B2 true JPH0691609B2 (en) 1994-11-14

Family

ID=17043085

Family Applications (1)

Application Number Title Priority Date Filing Date
JP60239328A Expired - Fee Related JPH0691609B2 (en) 1985-10-28 1985-10-28 Linear array sensor abnormality detection method

Country Status (1)

Country Link
JP (1) JPH0691609B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6480195A (en) * 1987-09-22 1989-03-27 Toshiba Corp Film picture reader

Also Published As

Publication number Publication date
JPS62100070A (en) 1987-05-09

Similar Documents

Publication Publication Date Title
EP1328903B1 (en) Method and apparatus for digital image defect correction and noise filtering
Tang et al. Cone beam volume CT image artifacts caused by defective cells in x‐ray flat panel imagers and the artifact removal using a wavelet‐analysis‐based algorithm
EP1291677A2 (en) Method for stitching partial radiation images to reconstruct a full image
US6801670B2 (en) Image reading apparatus and shading correction data acquiring method
JP4377571B2 (en) X-ray flat panel detector, X-ray image diagnostic apparatus, and X-ray image correction method
JPH0642882B2 (en) Desired image signal range determination method
JPH0877340A (en) Offset for scanning system and method for decision of gain correction
CN102970931B (en) For determining the method for the roomage response signature of the detector in computed radiography
WO2020073132A1 (en) Skewed x-ray detection apparatus and method for pipeline use
US6697663B1 (en) Method and apparatus for reducing noise artifacts in a diagnostic image
US5822392A (en) Multi-resolution detection for increasing in an x-ray imaging implementation of an object
JPH08201464A (en) Method for detecting S / N value of television video signal
JPH0691609B2 (en) Linear array sensor abnormality detection method
JP2012049947A5 (en)
CN110400358B (en) Annular artifact removing method for flat panel detector
US11769331B2 (en) Side view camera detecting wheels
US6455832B1 (en) Defect discriminating method and device for solid-state detector
CN111208491B (en) Method for eliminating miscellaneous points of high-resolution solid-state laser radar point cloud
CN110220475B (en) A Linear CCD 2D Variable Speed Imaging Method Based on Image Segmentation
JP2598636B2 (en) Image processing condition determination method
JP2000244825A (en) Method and device for judging defect of solid-state detector
JPS58163339A (en) Subtraction treatment of x-ray image
EP0929194B1 (en) Motion detection
JP3253990B2 (en) Multi-dimensional image measurement device
JP3305749B2 (en) Radiation image reader

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees