PCI遥感实验报告

1、PCI遥感实验报告地信1101 XXX2013年6月18日实验二几何校正准备数据：125-41.pix，125-42.pix （均包含7个波段，见实验一）靖县地形图.jpg。实验步骤：第一步、熟悉图像几何校正、镶嵌的基本原理第二步、图像几何校正运行 PCI，选择 GCPWorks 模块，在 Source of GCPs 选择 User EnteredCoordinates （用户输入投影坐标系统），点击Accept后，弹出校正模块:OptionsAccept =选择第一项加载需要校正的图像:FrocesiTig Requiraments* Full FrocesingCollect/rfivi

2、 e* GCPs only Hosai c Onlyhl鱼al Model* F alynomi alNoti4 训osHckin电 wd Dr thoEn 菖ine)C Thin Plate Spline Satellite Ortho Correct!onSource GCPsf (Jftoc&dfid Image H&rAcopy Hap on Bigi tiling T&bl C VectarsSelect Uncorrecttd Image Define Getreferencing Units Collaet GCFeSeiect Mos&ic AreaPer form Regis

3、trtion to Disk柯 GCPWorks SetupUssr jut ered Cacrdinates；|P Chip DtabiseCorrespanding Processing Steps厂ffi PCI GCPWor.,.T| Select Vikcorrtciftd Tin包官e | Dwfiri电 Gsor&ftr tncir(g Vni ts _| Collect GCFs | Select c Area | Frform Regi stratiort to Bi sk选择第二项，选择Other确定投影系统，点击Earth Model确定地球模型r柯 PCI GCPWor

4、Select Unccrrecta ImageFile OptionsBefiGeorUni tsDIo+-Genw 1 c Projeeti-orLS User PrsjetiLCCMCMERMSC Other GearferfKe UnitsROBRSOSGSGDOSIHSOM- Lambert Conform*丄 Conic-Miller Cylindrical-Mercator-Modif ied Stereographic Confor*al-Or thograph i c- Oblique Mercator (2 Points)- OhLique Mercaten? (Point

5、and Azimuth)- Folyconic- Polar- Stereographic-l?oh inson- R曰匚tified Skew Orthonorphic (Point and Azimuth)- S t ereograph ic- Stereographic Double-Sinusoidal-Space Oblique Mercat口rTran averse Iferu鼻 t or (Gmgs-Xr ueqmrVDG Van der GrintenCancelHelp | i电 Projection Definition注意输入6度带的中央经度与向东平移 500公里（500

6、000米）Earth ModelsDatws Ellipsoi dsEOEDEOEOEOEOEOEOEOOOOOOOOOOOE010EO11EO12EO13E014Clarke 1866*Clarke 1880Bessel 1S41Ne International 1967International 1S24WQS 72Everest (India 1830)WGS 66GRS 1980Airy 1830Everest (U. Malaysia and Singapore 194S)Modified AiryUQS 84Modifled FischeT I960Australian Natia

7、nal 196SKrassovsky 1940E015Hough i960T7* j 1- r 门产AcceptCancelHelp |选择第三项采集地面控制点。在采集地面控制点之前，利用看图软件打开扫 描后的电子地形图“靖县地形图jpg。分别在遥感图像和地形图中找到一个同名点， 然后在地形图中量算出该点的 大地坐标，精确到米，X坐标为6位（要去掉2位6度带的带号），丫坐标7位 再将坐标输入到GCP编辑窗口中，并点击Accept as GCP接受为一个控制点。 用同样的方法采集更多的地面控制点。 注意：前三个点不显示误差，从第四个点 开始才出现误差。一般要求选择15个以上控制点，并且各点的误差

8、1个像元， 将误差过大的点删除，直到误差1个像元为止（实验时误差可放宽到2 - 3个像元）。N| Define Georeferencing Units| Collect GCPs| Select Mosaic Area_| Perform Registrytion to Disktg PCI GCPWor.TSHelpFile Options_| Select Uncorrected ImageS432109876S52.;2 Ar7:0 69(-0,30,-0.62)(363000.0013 ;0.92(-0.38.0.84)(360000 OOi9 :0.52(0.27,0.45)(35

9、7000 00110 :0.93(0.0 66)(356000.00i11:0.80(-0,79,0,13)(359000.00112 :0.4S(0.32,-0.35)(367000.00i匸13 :0.74(0.22”-0.70)(36BOOO.00114 :1 14(7 93)(3S3000.00i15 ;0.37(0.08,-0.36)(373000 OQi*rpr1最后，关闭GCP编辑窗口，选择第五项（如果对整个图像进行几何校正，则第四项可以省略），执行几何校正闻 New File Type Selector点击OKGeo-Referencing Information 选择 Use

10、 bands and resolution。并点击 Create与Yes生产校正后的图像选择 Change ResolutionQu 亡 stiona x FIe弓e reviev ill iriformatiotk to ensure that it is correct. The fils x* J will be up to about 0. 79 Mbytes in si ze on 鱼 di sk with B3224. 66 MBytes apparently free.Dq you wish to crthi z匸直勺皿4 |选择Yes电 Disk to Disk Registr

11、ation回到图像校正窗口，建立原图像与校正后图像各通道的对应关系， 设置重采样方式（一般选择Cubic）与内存大少（一般300左右），点击Perform Registration执行校正下图为具有坐标系的靖县地形图:EC站县地形ImageryRGBjGeocode I|369731?258E 2941062. 622N242:1235:2212:3Whies AU 0| All 0n|时nrVectorsAU 0| All 0n|EnhancementsZoomQIQP二1He(ppdI369T3L.258E 2941062. 622N242:1235.2 212:3 |R: 1 G.2 B

12、:3PQ InugeWoriu V8.ZO (VDO)洁县瞬固pbcFil Edit yiew Icds ik&jify轲旨團0#|谓|叫旬|国珈目|/|/|好I ?F面我们利用加了地理坐标的地形图来对 125-42遥感影像进行几何校正:利用 GCPWorksGeomatica Toolbar V8.2.C51 GCPWorks SetupProcessing RequirementsG Full Processing Collect/review GCPs only 广 Mosaic OnlyMathematical ModelG PolynomialNone (Mosaicking and

13、 OrthoEngine) 厂 Thin Platc SplineSaLt ell it Ortho CorrectionSource of GCPsG Geocoded ImageHardcopy Map on Dig it izing Table ygctors User Entered Coordinates Chip DatabaseCorresponding Processing StepsSelect Uncorrected ImageSelect Georeferenced DatabaseCollect GCPsSelect Mosaic AreaColour Matching

14、Preregi strati on CheckingPerform Registration to DiskAccept ICancelHelp |S PCI GCPWork.HelpEil 皀 Optionsj Select Uncorrected ImageI Select Georeferenced Database| Collect GCPs_ Select Mosai c AreaColour Matching| Pre-registrati on Checking| Perform Registration to DiskFi e: C:Geomatica_V82user125-4

15、2.pixDatabase Channels12345678u 8u 8u 8u 8u8u8ul1imagewo:Exported imagevo:Transferred imagewo:Transferred imagewo:Transferred imagewo:Transferred imagewo:Transferred imagewo:Transferredfrom 1 from from from from from fromon111V82use11T R- _2j：|2 jJ：|3Car |巫劉Average ImageHelp |S PCI GCPWork.Voffneart

16、$j Uncorrected Image Overview控制点尽量多!tg GCP Selection an Ifl，在弹出的对话框中设置相关参数:o=III Iage TransferDEIC2打开125-4142.pix后，我们可以看到结果，如下所示Ac11AvailableStatus Run RunbgckClose HelpParameters to SetBEI-tt O； 0; 7844,_ ri17DEOir iooqT o,依次单击Status、Run即可。FILO.pixDEOCi PCI IaaceWorks V7.0.0: (TDO) 125-4142. ph点击GC

17、PWorks模块，默认设置，点击Accept确定，先后选择并加载未校正（需要镶嵌的图像）和参考图像，得到两景整图及其放大窗口图像共四个窗口 再选择Collect GCPs，分别在两景图像中初步选择重叠区域内的同名点， 并在放大窗口中精确定位。呵 GCPWorks SetupProcessing RequirementsG Full Processing Collect/review GCPs only 厂 Mosaic OnlyMathematical ModelG Polynomial None (Mosaicking and OrthoEngine) 厂 Thin Plat Spline

18、Satell1t Ortho CorrectionSource of GCPsW Geocoded ImageHardcopy Map on Digitizing Table( Vectors User Entered Coordinatesr Chip DatabaseCorresponding Processing Steps Select Uncorrected Image Select Georeferenced Database Collect GCPsSelect Mosaic AreaColour MatchingPre*regi strati on Checking Perfo

19、rm Registration to DiskAc ceptCancelHelp |5； PCI GCPWork.on11111Database ChannelsihESwTjl2 |：|23 I：目12345678u 8u 8u 8u 8u8u8ul1imagewo:Exported imagevo:Transferred imagewo:Transferred imagewo:Transferred imagewo:Transferred imagewo:Transferred imagewo:Transferredfrom 1 from from from from from fromF

20、:125-41pl25r41_5tl98125-41pl25r41_5t 125-41pl25r41_5t 125-41pl25r41_5t 125-41pl25r41_5t 125-41pl25r41_5t 125-41pl25r41_5tDatabase File SelectionLook in: | 詁 user二I NameDate modified *I 1CG2013/6/1215:KEi. SG2013/6/1215:|125.4i.pix2013/6/121%5!叵 125-41 42pix2013/6/18 22:11l125-42.pix2013/6/1216:0； |n

21、rFile name:|l2541.pix| Open |Ries of type: |AII Files CJ |CancelHelp|noose Directory | Choose Other |Average ImageClear | 瞌迫Load & CloseLoad |CloseHelp |PaGCPWork.电E PCI GCPWork.再选择Collect GCPs，分别在两景图像中初步选择重叠区域内的同名点,并在放大窗口中精确定位tg pa GCPWork.l3M 2*i X 、-JiwSPC15j GCP Selection and Editingle Reports E

22、levationGCP EditingF IAccept GCPUs*Check PtICkck GCPe： 0HJp|Modellit f* End3rd C 4th C* !EMS Error： 0.7Z0 60Accted GCFs： 10GCP 0rdring： W&rst |P* By IdId1 :Err口:r (Residual K.Y)(GeorefX.U心(-0.SO,0 27)(4191.125,2:0.70 (-0,42,-056)(4222.62S,3;0.61 (0.05f-061)(4256 G75,4:0.66 (0.66,0.0B)(4274 125,5 :1.

23、10 (1.09,-0.13)(435?.5G3,6：0.29 (-0,26, -009)(4409.500,7 :1.03 (-0.51,0.90)(4444 500,3 :0.59 (H -D.47)(3125.500.9 :1.23 (-0,87,0 86 )(339S.500,10 :0.S1 (0.50,-0.09)(3497 500,* 1用控制点尽量多选C pa fiCPWork.lWFile OptionsHelpSelect Vncorrecttd Tm幅圭Selectenced Databas&Colltct CCFiSelect Mos&i c AraCal our la

24、tchingFre-registration CheckingFerfcrm &egi strati on to Di sk+ lEtf* 2nd C* .f 5 thVertexEdi tingI2FNo.1069.49X I1T6.013AddMoveDelete |Insertini eK、：NwAcceptftd Gtoreferineed Vertices： 22l1:18:IT： 低 19： 20：6172,696S,7357,T962;84HBSOG,9024025691763712T65622421:18fi176J 22：】现15E沁p|呵 Colour,.点击MatchRe

25、gistr凰tiun Opti onsJfeiui siuiaiuuResampling Mode： r Nearest 厂 Bilinear i* Slibicj1- 8pti 16ptModel Ordr 1 1st 1 2nd 厂 Jz_d L 4th C* thBackground：Colour Matching： On OffBlend Iidth: (0?0CloseRegi str atioik Overvi ?aH*lp|选择Cubic,点击Registration Overview 查看预览效果点击 Un corrected Disk Registratio n 下的第 1、

26、2、3 行，再点击 Output FileImage Channels下的第1、2、3行。直到出现上面红框处的显示。如果点击错误，点击Reselect (重选)进行重新设定，Memory (MB)处改为500左右。点击Perform，出现下面的窗口。利用Image Works打开前三个波段镶嵌的图像如下:Help.fFite Edit View Tools Classify宝旦回也邑色刽回鱼u口 PCI ImageWorld V8.2.0： DO； 125-41-4.T |Fu11ks7 Window | 70:13230:3 |R：rG:2_B:3F面我们继续镶嵌第4、5、6波段，重复上面的

27、步骤，打开 GCPWorks选择125-41.pix文件点击第二项回73 mZl SKiia亠鼻砂i li* *上一火叱怎 Ungnvrtv SaWHofi电 GMiAlrmjnJ IrnjBjA 45-56-6|对0 Procrp Mcriitor,（_）| Performing Disk to Disk Registration Help25 65 B:6FOe Edit View Tools ClassifyBQ PCI ImageWprks V8.2.0: (VDO) 12 5Im agsry71173P 257L | 70:478:5 133:70:473:5133:6Enkance

28、iment s/|/l ：o：H-lIgliul old-lGraphics aii Off I All 血dJ Cursor|Di splay1T3FVectorsAll OffAllIIII1帀至此，七个波段中的六个已经镶嵌完成!第七波段部分镶嵌过程及结果如下:1234568u8u8u8u8u8ulFi e: C:GeomaticaDatabase Channelsimagewo:Exported imagevo:Transferred imagewo:Transferred imagewo:Transferred imagewo:Transferred imagewo:Transferr

29、edfrom 1 on F:125-41pl25r41_5tl98125-41pl25r41_5t125-41pl25r41_5t125-41pl25r41_5t125-41pl25r41_5t7 8u : iwo: Trainsf erredfrom 1 on F from 1 on F from 1 on F from 1 on F from 1 on Ff rom 1 on F:只剩下第七波段未镶嵌Car | DefaultAverage ImageLoad |Help |Load & CloseClse匚7H*y,回Linear P包鲫t$j Georeferencid Selecti

30、ont$j Uncorrected Image Oven/iewLinear | 1 ylloff J kB J BriPCl呵 Uncorrected Selection GCP Selection and EditinFile Reports ElevationGCP EditingAccept as GCPUse as Check PtDeleteModele 1st r 2nd C 3rd C 4th C !RMS Error： 0.510.3910.22(-0.14,-0.16)(4193.438,20.31(-0.30,-0.05)(4224.313,30.45(0.10,一043

31、)(4276.563,40.24(0.23,0.08)(2988.438,50.53(-0.44,0.29)(3513.438,60.36(-0.27,0.23)(2445.438,70.70(0.59,0.38)(3127.469,80.70(-0.66,-0.24)(-055)(4524.438,90.57(0.13z2618.313,100.89(0.76,0.47)(4918.688,4nrAccepted GCPs： 10GCP Ordering：( Worst By IdIdError (ResidualY)(GeorefCheck GCPs： 0Close|HJp|Colour.XDisk to Disk RegistrationUncorrected Image Channels5 8u: imagewo：Transferred from 1 on F:125-41pl25k41_5t6 Out put File Image Channels u: gepwork:Regi strati on Channel=5? file