dpcm.zip_DPCM_DPCM matlab_dcpm_dpcm matlab image_zip

function X_dpcm = dpcm(X,p) % DPCM Codage d'image DPCM (Differential Pulse Code Modulation) % % Utilisation : % X_DPCM = dpcm(X,predicteur) % % Variables d'entrée : % X : Image devant etre codée. % predicteur : Description du prédicteur. % % Variable de sortie : % X_DPCM : Image codée DPCM. % % % * Description du prédicteur : p = [A B C D] % ____________________ % | | | | % | A | B | C | % |______|______|______| % | | | | % | D | X | - | % |______|______|______| % | | | | % | - | - | - | % |______|______|______| % % * Exemples de prédicteurs (mode JPEG sans pertes) % % - p1 = [+1 0 0 0] A % - p2 = [0 +1 0 0] B % - p3 = [0 0 0 +1] D % - p4 = [-1 +1 0 +1] B+D-A % - p5 = [+1/2 -1/2 0 +1] D-(B-A)/2 % - p6 = [+1/2 +1 0 -1/2] B-(D-A)/2 % - p7 = [0 +1/2 0 +1/2] (B+D)/2 % % Voir aussi : IDPCM % AUTEUR : Alain Le Duff - ESEO <alain.le_duff@eseo.fr> % % VERSION : 23/04/04 T = size(X); T1 = T(1); T2 = T(2); Xx = zeros(T1+1,T2+2); X_point = Xx; Xx(2:T1+1,2:T2+1) = X; E_chap = Xx; X_chap = E_chap; for k = 2:T1+1 for l = 2:T2+1 X_point(k,l) = p(1)*X_chap(k-1,l-1) + p(2)*X_chap(k-1,l) + p(3)*X_chap(k-1,l+1) + p(4)*X_chap(k,l-1); E_chap(k,l) = round(Xx(k,l) - X_point(k,l)); X_chap(k,l) = X_point(k,l) + E_chap(k,l); end end X_dpcm = E_chap(2:T1+1,2:T2+1);