画个圆吧---C语言画圆

http://blog.chinaunix.net/uid-24774106-id-3047520.html

我知道C语言的能力很强，但是对于图形图像类的C语言实现很是好奇，觉得很难。我读书的时候，用C语言画过BMP格式的图像，因为没有任何的压缩，所以图像的数据量比较大，基本上就是RGB 字节流。今天看了Banu blog中的博文 Drawing circles，想起了以前用C语言进行BMP图像的处理。

     写这篇博客，是给自己枯燥的学习带来一些乐趣，C语言也可以干些比较有有趣的事情。代码基本上是照搬Banu，只是，文章画的图像是灰度图像，不够靓丽，我将代码改成了生成彩色图像。不管怎么说，光荣属于Banu，我只是仰望前辈的学习者。

    生成一个文件，后缀为ppm，Linux下的GIMP自然可以解析这种格式。ppm格式和BMP很像，基本属于傻瓜型图像格式，每个像素就是RGB三个字节表示，当然为了声明文件类型，会有文件头。看下WIKI中的介绍。    

Magic Number	Type	Encoding
P1	Portable bitmap	ASCII
P2	Portable graymap	ASCII
P3	Portable pixmap	ASCII
P4	Portable bitmap	Binary
P5	Portable graymap	Binary
P6	Portable pixmap	Binary

P3 

# The P3 means colors are in ASCII, then 3 columns and 2 rows, 

# then 255 for max color, then RGB triplets

 3 2

 255

 255 0 0       0 255 0        0 0 255 

255 255 0      255 255 255    0 0 0

我们看到了上面虚框中的文件内容，画出来就是下面这个彩色图像。当然，下面的图用一个比较大的正方形来表示一个像素点（255，0，0）表示红色，所以代表第一个像素点的正方形是红色的，依次类推。

OK，下面看下代码：

1. static void
   
2. image_save (const Image *image,
   
3.             const char *filename)
   
4. {
   
5.   FILE *out;
   
6. 
   
7.   out = fopen (filename, "wb");
   
8.   if (!out)
   
9.     return;
   
10. 
    
11.   fprintf (out, "P6\n");
    
12.   fprintf (out, "%zu %zu\n", image->width, image->height);
    
13.   fprintf (out, "255\n");
    
14. 
    
15.   fwrite (image->data, 1, image->width * image->height*3, out);
    
16. 
    
17.   fclose (out);
    
18. }

P6告诉文件解析程序（如GIMP）我是PPM格式的文件（Portable pixmap）。

第二行告诉文件的长和高。

下面是完整的代码，再次强调，荣耀和版权都属于Banu，我只是仰望Banu的学习者 。

1. #include <stdio.h>
   
2. #include <stdlib.h>
   
3. #include <string.h>
   
4. #include <math.h>
   
5. 
   
6. typedef struct {
   
7.   size_t width;
   
8.   size_t height;
   
9.   unsigned char *data;
   
10. } Image;
    
11. 
    
12. static Image *
    
13. image_new (size_t width,
    
14.            size_t height)
    
15. {
    
16.   Image *image;
    
17. 
    
18.   image = malloc (sizeof *image);
    
19.   image->width = width;
    
20.   image->height = height;
    
21.   image->data = malloc (width * height*3);
    
22. 
    
23.   return image;
    
24. }
    
25. 
    
26. static void
    
27. image_free (Image *image)
    
28. {
    
29.   free (image->data);
    
30.   free (image);
    
31. }
    
32. 
    
33. static void
    
34. image_fill (Image *image,
    
35.             unsigned char value)
    
36. {
    
37.   memset (image->data, value, image->width * image->height*3);
    
38. }
    
39. 
    
40. /**
    
41.  * image_set_pixel:
    
42.  *
    
43.  * Sets a pixel passed in signed (x, y) coordinates, where (0,0) is at
    
44.  * the center of the image.
    
45.  **/
    
46. static void
    
47. image_set_pixel (Image *image,
    
48.                  ssize_t x,
    
49.                  ssize_t y,
    
50.                  unsigned char Rvalue,
    
51.                  unsigned char Gvalue,
    
52.                  unsigned char Bvalue)
    
53. {
    
54.   size_t tx, ty;
    
55.   unsigned char *p;
    
56. 
    
57.   tx = (image->width / 2) + x;
    
58.   ty = (image->height / 2) + y;
    
59. 
    
60.   p = image->data + (ty * image->width*3) + tx*3;
    
61. 
    
62.   *p = Rvalue;
    
63.     *(p+1) = Gvalue;
    
64.     *(p+2) = Bvalue;
    
65. }
    
66. 
    
67. static void
    
68. image_save (const Image *image,
    
69.             const char *filename)
    
70. {
    
71.   FILE *out;
    
72. 
    
73.   out = fopen (filename, "wb");
    
74.   if (!out)
    
75.     return;
    
76. 
    
77.   fprintf (out, "P6\n");
    
78.   fprintf (out, "%zu %zu\n", image->width, image->height);
    
79.   fprintf (out, "255\n");
    
80. 
    
81.   fwrite (image->data, 1, image->width * image->height*3, out);
    
82. 
    
83.   fclose (out);
    
84. }
    
85. 
    
86. static void
    
87. draw_circle (Image *image,
    
88.              int radius,
    
89.              unsigned char Rvalue,
    
90.                          unsigned char Gvalue,
    
91.                          unsigned char Bvalue)
    
92. {
    
93.   int x, y;
    
94. 
    
95.   for (y = -radius; y <= radius; y++)
    
96.     for (x = -radius; x <= radius; x++)
    
97.       if ((x * x) + (y * y) <= (radius * radius))
    
98.         image_set_pixel (image, x, y, Rvalue,Gvalue,Bvalue);
    
99. }
    
100. 
     
101. int
     
102. main (int argc, char *argv[])
     
103. {
     
104.   Image *image;
     
105. 
     
106.   image = image_new (600, 600);
     
107. 
     
108.   image_fill (image, 0xff);
     
109.   draw_circle (image, 200, 0x00,0x00,0xff);
     
110.   image_save (image, "circle.ppm");
     
111. 
     
112.   image_free (image);
     
113. 
     
114.   return 0;
     
115. }

就可以在本地查看一个蓝色的圆图像。下图格式是jpg格式的，我用QQ截图工具生成的。

由于CU不支持ppm格式的图像，所以只能传jpg了。

参考文献：

1 Banu 博文 Drawing Circles

2 Wiki