KR20100077772A - 위성용 광학 센서의 보정방법 - Google Patents

위성용 광학 센서의 보정방법 Download PDF

Info

Publication number
KR20100077772A
KR20100077772A KR1020080135810A KR20080135810A KR20100077772A KR 20100077772 A KR20100077772 A KR 20100077772A KR 1020080135810 A KR1020080135810 A KR 1020080135810A KR 20080135810 A KR20080135810 A KR 20080135810A KR 20100077772 A KR20100077772 A KR 20100077772A
Authority
KR
South Korea
Prior art keywords
gain
optical sensor
light source
intensity
satellite
Prior art date
Application number
KR1020080135810A
Other languages
English (en)
Inventor
명환춘
윤형식
Original Assignee
한국항공우주연구원
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 한국항공우주연구원 filed Critical 한국항공우주연구원
Priority to KR1020080135810A priority Critical patent/KR20100077772A/ko
Publication of KR20100077772A publication Critical patent/KR20100077772A/ko

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/70—Circuitry for compensating brightness variation in the scene
    • H04N23/76—Circuitry for compensating brightness variation in the scene by influencing the image signals
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/80—Camera processing pipelines; Components thereof
    • H04N23/81—Camera processing pipelines; Components thereof for suppressing or minimising disturbance in the image signal generation

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

위성용 광학 센서의 보정방법을 개시한다. 상기 보정방법은 지구에 대한 영상 정보를 보정하는 방법에 있어서, (a) 상기 광학 센서를 통해 입력된 광원의 세기의 제곱이 잡음과 동일하다고 가정한 경우 상기 광학 센서의 디지털 출력값의 확률변수에서 양자화 잡음을 고려한 제1 확률변수의 평균값을 구하는 단계, (b) 상기 광학 센서의 비선형 이득, 픽셀 이득 및 광학계 이득의 곱(gnso)의 후보값을 설정하는 단계, (c) 상기 gnso의 후보값을 이용하여 상기 광학 센서의 입력 광원의 세기를 구하는 단계, (d) 상기 단계 (c)에서 구해진 입력 광원의 세기를 이용하여 상기 광학 센서를 통해 입력된 광원의 세기가 잡음과 동일하다고 가정한 경우 상기 광학 센서의 디지털 출력값의 확률변수에서 양자화 잡음을 고려한 제2 확률변수의 평균값을 구하는 단계, (e) 상기 제2 확률변수의 평균값을 이용하여 상기 광학 센서의 입력 광원의 세기를 구하는 단계를 포함하되, 상기 단계 (c) 및 단계 (e)의 상기 광학 센서의 입력 광원의 세기가 동일해질 때까지 상기 단계 (b) 내지 단계 (e)를 적어도 한번 반복한 후 이득값(gain)을 얻어 상기 영상 정보를 보정한다.
위성, 광학 센서, 잡음, 광원의 세기, 이득값

Description

위성용 광학 센서의 보정방법{Calibration method of optics-sensor for satellite}
본 발명은 광학 센서의 보정방법에 관한 것으로, 특히 위성용 광학 센서의 보정방법에 관한 것이다.
위성용 광학 센서는 위성(satellite)에 장착되어 지구의 영상을 관측하기 위한 광학 센서를 말한다. 이러한 위성용 광학 센서는 위성 운용 중 이득값(gain)의 변화(degradation)를 추정하여 지상에서 측정된 영상 정보로부터 원래의 광량을 도출함으로써 영상 정보를 보정해야 한다.
위성용 광학 센서의 보정방법은 크게 세 가지로 분류된다.
첫째는 절대 광원을 이용하여 위성용 광학 센서를 보정하는 것이다. 절대 광원으로는 크게 두 가지 방법이 이용된다. 하나는 위성 내에 절대 광원을 발생시키는 장치를 장착하는 것이고 다른 하나는 태양을 광원으로 이용하는 것이다. 절대 광원을 발생시키는 장치를 장착하는 것은 측정의 용이함으로 인하여 위성용 광학 센서의 운용에 있어 제약이 적다는 장점이 있지만, 부수적인 장치가 요구된다는 점에서 적용의 어려움이 있다. 하지만, 현재 많은 기상 관측용 망원경의 경우에 적외 선 센서의 보정을 위하여 위성망원경에 함께 장착된 흑체(black body) 망원경을 사용하고 있다. 위성용 광학 센서는 일반적으로 태양을 광원으로 이용하여 센서를 보정하고 있는 실정이다. 이 경우, 부수적인 장치 없이 안정적인 광원을 얻을 수 있다는 장점이 있지만, 태양을 관측하기 위하여 위성용 광학 센서의 운용에 부수적인 제약이 요구되며 태양을 광원으로 받아들이기 위하여 요구되어지는 확산기 (diffuser)의 성능 감환 현상이 단점이다. 이 방법은 MOS, SeaWiFS, MERIS 등에 적용되어 이용되고 있다.
둘째는 절대 광원과는 달리 센서의 상대적인 변화량을 관측하는 방법이다. 지표면의 일정한 반사광, 달의 반사광, 위성용 광학 센서에 장착된 LED와 램프(lamp) 등을 이용하여 센서의 이득 변화량을 감지하는 것이다. 센서의 상대적인 변화량을 관측하는 방법은 절대 광원을 사용하는 것에 비하여 운용과 장비에 있어서 비교적 수월한 장점이 있지만, 센서의 검출기 전체에 대하여 동일한(homogeneous) 입력을 가하기 어렵다는 문제점이 있다.
세째는 여러 위성들 간의 측정 데이터를 공유함으로써, 다른 위치와 시간에 측정된 동일한 데이터의 오차를 이용하여 센서를 보정하는 방법(intercalibration)이다. 하지만, 이는 많은 위성용 광학 센서 시스템간의 협력이 요구된다는 문제점이 있다.
본 발명이 이루고자 하는 기술적 과제는 상기 문제점들을 해결하기 위한 것으로 광원의 동일성(homogeneity)에 관한 제약성을 극복하고 위성용 광학 센서에 용이하게 적용할 수 있는 보정방법을 제공하는 데 있다.
본 발명의 일 양태에 따르면, 위성용 광학 센서에서 지구에 대한 영상 정보를 보정하는 방법에 있어서, (a) 상기 광학 센서를 통해 입력된 광원의 세기의 제곱이 잡음과 동일하다고 가정한 경우 상기 광학 센서의 디지털 출력값의 확률변수에서 양자화 잡음을 고려한 제1 확률변수의 평균값을 구하는 단계, (b) 상기 광학 센서의 비선형 이득, 픽셀 이득 및 광학계 이득의 곱(gnso)의 후보값을 설정하는 단계, (c) 상기 gnso의 후보값을 이용하여 상기 광학 센서의 입력 광원의 세기를 구하는 단계, (d) 상기 단계 (c)에서 구해진 입력 광원의 세기를 이용하여 상기 광학 센서를 통해 입력된 광원의 세기가 잡음과 동일하다고 가정한 경우 상기 광학 센서의 디지털 출력값의 확률변수에서 양자화 잡음을 고려한 제2 확률변수의 평균값을 구하는 단계, (e) 상기 제2 확률변수의 평균값을 이용하여 상기 광학 센서의 입력 광원의 세기를 구하는 단계를 포함하되, 상기 단계 (c) 및 단계 (e)의 상기 광학 센서의 입력 광원의 세기가 동일해질 때까지 상기 단계 (b) 내지 단계 (e)를 적어도 한번 반복한 후 이득값(gain)을 얻어 상기 영상 정보를 보정하는 위성용 광학 센서의 보정방법을 제공한다.
본 발명의 다른 양태에 따르면, 위성용 광학 센서의 보정방법을 수행하기 위해 디지털 처리 장치에 의해 실행될 수 있는 명령어들의 프로그램이 구현되어 있으며, 상기 디지털 처리 장치에 의해 판독될 수 있는 프로그램이 기록된 기록매체를 제공한다.
본 발명의 또 다른 양태에 따르면, 지구의 영상을 관측하는 위성용 광학 센서에 있어서, 상기 지구의 영상에 대한 광원을 입력받아 빛의 투과율이 서로 상이한 적어도 하나의 필터를 통하여 상기 광원을 통과시키는 대역 통과 필터부, 상기 광원을 전자의 형태인 영상 정보로 변환하는 검출부, 및 상기 대역 통과 필터부에 의해 임의로 만들어지는 인공 잡음인 △L이 △L=L2 및 △L=L(여기서, L은 광원의 세기임)인 두 가지 가정에 대하여 변형된 카이스퀘어 분포를 얻어, 상기 얻어진 변형된 카이스퀘어 분포에서 상기 두 가지 가정에 대하여 동일한 상기 광원의 세기가 구해질 때까지 반복함으로써 상기 위성용 광학 센서의 이득값을 얻어서 상기 지구의 영상에 대한 영상 정보를 보정하는 영상 처리부를 포함하는 위성용 광학 센서를 제공한다.
상기 대역 통과 필터부는 빛의 투과율은 서로 상이하고 동일한 주파수 대역의 빛을 통과시키는 적어도 하나의 필터를 포함할 수 있다.
상기 영상 처리부는 지상에서 상기 지구의 영상에 대한 영상 정보를 보정할 수 있다.
위성용 광학 센서에서 임의의 지구 영상만으로 광학 센서의 보정을 가능하게 하여 위성용 광학 센서의 운용의 자율성을 높이고 태양 빛에 대한 노출 위험도를 최소화한다. 입력 광원과 운용 방식에 대한 별도의 제약 조건 없이 위성용 광학 센서의 센서 이득값들을 필요할 때마다 용이하게 보정할 수 있다.
이하, 첨부된 도면을 참조하여 본 발명의 바람직한 실시예를 상세히 설명한다. 다음 실시예들은 본 발명의 이해를 돕기 위한 것이며 본 발명이 하기의 실시예에 한정되는 것은 아니다.
본 실시예에서 사용되는 '~부'라는 용어는 소프트웨어 또는 FPGA(field-programmable gate array) 또는 ASIC과 같은 하드웨어 구성요소를 의미하며, '~부'는 어떤 역할들을 수행한다. 그렇지만 '~부'는 소프트웨어 또는 하드웨어에 한정되는 의미는 아니다. '~부'는 어드레싱할 수 있는 저장 매체에 있도록 구성될 수도 있고 하나 또는 그 이상의 프로세서들을 재생시키도록 구성될 수도 있다. 따라서, 일 예로서 '~부'는 소프트웨어 구성요소들, 객체지향 소프트웨어 구성요소들, 클래스 구성요소들 및 태스크 구성요소들과 같은 구성요소들과, 프로세스들, 함수들, 속성들, 프로시저들, 서브루틴들, 프로그램 코드의 세그먼트들, 드라이버들, 펌웨어, 마이크로코드, 회로, 데이터,데이터베이스, 데이터 구조들, 테이블들, 어레이들, 및 변수들을 포함한다. 구성요소들과 '~부'들 안에서 제공되는 기능은 더 작은 수의 구성요소들 및 '~부'들로 결합되거나 추가적인 구성요소들과 '~부'들로 더 분 리될 수 있다. 뿐만 아니라, 구성요소들 및 '~부'들은 디바이스 또는 보안 멀티미디어카드 내의 하나 또는 그 이상의 CPU들을 재생시키도록 구현될 수도 있다.
도 1은 본 실시예에 따른 위성용 광학 센서를 나타내는 블록도이다.
도 1을 참조하면, 위성용 광학 센서(100)는 대역 통과 필터부(110), 광학부(120), 관측용 필터부(130), 검출부(140) 및 영상 처리부(150)를 포함한다.
대역 통과 필터부(110)는 빛의 투과율은 서로 상이하고 동일한 주파수 대역의 빛을 통과시키는 임의의 n개의 필터로 구성되며, 위성용 광학 센서(100)의 입력 광원에 대하여 보정을 위한 변형을 가한다. 대역 통과 필터부(110) 내 각각의 필터들의 출력값은 서로 상이한 투과율로 인하여, 원래 광원의 세기(L with upper bar)와는 다르게 　필터 투과율의 분포에 의지하여 광원(L)의 새로운 분포를 형성한다. 이러한 광원(L)의 새로운 분포는 최종적인 영상의 출력값들이 계산된 비-중심 카이스퀘어 분포(non-central Chi square distribution)를 형성하도록 만들며, 추출된 확률분포는 평균과 분산값들을 통하여 위성용 광학 센서(100)의 이득값들과 관련된 유용한 정보를 제공한다.
광학부(120)는 광학 반사 미러(Mirror) 또는 굴절 미러로 구성되며, 대역 통과 필터부(110)를 통과한 변형이 가해진 입력 광원을 입력받아, 광원의 형태인 광정보(photon)를 집적하여 증폭시킨다.
관측용 필터부(130)는 증폭된 광정보(photon)를 입력받아, 유용한 주파수 대역의 광정보들만을 선별한다.
검출부(140)는 CCD(charge-coupled device) 또는 CMOS로 구성되며, 유용한 주파수 대역의 광정보(photon)들을 전자(electron)로 변환한다.
영상 처리부(150)는 변환된 전자의 형태인 영상 정보를 입력받아 영상 정보에 대하여 보정을 수행하여 원하는 영상 정보를 획득하게 한다. 이러한 영상 처리부(150)는 위성용 광학 센서(100)와 일체이거나 지상에 있을 수 있다. 위성용 광학 센서(100)가 지상에 있는 경우는 위성에 장착된 위성용 광학 센서(100)의 검출부(140)로부터 입력 광원이 변환된 전자의 형태인 영상 정보를 입력받아 지상에서 보정을 수행하여 원하는 영상 정보를 획득케 할 수 있다.
본 실시예에서는 임의의 광원의 세기(L)에 대하여 센서 모델의 이득값들을 발견하여 위성용 광학 센서(100)에 의해 획득된 영상 정보의 보정에 활용한다.
도 2는 본 실시예에 따른 위성용 광학 센서의 보정방법을 설명하는 도면이다.
본 실시예에서는 비선형 센서 모델을 가정한다. 비선형 센서 모델에 따른 디지털 출력값(D)는 다음 수학식 1과 같다.
Figure 112008090001338-PAT00001
여기서, L(input radiance)은 입력 광원의 세기(빛의 세기)이며 단위는 W/m2/μm/sr(steradian), T는 노출 시간, gl는 선형 이득, gn은 비선형 이득, gs는 픽셀 이득 및 go는 광학계 이득을 나타낸다.
센서 모델에 있어 고려한 잡음의 종류는 암전류 옵셋(OD), 전기 옵셋(OE), shot 잡음(Ns), read-out 잡음(NR), 양자화 잡음(NQ) 등이다.
NDSR 및 NEQ는 각각 다음 수학식 2 및 수학식 3과 같다.
Figure 112008090001338-PAT00002
Figure 112008090001338-PAT00003
대역 통과 필터부(110)를 이용하여 입력 광원의 세기에 변형을 가하여 다음 수학식 4와 같이 입력 광원의 세기(L)에 잡음을 가정한다.
Figure 112008090001338-PAT00004
여기서,
Figure 112008090001338-PAT00005
는 입력 광원의 잡음을 가정했을 때의 입력 광원의 세기의 평균값,
Figure 112008090001338-PAT00006
는 대역 통과 필터부(110)에 의해 임의로 만들어지는 인공 잡음을 나타내며, 0을 평균으로 하는 가우시안(Gaussian) 확률변수로 가정된다.(도 2의 210 참조)
이에 따라 수학식 1의 디지털 출력값(D)는 다음 수학식 5와 같이 변형된 디지털 출력값(
Figure 112008090001338-PAT00007
)이 된다.
Figure 112008090001338-PAT00008
여기서, DN은 디지털 출력값(D)의 확률변수를 말하며, 다음 수학식 6과 같다.
Figure 112008090001338-PAT00009
여기서,
Figure 112008090001338-PAT00010
,
Figure 112008090001338-PAT00011
,
Figure 112008090001338-PAT00012
Figure 112008090001338-PAT00013
및
Figure 112008090001338-PAT00014
이다.
확률변수 DN의 확률밀도함수(probability density function: PDF)는 다음 수학식 7과 같다.
Figure 112008090001338-PAT00015
Figure 112008090001338-PAT00016
,
Figure 112008090001338-PAT00017
,
Figure 112008090001338-PAT00018
및
Figure 112008090001338-PAT00019
이다.
위성용 광학 센서의 디지털 출력값(D)의 신호 대 잡음비(signal-to-noise ratio: SNR)를 개선하기 위하여 평균화(averaging) 과정 또는 누적화(summing) 과정을 거친다.
디지털 출력값(D)의 확률변수(DN)에서 양자화 잡음(NQ)을 제외한 확률변수(XLI)는 다음 수학식 8과 같다.
Figure 112008090001338-PAT00020
여기서, k는 누적 개수, n은 평균화 개수 및
Figure 112008090001338-PAT00021
는 다음 수학식 9와 같다.
Figure 112008090001338-PAT00022
디지털 출력값(D)의 확률변수(DN)에서 양자화 잡음(NQ)을 제외한 확률변수(XLI)의 확률밀도함수는 다음 수학식 10과 같다.
Figure 112008090001338-PAT00023
,
Figure 112008090001338-PAT00024
,
Figure 112008090001338-PAT00025
및
Figure 112008090001338-PAT00026
이다.
디지털 출력값(D)의 확률변수(DN)에서 양자화 잡음(NQ)을 고려한 확률변수(XLIQ)의 확률밀도함수는 다음 수학식 11 및 수학식 12와 같다(도 2의 220 참조).
Figure 112008090001338-PAT00027
여기서, 상기 수학식 11에서 n=k=1이다.
Figure 112008090001338-PAT00028
여기서, 상기 수학식 12에서 k〉1이다.
디지털 출력값(D)의 확률변수(DN)에서 양자화 잡음(NQ)을 고려한 확률변수(XLIQ)의 평균(
Figure 112008090001338-PAT00029
) 및 분산(V[XLIQ])은 각각 수학식 13 및 수학식 14와 같다.
Figure 112008090001338-PAT00030
Figure 112008090001338-PAT00031
대역 통과 필터부(110)를 이용하여 γ=L2 및 γ=L인 두 가지 경우에 대하여 변형된 카이스퀘어 분포를 얻어, 얻어진 결과를 이용하여 반복하여 위성용 광학 센서의 이득값(gain)들을 획득한다.
도 3은 본 실시예에 따른 위성용 광학 센서의 보정방법을 나타내는 순서도이다.
도 3을 참조하면, 위성에 장착된 위성용 광학 센서(100)의 검출부(140)로 부터 입력 광원이 변환된 전자의 형태의 영상 정보를 이용하여 △L=L2 인 경우, XLIQ,1의 평균값을 계산한다(S110).
그런 다음, gnso의 후보값을 설정한다(S120). 예를 들어, 위성용 광학 센서의 보정은 궤도에서 주기적으로 수행되어지며 매 보정시 gnso의 후보값은 이전의 보정 과정에서 얻어진 gnso값을 gnso의 후보값으로 설정할 수 있다.
그런 다음, 설정된 gnso의 후보값을 이용하여 다음 수학식 15를 통해 입력 광원의 세기(L1)을 계산한다(S130).
Figure 112008090001338-PAT00032
여기서,
Figure 112008090001338-PAT00033
이다.
노출 시간(T)을 제외한 나머지 부분들은 보정하지 않은 정상적인 영상 측정 모드에서 서로 다른 노출 시간(T1, T2)을 이용하여 다음 수학식 16 및 수학식 17을 구한 후, 두 식의 연립방정식으로 풀어 계산한다.
Figure 112008090001338-PAT00034
Figure 112008090001338-PAT00035
그런 다음, 입력 광원의 세기(L1)를 이용하여 △L=L인 경우 XLIQ,2의 평균값을 계산한다(S140).
그런 다음, S140에서 계산한 XLIQ,2의 평균값을 이용하여 다음 수학식 18을 이용하여 입력 광원의 세기(L2)를 계산한다(S150).
Figure 112008090001338-PAT00036
그런 다음, 입력 광원의 세기(L1)와 입력 광원의 세기(L2)가 동일한지 판단한다(S160). 만약, 입력 광원의 세기(L1)와 입력 광원의 세기(L2)가 동일하다면 다음 수학식 19를 이용하여 선형 이득(glso)를 계산한다(S170). 여기서, 입력 광원의 세기(L1)와 입력 광원의 세기(L2)가 동일한 경우 입력 광원의 세기를 L이라 한다.
Figure 112008090001338-PAT00037
그런 다음, 입력 광원의 세기(L)가 0인 암전류 측정 모드를 통해 암전류 보정값(F'')을 다음 수학식 20을 통해 계산한다.
Figure 112008090001338-PAT00038
만약, S160에서 입력 광원의 세기(L1)와 입력 광원의 세기(L2)가 동일하지 않다면 다시 S120으로 돌아가 gnso의 후보값을 다시 설정한다. 그런 다음 다시 S130 내지 S150을 거친 후, 입력 광원의 세기(L1)와 입력 광원의 세기(L2)가 동일한지 판단한다(S160).
입력 광원의 세기(L1)와 입력 광원의 세기(L2)가 동일해질 때까지 S120 내지 S150을 반복한 후, S170을 거치고 상기 수학식 20을 통해 암전류 보정값(F'')을 계산한다(S180).
상술한 모든 기능은 상기 기능을 수행하도록 코딩된 소프트웨어나 프로그램 코드 등에 따른 마이크로프로세서, 제어기, 마이크로제어기, ASIC(Application Specific Integrated Circuit) 등과 같은 프로세서에 의해 수행될 수 있다. 상기 코드의 설계, 개발 및 구현은 본 발명의 설명에 기초하여 당업자에게 자명하다고 할 것이다.
이상 본 발명에 대하여 실시예를 참조하여 설명하였지만, 해당 기술 분야의 통상의 지식을 가진 자는 본 발명의 기술적 사상 및 영역으로부터 벗어나지 않는 범위 내에서 본 발명을 다양하게 수정 및 변경시켜 실시할 수 있음을 이해할 수 있을 것이다. 따라서 상술한 실시예에 한정되지 않고, 본 발명은 이하의 특허청구범위의 범위 내의 모든 실시예들을 포함한다고 할 것이다.
도 1은 본 실시예에 따른 위성용 광학 센서를 나타내는 블록도이다.
도 2는 본 실시예에 따른 위성용 광학 센서의 보정방법을 설명하는 도면이다.
도 3은 본 실시예에 따른 위성용 광학 센서의 보정방법을 나타내는 순서도이다.

Claims (11)

  1. 위성용 광학 센서에서 지구에 대한 영상 정보를 보정하는 방법에 있어서,
    (a) 상기 광학 센서를 통해 입력된 광원의 세기의 제곱이 잡음과 동일하다고 가정한 경우 상기 광학 센서의 디지털 출력값의 확률변수에서 양자화 잡음을 고려한 제1 확률변수의 평균값을 구하는 단계;
    (b) 상기 광학 센서의 비선형 이득, 픽셀 이득 및 광학계 이득의 곱(gnso)의 후보값을 설정하는 단계;
    (c) 상기 gnso의 후보값을 이용하여 상기 광학 센서의 입력 광원의 세기를 구하는 단계;
    (d) 상기 단계 (c)에서 구해진 입력 광원의 세기를 이용하여 상기 광학 센서를 통해 입력된 광원의 세기가 잡음과 동일하다고 가정한 경우 상기 광학 센서의 디지털 출력값의 확률변수에서 양자화 잡음을 고려한 제2 확률변수의 평균값을 구하는 단계;
    (e) 상기 제2 확률변수의 평균값을 이용하여 상기 광학 센서의 입력 광원의 세기를 구하는 단계를 포함하되,
    상기 단계 (c) 및 단계 (e)의 상기 광학 센서의 입력 광원의 세기가 동일해질 때까지 상기 단계 (b) 내지 단계 (e)를 적어도 한번 반복한 후 이득값(gain)을 얻어 상기 영상 정보를 보정하는 위성용 광학 센서의 보정방법.
  2. 제 1 항에 있어서,
    상기 단계 (a)에서 상기 제1 확률변수의 평균값은 하기 식
    Figure 112008090001338-PAT00039
    를 이용하여 구하는 것을 특징으로 하는 위성용 광학 센서의 보정방법.
    여기서,
    Figure 112008090001338-PAT00040
    ,
    Figure 112008090001338-PAT00041
    , gn은 비선형 이득, gs는 픽셀 이득 및 go는 광학계 이득이다.
  3. 제 1 항에 있어서,
    상기 단계 (c)에서 상기 광학 센서의 입력 광원의 세기는 하기 식
    Figure 112008090001338-PAT00042
    를 이용하여 구하는 것을 특징으로 하는 위성용 광학 센서의 보정방법.
    ,
    Figure 112008090001338-PAT00043
    ,
    Figure 112008090001338-PAT00044
    ,
    Figure 112008090001338-PAT00045
    , gn은 비선형 이득, gs는 픽셀 이득 및 go는 광학계 이득이다.
  4. 제 3 항에 있어서,
    상기 식 중 상기 노출 시간(T)을 제외한 나머지 구성부분들은 보정하지 않은 정상적인 영상 측정 모드에서 서로 다른 노출 시간(T1, T2)을 이용하여 하기 식 (E-1) 및 하기 식 (E-2)를 구한 후, 두 식의 연립방정식으로 풀어 구하는 것을 특징으로 하는 위성용 광학 센서의 보정방법.
    Figure 112008090001338-PAT00046
    (E-1)
    Figure 112008090001338-PAT00047
    (E-2)
    여기서, OD는 암전류 옵셋,
    Figure 112008090001338-PAT00048
    ,
    Figure 112008090001338-PAT00049
    , gn은 비선형 이득, gs는 픽셀 이득 및 go는 광학계 이득이다.
  5. 제 1 항에 있어서,
    상기 단계 (e)에서 상기 광학 센서의 입력 광원의 세기는 하기 식
    Figure 112008090001338-PAT00050
    를 이용하여 구하는 것을 특징으로 하는 위성용 광학 센서의 보정방법.
    여기서,
    Figure 112008090001338-PAT00051
    ,
    Figure 112008090001338-PAT00052
    ,
    Figure 112008090001338-PAT00053
    , gn은 비선형 이득, gs는 픽셀 이득 및 go는 광학계 이득이다.
  6. 제 1 항에 있어서,
    상기 이득값(gain)은 선형 이득, 픽셀 이득 및 광학계 이득의 곱인 선형 이득 glso 인 것을 특징으로 하는 위성용 광학 센서의 보정방법.
    여기서, gl는 선형 이득, gs는 픽셀 이득 및 go는 광학계 이득이다.
  7. 제 6 항에 있어서,
    상기 선형 이득 glso 은 하기 식
    Figure 112008090001338-PAT00054
    를 이용하여 구하는 것을 특징으로 하는 위성용 광학 센서의 보정방법.
    여기서,
    Figure 112008090001338-PAT00055
    ,
    Figure 112008090001338-PAT00056
    , gn은 비선형 이득, gs는 픽셀 이득 , go는 광학계 이득,
    Figure 112008090001338-PAT00057
    , gl는 선형 이득 및 OE 는 전기 옵셋이다.
  8. 제1항 내지 제7항 중 어느 한 항에 기재된 위성용 광학 센서의 보정방법을 수행하기 위해 디지털 처리 장치에 의해 실행될 수 있는 명령어들의 프로그램이 구현되어 있으며, 상기 디지털 처리 장치에 의해 판독될 수 있는 프로그램이 기록된 기록매체.
  9. 지구의 영상을 관측하는 위성용 광학 센서에 있어서,
    상기 지구의 영상에 대한 광원을 입력받아 빛의 투과율이 서로 상이한 적어도 하나의 필터를 통하여 상기 광원을 통과시키는 대역 통과 필터부;
    상기 광원을 전자의 형태인 영상 정보로 변환하는 검출부; 및
    상기 대역 통과 필터부에 의해 만들어지는 잡음인 △L이 △L=L2 및 △L=L(여기서, L은 광원의 세기임)인 두 가지 가정에 대하여 변형된 카이스퀘어 분포를 얻어, 상기 얻어진 변형된 카이스퀘어 분포에서 상기 두 가지 가정에 대하여 동일한 상기 광원의 세기가 구해질 때까지 반복함으로써 상기 위성용 광학 센서의 이득값을 얻어서 상기 지구의 영상에 대한 영상 정보를 보정하는 영상 처리부를 포함하는 위성용 광학 센서.
  10. 제 9 항에 있어서,
    상기 대역 통과 필터부는 빛의 투과율은 서로 상이하고 동일한 주파수 대역의 빛을 통과시키는 적어도 하나의 필터를 포함하는 것을 특징으로 하는 위성용 광학 센서.
  11. 제 9 항에 있어서,
    상기 영상 처리부는 지상에서 상기 지구의 영상에 대한 영상 정보를 보정하는 것을 특징으로 하는 위성용 광학 센서.
KR1020080135810A 2008-12-29 2008-12-29 위성용 광학 센서의 보정방법 KR20100077772A (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020080135810A KR20100077772A (ko) 2008-12-29 2008-12-29 위성용 광학 센서의 보정방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020080135810A KR20100077772A (ko) 2008-12-29 2008-12-29 위성용 광학 센서의 보정방법

Publications (1)

Publication Number Publication Date
KR20100077772A true KR20100077772A (ko) 2010-07-08

Family

ID=42639082

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020080135810A KR20100077772A (ko) 2008-12-29 2008-12-29 위성용 광학 센서의 보정방법

Country Status (1)

Country Link
KR (1) KR20100077772A (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116167445A (zh) * 2023-02-20 2023-05-26 北京百度网讯科技有限公司 量子测量模式的处理方法、装置及电子设备

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116167445A (zh) * 2023-02-20 2023-05-26 北京百度网讯科技有限公司 量子测量模式的处理方法、装置及电子设备
CN116167445B (zh) * 2023-02-20 2023-10-03 北京百度网讯科技有限公司 量子测量模式的处理方法、装置及电子设备

Similar Documents

Publication Publication Date Title
US7294817B2 (en) System and methods for determining nonuniformity correction parameters in detector-array imaging
CN105841815B (zh) 一种成像光谱仪ccd光谱图像暗电流校正方法
US9876968B2 (en) Drift correction method for infrared imaging device
US9058648B2 (en) Image acquisition for chemiluminescent samples
US7230741B2 (en) Optimum non-uniformity correction for imaging sensors
CN101889869B (zh) 摄像设备及其控制方法
Kobelski et al. Calibrating data from the Hinode/X-Ray Telescope and associated uncertainties
US10070076B2 (en) Drift correction method for infrared imaging device
US6219405B1 (en) Method and apparatus for sensing an image
KR20100077772A (ko) 위성용 광학 센서의 보정방법
Miville-Deschênes et al. Optimizing ISOCAM data processing using spatial redundancy
CN110807812B (zh) 一种基于先验噪声模型的数字图像传感器系统误差标定方法
CN114240801A (zh) 一种遥感图像非均匀校正方法
KR101637552B1 (ko) 렌즈에 의한 불균일 영상 보정 장치 및 그 방법
Christen et al. Fast conversion factor (gain) measurement of a CCD using images with vertical gradient
KR102286373B1 (ko) 디지털 이미저의 교정 방법
US9854138B2 (en) Fixed pattern noise reduction
US11736828B2 (en) Simultaneous and consistent handling of image data and associated noise model in image processing and image synthesis
RU2639680C2 (ru) Способ и система определения в реальном времени сигналов, подлежащих суммированию, среди множества принятых сигналов
Kerr The ISO definition of the dynamic range of a digital still camera
KR20090058299A (ko) 확률적 노출시간 조정에 의한 위성 광학 망원경의 센서보정방법
KR102013206B1 (ko) 보정 데이터의 생성 방법 및 이를 이용한 적외선 영상 처리 방법
CN107249429A (zh) 增益调整装置、增益调整程序、内窥镜以及内窥镜装置
US20060186338A1 (en) Method for improving measurement accuracy of infrared imaging radiometers
Mansouri et al. Systematic noise characterization of a CCD camera: application to a multispectral imaging system

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E601 Decision to refuse application