一種星載面陣ccd相機點擴散函數(shù)測量裝置及方法
【專利摘要】一種星載面陣CCD相機點擴散函數(shù)測量裝置及方法，裝置包括光源系統(tǒng)、靶標、平行光管、面陣CCD相機和圖像采集系統(tǒng)；靶標放置在平行光管的焦面上，平行光管與面陣CCD相機的光軸共軸放置，靶標被光源系統(tǒng)產(chǎn)生的均勻光照亮，經(jīng)過平行光管和面陣相機，成像在面陣CCD相機的焦面上，形成靶標的圖像，圖像采集系統(tǒng)對成像結(jié)果進行采集和處理。本發(fā)明在靶標設計、實驗方法和數(shù)據(jù)處理中考慮了相位對準，能夠保證PSF測量精度，填補了星載面陣CCD相機的PSF測試方法空白。
【專利說明】-種星載面陣CCD相機點擴散函數(shù)測量裝置及方法

【技術(shù)領域】
[0001] 本發(fā)明涉及點擴散函數(shù)的測量裝置及方法，特別是一種星載面陣CCD相機點擴散 函數(shù)的測量裝置及方法，適用于面陣凝視型CCD相機或CMOS相機的點擴散函數(shù)測試，屬于 航天光學遙感【技術(shù)領域】。

【背景技術(shù)】
[0002] 光學系統(tǒng)在理想狀態(tài)下，物空間一點發(fā)出的光能量在像空間也集中在一點上，但 是實際的光學系統(tǒng)成像時，由于衍射和像差以及其它因素的影響，物空間一點發(fā)出的光在 像空間是分布在一定的區(qū)域內(nèi)，其分布曲線稱為點擴散函數(shù)PSF。調(diào)制傳遞函數(shù)MTF是客 觀評價光學系統(tǒng)成像質(zhì)量的重要指標，目前國內(nèi)為了定量化評價相機成像質(zhì)量，發(fā)射之前 都進行實驗室的MTF測試，但MTF只是PSF在頻域的幅值信息，不包含相位信息。目前常用 的MTF測試方法包括以下三種，第一種是對比度矩形靶標法，可以用于面陣或推掃相機MTF 測試，黃巧林等人2006年在《航天返回與遙感》上發(fā)表的《航天光學遙感器MTF測試技術(shù)研 究》一文中給出兩種測試方法，一種是高對比度矩形靶標測試方法，該方法受采樣相位的影 響，測試誤差較大，另一種低頻靶標法對相位不敏感，但易受噪聲影響；第二種是斜刃邊靶 標法，適用于面陣CCD的MTF地面測試，張孝弘等人2006年在中國空間科學學會空間探測 專業(yè)委員會第十九次學術(shù)會議上發(fā)表的《面陣CCD相機的MTF測試技術(shù)》一文介紹了實驗 方法和數(shù)據(jù)處理算法，該處理方法增加了采樣點，但無法準確定位初始相位；第三種是點源 或狹縫測試法，適用于面陣CCD的MTF測試，吳海平碩士學位論文《點源法MTF測試技術(shù)研 究》中詳細介紹了點源法實驗系統(tǒng)，實驗方法是對點光源成像，該方法計算精度較高，但對 準直顯微物鏡的對準精度要求較高，必須保證點光源在一個像素內(nèi)成像。


【發(fā)明內(nèi)容】

[0003] 本發(fā)明解決的技術(shù)問題是：針對國內(nèi)現(xiàn)有星載面陣CCD相機地面只測MTF的現(xiàn)狀， 提出一種星載面陣CCD相機點擴散函數(shù)的實驗室測量裝置及方法，本發(fā)明通過靶標設計和 數(shù)據(jù)處理方法解決相位對準問題，保證面陣CCD相機PSF的測試準確性，最大程度上滿足了 面陣C⑶相機PSF測試的需求。
[0004] 本發(fā)明的技術(shù)方案是：
[0005] -種星載面陣CCD相機點擴散函數(shù)的測量裝置，其特征在于：包括光源系統(tǒng)、靶 標、平行光管、面陣CCD相機和圖像采集系統(tǒng)；
[0006] 靶標放置在平行光管的焦面上，平行光管與面陣CCD相機的光軸共軸放置；
[0007] 光源系統(tǒng)產(chǎn)生均勻光后照射在靶標上，然后經(jīng)過平行光管和面陣C⑶相機，成像 在面陣CCD相機的焦面上，形成靶標的圖像，圖像采集系統(tǒng)對面陣CCD相機的成像結(jié)果進行 采集和處理；
[0008] 靶標包括一百個星點、第一定標線和第二定標線，所述一百個星點分別為第一星 點至第一百星點；
[0009] 第一星點至第一百星點形狀和尺寸均相同，且均為正方形透光孔，正方形透光孔 的邊長為標準長度L的0. 9倍；所述L = (?·Μ?/〇 d，其中fMl為平行光管的焦距，feam為面 陣CCD相機的焦距，d為面陣CCD相機中探測器單元間距；第一定標線和第二定標線均為長 方形透光孔，長方形透光孔的寬度為標準長度L的0. 9倍，長度大于標準長度L的5倍；
[0010] 第一星點至第一百星點、第一定標線和第二定標線在靶標平面直角坐標系XOY中 的坐標分別為：
[0011] 第一定標線長方形結(jié)構(gòu)左上角的頂點坐標為（-20 · L，80 · L);第二定標線長方形 結(jié)構(gòu)左上角的頂點坐標為（80 · L, 200 · L);第一星點的坐標為（0,0);第i星點的坐標為 (a · 19. 9 · L，b · 19. 9 · L)，其中a為i除以10的余數(shù)部分減1，b為i除以10的整數(shù)部 分，l〈i〈100 ;第一百星點的坐標為（179. I *L, 179. 1噸）點；所述星點坐標為星點正方形結(jié) 構(gòu)左上角的坐標；
[0012] 所述靶標平面直角坐標系XOY的坐標原點0為第一星點正方形結(jié)構(gòu)位于左上角的 頂點，X軸的正方向為第一星點指向第二星點，Y軸的正方向為第一星點指向第十一星點；
[0013] 所述的靶標大小長度和寬度均大于300 · L。
[0014] 所述的靶標中的第一星點至第一百星點、第一定標線和第二定標線透光，其余部 分不透光，透光部分與不透光部分對比度不低于100:1。
[0015] 一種利用權(quán)利要求1所述測量裝置的面陣CCD相機點擴散函數(shù)的測量方法，步驟 如下：
[0016] (1)圖像采集系統(tǒng)采集靶標圖像的同時不斷調(diào)整靶標的位置，直至靶標的第一定 標線位于圖像矩陣同一列上，同時第二定標線位于圖像矩陣同一行上，固定此時靶標的位 置；
[0017] (2)設置面陣C⑶相機的曝光時間，使得第一星點至第一百星點在圖像中灰度值 的范圍為：[0. 5Κ，0. 6Κ];利用圖像采集系統(tǒng)采集得到P幀當前曝光時間下的靶標的圖像， 所述K為飽和灰度值，P為大于100的正整數(shù)；
[0018] (3)逐點對第X幀靶標圖像上的星點i的相鄰像素點進行采樣，具體為：
[0019] 令該星點i在圖像矩陣中的坐標為（x，y)，則對其相鄰像素點（x+n，y)和（x，y+n) 進行采樣，所述η為整數(shù)，且-N彡η彡N ;所述i為正整數(shù)，且1〈 = i〈 = 100 ;所述N為正 整數(shù)，且N彡10 ;
[0020] 采樣得到的序列長度為2N+1，其中，與星點i的y坐標相同的采樣點灰度值序列 為Xi-N，Xi-N+l，Xi-N+2，…，Xi，Xi+l，Xi+2，…，Xi+N-l，Xi+N，與生點?的X坐豐不相問的木樣點灰度值序 列為y^+i，yi^，…，yi，y i+i，yi+2,…，yi+N-i，yi+N ;若關于星點i點對稱的采樣點的灰度 值相問，即 Xi-I - Xj+I, Xj-2 - Xi+2> Xi-3 - Xi+3…，Xi-N - Xi+N，且 Yi-I - Yi+l) Yi-2 - 5^1+2) 5^1-3 - yi+3…，yg = yi+N，則選取該點序列為目標點序列，進入步驟（4);否則返回步驟（3)對下一 幀圖像進行處理；
[0021] (4)根據(jù)面陣C⑶相機的滿阱電子數(shù)M以及圖像量化位數(shù)bs將 步驟（3)中得到的目標點序列轉(zhuǎn)換為電子數(shù)序列，所述電子書序列為 …"匕，所述轉(zhuǎn)換過程具體由公式：
[0022]

【權(quán)利要求】
1. 一種星載面陣CCD相機點擴散函數(shù)的測量裝置，其特征在于：包括光源系統(tǒng)（103)、 靶標（104)、平行光管（105)、面陣CCD相機（106)和圖像采集系統(tǒng)（107); 靶標（104)放置在平行光管（105)的焦面上，平行光管（105)與面陣CCD相機（106) 的光軸共軸放置； 光源系統(tǒng)（103)產(chǎn)生均勻光后照射在靶標（104)上，然后經(jīng)過平行光管（105)和面陣CCD相機（106)，成像在面陣CCD相機（106)的焦面上，形成靶標（104)的圖像，圖像采集系 統(tǒng)（107)對面陣CCD相機（106)的成像結(jié)果進行采集和處理； 靶標（104)包括一百個星點、第一定標線（101)和第二定標線（102)，所述一百個星點 分別為第一星點（1)至第一百星點（100); 第一星點（1)至第一百星點（100)形狀和尺寸均相同，且均為正方形透光孔，正方形透 光孔的邊長為標準長度L的0. 9倍；所述L= (f^/Od，其中fMl為平行光管（105)的焦 距，f?！窞槊骊嘋⑶相機（106)的焦距，d為面陣C⑶相機（106)中探測器單元間距；第一定 標線（101)和第二定標線（102)均為長方形透光孔，長方形透光孔的寬度為標準長度L的 〇. 9倍，長度大于標準長度L的5倍； 第一星點（1)至第一百星點（100)、第一定標線（101)和第二定標線（102)在靶標 (104)平面直角坐標系XOY中的坐標分別為： 第一定標線（101)長方形結(jié)構(gòu)左上角的頂點坐標為（-20噸，80噸）；第二定標線（102) 長方形結(jié)構(gòu)左上角的頂點坐標為（80噸，200噸）；第一星點（1)的坐標為（0,0);第i星點 的坐標為（a· 19. 9 ·L，b· 19. 9 ·L)，其中a為i除以10的余數(shù)部分減1，b為i除以10 的整數(shù)部分，l〈i〈l〇〇 ;第一百星點（100)的坐標為（179. 1噸，179. 1噸）點；所述星點坐標 為星點正方形結(jié)構(gòu)左上角的坐標； 所述靶標（104)平面直角坐標系XOY的坐標原點0為第一星點（1)正方形結(jié)構(gòu)位于左 上角的頂點，X軸的正方向為第一星點（1)指向第二星點（2)，Υ軸的正方向為第一星點（1) 指向第十一星點（11); 所述的靶標（104)大小長度和寬度均大于300 ·L。
2. 根據(jù)權(quán)利要求1所述的一種面陣CCD相機點擴散函數(shù)的測量裝置，其特征在于：所 述的靶標（104)中的第一星點（1)至第一百星點（100)、第一定標線（101)和第二定標線 (102)透光，其余部分不透光，透光部分與不透光部分對比度不低于100:1。
3. -種利用權(quán)利要求1所述測量裝置的面陣CCD相機點擴散函數(shù)的測量方法，其特征 在于步驟如下： (1) 圖像采集系統(tǒng)（107)采集靶標（104)圖像的同時不斷調(diào)整靶標（104)的位置，直至 靶標（104)的第一定標線（101)位于圖像矩陣同一列上，同時第二定標線（102)位于圖像 矩陣同一行上，固定此時靶標（104)的位置； (2) 設置面陣CCD相機（106)的曝光時間，使得第一星點（1)至第一百星點（100)在圖 像中灰度值的范圍為：[〇. 5Κ，0. 6Κ];利用圖像采集系統(tǒng)（107)采集得到P幀當前曝光時間 下的靶標（104)的圖像，所述K為飽和灰度值，P為大于100的正整數(shù)； (3) 逐點對第X幀靶標（104)圖像上的星點i的相鄰像素點進行采樣，具體為： 令該星點i在圖像矩陣中的坐標為（X，y)，則對其相鄰像素點（x+n，y)和（X，y+n)進 行采樣，所述η為整數(shù)，且-N彡η彡N;所述i為正整數(shù)，且1〈 =i〈 = 100 ;所述N為正整 數(shù)，且N彡10; 采樣得到的序列長度為2N+1，其中，與星點i的y坐標相同的采樣點灰度值序列為Xi,XiL1,Xi_N+2，…，Xi,xi+1,xi+2，…，Xm,xi+N，與星點i的X坐標相同的采樣點灰度值序列為 yi_N，yi_N+i，yi-N+2,…，yi，yi+i，yi+2,…，yi+N-i，yi+N;若關于星點i點對稱的采樣點的灰度值相同， 艮PXi-I - Xi+1) Xi-2 - Xi+2, Xi-3 - Xi+3...，Xi-N - Xi+N?-S-Yi-I - Yi+1) Yi-2 -Ii+2,h- 5^+3...，Yi-N =yi+N，則選取該點序列為目標點序列，進入步驟（4);否則返回步驟（3)對下一幀圖像進行 處理； (4) 根據(jù)面陣CCD相機（106)的滿阱電子數(shù)M以及圖像量化位數(shù)bs將 步驟（3)中得到的目標點序列轉(zhuǎn)換為電子數(shù)序列，所述電子書序列為 -.C.，-C-…一-二，…，C ，=^^ XM2々s-1 給出，其中X為采樣點的灰度值，Xe為電子數(shù)； (5) 對步驟（4)中得到的電子數(shù)序列進行歸一化處理，得到歸一化電子數(shù)序列后利用 高斯函數(shù)Ν(μpμ2,σσ2)進行擬合，得到第X巾貞圖像高斯分布點擴散函數(shù)方差σi和 0 2 ； 所述歸一化電子數(shù)序列由公式： 一Λ -瓦 / +Λ 給出，式中f為歸一化的電子數(shù)序列，^ =Σ<; -"Λ. (6) 對當前曝光時間下P幀圖像中的每一幀圖像按照步驟（4)?步驟（5)的方法進行 處理，得到每一幀圖像的高斯分布點擴散函數(shù)方差，最終得到高斯分布點擴散函數(shù)方差序 列％:和，對％:和作直方圖分布統(tǒng)計，得到直方圖序列，并將該直方圖序列擬合 得到高斯分布函數(shù)Ν(μhl,μκ，Ohl,Oh2),N(0,0,Ohl,Oh2)即為面陣CCD相機（106)的點 擴散函數(shù)； (7) 利用步驟（6)中得到的面陣CCD相機（106)的點擴散函數(shù)對面陣CCD相機（106) 的成像數(shù)據(jù)進行圖像復原，提高遙感圖像質(zhì)量。
【文檔編號】G01M11/02GK104236859SQ201410484237
【公開日】2014年12月24日 申請日期:2014年9月19日 優(yōu)先權(quán)日:2014年9月19日
【發(fā)明者】劉薇, 高慧婷, 何紅艷, 鮑云飛, 邢坤, 曹世翔, 江澄, 王殿中, 李巖, 李方琦 申請人:北京空間機電研究所