• ArcGIS中配准栅格图像并且矢量化(转载)


    栅格图像的校正和坐标系确定

    启动ArcMap,新建一个新工程,右键Layers选择Add Data…添加TIF图像,将出现如下提示(如果提示无法加载rester data时请安装ArcGIS Desktop SP3补丁),单击Yes确定,加载图像后提示图像没有进行配准,确定然后配准图像。

    图像加载后即可看到图像内容,右键工具栏打开Georeferencing工具条,进行图像的配准工作,在配准之前最好先保存工程。

    File菜单下打开Map Properties编辑地图属性,Data Source Options可设置保存地图文件的相对路径和绝对路径。(这里选择相对路径以确保将工程复制到其他机器可用)。

    配准前要先读懂地图,望都县土地利用现状图采用1954北京坐标系,比例尺1:40000,查阅河北省地图发现望都县位于东经115度附近,那么按6度分带属于20带中央经线117度,按3度分带属于38带。从图框看到的公里数发现没有带号,应该是公里数。 

    这里只找了4个点进行配置(可以找更多的点),从左到右从下到上,逆时针编号为1234;在ArcMap中单击Georefercning工具条上的Add Control Ponit工具(先掉Auto Adjuest选项),添加4个点控制点。

    然后编辑Link Table中的4个控制点的代表的公里数,然后单击“Georeferecning下拉菜单的Auto Adjuest”图像即进行校正这时可看到参差值这里是0.00175Total RMS)非常小说明配准较为精确。单击Save按钮可将控制点信息保存到文件,单击Load按钮可从文件加载控制点坐标。

    给校准后的地图选择适合的坐标系,右键Layers打开Properties对话框属性对话框选择投影坐标系,(Prokected Coordinate Systems)展开Predefined/ Prokected Coordinate Systems/Gauss Kruger/Beijing 1954下找Beijing 1954 GK Zone 20坐标系(高斯克里克投影20带无带号),单击确定保存工程;这时配准工作即完成,在状态栏就可以看到正确的坐标单位了。 

    最后保存校正重新生成采样数据,单击“Georeferencing”工具条的“Rectify”菜单矫正并重采样栅格生成新的栅格文件; 

    打开“Save As”对话框设置保存参数和保存的校正栅格图像格式、保存位置等信息;其中“Resample Type”指定重新采样的类型,有三种选项Nearest Neighborfor discreate data)自然邻近内插(不连续数据)、Bilinear Interpolationfor continuous data)双线性内插(连续数据)、Cubic Convolutionfor continuous data)立方卷积内插(连续数据);这里选择了第二项,确定生成新的栅格数据,然后加载到ArcMap中进行下一步矢量化工作。

     

    栅格图像的矢量化工作

    单击ArcMap工具条上的ArcCatalog按钮打开ArcCatalog程序(ArcGIS的地理信息资源都这里完成创建、删除、复制等管理工作),在Catalog树下展开工程所在位置,鼠标右键菜单中选择New子菜单的Shapefile…新建一个Shape格式的地理要素文件(地理要素可存储为其他格式);

    在“Create New Shapefile”对话框中给新的要素命名,在“Feature Type”要素类型下拉列表框中选择创建要素的类型(一个Shape文件只能表示一种要素),如:Ponit点、Polyline多边型线、Polygon多边形面和MultiPointMultiPatch,这里我们创建一个“Ponit点”文件;在“Spatial Reference”框中没有指定坐标系,因此单击“Edit…”按钮给新建的要素类指定坐标系;

    单击“Edit…”按钮后打开“Spatial Reference Properties”对话框,然后单击“Select…”按钮选择一个合适的投影坐标系,这里我们选择“Beijing 1954 GK Zone 20”坐标系,最后确定返回完成Shapefile的创建;

    回到ArcMap中添加新建的Shapefile点文件到Layers下;

    将栅格图对准到矢量图上,先将两个图层打开,并且将georeferencing工具条调出来(在标题栏空白处右键即可找到)。然后在矢量图上分别在上下左右中间等各找一个参照点(以易于分辨为准)并且纪录她们的坐标值。如果栅格和矢量的坐标系统不一样,或者两个的单位不一样,要将栅格的坐标转一下。使二者一致。一般要把参考点的经纬度的度分秒形式坐标改成带小数的度的坐标。这几个点坐标形式的转化可通过自己编写一个小程序来实现。

    找好参照点(至少要五个)之后,在georeferencing中的layer选择栅格图层,点击add control point图标(一个一头红一头绿的双箭头图标)在栅格图层中找到你的参照点,单击之后再击右键选择inputX andY,在弹出的对话框中输入刚刚采集到的点的坐标即可,依次将各个参照点全部输入。即完成了配准的工作。

     

    配准前的准备工作:

        要配准的栅格和矢量图层必须具有相同的投应坐标系。如果不相同,要修改其中一个的投影坐标系。修改(添加)投影要在arccatlog里面进行。在目录中找到文件点击右键打开属性表,在Spatial Reference中点击编辑,打开对话框,点import导入参考图层即可加载参考图层的投影坐标系统。应用并确定即可。

    1 图幅拼接针对ARCGIS 9.0
    a)         保证待拼接图层的属性项名称和类型都相同,不相同则修改(修改方法见属性表操作);
    b)        本操作假定图层已做过接边,如没有接边请先在arcedit中使用edgematch工具接边;
    c)        ArcToolbox>>Data Management ToolsGeneral中找到Append工具,以待拼接图层的任意一个图层作为目标图层,把其它图层加入这个图层,Schema Type (Optional)选择TEST,然后点击OK按钮;
    d)        进入ArcInfo Workstation,在Arc提示符(Arc:)后依次输入&wo …aedisp 9999de allec …drawef arcsel allunsplit none命令消除伪结点;
    e)         这样目标图层就是拼接好的图层(但是仍然可能存在未消除的伪结点,有打断的线实体,可多使用unsplit none命令几次);
    2 图幅裁剪针对ARCGIS 9.0
    a)         ArcTools中选择Command Tools点击OK,出现Command Tools菜单;
    b)        在菜单中选择Analysis>>Map SheetsClip工具;
    c)        依次填写Clip featureInput coverClip coverOutput cover等参数,点击OK按钮即完成
    注意:
    如使用ArcToolbox中的Clip工具则ARCGIS有可能会自动改变属性某些属性项的类型。
    3 属性表操作针对ARCGIS 9.0
    ArcCatalogArcMapArcInfo Workstation中都可以对图层属性表进行操作,推荐使用ArcInfo Workstation
    a)         进入ArcInfo Workstation,在Arc提示符(Arc:)后输入tables进入属性表操作工具(提示符变成Tables:);
    b)        使用select …(如select terlk.aat)命令选择要编辑的属性表;
    c)        items …命令显示所有的属性字段定义,additem …增加属性项,dropitem …删除属性项,alter改变属性项;
    d)        cal …命令可以录入或修改属性项值;

  • 相关阅读:
    oracle 查看表空间使用率
    解决linux下vim中文乱码问题
    linux 时间同步
    oracle ho与mysql system命令
    mysql 重置root密码
    2020 10 26
    2020 10 24
    2020 10 23
    2020 10 22
    2020 10 21
  • 原文地址:https://www.cnblogs.com/zhangjun1130/p/1637400.html
Copyright © 2020-2023  润新知