PDA环境下GPS信号的接收和处理方法

　　Comml．CommPort=6

　　Comml．Settings=″57600，N，8．1″

　　Comml．PortOpen=True

　　strsmsdata=Comml．Input′去除串口Buffer中的暂存信息

　　Timerl_Timer

　　EndSub

　　PrivateSubTimerLTimer()

　　DimIngNextvbCrLfAsLong′记录下一个回车换行符位置

　　DimIngFirstGPGGAAsLong

　　DimMessageAsString

　　DimsDateAsString

　　Dimcn

　　strsmsdata=Comml．Input

　　IngFirstGPGGA=Instr(1，strsmsclata，″＄GPGGA″)

　　IflngFirstGPGGAThen

　　IngNextVbCrLf=lnstr(IngFirstGPGGA+1，

　　strsmsdata，vbCrLf)

　　Message=Mid(Strsmsdata，ingFirstGPGGA，lngNextvb

　　CrLf-ingFirstGPGGA)

　　ParseNMEAMessageMessage，dblXCoord，dblYCoord，sDate

　　ProcessInputdblXCoord,dblYCoord

　　EndTf

　　GPS坐标变换

　　GPS定位所得出的结果属于WGS84坐标，而在工程上实用的大多是国家坐标系，或者是独立坐标系。独立坐标系一般是在国家坐标系基础上形成的，因此，GPS定位结果的使用中就有与国家坐标系间的坐标转换问题。一般要通过两步转换：首先将WGS84的大地坐标(L，B)转换为对应于WGS84椭球的高斯平面坐标(X84、Y84)，然后再经过平面坐标转换，将高斯平面坐标(X84、Y84)强制附合到本地高斯平面坐标系统。

　　高斯换算

　　将GPS定位所得出的大地坐标(L，B)转换为高斯平面坐标(X84、Y84)。有关的推导过程较复杂，本文只给出结果：

　　式中，ι=(L—L。)／p，ιo为投影带中央子午线经度，ρ=206265s/rad；酉圈曲率半径椭球第一偏心率e=2α-α2；辅助变量t=tgB；辅助变量η=e′cosB；椭圆第二偏心率分别为参考椭球的长短半径）；扁率α=（a-b）/a；X为赤道至纬度为B的平行圈的子午线弧长，其计算公式为