最近著手把CSK移植到DSP中沽瞭，先看一些DSP中圖像處理的一些例子，第一件事當然就是怎么把圖像數(shù)據(jù)倒入CCS工程中了豆同，去年倒是用過一點CCS贷帮，再拿起來已經(jīng)忘得差不多了，這篇文章主要記錄一些學習的過程：

目錄

1. DSP導入圖像數(shù)據(jù)
2. 窗函數(shù)和加窗
3. 定點數(shù)和浮點數(shù)
4. 圖像的FFT運算
5. 圖像的IFFT運算

一. DSP導入圖像數(shù)據(jù)

搞了一下午大概可以了诱告，主要是如何導入數(shù)據(jù)撵枢，如何利用CCS的Image analyzer來做顯示民晒。

1. 準備dat形式的數(shù)據(jù)

CCS導入數(shù)據(jù)是dat格式的，利用matlab把數(shù)據(jù)保存到dat文件里锄禽，這個直接找的代碼：

clear
clc
I=imread('img.bmp');       //這里把圖片都進來就行
[s1,s2]=size(I);
Num=s1*s2/4;
fid=fopen('img.dat','w');      //“寫”形式打開
fprintf(fid,'1651 1 80000000 0 %X\n',s1*s2/4);   //信息頭
for i=1:4:(s1*s2-3)
    fprintf(fid,'0X% 02X% 02X% 02X% 02X\n',double(I(i+3:-1:i)));
end
fclose(fid);    //關閉

不用去糾結得到的dat的格式規(guī)則潜必，要用到的話再細究。

2. 準備內(nèi)存沃但，導入數(shù)據(jù)

load memory

load memory

memory browser

查看圖像灼伤。

參數(shù)設置

image

替代方案

最后也沒有發(fā)現(xiàn)到底是什么樣的原因颖侄，但是在李老師（李翔宇）的指導下，把圖片數(shù)據(jù)的首地址給定隆敢，最后給在的是0x00850000发皿，這樣就可以了崔慧，老師說也 不知道什么原因拂蝎，先放下吧，最起碼這樣可以進行一些簡單的算法仿真了惶室。
具體這樣：
unsigned char *img=(unsigned char *)0x0850000;
再load memory就沒有問題了温自，也可以進行操作。

二.窗函數(shù)的實現(xiàn)和導入

CSK的實現(xiàn)過程中要用到兩種窗函數(shù)皇钞，分別是高斯和漢明悼泌，這兩種窗函數(shù)可以利用matlab提前生成好，然后作為頭文件來導入到CCS工程中夹界。這個實現(xiàn)起來也不難馆里。

matlab代碼

%得到漢明窗和高斯窗的c代碼，不能直接寫入h文件，就先寫入txt再復制過去了鸠踪，主要是要中間的逗號丙者。
fid=fopen('cos.txt','wt');
cos_window = hann(64) * hann(64)';
cos_window=cos_window(:);
fprintf(fid,'cos_window_64_64={');
for i=1:size(cos_window)
    fprintf(fid,'%f,',cos_window(i));
end
fclose(fid);

這樣生成的是txt文件，主要是包含逗號就可以了营密，稍加修改就可以復制到CCS建立的頭文件中了械媒，數(shù)據(jù)類型選擇float類型。我試著導入了一個漢明窗评汰，64*64的纷捞，拉成一維數(shù)組后導入是這個樣子：

漢明窗

CCS中如果想看二維圖像長什么樣子的話還得把float轉(zhuǎn)換成8位int，挺麻煩的被去，所以就這樣看看吧主儡。

然后我們把窗函數(shù)加在圖像上看一看長什么樣子，我這里換了一副圖像编振，所有的操作都是針對6464的圖像來的缀辩。
導入的圖像為

原圖

img_with_window

for(i=0;i<4096;i++)
    {
        *(img_smooth++)=*(img_with_window++);  //轉(zhuǎn)換為uchar
    }

會莫名奇妙出現(xiàn)問題畅蹂，不僅img_smooth的結果不對健无，還會把原先img_with_window的數(shù)據(jù)給沖掉。數(shù)組地址我都是手動給定的液斜，保證不會沖突的累贤，也沒找到什么原因。（如有人知道還望指點Ｉ倨帷）

img_smooth

img_smooth就變成這樣了臼膏，這肯定是不對的。而且img_with_window就全變成0了示损。
后來換了一種方法用下標來做渗磅，就沒有問題了，真是神奇检访。

image_smooth

image_smooth

都是沒有問題的脆贵，這個smooth的命名的原因就是因為上次寫均值濾波時用的那個變量医清，還沒有改。

三.定點數(shù)和浮點數(shù)的區(qū)別

PC編程很少遇到這么細節(jié)的問題卖氨，但是DSP上就不同了会烙，以前只知道定點數(shù)需要定標负懦，浮點數(shù)是采用類似于科學計數(shù)法的一種方法，具體的細節(jié)就不清楚了柏腻，DSP還有定點和浮點之分密似，所以把這里的細節(jié)看了看，然后整理下葫盼。

定點數(shù)

所謂定點數(shù)就是約定小數(shù)點的位置是不變的残腌，這樣通常將定點數(shù)據(jù)表示成純小數(shù)或者純整數(shù)，純小數(shù)的話贫导，就把小數(shù)點固定在數(shù)值部分的最高位之前抛猫；純小數(shù)的話，則把小數(shù)點固定在數(shù)值部分的最后面孩灯，如圖：

小數(shù)點在機器中是不表示出來的闺金，是一個提前約定好的位置，對于一臺計算機峰档，一旦確定了小數(shù)點的位置败匹，就不再改變。
必須明確的一點就是計算機是無法計算小數(shù)的讥巡，所有的小數(shù)計算都是通整數(shù)計算完成的掀亩，這就導致事先約定小數(shù)點的位置尤為重要，這就是定標欢顷，在定點運算中槽棍，定標很重要。
為了簡單抬驴，我們用16位的來說吧炼七，也就是說參與運算的數(shù)都是16位的整型數(shù)，那么如何處理小數(shù)的關鍵就在于如何確定一個數(shù)據(jù)小數(shù)點的位置布持，這樣的話數(shù)據(jù)的精度和數(shù)據(jù)的范圍就是一對矛盾了豌拙，精度越高，范圍則會越小题暖。
數(shù)的標定有Q和S兩種表示方法：
對于16位的來說按傅，最高位為符號保留位，那么根據(jù)定標的不同芙委，有Q0到Q15共16中定標方式逞敷，精度依次增高狂秦，范圍依次降低灌侣。
例：
0010 0000 0000 0000b 表示8192 Q0定標
而同樣的一個數(shù)：Q15定標
0010 0000 0000 0000b 表示0.25
這個圖片不是特別清楚，可以看下：

這個算法是直接給出的蛾魄，也很簡單，實際上是一個移位（乘法就相當于移位）湿滓，注意調(diào)用之前x應該是歸一化到[-1,1)之間的滴须。

void fltoq15(float x[], short r[], short nx)
{
int i, a;
for(i = 0; i < nx; i++)
{
a = 32768 * x[i];
// saturate to 16.bit //
if (a>32767) a = 32767;
if (a<.32768) a = .32768;
r[i] = (short) a;
}
}

至于浮點數(shù)怎么判斷溢出和運算規(guī)則就暫且不細究了，我這次的應用只涉及到所用的算法是要求Q15格式叽奥，所以需要手動轉(zhuǎn)換扔水。

浮點數(shù)

正是由于定點數(shù)的表示方法太過僵硬，固定的小數(shù)點位置決定了固定位數(shù)的整數(shù)部分和小數(shù)部分朝氓，不利于表示特別大和特別小的數(shù)魔市，現(xiàn)在絕大多數(shù)的現(xiàn)代計算機系統(tǒng)都采納了浮點數(shù)表達實數(shù)。
浮點數(shù)用一個尾數(shù)赵哲，一個基數(shù)待德，一個指數(shù)來表示一個數(shù)。
比如123.45用十進制的科學技術法可以表示為：1.2345*10^2枫夺。其中1.2345就是尾數(shù)磅网，10是基數(shù)，2是指數(shù)筷屡，這樣可以更加靈活得表示更大范圍的實數(shù)涧偷。

IEEE浮點數(shù)

IEEE標準規(guī)定了兩種基本的浮點格式：單精度和雙精度，其中單精度具有23位有效數(shù)字毙死，總共占用32位燎潮，那么自然是有8位指數(shù)，1位是保留符號位扼倘，雙精度共64位确封，其中有52位尾數(shù)，11位指數(shù)再菊，1位是保留符號位爪喘。

單精度

雙精度

現(xiàn)代的X86架構早已摒棄定點數(shù)，在可以進行浮點數(shù)進行計算的計算環(huán)境中纠拔，可以自由得利用float表示小數(shù)進行計算秉剑。
細節(jié)不說了，一般不是很涉及底層的話也用不到這些知識稠诲。

四.圖像的FFT運算

我在另一篇文章里總結了侦鹏，二維的FFT可以轉(zhuǎn)換為一維的FFT進行計算诡曙，即先對每一列進行FFT運算，然后再第一次FFT運算的結果上進行列FFT變換略水，具體的證明去看那篇文章价卤，這樣二維的就可以用一維的來做，利用的是dsplib的fft32x32函數(shù)渊涝，這個函數(shù)可以快速計算FFT慎璧，先來看這個函數(shù)：

fft32x32

exe

用cmd來打開這個：

直接看截圖吧炸卑，這個是一位老師指導的，使用說明也給的很清楚煤傍，這樣之后就會生成一個文件盖文，這個文件可以是頭文件，也可以是c文件蚯姆，我是直接生成c文件五续，把系數(shù)復制到我的程序里了。

另外x[2nx]和y[2nx]都是32位的數(shù)據(jù)龄恋，虛部緊跟實部連續(xù)存儲疙驾，32位數(shù)據(jù)采用Q15定點，所以這里涉及到一個數(shù)據(jù)轉(zhuǎn)換郭毕，整個fft二維算法的大框架我是借鑒我們組另外一個老師的它碎，所以不便放出來，需要把數(shù)據(jù)轉(zhuǎn)換成Q15格式的显押，這里的轉(zhuǎn)換建議先歸一化到[-1.1),然后左移15位賦值給int型的扳肛，這樣就能得到Q15的數(shù)據(jù)。
nx是FFT變換的點數(shù)乘碑，無需多說挖息。
看下結果：對于一個64*64的圖像，我用matlab和ccs分別進行了計算兽肤，先歸一化套腹，然后去做FFT，這里只看實部：

matlab

ccs

我先把中間的結果輸出底洗，即每一列先做FFT的結果保存起來看一下腋么，我發(fā)現(xiàn)這個結果是在正確的，在matlab里和上面同樣做了差亥揖，兩位小數(shù)精度下是完全準確的珊擂。

所以說這個fft32x32的函數(shù)是不要求輸入亂序的，這個我翻手冊的時候也沒發(fā)現(xiàn)要求輸入亂序费变。我試了試暫時還沒有找到結果錯誤的原因摧扇。我也試了下看是不是亂序，以matlab_fft變換的第二列為標準挚歧，在ccs的結果中找最相似的一整列扛稽，實驗證明并不存在這樣的一整列。說明不是亂序的結果滑负。

搞了半個下午在张，還是沒有發(fā)現(xiàn)問題所在，我去干干別的再來搞這個矮慕，是在不行要去請教下老師帮匾！

18年1月19號上午更新：

剛才正在看別的東西，突然靈機一動痴鳄，想了想是不是內(nèi)存地址互相沖突呢瘟斜？于是趕緊打開CCS試了一下，第一次沒成功痪寻，趕緊檢查發(fā)現(xiàn)還是有一個變量的地址沒有完全指定螺句，導致兩個變量的地址是完全一樣的，我的程序里是xtemp和ytemp,指定了以后趕緊去看第二行第一個數(shù)橡类，果然就和以前的不同了壹蔓，導入matlab簡單畫了個圖，perfect猫态！

matlab和CCS傅里葉變換

誤差基本在兩個三個小數(shù)點之后佣蓉，可以接受！

五.圖像的IFFT運算

FFT成功之后亲雪，還是花了半個下午的時間把測試的代碼整理成函數(shù)放到頭文件里勇凭，因為代碼不多，就沒有新建源文件义辕，之所以花了半個下午是因為當時一個內(nèi)存死活寫不進去數(shù)虾标，對于一個int型的數(shù)據(jù)取float直接就是零，搞得我都想把電腦砸了灌砖，后來重啟了一下CCS竟然好了璧函，神奇的錯誤傀蚌，這也是現(xiàn)在不想做太多底層的原因，特別是和硬件相關的蘸吓。
這兩天爸媽來西安了善炫，中午剛送走，從車站趕回來剛趕上下午上班時間库继，實際上都困得不行了箩艺，就想去睡覺，到了辦公室把早上改的IFFT看了一下宪萄，其實IFFT就只需要把FFT的代碼中的dsp_fft32x32改成dsp_ifft32x32就可以了艺谆，其他的東西都不用改，因為其他得輔助代碼只是為了先行后列這樣操作罷了拜英。
但是算的結果確差了很多静汤，早上沒想到底是怎么回事，下午回來把數(shù)據(jù)導入到matlab里reshape顯示一看居凶，竟然是原圖撒妈，那么肯定是數(shù)據(jù)整個擴大了一個整倍數(shù)，找了一組對應的數(shù)字一除排监，4095.6狰右，簡直完美，這才想起逆傅里葉變換是有個系數(shù)的舆床，DSPlib中的ifft函數(shù)并沒有包含這個系數(shù)棋蚌，因為是行列做了兩次，所以應該是1/64*64=1/4096挨队，于是在最后除上這個數(shù)就ok了谷暮。
這樣這個CSK算法最難的地方就過去了。

2018/1/24更新

六：C674Xdsplib導入及使用

昨天碰見老師了盛垦，問我做的怎么樣了湿弦，交流了一下決定先用浮點的DSP來做，如果速度可以的話再移植到定點上去腾夯，所以前面用定點做的要改一改颊埃，于是又去下載了C674x的daplib的庫，這里面有一個dsp的函數(shù)是浮點的蝶俱，還有一些矩陣乘法數(shù)組乘法之類的東西班利，昨天一天都在鼓搗這個東西，因為網(wǎng)上的資料實在是太少榨呆，搞得有一陣子都想把電腦砸了罗标，不過趕在晚飯之前還是弄好了，不過就是這樣也很氣，感覺做了許多無用功闯割，沒人帶做東西真的累彻消。一晚上做夢也都是浮點定點這點破事，一天好像都在看內(nèi)存里存的什么東西宙拉。
不過抱怨歸抱怨宾尚，還是把這里總結下，如果有人能看到且要用到鼓黔，希望有用央勒。

全部復制過去缎谷，然后項目里吧include文件夾import進去，用的時候要包含相應的頭文件灶似，這個現(xiàn)在不用了不說了列林，很簡單。

重點說下dspC674x的dsplib的調(diào)用酪惭，這個方法可以用在許多l(xiāng)ib的導入上希痴，比如mathlib，DSPc674的dsplib需要用到mathlib春感，所以這兩個都導入了砌创，首先去下載兩個庫，TI的官網(wǎng)上就可以直接下鲫懒，exe格式的嫩实，安裝完成后是一個文件夾，長這樣：

unsigned char brev[64] = {
    0x0, 0x20, 0x10, 0x30, 0x8, 0x28, 0x18, 0x38,
    0x4, 0x24, 0x14, 0x34, 0xc, 0x2c, 0x1c, 0x3c,
    0x2, 0x22, 0x12, 0x32, 0xa, 0x2a, 0x1a, 0x3a,
    0x6, 0x26, 0x16, 0x36, 0xe, 0x2e, 0x1e, 0x3e,
    0x1, 0x21, 0x11, 0x31, 0x9, 0x29, 0x19, 0x39,
    0x5, 0x25, 0x15, 0x35, 0xd, 0x2d, 0x1d, 0x3d,
    0x3, 0x23, 0x13, 0x33, 0xb, 0x2b, 0x1b, 0x3b,
    0x7, 0x27, 0x17, 0x37, 0xf, 0x2f, 0x1f, 0x3f
};

另外有一個函數(shù)：

void tw_gen (float *w, int n)
{
    int i, j, k;
    double x_t, y_t, theta1, theta2, theta3;
    const double PI = 3.141592654;

    for (j = 1, k = 0; j <= n >> 2; j = j << 2)
    {
        for (i = 0; i < n >> 2; i += j)
        {
            theta1 = 2 * PI * i / n;
            x_t = cos (theta1);
            y_t = sin (theta1);
            w[k] = (float) x_t;
            w[k + 1] = (float) y_t;

            theta2 = 4 * PI * i / n;
            x_t = cos (theta2);
            y_t = sin (theta2);
            w[k + 2] = (float) x_t;
            w[k + 3] = (float) y_t;

            theta3 = 6 * PI * i / n;
            x_t = cos (theta3);
            y_t = sin (theta3);
            w[k + 4] = (float) x_t;
            w[k + 5] = (float) y_t;
            k += 6;
        }
    }
}

我想直接調(diào)用這個函數(shù)出不來输枯，覺得還不如自己生成數(shù)據(jù)拷貝進來方便议泵，F(xiàn)FT和IFFT用的是一套w，所以只需要生成一次復制進自己的頭文件或者源文件就行了桃熄。我是直接在VS里把這個函數(shù)復制進去先口，生成的數(shù)據(jù)拷貝到CCS的頭文件中，竟然出奇得順利就完成了瞳收。

未完待續(xù)---

2018.7.16： 這個暫時就終結了碉京，這個算法最終大部分的函數(shù)都寫完了但是沒有調(diào)通，主要的原因可能還是在我比較抗拒DSP吧螟深，所以做的成果都打包發(fā)給老師了谐宙，老師讓一個職工往下做了，希望可以有些不錯的結果界弧。