Discussion :

[nant44]

Bonjour,

Je souhaite faire interpolation/extrapolation un peu particulière d'une image et je me demande si cela est déjà fait par les spécialistes du traitement d'image.
Avant d'exposer le problème de 2D, je commence par le cas simple 1D:

Code :

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

1
2
x = 0.5:n-0.5;
f = [f(1), ..., f(n)];

je souhaite faire ceci :

Code :

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

1
2
y = 0:n;
g = [g(0),...,g(n)];

avec

Code :

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

g(k) = [f(k)+f(k+1)]/2 pour k = 1,..., n-1

et

Code :

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

1
2
[g(1)-g(0)]/[g(1)+g(0)]         =  [f(2)-f(1)]/[f(2)+f(1)]                 =>           g(0) = f(1)*[f(1)+f(2)]/[2*f(2)]
[g(n)-g(n-1)]/[g(n)+g(n-1)]   =  [f(n)-f(n-1)]/[f(n)+f(n-1)]           =>           g(n) = f(n)*[f(n)+f(n-1)]/[2*f(n-1)]

Sur un exemple numérique ceci donne:

Code :

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

1
2
f = [1, 1, 2, 5]
g = [1, 1, 1.5, 3.5, 8.75]

donc jusqu'ici c'est facile mais en passant en 2D les choses se compliquent !
Question, comment généraliser ce principe au cas 2D ??

Merci

[nant44]

Pour le moment je fait un traitement séparé sur les lignes et sur les colonnes et à la fin je fait une moyenne (voir le code ci-dessous), mais je pense que ce n'est pas optimal!

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
    function mtrxout = InterpExtrap(mtrxin)
        Arow = RowInterpExtrap(mtrxin);
        Acol = RowInterpExtrap(mtrxin.').';
        mtrxout_col = RowInterpExtrap(Arow.').';
        mtrxou_row = RowInterpExtrap(Acol);
        mtrxout = (mtrxout_col + mtrxou_row)/2;
        function grow = RowInterpExtrap(Ain)
            grow = zeros(size(Ain,1), size(Ain,2)+1);
            for row = 1:size(Ain,1)
                growtemp = 0.5*(Ain(row,1:end-1) + Ain(row,2:end));
                grow(row,:) = [Ain(row,1)*mean(Ain(row,1:2))/Ain(row,2), ...
                    growtemp,...
                    Ain(row,end)*mean(Ain(row,end-1:end))/Ain(row,end-1)];
            end
        end

Des idées ?