摘要：介紹一種基于ADSP21060和Virtex II的星載圖像處理系統。分析了圖像處理系統的功能和任務，給出了處理系統的硬件結構、FPGA的功能模塊、DSP的軟件框架和模塊。通過地面原理樣機開發，驗證了系統設計的正確性和高效性。
圖像處理系統多采用DSP陣列、DSP加FPGA/CPLD或單由FPGA/CPLD器件等方式構成。采用DSP陣列構成的圖像處理系統，其優點是處理功能可以通過軟件靈活修改，其缺點主要有功耗大、體積大、成本高；采用單由FPGA/CPLD等可編程器件構成的系統，可以避免DSP陣列系統功耗大、體積大的缺點，并可以根據系統要求，在硬件構造上作靈活的配置，但對于復雜的算法，其設計復雜度急劇上升。正因為如此，當前很多設計選擇DSP加FPGA/CPLD的結構來構造圖像處理系統。本空間太陽望遠鏡星載圖像處理系統，正是一種基于DSP和FPGA構造的圖像處理系統。本文討論了系統的功能任務、系統結構、FPGA設計的邏輯模塊、DSP的軟件結構和原理樣機的實驗結果。
1 圖像處理系統的功能和任務
1.1 系統的外圍接口
系統的外圍接口如圖1。處理系統接收來自CCD單元的圖像數據和同步信號，同時向CCD單元提供復位、曝光開始、數據讀出等控制信號；CEU是一個管理機，管理和控制CCD單元、圖像處理單元和壓縮單元，CEU單元向圖像處理系統發送各種命令，如CCD標定、正常的觀測測試及觀測數據處理所需要的數據，如觀測模式、時間等；系統將處理完成的數據送到壓縮單元進一步壓縮。
1.2 系統的圖像處理任務
圖像處理系統的任務是接收CCD傳來的圖像和同步信號，進行圖像積分、輻射校正、幾何校正、磁圖偏陣計算和圖像格式化，然后把數據送到壓縮單元進行壓縮。
(1)圖像積分：也就是圖像的多幀疊加，其目的是提高信噪比和消除圖像的時間相關性以壓縮數據量。空間太陽望遠鏡有三個觀測模式，分別是爆發模式、活動區模式和寧靜態模式。經計算，三種觀測模式下圖像積分的幀數分別是：4幀、24幀、240幀；每個模式又順次觀測六個分量，六個分量分別是：SV1、SV2、SQ1、SQ2、SU1、SU2。爆發模式的觀測時間是30s。30s的時間是這樣安排的：每個分量連續采集4幀，每秒一幀，然后1s進行偏振光學元件的切換，如此直到六個分量觀測完畢。活動區模式由6個連續的爆發周期構成，觀測時間是3min；而寧靜態模式又是10次活動區模式周期構成，觀測周期是30min。
(2)圖像的輻射校正：圖像由于受到大氣擾動、傳感器自身特性、宇宙射線等輻射的影響，圖像的質量會產生退化，必須進行輻射校正。主要考慮CCD的校正(暗流、偏置、平場)和條紋斑點的消除兩個方面。
(3)圖像的幾何校正：圖像的幾何校正，成像過程中引起的幾何畸變的校正；由于系統或隨機因素造成圖像產生幾何畸變，必須進行校正。考慮圖像漂移的校正，采用相關技術，求出活動圖像相對于參考圖像的的漂移量，再采用線性內插累加的方式，實現漂移前后圖像的積分問題。根據太陽米粒的存活壽命，幾何校正只對寧靜模式觀測進行，并且是以3min作為校正時間間隔。
(4)圖像的偏振計算：空間太陽望遠鏡主要進行太陽磁場測量，磁圖的偏振歸算也是重要任務之一。磁場圖像的歸算還可以是數據量壓縮1.5倍。
(5)圖像格式化：給圖像加上圖像頭，應包括：儀器號、時間、觀測模式、觀測分量等。
2 圖像處理系統設計
2.1 系統結構
系統的硬件結構框圖如圖2。FPGA中的硬件模塊，包括積分單元、多路開關單元、兩個SDRAM控制器和串口控制器。DSP的所有總線都伸入到FPGA內部。存儲器由兩個SDRAM單元和一個SRAM單元構成，SRAM只供DSP使用。系統接收CCD傳來的圖像，在積分單元進行圖像積分，積分就緒的圖像再由DSP進行輻射校正、幾何校正、偏振計算和圖像格式化等操作，最后通過DSP的片上串口將圖像送到壓縮單元。通過串口單元系統與管理機單元通信，接受管理和控制。由于圖像積分與DSP的數據處理、數據傳送是同時并行進行的，因此積分器和DSP必須采用相互獨立的存儲區。為了避免圖像在積分器和DSP的存儲區之間來回進行圖像傳送，在設計中采用多路開關進行存儲區乒乓方式切換。
圖2 系統硬件結構圖
2.2 FPGA邏輯模塊
FPGA內部設計的邏輯模塊主要包括圖像積分單元、通路選擇單元、異步通用串口(UART)單元、兩個SDRAM控制和其它附加邏輯。下面分別介紹各塊的結構，由于SDRAM控制器和串口單元(UART)均采用IP核實現，這里不作詳細介紹。
2.2.1 積分單元
積分單元主要完成如下功能：(1)按給定的時序接收CCD的圖像，并保存在寄存器中；(2)積分地址產生；(3)存儲器控制；(4)積分操作；(5)積分擴展位清零，當每個分量的第一幀圖像接收時，無須積分，但必須把存儲器的高位(積分擴展位)清零，其它情況不必清零；(6)每幀圖像積分結束，產生Irq信號中斷DSP，以便由DSP進行工作安排。
積分單元的邏輯框圖如圖3，其輸入輸出如下：(1)CCD輸入數據(CCD data)，同步信號——點同步(ps)、行同步(ls)、幀同步(fs)；(2)DSP數據、地址和控制總線，irq是積分單元發給DSP的中斷信號；(3)與存儲器連接的數據、地址和存儲器控制總線——data，mctrl，addr；(4)控制CCD的信號線——reset、start、read。
圖3充色的方框是DSP存儲器地址映射的寄存器，用于對積分單元進行控制和對DSP數據處理狀態進行標識；無充色的方框則是積分單元的硬件結構，從CCD送來的數據通過數據寄存器(OP1)，作為寄分的第一個操作數，從存儲器送來的數據通過數據寄存器(OP2)作為積分的第二個操作數；從CCD送來的同步信號，送到幀內地址發生器，作為地址產生信號，并生成存儲器需要的控制信號；CCD送來的OP1被寄存，同時產生地址從內存中取出要進行積分的操作數OP2，然后在點同步信號的下降沿，進行積分運算并存回存儲器。
圖3 積分單元結構
2.2.2 多路開關單元
多路開關單元如圖4所示。本單元的功能是完成積分總線INT bus和DSP bus與SDRAM1和SDRAM2乒乓連接結構的切換。切換包括三大總線：數據、地址和控制總線。切換通過一位切換控制信號實現，控制信號是DSP地址映射的一個一位寄存器，由DSP控制乒乓切換。
圖4 多路開關單元
2.3 DSP軟件結構
軟件部分包括系統軟件框架和各處理軟件模塊、中斷服務程序等；程序的流程和數據處理控制，都依賴于由DSP地址映射的寄存器的狀態。
2.3.1 DSP在FPGA中的存儲器映像寄存器及其作用
這一組寄存器都是在FPGA中實現的。作業(指令)采用16位寄存器，由于DSP地址空間較大，為了方便編程，寄存器不采用拼位的方法。系統的工作 受CEU的控制。CEU發送過來的指令格式如圖5(a)，其中DT表示CEU傳來的數據類型，00表示數據、01表示觀測指令、10表示CCD標定、11表示自測。接收時，若是命令則直接寫到JOB或放入命令隊列；若是數據，根據類型存放到相應的位置。Parameters是各種工作狀態參數，如待機/自測要求的參數、測試的要求(內容、測試選擇)、響應信息要求、觀測模式和次數、CCD標定要求的參數、標定的內容等，以便于DSP查看。
圖5(a) CEU命令格式 圖5(b) 系統作業格式
如圖5(b)，當CEU送來命令DT=01時，即為觀測命令。觀測作業的格式如：00表示無作業待機、01表示觀測指令、10表示CCD標定、JT=11表示自測；mode是當前觀測模式標志，001表示爆發模式、010表示活動區模式、100表示寧靜模式；Observ_Num是連續觀測同一種模式的次數標志。
(1)作業類型寄存器(JOB)：標識系統工作內容(作業類型)，如CCD標定、觀測、自測試等；(2)模式寄存器(MODE)：標識觀測的模式，爆發模式、活動區模式、寧靜模式；(3)觀測次數寄存器(OBSN)：用于保存一次作業的觀測次數；(4)寧靜態的觀測次數寄存器(RESTN)：只用于寧靜模式的觀測，當為寧靜態模式時，用它來標志3min觀測的次數；(5)幀計數器(FRMN)：用于計數爆發模式、活動區模式及寧靜模式下的前24幀圖像，以便用于高位清零信號生成和地址高位生成；(6)SW電路控制寄存器(SWCTRL)：控制多路開關的切換；(7)reset寄存器(FRST)：用于DSP控制FPGA中INT單元的時序邏輯復位，每次觀測結束，必須復位一次；(8)清零標識寄存器(ZERO)：由DSP根據FRMN的值修改它，以實現高位的清零任務；(9)分量寄存器(COMPO)：由DSP根據FRMN的值修改它，以實現高位地址的形成；(10)作業結束寄存器(END)：當一個作業結束時，用次寄存器進行標志，以便將程序的PC跳到作業檢查入口；(11)數據狀態寄存器(DAST)：DSP存儲區數據狀態標識：沒有數據、爆發模式的待處理數據、活動區模式的待處理數據、寧靜模式的待處理數據、寧靜模式處理好數據及待送出數據等。DSP根據數據狀態和觀測模式，可以進行相應的處理；(12)相關數據狀態(CORR)：標識相關運算用的數據的狀態。
2.3.2 DSP的存儲器空間分配
DSP的SYSCON寄存器的15～12位(MSIZE)，可以設置以MS0～MS3作為片選信號時，每個分組的空間大小。其公式為：MSIZE=log2(塊大小)-13；取每個塊的大小為16MB，因此MSIZE=1011。
·片內空間：0000，0000H～0007，FFFFH；
·MS0——數據空間：0040，0000H～13F，FFFFH；數據是倒序放置的，最后一個單元放第一個數據；其中：
SQa：40，0000～5F，FFFF SQb：60，0000～7F，FFFF
SUa：80，0000～9F，FFFF SUb：a0，0000～bF，FFFF SVa：c0，0000～dF，FFFF SVb：e0，0000～fF，FFFF
Rest2m：100，0000～11F，FFFF unused：120，0000～13F，FFFF
·MS1——SRAM空間：140，0000H～23F，FFFFH；
·MS2——CCD參數：240，0000H～33F，FFFFH；放置暗流、偏置和平場圖像3幀2048×1024圖像；
·MS3——未使用：340，0000H～43F，FFFFH；
·未分組空間：4040，0000～FFFF，FFFF所有存儲器映射的寄存器。
2.3.3 DSP的軟件結構
DSP的主程序框架如圖6。系統初始化后，進行作業查詢，根據查詢結果，進行相應的處理任務(CCD標定、自測試、觀測)。圖6中僅以觀測作業為例進行說明。進入觀測后，首先由DSP對積分單元作相應的設置，設置完成后，積分器開始單獨地積分工作；DSP接著查詢本次觀測是否結束，若結束回到作業查詢，否則判斷存儲器中數據的 狀態，根據數據狀態，進行相應的處理；處理就緒的數據根據FLAG0的標志，可以向壓縮單元傳送數據。三種模式下的圖像處理程序流圖如圖7，圖7中ccdst( )是圖像輻射標定子程序，format()是格式化子程序，內插運算采用線性內插方法，相關運算采用FFT和逆FFT方法。
圖6 DSP主程序框架
圖7 圖像處理程序模塊流圖
圖8 積分單元中斷服務程序流圖
除了主程序和數據處理程序外，DSP還有一個重要的中斷子服務，就是積分單元每幀圖像積分完畢的中斷服務程序，其流程框圖如圖8。它根據積分的圖像幀數和觀測模式，判斷觀測是否結束，根據不同的情況設置積分器中相應的寄存器，這些寄存器控制著積分器的工作行為和DSP的數據處理方法。圖8中參考圖像和活動圖像的獲取是在采集到的2048×1024圖像的中心，獲取256×256的子塊，也就是說，相關運算是在兩幀256×256的圖像之間進行的；對于圖像傳送操作，考慮到寧靜模式的數據處理，要求的存儲器比爆發和活動區模式要大一倍。為了不把存儲器設計得太大，對于寧靜模式，不再采用乒乓存儲器切換的結構，而是在數據就緒時，利用光學儀器調整的1s時間，將積分就緒的圖像傳到DSP的工作區內。
3 樣機研制和實驗結果
原理樣機的研制過程中，DSP器件是AD公司ADSP21060，開發平臺采用VisualDSP2.0。FPGA器件選用Xilinx公司100萬門的Virtex II系列器件，開發平臺采用ISE6.2。已完成的地面樣機達到了合同的要求。用(接上頁) DSP硬件仿真器對樣機作了軟硬件測試，首先對各接口進行了測試，在保證數據可以正確地接收和發送后，測試系統的功能。這些測試均在一個專門開發的地面支持系統測試平臺上進行。經過單元測試，再與其它單元接在一起進行系統級的聯調和系統級測試。結果表明達到設計目標，可以按要求完成系統的各項功能。
參考文獻
[1] 中科院國家天文臺.德國馬普高層大氣研究所. SSTA相報告[R].北京：國家天文臺，1997；1(2)：1～30.
[2] 中科院國家天文臺，德國航空宇航研究所.SST評估研究報告[R].北京：國家天文臺，1997；30～36.
[3] 王宇舟.單片FPGA圖像預處理系統的設計與實現.微型機與應用，2004；23(4)
[4] 李廣軍，孟憲元．可編程ASIC設計及應用.成都：電子科技大學，2000
[5] 北京理工大學ASIC研究所.VHDL語言100例詳解[M].北京：清華大學出版社，1999
[6] 劉書明，羅軍輝. ADSP SHARC系列DSP應用系統設計[M]. 北京：電子工業出版社，2003
