用RM软件推算方格网点的经纬度坐标.doc

用RM软件推算方格网点的经纬度坐标 在进行小流域规划的时候，有些时候提供的扫描图上面的四个角点并没有提供经纬度，如果此时要求我们对这幅图进行投影变换的话，需要计算一些方格网点的经纬度。因为RM软件在进行投影变换的时候需要一些TIC点，而这些点的坐标必须是经纬度的。下面就简要的说一下，在RM软件中如何根据方格网点的坐标推算经纬度坐标。 如果底图上面有公里网格，根据这些公里网格上面标注的数字可以推算其采用的投影方式，对于小流域常见的地图来说，我国常采用高斯克吕格投影，一般在图上应该能看到其投影代号。比如图上坐标为21463km，这样可以看出这幅图位于高斯克吕格投影的21带，对于具体的小流域，我们知道其所在地的大致经纬度，这样就可以判断是采用的3度带还是6度带。从而确定其中央经线。比如对于上面的那个坐标，我们知道该地区的经度大致为东经122度。这样就可以推断其采用的是高斯克吕格投影，位于六度带，中央经线采用123度。从21463km中的4可以推断，其偏移距离应该为500000m。有的坐标会是“带号4实际坐标值”、“带号5实际坐标值”、“带号6实际坐标值”即21463000m、21563000m或者21663000m。通过4、5、6这些数都可以推断其坐标东移距离为500000m。我们这次的底图上面只有公里网格，上面的横坐标标有93、94、95等数字，由于底图是在MAPGIS下面做的，开始怀疑这些坐标标的有问题，后来到宏图打印社察看，这些坐标前面少标了一个6，其正确的x坐标应该是693000m、694000m、695000m等，而Y坐标同样少标注了52，其真实坐标应该是52930005304000m。由于图上的公里方格网是没有问题的，因此，可以得出每一个方格网节点的x、y坐标。经过推 算，结合本地区的大致经纬度，加上他们提供的一些信息，确定该底图采用的投影为：高斯克吕格投影，中央经线为120度，椭球体采用1980西安坐标系。1 知道各网格节点的坐标后，可以读出这些点的坐标将其输入到EXCEL表格里面，如下图2 输入完这些坐标后，将其另存为DBFIV格式，便于RM软件调用3 下面就开始将这些点的坐标转换成经纬度坐标：首先先新建一个工程，在工程里面新建或随便导入一个图层，打开其属性表4 要是不打开一个属性表的话，上面“数据表”的菜单不会显示出来。5 数据表/获取外部数据/链接表，出现下面的对话框6 数据类型选择DBASE IV，单击“下一步”7 在弹出的对话框里面单击“浏览”，选择存放坐标值的文件8 这样就可以把一个外部表链接到RM软件中了，由于表中还有一些数据的XY坐标值为0，我们下面的操作只是对于坐标值不为0的记录进行操作，因此需要将其筛选出来。9 这样就将所要记录筛选出来了。10 数据表/生成点文件11 设置横坐标字段和纵坐标字段，以及坐标单位和点图层名称后，就可以生成这些点了。如下图所示12 打开坐标点图层的属性对话框，更该其投影信息13 设置如下投影参数14 更改图层比例，改为1：1000015 确定后，保存坐标点文件，重新打开坐标点文件，可以设置经纬度显示坐标值16 这样就可以查看每一个点的经纬度了17 新建x1和x2字段，用来存放经纬度坐标，如果现在对字段进行负值的话，得到的还是平面坐标而不是经纬度坐标，如下图所示18 为了得到这些点的经纬度坐标，并把他们自动填写到x1和y1字段中间，必须将其投影转换成经纬度投影方式，由于图层现在已经有投影信息了，因此可以依据投影转换公式很容易对其进行坐标转换19 设置经纬度投影，设置椭球体后，确定就可以进行转化了20 转换后的图层如下图所示21 打开投影变化后的图层的属性表22 确定后，就计算出了每一个点的经度，结果如下图，同样的方法可以计算其Y坐标这样就把图上方格点的经纬度坐标计算出来了，这样对扫描的栅格图像进行投影转换的时候，编辑TIC点的时候，需要用那个点的经纬度，就可以依据其XY坐标在表里面查询对应的经纬度坐标。说明：如果只是要使栅格图层有坐标的话，其实可以依据公里网格的坐标值，对其进行几何校正，几何校正后的图的坐标值也是完全正确的坐标值，几何校正后的图层并没有投影信息，因此打开几何校正后的图，在图层属性里面给其填上对应的投影信息，这样得到的这幅图和把扫描图进行投影变换后的效果是一样的。由于本小流域规划的时候，开始没有进行扫描图的投影变换