PCI （PCI Geomatica）裸数据配准 与 ArcGIS配准比较

这里我们把"纠正"与"配准"当成同一个东西.

PCI对一幅裸数据图像进行配准的话，按照一般教程的步骤是：

为裸数据图像选择投影

-->根据裸数据的头文件的内容输入GCP（地面控制点）（这里一般为四个边角点）及其对应大地坐标

-->然后保存为一个新的文件（PCI9.0.0貌似要另外新建给它来保存，而且新建的文件要有之前选择的投影），这个过程成为粗纠正

-->然后是精纠正（配准），步骤是选择相同地物点，多选几个满足对应配准的需求数

-->最后保存最终的图像，完成。

  

这里我们会发现，在精纠正完了之后，我们在粗纠正里面输入的四个边角点的坐标已经发生变化了（不管多少，反正就变了），那么我们之前做的那步粗纠正是不是就没用了？不是的，它可以为我们选定大致的范围，便于我们进行下一步的精纠正。

这个说法貌似很合理，但是不然，这是用了ArcGIS中的配准后才想到的，粗纠正在ArcGIS里基本没太大必要，主要是因为输入裸数据的头文件的那四个边角点比较麻烦，而且容易出错，既然能在不进行粗纠正的前提下找到相同的地物进行配准的话，就可以直接这样做而略过输入四个边角点的步骤。

以TM图像来说，它的头文件（见文章最后附录）已经指定了椭球体，投影方式，也就是说，我们的粗纠正就是把原图（裸数据）的四个边角点与其投影平面上的坐标（x，y）分别对应上，即,一种平面与平面间的对应关系,需要插值,不涉及到与空间坐标的对应;同样的,精纠正也是这样子,只是精纠正的地面控制点是仔细来选,多选几个,那么也就没必要之前的粗纠正了,唯一的区别是,粗纠正后会保存文件一次,也就产生了插值,与直接就进行精纠正比较,可能会因为这个"插值"而导致图像最终配准结果的微小差别.TM的头文件给定的椭球体和投影信息并不是最准确的(否则也不会需要我们去精纠正,直接粗纠正就OK了),既然连TM官方放出的头文件数据也是不够准确的,那么被我们用于精纠正的参考图是怎么确定它自己本身的正确性的呢?猜想是用GPS对某些明显地物进行定位,然后匹配到地图上.

  

这里顺便提一个在ArcGIS里面比较简单准确的配准方法，我们知道ArcGIS里可以在Effect工具里面设置图层的透明度，而且ArcGIS中Georeferencing工具的配准是实时性的，多加一个点图像也跟着变，所以利用这两个工具进行配准会快速很多，效果又直观。师兄nice啊！

 

TM某幅图像的头文件

PRODUCT =99094126-01

WRS =122/04400

ACQUISITION DATE =19981222

SATELLITE =L5

INSTRUMENT =TM10

PRODUCT TYPE =ORBIT ORIENTED

PRODUCT SIZE =FULL SCENE

TYPE OF GEODETIC PROCESSING =SYSTEMATIC

RESAMPLING =D6

RAD GAINS/BIASES =

    00.00398/-0.0100

    00.00964/-0.0232

    00.00540/-0.0078

    00.01043/-0.0193

    00.00235/-0.0080

    00.00000/0.00000

    00.00154/-0.0040

TAPE SPANNING FLAG=1/1

START LINE #= 1

LINES PER VOL=40096

ORIENTATION = -7.90

PROJECTION =TM USGS PROJECTION # = 9 USGS

MAP ZONE = 100 USGS

PROJECTION PARAMETERS =

0.637824500000000D+07 0.669342748981920D-02 0.100000000000000D+01 0.000000000000000D+00 0.111000000009735D+09

0.000000000000000D+00 0.500000000000000D+06 0.000000000000000D+00 0.000000000000000D+00 0.000000000000000D+00

0.000000000000000D+00 0.000000000000000D+00 0.000000000000000D+00 0.000000000000000D+00 0.000000000000000D+00

EARTH ELLIPSOID =KRASOVSKY 1940

SEMI-MAJOR AXIS =6378245.000

SEMI-MINOR AXIS =6356863.000

PIXEL SIZE =30.00

PIXELS PER LINE= 6920

LINES PER IMAGE= 5728

UL     1125416.5166E 235948.8775N

    693826.876 2656304.236

UR     1143956.9507E 234422.2642N

    873930.222 2631305.562

LR     1141734.6921E 221334.0455N

    839637.283 2462578.052

LL     1123303.2182E 222839.5917N     659623.082 2487564.352

BANDS PRESENT =1234567

BLOCKING FACTOR = 1

RECORD LENGTH = 6920

SUN ELEVATION =37 SUN AZIMUTH =149

CENTER     1133556.9270E 230647.7510N

    766316.067 2559498.887

    3338 2865

OFFSET= 168 REVB