 |
|
|
Elektronika.lt portalo forumas
Jūs esate neprisijungęs lankytojas. Norint dalyvauti diskusijose, būtina užsiregistruoti ir prisijungti prie forumo.
Prisijungę galėsite kurti naujas temas, atsakyti į kitų užduotus klausimus, balsuoti forumo apklausose.
Administracija pasilieka teisę pašalinti pasisakymus bei dalyvius,
kurie nesilaiko forumo taisyklių.
Pastebėjus nusižengimus, prašome pranešti.
Dabar yra 2025 05 13, 10:43. Visos datos yra GMT + 2 valandos.
|
|
|
 |
Forumas » Programavimas » Padekit su matlab
|
Jūs negalite rašyti naujų pranešimų į šį forumą
Jūs negalite atsakinėti į pranešimus šiame forume
Jūs negalite redaguoti savo pranešimų šiame forume
Jūs negalite ištrinti savo pranešimų šiame forume
Jūs negalite dalyvauti apklausose šiame forume
|
|
|
|
|
 |
Padekit su matlab |
 Parašytas: 2009 05 22, 17:19 |
  |
|
|
|
Sveiki, su matlab bandau atsatyt paveiksliuka tik is fazes, jos informacija saugodamas kaip 1bita. po 100 iteraciju pradinis paveiksliukas turi atsistatyti, taciau ne pas mane:( gal galit pasakyti kur yra problema, busiu labai dekingas! Stai kodas: (paveiksliuka patikrinimui isikelkit koki norit, bet 64x64 didumo tik)
pav = imread('logo.jpg');
pav=rgb2gray(pav);
l=256;
[m,n]=size(pav);
len=m*n;
image2=reshape(pav,len,1);
h = hist(image2,[0:l-1]);
%h=imhist(pav);
save h;
xform=[
1 0 0
0 1 0
64 64 1];
tform_translate=maketform('affine',xform);
pav2 = imtransform(pav, tform_translate,...
'XData', [1 (size(pav,2)+xform(3,1))],...
'YData', [1 (size(pav,1)+xform(3,2))]);
figure, imshow(pav2)
pavfft=fft2(pav2);
phi=angle(pavfft);
ampl=abs(pavfft);
ampl( =4760;
phi(phi<0)=0;
phi(phi>0)=1;
save phi;
phi2=phi;
phi2(phi==0)=-pi/2;
phi2(phi==1)=pi/2;
ipav= ampl.*exp(phi2*sqrt(-1));
ipavfft=uint8(abs(ifft2(ipav)));
his_n=ipavfft(65:128,65:128);
his=histeq(his_n,h);
xform=[
1 0 0
0 1 0
64 64 1];
tform_translate=maketform('affine',xform);
his2 = imtransform(his, tform_translate,...
'XData', [1 (size(his,2)+xform(3,1))],...
'YData', [1 (size(his,1)+xform(3,2))]);
his3=his2;
figure, imshow(his2)
iter=0;
while iter<10
pavfft2=fft2(his3);
phi3=angle(pavfft2);
ampl(=5500;
[M,N] = size(phi3);
for i = 1:1:M,
for j = 1:1:N,
if phi3(i,j)>0 & phi(i,j)==0
phi3(i,j) = -pi/2 ; end
if phi3(i,j)<0 & phi(i,j)==1
phi3(i,j)=pi/2;end
end
end
ipavfft2=uint8(abs(ifft2(ampl.*exp(phi3*sqrt(-1)))));
his_n2=ipavfft2(65:128,65:128);
his2=histeq(his_n2,h);
xform=[
1 0 0
0 1 0
64 64 1];
tform_translate=maketform('affine',xform);
his3 = imtransform(his2, tform_translate,...
'XData', [1 (size(his2,2)+xform(3,1))],...
'YData', [1 (size(his2,1)+xform(3,2))]);
iter=iter+1;
end
figure, imshow(his3); |
|
|
|
|
|
|
|
Padekit su matlab |
| Parašytas: 2009 05 26, 16:24 |
|
|
|
|
| arba bent jau zinote kaip atsatyti paveiksliuka kai amplitude yra konstanta, prilyginta kokiam nors skaiciui, tarkim visos matricos vertes lygios 1000, nes butent kai amplitude prilyginu kokiam skaiciui, nebeastato paveiksliuko, busiu labai dekingas! |
|
|
|
|
|
|
|
Padekit su matlab |
| Parašytas: 2009 06 02, 01:48 |
|
|
|
|
pats po truputi arteju prie tikslo, bet gal dar pasiziuresit ir suprasit kame likusi problema, nes dar nepilnai veikia kazkas cia, paziurekit jei galit:(
pav = imread('logo.jpg');
pave=rgb2gray(pav);
%h=imhist(pave,64);
h=imhist(pave);
%pav=double(pave);
save h;
xform=[
1 0 0
0 1 0
64 64 1];
tform_translate=maketform('affine',xform);
pav2 = imtransform(pave, tform_translate,...
'XData', [1 (size(pave,2)+xform(3,1))],...
'YData', [1 (size(pave,1)+xform(3,2))]);
figure, imshow(pave)
pav2=double(pav2);
pavfft=fft2(pav2);
phi=angle(pavfft);
ampl=abs(pavfft);
ampl(:)=5000;
phi(phi<0)=0;
phi(phi>0)=1;
save phi;
phi2=phi;
phi2(phi==0)=-pi/2;
phi2(phi==1)=pi/2;
pavi= ampl.*exp(1i*phi2);
pavifft=uint8(abs(ifft2(pavi)));
his_n=pavifft(65:128,65:128);
his=histeq(his_n,h);
tform_translate=maketform('affine',xform);
his2 = imtransform(his, tform_translate,...
'XData', [1 (size(his,2)+xform(3,1))],...
'YData', [1 (size(his,1)+xform(3,2))]);
his3=his2;
%paviks=fft2(his2);
%ampl2=abs(paviks);
figure, imshow(his2)
for iter=1:100
pavfft2=fft2(his3);
phi3=angle(pavfft2);
%ampl(:)=800;
ampl2=abs(pavfft2);
[M,N] = size(phi3);
for k = 1:1:M,
for l = 1:1:N,
if phi3(k,l)>0 && phi(k,l)==0
phi3(k,l) = -pi/2 ; end
if phi3(k,l)<0 && phi(k,l)==1
phi3(k,l)=pi/2; end
end
end
pavi2=ampl2.*exp(1i*phi3);
pavifft2=uint8(abs(ifft2(pavi2)));
his_n2=pavifft2(65:128,65:128);
his2=histeq(his_n2,h);
tform_translate=maketform('affine',xform);
his3 = imtransform(his2, tform_translate,...
'XData', [1 (size(his2,2)+xform(3,1))],...
'YData', [1 (size(his2,1)+xform(3,2))]);
iter=iter+1;
end
figure, imshow(his3); |
|
|
|
|
|
|
  |
|
|
|
 |
Google paieška forume |
|
|
|
Naujos temos forume |
|
|
FS25 Tractors
Farming Simulator 25 Mods,
FS25 Maps,
FS25 Trucks |
|
ETS2 Mods
ETS2 Trucks,
ETS2 Bus,
Euro Truck Simulator 2 Mods
|
|
FS22 Tractors
Farming Simulator 22 Mods,
FS22 Maps,
FS25 Mods |
|
VAT calculator
VAT number check,
What is VAT,
How much is VAT |
|
FS25 Mods
FS25 Harvesters,
FS25 Tractors Mods,
FS25 Maps Mods |
|
Dantų protezavimas
All on 4 implantai,
Endodontija mikroskopu,
Dantų implantacija |
|
FS25 Mods
FS25 Maps,
FS25 Cheats,
FS25 Install Mods |
|
FS25 Mods
Farming Simulator 25 Mods,
FS25 Maps |
|
ATS Trailers
American Truck Simulator Mods,
ATS Trucks,
ATS Maps |
|
Football Training Kit
Football Training Equipment,
Football Skills,
Football Training |
|
|
 |
