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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 波变换 > 查看源码
小波变换Doc.cpp
资源名称：imagefusionwavelettransform.rar　[点击查看]
上传用户：aqingfeng
上传日期：2014-03-25
资源大小：1839k
文件大小：38k
源码类别：

波变换

开发平台：

Visual C++

小波变换Doc.cpp：源码内容
  1. // 小波变换Doc.cpp : implementation of the CMyDoc class
  2. //
  4. #include "stdafx.h"
  5. #include "小波变换.h"
  6. #include "WvltTrans.h"
  7. #include "小波变换Doc.h"
  8. #include "Diproc.h"
  9. #include "WFilter.h"
  10. #include "math.h"
  13. #ifdef _DEBUG
  14. #define new DEBUG_NEW
  15. #undef THIS_FILE
  16. static char THIS_FILE[] = __FILE__;
  17. #endif
  19. /////////////////////////////////////////////////////////////////////////////
  20. // CMyDoc
  22. IMPLEMENT_DYNCREATE(CMyDoc, CDocument)
  24. BEGIN_MESSAGE_MAP(CMyDoc, CDocument)
  25. //{{AFX_MSG_MAP(CMyDoc)
  26. ON_COMMAND(ID_FILE_REOPEN, OnFileReopen)
  28. ON_COMMAND(IDM_WVLT_HORTRANS, OnWvltHortrans)
  29. ON_COMMAND(IDM_WVLT_VERTRANS, OnWvltVertrans)
  30. ON_COMMAND(IDM_WVLT_ONCE, OnWvltTransOnce)
  31. ON_COMMAND(IDM_WVLT_TWICE, OnWvltTransTwice)
  32. ON_COMMAND(IDM_WVLT_TRBL, OnWvltTransTrbl)
  34. ON_COMMAND(IDM_DIPROC_FUSION, OnDiprocFusion)
  35. ON_COMMAND(IDM_DIPROC_REVER, OnDiprocRever)
  36. ON_COMMAND(IDM_FILTER_BLUR, OnFilterBlur)
  37. ON_COMMAND(IDM_FILTER_SHARPNESS, OnFilterSharpness)
  38. ON_COMMAND(IDM_FILTER_BLUR2, OnFilterBlur2)
  39. ON_COMMAND(IDM_FILTER_SHARPNESS2, OnFilterSharpness2)
  41.    
  42. //}}AFX_MSG_MAP
  43. END_MESSAGE_MAP()
  45. /////////////////////////////////////////////////////////////////////////////
  46. // CMyDoc construction/destruction
  48. CMyDoc::CMyDoc()
  49. {
  50. // TODO: add one-time construction code here
  51. m_fDirty = FALSE;
  52. m_pBitmap = NULL ;
  53. m_pTransfered = NULL;
  54. m_bTribl = FALSE; m_bOnce = FALSE; m_bTwice = FALSE; m_bFilter = FALSE;
  56. m_hDIB = NULL;
  57. m_palDIB = NULL;
  58. m_sizeDoc = CSize(1,1);
  59. m_hMulDIB= NULL;
  60. }
  62. CMyDoc::~CMyDoc()
  63. {
  64. if( m_pBitmap != NULL)
  65. free( m_pBitmap);
  66. if( m_pTransfered != NULL)
  67. free( m_pTransfered );
  68. // 判断DIB对象是否存在
  69. if (m_hDIB != NULL)
  70. {
  71. // 清除DIB对象
  72. ::GlobalFree((HGLOBAL) m_hDIB);
  73. }
  75. // 判断调色板是否存在
  76. if (m_palDIB != NULL)
  77. {
  78. // 清除调色板
  79. delete m_palDIB;
  80. }
  81. }
  84. /////////////////////////////////////////////////////////////////////////////
  85. // CMyDoc serialization
  87. void CMyDoc::Serialize(CArchive& ar)
  88. {
  89. if (ar.IsStoring())
  90. {
  91. // TODO: add storing code here
  92. }
  93. else
  94. {
  95. // TODO: add loading code here
  96. }
  97. }
  99. /////////////////////////////////////////////////////////////////////////////
  100. // CMyDoc diagnostics
  102. #ifdef _DEBUG
  103. void CMyDoc::AssertValid() const
  104. {
  105. CDocument::AssertValid();
  106. }
  108. void CMyDoc::Dump(CDumpContext& dc) const
  109. {
  110. CDocument::Dump(dc);
  111. }
  112. #endif //_DEBUG
  114. /////////////////////////////////////////////////////////////////////////////
  115. // CMyDoc commands
  116. BOOL CMyDoc::OnNewDocument()//建立新文档
  117. {
  118. if (!CDocument::OnNewDocument())
  119. return FALSE;
  121. // TODO: add reinitialization code here
  122. // (SDI documents will reuse this document)
  124. return TRUE;
  125. }
  126. BOOL CMyDoc::OnOpenDocument(LPCTSTR lpszPathName)//打开已有文档，lpszPathName是待打开文件的文件名 
  127. {
  128. CFile file;
  129. CFileException fe;
  130. // 打开文件
  131. if (!file.Open(lpszPathName, CFile::modeRead | CFile::shareDenyWrite, &fe))
  132. {
  133. // 失败
  134. ReportSaveLoadException(lpszPathName, &fe,FALSE, AFX_IDP_FAILED_TO_OPEN_DOC);
  135. // 返回FALSE
  136. return FALSE;
  137. }
  138. DeleteContents();
  139. // 更改光标形状
  140. BeginWaitCursor();
  141. // 尝试调用ReadDIBFile()读取图像
  142. TRY
  143. {
  144. m_hDIB = ::ReadDIBFile(file);
  145. }
  146. CATCH (CFileException, eLoad)
  147. {
  148. // 读取失败
  149. file.Abort();
  150. // 恢复光标形状
  151. EndWaitCursor();
  152. // 报告失败
  153. ReportSaveLoadException(lpszPathName, eLoad,FALSE, AFX_IDP_FAILED_TO_OPEN_DOC);
  154. // 设置DIB为空
  155. m_hDIB = NULL;
  156. // 返回FALSE
  157. return FALSE;
  158. }
  159. END_CATCH
  160. // 初始化DIB
  161. InitDIBData();
  162. // 恢复光标形状
  163. EndWaitCursor();
  164. // 判断读取文件是否成功
  165. if (m_hDIB == NULL)
  166. {
  167. // 失败，可能非BMP格式
  168. CString strMsg;
  169. strMsg = "读取图像时出错！可能是不支持该类型的图像文件！";
  170. // 提示出错
  171. MessageBox(NULL, strMsg, "系统提示", MB_ICONINFORMATION | MB_OK);
  172. // 返回FALSE
  173. return FALSE;
  174. }
  175. // 设置文件名称
  176. SetPathName(lpszPathName);
  177. // 初始化胀标记为FALSE
  178. SetModifiedFlag(FALSE);
  179. // 返回TRUE
  180. return TRUE;
  181. }
  184.  void CMyDoc::OnFileReopen() 
  185. {
  186.   // 重新打开图像，放弃所有修改
  187.  
  188. // 判断当前图像是否已经被改动
  189. if (IsModified())
  190. {
  191. // 提示用户该操作将丢失所有当前的修改
  192. if (MessageBox(NULL, "重新打开图像将丢失所有改动！是否继续？", "系统提示", MB_ICONQUESTION | MB_YESNO) == IDNO)
  193. {
  194. // 用户取消操作，直接返回
  195. return;
  196. }
  197. }
  199. CFile file;
  200.   CFileException fe;
  201.  
  202.   CString strPathName;
  203.   
  204. // 获取当前文件路径
  205.   strPathName = GetPathName();
  206.  
  207.   // 重新打开文件
  208.   if (!file.Open(strPathName, CFile::modeRead | CFile::shareDenyWrite, &fe))
  209.   {
  210.   // 失败
  211.   ReportSaveLoadException(strPathName, &fe,
  212.   FALSE, AFX_IDP_FAILED_TO_OPEN_DOC);
  213.  
  214.   // 返回
  215.   return;
  216.   }
  217.  
  218.   // 更改光标形状
  219.   BeginWaitCursor();
  220.  
  221.   // 尝试调用ReadDIBFile()读取图像
  222.   TRY
  223.   {
  224.   m_hDIB = ::ReadDIBFile(file);
  225.   }
  226.   CATCH (CFileException, eLoad)
  227.   {
  228.   // 读取失败
  229.   file.Abort();
  230.  
  231.   // 恢复光标形状
  232.   EndWaitCursor();
  233.  
  234.   // 报告失败
  235.   ReportSaveLoadException(strPathName, eLoad,
  236.   FALSE, AFX_IDP_FAILED_TO_OPEN_DOC);
  237.  
  238.   // 设置DIB为空
  239.   m_hDIB = NULL;
  240.  
  241.   // 返回
  242.   return;
  243.   }
  244.   END_CATCH
  245.  
  246.   // 初始化DIB
  247.   InitDIBData();
  248.  
  249.   // 判断读取文件是否成功
  250.   if (m_hDIB == NULL)
  251.   {
  252.   // 失败，可能非BMP格式
  253.   CString strMsg;
  254.   strMsg = "读取图像时出错！可能是不支持该类型的图像文件！";
  255.  
  256.   // 提示出错
  257.   MessageBox(NULL, strMsg, "系统提示", MB_ICONINFORMATION | MB_OK);
  258.   
  259.   // 恢复光标形状
  260.   EndWaitCursor();
  261.  
  262.   // 返回
  263.   return;
  264.   }
  265.  
  266.   // 初始化胀标记为FALSE
  267.   SetModifiedFlag(FALSE);
  268.  
  269.   // 刷新
  270.   UpdateAllViews(NULL);
  271.   
  272.   // 恢复光标形状
  273.   EndWaitCursor();
  275.   // 返回
  276.   return;
  277.  
  278. }
  282. void CMyDoc::InitDIBData()
  283. {
  284. // 初始化DIB对象
  285. // 判断调色板是否为空
  286. if (m_palDIB != NULL)
  287. {
  288. // 删除调色板对象
  289. delete m_palDIB;
  290. // 重置调色板为空
  291. m_palDIB = NULL;
  292. }
  293. // 如果DIB对象为空，直接返回
  294. if (m_hDIB == NULL)
  295. {
  296. // 返回
  297. return;
  298. }
  299. LPSTR lpDIB = (LPSTR) ::GlobalLock((HGLOBAL) m_hDIB);
  300. // 判断图像是否过大
  301. if (::DIBWidth(lpDIB) > INT_MAX ||::DIBHeight(lpDIB) > INT_MAX)
  302. {
  303. ::GlobalUnlock((HGLOBAL) m_hDIB);
  304. // 释放DIB对象
  305. ::GlobalFree((HGLOBAL) m_hDIB);
  306. // 设置DIB为空
  307. m_hDIB = NULL;
  308. CString strMsg;
  309. strMsg = "BMP图像太大！";
  310. // 提示用户
  311. MessageBox(NULL, strMsg, "系统提示", MB_ICONINFORMATION | MB_OK);
  312. // 返回
  313. return;
  314. }
  315. // 设置文档大小
  316. m_sizeDoc = CSize((int) ::DIBWidth(lpDIB), (int) ::DIBHeight(lpDIB));
  317. ::GlobalUnlock((HGLOBAL) m_hDIB);
  318. // 创建新调色板
  319. m_palDIB = new CPalette;
  320. // 判断是否创建成功
  321. if (m_palDIB == NULL)
  322. {
  323. // 失败，可能是内存不足
  324. ::GlobalFree((HGLOBAL) m_hDIB);
  325. // 设置DIB对象为空
  326. m_hDIB = NULL;
  327. // 返回
  328. return;
  329. }
  330. // 调用CreateDIBPalette来创建调色板
  331. if (::CreateDIBPalette(m_hDIB, m_palDIB) == NULL)
  332. {
  333. // 返回空，可能该DIB对象没有调色板
  334. // 删除
  335. delete m_palDIB;
  336. // 设置为空
  337. m_palDIB = NULL;
  338. // 返回
  339. return;
  340. }
  342. }
  343. BOOL CMyDoc::OnSaveDocument(LPCTSTR lpszPathName) 
  344. {
  345. CFile file;
  346. CFileException fe;
  347. // 打开文件
  348. if (!file.Open(lpszPathName, CFile::modeCreate |
  349.   CFile::modeReadWrite | CFile::shareExclusive, &fe))
  350. {
  351. // 失败
  352. ReportSaveLoadException(lpszPathName, &fe,TRUE, AFX_IDP_INVALID_FILENAME);
  353. // 返回FALSE
  354. return FALSE;
  355. }
  356. // 尝试调用SaveDIB保存图像
  357. BOOL bSuccess = FALSE;
  358. TRY
  359. {
  360. // 更改光标形状
  361. BeginWaitCursor();
  362. // 尝试保存图像
  363. bSuccess = ::SaveDIB(m_hDIB, file);
  364. // 关闭文件
  365. file.Close();
  366. }
  367. CATCH (CException, eSave)
  368. {
  369. // 失败
  370. file.Abort();
  371. // 恢复光标形状
  372. EndWaitCursor();
  373. ReportSaveLoadException(lpszPathName, eSave,TRUE, AFX_IDP_FAILED_TO_SAVE_DOC);
  374. // 返回FALSE
  375. return FALSE;
  376. }
  377. END_CATCH
  378. // 恢复光标形状
  379. EndWaitCursor();
  380. // 重置胀标记为FALSE
  381. SetModifiedFlag(FALSE);
  382. if (!bSuccess)
  383. {
  384. // 保存失败，可能是其它格式的DIB，可以读取但是不能保存
  385. // 或者是SaveDIB函数有误
  386. CString strMsg;
  387. strMsg = "无法保存BMP图像！";
  388. // 提示出错
  389. MessageBox(NULL, strMsg, "系统提示", MB_ICONINFORMATION | MB_OK);
  390. }
  391. return bSuccess;
  392. }
  396. void CMyDoc::OnWvltHortrans() 
  397. {
  398. // TODO: Add your command handler code here
  399. //读取数字图像的文件头，获取图像的属性参数
  400. LPSTR lpDIB;
  401. // 指向DIB象素指针
  402. LPSTR    lpDIBBits;
  403. // 锁定DIB
  404. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  405. // 找到DIB图像象素起始位置
  406. lpDIBBits = ::FindDIBBits(lpDIB);
  407. unsigned long biWidth =::DIBWidth(lpDIB);
  408. unsigned long biHeight = ::DIBHeight(lpDIB);
  409. unsigned long biAlign = (biWidth*3+3)/4 *4;
  410. unsigned long bmSize = biHeight * biAlign;
  411. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  413. //图像矩阵坐标与像素数值
  414. int x,y,cur;
  415. unsigned char tempR, tempG, tempB;
  416. short **spOriginData, **spTransData0;
  417. float fTempBufforDisp;
  418. short MaxPixVal,MinPixVal,Diff;
  419. //分配数据空间
  420. spOriginData = new short* [biHeight];
  421. spTransData0 = new short* [biHeight];
  422. //
  423. for(int i = 0; i < biHeight; i ++)
  424. {
  425. spOriginData[i] = new short [biWidth];
  426. spTransData0[i] = new short [biWidth];
  427. }
  428. //创建图像小波变换类
  429. CWvltTrans *pTrans;
  430. //从设备缓存中获得原始图像数据
  431. for(y=0; y<(int)biHeight; y++)
  432. {
  433. for( x=0; x<(int)biWidth; x++)
  434. {
  435. cur = y*biAlign+3*x;
  436. tempB=lpData[cur];
  437. tempG=lpData[cur+1];
  438. tempR=lpData[cur+2];
  439. spOriginData[biHeight-1-y][x]=(short)(0.3*tempR+0.59*tempG+0.11*tempB);  //no problem
  440. }
  441. }
  442. //完成图像水平方向的一次小波变换
  443. pTrans->Hori_Transform(spOriginData,spTransData0,biHeight,biWidth/2,1);
  444. //屏蔽图像复原操作标志
  445. m_bFilter = TRUE;
  446. m_bFilter = m_bFilter & ~m_bOnce & ~m_bTwice & ~m_bTribl;
  447. MaxPixVal=spTransData0[0][0];
  448. MinPixVal=spTransData0[0][0];
  449. //得到图像小波系数的极大值与极小值
  450. for( y=0; y<(int)biHeight; y++)
  451. {
  452. for( x=0; x<(int)biWidth; x++)
  453. {
  454. if(MaxPixVal<spTransData0[y][x])
  455. MaxPixVal=spTransData0[y][x];
  456. if(MinPixVal>spTransData0[y][x])
  457. MinPixVal=spTransData0[y][x];
  458. }
  459. }
  460. //计算获取图像小波系数的极值差
  461. Diff=MaxPixVal-MinPixVal;
  462. //小波经过处理后，放入显示缓存中
  463. for(y=0; y<(int)biHeight; y++)
  464. {
  465. for(x=0; x<(int)biWidth; x++)
  466. {
  467. //因为小波变换后的小波系数有可能超过255甚至更多，那么就将
  468. //小波系数的范围映射到0~255区间内，以后出现类似的处理，目的都是一样的
  469. fTempBufforDisp=spTransData0[biHeight-1-y][x];
  470. fTempBufforDisp-=MinPixVal;
  471. fTempBufforDisp*=255;
  472. fTempBufforDisp/=Diff;
  473. cur= y*biAlign+3*x; //current pixel
  474. lpData[cur] = (unsigned char)fTempBufforDisp;
  475. lpData[cur+1]= (unsigned char)fTempBufforDisp;
  476. lpData[cur+2]= (unsigned char)fTempBufforDisp;
  477. }
  478. }
  479. //小波图像水平方向小波变换
  480. SetModifiedFlag(TRUE);
  481. UpdateAllViews(NULL);
  482. ::GlobalUnlock((HGLOBAL)GetHDIB());
  483. //删除临时的数据空间
  484. delete spOriginData;
  485. delete spTransData0;
  486. }
  488. void CMyDoc::OnWvltVertrans() 
  489. {
  490. // TODO: Add your command handler code here
  491. //读取数字图像的文件头，获取图像的属性参数
  492. LPSTR lpDIB;
  493. // 指向DIB象素指针
  494. LPSTR    lpDIBBits;
  495. // 锁定DIB
  496. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  497. // 找到DIB图像象素起始位置
  498. lpDIBBits = ::FindDIBBits(lpDIB);
  499. unsigned long biWidth =::DIBWidth(lpDIB);
  500. unsigned long biHeight = ::DIBHeight(lpDIB);
  501. unsigned long biAlign = (biWidth*3+3)/4 *4;
  502. unsigned long bmSize = biHeight * biAlign;
  503. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  506. //图像矩阵坐标与像素数据
  507. int x,y,cur;
  508. unsigned char tempR, tempG, tempB;
  509. float fTempBufforDisp;
  510. short MaxPixVal,MinPixVal,Diff;
  511. short **spOriginData, **spTransData1;
  512. //分配图像小波变换所需的数据空间
  513. spOriginData = new short* [biHeight];
  514. spTransData1 = new short* [biHeight];
  516. for(int i = 0; i < biHeight; i ++)
  517. {
  518. spOriginData[i] = new short [biWidth];
  519. spTransData1[i] = new short [biWidth];
  520. }
  521. //创建图像小波变换类
  522. CWvltTrans *pTrans;
  523. //从设备缓存中获取原始图像数据
  524. for(y=0; y<(int)biHeight; y++)
  525. {
  526. for( x=0; x<(int)biWidth; x++)
  527. {
  528. cur = y*biAlign+3*x;
  529. tempB=lpData[cur];
  530. tempG=lpData[cur+1];
  531. tempR=lpData[cur+2];
  532. spOriginData[biHeight-1-y][x]=(short)(0.3*tempR+0.59*tempG+0.11*tempB);  //no problem
  533. }
  534. }
  535. //完成一次竖直方向的图像小波变换
  536. pTrans->Vert_Transform(spOriginData,spTransData1,biHeight/2,biWidth,1);
  537. //屏蔽图像复原操作标志
  538. m_bFilter = TRUE;
  539. m_bFilter = m_bFilter & ~m_bOnce & ~m_bTwice & ~m_bTribl;
  540. MaxPixVal=spTransData1[0][0];
  541. MinPixVal=spTransData1[0][0];
  542. //计算得到图像小波系数的极大值和极小值
  543. for( y=0; y<(int)biHeight; y++)
  544. {
  545. for( x=0; x<(int)biWidth; x++)
  546. {
  547. if(MaxPixVal<spTransData1[y][x])
  548. MaxPixVal=spTransData1[y][x];
  549. if(MinPixVal>spTransData1[y][x])
  550. MinPixVal=spTransData1[y][x];
  551. }
  552. }
  553. //计算小波系数的极值差
  554. Diff=MaxPixVal-MinPixVal;
  555. //将小波系数处理后，放入显示缓存中
  556. for(y=0; y<(int)biHeight; y++)
  557. {
  558. for(x=0; x<(int)biWidth; x++)
  559. {
  560. //因为小波变换后的小波系数有可能超过255甚至更多，那么就将
  561. //小波系数的范围映射到0~255区间内，以后出现类似的处理，目的都是一样的
  562. fTempBufforDisp=spTransData1[biHeight-1-y][x];
  563. fTempBufforDisp-=MinPixVal;
  564. fTempBufforDisp*=255;
  565. fTempBufforDisp/=Diff;
  566. cur= y*biAlign+3*x; //current pixel
  567. lpData[cur] = (unsigned char)fTempBufforDisp;
  568. lpData[cur+1]= (unsigned char)fTempBufforDisp;
  569. lpData[cur+2]= (unsigned char)fTempBufforDisp;
  570. }
  571. }
  572. SetModifiedFlag(TRUE);
  573. //显示图像的小波系数
  574. UpdateAllViews(NULL);
  575. ::GlobalUnlock((HGLOBAL)GetHDIB());
  576. //删除临时的数据空间
  577. delete spOriginData;
  578. delete spTransData1;
  580. }
  582. void CMyDoc::OnWvltTransOnce() 
  583. {
  584. // TODO: Add your command handler code here
  585. //读取数字图像的文件头，获取图像的属性参数
  586. LPSTR lpDIB;
  587. // 指向DIB象素指针
  588. LPSTR    lpDIBBits;
  589. // 锁定DIB
  590. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  591. // 找到DIB图像象素起始位置
  592. lpDIBBits = ::FindDIBBits(lpDIB);
  593. unsigned long biWidth =::DIBWidth(lpDIB);
  594. unsigned long biHeight = ::DIBHeight(lpDIB);
  595. unsigned long biAlign = (biWidth*3+3)/4 *4;
  596. unsigned long bmSize = biHeight * biAlign;
  597. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  600. //图像矩阵坐标与像素数据
  601. int x,y,cur;
  602. unsigned char tempR, tempG, tempB;
  603. float fTempBufforDisp;
  604. short MaxPixVal,MinPixVal,Diff;
  605. short **spOriginData, **spTransData0, **spTransData1;
  606. //分配图像小波变换所用的数据空间
  607. spOriginData = new short* [biHeight];
  608. spTransData0 = new short* [biHeight];
  609. spTransData1 = new short* [biHeight];
  610. m_WvltCoeff = new short * [biHeight];
  612. for(int i = 0; i < biHeight; i ++)
  613. {
  614. spOriginData[i] = new short [biWidth];
  615. spTransData0[i] = new short [biWidth];
  616. spTransData1[i] = new short [biWidth];
  617. m_WvltCoeff[i] = new short [biWidth];
  618. }
  619. //创建图像小波变换类
  620. CWvltTrans *pTrans;
  621. //从设备缓存中获取原始图像数据
  622. for(y=0; y<(int)biHeight; y++)
  623. {
  624. for( x=0; x<(int)biWidth; x++)
  625. {
  626. cur = y*biAlign+3*x;
  627. tempB=lpData[cur];
  628. tempG=lpData[cur+1];
  629. tempR=lpData[cur+2];
  630. spOriginData[biHeight-1-y][x]=(short)(0.3*tempR+0.59*tempG+0.11*tempB);  //no problem
  631. }
  632. }
  633. //完成一次图像小波变换
  634. pTrans->DWT_Once(spOriginData,spTransData0,spTransData1,biHeight,biHeight/2,biWidth,biWidth/2,1,1.414);
  635. //允许图像复原操作标志
  636. m_bOnce = TRUE;
  637. m_bOnce = m_bOnce & ~m_bTwice & ~m_bTribl & ~m_bFilter;
  638. MaxPixVal=spTransData1[0][0];
  639. MinPixVal=spTransData1[0][0];
  640. for( y=0; y<(int)biHeight; y++)
  641. {
  642. for( x=0; x<(int)biWidth; x++)
  643. {
  644. if(MaxPixVal<spTransData1[y][x])
  645. MaxPixVal=spTransData1[y][x];
  646. if(MinPixVal>spTransData1[y][x])
  647. MinPixVal=spTransData1[y][x];
  648. m_WvltCoeff[y][x] = spTransData1[y][x];
  649. }
  650. }
  651. Diff=MaxPixVal-MinPixVal;
  652. for(y=0; y<(int)biHeight; y++)
  653. {
  654. for(x=0; x<(int)biWidth; x++)
  655. {
  656. //因为小波变换后的小波系数有可能超过255甚至更多，那么就将
  657. //小波系数的范围映射到0~255区间内，以后出现类似的处理，目的都是一样的
  658. fTempBufforDisp=spTransData1[biHeight-1-y][x];
  659. fTempBufforDisp-=MinPixVal;
  660. fTempBufforDisp*=255;
  661. fTempBufforDisp/=Diff;
  662. cur= y*biAlign+3*x; //current pixel
  663. lpData[cur] = (unsigned char)fTempBufforDisp;
  664. lpData[cur+1]= (unsigned char)fTempBufforDisp;
  665. lpData[cur+2]= (unsigned char)fTempBufforDisp;
  666. }
  667. }
  668. SetModifiedFlag(TRUE);
  669. //显示图像的小波系数
  670. UpdateAllViews(NULL);
  671. ::GlobalUnlock((HGLOBAL)GetHDIB());
  672. //删除临时的数据空间
  673. delete spOriginData;
  674. delete spTransData0;
  675. delete spTransData1;
  676. }
  678. void CMyDoc::OnWvltTransTwice() 
  679. {
  680. // TODO: Add your command handler code here
  681. //读取数字图像的文件头，获取图像的属性参数
  682. LPSTR lpDIB;
  683. // 指向DIB象素指针
  684. LPSTR    lpDIBBits;
  685. // 锁定DIB
  686. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  687. // 找到DIB图像象素起始位置
  688. lpDIBBits = ::FindDIBBits(lpDIB);
  689. unsigned long biWidth =::DIBWidth(lpDIB);
  690. unsigned long biHeight = ::DIBHeight(lpDIB);
  691. unsigned long biAlign = (biWidth*3+3)/4 *4;
  692. unsigned long bmSize = biHeight * biAlign;
  693. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  696. //图像矩阵坐标与像素数值
  697. int x,y,cur;
  698. unsigned char tempR, tempG, tempB;
  699. float fTempBufforDisp;
  700. short MaxPixVal,MinPixVal,Diff;
  701. short **spOriginData, **spTransData0, **spTransData1;
  702. //分配数据空间
  703. spOriginData = new short* [biHeight];
  704. spTransData0 = new short* [biHeight];
  705. spTransData1 = new short* [biHeight];
  706. m_WvltCoeff = new short * [biHeight];
  708. for(int i = 0; i < biHeight; i ++)
  709. {
  710. spOriginData[i] = new short [biWidth];
  711. spTransData0[i] = new short [biWidth];
  712. spTransData1[i] = new short [biWidth];
  713. m_WvltCoeff[i] = new short [biWidth];
  714. }
  715. //创建图像小波变化类
  716. CWvltTrans *pTrans;
  717. //从设备缓存中获取原始图像数据
  718. for(y=0; y<(int)biHeight; y++)
  719. {
  720. for( x=0; x<(int)biWidth; x++)
  721. {
  722. cur = y*biAlign+3*x;
  723. tempB=lpData[cur];
  724. tempG=lpData[cur+1];
  725. tempR=lpData[cur+2];
  726. spOriginData[biHeight-1-y][x]=(short)(0.3*tempR+0.59*tempG+0.11*tempB);  //no problem
  727. }
  728. }
  729. //完成图像的两次小波变换
  730. pTrans->DWT_TwoLayers(spOriginData,spTransData0,spTransData1,biHeight,biHeight/2,biWidth,biWidth/2,2,1.414);
  731. //允许图像复原操作标志
  732. m_bTwice = TRUE;
  733. m_bTwice = m_bTwice & ~m_bTribl & ~m_bOnce & ~m_bFilter;
  734. MaxPixVal=spTransData1[0][0];
  735. MinPixVal=spTransData1[0][0];
  736. //计算得到图像小波系数的极大值与极小值
  737. for( y=0; y<(int)biHeight; y++)
  738. {
  739. for( x=0; x<(int)biWidth; x++)
  740. {
  741. if(MaxPixVal<spTransData1[y][x])
  742. MaxPixVal=spTransData1[y][x];
  743. if(MinPixVal>spTransData1[y][x])
  744. MinPixVal=spTransData1[y][x];
  745. m_WvltCoeff[y][x] = spTransData1[y][x];
  746. }
  747. }
  748. //计算获得小波系数的极值差
  749. Diff=MaxPixVal-MinPixVal;
  750. //小波系数经过处理后，放入显示缓存中
  751. for(y=0; y<(int)biHeight; y++)
  752. {
  753. for(x=0; x<(int)biWidth; x++)
  754. {
  755. //因为小波变换后的小波系数有可能超过255甚至更多，那么就将
  756. //小波系数的范围映射到0~255区间内，以后出现类似的处理，目的都是一样的
  757. fTempBufforDisp=spTransData1[biHeight-1-y][x];
  758. fTempBufforDisp-=MinPixVal;
  759. fTempBufforDisp*=255;
  760. fTempBufforDisp/=Diff;
  761. cur= y*biAlign+3*x; //current pixel
  762. lpData[cur] = (unsigned char)fTempBufforDisp;
  763. lpData[cur+1]= (unsigned char)fTempBufforDisp;
  764. lpData[cur+2]= (unsigned char)fTempBufforDisp;
  765. }
  766. }
  767. SetModifiedFlag(TRUE);
  768. //显示图像
  769. UpdateAllViews(NULL);
  770. ::GlobalUnlock((HGLOBAL)GetHDIB());
  771. //删除临时的数据空间
  772. delete spOriginData;
  773. delete spTransData0;
  774. delete spTransData1;
  775. }
  777. void CMyDoc::OnWvltTransTrbl() 
  778. {
  779. // TODO: Add your command handler code here
  780. //读取数字图像的文件头，获取图像的属性参数
  781. LPSTR lpDIB;
  782. // 指向DIB象素指针
  783. LPSTR    lpDIBBits;
  784. // 锁定DIB
  785. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  786. // 找到DIB图像象素起始位置
  787. lpDIBBits = ::FindDIBBits(lpDIB);
  788. unsigned long biWidth =::DIBWidth(lpDIB);
  789. unsigned long biHeight = ::DIBHeight(lpDIB);
  790. unsigned long biAlign = (biWidth*3+3)/4 *4;
  791. unsigned long bmSize = biHeight * biAlign;
  792. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  794. //图像矩阵坐标与像素数据
  795. int x,y,cur;
  796. unsigned char tempR, tempG, tempB;
  797. float fTempBufforDisp;
  798. short MaxPixVal,MinPixVal,Diff;
  799. short **spOriginData, **spTransData0, **spTransData1;
  800. //分配图像小波变换的数据内存空间
  801. spOriginData = new short* [biHeight];
  802. spTransData0 = new short* [biHeight];
  803. spTransData1 = new short* [biHeight];
  804. m_WvltCoeff = new short * [biHeight];
  805. for(int i = 0; i < biHeight; i ++)
  806. {
  807. spOriginData[i] = new short [biWidth];
  808. spTransData0[i] = new short [biWidth];
  809. spTransData1[i] = new short [biWidth];
  810. m_WvltCoeff[i] = new short [biWidth];
  811. }
  812. //创建图像小波类
  813. CWvltTrans *pTrans;
  814. //从设备的图像缓存中获取原始图像的数据
  815. for(y=0; y<(int)biHeight; y++)
  816. {
  817. for( x=0; x<(int)biWidth; x++)
  818. {
  819. cur = y*biAlign+3*x;
  820. tempB=lpData[cur];
  821. tempG=lpData[cur+1];
  822. tempR=lpData[cur+2];
  823. spOriginData[biHeight-1-y][x]=(short)(0.3*tempR+0.59*tempG+0.11*tempB);  //no problem
  824. }
  825. }
  826. //完成图像的三次小波变换
  827. pTrans->DWT_TriLayers(spOriginData,spTransData0,spTransData1,biHeight,biHeight/2,biWidth,biWidth/2,3,1.414);
  828. //允许图像复员操作标志
  829. m_bTribl = TRUE;
  830. m_bTribl = m_bTribl & ~m_bTwice & ~m_bOnce & ~m_bFilter;
  831. MaxPixVal=spTransData1[0][0];
  832. MinPixVal=spTransData1[0][0];
  833. //得到小波系数的极大值和极小值
  834. for( y=0; y<(int)biHeight; y++)
  835. {
  836. for( x=0; x<(int)biWidth; x++)
  837. {
  838. if(MaxPixVal<spTransData1[y][x])
  839. MaxPixVal=spTransData1[y][x];
  840. if(MinPixVal>spTransData1[y][x])
  841. MinPixVal=spTransData1[y][x];
  842. m_WvltCoeff[y][x] = spTransData1[y][x];
  843. }
  844. }
  845. //计算出小波系数的极值差
  846. Diff=MaxPixVal-MinPixVal;
  847. //将图像的小波数据处理后放入显示缓存中
  848. for(y=0; y<(int)biHeight; y++)
  849. {
  850. for(x=0; x<(int)biWidth; x++)
  851. {
  852. //因为小波变换后的小波系数有可能超过255甚至更多，那么就将
  853. //小波系数的范围映射到0~255区间内，以后出现类似的处理，目的都是一样的
  854. fTempBufforDisp=spTransData1[biHeight-1-y][x];
  855. fTempBufforDisp-=MinPixVal;
  856. fTempBufforDisp*=255;
  857. fTempBufforDisp/=(float)Diff;
  858. cur= y*biAlign+3*x; //current pixel
  859. lpData[cur] = (unsigned char)fTempBufforDisp;
  860. lpData[cur+1]= (unsigned char)fTempBufforDisp;
  861. lpData[cur+2]= (unsigned char)fTempBufforDisp;
  862. }
  863. }
  864. SetModifiedFlag(TRUE);
  865. //显示图像
  866. UpdateAllViews(NULL);
  867. ::GlobalUnlock((HGLOBAL)GetHDIB());
  868. //删除临时的数据空间
  869. delete spOriginData;
  870. delete spTransData0;
  871. delete spTransData1;
  872. }
  876. void CMyDoc::OnDiprocFusion() 
  877. {
  878. // TODO: Add your command handler code here
  879. //读取数字图像的文件头，获取图像的属性参数
  881. LPSTR lpDIB;
  882. // 指向DIB象素指针
  883. LPSTR    lpDIBBits;
  884. // 锁定DIB
  885. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  886. // 找到DIB图像象素起始位置
  887. lpDIBBits = ::FindDIBBits(lpDIB);
  888. unsigned long biWidth =::DIBWidth(lpDIB);
  889. unsigned long biHeight = ::DIBHeight(lpDIB);
  890. unsigned long biAlign = (biWidth*3+3)/4 *4;
  891. unsigned long bmSize = biHeight * biAlign;
  892. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  895. //定义图像的内存数据空间指针
  896. short **pData, **pDataFusion;
  897. CString cstrFileName = m_strFileFusion;
  898. //图像矩阵的坐标与像素数据
  899. int x, y, cur;
  900. int tempR, tempG, tempB;
  901. //分配图像的数据空间
  902. pData = new short * [biHeight];
  903. pDataFusion = new short * [biHeight];
  904. for(int i = 0; i < biWidth; i ++)
  905. {
  906. pDataFusion[i] = new short [biWidth];
  907. pData[i] = new short [biWidth];
  908. }
  909. //从设备的显示缓存中获取原始图像的数据
  910. for(y = 0; y < (int) biHeight; y ++)
  911. {
  912. for(x = 0; x < (int) biWidth; x ++)
  913. {
  914. cur = y * biAlign + 3 * x;
  915. tempB = lpData[cur];
  916. tempG = lpData[cur + 1];
  917. tempR = lpData[cur + 2];
  918. pData[biHeight - 1 - y][x] = 0.3 * tempR + 0.59 * tempG + 0.11 * tempB;
  919. }
  920. }
  921. //释放设备占用的显示缓存
  922. m_hDIB = NULL;
  923. ::GlobalUnlock((HGLOBAL)GetHDIB());
  924. //打开另外一幅用于融合的图像文件
  925. MessageBoxA(NULL,_T("请在WvltDip\Fusion目录中打开另外一幅图像，进行图像融合"),_T("Message"),MB_ICONEXCLAMATION|MB_OK);
  926. while(!abs((int)strcmp(m_strFileFusion, cstrFileName)))
  927. {
  928. CFileDialog dlg( TRUE,NULL,NULL,
  929.                  OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,
  930.                  _T("Bitmap (*.BMP)|*.BMP|"));
  931. if( dlg.DoModal()==IDOK )
  932. cstrFileName = dlg.GetPathName();
  933. //如果文件重名，提示重新打开文件
  934.     
  935. if(!abs((int)strcmp(m_strFileFusion, cstrFileName))) 
  936. MessageBoxA(NULL,_T("文件"+m_strFileFusion+"已打开，请打开该目录中另外一个文件"),_T("Message"),MB_ICONEXCLAMATION|MB_OK);
  937. else
  938. MessageBoxA(NULL,_T("用于图像融合的文件是：n"+m_strFileFusion+"n"+cstrFileName),_T("Message"),MB_ICONEXCLAMATION|MB_OK);
  939. }
  940. CFile file;
  941.     CFileException fe;
  942. if (!file.Open(cstrFileName, CFile::modeRead | CFile::shareDenyWrite, &fe))
  943.   {
  944.   // 失败
  945.   ReportSaveLoadException(cstrFileName, &fe,
  946.   FALSE, AFX_IDP_FAILED_TO_OPEN_DOC);
  947.   }
  948.  
  949.   // 更改光标形状
  950.   BeginWaitCursor();
  951.  
  952.   // 尝试调用ReadDIBFile()读取图像,读取位图数据，载入显示缓存中
  953.   TRY
  954.   {
  955.   m_hDIB = ::ReadDIBFile(file);
  956.   }
  957.   CATCH (CFileException, eLoad)
  958.   {
  959.   // 读取失败
  960.   file.Abort();
  961.  
  962.   // 恢复光标形状
  963.   EndWaitCursor();
  964.  
  965.   // 报告失败
  966.   ReportSaveLoadException(cstrFileName, eLoad,
  967.   FALSE, AFX_IDP_FAILED_TO_OPEN_DOC);
  968.  
  969.   // 设置DIB为空
  970.   m_hDIB = NULL;
  971.  
  972.   }
  973.   END_CATCH
  974.  
  975.   // 初始化DIB
  976.   InitDIBData();
  977.  
  978.   // 判断读取文件是否成功
  979.   if (m_hDIB == NULL)
  980.   {
  981.   // 失败，可能非BMP格式
  982.   CString strMsg;
  983.   strMsg = "读取图像时出错！可能是不支持该类型的图像文件！";
  984.  
  985.   // 提示出错
  986.   MessageBox(NULL, strMsg, "系统提示", MB_ICONINFORMATION | MB_OK);
  987.   
  988.   // 恢复光标形状
  989.   EndWaitCursor();
  990.   }
  991.  
  992.   // 初始化胀标记为FALSE
  993.   SetModifiedFlag(FALSE);
  994.  
  995.   //显示新打开的图像
  996.   UpdateAllViews(NULL);
  997.   
  999. // 锁定DIB
  1000. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  1001. // 找到DIB图像象素起始位置
  1002. lpDIBBits = ::FindDIBBits(lpDIB);
  1003. //读取另外一幅图像的数据
  1004. unsigned char *lpImgData = (unsigned char*)lpDIBBits ;
  1005. //从设备的显示缓存中获取原始图像的数据
  1006. for(y = 0; y < biHeight; y ++)
  1007. {
  1008. for (x = 0; x < biWidth; x++)
  1009. {
  1010. cur = y * biAlign + 3 * x;
  1011. tempB = lpImgData[cur];
  1012. tempG = lpImgData[cur + 1];
  1013. tempR = lpImgData[cur + 2];
  1014. pDataFusion[biHeight - 1 - y][x] = 0.3 * tempR + 0.59 * tempG + 0.11 * tempB;
  1015. }
  1016. }
  1018. //创建用于图像融合的Diproc类
  1019. CDiproc *pDIP;
  1020. pDIP->DIP_ImageFusion(pData, pDataFusion, biHeight, biWidth);
  1021. //将图像数据放入显示缓存中
  1022. for(y=0; y<(int)biHeight; y++)
  1023. {
  1024. for(x=0; x<(int)biWidth; x++)
  1025. {
  1026. cur= y*biAlign+3*x; //current pixel
  1027. lpImgData[cur] = (unsigned char)pDataFusion[biHeight - 1- y][x];
  1028. lpImgData[cur+1]= (unsigned char)pDataFusion[biHeight - 1 - y][x];
  1029. lpImgData[cur+2]= (unsigned char)pDataFusion[biHeight - 1 - y][x];
  1030. }
  1031. }
  1032. SetModifiedFlag(TRUE);
  1033. //显示图像
  1034. UpdateAllViews(NULL);
  1035. ::GlobalUnlock((HGLOBAL)GetHDIB());
  1036. //删除临时的数据内存空间
  1037. delete pData;
  1038. delete pDataFusion;
  1039. }
  1040. void CMyDoc::OnDiprocRever() 
  1041. {
  1042. // TODO: Add your command handler code here
  1043. //读取数字图像的文件头，获取图像的属性参数
  1045. LPSTR lpDIB;
  1046. // 指向DIB象素指针
  1047. LPSTR    lpDIBBits;
  1048. // 锁定DIB
  1049. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  1050. // 找到DIB图像象素起始位置
  1051. lpDIBBits = ::FindDIBBits(lpDIB);
  1052. unsigned long biWidth =::DIBWidth(lpDIB);
  1053. unsigned long biHeight = ::DIBHeight(lpDIB);
  1054. unsigned long biAlign = (biWidth*3+3)/4 *4;
  1055. unsigned long bmSize = biHeight * biAlign;
  1056. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  1059. //图像矩阵的坐标
  1060. int x,y,cur;
  1061. int tempR, tempG, tempB;
  1062. //图像数据与图像处理数据的存放空间
  1063. short **spOriginData, **spTransData0,**spTransData1;
  1064. spOriginData = new short* [biHeight];
  1065. spTransData0 = new short* [biHeight];
  1066. spTransData1 = new short* [biHeight];
  1067. //分配数据的内存空间
  1068. for(int i = 0; i < biHeight; i ++)
  1069. {
  1070. spOriginData[i] = new short [biWidth];
  1071. spTransData0[i] = new short [biWidth];
  1072. spTransData1[i] = new short [biWidth];
  1073. }
  1074. //创建图像小波变换类
  1075. CWvltTrans *pTrans;
  1076. //从设备缓存中获取原始图像数据
  1077. for(y=0; y<(int)biHeight; y++)
  1078. {
  1079. for( x=0; x<(int)biWidth; x++)
  1080. {
  1081. cur = y*biAlign+3*x;
  1082. tempB=lpData[cur];
  1083. tempG=lpData[cur+1];
  1084. tempR=lpData[cur+2];
  1085. spOriginData[biHeight-1-y][x]=(short)(0.3*tempR+0.59*tempG+0.11*tempB);  //no problem
  1086. }
  1087. }
  1089. //创建数字处理类
  1090. CDiproc *pDip;
  1091. //如果图像只经过一次小波变换，则进行相应的逆变换
  1092. if(m_bOnce)
  1093. { //完成一次图像小波变换
  1094.      pDip->DIP_WvltRevers(spOriginData, spTransData0, m_WvltCoeff,biHeight,biHeight/2,biWidth,biWidth/2,1,1.414);
  1095. m_bOnce = FALSE;
  1096. delete m_WvltCoeff;
  1097. }
  1098. //如果图像只经过两次小波变换，则进行相应的逆变换
  1099. else if(m_bTwice)
  1100. {
  1101. pDip->DIP_WvltRevers(spOriginData, spTransData0, m_WvltCoeff,biHeight ,biHeight/2,biWidth,biWidth/2,2,1.414);
  1102. m_bTwice = FALSE;
  1103. delete m_WvltCoeff;
  1104. }
  1105. //如果图像只经过三次小波变换，则进行相应的逆变换
  1106. else if(m_bTribl)
  1107. {
  1108. pDip->DIP_WvltRevers(spOriginData, spTransData0, m_WvltCoeff, biHeight, biHeight/2,biWidth,biWidth/2,3,1.414);
  1109. m_bTribl = FALSE;
  1110. delete m_WvltCoeff;
  1111. }
  1112. //判断是否有小波系数可以进行图像的复原
  1113. else if(m_bFilter)
  1114. {
  1115. MessageBoxA(NULL,_T("没有找到可用小波系数，或者是未做小波变换"),_T("Message"),MB_ICONEXCLAMATION|MB_OK);
  1116. m_bFilter = FALSE;
  1117. }
  1118. //将图像数据放入显示缓存中
  1119. for(y=0; y<(int)biHeight; y++)
  1120. {
  1121. for(x=0; x<(int)biWidth; x++)
  1122. {
  1123. cur= y*biAlign+3*x; //current pixel
  1124. lpData[cur] = (unsigned char)spOriginData[biHeight - 1- y][x];
  1125. lpData[cur+1]= (unsigned char)spOriginData[biHeight - 1 - y][x];
  1126. lpData[cur+2]= (unsigned char)spOriginData[biHeight - 1 - y][x];
  1127. }
  1128. }
  1129. SetModifiedFlag(TRUE);
  1130. //显示图像
  1131. UpdateAllViews(NULL);
  1132. ::GlobalUnlock((HGLOBAL)GetHDIB());
  1133. //删除临时的数据内存空间
  1134. delete spOriginData;
  1135. delete spTransData0;
  1136. //delete spTransData1;
  1137. }
  1139. void CMyDoc::OnFilterBlur() 
  1140. {
  1141. // TODO: Add your command handler code here
  1142. //读取数字图像的文件头，获取图像的属性参数
  1144. LPSTR lpDIB;
  1145. // 指向DIB象素指针
  1146. LPSTR    lpDIBBits;
  1147. // 锁定DIB
  1148. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  1149. // 找到DIB图像象素起始位置
  1150. lpDIBBits = ::FindDIBBits(lpDIB);
  1151. unsigned long biWidth =::DIBWidth(lpDIB);
  1152. unsigned long biHeight = ::DIBHeight(lpDIB);
  1153. unsigned long biAlign = (biWidth*3+3)/4 *4;
  1154. unsigned long bmSize = biHeight * biAlign;
  1155. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  1158. //图像小波系数的低通滤波层数为3层
  1159. int MaxLayer = 3;
  1160. short **ImgData;
  1161. float fTempBufforDisp;
  1162. //图像矩阵坐标与图像数据
  1163. int tempR, tempG, tempB, x, y, cur;
  1164. //分配图像数据的内存空间
  1165. ImgData = new short * [biHeight];
  1166. for(int i = 0; i < biHeight; i ++)
  1167. {
  1168. ImgData[i] = new short [biWidth];
  1169. }
  1170. //获取显示缓存中的原始图像数据
  1171. for(y=0; y<(int)biHeight; y++)
  1172. {
  1173. for(x=0; x<(int)biWidth; x++)
  1174. {
  1175. cur = y*biAlign+3*x;
  1176. tempB=lpData[cur];
  1177. tempG=lpData[cur+1];
  1178. tempR=lpData[cur+2];
  1179. ImgData[biHeight-1-y][x]=(short)(0.3*tempR+0.59*tempG+0.11*tempB);  //no problem
  1180. }
  1181. }
  1182. //创建小波滤波器类
  1183. CWFilter Filter;
  1184. //小波低通滤波处理
  1185. Filter.LPass_Filter(ImgData, biHeight, biWidth, MaxLayer, 1);
  1186. //屏蔽图像复原操作标志
  1187. m_bFilter = TRUE;
  1188. //将处理后的图像数据放入显示缓存中
  1189. for(y=0; y<(int)biHeight; y++)
  1190. {
  1191. for(x=0; x<(int)biWidth; x++)
  1192. {
  1193. cur= y*biAlign+3*x; //当前像素的位置
  1194. lpData[cur] = (unsigned char)Filter.m_FilterData[biHeight - 1- y][x];
  1195. lpData[cur+1]= (unsigned char)Filter.m_FilterData[biHeight - 1 - y][x];
  1196. lpData[cur+2]= (unsigned char)Filter.m_FilterData[biHeight - 1 - y][x];
  1197. }
  1198. }
  1199. SetModifiedFlag(TRUE);
  1200. //显示图像
  1201. UpdateAllViews(NULL);
  1202. //删除临时的图像数据空间
  1203. ::GlobalUnlock((HGLOBAL)GetHDIB());
  1204. //删除临时的图像数据空间
  1205. delete ImgData;
  1206. }
  1208. void CMyDoc::OnFilterBlur2() 
  1209. {
  1210. // TODO: Add your command handler code here
  1211. //读取数字图像的文件头，获取图像的属性参数
  1213. LPSTR lpDIB;
  1214. // 指向DIB象素指针
  1215. LPSTR    lpDIBBits;
  1216. // 锁定DIB
  1217. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  1218. // 找到DIB图像象素起始位置
  1219. lpDIBBits = ::FindDIBBits(lpDIB);
  1220. unsigned long biWidth =::DIBWidth(lpDIB);
  1221. unsigned long biHeight = ::DIBHeight(lpDIB);
  1222. unsigned long biAlign = (biWidth*3+3)/4 *4;
  1223. unsigned long bmSize = biHeight * biAlign;
  1224. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  1227. //图像小波系数的低通滤波层数为3层
  1228. int MaxLayer = 3;
  1229. short **ImgData;
  1230. float fTempBufforDisp;
  1231. //图像矩阵坐标与图像数据
  1232. int tempR, tempG, tempB, x, y, cur;
  1233. //分配图像数据的内存空间
  1234. ImgData = new short * [biHeight];
  1235. for(int i = 0; i < biHeight; i ++)
  1236. {
  1237. ImgData[i] = new short [biWidth];
  1238. }
  1239. //获取显示缓存中的原始图像数据
  1240. for(y=0; y<(int)biHeight; y++)
  1241. {
  1242. for(x=0; x<(int)biWidth; x++)
  1243. {
  1244. cur = y*biAlign+3*x;
  1245. tempB=lpData[cur];
  1246. tempG=lpData[cur+1];
  1247. tempR=lpData[cur+2];
  1248. ImgData[biHeight-1-y][x]=(short)(0.3*tempR+0.59*tempG+0.11*tempB);  //no problem
  1249. }
  1250. }
  1251. //创建小波滤波器类
  1252. CWFilter Filter;
  1253. //小波低通滤波处理
  1254. Filter.LPass_Filter(ImgData, biHeight, biWidth, MaxLayer, 2);
  1255. //屏蔽图像复原操作标志
  1256. m_bFilter = TRUE;
  1257. //将处理后的图像数据放入显示缓存中
  1258. for(y=0; y<(int)biHeight; y++)
  1259. {
  1260. for(x=0; x<(int)biWidth; x++)
  1261. {
  1262. cur= y*biAlign+3*x; //当前像素的位置
  1263. lpData[cur] = (unsigned char)Filter.m_FilterData[biHeight - 1- y][x];
  1264. lpData[cur+1]= (unsigned char)Filter.m_FilterData[biHeight - 1 - y][x];
  1265. lpData[cur+2]= (unsigned char)Filter.m_FilterData[biHeight - 1 - y][x];
  1266. }
  1267. }
  1268. SetModifiedFlag(TRUE);
  1269. //显示图像
  1270. UpdateAllViews(NULL);
  1271. //删除临时的图像数据空间
  1272. ::GlobalUnlock((HGLOBAL)GetHDIB());
  1273. delete ImgData;
  1275. }
  1277. void CMyDoc::OnFilterSharpness() 
  1278. {
  1279. // TODO: Add your command handler code here
  1280. //读取数字图像的文件头，获取图像的属性参数
  1282. LPSTR lpDIB;
  1283. // 指向DIB象素指针
  1284. LPSTR    lpDIBBits;
  1285. // 锁定DIB
  1286. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  1287. // 找到DIB图像象素起始位置
  1288. lpDIBBits = ::FindDIBBits(lpDIB);
  1289. unsigned long biWidth =::DIBWidth(lpDIB);
  1290. unsigned long biHeight = ::DIBHeight(lpDIB);
  1291. unsigned long biAlign = (biWidth*3+3)/4 *4;
  1292. unsigned long bmSize = biHeight * biAlign;
  1293. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  1295. //图像小波系数的高通滤波层数为3层
  1296. int MaxLayer = 3;
  1297. short **ImgData;
  1298. float fTempBufforDisp;
  1299. //图像矩阵坐标与图像数据
  1300. int tempR, tempG, tempB, x, y, cur;
  1301. //分配图像数据的内存空间
  1302. ImgData = new short * [biHeight];
  1303. for(int i = 0; i < biHeight; i ++)
  1304. {
  1305. ImgData[i] = new short [biWidth];
  1306. }
  1307. //获取显示缓存中的原始图像数据
  1308. for(y=0; y<(int)biHeight; y++)
  1309. {
  1310. for(x=0; x<(int)biWidth; x++)
  1311. {
  1312. cur = y*biAlign+3*x;
  1313. tempB=lpData[cur];
  1314. tempG=lpData[cur+1];
  1315. tempR=lpData[cur+2];
  1316. ImgData[biHeight-1-y][x]=(short)(0.3*tempR+0.59*tempG+0.11*tempB);  //no problem
  1317. }
  1318. }
  1319. //创建小波滤波器类
  1320. CWFilter Filter;
  1321. //小波高通滤波处理
  1322. Filter.HPass_Filter(ImgData, biHeight, biWidth, MaxLayer);
  1323. //屏蔽图像复原操作标志
  1324. m_bFilter = TRUE;
  1325. //将处理后的图像数据放入显示缓存中
  1326. for(y=0; y<(int)biHeight; y++)
  1327. {
  1328. for(x=0; x<(int)biWidth; x++)
  1329. {
  1330. cur= y*biAlign+3*x; //当前像素的位置
  1331. lpData[cur] = (unsigned char)Filter.m_FilterData[biHeight - 1- y][x];
  1332. lpData[cur+1]= (unsigned char)Filter.m_FilterData[biHeight - 1 - y][x];
  1333. lpData[cur+2]= (unsigned char)Filter.m_FilterData[biHeight - 1 - y][x];
  1334. }
  1335. }
  1336. SetModifiedFlag(TRUE);
  1337. //显示图像
  1338. UpdateAllViews(NULL);
  1339. //删除临时的图像数据空间
  1340. ::GlobalUnlock((HGLOBAL)GetHDIB());
  1341. //删除临时的图像数据空间
  1342. delete ImgData;
  1343. }
  1349. void CMyDoc::OnFilterSharpness2() 
  1350. {
  1351. // TODO: Add your command handler code here
  1352. //读取数字图像的文件头，获取图像的属性参数
  1355. LPSTR lpDIB;
  1356. // 指向DIB象素指针
  1357. LPSTR    lpDIBBits;
  1358. // 锁定DIB
  1359. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)GetHDIB());
  1360. // 找到DIB图像象素起始位置
  1361. lpDIBBits = ::FindDIBBits(lpDIB);
  1362. unsigned long biWidth =::DIBWidth(lpDIB);
  1363. unsigned long biHeight = ::DIBHeight(lpDIB);
  1364. unsigned long biAlign = (biWidth*3+3)/4 *4;
  1365. unsigned long bmSize = biHeight * biAlign;
  1366. unsigned char *lpData = (unsigned char*)lpDIBBits ;
  1370.     //图像小波系数的高通滤波层数为3层
  1371. int MaxLayer = 3;
  1372. short **ImgData;
  1373. float fTempBufforDisp;
  1374. //图像矩阵坐标与图像数据
  1375. int tempR, tempG, tempB, x, y, cur;
  1376. //分配图像数据的内存空间
  1377. ImgData = new short * [biHeight];
  1378. for(int i = 0; i < biHeight; i ++)
  1379. {
  1380. ImgData[i] = new short [biWidth];
  1381. }
  1382. //获取显示缓存中的原始图像数据
  1383. for(y=0; y<(int)biHeight; y++)
  1384. {
  1385. for(x=0; x<(int)biWidth; x++)
  1386. {
  1387. cur = y*biAlign+3*x;
  1388. tempB=lpData[cur];
  1389. tempG=lpData[cur+1];
  1390. tempR=lpData[cur+2];
  1391. ImgData[biHeight-1-y][x]=(short)(0.3*tempR+0.59*tempG+0.11*tempB);  //no problem
  1392. }
  1393. }
  1394. //创建小波滤波器类
  1395. CWFilter Filter;
  1396. //小波高通滤波处理
  1397. Filter.HPass_Filter2(ImgData, biHeight, biWidth, MaxLayer);
  1398. //屏蔽图像复原操作标志
  1399. m_bFilter = TRUE;
  1400. //将处理后的图像数据放入显示缓存中
  1401. for(y=0; y<(int)biHeight; y++)
  1402. {
  1403. for(x=0; x<(int)biWidth; x++)
  1404. {
  1405. cur= y*biAlign+3*x; //当前像素的位置
  1406. lpData[cur] = (unsigned char)Filter.m_FilterData[biHeight - 1- y][x];
  1407. lpData[cur+1]= (unsigned char)Filter.m_FilterData[biHeight - 1 - y][x];
  1408. lpData[cur+2]= (unsigned char)Filter.m_FilterData[biHeight - 1 - y][x];
  1409. }
  1410. }
  1411. SetModifiedFlag(TRUE);
  1412. //显示图像
  1413. UpdateAllViews(NULL);
  1414. //删除临时的图像数据空间
  1415. ::GlobalUnlock((HGLOBAL)GetHDIB());
  1416. delete ImgData;
  1417. }
  1422. void CMyDoc::Openanother()
  1423. {CString cstrFileName = m_strFileFusion;
  1425. //释放设备占用的显示缓存
  1427. ::GlobalUnlock((HGLOBAL)GetHDIB());
  1428. //打开另外一幅用于融合的图像文件
  1429. MessageBoxA(NULL,_T("请在WvltDip\Fusion目录中打开另外一幅图像，进行图像融合"),_T("Message"),MB_ICONEXCLAMATION|MB_OK);
  1430. while(!abs((int)strcmp(m_strFileFusion, cstrFileName)))
  1431. {
  1432. CFileDialog dlg( TRUE,NULL,NULL,
  1433.                  OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,
  1434.                  _T("Bitmap (*.BMP)|*.BMP|"));
  1435. if( dlg.DoModal()==IDOK )
  1436. cstrFileName = dlg.GetPathName();
  1437. //如果文件重名，提示重新打开文件
  1438.     
  1439. if(!abs((int)strcmp(m_strFileFusion, cstrFileName))) 
  1440. MessageBoxA(NULL,_T("文件"+m_strFileFusion+"已打开，请打开该目录中另外一个文件"),_T("Message"),MB_ICONEXCLAMATION|MB_OK);
  1441. else
  1442. MessageBoxA(NULL,_T("用于图像融合的文件是：n"+m_strFileFusion+"n"+cstrFileName),_T("Message"),MB_ICONEXCLAMATION|MB_OK);
  1443. }
  1444. CFile file;
  1445.     CFileException fe;
  1446. if (!file.Open(cstrFileName, CFile::modeRead | CFile::shareDenyWrite, &fe))
  1447.   {
  1448.   // 失败
  1449.   ReportSaveLoadException(cstrFileName, &fe,
  1450.   FALSE, AFX_IDP_FAILED_TO_OPEN_DOC);
  1451.   }
  1452.  
  1453.   // 更改光标形状
  1454.   BeginWaitCursor();
  1455.  
  1456.   // 尝试调用ReadDIBFile()读取图像,读取位图数据，载入显示缓存中
  1457.   TRY
  1458.   {
  1459.   m_hMulDIB = ::ReadDIBFile(file);
  1460.   }
  1461.   CATCH (CFileException, eLoad)
  1462.   {
  1463.   // 读取失败
  1464.   file.Abort();
  1465.  
  1466.   // 恢复光标形状
  1467.   EndWaitCursor();
  1468.  
  1469.   // 报告失败
  1470.         ReportSaveLoadException(cstrFileName, eLoad,
  1471.   FALSE, AFX_IDP_FAILED_TO_OPEN_DOC);
  1472.  
  1473.   // 设置DIB为空
  1474.   m_hMulDIB = NULL;
  1475.  
  1476.   }
  1477.   END_CATCH
  1478.  
  1479.   // 初始化DIB
  1480.   InitDIBData();
  1481.  
  1482.   // 判断读取文件是否成功
  1483.   if (m_hMulDIB == NULL)
  1484.   {
  1485.   // 失败，可能非BMP格式
  1486.   CString strMsg;
  1487.   strMsg = "读取图像时出错！可能是不支持该类型的图像文件！";
  1488.  
  1489.   // 提示出错
  1490.   MessageBox(NULL, strMsg, "系统提示", MB_ICONINFORMATION | MB_OK);
  1491.   
  1492.   // 恢复光标形状
  1493.   EndWaitCursor();
  1494.   }
  1495. }