[CichyGlosZTyluGlowy]

#matlab #matematyka Mam taki skrypt http://pastebin.com/DDkgjW5p który ma za zadanie wyszukać krawędzie na obrazku. Ignoruję w nim skrajne piksele, ale żeby rozmiar grafiki się zgadzał jako "podkład" do edytowania kopiuję czarno-białe zdjęcie
Masked = Gray;
Skrypt zdaje się jako-tako działać, jeśli właśnie to jest w podkładzie, ale kiedy zmieniam tę linijkę na zeros(dim), GrayFFT, albo ones(dim), dostaję w wyniku lekko zaciemniony oryginał. Dlaczego?

[CichyGlosZTyluGlowy]

@CichyGlosZTyluGlowy: wersja z GrayFFT jako podkładem

[CichyGlosZTyluGlowy]

@CichyGlosZTyluGlowy: i zerami

[E.....L]

@CichyGlosZTyluGlowy: a ta maska do konwolucji na pewno jest dobra?

[CichyGlosZTyluGlowy]

@EssePL: wziąłem ją z jakiegoś skryptu w necie http://gdr.geekhood.net/gdrwpl/heavy/studia/Lecture2007%238__Splot_Korelacja.pdf

[E.....L]

@CichyGlosZTyluGlowy: jak wrócę do domu mogę na to zerknąć - ta maska mi niespecjalnie się podoba, a nie mam teraz ani jak zerknąć do książki, ani jak matlaba odpalić ;)

[CichyGlosZTyluGlowy]

@EssePL: byłbym wdzięczny :)

[prusi]

maska jest ok.
@CichyGlosZTyluGlowy: nie za bardzo rozumiem co robi ta linijka:

Masked(row,col) = sum(sum(sample));

[CichyGlosZTyluGlowy]

@prusi: sumuje wszystkie elementy macierzy- pojedyncze sum daje wektor sum kolumn.

[prusi]

@CichyGlosZTyluGlowy: no ok, ale po co? :)

[CichyGlosZTyluGlowy]

@prusi: to ma być suma ważona pikseli otaczających ten na którym akurat pracuję. Opierałem się na tym slajdzie:

[CichyGlosZTyluGlowy]

@CichyGlosZTyluGlowy: no chyba, że źle to zrozumiałem...

[prusi]

@CichyGlosZTyluGlowy: dobrze rozumiesz, ale to jest bez użycia ftt, jak wywalisz transformację na ftt i z powrotem to w ten sposób możesz to zrobić. z ftt splot wygląda trochę inaczej.

[CichyGlosZTyluGlowy]

@prusi: a jak by to wyglądało z fft?

[prusi]

@CichyGlosZTyluGlowy: zaczekaj 5 min, zaraz sprawdzę u siebie czy sam nic nie poknociłem.

[CichyGlosZTyluGlowy]

@prusi: faktycznie bez fft działa. Tylko dlaczego facet opowiadał nam o tym na wykładzie o fft?

[prusi]

@CichyGlosZTyluGlowy: bo można splot realizować w dziedzinie częstotliwości, wtedy to jest zwykłe mnożenie (bez sumowania), czyli maskę i próbkę poddajesz ftt, mnożysz i masz wynik, który poddajesz odwrotnej transformacji. w ten sposób też powinno działać. powinno się chyba jeszcze realizować jakieś nakładanie się tych sampli, ale nie wiem czy przy wykrywaniu krawędzi to ma jakikolwiek sens :)

[CichyGlosZTyluGlowy]

@prusi: do takiego splotu Maska i Obrazek muszą być chyba uzupełnione zerami, żeby miały wymiary sumy ich długości-1?

[CichyGlosZTyluGlowy]

@CichyGlosZTyluGlowy: a nie, dobra nie doczytałem, że chodzi o próbkę :D

[prusi]

@CichyGlosZTyluGlowy: właściwie można bez próbkowania, też powinno działać. popróbuj i daj znać czy się udało, bo Octave mi wariuje.