湿地调查-CBERS影像几何校正流程规范.docVIP

CBERS影像几何校正流程规范 影像准备 数据源 以2008年7-9月份中巴地球资源卫星(CBERS-02B和CBERS-02)CCD遥感影像为主，缺少部分由国产环境减灾小卫星星座A星和B星CCD影像补充，或由Landsat TM、ETM数据等补充。 数据质量检查 影像各波段之间配准无问题(小于0.5个像元)，影像不存在明显色差，影像数据不缺少波段，影像无明显条带噪声和掉线情况。 单景影像平均云量小于10%，湿地上空无云覆盖；受人为干扰影响比较小的不易发生变化的区域，可适当放宽；云量没有明显降低影像质量，影响后续的影像解译判读。 影像合成 利用CBERS-CCD影像的前四个波段进行合成(第5波段是全色波段，不予考虑)，通过ERDAS软件实现(软件版本8.7以上)。具体流程如下： 解压CBERS-CCD的2级压缩包文件，得到5个Geotiff文件和一个xml头文件(见图1.1)； 图1.1 CBERS数据解压缩 双击打开xml头文件，读取影像的成像时间，记录下来，见图1.2； 图1.2 从CBERS头文件中读取影像成像时间 打开ERDAS，点击Interpreter模块，再依次点击Utility、Layer Stack，然后点击打开文件夹图标(见图1.3)，从数据存放路径中找到相应的解压缩后CBERS影像： 图1.3 进入LayerStack对话框 选择文件类型(Files of type)，改为TIFF，选择第一个文件“369-49-L20000157653-1.TIF”,然后点击OK。 图1.4 选择波段数据 点击Add加入1波段“369-49-L20000157653-1.TIF”数据，重复图1.4步骤，依次顺序加入波段2、3、4的tiff文件；勾上“Ignore Zero in Stats”选项，然后在Output File(*.img)输出框中选择文件存放目录，给波段合成文件命名，最后点击OK(见图1.5)。 附文件命名原则：影像所在区域名(拼音)-卫星名-传感器名-行号-列号-成像时间，如示例文件为CBERS-02B星2008年9月6日成像影像，影像所在区域为白城，因此命名为“baicheng-cbers02b-ccd-369-49-080906”，影像成像时间的获得见图1.2。 图1.5 波段合成 几何校正 在ERDAS软件下，以1：5万国家标准分幅地形图(西安80坐标系)为参考图像，逐一对每景CBERS-CCD影像进行几何精校正，以满足1：5万到1：10万的制图精度。 启动ERDAS，点击面板上Viewer模块，在View1中打开需要校正的遥感影像，再次点击Viewer模块，在View2中打开参考地形图(见图2.1)； 图2.1 打开影像和地形图 在打开cbers影像的View1选择Raster菜单下Geometric Correction并单击，弹出Set Geometric Model对话框(见图2.2)； 图2.2 进入Geomtric Correction 在对话框Set Geometric Model中选择多项式模型Polynominal，然后点击OK(见图2.3)； 图2.3 选择几何校正模型 不改动Polynominal Order的次数1(为了后续选点方便，控制点选足6个后，在把多项式次数改为2)，点击Close (见图2.4)； 图2.4 选择多项式次数 在出现对话框中点击OK(见图2.5)； 图2.5 选择控制点来源 用鼠标点击打开地形图的View2作为参考影像 (见图2.6)； 图2.6 选择控制点来源 在弹出对话框里点击OK，进入选择控制点界面 (见图2.7)； 图2.7 选择控制点界面 在View2菜单的Utility下选择Selector Properties，在弹出的对话框中选择颜色，在GCP Tool内Color下方空白格内单击右键，选择颜色 (见图2.8)； 图2.8 设置放大框颜色和控制点颜色 现在开始选择同名控制点 (Ground Control Point, GCP)，在影像上和地形图上选择对应点。点击上方Geo Correction Tools中圆形中间带十字图标，在View1中选择控制点，再次点击十字图标，在View2中选择对应的控制点 (见图2.9)。 附控制点选择原则：a) 由于影像的成像时间与地形图成图时间相差在20年以上，控制点要以不易随时间变化的标志物为主(人造物体)，例如：桥梁、大型道路交叉口、人工沟渠和道路交叉口、人工堤坝交叉口等，不建议选取水体、农田、城市、居民点等易随时间变化的地理标志；b) 控制点分布要均匀，保证影像的四角和四边都有点分布；c) 一景影像的控制点在15-20个左右，保证均方根误差小于20米，控制在CCD影像一个像元以内。 图2.9