interpolation et extrapolation d'image

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:

---

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

---

je souhaite faire ceci :

Code:

---

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

---

avec

Code:

---

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

---

et

Code:

---

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:

---

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

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:

---

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 ?