数字图像处理实验报告(4)

a.at(j,k) = floor(t.at(j,k)); } }

Mat proj1(a+ctrldx+1);

for(int j=0; j

for(int k=0; k(j,k);

proj1.at(j,k) = (float)proj.at(data-1,0); } }

Mat proj2(a+ctrldx);

for(int j=0; j

for(int k=0; k(j,k);

proj2.at(j,k) = (float)proj.at(data-1,0); } }

Mat temp = (t-a).mul(proj1) + (a+1-t).mul(proj2); img = img + temp; }

img = img*M_PI/(2*360);

normalize(img, img, 0, 1, CV_MINMAX);

imshow(\ waitKey(0); }

16