| 查看: 301 | 回复: 0 | ||
[求助]
分水岭算法中计算输入图像的梯度这段代码是什么意思
|
|
先上代码，代码是我从网上下的,用MFC写的 void CMyImageDBDoc::GetGradient(BYTE* image, INT width, INT height , FLOAT* deltar, FLOAT* deltasita) //得到输入图像的梯度； { //下面计算各像素在水平和垂直方向上的梯度,边缘点梯度计为0； INT* deltaxarr; INT* deltayarr; INT grawidth = width; INT graheight = height; INT deltacount = grawidth * graheight; deltaxarr = new INT[deltacount]; deltayarr = new INT[deltacount]; //暂不计算边缘点； for (INT y=1; y for (INT x=1; x INT inarrpos = ((y)*width + (x))*3 + 1;//在输入块中的位置； INT deltaarrpos = y*grawidth + x;//在梯度数组中的位置； //卷积计算； deltaxarr[deltaarrpos] = (INT) ( ( image[((y-1)*width + (x+1))*3 + 1] //右上 + image[((y)*width + (x+1))*3 + 1] //右 + image[((y+1)*width + (x+1))*3 + 1] //右下 - image[((y-1)*width + (x-1))*3 + 1] //左上 - image[((y)*width + (x-1))*3 + 1] //左 - image[((y+1)*width + (x-1))*3 + 1] ) / 3 );//左下 deltayarr[deltaarrpos] = (INT) ( ( image[((y-1)*width + (x+1))*3 + 1] //右上 + image[((y-1)*width + (x))*3 + 1] //上 + image[((y-1)*width + (x-1))*3 + 1] //左上 - image[((y+1)*width + (x-1))*3 + 1] //左下 - image[((y+1)*width + (x))*3 + 1] //下 - image[((y+1)*width + (x+1))*3 + 1]) / 3 );//右下 } } //边缘赋为其内侧点的值； for (y=0; y INT x1 = 0; INT pos1 = y*grawidth + x1; deltaxarr[pos1] = deltaxarr[pos1+1]; deltayarr[pos1] = deltayarr[pos1+1]; INT x2 = grawidth-1; INT pos2 = y*grawidth + x2; deltaxarr[pos2] = deltaxarr[pos2-1]; deltayarr[pos2] = deltayarr[pos2-1]; } for (INT x=0; x INT y1 = 0; INT pos1 = x; INT inner = x + grawidth;//下一行； deltaxarr[pos1] = deltaxarr[inner]; deltayarr[pos1] = deltayarr[inner]; INT y2 = graheight-1; INT pos2 = y2*grawidth + x; inner = pos2 - grawidth;//上一行； deltaxarr[pos2] = deltaxarr[inner]; deltayarr[pos2] = deltayarr[inner]; } for (y=0; y for (x=0; x INT temppos = y*grawidth + x; if ( (deltaxarr[temppos])==0 ) { if (deltayarr[temppos]!=0) { deltasita[temppos] = 0;//水平方向; deltar[temppos] = (FLOAT) abs(deltayarr[temppos]); } else { deltasita[temppos] = -1;//无确定方向; deltar[temppos] = (FLOAT) abs(deltayarr[temppos]); } continue; } deltasita[temppos] = (FLOAT) ( atan( (FLOAT)deltayarr[temppos] / (FLOAT)deltaxarr[temppos] ) + PI/2. ); deltar[temppos] = (FLOAT) sqrt((DOUBLE) ( deltayarr[temppos]*deltayarr[temppos] + deltaxarr[temppos]*deltaxarr[temppos] ) ); } } delete [] deltaxarr; deltaxarr = NULL; //删除水平和垂直梯度数组； delete [] deltayarr; deltayarr = NULL; } 我想问的是， INT inarrpos = ((y)*width + (x))*3 + 1;//在输入块中的位置； 这段代码是什么意思呢，为什么要乘以3再加1，然后还有下面这整个循环啥意思，我也不明白，请懂的人详细给我说说这段代码 for (y=0; y for (x=0; x INT temppos = y*grawidth + x; if ( (deltaxarr[temppos])==0 ) { if (deltayarr[temppos]!=0) { deltasita[temppos] = 0;//水平方向; deltar[temppos] = (FLOAT) abs(deltayarr[temppos]); } else { deltasita[temppos] = -1;//无确定方向; deltar[temppos] = (FLOAT) abs(deltayarr[temppos]); } continue; } deltasita[temppos] = (FLOAT) ( atan( (FLOAT)deltayarr[temppos] / (FLOAT)deltaxarr[temppos] ) + PI/2. ); deltar[temppos] = (FLOAT) sqrt((DOUBLE) ( deltayarr[temppos]*deltayarr[temppos] + deltaxarr[temppos]*deltaxarr[temppos] ) ); } } [ 来自科研家族 哲学人生 ] |
回复此楼» 猜你喜欢
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有11人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
找到一些相关的精华帖子，希望有用哦~
SEM图片平均晶粒尺寸的计算
已经有7人回复
- 论文修回后编辑给出如下意见，有些看不懂什么意思，指高人指点
已经有25人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 求助hydridosilanes单词的意思？
已经有1人回复
- 服务器上计算，图形界面及图像处理问题思考
已经有19人回复
- 关于Dictra计算的问题（主要是参数调整，数据库数据添加）
已经有20人回复
- 求：如何用简便方法把图像的边缘进行清晰处理
已经有16人回复
- 【求助】PBE0 计算Si、Ge能带结构
已经有9人回复
- 【求助】哪位高手帮我翻译下这段英文算法的说明【已解决】
已经有5人回复
- 【资源】提高优化计算系统程序的应用软件！
已经有4人回复
- 【求助】算法和计算复杂性，就业前途如何？
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
30