适普卫星影像处理方案

1、适普公司卫星影像处理解决方案随着分辩率和影像质量的不断提高，高分辩率卫星影像在测绘及相关行业的应用已越来越广 泛。适普公司的VirtuoZo Sat 已很好的实现了关于高分辩率卫星影像处理的支持。一、软件特点1、处理的卫星影像类型全面影像类型影像类型全色影像星下点分辩率SPOT5IRS-P5HRS10X5米HRG : 2.5 米2.5 米ALOSPRISM2.5 米EROS-A1.9 米EROS-B0.7 米OrbView1 米IKONOS0.82 米Geoeye-10.41 米QuickBir0.61 米WorlView0.45 米TH2.5米全色+10米多光谱2、空三作业方便支持半自动和全

2、人工两种方式添加点。3、处理精度高采用鉴于RFM 的区域网平差方案。使用 RPC 参数（直接读取已提供的历和姿态角进行解算）和仿射变幻参数进行定向，只需少量的控制点（单幅影像达到较高的定向精度 4、支持的坐标系全面RPC 参数或利用卫星星4 至 6 个）便能目前支持的坐标系囊括坐标系。WGS84 北京 54 坐标系、西安 80 坐标系、中国 2000 坐标系及自定义二、数据准备卫星影像的原始数据囊括影像文件和PRC参数文件，不同的传感器的PRC文件名也有差异。影像类型Spot5IKONOS,GeoEye-1，P5RPC文件METAATA.IM*_rpc.txt匚 INAGERY.TIF MET

3、AAIA,M_ BAWA.RPCh- BariA.tif样例（必备文件组成）QuickBri WorlViewOrbView-*.rpb-*.pvl08JAN06230iS-PlES-00572S933010_01_P001.ATT _08JANO&023015*P1ES-&0572S?330LO_01_F01EPHOBJANO&023O15-PlBS-C5728933aWl_PM)l.GEQ 08JAN06023015-PlES-a0572S933010_01_i.JM _08JAN06023O15-PlBS-5728932010_01_P0OLRPB 曲08JANO&023015-P1BS

4、-005728933010_01_P001STE 固O8JANOfr023O15*PlBS-(5728933Ol_Ol_P0OL71F 08JANa&023015-PlBS-&0(XOOg45271A5200001001010_()06ia590.tif L. 3VO51023PO000945271AS200001001010_00&l&?BO,pvlALOS-_RPC-*（PRCIM G-03-ALPS MB044052 925-O1B1R文件分块存储， 由 summary.txt 像时请选择summary.txt）IM G-04-ALPSMB044052925-01B1_BIMG -05-

5、ALPS MB044052 925-01Bl_B IMG- 06-ALPtMB0440525-01Bl_B LE- ALPSMB44052925-01B1_BR PC-03-ALPS M &CH4O52925-01Bl_B RPC-M-ALPEMB(M4n52925-OlBL_B RPC- 05-ALPSMBQ4452925 O1B1_BRPC-06-ALPS M &04405 292 5-01B L_ &RPC-AI_PSMB04405225-01B1_BSLinnmary.txtSU-ALFSMB04452925-01E1_B TRL-ALPSMB04405395-01B1_B VOL-A

6、LPSME4405292501B1 BEROSTH*.rpc*.passSXZ+GFB+GP*.rpc.MBTVe2161&63-lAUbrtkt；f _MBT12161S3.pa 匚MaTl.e2161063.fpcB_SXZ_1_O8 14-136 .RAWB_XZ_l_08-14-134.rptB_SXZ_l_08-14-l36.)cmlB_GP_08_014_136.中cB_GP_O0_O14_136.T1FB_GP_0S_014_136.xml B_GFB_OB_014_136.rpcB.GFe0S_O14_136fTIF B G FB 08 014 136jcmlRPCLINE_O

7、rF: 4-5985,82 pixels+5952.4G pixe sLATOFF: +30,57563548 egreesLCNG_OrF: 4-101 40477032 ecrees HEIGHTOFF: +700.000metersLINESCALE:SAMPSCALE：+9231口 5 pixels+9BS6.S1 pixelsLATSCALE:+0.18000000 egreesLCNGSCALE:+0.19500000 egreesLINENUKCOEFFlHEIGHT SCALE: 4-500.000 metersLINENUMCOEFFS+1.958730240027644E-

8、04-N 【NU匸C 【3611667 630373 E-01一&LINE-NUK-COEFF_4307726057910314E-01 -LINENUHC0EFF52.18350596C8S0724E-03L INE【NU匸 8EFF【5+5 3 6599508561947-02 LINENUMCOEFFj12167121516405E-3 LINENUKCOEFF8+2.86S9113993S9206E-02LINE-NUK_COEFFg+6.69522S93757319E-03十6.61426909057129&E-02 LINE NUM COEFF 10:LINE_NUH_COEFf_

9、11:-2 103313692078135E-04LmEjjU 匸 COEFFp :-5 *913 63 4632529452-06LINENUMCOEFF13:LINE-K(JK-COEFF-i5:-1.435521O637B1549E-0S-1.57807590026760-04-13S1651962131069E-04LINENUMCOEFFl6:+7 640599772S565E-04 LINENU 匸 C OEFFl?：-6.867529676123316E-O5LINE-NUK-COEFF_18:-221083294596713OE-OSLINE_NUM_COEFF_19+1 40

10、405217287949BE-05 LINE 0: -866B49B305394E-07N-EF-1: 41.: OOOOOOOOOOOOOOE+OOE N一马.329588103769903E-02LINEEMCOEFF_4:LINEEZCOEFFS:1*037352298081656E-1-3.422620810240033E-02+3LINEENCOEFFG:LINE3044483466701602E-04ENCOEFF7:-LTNEE:N8EFFB:+2.40043S758372939E-04LINEENCOErF1 311?81726180937E-0S+1T:TNEENCOEFIC

11、I727701314203965E-04-LTNEENCOEFF8.926932463381600E-0411LINEENCOEFFiS LINEENCOEFF:L3 LINEENCOEFFH+9,2683539552905-05-2072Se5532705432E-07-1.12193423S451222E-07-l,411625901401355E-06-LINEENCOEFFISLINE【E 匸COEFF1.7217943239252lEOG16LINEENCOEFFl?+1.012777647990264E-06 LINKEENCOEFFIS908842672488080E-07 LI

12、NEE 匸C0EFF【19L+.454120828626039E-0S JEEN3： OEFF)0+1.703707091275895E-07-1soPNUHCOErri：匸+3,3800156387142-07-0EFF 【2:COEFF3 :+1*458288823041920E-03+7988261872506678E-1-1.952847294878630E-01三、处理流程1所示建立工程 数据准备4生产图1卫星影像处理流程图1、新建工程线相关参数；2、数据准备EMOM和等高RPC其原理如下：Wallis果；旋转：所选旋转角度的旋转方向为沿顺时针方向，设置此系数的主要目的在转（SPOT

13、5-HRS异轨像关于），关于于同轨立体像关于，左右影像应902所示:左右影像同时沿Q 顺时軒方向脸转Q图 2 同轨影像旋转 90 度原理示意图说明：卫星轨道方向一般为南北方向，因此同轨立体像关于为上下重叠关系，这不利于体观测，立体像关于左右影像同时沿顺时针旋转90度，上下重叠关系变成了左右重叠关系X丫致，使地物各个方向的摄影比例尺一致，从而比较有真实感。SPOT5-HRS200%3所示：5m10m，因扫荷方向-行方向飞行向方1个像素10m2个像素10m,向则每个像素 5m图 3 影像拉伸原理图SP0T5-HRS 同轨立体像关于，旋转角度为 90 度，X 方向缩放和Y 100% 200%4 所示

14、。十扫描方向分浙率 10m图 4 SPOT5-HRS像预处理示意图3、空三定向量测控制点的像点坐标（无控制点时添加一定数量的连接点）进行绝关于定向，解算仿纠正系数。4、4 产品生产利用IGS 进行数字化测图，并且生成 EM、等高线和 OM。四、ALOS数据案例1、新建工程设置投影坐标系，根据控制点的投影坐标来设置，与控制点一致。2、引入影像选择影像类型为 ALOS vz tif 均可，ALOSsummary.txt文件引入输出目录：|H： ataSATALOS-TestImages输入影像測览rirairgi影像类型PRISM 二式VZ)RPC 错I I 移I 除一切I任务勺All Files

15、 (.)时紛间3012-11-29 141/2 a Elinifnary. tst 开始处理9 2012-11-2914影豫连播0 2012-11-29 3012-11-29 1459:1859:35生咸祇文件0 度：K 方向缩锁2012-11-29142012-11-29 14紛S012-L1-29 14生成砂 t 文件完成处連2/2 - Etimmary. titt ：幵始处理2012-11-291500:272012-11-2915翱憬连接.生咸眈文件0 片方向缩锁o 2012-11-29 152012-11-29 151:1601:16生咸Pt 文件* ” 任务列康一.处理消息完成处理

16、一切廳像！处理 |退出、设置影像列表添加影像，根据引入影像时设置的输出影像类型来选择。、设置模型列表Hi2 钉 4 韜 刍辰阴aanitae EJWIS EZAMCS UKt- *wIt闵谢匚尽 t_； WUta *I MB.p RIMff IiTMM、定向添加控制点，多条带的话需要在条带间添加连接点。接点，进行疋向修正 rpc 参数。该数据没有控制点，均匀添加了连察看定向报告，像方小于半个像元，物方控制点的精度要求参照国家规矩。ALOS-Test.rep 记事本1 口 II回文件旧巒旧柜式(O)亘看 M 箱旳(H)験!評烁 加*2工工辔H: at aSATALOS -Test 0r i en

17、t at i onALO S-Test. pro ALOS-Test unknown工用称单位权中误差二0.138 (像元) 像点坐标残差;其 口 aeH型2 10 0 0 0J连 口点M点 点其 口 定 点模型 15011 - 951195110002-0.2200, 0000. 4740.000951100030.009-0. 000-0. 01B-o. ooo150110002-0.043-0. 0000. 0910. 000150110003-0.076-0. 0000.1.64-0. 000150110004-0.035-0. 0000.1.S4-0. 0001501100050.1

18、150. 000-0. 247-0. 000150110006u. 142-0. 000-0. 3040.0006、模型裁切，一般选择“原始影像裁切”7、定义裁切方式L4Wg C W jOMt&裁切后就可以用模型进行测图，匹配EM等。H:ataSAnALOS-TestMoelsG48R)35026雋（臼童看（V）工具（T;帮旳H）包含到库中-tt窖称Images修改日期2012/11/915:24文衣谊问的链置k ProuctG43F03502S1 ProGeosys.ge2012/11/2915:24文磁2012/11/25 15：24MR :3012/11/9 15:12G8S9、生产 E

19、M EM，速度较慢。SM，是一个立体像关于匹配一整块另一种是模型匹配 EM，裁切模型后，关于小块的模型进行匹配EM，速度快。EM 匹配完后进行编辑，出 EM 成果。10、生产正射影像按默认纠正方法纠正即可。可以添加修改任务11、可将 EM 叠加正射影像察看地形是否匹配。用【生产工具】打开em，然后加载正射影像。EM】菜单下的【EM 三维显示2、选择文件菜单下的 “同轨影像拼合五、WorlView数据案例WorlVIew 数据的一景影像通常被裁切成多块小影像存储，在定向添加连接点时不太方便，因此关于于裁切的过的worlView 影像的处理，首先要进行同轨影像拼合处理。1、新建工程设置投影坐标系，

20、根据控制点的投影坐标来设置，与控制点一致。添加拼接任务，此处的处理设置与引入影像是相同的作用，故做了影像拼接就不必再常规流程中的“引入影像”一步。添加任务后，开始处理|1001 二右影像：|1002丁5、定向添加控制点，多条带的话需要在条带间添加连接点。该数据没有控制点，均匀添加了 连接点，进行定向修正 rpc 参数。beijing.rep - 记事本文件(F)耨(E)稻式(0)鶴助H)型接2 16 0 06髻蒼HE点连HE点苴 小S 苴 定 向 占总其籍查&单位权中误差 0.836 (像元) 像点坐标残差模型；21001 - 21002210020001-0. 534-0. 3200. 60

21、50. 220210020003-0. 249-0,1470. 2790.108210020004-0.562-0. 3366 6440.2342100200050.0270+016-0. 038-0. 0122100200066 5120. 904-0. 704-0. 628 |210020007-0.194-0.1180. 2130. 079定向控制点的均方根中误差和残差;工程中没有定向控制点工程中役有检查控制点6、模型裁切，worlview 射重采样”。立体像关于:21001 - 21002左爭橡：H：mih&rkgjcamplettarlVij ingMin iges右毙橡：H： mz

22、h4ngExajnpleWarlinvJQbijages EMHUL右片仿：核线采样方法:II右片仿射重采样口II确定取消7、设置裁切方法&裁切后就可以用模型进行测图，匹配EM等。(E)査看w工焙到库中T并且享-名称劃录新建却1 决*修改日期Images啲位置.Prouct,.K5OGO96O35.mr,: ProGeosys-2012/12/3 16:052012/12/3 16:052012/11/16 2012/11/16 文 怯 文磁MR 文 文住9、生产 EM EM，速度较慢。SM，是一个立体像关于匹配一整块另一种是模型匹配 EM，裁切模型后，关于小块的模型进行匹配EM，速度快。EM

23、匹配完后进行编辑，出EM成果。10、生产正射影像按默认纠正方法纠正即可。可以添加修改任务11、可将 EM 叠加正射影像察看地形是否匹配。用【生产工具】打开em，然后加载正射影像。EM】菜单下的【EM 三维显示六、TH（ Z3）数据案例1、新建工程设置投影坐标系，根据控制点的投影坐标来设置，与控制点一致。2、引入影像选择影像类型为天绘一号，输出格式为vz或tif均可I处理类盅转角莎向缩放丫方向型 号度放曲 - 刀 0F珂 EV 任务列夷处理消息1J J_J处理I退出、设置影像列表添加影像，根据引入影像时设置的输出影像类型来选择。影橡塑塑I 斑二号二 J输出格式|TIFFC*. tif) 211

24、r 厂ETC错误 機索|移除一切I、设置三线阵列表S xZ n)YiswwM-a rfnrxv rufMi m陛 i njuu科凹祸B 已旳屏 I-11TWL L_P2310?riwL.ij空niQM .THlL-n._P31ff?!口 曰 LJ MSP阳：叮令 U1H 14m H:rKHOLJ1耳】LQ：nj郴1闻|加乍帶EHM时心外芒ai細|LL悴3.TW2 |L?祥常？田兰创田兰创盘制 观忌界anwmcp 易刃葩曲a *、设置控制点交换5交换V/Z交换VZT販*1 添加X 删除 固编辑 分组EI 保存 钉退出6、疋向添加控制点，多条带的话需要在条带间添加连接点。该数据没有控制点，均匀添加了连接点，进行定向修正 rpc 参数。活=察看定向报告，像方小于半个像元，物方控制点的精度要求参照国家规矩。U&.rep -记!E）Q）直百凹卒助H）口區1工 PI煙1E:Testl 360ri ent at i onl 3 &Bpro 136称丄i.mknnwn称用占Joo0301向9查9