0

我正在尝试使用图像 j 实现中值滤波。我在填充零时遇到问题,因为它在图片的底部和最左侧添加了额外的零。

这是我到目前为止所做的,如果你们能帮助我:

Dialog.create("9x9 median filtering");
Dialog.addMessage("9x9 median filtering");
Dialog.show();
setBatchMode(true);

median_filter_9();

setBatchMode("exit and display");

// Produce the 9x9 median image
function median_filter_9() 
{
 width = getWidth();
 height= getHeight();

 //if you want to apply this median filter to 16bit 
 depth = bitDepth();
 nBin= pow(2, depth);
 //nBin hold max gray intensity value 
 filteHisto = newArray(nBin);
 //filteHisto = newArray(255);

 fiveBYFive = newArray(81);
 //this is what i used for middle position of array to get median
 middlePos = round(81/2);


//-3, -3 will get you position 0,0 of a 9x9 matrix if you start in the middle
 for(j=-2;j<width-2;j++){
  for(i=-2;i<height-2;i++){ 
   z=0;
   for(r=0;r<9;r++){
    for(c=0;c<9;c++){
     //Extend outside image boundaries using zero padding.
    //error here: adds extra to bottom and farleft of picture
     if(j+r<0||j+r>=width||i+c<0||i+c>=height){
      fiveBYFive[z]=0;
      z++;
     }else{  
      v = getPixel(j+r,i+c);
      fiveBYFive[z]= v;
      z++;
     }
    }
   }
   //sort the array to find median
   Array.sort(fiveBYFive);
   median = fiveBYFive[middlePos];
   setPixel(j, i, median);
  } 
  updateDisplay();
 }



}
4

1 回答 1

0

您在图像边缘看到的一个问题是因为您正在用零填充 9x9 窗口,但您仍然将中值作为 81 项目窗口的中间值。

因此,例如,在图像的第一列中,您对至少 36 个元素进行零填充(顶部和底部更多),这意味着您只需要在图像中再找到 4 或 5 个零像素即可使中值元素零。

最简单的解决方法是根据添加的零数量向上调整中值元素的索引(在每次迭代中初始化为 81/2),或者只计算使用了多少非零像素,然后在中间找到中值排序数组中的范围(考虑排序顺序)。

通过这种方式,您可以获取找到的实际像素的中值并忽略填充的零。

可能,您错过了将代码从原始 5x5 更改为 9x9,因为开始/结束索引在任何情况下都是错误的,应该是

   for(j=-4;j<width;j++){
     for(i=-4;i<height;i++){ 

后面另一个可能的混淆来源是这条线,看起来你混淆了宽度和高度

if(j+r<0||j+r>=width||i+c<0||i+c>=height)

如果j是列索引并且i是行索引,应该是

if(j+c<0||j+c>=width||i+r<0||i+r>=height)

尽管对于方形窗口,这在实践中实际上并没有任何区别。

于 2013-06-04T04:56:53.707 回答