如何用简单的算法生成一个类似『光盘』的彩色圆形图片？

大概就像下图所示的样子

回复内容：

补充一下MATLAB的代码：

<span class="n">t</span> <span class="p">=</span> <span class="p">(</span><span class="mi">0</span><span class="p">:.</span><span class="mi">02</span><span class="p">:</span><span class="mi">2</span><span class="p">)</span><span class="o">*</span><span class="nb">pi</span><span class="p">;</span>
<span class="n">r</span> <span class="p">=</span> <span class="mi">0</span><span class="p">:.</span><span class="mi">02</span><span class="p">:</span><span class="mi">1</span><span class="p">;</span>
<span class="n">pcolor</span><span class="p">(</span><span class="nb">cos</span><span class="p">(</span><span class="n">t</span><span class="p">)</span><span class="o">'*</span><span class="n">r</span><span class="p">,</span><span class="nb">sin</span><span class="p">(</span><span class="n">t</span><span class="p">)</span><span class="o">'*</span><span class="n">r</span><span class="p">,</span><span class="n">t</span><span class="o">'*</span><span class="p">(</span><span class="n">r</span><span class="o">==</span><span class="n">r</span><span class="p">))</span>
<span class="n">colormap</span><span class="p">(</span><span class="n">hsv</span><span class="p">(</span><span class="mi">256</span><span class="p">)),</span> <span class="n">shading</span> <span class="n">interp</span><span class="p">,</span> <span class="n">axis</span> <span class="n">image</span> <span class="n">off</span>

...被抢先了...

基本上就是在这个圆上颜色的色相(Hue)，只与该点与圆心连线的夹角有关。知道这一点你就做一幅图像，然后遍历每一个点，把该点的坐标值转换为极坐标，色相值就是极坐标角度。

用 matplotlib 实现的话有个小技巧，把imshow的cmap改为hsv就可以直接按照色相来画。

<span class="c"># -*- coding:utf-8 -*-</span>
<span class="kn">from</span> <span class="nn">pylab</span> <span class="kn">import</span> <span class="o">*</span>

<span class="n">center</span> <span class="o">=</span> <span class="p">(</span><span class="mi">250</span><span class="p">,</span> <span class="mi">250</span><span class="p">)</span>
<span class="n">radius</span> <span class="o">=</span> <span class="mi">250</span>

<span class="n">img</span> <span class="o">=</span> <span class="n">zeros</span><span class="p">((</span><span class="mi">500</span><span class="p">,</span><span class="mi">500</span><span class="p">))</span>

<span class="k">for</span> <span class="n">i</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="mi">500</span><span class="p">):</span>
    <span class="k">for</span> <span class="n">j</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="mi">500</span><span class="p">):</span>
        <span class="n">x</span> <span class="o">=</span> <span class="n">i</span> <span class="o">-</span> <span class="n">center</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">*</span> <span class="mf">1.0</span>
        <span class="n">y</span> <span class="o">=</span> <span class="n">j</span> <span class="o">-</span> <span class="n">center</span><span class="p">[</span><span class="mi">1</span><span class="p">]</span> <span class="o">*</span> <span class="mf">1.0</span>
        <span class="k">if</span> <span class="n">x</span><span class="o">**</span><span class="mi">2</span><span class="o">+</span><span class="n">y</span><span class="o">**</span><span class="mi">2</span> <span class="o"><</span> <span class="n">radius</span><span class="o">**</span><span class="mi">2</span><span class="p">:</span>
            <span class="k">if</span> <span class="n">x</span> <span class="o">></span> <span class="mi">0</span><span class="p">:</span>
                <span class="n">img</span><span class="p">[</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">]</span> <span class="o">=</span> <span class="n">arctan</span><span class="p">(</span><span class="n">y</span><span class="o">/</span><span class="n">x</span><span class="p">)</span>
            <span class="k">elif</span> <span class="n">x</span><span class="o"><</span><span class="mi">0</span> <span class="ow">and</span> <span class="n">y</span><span class="o">>=</span><span class="mi">0</span><span class="p">:</span>
                <span class="n">img</span><span class="p">[</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">]</span> <span class="o">=</span> <span class="n">arctan</span><span class="p">(</span><span class="n">y</span><span class="o">/</span><span class="n">x</span><span class="p">)</span> <span class="o">+</span> <span class="n">pi</span>
            <span class="k">elif</span> <span class="n">x</span><span class="o"><</span><span class="mi">0</span> <span class="ow">and</span> <span class="n">y</span><span class="o"><</span><span class="mi">0</span><span class="p">:</span>
                <span class="n">img</span><span class="p">[</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">]</span> <span class="o">=</span> <span class="n">arctan</span><span class="p">(</span><span class="n">y</span><span class="o">/</span><span class="n">x</span><span class="p">)</span> <span class="o">-</span> <span class="n">pi</span>
            <span class="k">elif</span> <span class="n">x</span><span class="o">==</span><span class="mi">0</span> <span class="ow">and</span> <span class="n">y</span><span class="o">></span><span class="mi">0</span><span class="p">:</span>
                <span class="n">img</span><span class="p">[</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">]</span> <span class="o">=</span> <span class="n">pi</span> <span class="o">/</span> <span class="mi">2</span>
            <span class="k">elif</span> <span class="n">x</span><span class="o">==</span><span class="mi">0</span> <span class="ow">and</span> <span class="n">y</span><span class="o"><</span><span class="mi">0</span><span class="p">:</span>
                <span class="n">img</span><span class="p">[</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">]</span> <span class="o">=</span> <span class="n">pi</span> <span class="o">/</span> <span class="o">-</span><span class="mi">2</span>
            <span class="k">elif</span> <span class="n">x</span><span class="o">==</span><span class="mi">0</span> <span class="ow">and</span> <span class="n">y</span><span class="o">==</span><span class="mi">0</span><span class="p">:</span>
                <span class="n">img</span><span class="p">[</span><span class="n">i</span><span class="p">,</span><span class="n">j</span><span class="p">]</span> <span class="o">=</span> <span class="mf">0.0</span>

<span class="k">print</span> <span class="n">img</span>
<span class="n">imshow</span><span class="p">(</span><span class="n">img</span><span class="p">,</span> <span class="n">cmap</span><span class="o">=</span><span class="n">cm</span><span class="o">.</span><span class="n">hsv</span><span class="p">)</span>
<span class="n">show</span><span class="p">()</span>

改编自@冯昱尧的代码

<span class="kn">from</span> <span class="nn">pylab</span> <span class="kn">import</span> <span class="o">*</span>
<span class="n">n</span><span class="o">=</span><span class="mi">300</span>
<span class="n">img</span><span class="o">=</span><span class="p">[[</span><span class="n">arctan2</span><span class="p">(</span><span class="n">x</span><span class="p">,</span><span class="n">y</span><span class="p">)</span> <span class="k">if</span> <span class="n">x</span><span class="o">*</span><span class="n">x</span><span class="o">+</span><span class="n">y</span><span class="o">*</span><span class="n">y</span><span class="o"><</span><span class="n">n</span><span class="o">*</span><span class="n">n</span> <span class="k">else</span> <span class="mi">0</span> <span class="k">for</span> <span class="n">y</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="o">-</span><span class="n">n</span><span class="p">,</span><span class="n">n</span><span class="p">)]</span> <span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="o">-</span><span class="n">n</span><span class="p">,</span><span class="n">n</span><span class="p">)]</span>
<span class="n">imshow</span><span class="p">(</span><span class="n">img</span><span class="p">,</span> <span class="n">cmap</span><span class="o">=</span><span class="n">cm</span><span class="o">.</span><span class="n">hsv</span><span class="p">)</span>
<span class="n">show</span><span class="p">()</span>

如果根据这幅图，“光盘”的意思就是每一条半径的颜色在色相Hue上渐变。
关于色相是什么（反正这里说的不是可以用来牺牲的那个色相了），可以参看Wiki：HSL和HSV色彩空间，那个H就是色相Hue。
有点像下面这幅图的横截面（摘自Wiki）
只不过明度Chroma没有变化。
这样的话，蛮简单的啊。对每个像素算出它所在半径的角度 0~2π 然后映射到Hue的值域上就是了。
具体实现的话，话说可以用OpenCV吗？反正一些图像库里面是有HSR到RGB的转换函数的。

（我一直都不明白所谓邀请是怎么回事，以及为什么大家都爱说谢邀） 补充一下Mathematica代码

n=300;
c=Hue[i/n];
Graphics@Table[{c,EdgeForm@c,Disk[{0,0},1,2Pi/n{i-1,i}]},{i,n}]

期待mma的代码，估计一条语句搞定