Sat, 20 Oct 2007 14:33:00 GMT

摘要: 该算法有几个可学习的地方�Q�（1�Q�正�?思想�Q?�Q�对边界条�g的处理（3�Q�数据结构的选择code:sweep.h #ifndef SWEEP_H#define SWEEP_Hstruct Edge { int nxty; int curx; ... 阅读全文

Thu, 11 Oct 2007 03:57:00 GMT

昨晚��d��书馆看了《计��机囑�Ş学——OpenGL实现》关于Bresenham��法的另一�U�推导方式�?br>Bresenham最�_�֦�之处在于通过方程变换�Q�然后得到�P代方�E�，从而消除了��点�q�算�?br>
下面��单写写自己对中点法推导的理解�Q?br>
讎ͼ�W = bx - ax, H = by - ay
         所�?�Q�ax, ay�Q�和�Q�bx, by�Q�的理想直线为：
         -W*(y-ay) + H*(x-ax) = 0

讎ͼ�函数 f(x, y) = -2*W*(y-ay) + 2*H*(x-ax);
         f(x,y)有如下性质�Q?br>         f(x, y) < 0, 那么(x, y)在直�U�上�?br>         f(x, y) > 0, 那么(x, y)在直�U�下�?br>
现考虑点L(Px+1, Py), 点U(Px+1, Py+1), 则LU中点M(Px+1, Py+1/2) 有：
         如果f(Mx, My) < 0, 则M在理想直�U�上方，所以选择L
         如果f(Mx, My) > 0, 则M在理想直�U�下方，所以选择U
则：
         f(Mx,My) = -2*w*(Py+1/2-ay) + 2*H*(Px+1-ax)
�?x从Px+1�U�d��到Px+2�Ӟ�� 考虑f变化M'和M'':
         M'�Q�在前一步没有增加y, M' = (Px+2, Py+1/2)
         M''�Q�在前一步增加了y, M' = (Px+2, Py+3/2)
对于 M'�Q?br>         f(M'x, M'y) = -2*w*(Py+1/2-ay) + 2*H*(Px+2-ax) = f(Mx, My) + 2 * H
对于 M''�Q?br>         f(M''x, M''y) = -2*w*(Py+3/2-ay) + 2*H*(Px+2-ax) = f(Mx, My) - 2 * (W-H)
所�?br>         对于下一�?#8220;��试�?#8221;都有一个常数增量：前一�ơ没有增加y�Q�增量�ؓ2*H�Q�如果增加了y�Q�则增量�?2*(W-H)

对于初始条�g�Q�x = ax, y = ay
         M = (ax+1, ay+1/2);
         f(Mx, My) = -2*W*(ay+1/2-ay) + 2*H(ax+1-ax) = 2*H-W

Code:

#include <stdlib.h>
#include <math.h>
#include <GL/glut.h>

void myInit() {
    glClearColor(1.0, 1.0, 1.0, 0.0);
    glColor3f(0.0, 0.0, 0.0);
    //glPointSize(2.0);
    glMatrixMode(GL_PROJECTION);
    glLoadIdentity();
    gluOrtho2D(0.0, 640.0, 0.0, 480.0);
}

void setPixel(int x, int y) {
    glBegin(GL_POINTS);
    glVertex2i(x, y);
    glEnd();

}

void lineBres(int xs, int ys, int xe, int ye) {
    int W = xe - xs, H = ye - ys, f = 2 * H - W, tH = 2 * H, tHW = 2 * (H - W);
    int x, y;
    if (xs > xe) {
        x = xe;
        y = ye;
        xe = xs;
    } else {
        x = xs;
        y = ys;
    }
    while (x <= xe) {
        setPixel(x, y);
        x++;
        if (f<0) {
            f += tH;
        } else {
            y++;
            f += tHW;
        }
    }
}

void myDisplay() {
    glClear(GL_COLOR_BUFFER_BIT);
    lineBres(20, 10, 300, 180);
    glFlush();
}

int main(int argc, char **argv) {
    glutInit(&argc, argv);
    glutInitDisplayMode(GLUT_SINGLE|GLUT_RGB);
    glutInitWindowSize(640, 480);
    glutInitWindowPosition (100, 150);
    glutCreateWindow("Bresenham�ȝ��");
    glutDisplayFunc(myDisplay);
    myInit();
    glutMainLoop();
    return 0;
}

�?/a> 2007-10-11 11:57 发表评论

亚洲黄色免费电影,另类尿喷潮videofree,亚洲伊人网站