这位同学真有缘，我今天正打算研究这段代码是什么鸟意思。想知道具体回答qq上解答 case 'motion' % Motion filter uses bilinear interpolation关键思想是双线性差分？ len = max(1,p2);输入的模糊核的长度必须大于1 half = (len-1)/2;% rotate half length around center求出模糊核一半的长度 phi = mod(p3,180)/180*pi;这里的输出角度全部都为0到180度，这也是后面sy求解中的不使用ysign的缘由 cosphi = cos(phi); sinphi = sin(phi); xsign = sign(cosphi);确定网格中x轴的生成方向 linewdt = 1;设置的模糊核的宽度，为什么要设置为1？ % define mesh for the half matrix, eps takes care of the right size其中eps用于保证mesh有正确的size？Eps是matlab中能计算的最小的数，比这个数小的数就认为是0. % for 0 & 90 rotation sx = fix(half*cosphi + linewdt*xsign - len*eps);fix这个函数是将小数向0方向舍入，如果half*cosphi包含小数则会因为向零舍入造成矩阵的size不能包含全部的长度，所以通过linewdt*xsign进行修正。 sy = fix(half*sinphi + linewdt - len*eps);与sx类似。由于eps的实际值很小，即使乘以len也不存在较大的值，所以我认为其影响只在某些特定的条件下产生：当sinphi=0时如果没有eps的存在，则sy=1，通过矩阵变换最终的模糊核有3行而不是一行，显然一行的情况更合乎逻辑，所以使用eps就是为了这点，但是为何乘以len还不得而知。[x y] = meshgrid([0:xsign:sx],[0:sy]);mesh这条命令在这个应用里面起到了两层作用：1.产生一个能够容纳线段的矩阵，也就是说这个矩阵有合适的size；2.矩阵的元素是坐标！所以可以通过矩阵元素进行坐标方面的运算。 % define shortest distance from a pixel to the rotated line dist2line = (y*cosphi-x*sinphi);% distance perpendicular to the line这条语句就使用了mesh产生的坐标性质，计算这些离散的点距离我们实际的这条线段的距离，直观地认为在直线上的点最好，如果没有直线上的点的话，我们可以选取离直线近的点来充当直线成分。后面就是根据这些距离与linewdt的大小比较获取的。 rad = sqrt(x.^2 + y.^2);计算mesh区域中离散点到原点的距离，主要用于区别超出线段长度的点。这一点也是利用mesh出来的矩阵的坐标特性。 % find points beyond the line's end-point but within the line widthlastpix = find((rad >= half)&(abs(dist2line)<=linewdt));找出距离原点超出half同时距离直线距离在linewdt宽度内的离散点坐标。这些点留作进一步判断，判断其是否还能用作表示线段的点。 %distance to the line's end-point parallel to the line x2lastpix = half - abs((x(lastpix) + dist2line(lastpix)*sinphi)/cosphi);dist2line(lastpix) = sqrt(dist2line(lastpix).^2 + x2lastpix.^2);该句连同上句是计算lastpix中的点到half末端的距离。后面的判断机制就是看这个距离是否在linewdt之内，在内则可以将该点算作线段内部，不在则将其置为0.dist2line = linewdt + eps - abs(dist2line); dist2line(dist2line<0) = 0;% zero out anything beyond line width % unfold half-matrix to the full size h = rot90(dist2line,2);后面这些就是将矩阵进行翻转等，补充一个完整的拥有len长度的线段。 h(end+[1:end]-1,end+[1:end]-1) = dist2line; h = h./(sum(h(:)) + eps*len*len); if cosphi>0, h = flipud(h); end 我只能按照代码进行解释，如果还有看不懂的，自己举一个实例就能搞懂，这种算法的原理我还找不到相关文献有一个理论解释，如果你能找到能告诉我吗？谢谢了。
|
