三次线性插值//at3是什么意思。。。

oirrm

1. 
   
2. 
   
3. #include "stdafx.h"
   
4. 
   
5. ;
   
6. #define BM 0xff         //255  
   
7. #define N 0x1000        //4096  
   
8. #define B 0x100        //256  
   
9. 
   
10. int b00, b10, b01, b11;
    
11. int bx0, bx1, by0, by1, bz0, bz1;
    
12. float rx0, rx1, ry0, ry1, rz0, rz1;
    
13. float sx,sy,sz;
    
14. float t,u,v,a,b,c,d;
    
15. float* q;
    
16. float vec[3];
    
17. static int p[ B + B + 2];  
    
18. static float g3[B + B + 2][3];  
    
19. static float g2[B + B + 2][2];  
    
20. static float g1[B + B +2];  
    
21. int i,j;
    
22. 
    
23. #define lerp(t, a, b) ( a + t * (b - a) )  
    
24. #define s_curve(t) ( t * t * (3. - 2. * t) )  
    
25. #define setup(i,b0,b1,r0,r1)\
    
26.     t = vec[i] + N;\
    
27.     b0 = ((int)t) & BM;\
    
28.     b1 = (b0+1) & BM;\
    
29.     r0 = t - (int)t;\
    
30.     r1 = r0 - 1.;
    
31. #define at3(rx,ry,rz) ( rx * q[0] + ry * q[1] + rz * q[2] )
    
32. 
    
33. 
    
34. 
    
35. int _tmain(int argc, _TCHAR* argv[])
    
36. {
    
37.         setup(0, bx0, bx1, rx0, rx1); //0,1,0,-1
    
38.     setup(1, by0, by1, ry0, ry1); //0,1,0,-1
    
39.     setup(2, bz0, bz1, rz0, rz1); //0,1,0,-1
    
40.        
    
41.         i = p[bx0];  
    
42.     j = p[bx1];  
    
43.   
    
44.     b00 = p[i+by0];//得到四个随机值  
    
45.     b10 = p[j+by0];  
    
46.     b01 = p[i+by1];  
    
47.     b11 = p[j+by1];  
    
48.   
    
49.         sx = s_curve(rx0);  
    
50.     sy = s_curve(ry0);  
    
51.     sz = s_curve(rz0);
    
52. 
    
53.         q = g3[b00 + bz0];u = at3(rx0, ry0, rz0);//将(bx0,by0,bz0)处的梯度与中间点到该点的向量(rx0,ry0,rz0)相乘  
    
54.     q = g3[b10 + bz0]; v = at3(rx1, ry0, rz0);//将(bx1,by0,bz0)处的梯度与中间点到该点的向量(rx1,ry0,rz0)相乘  
    
55.     a = lerp(sx, u, v);  
    
56.   
    
57.     q = g3[b01 + bz0]; u = at3(rx0, ry1, rz0);//将(bx0,by1,bz0)处的梯度与中间点到该点的向量(rx0,ry1,rz0)相乘  
    
58.     q = g3[b11 + bz0];v = at3(rx1, ry1, rz0);//将(bx1,by1,bz0)处的梯度与中间点到该点的向量(rx1,ry1,rz0)相乘  
    
59.     b = lerp(sx, u, v);  
    
60.   
    
61.     c = lerp(sy, a, b);
    
62. 
    
63.         q = g3[b00 + bz1];u = at3(rx0, ry0, rz1);//将(bx0,by0,bz1)处的梯度与中间点到该点的向量(rx0,ry0,rz1)相乘  
    
64.     q = g3[b10 + bz1];v = at3(rx1, ry0, rz1);//将(bx1,by0,bz1)处的梯度与中间点到该点的向量(rx1,ry0,rz1)相乘  
    
65.     a = lerp(sx, u, v);  
    
66.   
    
67.     q = g3[b01 + bz1];u = at3(rx0, ry1, rz1);//将(bx0,by1,bz1)处的梯度与中间点到该点的向量(rx0,ry1,rz1)相乘  
    
68.     q = g3[b11 + bz0];v = at3(rx1, ry1, rz1);//将(bx1,by1,bz1)处的梯度与中间点到该点的向量(rx1,ry1,rz1)相乘  
    
69.     b = lerp(sx, u, v);  
    
70.   
    
71.     d = lerp(sy, a, b);  
    
72.   
    
73.     lerp(sz, c, d);//进行三次线性插值
    
74.         getchar();
    
75.         return 0;
    
76. }
    

复制代码