- 相關(guān)推薦
DSP技術(shù)實(shí)習報告
一 、實(shí)習題目
1、 灰度線(xiàn)性變換
2、 灰度的對數變換
3、 鋸齒波變換
4、 圖像剪取
二 、實(shí)習目的：
1、 鞏固和深化數字圖像處理技術(shù)所涉及的數理基礎、基本算法和各種圖像處理技術(shù)方法，學(xué)習和掌握圖像變換。
2、 對學(xué)習圖像處理的基礎知識對其應用工程實(shí)踐有一定的認識，提高學(xué)生對應用軟件的使用能力。
3、 通過(guò)理論聯(lián)系實(shí)際，綜合運用所學(xué)知識，提高學(xué)生獨立分析和解決實(shí)際問(wèn)題的能力，增強學(xué)生的工程意識，打好專(zhuān)業(yè)基礎
三、實(shí)習要求:
1、能夠根據設計題目要求查閱檢索有關(guān)的文獻資料，結合題目選學(xué)有關(guān)參考書(shū);
2、熟悉計算機圖像處理的設計方法;
3、熟悉圖像灰度變換程序的設計方法;
4、掌握圖像灰度變換的仿真方法;
5、完成圖像的灰度變換。
四、系統原理描述:
灰度變換將輸入圖像映射為輸出圖像，輸出圖像每個(gè)像素點(diǎn)的灰度值僅由對應的輸入像素點(diǎn)的值決定。它常用于改變圖像的灰度范圍及分布，是圖像數字化及圖像顯示的重要工具�；叶茸儞Q因其作用性質(zhì)有時(shí)也被稱(chēng)為對比度增強、對比度拉伸或點(diǎn)運算，稱(chēng)之為灰度變換。
灰度變換實(shí)際上是灰度到灰度的映射過(guò)程。設輸入圖像為A(x，y)輸出圖像為B(x，y)，則灰度變換可表示為：
即灰度變換完全由灰度映射函數決定。顯然灰度變換不會(huì )改變圖像內像素點(diǎn)之間的空間關(guān)系。
1. 灰度線(xiàn)性變換
設圖像灰度值f(x，y)的可能值域為D，但在一定條件下，使得其取值范圍縮小為且如圖 7-11所示，這種狀態(tài)常出現于連續圖f(x，y)值的動(dòng)態(tài)范圍小，或攝影曝光不足中，因而其對應的直方圖P (D)如圖所示，集中于某個(gè)較小的灰度區間內。但是通常人們希望灰度直方圖在全部灰度區域內均勻分布，最簡(jiǎn)單地是把低反差圖像進(jìn)行灰度變換得到高反差圖像，即線(xiàn)性變換，如下式：
該式可以使變換后的灰度的取值范圍擴大到
一般來(lái)說(shuō)，只有當兩個(gè)相鄰像素的灰度值(亮度值)相差到一定程度時(shí)，人的視覺(jué)才能分辨。若灰度值D僅在較小區間內時(shí)，則人眼可分辨的亮度差的總級數則亦很少，從而造成目標圖像灰度值與背景灰度值相接近，人眼而無(wú)法分辨檢出。而對其進(jìn)行的變換后，則如上圖可使變得更大，從而使得視覺(jué)上對變換后圖像能夠分辨的亮度差的總級數增加，造成目標圖像與背景間亮度差異的加大，使原先無(wú)法被人眼檢出的目標圖像也能檢出，且變換后圖像清晰度也大大提高。
2. 灰度的對數變換
對數變換常用來(lái)擴展低值灰度，而壓縮高值灰度，這樣可以使低值灰度的圖像細節更容易看清楚，對數變換的表達式為：
其中K為一根據圖像內容來(lái)確定的常數因子。
3. 鋸齒波變換
這種方法常用來(lái)在動(dòng)態(tài)范圍較小的顯示器上顯示動(dòng)態(tài)范圍較大的圖像，或者用來(lái)發(fā)現圖像中基本輪廓，變換曲線(xiàn)如圖所示。
4. 圖像剪取
對于灰度圖像f，設定一個(gè)閾值t，若f(i，j) t，則f(i，j)的值保持不變。處理后的圖像仍然是灰度圖像，只不過(guò)閾值t以下的像素被置0。
五、方案論證及系統設計
六、系統程序的設計
1、灰度線(xiàn)性變換
#include
#include
#include "imagelib.h"
//添加圖像和視頻處理庫頭文件
#include "imagesample.h"
//添加圖像例子頭文件，
//包含原始的圖像文件goldhill[128][128]
#define MAX_PIXEL_VALUE 256
//包含原始圖像像素的最大灰度值
#define WIDTH 128
#define HEIGHT 128
//原始圖像行和列的像素數
int my_result[256];
//保存自編函數得到的灰度圖結果的數組
short histogram_output[256];
//保存調用庫函數得到的灰度圖結果的數組
short my_threshold_val=175;
//short my_threshold_val=127;
//閾值的設定
void my_histogram();
//聲明自編灰度圖算法顯示函數
void main( )
{
int i;
int my_size;
short *input;
my_size = WIDTH*HEIGHT;
//計算所有像素點(diǎn)數
input = &goldhill[0][0];
//獲取原始圖像的首地址
for(i= 0;i< MAX_PIXEL_VALUE;i++) histogram_output[i]= 0;
IMG_histogram(input，histogram_output，my_size);
//調用image.lib庫函數實(shí)現灰度圖顯示
my_histogram();
//自編C語(yǔ)言算法程序實(shí)現灰度圖顯示
IMG_threshold( &goldhill[0][0]，&goldhill[0][0]，WIDTH ，HEIGHT， my_threshold_val );
//調用image.lib庫函數實(shí)現閾值處理與顯示
while(1);
}
void my_histogram()
{
int za=10，zb=235，i，j，z1=0，z2=255，zzba;
short (*p)[128];
//定義了一個(gè)指針變量p，
//指向了一個(gè)包含128個(gè)元素的一維數組，
//即行指針變量p
p=goldhill;
zzba=(z2-z1)/(float)(zb-za);
//獲取原始圖像的首地址
// for(i=0;i<256;i++) my_result[i]=0;
for (i=0;i<128;i++)
{ for (j=0;j<128;j++)
{
if(*(*(p+i)+j)<=za)*(*(p+i)+j)=z1;
else{
if(*(*(p+i)+j)>=zb)*(*(p+i)+j)=z2;
else *(*(p+i)+j)=zzba*(*(*(p+i)+j)-za)+z1+0.5;
}
// k=*(*(p+i)+j);
// my_result[k]++;
}
}
//自編灰度直方圖算法顯示函數
}
2、灰度的對數變換
#include
#include
#include "imagelib.h"
//添加圖像和視頻處理庫頭文件
#include "imagesample.h"
//添加圖像例子頭文件，
//包含原始的圖像文件goldhill[128][128]
#define MAX_PIXEL_VALUE 256
//包含原始圖像像素的最大灰度值
#define WIDTH 128
#define HEIGHT 128
#define ZERO 0
//原始圖像行和列的像素數
int my_result[256];
//保存自編函數得到的灰度直方圖結果的數組
short histogram_output[256];
//保存調用庫函數得到的灰度直方圖結果的數組
short my_threshold_val=175;
//short my_threshold_val=127;
//閾值的設定
void my_histogram();
//聲明自編灰度直方圖算法顯示函數
void main( )
{
int i;
int my_size;
short *input;
my_size = WIDTH*HEIGHT;
//計算所有像素點(diǎn)數
input = &goldhill[0][0];
//獲取原始圖像的首地址
for(i= 0;i< MAX_PIXEL_VALUE;i++) histogram_output[i]= 0;
IMG_histogram(input，histogram_output，my_size);
//調用image.lib庫函數實(shí)現灰度直方圖顯示
my_histogram();
//自編C語(yǔ)言算法程序實(shí)現灰度直方圖顯示
IMG_threshold( &goldhill[0][0]，&goldhill[0][0]，WIDTH ，HEIGHT， my_threshold_val );
//調用image.lib庫函數實(shí)現閾值處理與顯示
while(1);
}
void my_histogram()
{
int i，j;
double za=2，zb=128，z1=0，z2=255，fij，gij，zzlba;
short (*p)[128];
//定義了一個(gè)指針變量p，
//指向了一個(gè)包含128個(gè)元素的一維數組，
//即行指針變量p
p=goldhill;
zzlba=(z2-z1)/(double)(log(zb)-log(za));
//獲取原始圖像的首地址
// for(i=0;i<256;i++) my_result[i]=0;
//初始化灰度直方圖結果數組
for (i=0;i<128;i++)
{ for (j=0;j<128;j++)
/* log transform.*/
{
fij=*(*(p+i)+j);
if(fij==0)fij=0.001;
gij=zzlba*(log(fij)-log(za))+z1;
*(*(p+i)+j)=gij+0.5;
}
}
// k=*(*(p+i)+j);
// my_result[k]++;
}
3、鋸齒波變換
#include
#include
#include "imagelib.h"
//添加圖像和視頻處理庫頭文件
#include "imagesample.h"
//添加圖像例子頭文件，
//包含原始的圖像文件goldhill[128][128]
#define MAX_PIXEL_VALUE 256
//包含原始圖像像素的最大灰度值
#define WIDTH 128
#define HEIGHT 128
//原始圖像行和列的像素數
int my_result[256];
//保存自編函數得到的灰度直方圖結果的數組
short histogram_output[256];
//保存調用庫函數得到的灰度直方圖結果的數組
short my_threshold_val=175;
//short my_threshold_val=127;
//閾值的設定
void my_histogram();
//聲明自編灰度直方圖算法顯示函數
void main( )
{
int i;
int my_size;
short *input;
my_size = WIDTH*HEIGHT;
//計算所有像素點(diǎn)數
input = &goldhill[0][0];
//獲取原始圖像的首地址
for(i= 0;i< MAX_PIXEL_VALUE;i++) histogram_output[i]= 0;
IMG_histogram(input，histogram_output，my_size);
//調用image.lib庫函數實(shí)現灰度直方圖顯示
my_histogram();
//自編C語(yǔ)言算法程序實(shí)現灰度直方圖顯示
IMG_threshold(&goldhill[0][0]，&goldhill[0][0]，WIDTH，HEIGHT， my_threshold_val );
//調用image.lib庫函數實(shí)現閾值處理與顯示
while(1);
}
void my_histogram()
{
int i，j;
float lwidth=128，ratio=128，gg=235;
short (*p)[128];
//定義了一個(gè)指針變量p，
//指向了一個(gè)包含128個(gè)元素的一維數組，
//即行指針變量p
p=goldhill;
//獲取原始圖像的首地址
// for(i=0;i<256;i++) my_result[i]=0;
//初始化灰度直方圖結果數組
for (i=0;i<128;i++)
{ for (j=0;j<128;j++)
{
gg=*(*(p+i)+j);
while(gg>lwidth) gg=gg-lwidth;
*(*(p+i)+j)=gg*ratio+0.5;
}
// k=*(*(p+i)+j);
// my_result[k]++;
}
//自編灰度直方圖算法顯示函數
}
4、圖像剪取
#include
#include
#include "imagelib.h"
//添加圖像和視頻處理庫頭文件
#include "imagesample.h"
//添加圖像例子頭文件，
//包含原始的圖像文件goldhill[128][128]
#define MAX_PIXEL_VALUE 256
//包含原始圖像像素的最大灰度值
#define WIDTH 128
#define HEIGHT 128
//原始圖像行和列的像素數
#define ZERO 0
int my_result[256];
//保存自編函數得到的灰度直方圖結果的數組
short histogram_output[256];
//保存調用庫函數得到的灰度直方圖結果的數組
short my_threshold_val=175;
//short my_threshold_val=127;
//閾值的設定
void my_histogram();
//聲明自編灰度直方圖算法顯示函數
void main( )
{
int i;
int my_size;
short *input;
my_size = WIDTH*HEIGHT;
//計算所有像素點(diǎn)數
input = &goldhill[0][0];
//獲取原始圖像的首地址
for(i= 0;i< MAX_PIXEL_VALUE;i++) histogram_output[i]= 0;
IMG_histogram(input，histogram_output，my_size);
//調用image.lib庫函數實(shí)現灰度直方圖顯示
my_histogram();
//自編C語(yǔ)言算法程序實(shí)現灰度直方圖顯示
IMG_threshold(&goldhill[0][0]，&goldhill[0][0]，WIDTH，HEIGHT，my_threshold_val );
//調用image.lib庫函數實(shí)現閾值處理與顯示
while(1);
} void my_histogram()
{
int i，j，threshold=0;
short (*p)[128];
//定義了一個(gè)指針變量p，
//指向了一個(gè)包含128個(gè)元素的一維數組，
//即行指針變量p
p=goldhill;
//獲取原始圖像的首地址
// for(i=0;i<256;i++) my_result[i]=0;
//初始化灰度結果數組
for (i=0;i<128;i++)
{ for (j=0;j<128;j++)
{
*(*(p+i)+j)=(*(*(p+i)+j)
// k=*(*(p+i)+j);
// my_result[k]++;
}
}
//自編灰度算法顯示函數
}
七、調試方法及性能分析與實(shí)驗
通過(guò)灰度線(xiàn)性變換程序將一副圖像進(jìn)行線(xiàn)性處理，線(xiàn)性變換后得到的圖像如圖(一)。通過(guò)灰度對數變換程序將圖像進(jìn)行對數處理，得到對數變換后的圖像如圖(二)。通過(guò)鋸齒波程序將圖像進(jìn)行處理 ，得到處理后的圖像如圖(三)。通過(guò)圖像剪取程序將圖像進(jìn)行處理，得到處理后的圖像如圖(四)。
圖一 線(xiàn)性變換
圖二 對數變換
圖三 鋸齒波變換
圖四 圖像剪取
八、實(shí)習總結：
基于CCS圖像處理程序調試與設計主要是在C語(yǔ)言的環(huán)境下與圖象處理的基本算法相結合。從調出圖象到編寫(xiě)C程序的圖象處理算法，最后到輸出處理后的圖象，形成一個(gè)一體化的模式，從而達到圖象處理的應用效果。通過(guò)算法和具體程序的結合，從而達到我們所要的結果。
九、實(shí)習體會(huì )：
經(jīng)過(guò)兩周的實(shí)習，我在書(shū)本中，上課時(shí)學(xué)習的知識得到了靈活運用，尤其通過(guò)對于變換前后的圖像對比，我更清楚地看到各種處理方法對于圖像的作用，本次實(shí)驗也讓我對CCS軟件有了更多的了解和認識，為進(jìn)一步學(xué)習數圖知識打下了良好的基礎。
十、參考文獻
1.TMS320C55x DSP原理及應用 汪春梅 孫洪波 編著(zhù) 電子工業(yè)出版社
2..DSP原理與應用實(shí)驗 姜陽(yáng) 周錫青 張雪英 編著(zhù) 西安電子科技大學(xué)出版社
【DSP技術(shù)實(shí)習報告】相關(guān)文章：
DSP技術(shù)的發(fā)展及應用08-01
DSP技術(shù)推動(dòng)智能物聯(lián)網(wǎng)發(fā)展09-28
基于幀的實(shí)時(shí)話(huà)音變換技術(shù)與DSP實(shí)現07-09
一種基于CAN總線(xiàn)的DSP程序加載技術(shù)06-21