GIS基础

最近在摸索GIS方面的东西，感觉此网站不错
http://www.gaoang.com/

用过Google Maps的朋友都会认为它很Cool。其实创建一个类似的地图服务站点，在技术上难度并不大，已经有较为成熟的开源解决方案。但是能像Google那样有原始影像数据和庞大的集群支持，不是一般站点可以做到的。

这次把Google Maps搬出来，是想分享一下比较有意思的Google坐标换算方式。Google Maps使用Keyhole公司提供的卫星影像，全世界众多大城市的卫星图片可以达到1m以下的精度。Google Maps的卫星影像使用等角正切圆柱的墨卡托(Mercator)投影，被预先处理成按照不同精度划分的影象金字塔，提供了0-17共18级的缩放等级并进行四叉树编码，每张卫星图都由切片成256*256像素大小的影像组成，根据用户请求的位置拼接组成大的影像图。

Google Maps的坐标使用QRTS这四个字母进行编码，如右图所示，一张图片被分成不同字母标识的4块，根据请求的位置逐级细化，最后达到所需精度的卫星图片。 如果知道特定位置的经纬度，换算成Google Maps坐标并不复杂。首先将经纬度转换为弧度，规格化后使得变化区间在0－1之内，然后根据需要的图像精度级数确定迭代次数，进而计算出经纬度对应的 Google坐标字符串。

按照计算公式，北京的经纬度以天安门作为地理标准坐标，位于北纬39度54分27秒，东经116度23分17秒。换算后得到北京的Google坐标字符串是trstrqqrqssttttrrrstq，具体对应的卫星影像地址为：http://kh.google.com/kh?v=3&t=trstrqqrqssttt 网上提供了C和Python写的换算程序，感兴趣的朋友可以一窥究竟，方便的把经纬度换算成Google卫星影像坐标。