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

看板C_and_CPP作者 (髍鵀)時間15年前 (2009/04/09 09:50), 編輯推噓4(409)
留言13則, 3人參與, 最新討論串1/1
各位板友好 小弟最近開始學習影像處理 目前遇到了一些問題想請教各位 以下有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

04/09 10:55, , 1F
給公式 可以中斷偵錯看
04/09 10:55, 1F

04/09 11:22, , 2F
原本的雜訊是某種機率下會使pixel level跳離原始得值.
04/09 11:22, 2F

04/09 11:23, , 3F
不同雜訊有不同機率分布, 也有不同的雜訊強度(如salt &
04/09 11:23, 3F

04/09 11:23, , 4F
papper), 像一樓説的, 把公式拿出來對著code看一下:)
04/09 11:23, 4F

04/09 11:25, , 5F
提供一個應該是小弟以前的實作方式, 看loop裡要幹麻:
04/09 11:25, 5F

04/09 11:27, , 6F
1.loop每個pixel; 2.計算亂數決定要不要change這個pixel
04/09 11:27, 6F

04/09 11:27, , 7F
3.(可能要再跑亂數)決定這個值要+-或改成多少; 4.寫回去
04/09 11:27, 7F

04/09 11:48, , 8F
喵的又記錯了Orz 是先依照noise公式產生了這個pixel從
04/09 11:48, 8F

04/09 11:49, , 9F
原值變化從-255~+255的機率密度函數分布結果, 然後才跑
04/09 11:49, 9F

04/09 11:50, , 10F
loop對每個pixel算random值, 查表決定level要怎麼變化:)
04/09 11:50, 10F

04/09 11:51, , 11F
記得Prentice Hall出板的, Gonzalez and Woods寫的這本
04/09 11:51, 11F

04/09 11:52, , 12F
Digital Image Processing這本書裡有介紹要怎麼計算^^
04/09 11:52, 12F

04/09 12:56, , 13F
感謝p大和V大的指點 我再努力研究看看 非常感謝
04/09 12:56, 13F
文章代碼(AID): #19tLHFp8 (C_and_CPP)