[問題] 影像處理noise的問題

[yanqinru (髍鵀)]

各位板友好

小弟最近開始學習影像處理

目前遇到了一些問題想請教各位

以下有Impulse、Gaussuan、Uniform三個雜訊程式段

想請問三個程式中for迴圈裡大概是在做什麼?

因為一直無法將理論和程式結合在一起

又不確定自己的想法到底有沒有錯誤

所以想問看看

不好意思麻煩各位了

(Impulse Noise)
prob = StrToInt(Form11->Edit1->Text);
level = StrToInt(Form11->Edit2->Text);
randomize();
for (int y = 0;y < nBitmap->Height;y++)
{
aptr = (Byte*)nBitmap->ScanLine[y];
for (int x = 0;x < nBitmap->Width;x++)
if (random(100) <= prob)
{
pix1 = cp * x;
aptr[pix1] = (Byte)(aptr[pix1]+level);
aptr[pix1+1] = (Byte)(aptr[pix1+1]+level);
aptr[pix1+2] = (Byte)(aptr[pix1+2]+level);
}
}



(Gaussian Noise)
mean = StrToInt(Form11->Edit1->Text);
sd = StrToFloat(Form11->Edit2->Text);
//ShowMessage("Mean = "+IntToStr(mean)+",\nSD ="+FloatToStr(sd)+".");
randomize();
for (int y = 0;y < nBitmap->Height;y++)
{
aptr = (Byte*)nBitmap->ScanLine[y];
for (int x = 0;x < nBitmap->Width;x++)
{
total = 0.0;
for (int k = 0;k < 10;k++)
{
ran = (float)(random(100)+1) / 100;
total += (float)sqrt12 * sd * (ran - 0.5);
}
total /= (float)sqrt10;
total += (float)mean;
pix1 = cp * x;
sum1 = total+(float)aptr[pix1];
sum2 = total+(float)aptr[pix1+1];
sum3 = total+(float)aptr[pix1+2];
if (sum1 < 0.0) sum1 = 0.0;
if (sum2 < 0.0) sum2 = 0.0;
if (sum3 < 0.0) sum3 = 0.0;
aptr[pix1] = (Byte)(sum1);
aptr[pix1+1] = (Byte)(sum2);
aptr[pix1+2] = (Byte)(sum3);
}
}




(Uniform Noise)
mean = StrToInt(Form11->Edit1->Text);
sd = StrToFloat(Form11->Edit2->Text);
randomize();
for (int y = 0;y < nBitmap->Height;y++)
{
aptr = (Byte*)nBitmap->ScanLine[y];
for (int x = 0;x < nBitmap->Width;x++)
{
ran = (float)(random(100)+1) / 100;
total = ran - 0.5;
total *= (float)sqrt12;
total *= sd;
total += (float)mean;
pix1 = cp * x;
sum1 = total+(float)aptr[pix1];
sum2 = total+(float)aptr[pix1+1];
sum3 = total+(float)aptr[pix1+2];
if (sum1 < 0.0) sum1 = 0.0;
if (sum2 < 0.0) sum2 = 0.0;
if (sum3 < 0.0) sum3 = 0.0;
aptr[pix1] = (Byte)(sum1);
aptr[pix1+1] = (Byte)(sum2);
aptr[pix1+2] = (Byte)(sum3);
}
}

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 203.64.186.243
推 pinchieh:給公式 可以中斷偵錯看
→ VictorTom:原本的雜訊是某種機率下會使pixel level跳離原始得值.
→ VictorTom:不同雜訊有不同機率分布, 也有不同的雜訊強度(如salt &
→ VictorTom:papper), 像一樓説的, 把公式拿出來對著code看一下:)
推 VictorTom:提供一個應該是小弟以前的實作方式, 看loop裡要幹麻:
→ VictorTom:1.loop每個pixel; 2.計算亂數決定要不要change這個pixel
→ VictorTom:3.(可能要再跑亂數)決定這個值要+-或改成多少; 4.寫回去
推 VictorTom:喵的又記錯了Orz 是先依照noise公式產生了這個pixel從
→ VictorTom:原值變化從-255~+255的機率密度函數分布結果, 然後才跑
→ VictorTom:loop對每個pixel算random值, 查表決定level要怎麼變化:)
推 VictorTom:記得Prentice Hall出板的, Gonzalez and Woods寫的這本
→ VictorTom:Digital Image Processing這本書裡有介紹要怎麼計算^^
→ yanqinru:感謝p大和V大的指點 我再努力研究看看  非常感謝