毛詩(shī)漫 李敏

（1.河南省測(cè)繪工程院，河南 鄭州 450000；2.河南省基礎(chǔ)地理信息中心，河南 鄭州 450000）

PCI GXL在地理國(guó)情普查標(biāo)準(zhǔn)時(shí)點(diǎn)核準(zhǔn)中的應(yīng)用

毛詩(shī)漫1李敏2

（1.河南省測(cè)繪工程院，河南 鄭州 450000；2.河南省基礎(chǔ)地理信息中心，河南 鄭州 450000）

本文介紹了為什么使用GXL進(jìn)行影像糾正以及糾正的流程和方法，并通過(guò)GXL在第一次全國(guó)地理國(guó)情普查標(biāo)準(zhǔn)時(shí)點(diǎn)核準(zhǔn)中的應(yīng)用對(duì)糾正精度進(jìn)行了分析,對(duì)出現(xiàn)的問(wèn)題給出了解決方法，并列出了實(shí)際應(yīng)用的一些方法。

PCI GXL；地理國(guó)情普查；衛(wèi)星遙感影像；標(biāo)準(zhǔn)時(shí)點(diǎn)核準(zhǔn)

1 有關(guān)GXL的說(shuō)明

1.1 GXL自動(dòng)化處理系統(tǒng)

GXL（Geoimaging Accelerator，GXL）是PCI公司面向海量影像自動(dòng)化生產(chǎn)提出的新一代解決方案產(chǎn)品，主要用于衛(wèi)星遙感影像數(shù)據(jù)的自動(dòng)化生產(chǎn)。在保證用戶處理精度的情況下，大大提高數(shù)據(jù)的處理速度和處理效率。

PCI地理成像加速器（GXL）將高性能的計(jì)算能力與PCI Geomatics的專業(yè)處理技術(shù)相結(jié)合，使數(shù)據(jù)處理在速度和效率上都得到了顯著提高。通過(guò)定制自動(dòng)化的GXL工作流，能夠大大提升傳統(tǒng)的影像處理能力，從而提高處理衛(wèi)星影像的實(shí)力。

1.2 GXL的系統(tǒng)架構(gòu)

GXL系統(tǒng)根據(jù)客戶處理數(shù)據(jù)的規(guī)模和處理速度需求設(shè)計(jì)系統(tǒng)處理節(jié)點(diǎn)數(shù)目。通常系統(tǒng)設(shè)計(jì)中，GXL處理服務(wù)器包括日常工作響應(yīng)的常規(guī)服務(wù)器和緊急情況下啟動(dòng)的備份服務(wù)器兩部分。系統(tǒng)中大數(shù)據(jù)量的自動(dòng)化處理在GXL服務(wù)器中完成，數(shù)據(jù)的質(zhì)量檢查在工作站上實(shí)現(xiàn)。

GXL系統(tǒng)是分布式的處理系統(tǒng)，系統(tǒng)管理員可以通過(guò)GXL管理工作站執(zhí)行處理任務(wù)、處理任務(wù)監(jiān)控和終止處理任務(wù)。GXL系統(tǒng)也是一個(gè)開(kāi)放的體系架構(gòu)，可以和用戶的數(shù)據(jù)存儲(chǔ)系統(tǒng)結(jié)合搭建形成集數(shù)據(jù)存儲(chǔ)、數(shù)據(jù)處理和分發(fā)為一體的完整系統(tǒng)。GXL使用中心數(shù)據(jù)庫(kù)儲(chǔ)存所有作業(yè)的歷史記錄和狀態(tài)。

1.3 GXL設(shè)計(jì)的優(yōu)點(diǎn)

GXL將專業(yè)的地理圖像引入開(kāi)放式體系結(jié)構(gòu)及標(biāo)準(zhǔn)化的硬件組件中，根據(jù)用戶數(shù)據(jù)處理量需求，同時(shí)支持多個(gè)CPU、CPU多核以及多GPU的數(shù)據(jù)并行處理和分布式處理。

通過(guò)靈活穩(wěn)定的PCI ProLines技術(shù)，將PCI傳統(tǒng)的優(yōu)秀圖像處理算法進(jìn)行了流程化封裝，用戶可以根據(jù)工作需要自行選擇所需的功能模塊搭建個(gè)性化的ProLines。同時(shí)兼容Windows和Linux操作系統(tǒng)。

PCI GXL設(shè)計(jì)的優(yōu)點(diǎn)：

（1）采用GPU并行處理來(lái)加速圖像處理速度；

（2）支持CPU多核并行計(jì)算；

（3）各模塊可隨意組合應(yīng)用靈活易用；

（4）支持用戶定制開(kāi)發(fā)；

（5）基于Web的用戶界面，支持遠(yuǎn)程訪問(wèn)，作業(yè)模塊啟動(dòng)時(shí)間的預(yù)定；

（6）全自動(dòng)化處理，可人工干預(yù)；

（7）硬件選擇多樣，支持筆記本電腦、臺(tái)式機(jī)、刀片機(jī)群；

（8）支持當(dāng)前主流的所有數(shù)據(jù)格式，支持國(guó)產(chǎn)衛(wèi)星數(shù)據(jù)；

（9）提供顯示處理效果和處理結(jié)果報(bào)告。

當(dāng)今可以糾正影像的軟件也不少，各有特點(diǎn)，以下列出幾種軟件的對(duì)比表（表1）。

表1 GXL與其他軟件對(duì)比表

2 GXL在時(shí)點(diǎn)核準(zhǔn)中的技術(shù)方案

2.1 原始數(shù)據(jù)

原始數(shù)據(jù)是ZY02C、ZY3、KOMPSAT-3、Geoeye-1等的衛(wèi)星影像資料。

2.2 已有圖件資料

（1）2014年制作的地面分辨率為0.5m的河南省第一次全國(guó)地理國(guó)情普查正射影像（DOM）。

（2）2011年1∶50000DEM成果(河南省測(cè)繪地理信息局提供)。

（3）河南省第一次全國(guó)地理國(guó)情普查控制點(diǎn)庫(kù)。

2.3 數(shù)學(xué)基礎(chǔ)

(1)平面坐標(biāo)系統(tǒng)[1]

大地基準(zhǔn)：2000國(guó)家大地坐標(biāo)系。

投影方式：高斯―克呂格投影。

分帶方式：6°分帶。

(2)高程基準(zhǔn)[1]

1985國(guó)家高程基準(zhǔn)。

(3)成圖地面分辨率

DOM地面分辨率：ZY3為2m、GF1為2m、TH為2m、ZY2為 2.3m、GF2為 0.8m、Geoeye-1為 0.4m、KOMPSAT-3為1m。

2.4 資源投入概述

影像處理軟件: PCI GXL；

投入人員：1名；

投入設(shè)備：1臺(tái)服務(wù)器。

圖1 GXL生產(chǎn)處理流程

2.5 GXL生產(chǎn)處理流程

GXL采用的是Web風(fēng)格的后臺(tái)處理操控界面，操控界面簡(jiǎn)潔易懂，處理流程如圖1:

2.5.1 影像導(dǎo)入

在GXL Web界面中啟動(dòng)影像導(dǎo)入功能，輸入原始影像存放位置和影像輸出目錄。提交新作業(yè)后系統(tǒng)同時(shí)進(jìn)行所有影像的讀入操作，讀入的同時(shí)自動(dòng)完成分塊影像的合并、波段組合等工作，轉(zhuǎn)換成PIX格式，采用GDB技術(shù)進(jìn)行統(tǒng)一管理，方便后期應(yīng)用。

GXL數(shù)據(jù)讀入功能模塊支持當(dāng)前主流的所有數(shù)據(jù)格式，支持應(yīng)用RPC模型和嚴(yán)格軌道模型方式讀入數(shù)據(jù)。支持國(guó)產(chǎn)衛(wèi)星數(shù)據(jù)，如ZY02C/GF1/ TH-01等。

2.5.2 自動(dòng)采集控制點(diǎn)

GXL提供了根據(jù)參考影像（DOM）、參考矢量、控制點(diǎn)庫(kù)和GCPs/TPs文本文件等幾種控制點(diǎn)來(lái)源，以滿足不同的控制資料情況。

在本區(qū)地理國(guó)情普查標(biāo)準(zhǔn)時(shí)點(diǎn)核準(zhǔn)中有地面分辨率為0.5m的河南省第一次全國(guó)地理國(guó)情普查正射影像（DOM）、河南省第一次全國(guó)地理國(guó)情普查控制點(diǎn)庫(kù)。適合采用參考影像（DOM），也可導(dǎo)入控制點(diǎn)庫(kù)。

對(duì)于大多數(shù)現(xiàn)代的傳感器所分配的RPC是非常良好準(zhǔn)確的，因此在控制點(diǎn)采集的時(shí)候使用有理函數(shù)模型是最穩(wěn)固的模型。GCP的個(gè)數(shù)根據(jù)需要來(lái)定，如果RPC的精度滿足要求的話，不需要采集太多GCP，只要分布均勻即可。搜索半徑不可過(guò)大或過(guò)小，過(guò)大或過(guò)小可能會(huì)匹配不到點(diǎn)，一般搜索半徑是原始影像和參考影像的誤差與參考影像的分辨率的比值的1.5倍或2倍。

2.5.3 區(qū)域網(wǎng)平差

平差是通過(guò)在重疊區(qū)域采集連接點(diǎn)（tie point），調(diào)整相鄰影像的接邊，進(jìn)而提高區(qū)域的整體精度。平差模塊根據(jù)用戶指定的搜索策略、搜索半徑（搜索半徑和采集控制點(diǎn)時(shí)的搜索半徑的方法相同）、連接點(diǎn)個(gè)數(shù)進(jìn)行自動(dòng)采集，采集后的連接點(diǎn)可以查看其分布和誤差，還可對(duì)連接點(diǎn)進(jìn)行編輯修改。區(qū)域網(wǎng)平差任務(wù)提交后，系統(tǒng)會(huì)自動(dòng)查找影像重疊像對(duì)，并按設(shè)定參數(shù)采集連接點(diǎn)，如果提交影像較多，GXL將自動(dòng)將一個(gè)大區(qū)域分成若干個(gè)小區(qū)域,同時(shí)進(jìn)行連接點(diǎn)采集，采集完成后將小區(qū)域重新合并并進(jìn)行平差，這樣既能保證成果精度，又能大幅提高生產(chǎn)效率。平差的結(jié)果主要體現(xiàn)在對(duì)RPC的改善上。

2.5.4 全色正射校正

PCI GXL既可使用嚴(yán)格模型（也稱Toutin模型）也可使用RPC模型技術(shù)進(jìn)行衛(wèi)星影像正射校正。利用采集的足夠且符合要求的GCPs/TPs、足夠精度的DEM數(shù)據(jù)，PCI可以生產(chǎn)高達(dá)亞像元級(jí)精度的正射產(chǎn)品。軟件在進(jìn)行全色校正的同時(shí)會(huì)掃描目錄中的多光譜影像進(jìn)行匹配和校正。

2.5.5 銳化融合

PANSHARP是一個(gè)自動(dòng)融合工具，該算法是一種基于最小二乘法在原始的多光譜、全色影像間達(dá)到最佳近似灰度值關(guān)系[2]，從而達(dá)到最佳色彩組合的融合方法。

2.5.6 鑲嵌分幅

后期如果需要分幅圖時(shí)可在GXL中鑲嵌，用戶必須先進(jìn)行鑲嵌預(yù)處理以生成包含勻色效果和接邊線的XML文件，然后才能進(jìn)行鑲嵌，影像間由于不同的光照和視角，同一地物在重疊影像中看起來(lái)存在差異，得不到理想的效果，在鑲嵌預(yù)處理的參數(shù)設(shè)置的勻色方法新增加了“捆綁”，使用“捆綁”的方法勻色使影像最大化地消除差異。用戶也可以對(duì)鑲嵌預(yù)處理自動(dòng)采集的接邊線進(jìn)行手動(dòng)編輯，并應(yīng)用到影像正式鑲嵌結(jié)果中（圖2）。

圖2 GXL自動(dòng)鑲嵌效果圖

鑲嵌結(jié)果可以直接根據(jù)文本坐標(biāo)或者矢量shp數(shù)據(jù)等方式進(jìn)行分幅裁切。根據(jù)需要選擇輸出格式。

3 遇到的問(wèn)題及解決方法

以上流程適合大部分的衛(wèi)星影像，在作業(yè)過(guò)程中控制點(diǎn)和平差點(diǎn)無(wú)需手工調(diào)整基本都能滿足精度要求。但是有些影像的RPC精度不高或原始影像有局部變形，則需要調(diào)整一下方法，如有些TH1A、GF2的RPC的精度不能滿足作業(yè)精度，ZY02C的影像內(nèi)部有的局部變形，以下是針對(duì)這些情況的處理方法。

3.1 TH1A的處理方法

由于TH1A的數(shù)據(jù)自帶RPC的精度不能滿足要求，如果利用自帶RPC會(huì)出現(xiàn)大部分地方匹配不上點(diǎn)且校正后的影像和參考影像誤差過(guò)大，所以在處理過(guò)程中就不能再用它自帶的RPC來(lái)計(jì)算，而是要從GCP中來(lái)反算一個(gè)RPC。影像導(dǎo)入之后在Focus中打開(kāi)衛(wèi)星影像刪除模型段（Rational Functions Model Segment），刪除之后再進(jìn)行控制點(diǎn)的自動(dòng)采集，在控制點(diǎn)采集參數(shù)中要選上從GCP中計(jì)算RPC而且控制點(diǎn)個(gè)數(shù)也要多，通常設(shè)1000個(gè)或更多。控制點(diǎn)采集完成后進(jìn)入桌面版（Geomatica）中進(jìn)行手工干預(yù)控制點(diǎn)。在桌面版中選擇Optical Satellite Modelling和Rational Function(compute from GCPS)，導(dǎo)入自動(dòng)采集的控制點(diǎn)修改或刪除誤差大的點(diǎn)，控制點(diǎn)調(diào)整完成就可以直接糾正輸出正射影像，多光譜影像采用同樣的方法做一遍。

3.2 GF2的處理方法

GF2相對(duì)于TH1A要簡(jiǎn)單多了，GF2只是出現(xiàn)一些局部的套合不太好，只需要在自動(dòng)采集控制點(diǎn)的參數(shù)設(shè)置里把數(shù)學(xué)模型的次數(shù)選為2階，也要從GCP中計(jì)算RPC，除這兩處的設(shè)置不同外，其他均和正常流程相同即可。

3.3 ZY02C的處理方法

ZY02C的衛(wèi)星影像和TH1A、GF2的不太相同，ZY02C_HR1左邊的原始影像本身變形，內(nèi)部局部套合不好，出現(xiàn)原始影像變形的情況下用一般的方法是糾不過(guò)去的，碰到這種情況只有自動(dòng)和手動(dòng)相結(jié)合才能使糾正的影像和參考影像套合。先使用正常的流程做一遍，如果生成的正射不能滿足精度要求，利用生成的正射影像再自動(dòng)采集一遍控制點(diǎn)，控制點(diǎn)采集結(jié)束進(jìn)入桌面版新建工程時(shí)要選擇用多項(xiàng)式（Polynomial）,把正射采集的控制點(diǎn)導(dǎo)入手動(dòng)修改誤差大的點(diǎn)，修改完成后再在此工程下自動(dòng)采集控制點(diǎn)，控制點(diǎn)個(gè)數(shù)要盡量多一些。把采集完的控制點(diǎn)再導(dǎo)入波段樣條（Thin Plate Spline）的工程中進(jìn)行正射糾正，如果生成的正射仍有套合不住的地方可返回工程在變形處手動(dòng)加點(diǎn)，經(jīng)過(guò)反復(fù)加點(diǎn)，直到結(jié)果滿足精度要求。

4 質(zhì)量及效率

4.1 平面精度檢查

方法一：將參考影像與正射成果疊加，利用“卷簾”目視檢查（圖3）。

方法二：打開(kāi)采集的控制點(diǎn)檢查控制點(diǎn)的誤差是否滿足精度要求（圖4）。

圖3 正射影像與參考影像疊加“卷簾”

圖4 控制點(diǎn)誤差表

表中的單位為m，從表中可以看出X的中誤差為1.33m、Y的中誤差為1.73m，控制點(diǎn)最大誤差3.9m，地理國(guó)情普查標(biāo)準(zhǔn)時(shí)點(diǎn)核準(zhǔn)中的精度要求是中誤差不大于4m，最大誤差不大于兩倍的中誤差，表中控制點(diǎn)的誤差完全滿足精度要求[3]。

4.2 GXL的生產(chǎn)效率

下面的工程中共16景影像，時(shí)間統(tǒng)計(jì)是采用單節(jié)點(diǎn)CPU處理效率，若采用多節(jié)點(diǎn)處理效率可成倍提高。

表2 耗時(shí)統(tǒng)計(jì)表

4.3 融合后效果

PANSHARP是一個(gè)自動(dòng)影像融合工具，該算法是一種基于最小二乘法在原始的多光譜、全色影像間達(dá)到最佳近似灰度值關(guān)系，這對(duì)地物真實(shí)色彩的恢復(fù)和解譯精度的提高具有重要的意義。融合狀況見(jiàn)圖5。

4.4 藍(lán)波段的模擬

圖5 融合前后多光譜影像對(duì)比

ZY02C多光譜的原始影像缺少藍(lán)波段，直接融合后的影像因缺少藍(lán)波段不易詳細(xì)區(qū)分地物，后期為作業(yè)方便可轉(zhuǎn)自然色或進(jìn)行藍(lán)波段的模擬，PCI GXL桌面版中的藍(lán)波段模擬完成后的結(jié)果為含有藍(lán)波段的多光譜成果，圖6為運(yùn)行界面。Model Text中的算法如下：

圖6 藍(lán)波段模擬

%1=%1

%2=%2%3=13.9054+0.873393*%1-0.138254*%2-0.100616*%3

圖7分別是融合后影像、PCI GXL桌面版的藍(lán)波段模擬和用erdas轉(zhuǎn)的自然色。

圖7 藍(lán)波段模擬效果圖

5 結(jié)論與建議

對(duì)于國(guó)外或國(guó)產(chǎn)的衛(wèi)星影像數(shù)據(jù)，采用GXL區(qū)域網(wǎng)平差的方法，不但可以滿足地理國(guó)情普查標(biāo)準(zhǔn)時(shí)點(diǎn)核準(zhǔn)的精度要求，而且可以大幅提高生產(chǎn)效率。采用GXL自動(dòng)化生產(chǎn)系統(tǒng)，明顯提高了生產(chǎn)效率，降低勞動(dòng)強(qiáng)度，節(jié)省了人力。采用一個(gè)處理節(jié)點(diǎn)，處理100景的衛(wèi)星影像數(shù)據(jù)，所有流程需要1.5天就可以完成，而兩個(gè)熟練的作業(yè)員處理完同樣的數(shù)據(jù)需要半個(gè)月左右，效率提高至少10倍。

對(duì)于地理國(guó)情普查標(biāo)準(zhǔn)時(shí)點(diǎn)核準(zhǔn)等需要處理海量影像數(shù)據(jù)的生產(chǎn)項(xiàng)目，最好采用GXL等自動(dòng)化并行處理系統(tǒng)。

[1]國(guó)務(wù)院第一次全國(guó)地理國(guó)情普查領(lǐng)導(dǎo)小組辦公室.數(shù)字正射影像生產(chǎn)技術(shù)規(guī)定：GDPJ 05—— 2013[S].北京：國(guó)務(wù)院第一次全國(guó)地理國(guó)情普查領(lǐng)導(dǎo)小組辦公室，2013.

[2]國(guó)務(wù)院第一次全國(guó)地理國(guó)情普查領(lǐng)導(dǎo)小組辦公室. 遙感影像解譯樣本數(shù)據(jù)技術(shù)規(guī)定：GDPJ 06——2013[S].北京：國(guó)務(wù)院第一次全國(guó)地理國(guó)情普查領(lǐng)導(dǎo)小組辦公室，2013.

毛詩(shī)漫 (1980― )，女，漢族，工程師，主要從事遙感和地圖制圖工作。E -mail:237966883@qq.com