| 2 | 1/1 | 返回列表 |
| 查看: 1425 | 回復(fù): 1 | ||
彩虹豬新蟲 (小有名氣)
|
[求助]
ssr單尺度算法
|
|
求各位大神指點:謝謝 根據(jù)Retinex圖像增強公式: S(x,y) = R(x,y)×L(x,y)??S為原圖 R為要求的入射分量圖 L為亮度圖 取對數(shù) log(S(x,y)) = log(R(x,y)) + log(L(x,y)); 定義一個高斯函數(shù) G(x,y)=Ke^(-(x*x+y*y)/c*c)這里的c是卷積尺度c=80 并且一般設(shè)L(x,y)=G(x,y)*S(x,y);這里的*是卷積 然后得到增強后的圖像log(R(x,y))= log(S(x,y))??+??log(G(x,y)*S(x,y) ); 我自己編寫的程序如下:可是現(xiàn)實都不對啊 為什么呢 要不全黑 要么全白的 有沒有高手指點一下 困擾很久了 clear all;clc; I=imread('bs.jpg'); R=I(:,:,1);%獲取R通道 G=I(:,:,2);%獲取G通道 B=I(:,:,3);%獲取B通道 [m,n]=size(R);%獲取圖像大小 Rd=double(R);%R將整形轉(zhuǎn)換為浮點型 Gd=double(G);%G將整形轉(zhuǎn)換為浮點型 Bd=double(B);%B將整形轉(zhuǎn)換為浮點型 Rlog=log(Rd+1);%R將浮點型轉(zhuǎn)換為對數(shù) Glog=log(Gd+1);%G將浮點型轉(zhuǎn)換為對數(shù) Blog=log(Bd+1);%B將浮點型轉(zhuǎn)換為對數(shù) fftr=fft2(Rd);%R二維離散傅立葉變換 fftg=fft2(Gd);%G二維離散傅立葉變換 fftb=fft2(Bd);%B二維離散傅立葉變換 %建立高斯函數(shù)進行卷積運算 c=80;%C是高斯環(huán)繞尺度,λ是一個尺度 s=zeros(m,n); for i=1:m ? ? for j=1:n ? ?? ???s(i,j)=exp(-(i*i+j*j)/c*c); ? ? end end K=s/(sum(sum(s)));%求歸一化K g=K.*s;%高斯函數(shù) g=double(g); g=fft2(g,size(g,1),size(g,2));%二維離散傅立葉變換 %Rc=conv(g,fftr);%R進行卷積運算,獲得低頻圖像 %Gc=conv(g,fftg);%G進行卷積運算 %Bc=conv(g,fftb);%B進行卷積運算 Rc=g.*fftr;%R進行卷積運算,獲得低頻圖像 Gc=g.*fftg;%G進行卷積運算 Bc=g.*fftb;%B進行卷積運算 iffRc=ifft2(Rc);%R傅里葉反轉(zhuǎn) iffGc=ifft2(Gc);%G傅里葉反轉(zhuǎn) iffBc=ifft2(Bc);%B傅里葉反轉(zhuǎn) Rclog=log(iffRc+1);%轉(zhuǎn)換到對數(shù)域 Gclog=log(iffGc+1);%轉(zhuǎn)換到對數(shù)域 Bclog=log(iffBc+1);%轉(zhuǎn)換到對數(shù)域 gR=Rlog-Rclog;%用原圖像減去地頻圖像 gG=Glog-Gclog;%用原圖像減去地頻圖像 gB=Blog-Bclog;%用原圖像減去地頻圖像 expR=exp(gR);%取反對數(shù),得到增強的圖像 expG=exp(gG);%取反對數(shù),得到增強的圖像 expB=exp(gB);%取反對數(shù),得到增強的圖像 %圖像灰度拉伸 %此時的expR,G,B值的范圍并不是0–255,所以還需要進行線性拉伸并轉(zhuǎn)換成相應(yīng)的格式輸出顯示 MIN = min(min(expR));%min(expR)求出每一列中最小的值min(min(expR))求出整個矩陣中最小的值 MAX = max(max(expR));%max(expR)求出每一列中最大的值max(max(expR))求出整個矩陣中最大的值 expR = (expR - MIN)/(MAX - MIN);%灰度變換增強 expR = adapthisteq(expR);%通過使用對比度限制的自適應(yīng)直方圖均衡化來增強灰度圖像的對比度 MIN1 = min(min(expG));%min(expG)求出每一列中最小的值min(min(expG))求出整個矩陣中最小的值 MAX1 = max(max(expG));%max(expG)求出每一列中最大的值max(max(expG))求出整個矩陣中最大的值 expG = (expG - MIN1)/(MAX1 - MIN1);%灰度變換增強 expG = adapthisteq(expG);%通過使用對比度限制的自適應(yīng)直方圖均衡化來增強灰度圖像的對比度 MIN2 = min(min(expB));%min(expB)求出每一列中最小的值min(min(expB))求出整個矩陣中最小的值 MAX2 = max(max(expB));%max(expB)求出每一列中最大的值max(max(expB))求出整個矩陣中最大的值 expB = (expB - MIN2)/(MAX2 - MIN2);%灰度變換增強 e 發(fā)自小木蟲Android客戶端 |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 5/250 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +9 | rare12345 2026-03-18 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 【同濟軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 9/450 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|