Discussion :

[defqoon]

Bonjour,

Dans le cadre de mon stage, j'analyse des cartes de variation de surface océanique. Je dois filtré spatialement les données pour retirer les signaux de grandes longueur d'onde, supérieurs à 20° de longitude par 10° de latitude.
J'ai des bases très faibles en traitement du signal. Après maintes recherches, je n'arrive pas à m'en sortir, entre les fonctions imfilter, filter2... et surtout comment créer un filtre passe haut qui corresponde à mes besoins.
Si quelqu'un à une piste, un lien qui pourrait m'éclaircir je suis preneur!

Merci

[EDIT] Je suis tombé la dessus mais je ne vois toujours pas comment m'en sortir:
http://www.cs.uregina.ca/Links/class...ab5/index.html

[FLB]

Salut,
as tu la signal processing toolbox à ta disposition?
Elle possède de nombreuses méthodes qui te simplifie la vie pour créer des filtres, je te conseil de regarder par exemple butter.
Sinon, la lecture de cet article de la faq te donnera une idée de la création de filtres simple avec les commandes de base sans toolbox :
http://matlab.developpez.com/faq/?pa...res_classiques.
N'hésite pas à reproduire les exemples pour comprendre le fonctionnement du filtrage.

[defqoon]

Merci pour tes informations mais cela ne m'a pas trop aidé! A part pour la connaissance pure bien sur. Voici mon code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
function tt = filter_SSH(data,lon_size,lat_size)
%filter SSH data
 
m=mean(data(:));
ee = data - m; %remove mean
dd = fft2(ee); %fft 2D
 
ii = fftshift(dd.*conj(dd));
uu = fftshift(dd);
 
 
s=size(ii);
hh = zeros(s(1),s(2));
 
 
for i=1:s(1)
    for j=1:s(2)
        if (i-s(1)/2)^2/lon_size+(j-s(2)/2)^2/lat_size < 1 %remove all the low frequency
            hh(i,j)=uu(i,j);
            ii(i,j)=0;
%             uu(i,j)=0;
        end
    end
end
 
% figure
% surface(1:s(2),1:s(1),hh)
% hold on
% plot(s(2)/2,s(1)/2,'kx')
 
D=hh;
 
% figure
tt = ifft2(ifftshift(D),'symmetric');

une idée ?