matlab练习程序（渲染三原色）

想成功必须强大

这里我用的空间是x向右为正，y向下为正，z向屏幕里面为正。相当于标准右手系绕x轴旋转了180度。

将三个点光源放在

r = [0.3,0,0.5];
g = [0.3,-0.5*cos(pi/6),-0.5*sin(pi/6)];
b = [0.3,0.5*cos(pi/6),-0.5*sin(pi/6)];

这三个位置上，向四周发射光线，取光线到y-z平面的模的倒数作为光的强度。

图像如下：

程序如下：

1. clear all;
   
2. 
   
3. close all;
   
4. clc;
   
5. 
   
6. r = [0.3,0,0.5];
   
7. g =
   
8. [0.3,-0.5*cos(pi/6),-0.5*sin(pi/6)];
   
9. b =
   
10. [0.3,0.5*cos(pi/6),-0.5*sin(pi/6)];
    
11. imgr1=[];
    
12. imgr2=[];
    
13. imgg1=[];
    
14. imgg2=[];
    
15. imgb1=[];
    
16. imgb2=[];
    
17. 
    
18. for y=-1:0.004:1
    
19. imgr1=[];
    
20. imgg1=[];
    
21. imgb1=[];  
    
22. for z=-1:0.004:1      
    
23.        pm=[0,y,z];         
    
24.        lightr=pm-r;
    
25.        lightg=pm-g;
    
26.        lightb=pm-b;
    
27.        imgr=1/norm(lightr);
    
28.        imgg=1/norm(lightg);
    
29.        imgb=1/norm(lightb);
    
30.       
    
31.        imgr1=[imgr1 imgr];
    
32.        imgg1=[imgg1 imgg];      
    
33.        imgb1=[imgb1 imgb];           
    
34.   end
    
35.    imgr2=[imgr2;imgr1];
    
36.    imgg2=[imgg2;imgg1];  
    
37.    imgb2=[imgb2;imgb1];
    
38. end[h w]=size(imgr2);
    
39. 
    
40. img=zeros(h,w,3);
    
41. img(:,:,1) = mat2gray(imgr2);
    
42. img(:,:,2) = mat2gray(imgg2);
    
43. img(:,:,3) = mat2gray(imgb2);
    
44. imshow(img)
    
45. imwrite(img,'img.png')

复制代码