Filtro Geométrico Móvel Médio

Eu preciso testar algumas técnicas básicas de processamento de imagem no Matlab. Eu preciso testar e comparar especialmente dois tipos de filtros: filtro médio e filtro mediano. Para alisar a imagem usando a mediana de filtragem, existe uma grande função medfilt2 da caixa de ferramentas de processamento de imagem. Existe alguma função semelhante para filtro médio Ou como usar a função filter2 para criar o filtro médio Uma das coisas mais importantes para mim é ter a possibilidade de configurar o raio do filtro. I. e. Para o filtro mediano, se eu quiser o raio de 3 x 3 (máscara), eu uso apenas eu gostaria de conseguir algo semelhante para o filtro médio. Perguntei Nov 15 09 at 16:12 Eu vejo boas respostas já foram dadas, mas eu pensei que poderia ser bom para dar apenas uma maneira de realizar filtragem média no MATLAB sem usar funções especiais ou caixas de ferramentas. Isso também é muito bom para entender exatamente como o processo funciona, uma vez que você é obrigado a definir explicitamente o kernel de convolução. O kernel médio do filtro é, felizmente, muito fácil: Note que para imagens em cores você teria que aplicar isso a cada um dos canais na imagem. User8264: Eu não tenho acesso ao livro agora, mas geralmente o kernel gaussiano fornece um efeito de suavização mais suave e tende a preservar as bordas melhor do que um filtro médio do mesmo tamanho. Pense na resposta de freqüência do filtro passa-baixo em ambos os casos. Aqui está uma página com uma boa explicação: homepages. inf. ed. ac. ukrbfHIPR2gsmooth. htm ndash Amro Aug 1 14 às 9: 48Eu preciso projetar um filtro de média móvel que tem uma freqüência de corte de 7,8 Hz. Eu usei filtros de média móvel antes, mas até onde eu estou ciente, o único parâmetro que pode ser alimentado é o número de pontos a serem calculados. Como isso pode se relacionar com uma freqüência de corte O inverso de 7,8 Hz é de 130 ms, e Im trabalhando com dados que são amostrados a 1000 Hz. Isso implica que eu deveria estar usando um tamanho de janela de filtro média móvel de 130 amostras, ou há algo mais que estou faltando aqui pediu Jul 18 13 at 9:52 O filtro de média móvel é o filtro usado no domínio do tempo para remover O ruído adicionado e também para o propósito de suavização, mas se você usar o mesmo filtro de média móvel no domínio da freqüência para a separação de freqüência, o desempenho será pior. Então, nesse caso, use filtros de domínio de freqüência O filtro de média móvel (por vezes conhecido coloquialmente como um filtro de caixa) tem uma resposta de impulso retangular: Ou, declarado de forma diferente: Lembrando que uma resposta em freqüência de sistemas de tempo discreto É igual à transformada de Fourier de tempo discreto da sua resposta de impulso, podemos calculá-la da seguinte forma: O que mais interessou para o seu caso é a resposta de magnitude do filtro, H (ômega). Usando algumas manipulações simples, podemos obter isso em uma forma mais fácil de compreender: Isso pode não parecer mais fácil de entender. No entanto, devido à identidade Eulers. Lembre-se que: Portanto, podemos escrever o acima como: Como eu disse antes, o que você está realmente preocupado com a magnitude da resposta de freqüência. Assim, podemos tomar a magnitude do acima para simplificá-lo ainda mais: Nota: Nós somos capazes de soltar os termos exponenciais, porque eles não afetam a magnitude do resultado e 1 para todos os valores de ômega. Como xy xy para quaisquer dois números finitos x e y, podemos concluir que a presença dos termos exponenciais não afeta a resposta da magnitude global (em vez disso, eles afetam a resposta da fase do sistema). A função resultante dentro dos parênteses de magnitude é uma forma de um kernel de Dirichlet. É chamado às vezes uma função periódica de sinc, porque se assemelha à função do sinc um tanto na aparência, mas é periódica preferivelmente. De qualquer forma, uma vez que a definição de freqüência de corte é um pouco underspecified (-3 dB ponto -6 dB ponto primeiro sidelobe nulo), você pode usar a equação acima para resolver o que você precisa. Especificamente, você pode fazer o seguinte: Definir H (omega) para o valor correspondente à resposta do filtro que você deseja na freqüência de corte. Defina ômega igual à freqüência de corte. Para mapear uma freqüência de tempo contínuo para o domínio de tempo discreto, lembre-se que omega 2pi frac, onde fs é sua taxa de amostragem. Encontre o valor de N que lhe dá o melhor acordo entre os lados esquerdo e direito da equação. Isso deve ser o comprimento de sua média móvel. Se N é o comprimento da média móvel, então uma frequência de corte aproximada F (válida para N gt 2) na frequência normalizada Fffs é: O inverso disso é Esta fórmula é assintoticamente correta para N grande e tem cerca de 2 erro Para N2, e menos de 0,5 para N4. P. S. Depois de dois anos, aqui finalmente qual foi a abordagem seguida. O resultado foi baseado na aproximação do espectro de amplitude da MA em torno de f0 como uma parábola (série de 2ª ordem) de acordo com MA (Omega) aproximadamente 1 (frac-fra) Omega2 que pode ser feita mais exata perto do cruzamento zero de MA (Omega) Frac por multiplicação de Omega por um coeficiente de obtenção de MA (Omega) aprox. 10.907523 (frac - frac) Omega2 A solução de MA (Omega) - frac 0 dá os resultados acima, onde 2pi F Omega. Tudo o que acima se refere à freqüência de corte -3dB, o sujeito deste post. Às vezes, porém, é interessante obter um perfil de atenuação em banda de parada que é comparável ao de um filtro passa-baixo IIR de 1ª ordem (LPF de um pólo) com uma determinada freqüência de corte -3dB (tal LPF é também chamado integrador com vazamento, Tendo um pólo não exatamente em DC, mas próximo a ele). De facto, tanto a MA como a Ia ordem IIR LPF têm uma inclinação de 20dBdecade na banda de paragem (é necessário um N maior do que o utilizado na figura, N32, para ver isto), mas enquanto MA tem nulos espectricos em FkN e um 1f evelope, o filtro IIR só tem um perfil 1f. Se se deseja obter um filtro MA com capacidades semelhantes de filtragem de ruído como este filtro IIR, e corresponder às frequências de corte 3dB para ser o mesmo, ao comparar os dois espectros, ele perceberá que a ondulação da banda de parada do filtro MA acaba 3dB abaixo do filtro IIR. Para obter a mesma ondulação de banda de parada (ou seja, a mesma atenuação de potência de ruído) como o filtro IIR as fórmulas podem ser modificadas da seguinte forma: Eu encontrei de volta o script Mathematica onde eu calculou o corte para vários filtros, incluindo o MA. O resultado foi baseado na aproximação do espectro MA em torno de f0 como uma parábola de acordo com MA (Omega) Sin (OmegaN2) Sin (Omega2) Omega 2piF MA (F) aproximadamente N16F2 (N-N3) pi2. E derivando o cruzamento com 1sqrt de lá. Ndash Massimo Jan 17 16 às 2: 08Qual é a diferença entre a média móvel ea média móvel ponderada A média móvel de 5 períodos, com base nos preços acima, seria calculada usando a seguinte fórmula: Com base na equação acima, o preço médio sobre a média móvel Período mencionado acima foi de 90,66. Usando médias móveis é um método eficaz para eliminar flutuações de preços fortes. A principal limitação é que os pontos de dados de dados mais antigos não são ponderados de forma diferente dos pontos de dados próximos ao início do conjunto de dados. É aqui que as médias móveis ponderadas entram em jogo. As médias ponderadas atribuem uma ponderação mais pesada a pontos de dados mais atuais, uma vez que são mais relevantes do que pontos de dados no passado distante. A soma da ponderação deve somar 1 (ou 100). No caso da média móvel simples, as ponderações são distribuídas igualmente, razão pela qual não são mostradas na tabela acima. Preço de Fechamento da AAPL