24小时热门版块排行榜

返回列表

xiashaoyan

新虫 (正式写手)

应助: 10 (幼儿园)
金币: 2574.9
散金: 30
红花: 10
帖子: 498
在线: 184.1小时
虫号: 2507802
注册: 2013-06-14
专业: 计算机软件

[求助] cany算子在局部非极大值抑制的计算过程的问题！！

在cany算子中局部非极大值抑制的计算过程中，要进行对插值，计算的程序如下：
（程序来源于http://blog.csdn.net/xiajun07061225/article/details/6926108）

%局部非极大值抑制
K = d;%记录进行非极大值抑制后的梯度
%设置图像边缘为不可能的边缘点
for j = 1:width
K(1,j) = 0;
end
for j = 1:width
K(height,j) = 0;
end
for i = 2:width - 1
K(i,1) = 0;
end
for i = 2:width - 1
K(i,width) = 0;
end

for i = 2:height - 1
for j = 2:width - 1
      %当前像素点的梯度值为0，则一定不是边缘点
      if d(i,j) == 0
         K(i,j) = 0;
      else
         gradX = dx(i,j);%当前点x方向导数
         gradY = dy(i,j);%当前点y方向导数
         gradTemp = d(i,j);%当前点梯度
         %如果Y方向幅度值较大
         if abs(gradY) > abs(gradX)
            weight = abs(gradX) / abs(gradY);%权重
            grad2 = d(i - 1,j);
            grad4 = d(i + 1,j);
            %如果x、y方向导数符号相同
            %像素点位置关系
            %g1 g2
            % C
            % g4 g3
            if gradX * gradY > 0
                  grad1 = d(i - 1,j - 1);
                  grad3 = d(i + 1,j + 1);
            else
                  %如果x、y方向导数符号反
                  %像素点位置关系
                  % g2 g1
                  % C
                  %g3 g4
                  grad1 = d(i - 1,j + 1);
                  grad3 = d(i + 1,j - 1);
            end
         %如果X方向幅度值较大
         else
            weight = abs(gradY) / abs(gradX);%权重
            grad2 = d(i,j - 1);
            grad4 = d(i,j + 1);
            %如果x、y方向导数符号相同
            %像素点位置关系
            %g3
            %g4 C g2
            %    g1
            if gradX * gradY > 0
                  grad1 = d(i + 1,j + 1);
                  grad3 = d(i - 1,j - 1);
            else
                  %如果x、y方向导数符号反
                  %像素点位置关系
                  %    g1
                  %g4 C g2
                  %g3
                  grad1 = d(i - 1,j + 1);
                  grad3 = d(i + 1,j - 1);
            end
         end
         %利用grad1-grad4对梯度进行插值
         gradTemp1 = weight * grad1 + (1 - weight) * grad2;
         gradTemp2 = weight * grad3 + (1 - weight) * grad4;
         %当前像素的梯度是局部的最大值，可能是边缘点
         if gradTemp >= gradTemp1 && gradTemp >= gradTemp2
            K(i,j) = gradTemp;
         else
            %不可能是边缘点
            K(i,j) = 0;
         end
      end
end
end
我捉摸了一下午，感觉考虑两个条件，abs(dy)与abs(dx)的大小，以及dy与dx乘积的符号，感觉理解的不是很透彻，于是我画了个图，如见下：局部非极大值抑制 .jpg
画完了我就困惑了，因为图中第3个图，g1,g2不在一条直线上，怎么插值呢？所以想是不是我认为的x,y轴相反了呢？还是我画的梯度方向本身就不对呢？我查了一下，x,y轴在图像中的方向的确像我这样画的，但是我看了一个帖子http://bbs.csdn.net/topics/370004267中给了一个图，见图2，帖子里面说gx*gy>0是因为图像坐标系中y轴是朝下的，我就更加困惑了，这个帖子里面定义x,y轴方向与一般图像处理教材的方向恰好颠倒了，后来我也试了x,y轴颠倒过来，画了四种类似于上面第一个图局部非极大值抑制 .jpg，但是还是出现，g1,g2不在同一条直线上，接着我仔细看那程序，发现，考虑，abs(gradY) > abs(gradX) 与abs(gradX) > abs(gradY) 情况下，在考虑相同的gradX * gradY > 0下，grad1 与grad1 恰好颠倒过来，但是如果这样，g1,g2是没法再同一条直线上  ，所以感觉很困惑，求助各位研究算法和图像处理的前辈们！！！

局部非极大值抑制 .jpg

2.png

回复此楼

» 猜你喜欢

求助中医院临床中医中标的青年标书学习，谢谢已经有6人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有259人回复
物理学I论文润色/翻译怎么收费? 已经有209人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复
25届硕士申请，材料计算方向已经有2人回复
课题组招1-2名硕士已经有1人回复
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上已经有0人回复
斯坦福大学与Elsevier世界前2%学者新排名查询已经有1人回复
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有19人回复
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计已经有1人回复
apl投稿状态已经有2人回复

1楼 2015-01-24 23:32:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 xiashaoyan 的主题更新

返回列表