矢量地图坐标转换方法及工具.docVIP

矢量地图坐标转换方法及工具 在MapInfo中，通常采用大地坐标的北京54标准。与GPS采用的WGS84在同一点上相差经纬度相差有10－20分（仅仅从地图在MAPINFO的显示数据来说没涉及到NumericCoordSys的坐标系统）。虽然现有GIS平台中都预定义有上百个基准面供用户选用，但均没有我们国家的基准面定义。假如精度要求不高，可利用前苏联的Pulkovo 1942基准面(Mapinfo中代号为1001)代替北京54坐标系；假如精度要求较高，如土地利用、海域使用、城市基建等GIS系统，则需要自定义基准面。实际工作中一般都根据工作区内已知的北京54坐标控制点计算转换参数，如果工作区内有足够多的已知北京54与WGS84坐标控制点，可直接计算坐标转换的7参数或3参数；当工作区内有3个已知北京54与WGS84坐标控制点时，可用下式计算WGS84到北京54坐标的转换参数(A、B、C、D、E、F)：x54 = AX84 + BY84 + C，y54 = DX84 + EY84 + F，多余一点用作检验；在只有一个已知控制点的情况下(往往如此)，用已知点的北京54与WGS84坐标之差作为平移参数，当工作区范围不大时精度也足够了。 一般使用3参数，将3个控制点的数据代入方程，求解出参数，可以利用Excel求解，如图 在黄色区域分别输入3个点的不同坐标系的数据，通过消元法，在红色区域会得出a、b、c、d、e、f参数的值，可以在下图的黄色区域输入数据进行验证，检验误差 为了实现地图数据的转换，我用MapBasic编写了一个程序。该程序在MapInfo程序界面中添加了一个转换的菜单，主要包括2部分：一是参数设置，用于设置并保存用于转换的3参数；二是执行转换，选择需要进行转换的表并转换。如图： 参数设置内容如图： Ｄｅｆａｕｌｔ按钮可以将３参数恢复到缺省值，点击ＯＫ按钮会将输入的参数设置进系统，并保存为config.ini文件。 转换内容如图： 可以选择需要转换的表，如果点击全选则自动全部选中所有表，再次点击则取消全部选中。 该程序具有记录操作的功能，将主要的操作记录到ｌｏｇ文件中。 附件：源程序 ***** 坐标转换 ***** Created by 孙国文,江苏电信徐州分公司无线中心 ***** Jan 28, 2005 include mapbasic.def declare sub main declare sub RCV declare sub CB declare sub HB declare sub editlayer declare sub minmap declare sub restoremap declare sub checktab declare sub objmove declare sub deltab declare sub YesButton declare sub NoButton declare sub EndPro declare sub About declare sub provision declare sub objmove_main declare sub getconffromfile declare sub setdefaultconfig declare Sub getres declare Sub changestatemulti declare Sub changestatecheck dim r,v,i,j,h,poly,pnode,cnt as integer dim z as Smallint dim numtab,tabtemp,tabnum,numlayer,objoff,numwin,mx,my,arccnt as integer dim ptcnt,plcnt,lcnt,regcnt,filecnt,rctcnt,rrtcnt,elpcnt,txcnt,ewindow as integer dim sobj as object dim nx,ny,xoff,yoff,go1,go2,go3,go4,tsp,ta,tlx,tly,gx1,gy1,gx2,gy2 as float dim tabstr(1) as string dim tabname,runstr,ts,elayer,ct as string dim tabobj as alias dim layerfound,cancel as logical dim filenum(1) as integer dim tf as font dim tj,tar as smallint dim tjs as string 存储参数 dim da,db,dc,dd,de,df as float dim a,b,c,d,e