Giải thuật Xiaolin Wu
![](http://upload.wikimedia.org/wikipedia/commons/thumb/2/27/XiaolinWuLine.png/220px-XiaolinWuLine.png)
Giải thuật vẽ đoạn thẳng Xiaolin Wu, tiếng Anh: XiaolinWu's line algorithm là giải thuật vẽ đường thẳng khử răng cưa, được giới thiệu lần đầu tiên trên bài báo An Efficient Antialiasing Technique vào tháng 7 năm 1991 trên tờ báo Computer Graphics, cũng như trên bài báo Fast Antialiasing vào tháng 6 năm 1992 trên tờ Dr. Dobb's Journal.
Giải thuật Bresenham vẽ đoạn thẳng vẽ đường thẳng rất nhanh, tuy nhiên lại không có chức năng khử răng cưa. Hơn nữa, giải thuật không kiểm soát được trường hợp điểm cuối của đoạn thẳng không nằm trên một điểm có tọa độ nguyên. Các phương pháp khử răng cưa thường tốn rất nhiều thời gian xử lý, nhưng giải thuật của Xiaolin Wu thì rất nhanh (mặc dù vẫn chậm hơn giải thuật của Bresenham). Thuật toán cơ bản của giải thuật là vẽ các cặp điểm gần nhau hai bên đoạn thẳng và tô màu dựa trên độ ưu tiên. Hai điểm ở hai đầu đoạn thẳng được kiểm soát riêng. Đoạn thẳng ngắn hơn 1 pixel sẽ được xử lý trong trường hợp riêng.
Một giải thuật mở rộng để vẽ đường tròn được giới thiệu bởi Xiaolin Wu trên tạp chí Graphics Gems II. Tương tự như giải thuật vẽ đoạn, giải thuật này dùng để thay thế giải thuật vẽ đường tròn của Bresenham.
function plot(x, y, c) is plot the pixel at (x, y) with brightness c (where 0 ≤ c ≤ 1) // integer part of x function ipart(x) is return floor(x) function round(x) is return ipart(x + 0.5) // fractional part of x function fpart(x) is return x - floor(x) function rfpart(x) is return 1 - fpart(x) function drawLine(x0,y0,x1,y1) is boolean steep := abs(y1 - y0) > abs(x1 - x0) if steep then swap(x0, y0) swap(x1, y1) end if if x0 > x1 then swap(x0, x1) swap(y0, y1) end if dx := x1 - x0 dy := y1 - y0 gradient := dy / dx if dx == 0.0 then gradient := 1.0 end if // handle first endpoint xend := round(x0) yend := y0 + gradient * (xend - x0) xgap := rfpart(x0 + 0.5) xpxl1 := xend // this will be used in the main loop ypxl1 := ipart(yend) if steep then plot(ypxl1, xpxl1, rfpart(yend) * xgap) plot(ypxl1+1, xpxl1, fpart(yend) * xgap) else plot(xpxl1, ypxl1 , rfpart(yend) * xgap) plot(xpxl1, ypxl1+1, fpart(yend) * xgap) end if intery := yend + gradient // first y-intersection for the main loop // handle second endpoint xend := round(x1) yend := y1 + gradient * (xend - x1) xgap := fpart(x1 + 0.5) xpxl2 := xend //this will be used in the main loop ypxl2 := ipart(yend) if steep then plot(ypxl2 , xpxl2, rfpart(yend) * xgap) plot(ypxl2+1, xpxl2, fpart(yend) * xgap) else plot(xpxl2, ypxl2, rfpart(yend) * xgap) plot(xpxl2, ypxl2+1, fpart(yend) * xgap) end if // main loop if steep then for x from xpxl1 + 1 to xpxl2 - 1 do begin plot(ipart(intery) , x, rfpart(intery)) plot(ipart(intery)+1, x, fpart(intery)) intery := intery + gradient end else for x from xpxl1 + 1 to xpxl2 - 1 do begin plot(x, ipart(intery), rfpart(intery)) plot(x, ipart(intery)+1, fpart(intery)) intery := intery + gradient end end if end function
Ghi chú: Nếu trong lần chạy đầu tiên điều kiện abs(dx) < abs(dy) là đúng, quá trình vẽ sẽ thực hiện với x và y ngược nhau.
Tham khảo
- Abrash, Michael (1992). “Fast Antialiasing (Column)”. Dr. Dobb's Journal. 17 (6): 139(7).
- Wu, Xiaolin (1991). “An efficient antialiasing technique”. Computer Graphics. 25 (4): 143–152. doi:10.1145/127719.122734. ISBN 0-89791-436-8.
- Wu, Xiaolin (1991). “Fast Anti-Aliased Circle Generation”. Trong James Arvo (biên tập). Graphics Gems II. San Francisco: Morgan Kaufmann. tr. 446–450. ISBN 0-12-064480-0.
Liên kết ngoài
- Xiaolin Wu's homepage