通过MATLAB实现基于冠状动脉造影图像序列的心脏及血管的三维光流运动估计

Good

1.仿真预览

 

2.部分核心代码

1. function [us,vs] = func_guangliu(images01,images02,areas);
3. imag01  = images01;
4. imag02  = images02;
5. %获取图像对边界进行调整
6. [r,c]   = size(images01);  
7. r_5     = r-5;
8. c_5     = c-5;
9. z       = zeros(r,c);
10. v1      = z;
11. v2      = z;
12. %初始化
13. gradt  = v1;
14. gradx  = v1;
15. grady  = v1;
16. alpha2  = 625;
17. times   = 20;
19. %计算梯度
20. gradt(5:r_5,5:c_5) = ( imag02(6:r_5+1,6:c_5+1)-imag01(6:r_5+1,6:c_5+1) + imag02(6:r_5+1,5:c_5)   - imag01(6:r_5+1,5:c_5)   + imag02(5:r_5  ,6:c_5+1) - imag01(5:r_5  ,6:c_5+1) + imag02(5:r_5  ,5:c_5)   - imag01(5:r_5 , 5:c_5))/4;
21. grady(5:r_5,5:c_5) = ( imag02(6:r_5+1,6:c_5+1)-imag02(5:r_5  ,6:c_5+1) + imag02(6:r_5+1,5:c_5)   - imag02(5:r_5  ,5:c_5)   + imag01(6:r_5+1,6:c_5+1) - imag01(5:r_5  ,6:c_5+1) + imag01(6:r_5+1,5:c_5)   - imag01(5:r_5 , 5:c_5))/4;
22. gradx(5:r_5,5:c_5) = ( imag02(6:r_5+1,6:c_5+1)-imag02(6:r_5+1,5:c_5)   + imag02(5:r_5  ,6:c_5+1) - imag02(5:r_5  ,5:c_5)   + imag01(6:r_5+1,6:c_5+1) - imag01(6:r_5+1,5:c_5)   + imag01(5:r_5  ,6:c_5+1) - imag01(5:r_5 , 5:c_5))/4;
28. for i=1:times
29.    delta =(gradx.*v1+grady.*v2+gradt)./(alpha2+gradx.^2+grady.^2);
30.    v1    = v1-gradx.*delta;
31.    v2    = v2-grady.*delta;
32. end;
34. u              = z;
35. u(5:r_5,5:c_5) = v1(5:r_5,5:c_5);
36. v              = z;
37. v(5:r_5,5:c_5) = v2(5:r_5,5:c_5);
40. xskip          = round(r/areas);
41. [us1,vs1]  = size(u(1:xskip:r,1:xskip:c));
42. us             = zeros(us1,vs1);
43. vs             = us;
45. N              = xskip^2;
47. for i=1:us1-1
48.   for j=1:vs1-1
49.       hk      = i*xskip-xskip+1;
50.       hl      = i*xskip;
51.       wk      = j*xskip-xskip+1;
52.       wl      = j*xskip;
53.       us(i,j) = sum(sum(u(hk:hl,wk:wl)))/N;
54.       vs(i,j) = sum(sum(v(hk:hl,wk:wl)))/N;
55.    end;
56. end;

复制代码

C31

来源：https://blog.csdn.net/ccsss22/article/details/123540836
免责声明：如果侵犯了您的权益，请联系站长，我们会及时删除侵权内容，谢谢合作！