偏振度圖像獲取(課設(shè))

PAGEPAGE17課程設(shè)計(jì)題目：偏振度圖像的獲取學(xué)院：電子工程與自動(dòng)化學(xué)院專業(yè)：光信息科學(xué)與技術(shù)學(xué)生姓名：學(xué)號(hào)：指導(dǎo)教師：2011年12內(nèi)容摘要： 21、設(shè)計(jì)任務(wù)及要求 32、實(shí)驗(yàn)原理 33、方案論證和選擇 43.1方案一：三個(gè)不同角度測(cè)量線偏振光強(qiáng)法 43.2方案二：四個(gè)不同角度測(cè)量線偏振光強(qiáng)法 43.3方案比較與選擇 54、光路設(shè)計(jì)及器件選擇 65、實(shí)驗(yàn)原理及過(guò)程 75.1光路圖 75.2實(shí)驗(yàn)工作原理 75.3實(shí)驗(yàn)步驟 75.4實(shí)驗(yàn)記錄數(shù)據(jù) 86、圖像數(shù)據(jù)處理及分析 96.1編程計(jì)算過(guò)程 96.1.1編程工具的選擇 96.1.2圖像數(shù)據(jù)處理算法 96.2合成的強(qiáng)度圖像 106.3合成的偏振度圖像 116.4數(shù)據(jù)分析 127、課題評(píng)價(jià) 127.1總結(jié)設(shè)計(jì)的特點(diǎn) 127.2方案的優(yōu)缺點(diǎn) 137.3課題核心價(jià)值 138、總結(jié)體會(huì) 139、參與文獻(xiàn) 13附錄1：實(shí)驗(yàn)器材 14附錄2：程序代碼段 14附錄3：實(shí)驗(yàn)光路圖 18附錄4：VC++編寫(xiě)軟件界面圖 19內(nèi)容摘要：由Stokes矢量計(jì)算偏振度原理，通過(guò)使白光光源照射在電路板表面，反射光通過(guò)偏振片后由CCD探測(cè)同一目標(biāo)在不同偏振方向的反射光采集到幾組圖像數(shù)據(jù)，最終可合成偏振度圖像，結(jié)果表明偏振度圖像包含有許多普通圖像所不能體現(xiàn)出的信息，利用這些信息可以識(shí)別區(qū)分在普通圖像中難以辨別的物體，可應(yīng)用于軍事目標(biāo)識(shí)別與遙感目標(biāo)識(shí)別領(lǐng)域。關(guān)鍵字：偏振度Stokes參量 目標(biāo)識(shí)別偏振度成像VC++編程1、設(shè)計(jì)任務(wù)及要求： 任務(wù)：設(shè)計(jì)實(shí)驗(yàn)光路，通過(guò)CCD成像器件采集目標(biāo)物體的偏振信息圖像數(shù)據(jù)并使用計(jì)算機(jī)編程做數(shù)據(jù)處理最終合成偏振度圖像。要求：獲取自然光（白光）照射條件下目標(biāo)偏振圖像，并計(jì)算合成偏振度圖像，并對(duì)圖像數(shù)據(jù)結(jié)果進(jìn)行必要的分析。對(duì)所選用設(shè)備、確定的方案給出技術(shù)分析。2、實(shí)驗(yàn)原理： 由于電磁波經(jīng)物體表面反射輻射后，根據(jù)表面的結(jié)構(gòu)、紋理以及入射的角度，其偏振狀態(tài)將發(fā)生改變，會(huì)使物體表面某些信息得到增強(qiáng)，從而可以更有效地鑒別物體。所以，通過(guò)實(shí)驗(yàn)手段獲得的偏振圖像可以表征一些強(qiáng)度測(cè)量很難表征的信息，如表面的結(jié)構(gòu)、粗糙程度等。 偏振度：是度量電磁波中偏振程度的參數(shù)，為偏振光在總光強(qiáng)中所占的比例。一般而言，偏振度在0(自然光)與1(全偏振)間變化。 斯托克斯參量（Stokesparameters）：是光學(xué)偏振態(tài)的數(shù)學(xué)表示。斯托克斯引入四個(gè)量作為表示偏振態(tài)的參量，即：

S0=〈|Ex|2〉+〈|Ey|2〉 ……（1-1）

S1=〈|Ex|2〉-〈|Ey|2〉 ……（1-2）

S2=〈2ExEycosδ〉 ……（1-3）

S3=〈2ExEysinδ〉 ……（1-4）

其中Ex，Ey是光的電矢量E在x、y方向上的振幅，δ是二者的相位差。這四個(gè)參量可以表示自然光及各種偏振光。一個(gè)給定的斯托克斯參量可看做一個(gè)矢量。Stokes參量中的S0、S1、S2、S3通常表示成I、Q、U、V。其中I表示光波的總強(qiáng)度。Q表示x方向與y方向上的線偏振光的強(qiáng)度差。U表示方向與方向上的線偏振光的強(qiáng)度差。V表示右旋不是左旋圓偏振分量占優(yōu)勢(shì)。在自然界大氣背景及目標(biāo)物對(duì)太陽(yáng)入射的偏振效應(yīng)中，圓偏振的分量極少，相對(duì)于儀器的誤差來(lái)說(shuō)可以忽略，通常假定V=0。因而，要完全確定一束光線的偏振狀態(tài)，還需要三個(gè)獨(dú)立數(shù)據(jù)來(lái)確定I、Q、U這三個(gè)參量。 得到Stokes矢量S=[IQUV]后即可通過(guò)以下公式計(jì)算出該光束的偏振度為： ，(V=0) ……(1-5)3、方案論證和選擇：3.1方案一：三個(gè)不同角度測(cè)量線偏振光強(qiáng)法 由上述原理有：I=S0=〈|Ex|2〉+〈|Ey|2〉 ……（2-1）

Q=S1=〈|Ex|2〉-〈|Ey|2〉 ……（2-2）

U=S2=〈2ExEycosδ〉 ……（2-3） 在任一xoy平面，在與X軸的夾角為的方向上進(jìn)行觀測(cè)所得到的光強(qiáng)可以用下式表述： ……（2-4）或者： ……（2-5） 這樣只要測(cè)出三個(gè)不同角度的線偏振光的光強(qiáng)I，就可以聯(lián)立求出Stokes的I、Q。U的值?，F(xiàn)取0o、60o、130o三個(gè)角度獲得I、Q、U值如下： ……（2-6） ……（2-7） ……（2-8） 這樣通過(guò)實(shí)驗(yàn)使用CCD成像器件測(cè)出0o、60o、130o三個(gè)角度下的灰度圖像即可沒(méi)得三個(gè)方向的光強(qiáng)I1、I2、I3的圖像。代入以上（2-6）（2-7）（2-8）三式求出I、Q、U，然后再代入式(1-5)即可求得偏振度值P。3.2方案二：四個(gè)不同角度測(cè)量線偏振光強(qiáng)法 按式（1-1）（1-2）（1-3）（1-4）可得： ……（2-9） 其中I為光的總強(qiáng)度，I-1、I0、I1、I2分別表示放置在光傳播方向上一理想偏振片透光軸與參考方向呈-45o、0o、45o、90o方向上的線偏振光強(qiáng)，Ir和Il表示右旋（r）和左旋（l）圓偏振光強(qiáng).線偏振光的Stokes矢量為[1100]，自然光的Stokes矢量為[1000]。 這樣通過(guò)實(shí)驗(yàn)使用CCD成像器件測(cè)出-45o、0o、45o、90o四個(gè)角度下的灰度圖像即可沒(méi)得四個(gè)方向的光強(qiáng)I-1、I0、I1、I2的圖像。代入上式（2-9）求出I、Q、U，然后再代入式(1-5)即可求得偏振度值P。3.3方案比較與選擇 兩個(gè)方案相比看出，在實(shí)驗(yàn)步驟上兩者大同小異，主要區(qū)別在于數(shù)據(jù)采集和數(shù)據(jù)處理過(guò)程。使用方案一則只需使檢偏器分別旋轉(zhuǎn)3個(gè)角度拍取圖像即可，由于在實(shí)驗(yàn)過(guò)程中，去旋轉(zhuǎn)檢偏器多多少少會(huì)使鏡片產(chǎn)生位移，這樣可能會(huì)引起拍攝圖像產(chǎn)生位移。3個(gè)角度方法，即只需旋轉(zhuǎn)兩次，相比4個(gè)角度方法少轉(zhuǎn)一次，從而減少由于鏡片位移引起的誤差。 使用方案二則需要測(cè)量4組數(shù)據(jù)，相比方案一多測(cè)1組。這樣360o劃分為45o，而方案二劃分為60o，所以相對(duì)于方案二，方案一計(jì)算測(cè)量的分辨率要高，從而計(jì)算誤差比方案一的要小。 綜上所述，最終選用方案二。4、光路設(shè)計(jì)及器件選擇 由方案二原理分析得，所需器件有：白光光源、被測(cè)目標(biāo)物（實(shí)驗(yàn)所用電路板）、偏振片（檢偏器）、成像透鏡、CCD成像探測(cè)器、計(jì)算機(jī)、固定支架等。 由于白光光源的光譜分布最接近于自然光的光譜分布，所以其光束各方向偏振強(qiáng)度均勻，比較符合實(shí)驗(yàn)原理要求。當(dāng)白光照射到被測(cè)物時(shí)，其反射的光偏振方向主要就取決于目標(biāo)物表面性質(zhì)的分布。 被測(cè)目標(biāo)物的選擇：由于要獲取的圖像為偏振度圖像，則圖像應(yīng)該具有明顯的輪廓和層次感，這樣才能體現(xiàn)出由于被測(cè)物表面性質(zhì)分布特征而反映出的物體表面某些信息。所以要選擇的被測(cè)物體表面性質(zhì)、結(jié)構(gòu)、粗糙程度等應(yīng)有明顯的變化。實(shí)驗(yàn)初步使用樹(shù)葉，用偏振度圖像反映樹(shù)葉的脈絡(luò)。另外也可以選擇表面物質(zhì)有性質(zhì)有明顯變化的物體，如電路板。實(shí)驗(yàn)最終選用性質(zhì)變化明顯的電路板。 偏振片：其對(duì)入射光具有遮蔽和透過(guò)的功能，可使縱向光或橫向光一種透過(guò)，一種遮蔽。在實(shí)驗(yàn)中，其起著關(guān)鍵作用。實(shí)驗(yàn)通過(guò)旋轉(zhuǎn)偏振片于不同方向可以讓不同方面的線偏振光束透過(guò)成為線偏振圖像。 成像透鏡：主要是作為CCD成像采集的輔助作用，其將被測(cè)物體成像于CCD探測(cè)面上。為了保證物體表面層次成像有好效果，應(yīng)使用光學(xué)傳遞函數(shù)MTF值在中低頻的強(qiáng)度曲線呈平緩的成像透鏡器件。 CCD成像器件：其主要用于光線偏振圖像的光強(qiáng)度直接采集，直接關(guān)系的源數(shù)據(jù)信息及誤差，所應(yīng)盡可以選擇分辨率高和暗電流小的CCD器件。 計(jì)算機(jī)：主要用于采集的圖像數(shù)據(jù)處理。實(shí)際編程計(jì)算時(shí)可使用MATLAB或VC++6.0等編程開(kāi)發(fā)軟件。為了保證數(shù)據(jù)處理的速度，應(yīng)先用內(nèi)存大、運(yùn)算速度高的計(jì)算機(jī)。5、實(shí)驗(yàn)原理及過(guò)程5.1光路圖：白光光源白光光源被測(cè)目標(biāo)偏振片成像透鏡CCD成像器件（約100o）圖1實(shí)驗(yàn)光路圖5.2實(shí)驗(yàn)工作原理： 如圖1實(shí)驗(yàn)光路圖所示，白光光束照射到物體表面后CCD探測(cè)方向的反射光線與入射光線之間的夾角應(yīng)大于或等于90度（實(shí)際選用約100度），這樣反射光的偏振信息才明顯。反射光束經(jīng)過(guò)偏振片后透過(guò)的只有線偏振光束，這樣得到的則為線偏振圖像。通過(guò)成像透鏡后，CCD探測(cè)采集到線偏振圖像數(shù)據(jù)。 實(shí)驗(yàn)中，通過(guò)旋轉(zhuǎn)偏振片分別為-45o、0o、45o、90o四個(gè)角度方向則可以采集到相應(yīng)方向上的線偏振圖像。通過(guò)數(shù)據(jù)處理，將采集到的灰度強(qiáng)度數(shù)據(jù)代入式（2-9）求出I、Q、U，然后再代入式(1-5)即可求得偏振度值P。最終將所有計(jì)算結(jié)果的像素合成的圖像即為偏振度圖像。5.3實(shí)驗(yàn)步驟：擺光路器件：按圖1實(shí)驗(yàn)光路圖，除偏振片外，按光路順序擺放好各器件，其中將電路板（目標(biāo)物）以白紙為背景固定在目標(biāo)靶上。調(diào)焦：在光軸方向上移動(dòng)成像透鏡，觀察采集圖像的，使圖像達(dá)最清晰程度，并固定好所有器件。采集普通圖像：如圖1，在成像透鏡前加入偏振片并旋轉(zhuǎn)偏振片角度使采集在的圖像呈最明亮的程度，記錄該圖像為普通圖像。采集線偏振圖像：按偏振片上的刻度值，分別將偏振片與參考刻度成：-45o、0o、45o、90o夾角，采集記錄相應(yīng)的圖像：I-1、I0、I1、I2。數(shù)據(jù)處理：將步驟4）采集到的I-1、I0、I1、I2圖像灰度值代入式（2-9）求出I、Q、U，然后再代入式(1-5)即可求得偏振度值P。最終將所有計(jì)算結(jié)果的像素合成的圖像即為偏振度圖像。5.4實(shí)驗(yàn)記錄數(shù)據(jù)：a.a.-45o方向線偏振圖像b.0o方向線偏振圖像c.45o方向線偏振圖像d.90o方向線偏振圖像圖2原始采集圖像6、圖像數(shù)據(jù)處理及分析6.1編程計(jì)算過(guò)程：6.1.1編程工具的選擇： 在圖像處理編程中，較好常用的是MATLAB軟件編程，其特點(diǎn)是方便程序設(shè)計(jì)，許多圖像處理的函數(shù)如：圖像文件輸入輸出（imread()、imwrite()）、圖像轉(zhuǎn)換（rgb2gray()）、圖像灰度增強(qiáng)運(yùn)算、數(shù)字圖像空域?yàn)V波等已經(jīng)包含在數(shù)據(jù)庫(kù)中，用戶只需調(diào)用這些現(xiàn)成的函數(shù)來(lái)處理相應(yīng)的數(shù)據(jù)及可快捷完成實(shí)現(xiàn)圖像的處理。其主要缺點(diǎn)是，一般設(shè)計(jì)的MATLAB程序只能運(yùn)行在MATLAB軟件的基礎(chǔ)上，即程序只能在打開(kāi)了MATLAB軟件通過(guò)它來(lái)調(diào)出程序才能正常運(yùn)行，程序獨(dú)立性較差，工作平臺(tái)受到了限制。 若在沒(méi)有涉及到太多的圖像數(shù)據(jù)處理，如：圖像增強(qiáng)，圖像濾波，圖像復(fù)原一的些復(fù)雜的運(yùn)算，使用VC++6.0通用的軟件編程不但可以完成簡(jiǎn)單的圖像處理功能，還解決了MATLAB中程序獨(dú)立性差的缺點(diǎn)。因?yàn)閂C++編程及為面向windows系統(tǒng)的編程，其編寫(xiě)完成的可執(zhí)行文件可以直接在操作系統(tǒng)上運(yùn)行。這樣，通過(guò)完成軟件成品的編程，在實(shí)驗(yàn)數(shù)據(jù)處理中很方便快捷的完成處理得到運(yùn)算結(jié)果。 相比上述兩種編程工具，雖然MATLAB方便設(shè)計(jì)，但不方便實(shí)驗(yàn)中使用來(lái)做數(shù)據(jù)處理，而使用VC++6.0雖然工作量較大，但更有利于實(shí)驗(yàn)數(shù)據(jù)處理使用方便，所以最終選用VC++6.0編程工具。6.1.2圖像數(shù)據(jù)處理算法： 圖像數(shù)據(jù)的處理主要是將所采集到的4個(gè)角度方向的線偏振圖像數(shù)據(jù)I-1、I0、I1、I2圖像灰度值代入式（2-9）求出I、Q、U，然后再代入式(1-5)即可求得偏振度值P。最終將所有計(jì)算結(jié)果的像素合成的圖像即為偏振度圖像。開(kāi)始i=0i%3=0?計(jì)算新偏振度值及強(qiáng)度所有參與運(yùn)算的數(shù)據(jù)塊指針加1i++i<圖像數(shù)據(jù)字節(jié)數(shù)結(jié)束偏振度及強(qiáng)度數(shù)據(jù)塊賦值圖3偏振度算法流程圖 通過(guò)讀取采集的BMP位圖文件，將四個(gè)角度的線偏振圖像數(shù)據(jù)保存在頭指針?lè)謩e為：lpImgDataN45，lpImgData0，開(kāi)始i=0i%3=0?計(jì)算新偏振度值及強(qiáng)度所有參與運(yùn)算的數(shù)據(jù)塊指針加1i++i<圖像數(shù)據(jù)字節(jié)數(shù)結(jié)束偏振度及強(qiáng)度數(shù)據(jù)塊賦值圖3偏振度算法流程圖如附錄2中代碼段所示，數(shù)據(jù)處理通過(guò)使用for循環(huán)語(yǔ)句逐個(gè)查找源數(shù)據(jù)的4個(gè)數(shù)據(jù)塊的元素來(lái)計(jì)算出偏振度對(duì)應(yīng)的灰度并存到計(jì)算結(jié)果數(shù)據(jù)塊中。最終發(fā)送PAINT消息，由OnPaint消息響應(yīng)函數(shù)來(lái)顯示計(jì)算結(jié)果的偏振度圖像。6.2合成的強(qiáng)度圖像：圖4強(qiáng)度圖像6.3合成的偏振度圖像：圖5偏振度圖像6.4數(shù)據(jù)分析： 如圖4的合成偏振度圖像所示，從中可看出圖像主要以被腐蝕了表面銅的塑料板電路為背景，可明顯的看出各不同器件的輪廓。電路板中包含有如下器件：直插電阻、二極管、三極管、電解電容。其均由不同的材料組成，如下：a.強(qiáng)度圖像部分區(qū)域b.偏振度圖像部分區(qū)域圖6強(qiáng)度圖像與偏振度圖像對(duì)比 電路板材料：合成塑料板； 電阻表面：碳膜金屬膜氧化膜；二極管：硅及一些摻雜物；三極管：硅、塑料、環(huán)氧樹(shù)脂；電解電容：通常是由金屬箔和絕緣氧化層（氧化鋁/鉭五氧化物）。a.強(qiáng)度圖像部分區(qū)域b.偏振度圖像部分區(qū)域圖6強(qiáng)度圖像與偏振度圖像對(duì)比 從圖6的亮暗對(duì)比程序看出，從以上同一區(qū)域的不同圖像看出，在強(qiáng)度圖像中分辨不明顯的三極管和電解電容在偏振度圖像中得到明顯區(qū)分。偏振圖像中三極管偏振度最高，其次是電阻和電解電容的偏振度相當(dāng)，而二極管的偏振度較低；其中偏振度表現(xiàn)最低的則為電路塑料板。 從偏振圖像中各元器件偏振度的對(duì)比可明顯發(fā)現(xiàn)，偏振度圖像表征了物體材料的物理化學(xué)特性。通過(guò)這點(diǎn)可說(shuō)明，偏振成像探測(cè)技術(shù)在云和大氣溶膠的探測(cè)、地質(zhì)勘探、土壤分析、環(huán)境監(jiān)測(cè)、材料分析等領(lǐng)域具有非常長(zhǎng)遠(yuǎn)的應(yīng)用價(jià)值。7、課題評(píng)價(jià)7.1總結(jié)設(shè)計(jì)的特點(diǎn)：該圖像偏振度的獲取方法是較為常用的的方法，通過(guò)使用偏振處過(guò)濾光束實(shí)現(xiàn)線偏振圖像的獲取，最后將4個(gè)角度的線偏振圖像通過(guò)公式計(jì)算得偏振度圖像。與通常方案相比，本方案的突出特點(diǎn)在于圖像數(shù)據(jù)處理環(huán)節(jié)。本方案采用VC++編程工具，實(shí)現(xiàn)獨(dú)立運(yùn)行于Windows操作系統(tǒng)上的圖像處理。7.2方案的優(yōu)缺點(diǎn): 方案優(yōu)點(diǎn)：使用VC++界面控制，在實(shí)驗(yàn)數(shù)據(jù)處理中方便快捷，具有實(shí)際應(yīng)用價(jià)值，在VC++上的設(shè)計(jì)比較靈活，可根據(jù)實(shí)驗(yàn)需要添加設(shè)計(jì)新功能。 方案缺點(diǎn)：未能達(dá)到實(shí)時(shí)性采集獲取偏振度圖像；實(shí)驗(yàn)裝置簡(jiǎn)單，測(cè)量數(shù)據(jù)容易外面干擾產(chǎn)生較大誤差。7.3課題核心價(jià)值： 本次課題主要在光偏振度的研究基礎(chǔ)上設(shè)計(jì)方案采集合成偏振度圖像，通過(guò)實(shí)驗(yàn)數(shù)據(jù)的處理比較，論證偏振度圖像在研究物體成像中所表征一些偏振信息在物體識(shí)別中所起到的突出作用。通過(guò)課題認(rèn)識(shí)到偏振成像探測(cè)技術(shù)在大氣溶膠的探測(cè)、地質(zhì)勘探、土壤分析、環(huán)境監(jiān)測(cè)、材料分析等領(lǐng)域具有非常長(zhǎng)遠(yuǎn)的應(yīng)用價(jià)值。8、總結(jié)體會(huì)： 通過(guò)本次課程設(shè)計(jì)，本人對(duì)光學(xué)技術(shù)中的偏振度圖像有了進(jìn)一步了解，認(rèn)識(shí)到了偏振成像探測(cè)技術(shù)在實(shí)際自然界研究起著重要的作用。在實(shí)驗(yàn)數(shù)據(jù)采集中，了解到了各器件的使用細(xì)節(jié)，提高了動(dòng)手能力。通過(guò)采用VC++編程對(duì)Windows程序設(shè)計(jì)和調(diào)試有了初步認(rèn)識(shí)，積累了初步經(jīng)驗(yàn)，培養(yǎng)了能通過(guò)使用計(jì)算機(jī)處理能力了對(duì)實(shí)驗(yàn)數(shù)據(jù)處理的能力。能使用所學(xué)知識(shí)解決現(xiàn)實(shí)中所遇到的問(wèn)題是一件非常有意義的事，通過(guò)本次課設(shè)也進(jìn)一步鞏固了大學(xué)期間所學(xué)的知識(shí)。 9、參與文獻(xiàn)[1]唐若愚,于國(guó)萍,王曉峰.自然光照下偏振度圖像的獲取方法[J].武漢大學(xué)學(xué)報(bào),2006,52(1),59~60.[2] 于國(guó)萍,王曉峰.偏振圖像的研究[J].物理實(shí)驗(yàn),2009,29(7),38~42.[3]曾延安,馬娟,常大定.基于VC++的偏振度圖像獲取[J].微計(jì)算機(jī)信息,2007,23(8),304~306.[4]任哲.Windows程序設(shè)計(jì)技術(shù)基礎(chǔ)：MFC與.NET.窗口軟件,Windows.北京：機(jī)械工業(yè)出版社，2009附錄1：實(shí)驗(yàn)器材實(shí)驗(yàn)器材數(shù)量備注白光光源1實(shí)驗(yàn)室中的白光點(diǎn)光源偏振片1成像透鏡1CCD1計(jì)算機(jī)1支架若干被測(cè)目標(biāo)若干樹(shù)葉、電路板附錄2：程序代碼段：BMP文件的讀取與保存：以下為BMP文件的讀取程序段：BYTE*lptImgData; //內(nèi)存數(shù)據(jù)塊指針LPCTSTRlpszFilter="BMPFiles(*.bmp)|*.bmp|任何文件|*.*||"; //文件過(guò)濾設(shè)置CFileDialogfiledlg(TRUE,NULL,NULL,OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,lpszFilter,NULL); //定義文件對(duì)話框?qū)ο驝Stringfilename;CFilefile;BYTE*lptImgData;if(filedlg.DoModal()==IDOK) //調(diào)出文件對(duì)話框{ filename=filedlg.GetPathName(); //獲取文件名 //打開(kāi)文件 file.Open(filename,CFile::modeRead|CFile::shareDenyNone,NULL); //開(kāi)辟內(nèi)存 lptImgData=(BYTE*)HeapAlloc(GetProcessHeap(),HEAP_ZERO_MEMORY,m_bmih.biSizeImage); file.Read(lptImgData,m_bmih.biSizeImage); //將文件數(shù)據(jù)讀入內(nèi)存 file.Close(); //關(guān)閉文件}以下為BMP文件的保存程序段：LPCTSTRlpszFilter="BMPFiles(*.bmp)|*.bmp|任何文件|*.*||"; //文件過(guò)濾設(shè)置CFileDialogfiledlg(FALSE,lpszFilter,"PolarImg.bmp",OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,lpszFilter,NULL); //定義文件對(duì)話框?qū)ο驝Stringfilename;if(IDOK==filedlg.DoModal()) //調(diào)出文件對(duì)話框{CFilefile(filedlg.GetPathName(),CFile::modeCreate|CFile::modeReadWrite); filename=file.GetFilePath(); //獲取文件 //將內(nèi)存數(shù)據(jù)寫(xiě)入磁盤(pán)文件 file.Write(m_lpImgDataResult,m_bmih.biSizeImage); file.Close(); //關(guān)閉文件}以下為源圖像數(shù)據(jù)計(jì)算偏振度圖像的程序段：intI,Q,U,P; //定義公式參數(shù)I,Q,U,PDWORDi; BYTE*ptImg0,*ptImg45,*ptImg90,*ptImgN45,*ptImgResult;//定義各數(shù)據(jù)塊指針ptImg0=m_lpImgData0; //局部指針與全局指針關(guān)聯(lián)ptImg45=m_lpImgData45;ptImg90=m_lpImgData90;ptImgN45=m_lpImgDataN45;m_lpImgDataResult=(BYTE*)HeapAlloc(GetProcessHeap(),HEAP_ZERO_MEMORY,m_bmih.biSizeImage); //開(kāi)辟存放計(jì)算結(jié)果數(shù)據(jù)塊ptImgResult=m_lpImgDataResult;for(i=0;i<m_bmih.biSizeImage;i++){//取源數(shù)據(jù)塊中被3整除的元素計(jì)算偏振度 if(i%3==0) { I=*ptImg0+*ptImg90;