在我们的GIS开发中,为了更好地表达和展现空间数据,经常要进行符号的绘制,圆在地图上常常作为定位分级符号,其圆心定位在要素的中心位置,圆的面积表示相应的数量指标。
在计算机中,圆是由多边形逼近的,当多边形的边数多到人眼不能分辨的地步时,就显示成了圆形。然而,变数过多的话,可能会影响绘制的效果,边数过少就不像圆。所以一般的做法是考虑多边形上相邻三点构成的三角形看上去是否是一条直线。
具体的推导过程如下:
设圆的半径为R
本人是在GDI环境下绘制的,代码如下:
void CGeoSymbol::DrawCycle(CDC* pDC, int x0, int y0, double radis, double ang1, double ang2){ //先将初始角度和终止角度转化为弧度制单位 double ag1 = 0.01745*ang1; double ag2 = 0.01745*ang2; //计算逼近圆的多边形的边数和角度增量值 int N = (int)(fabs(ag2-ag1)/acos(1-0.015/radis)) + 1; double da = fabs(ag2-ag1)/N; //开始循环计算每个顶点的坐标 int i = 0; while (i <= N) { double w = ag1+i*da; double x = x0 + radis*cos(w); double y = y0 + radis*sin(w); if (i<=0) { pDC->MoveTo(x,y); } else { pDC->LineTo(x,y); } i ++; }}
但是绘制的效果也有一个问题,就是符号精度不高,因为GDI本身的绘制函数都只支持整型坐标值,当传入浮点型坐标时,也是先转化为整型坐标再绘制的,所以可以考虑用GDI+或者DX、OPENGL之类的绘图接口。