基本原理：
假定直线斜率0<k<1，并且直线上当前已确定一个像素点为P(xp, yp)，则下一个与理想直线最接近的像素点只能是P点正右方的点P1或右上方的点P2两者之一。设M为线段P1P2的中点，Q为理想直线与线段P1P2的交点。现需要确定下一个像素点。
若M在Q上方，则P1离直线更近，应取P1为下一个像素；
若M在Q下方，则P2离直线更近，应取P2为下一个像素；
若M与Q重合，则P1或P2任取一点。
这种以中点M作为判别标志的方法即为中点画线法。
假设直线段的起点为(x1, y1)，终点为(x2, y2)。设直线方程为：F (x, y)= ax + by + c ＝0
其中：a = y1 - y2，b = x2 - x1，c = x1 y2 - x2 y1
F(x,y)= 0时，则点在直线上；
当F(x,y)> 0时，则点在直线上方；
当F (x,y)< 0时，则点在直线下方。
要判断M在Q的上方还是下方，只需把M代入F(x, y)，并判断它的符号
构造判别式：
d＝F(M)=F(xp+1, yp+0.5)=a(xp+1)+b(yp+0.5)+c
若d<0，M在直线（Q点）下方，则取右上方的P2点作为下一个像素；
若d>0，M在直线（Q点）上方，则取正右方的P1点作为下一个像素；
若d=0，选P1点或P2点均可，约定取正右方P1点。
d是xp和yp的线性函数，可采用增量计算，以提高运算效率。
当d≥0时，取正右方像素P1，则再下一个像素的判别式为：
d1 = F(xp+2, yp+0.5)= a(xp+2)+b(yp+0.5)+c
= a(xp+1)+b(yp+0.5)+c+a = d + a
所以此时 d 的增量为 a
当d<0时，取右上方像素P2，则再下一个像素的判别式为：
d2 = F(xp+2, yp+1.5)=a(xp+2)+b(yp+1.5)+c
= a(xp+1)+b(yp+0.5)+c+a+b = d + a + b
所以此时 d 的增量为 a+b
d 的初始值问题：
假设直线的第一个像素为左端点(x1, y1)，则相应的判别式为
d0 = F(x1+1, y1+0.5)=a(x1+1)+b(y1+0.5)+c
= (ax1+by1+c)+a+0.5b = F(x1, y1)+a+0.5b
所以左端点(x1, y1)在直线上
所以 F (x1, y1) = 0，即 d 的初始值 d0 = a+0.5b
中点画线法的算法步骤（0≤k≤1）：
1. 输入直线段的两个端点P1(x1, y1)和P2(x2, y2)。
2. 初始化：a，b，d=a+0.5b，x=x1，y=y1，画点(x, y)。
3. 若x<x2，则执行下列各步，否则算法结束。
4. 判断d的符号；若d<0，则(x, y)更新为(x+1, y+1)，d更新为d+a+b；否则(x, y)更新为(x+1, y)，d更新为d+a。
5. 画点(x, y)，返回3。
改进：用2d代替d的中点画线法的算法步骤（0≤k≤1 ）
1. 输入直线段的两个端点P1(x1, y1)和P2(x2, y2)。
2. 初始化：a，b，d=2a+b，x=x1，y=y1，画点(x, y)。
3. 若x<x2，则执行下列各步，否则算法结束。
4. 判断d的符号：若d<0，则(x, y)更新为(x+1, y+1)，d更新为d+2a+2b；否则(x, y)更新为(x+1, y)，d更新为d+2a。
5. 画点(x, y)，返回3。
例：用中点画线法画直线段 P1(0, 0)—P2(5, 2)
a=y1-y2=-2 b=x2-x1=5
d0=2a+b=1 d1=2a=-4 d2=2(a+b)=6
代码（0<k<1）:
void MPLine(int x1,int y1,int x2,int y2,int color){
int x,y,a,b,c,d,d1,d2;
a=y1-y2;b=x2-x1;
y=y1;
d=2*a+b;d1=2*a;d2=2*(a+b);
putpixel(x,y,color);
k=1.0*(y2-y1)/(x2-x1);
for(x=x1;x<=x2;x++){
if(d<0){
y++;
d+=d2;
}else{
d+=d1;
}
putpixel(x,y,color);
}
}
