DSP图像处理算法的实现.ppt

DSP图像处理算法的

实现姓名：王淑阁学号：DSP图像处理的原理介绍1.离散卷积：a.公式：b.算法实现：使用一个模板mask(大小可以是3×3，5×5或7×7等)，依次覆盖在要处理图像的各个像素上，模板中心与该像素重合；模板上的各个元素均有其自身的权值，将模板各元素权值与覆盖在其下的像素值分别相乘后求和，然后进行相应的处理，最后将此所得结果赋给当前处理的像素。

5.案例分3.中值滤波中值滤波是一种非线性滤波，对像素领域内灰度的中值代替该像素的值，即：其中：

?第l阶段：直接根据需要用高级C语言实现DSP功能，测试代码的正确性。然后，移植到C6000平台，利用C6000开发环境Profile测试程序的运行时间。若不满足要求，则进入下一阶段。?第2阶段：利用C6000提供的优化方式和其他各种优化技巧，如使用：编译器选项、内联函数、字访问短型数据、软件流水、循环展开等，优化C语言代码。如果还不能满足要求，则进入第3阶段。?第3阶段：将C语言代码中耗时最长的部分抽取出来，用线性汇编语言重写，用汇编优化器进行优化。使用profile确定这段代码是否需要进一步优化。

举例：33×33相关运算的算法实现mask的设置如下：的设置如下：该mask的作用是实现图像边缘的检测。

?输入图像：一幅8位无符号图像(8行8列)，黑色部分的像素值为0，白色部分的像素值为255；对它进行相关运算，要求输出图像为32位的图像。输入图像

标准C代码的部分程序：//processingtheimagefor(irows=1;irowsaSz.height+1;irows++){for(icols=1;icols=aSz.width+1;icols++){sum=0;for(j=0;j3;j++){for(i=0;i3;i++){sum+=temp_ptr-pData[(irows+j-1)*temp_aSz.width+icols-1+i]*mask[j*3+i];}}outptr-pData[(irows-1)*aSz.width+(icols-1)]=sum;}}

利用循环展开和软件流水的程序优化方法实现：//processingtheimagefor(irows=1;irowsaSz.height+1;irows++){for(icols=1;icols=aSz.width+1;icols++){sum=0;sum+=temp_ptr-pData[(irows-1)*temp_aSz.width+icols-1]*mask[0];sum+=temp_ptr-pData[(irows-1)*temp_aSz.width+icols]*mask[1];sum+=temp_ptr-pData[(irows-1)*temp_aSz.width+icols+1]*mask[2];sum+=temp_ptr-pData[(irows)*temp_aSz.width+icols-1]*mask[3];sum+=temp_ptr-pData[(irows)*temp_aSz.width+icols]*mask[4];sum+=temp_ptr-pData[(irows)*temp_aSz.width+icols+1]*mask[5];sum+=temp_ptr-pData[(irows+1)*temp_aSz.width+icols-1]*mask[6];sum+=temp_ptr-pData[(irows+1)*temp_aSz.width+icols]*mask[7];sum+=temp_ptr-pData[(irows+1)*temp_aSz.width+icols+1]*mask[8];}outptr-pData[(irows-1)*aSz.width+(icols-1)]=sum;

?3×3相关运算处理后图像：?程序中要求输出图像为32位，但是在本次设计中，图像都是用8位来显示的，因此实验结果只是黑白图像。

验证结果的正确性以下情况，output(x,y)=0

以下情况，output(x,y)=65025

以下情况，output(x,y)=130050

以下情况，output(x,y)=195075

?通过以上的计算可知，得到的实验结果是正确的。

浦东机场钢结构合同结构特点1.采用施工总承包