地理坐标与屏幕坐标转换（VC++）

在底层GIS开发中，我们经常会碰到空间数据可视化问题，当我们拿到地理数据之后，我们该怎么把它展现在我们的计算机屏幕中呢，这或许是一个比较简单的问题。

我们知道，地理坐标系是我们用户自定义的坐标系，而显示在计算机屏幕上的坐标是基于屏幕坐标来绘制的，那么改如何转化呢？

首先，我们来熟悉几个概念：

地理坐标系，通常原点在左下角，X轴的方向是向右递增，Y轴的方向是向上递增。

屏幕坐标系，通常它的原点在左上角，X轴的方向向右递增，Y轴的方向是向下递增的。正好与地理坐标系相反。

 

第一种方法：

中心点配对法，将视图的中心和地理空间的中心点匹配，然后给定一个地理坐标根据该坐标和视野中心点的偏移和比例尺来进行计算就可以获得屏幕坐标，反之也是如此；

 

关于这种方法，可以参考其他的资料。

其实这种方法也很简单，像ArcGIS等GIS软件就是基于这种转换方式进行的。其代码如下：

地理坐标转屏幕坐标

	OGRPoint ptCenter = GetMapCenter();		//地图中心点
	POINT pt;
	pt.x = draw->m_ScreenWidth/2 + (long)((poPoint->getX()-ptCenter.getX())*draw->m_Scale+0.5);
	pt.y = draw->m_ScreenHeigh/2 - (long)((poPoint->getY()-ptCenter.getY())*draw->m_Scale+0.5);
	return pt;

 

屏幕坐标转地理坐标

	OGRPoint ptCenter = GetMapCenter();		//地图中心点
	POINT pt;
	pt.x = draw->m_ScreenWidth/2 + (long)((poPoint->getX()-ptCenter.getX())*draw->m_Scale+0.5);
	pt.y = draw->m_ScreenHeigh/2 - (long)((poPoint->getY()-ptCenter.getY())*draw->m_Scale+0.5);
	return pt;

下过如下图所示：

第二种方法：利用相似比（完整显示地图）

 

这幅图是网上的一副图片，我希望没有版权问题。

通常情况下，我们取X和Y方向比例因子较大的作为比例因子，即

double xScale = rect.Width()/width;		//x方向的比例系数
				double yScale = rect.Height()/height;	//y方向的比例系数
				double scale = max(xScale,yScale);

从而地理坐标转换为屏幕坐标的代码如下：

pts[index].x = (int)((poPoint->getX()-minx)*scale);
 				pts[index].y = (int)(rect.Height()-(poPoint->getY()-miny)*scale);

在GDI设备描述表（device context,DC）包含了特定的显示设备的信息、图像的显示属性及特定的显示窗口等信息，用GDI编程的第一步必须获取当前的设备描述表的句柄（hDC），然后将该句柄作为参数调用GDI提供的API函数。我们在绘图时，通常是各个地理实体分别绘制，以线串为例：

//以下代码是根据直线的方法绘制
			for (int index = 0; index < numPoints - 1; index ++)
			{
				dcMem.MoveTo(pts[index]);
				dcMem.LineTo(pts[index + 1]);
			}

这只是在这方面的初步探索，还有就是现实效率不高，当元素特别多时，就显得非常慢，在以后会逐步加入空间索引机制，加快图形搜索和显示效率，下面是运行结果的截图，