|
|
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 2024 04 19, 10:21. 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 |
|
|
FS 22 Tractors
Farming Simulator 19 Mods,
FS 22 Maps,
FS22 Mods |
|
ETS2 Mods
ETS2 Trucks,
ETS2 Bus,
Euro Truck Simulator 2 Mods
|
|
FS22 Tractors
Farming Simulator 22 Mods,
FS22 Maps,
FS22 Trucks |
|
VAT calculator
VAT number check,
What is VAT,
How much is VAT |
|
Paskola internetu
Vartojimo paskola,
paskola automobiliui,
paskola būsto remontui |
|
Thermal monocular
Thermal vision camera,
Night vision ar scope,
Night vision spotting scope |
|
FS22 Mods
FS22 Harvesters,
FS22 Tractors Mods,
FS22 Maps Mods |
|
FS22 Mods
FS22 Maps,
FS22 Harvesters,
FS22 Tractors |
|
Dantų protezavimas
All on 4 implantai,
Endodontija mikroskopu,
Dantų implantacija |
|
Sims 4 Mods
Sims 4 CC Clothes,
Sims 4 Hair CC,
Sims 4 Skill Cheat |
|
Optic sight
Binoculars for hunting elk,
Best compact binoculars,
Riflescope hunting |
|
|
|