作者yanqinru (髍鵀)
看板C_and_CPP
標題[問題] 影像處理noise的問題
時間Thu Apr 9 09:50:06 2009
各位板友好
小弟最近開始學習影像處理
目前遇到了一些問題想請教各位
以下有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
1F：推 pinchieh:給公式 可以中斷偵錯看 04/09 10:55
2F：→ VictorTom:原本的雜訊是某種機率下會使pixel level跳離原始得值. 04/09 11:22
3F：→ VictorTom:不同雜訊有不同機率分布, 也有不同的雜訊強度(如salt & 04/09 11:23
4F：→ VictorTom:papper), 像一樓説的, 把公式拿出來對著code看一下:) 04/09 11:23
5F：推 VictorTom:提供一個應該是小弟以前的實作方式, 看loop裡要幹麻: 04/09 11:25
6F：→ VictorTom:1.loop每個pixel; 2.計算亂數決定要不要change這個pixel 04/09 11:27
7F：→ VictorTom:3.(可能要再跑亂數)決定這個值要+-或改成多少; 4.寫回去 04/09 11:27
8F：推 VictorTom:喵的又記錯了Orz 是先依照noise公式產生了這個pixel從 04/09 11:48
9F：→ VictorTom:原值變化從-255~+255的機率密度函數分布結果, 然後才跑 04/09 11:49
10F：→ VictorTom:loop對每個pixel算random值, 查表決定level要怎麼變化:) 04/09 11:50
11F：推 VictorTom:記得Prentice Hall出板的, Gonzalez and Woods寫的這本 04/09 11:51
12F：→ VictorTom:Digital Image Processing這本書裡有介紹要怎麼計算^^ 04/09 11:52
13F：→ yanqinru:感謝p大和V大的指點 我再努力研究看看 非常感謝 04/09 12:56