《遥感地学分析-图像解译(pci几何纠正+arcinfoworkstation目视解译+pci监督分类)【beta1.0】》由会员分享，可在线阅读，更多相关《遥感地学分析-图像解译(pci几何纠正+arcinfoworkstation目视解译+pci监督分类)【beta1.0】（41页珍藏版）》请在金锄头文库上搜索。
1、几何校正 图像解译 图像监督分类 附录By fhys（不少错漏，多指正）大致流程：几何校正几何粗 校正重投影裁剪部分区域几何精校正图像解译统一投影信息创建主控层创建解译图层添加解译属性字段解译各土地利用类型建立拓扑修改悬挂弧和 lab 点AML 批处理实用Arcmap 中成图ArcMap 中 coverage 文件与 shp 文件的转换图像监督分类附录相关概念coverage 格式 文 件概念如何在 win7 打开旧版 help 文件Arc 命令 &watchArcedit 命令 de【drawenvironment】Arc 命令 createArcedit 命令 sel【select】 Ar
2、cedit 命令putArc 命令 buildArc 命令 ab【label】Arc 命令 additem属性表的相关操作命令Arc 命令 CLEANCLEAN 与 build 区别AML 批处理实用nice 的命令一、 几何校正【Geometric correction】 本次所用软件：PCI Geomatica V9.0.0。 操作过程中有不懂的，可以查看每个操作框的“help”获取本机帮助，或者在 PCI 的安装目录的“manuals”文件夹找到对应的 PDF 帮助文档，如果我的PCI 是装在 “C:Geomatica_V90”下，那么对应的 PDF 帮助文档在“C:Geomatica_
3、V90manuals ”中，看帮助文档很直接，比在网上搜索教程详细很多。 这里把“配准”和“校正”认为是同一个概念，不做区分。 这里的“分辨率”指的是地物分辨率（如一个像元代表地物 30*30m 的分辨率就为 30m） ，而不是指图像的行列数的相乘 先列出原图与粗校正后的图像比较（左原，右后） ，便于后面清楚的步骤：#几何粗校正利用 PCI Geomatica 中的 GCPWorks 模块，进行几何粗校正。进入 GCPWorks 模块，在该菜单下，选择适当的选项，然后 Accept 继续:几何粗校正菜单 应选 备注Processing Requirements Full Processing
4、可选的校正过程（分别是）：1、 包括校正所有步骤2、 只是选择或查看 GCP（地面控制点）3、 只是拼接图像Mathematical Model Polynomial 配准方式：1、 多项式的方式配准2、 单纯的图像拼接？3、 薄板样条的方式配准Source of GCPs User Entered Coordinates地面控制点（GCP）的来源：1、 用一幅已配准的图像来对比选点2、 使用数字化仪（有点古老）3、 用矢量图中的点4、 用户自定义输入控制点（输入点的行列位置及其对应的投影坐标位置）5、 ？在下一个对话框中，点击 Select Uncorrected Image，选择待校正的图
5、像，即给定的 gz.pix，选定需要显示的三个波段（R5,G4,B3） ，合成后即为标准假彩色图像。接着，点击 Define Georeferencing Units，进入对话框，在原本显示“Meter”的下拉列表框中选择“Other” ，进入对话框，选择 TM-Transverse Hercator 横轴墨卡托投影。在 Earth Model 中，选择 TM-Transverse Hercator 横轴墨卡托投影方式。接着，设定投影参数，在 Longitude 中设定 True origin 的值为 111d0000.0000“E（TM 投影中，中央经线的位置，因为广州所跨经度约在 112-
6、113E 之间，其最接近的中央经线是 111E）,在 False 中设定 Easting 的值为 500000.000（在 TM 投影中，坐标原点西移 500km，为了使其 x、y 坐标均为正） ，如图：点击 Accept 进入下一个对话框，定义椭球类型，选择 E015-Krassovsky1940（E015 的 E 代表 Ellipsoids，椭球体的意思） ，如图：在投影信息显示框点击 Accept，确定设置，回到 GCPWorks 模块的对话框。接着，点击 Collect GCPs，进行控制点选择。由于是粗校正，所以只确定四个角点的大地坐标以及其相对应的行列号。方法可以是: 直接在主窗口
7、中输入 E、N 和对应的 P、L 然后点 Accept as GCP (用这个比较方便)也可以先写一个 txt 文件，然后一次性导入（Fileload GCP file） 具体是先把 TM 图像对应的头文件的四个边角点的投影的 X、Y 值按照容易理解的方式放在 Excel 表格中然后复制到 txt 文本中，txt 中的格式可以参见 PCI 图框右下角的 Help 获取最详细的帮助文档，这里用的是 IPLXY 格式（即序号 列数 行数 投影的X 值 投影的 Y 值，中间用空格或 tab 符隔开）接着导入到 PCI 中选择 IPLXY 格式，点 Apply Format点 Accept，返回到 G
8、CP 主窗口再点 Close，返回到 GCPWorks 主窗口最后是输出保存粗校正的图像 Perform Registration to Disk，但要先进行下面步骤：由于在 PCI Geomatica V9.0.0 中，需要新建一个含相同投影的（TM、E015、111E 中央经线、500km 西移） 、有足够图幅大小的图像来保存这个粗校正后的图像，才能保存粗校正后的图像（BigBigBug） 。所以打开 ImageWorks，随便选一个图像或者 User Defined，波段任意（主要是为了新建图像，不是使用那个打开的图） ，然后 accept & load（如果是User Defined
9、就点 Accept）接着，打开菜单栏 FileNew，进入新建文件框，选择 pix 格式，OK ，进入新建文件的详细设置（这里有些需要特别注意，最好按后面的顺序，便于按照这里进行说明）：7 个波段（因为被校正的图有 7 个波段）选择“Use bounds and resolution”（用边界【左上角和右上角】和分辨率（30*30m 代表一个像元）来决定行列数，所以不用设置行列数，它会随着下面的设置而自定改变）选择投影，跟前面一样（TM、E015、111E 中央经线、500km 西移）设置图像的左上角和右下角（因为只需要这两个点就可以确定图幅的大小） ，左上角对应 TM 图像对应头文件中 Xm
10、in、Ymax，右下角对应头文件的Xmax、Ymin到这步，Pixel Size 会变成其他数字，把它设为 30*30（这时会发现右下角的 X、Y 值会有轻微改变，这是因为像元是整个的，不能只有半个像元，为了刚好能把图像“切成”整数个像元，软件会自动设置 UL 和 LR 的与对应 E 和 N 的差为分辨率的倍数( 这里为 30)，原则是至少保留边角点在生成的图像中） 。到这里，可以看到图像的行列数已经变成比未校正的图像的行列数大一些（这是因为粗校正后，图像会按照对应投影的 X、Y 值变得倾斜，导致图幅要更大才【行列数更多】能容得下校正后的图像，所以不用改成与原图一样的行列数）最后点 Creat
11、e 生成新建文件新建文件完毕后，回到粗校正的输入文件的窗口中，选“Default” （自动按顺序对应）重采样（Resampling Mode）默认最近邻域插值（Nearest）Model 默认background 默认Memory（缓存） ，可以设成最大，可以少加快生成速度最后 Perform Registration 输出最终结果到文件中打开看看是否正确、#重投影 之所以要重投影是因为，精校正中的参考图（俗称“准确的图” ）的投影不是 TM E015 （横轴墨卡托投影【也即高斯克里格投影】 ，克拉索夫斯基1940 椭球体） ，而是 ACEA E015（Albers 投影【正轴等积割圆锥投影】
12、 ） ，要统一地理基准面和投影系统才能准确进行精校正。 其实也可以把“准确的图”转为 TM E015，只要统一就行。 在我国，一般 Albers 投影的中心经度是 105E（我国区域的中心大致在 115E） ，标准纬线是 25N 和 47N，原点经纬度（105E ，0） ，偏移量为 0【？】 。中国地图一般采用这个方式的投影，但在各省里，可以根据各自地理位置设置自定的参数（自定义中央经线、标准纬线、坐标原点） ，这里采用对广州地区范围的参数设置会更准确，但统一起见，直接采用全国地图的设置方式，也就是下面截图的设置。 这时的 Destination File 保存的文件可以直接指定，而不像之前那
13、样要新建了才能用来保存。 最后在 Reprojection 窗口中设置投影和分辨率时，步骤跟前面基本一样：选投影选择 Use bounds and resolution 方式 设分辨率为 30*30，提示更改时与前面一样，也是更改 change lower right 不用修改左上角点和右下角点，会自动随分辨率而变。全选到目标图层 destination layersselectallreproject 输出图像。 得到结果会发现丢失了部分图像数据，如下面的图，这大概是因为 PCI 只把原图像四个边角点（不是卫星图的轮廓边角点，而是原图像显示范围的四个边角点）转换投影后的坐标作为转投影后的图像
14、的界限，而当转投影后图像上下边界的一些地区超出了该范围后就变成丢失了。 但这里需要精校正的只是其中部分区域，所以丢失一点边界的信息也不影响效果。若一定要保留所有图像信息，可以把对应的 UL 和 LR 设置到合适的值。 另外，要查看一幅图像的坐标到底有没有像 TM（横轴墨卡托）投影那样西移了 500km 的话，要在上面的 Reprojection 中选 other，接着选一个投影，接着 accept 才会看到这幅图默认是否偏移了。#裁剪部分区域 之所以要裁剪只是因为“准确的图”只有一小部分，所以为了保证图像的精校正精度，最好是：“准确的图”的范围=“不准确的图”的范围。先看看效果（左为“准确的图
15、” ，右为对重投影后的大图“裁剪后的图” ）：步骤：依然在 Utility 的 tools 中选择 Subset，进入 Subset 主窗口输入输出文件的位置和名字（可以直接指定文件名，而不用新建文件才能保存）点 preview 预览，拉伸框框到合适区域（这个预览图看起来太小了）然后 select all，并 add 到右边去最后输出图像#几何精校正 类似于粗校正，也是在 GCPWorks 中，进入 GCPWorks 主窗口，在 Progressing Requirements 中选择 Full Processing，在 Mathematical Model 选择 Thin Plate Spline（薄板样条函数） 【也可以选择 Polynomial 多项式方式，各有优缺点】 ，accept 进入下一步。选择未精校正图像，这里是 gz_new1reprojSubset.pix，选 5 4 3 作为显示波段，load and close回到校正流程的主窗口，接着选择“准确的图”m4
