开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > matlab例程 > 查看源码
main44gray.m
资源名称：MPEG-2(MATLAB).rar　[点击查看]
上传用户：cxsjwj
上传日期：2022-08-09
资源大小：34k
文件大小：6k
源码类别：

matlab例程

开发平台：

Matlab

main44gray.m：源码内容
  2. close all
  3. clc
  4. tic
  6. gray
  7. prinf(resize,long,wide)
  8. startframe
  9. endframe
  11. aviobj = avifile('2.avi','fps',25,'COMPRESSION','None');
  12. %motion Vectors
  13. motionVect=zeros(2,256);
  14. fid=fopen('Qtable2.txt','r');
  15. %array is quantization matrix
  16. array=fscanf(fid,'%e',[8,inf]);
  17. %array=2*array;
  18. array=array;
  20. i = 0;
  21. %flag to tell if B frame was predicted from I or P frame flag=1 P frame
  22. %flag=0 I frame
  23. flag=2;
  24. %GOP is of 10 pictures IBBPBBPBBI
  25. mov=aviread('news.avi');
  26. movinfo=aviinfo('news.avi');
  27. noframe=movinfo.NumFrames;
  28. %mbSize macroblock size used for motion estimation
  29. mbSize = 8;
  30. %p for search area
  31. p = 7;
  34. global I1d P4d P7d I10d B2d B3d B5d B6d B8d B9d
  35. for i=startframe:10:endframe
  36.     i
  37. framedata=aviread('news.avi',i);
  38. I1=frame2im(framedata);
  39. if resize==1
  40.     I1=imresize(I1,[long,wide]);
  41.     end
  42. I1=rgb2gray(I1);
  44. I1=double(I1);
  45. JQ=forwardDCT(I1,array);
  46. bufferI1 = inverseDCT(JQ,array);
  47. c= makeLayers(JQ);
  49. %%calling decoder
  50. identifier=1;
  51. decoder(c,motionVect,identifier,flag);
  52. %%%this loop is for accessing frame from actual movie
  54.     if i~=1
  55.         j=i-round(i/10);
  56.         [B2,B3,P4,B5,B6,P7,B8,B9,I10]=GOP(j,resize,long,wide);
  57.     elseif i==1
  58.         j=i;
  59.         [B2,B3,P4,B5,B6,P7,B8,B9,I10]=GOP(j,resize,long,wide);
  60.     end
  61.         %processing on P4 frame
  62.         [bufferP4,streamP4,motionVect]=compensatedFrame(P4,bufferI1,mbSize,p,array); 
  63.         c= makeLayers(streamP4);
  65.         %calling decoder
  66.         identifier=4;
  67.         %decoder(diffDC,valVectorBase,lenVectorBase,valVectorEnhan,lenVectorEnhan,identifier,flag,motionVect);
  68.         decoder(c,motionVect,identifier,flag);
  72. % % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%        
  73.         
  74.         %for B2    
  75.         [B2t ,motionVect, flag] = bFrameProc(B2,bufferI1,bufferP4,0.5*array,mbSize,p);
  76.         c= makeLayers(B2t);
  77.         
  78.         %calling decoder
  79.         identifier=2;
  80.         decoder(c,motionVect,identifier,flag);
  81.         
  82.         %for B3    
  83.         [B3t ,motionVect ,flag] = bFrameProc(B3,bufferI1,bufferP4,0.5*array,mbSize,p);
  84.         c= makeLayers(B3t);
  85.         
  86.         %calling decoder
  87.         identifier=3;
  88.         decoder(c,motionVect,identifier,flag);
  89.  
  90.         %predicting p7 from p4
  91.         [bufferP7,streamP7,motionVect]=compensatedFrame(P7,bufferP4,mbSize,p,array);
  92.         c= makeLayers(streamP7);
  93.         
  94.         %calling decoder
  95.         identifier=7;
  96.         decoder(c,motionVect,identifier,flag);
  97. %          
  98.         %for B5    
  99.         [B5t motionVect flag] = bFrameProc(B5,bufferP4,bufferP7,0.5*array,mbSize,p);
  100.         c= makeLayers(B5t);
  101.         
  102.         %calling decoder
  103.         identifier=5;
  104.         decoder(c,motionVect,identifier,flag);
  105.    
  106.         %for B6    
  107.         [B6t motionVect flag] = bFrameProc(B6,bufferP4,bufferP7,0.5*array,mbSize,p);
  108.         c= makeLayers(B6t);
  109.         
  110.         %calling decoder
  111.         identifier=6;
  112.         decoder(c,motionVect,identifier,flag);
  113.  
  114.         %processing on I10 frame
  115.         streamI10=forwardDCT(I10,array);
  116.         bufferI10 = inverseDCT(streamI10,array);
  117.         c= makeLayers(streamI10);
  118.         %%calling decoder
  119.         identifier=10;
  120.         decoder(c,motionVect,identifier,flag);
  121. %         
  122.         %for B8    
  123.         [B8t motionVect flag] = bFrameProc(B8,bufferP7,bufferI10,0.5*array,mbSize,p);
  124.         c= makeLayers(B8t);
  125.         
  126.         %calling decoder
  127.         identifier=8;
  128.         decoder(c,motionVect,identifier,flag);
  130.         %for B9    
  131.         [B9t motionVect flag] = bFrameProc(B9,bufferP7,bufferI10,0.5*array,mbSize,p);
  132.         c= makeLayers(B9t);
  133.         
  134.         %calling decoder
  135.         identifier=9;
  136.         decoder(c,motionVect,identifier,flag);
  137.     
  138. %         
  139. %         
  140.         
  141. % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  145.         
  146.     
  148.         if i==1
  149.                 imshow(uint8(I1d))
  150.                 aviobj = addframe(aviobj,getframe);
  151.                 clf;
  152.                 
  153.                 imshow(uint8(B2d))
  154.                 aviobj = addframe(aviobj,getframe);
  155.                 clf;
  156.                 
  157.                 imshow(uint8(B3d))
  158.                 aviobj = addframe(aviobj,getframe);
  159.                 clf;
  160.                 
  161.                 imshow(uint8(P4d))
  162.                 aviobj = addframe(aviobj,getframe);
  163.                 clf;
  164.                 
  165.                 imshow(uint8(B5d))
  166.                 aviobj = addframe(aviobj,getframe);
  167.                 clf;
  168.                 
  169.                 imshow(uint8(B6d))
  170.                 aviobj = addframe(aviobj,getframe);
  171.                 clf;
  172.                 
  173.                 imshow(uint8(P7d))
  174.                 aviobj = addframe(aviobj,getframe);
  175.                 clf;
  176.                 
  177.                 imshow(uint8(B8d))
  178.                 aviobj = addframe(aviobj,getframe);
  179.                 clf;
  180.                 
  181.                 imshow(uint8(B9d))
  182.                 aviobj = addframe(aviobj,getframe);
  183.                 clf;
  184.                 
  185.                 imshow(uint8(I10d))
  186.                 aviobj = addframe(aviobj,getframe); 
  187.                 clf;
  188.             else
  190.                 
  191.                 imshow(uint8(B2d))
  192.                 aviobj = addframe(aviobj,getframe);
  193.                 clf;
  194.                 
  195.                 imshow(uint8(B3d))
  196.                 aviobj = addframe(aviobj,getframe);
  197.                 clf;
  198.                 
  199.                 imshow(uint8(P4d))
  200.                 aviobj = addframe(aviobj,getframe);
  201.                 clf;
  202.                 
  203.                 imshow(uint8(B5d))
  204.                 aviobj = addframe(aviobj,getframe);
  205.                 clf;
  206.                 
  207.                 imshow(uint8(B6d))
  208.                 aviobj = addframe(aviobj,getframe);
  209.                 clf;
  210.                 
  211.                 imshow(uint8(P7d))
  212.                 aviobj = addframe(aviobj,getframe);
  213.                 clf;
  214.                 
  215.                 imshow(uint8(B8d))
  216.                 aviobj = addframe(aviobj,getframe);
  217.                 clf;
  218.                 
  219.                 imshow(uint8(B9d))
  220.                 aviobj = addframe(aviobj,getframe);
  221.                 clf;
  222.                 
  223.                 imshow(uint8(I10d))
  224.                 aviobj = addframe(aviobj,getframe);
  225.                 clf;
  226.             end
  228.            
  230. end
  231. aviobj = close(aviobj);
  232. toc