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

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

波变换

开发平台：

Visual C++

小波变换View.cpp：源码内容
  1. // 小波变换View.cpp : implementation of the CMyView class
  2. //
  4. #include "stdafx.h"
  5. #include "小波变换.h"
  7. #include "小波变换Doc.h"
  8. #include "小波变换View.h"
  9. #include "dibapi.h"
  10. #include "math.h"
  11. #include "waveletlevel.h"
  12. #include "WaveletCoi.h"
  13. #include "Daubechiesdlg.h"
  14. #include "WaveTranform.h"
  15. #include "HebingDlg.h"
  16. #include "Symletsdlg.h"
  18. #ifdef _DEBUG
  19. #define new DEBUG_NEW
  20. #undef THIS_FILE
  21. static char THIS_FILE[] = __FILE__;
  22. #endif
  24. /////////////////////////////////////////////////////////////////////////////
  25. // CMyView
  27. IMPLEMENT_DYNCREATE(CMyView, CView)
  29. BEGIN_MESSAGE_MAP(CMyView, CView)
  30. //{{AFX_MSG_MAP(CMyView)
  31. ON_COMMAND(ID_COMPRESSION, OnCompression)
  32. ON_COMMAND(ID_DISCOMPRESSION, OnDiscompression)
  33. ON_COMMAND(ID_WAVELETCONSTRUCTION, OnWaveletconstruction)
  34. ON_COMMAND(IDM_DIPROC_FUSION2, OnDiprocFusion2)
  35. ON_COMMAND(IDM_DIPROC_FUSION3, OnDiprocFusion3)
  36. ON_COMMAND(IDM_DIPROC_FUSION4, OnDiprocFusion4)
  37.  ON_COMMAND(IDM_IMAGE_HEBING, OnImageHebing)
  38.  ON_COMMAND(ID_DETECT_TEMPLATE, OnDetectTemplate)
  39. //}}AFX_MSG_MAP
  40. // Standard printing commands
  41. ON_COMMAND(ID_FILE_PRINT, CView::OnFilePrint)
  42. ON_COMMAND(ID_FILE_PRINT_DIRECT, CView::OnFilePrint)
  43. ON_COMMAND(ID_FILE_PRINT_PREVIEW, CView::OnFilePrintPreview)
  44. END_MESSAGE_MAP()
  46. /////////////////////////////////////////////////////////////////////////////
  47. // CMyView construction/destruction
  49. CMyView::CMyView()
  50. {
  51. // TODO: add construction code here
  52. m_nLevel=1;
  53. FWidth=0;
  54. wavefunction=-1;
  55. }
  57. CMyView::~CMyView()
  58. {
  59. }
  61. BOOL CMyView::PreCreateWindow(CREATESTRUCT& cs)
  62. {
  63. // TODO: Modify the Window class or styles here by modifying
  64. //  the CREATESTRUCT cs
  66. return CView::PreCreateWindow(cs);
  67. }
  69. /////////////////////////////////////////////////////////////////////////////
  70. // CMyView drawing
  72. void CMyView::OnDraw(CDC* pDC)
  73. {
  74. BeginWaitCursor();
  75. CMyDoc* pDoc = GetDocument();
  76. ASSERT_VALID(pDoc);
  77. // TODO: add draw code for native data here
  78. // 获取DIB
  79. HDIB hDIB = pDoc->GetHDIB();
  80.     HDIB m_hMulDIB=pDoc->m_hMulDIB;
  82. // 判断DIB是否为空
  83. if (hDIB != NULL)
  84. {
  85. LPSTR lpDIB = (LPSTR) ::GlobalLock((HGLOBAL) hDIB);
  86. // 获取DIB宽度
  87. int cxDIB = (int) ::DIBWidth(lpDIB);
  88. // 获取DIB高度
  89. int cyDIB = (int) ::DIBHeight(lpDIB);
  90. ::GlobalUnlock((HGLOBAL) hDIB);
  91. CRect rcDIB;
  92. rcDIB.top = rcDIB.left = 0;
  93. rcDIB.right = cxDIB;
  94. rcDIB.bottom = cyDIB;
  95. // 输出DIB
  96. ::PaintDIB(pDC->m_hDC, &rcDIB, pDoc->GetHDIB(),
  97. &rcDIB, pDoc->GetDocPalette());
  98. }
  100. if (m_hMulDIB!= NULL)
  101. {
  102. LPSTR lpDIB = (LPSTR) ::GlobalLock((HGLOBAL)m_hMulDIB);
  103. // 获取DIB宽度
  104. int cxDIB = (int) ::DIBWidth(lpDIB);
  105. // 获取DIB高度
  106. int cyDIB = (int) ::DIBHeight(lpDIB);
  107. ::GlobalUnlock((HGLOBAL)m_hMulDIB);
  108. CRect rcDIB;
  109. rcDIB.top =0;
  110. rcDIB.left =600;
  111. rcDIB.right = cxDIB;
  112. rcDIB.bottom = cyDIB;
  113. // 输出DIB
  114. ::PaintDIB(pDC->m_hDC, &rcDIB, pDoc->m_hMulDIB,
  115. &rcDIB, pDoc->GetDocPalette());
  116. }
  117. // 恢复正常光标
  119. }
  121. /////////////////////////////////////////////////////////////////////////////
  122. // CMyView printing
  124. BOOL CMyView::OnPreparePrinting(CPrintInfo* pInfo)
  125. {
  126. // default preparation
  127. return DoPreparePrinting(pInfo);
  128. }
  130. void CMyView::OnBeginPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
  131. {
  132. // TODO: add extra initialization before printing
  133. }
  135. void CMyView::OnEndPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
  136. {
  137. // TODO: add cleanup after printing
  138. }
  140. /////////////////////////////////////////////////////////////////////////////
  141. // CMyView diagnostics
  143. #ifdef _DEBUG
  144. void CMyView::AssertValid() const
  145. {
  146. CView::AssertValid();
  147. }
  149. void CMyView::Dump(CDumpContext& dc) const
  150. {
  151. CView::Dump(dc);
  152. }
  154. CMyDoc* CMyView::GetDocument() // non-debug version is inline
  155. {
  156. ASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CMyDoc)));
  157. return (CMyDoc*)m_pDocument;
  158. }
  159. #endif //_DEBUG
  161. /////////////////////////////////////////////////////////////////////////////
  162. // CMyView message handlers
  164. void CMyView::OnWaveletconstruction() 
  165. {
  166. // TODO: Add your command handler code here
  167. //所谓的构建，其实是已经将数据存入了数组里。
  168. BeginWaitCursor();
  169. CMyDoc* pDoc = GetDocument();
  170. Ld = new double[20];  //分解尺度函数
  171. Hd = new double[20];  //分解母函数
  172. Lr = new double[20];  //重建尺度函数
  173. Hr = new double[20];  //重建母函数
  174. CWaveletCoi waveletcoi;
  175. if(waveletcoi.DoModal()==IDOK)
  176. {
  177. FWidth=waveletcoi.m_waveletlength;
  178. wavefunction=waveletcoi.m_waveletfunction;
  179. switch(wavefunction)
  180. {
  181. case 0://IDC_SHANNON
  182. {
  183. Ld[0]=1/sqrt(2);
  184. Ld[1]=1/sqrt(2);
  185. break;
  186. }
  187. case 1://IDC_DAUBECHIES
  188. {
  189. switch(FWidth)
  190. {
  191. case 2:
  192. Ld[0]=1/sqrt(2); Ld[1]=1/sqrt(2);
  193. m_wavetransform.m_preoffset=0;
  194. m_wavetransform.m_aftoffset=1;
  195. break;
  196. case 4:
  197. Ld[0]=0.4829629131445341; Ld[1]=0.8365163037378077; Ld[2]=0.2241438680420134;
  198. Ld[3]=-0.1294095225512603;
  199. m_wavetransform.m_preoffset=0;
  200. m_wavetransform.m_aftoffset=3;
  201. break;
  202. case 6:
  203. Ld[0]=0.3326705529500825; Ld[1]=0.8068915093110924; Ld[2]=0.4598775021184914;
  204. Ld[3]=-0.1350110200102546; Ld[4]=-0.0854412738820267; Ld[5]=0.035226218857095;
  205. m_wavetransform.m_preoffset=1;
  206. m_wavetransform.m_aftoffset=4;
  207. break;
  208. case 8:
  209. Ld[0]=0.230377813309; Ld[1]=0.714846570553; Ld[2]=0.630880767930;
  210. Ld[3]=-0.027983769417; Ld[4]=-0.187034811719; Ld[5]=0.030841381836;
  211. Ld[6]=0.032883011667; Ld[7]=-0.010597401785;
  212. m_wavetransform.m_preoffset=1;
  213. m_wavetransform.m_aftoffset=6;
  214. break;
  215. case 10:
  216. Ld[0]=0.160102397974; Ld[1]=0.603829269797; Ld[2]=0.724308528438; 
  217. Ld[3]=0.138428145901; Ld[4]=-0.242294887066; Ld[5]=-0.032244869585;
  218. Ld[6]=0.077571493840; Ld[7]=-0.006241490213; Ld[8]=-0.012580751999;
  219. Ld[9]=0.003335725285;
  220. m_wavetransform.m_preoffset=1;
  221. m_wavetransform.m_aftoffset=8;
  222. break;
  223. case 12:
  224. Ld[0]=0.111540743350;   Ld[1]=0.494623890398; Ld[2]=0.751133908021; 
  225. Ld[3]=0.315250351709; Ld[4]=-0.226264693965; Ld[5]=-0.129766867567;
  226. Ld[6]=0.097501605587; Ld[7]=0.027522865530; Ld[8]=-0.031582039318;
  227. Ld[9]=0.000553842201; Ld[10]=0.004777257511; Ld[11]=-0.001077301085;
  228. m_wavetransform.m_preoffset=1;
  229. m_wavetransform.m_aftoffset=10;
  230. break;
  231. case 14:
  232. Ld[0]=0.077852054085; Ld[1]=0.396539319482; Ld[2]=0.729132090846;
  233. Ld[3]=0.469782287405; Ld[4]=-0.143906003929; Ld[5]=-0.224036184994;
  234. Ld[6]=0.071309219267; Ld[7]=0.080612609151; Ld[8]=-0.038029936935;
  235. Ld[9]=-0.016574541631; Ld[10]=0.012550998556; Ld[11]=0.000429577973;
  236. Ld[12]=-0.001801640704; Ld[13]=0.000353713800;
  237. m_wavetransform.m_preoffset=1;
  238. m_wavetransform.m_aftoffset=12;
  239. break;
  240. case 16:
  241. Ld[0]=0.054415842243; Ld[1]=0.312871590914; Ld[2]=0.675630736297;
  242. Ld[3]=0.585354683654; Ld[4]=-0.015829105256; Ld[5]=-0.284015542962;
  243. Ld[6]=0.000472484574; Ld[7]=0.128747426620; Ld[8]=-0.017369301002;
  244. Ld[9]=-0.044088253931; Ld[10]=0.013981027917; Ld[11]=0.008746094047;
  245. Ld[12]=-0.004870352993; Ld[13]=-0.000391740373; Ld[14]=0.000675449406;
  246. Ld[15]=-0.000117476784;
  247. m_wavetransform.m_preoffset=1;
  248. m_wavetransform.m_aftoffset=14;
  249. break;
  250. case 18:
  251. Ld[0]=0.038077947364; Ld[1]=0.243834674613; Ld[2]=0.604823123690;
  252. Ld[3]=0.657288078051; Ld[4]=0.133197385825; Ld[5]=-0.293273783279;
  253. Ld[6]=-0.096840783223; Ld[7]=0.148540749338; Ld[8]=0.030725681479;
  254. Ld[9]=-0.067632829061; Ld[10]=0.000250947115; Ld[11]=0.022361662124;
  255. Ld[12]=-0.004723204758; Ld[13]=-0.004281503682; Ld[14]=0.001847646883;
  256. Ld[15]=0.000230385764; Ld[16]=-0.000251963189; Ld[17]=0.000039347320;
  257. m_wavetransform.m_preoffset=2;
  258. m_wavetransform.m_aftoffset=15;
  259. break;
  260. case 20:
  261. Ld[0]=0.026670057901; Ld[1]=0.188176800078; Ld[2]=0.527201188932;
  262. Ld[3]=0.688459039454; Ld[4]=0.281172343661; Ld[5]=-0.249846424327;
  263. Ld[6]=-0.195946274377; Ld[7]=0.127369340336; Ld[8]=0.093057364604;
  264. Ld[9]=-0.071394147166; Ld[10]=-0.029457536822; Ld[11]=0.033212674059;
  265. Ld[12]=0.003606553567; Ld[13]=-0.010733175483; Ld[14]=0.001395351747;
  266. Ld[15]=0.001992405295; Ld[16]=-0.000685856695; Ld[17]=-0.000116466855;
  267. Ld[18]=0.000093588670; Ld[19]=-0.000013264203;
  268. m_wavetransform.m_preoffset=2;
  269. m_wavetransform.m_aftoffset=17;
  270. break;
  271. }
  272. }
  273. break;
  274. case 2://IDC_MORLET
  275. break;
  276. case 3://IDC_MEXICANHAT
  277. break;
  278. case 4://IDC_MEYER
  279. break;
  280. case 5://IDC_BSPLINE
  281. break;
  282. case 6://IDC_SYMLETS
  283. {
  284. switch(FWidth)
  285. {
  286. case 8:
  287. Ld[0]=-0.107148901418/sqrt(2); Ld[1]=-0.041910965126/sqrt(2); Ld[2]=0.703739068656/sqrt(2);
  288. Ld[3]=1.136658243409/sqrt(2); Ld[4]=0.421234534204/sqrt(2); Ld[5]=-0.140317624179/sqrt(2);
  289. Ld[6]=-0.017824701442/sqrt(2); Ld[7]=0.045570345896/sqrt(2);
  290. m_wavetransform.m_preoffset=3;
  291. m_wavetransform.m_aftoffset=4;
  292. break;
  293. case 10:
  294. Ld[0]=0.038654795955/sqrt(2); Ld[1]=0.041746864422/sqrt(2); Ld[2]=-0.055344186117/sqrt(2);
  295. Ld[3]=0.281990696854/sqrt(2); Ld[4]=1.023052966894/sqrt(2); Ld[5]=0.896581648380/sqrt(2);
  296. Ld[6]=0.023478923136/sqrt(2); Ld[7]=-0.247951362613/sqrt(2); Ld[8]=-0.029842499869/sqrt(2);
  297. Ld[9]=0.027632152958/sqrt(2);
  298. m_wavetransform.m_preoffset=4;
  299. m_wavetransform.m_aftoffset=5;
  300. break;
  301. case 12:
  302. Ld[0]=0.021784700327/sqrt(2); Ld[1]=0.004936612372/sqrt(2); Ld[2]=-0.166863215412/sqrt(2);
  303. Ld[3]=-0.068323121587/sqrt(2); Ld[4]=0.0694457972958/sqrt(2); Ld[5]=1.113892783926/sqrt(2);
  304. Ld[6]=0.477904371333/sqrt(2); Ld[7]=-0.102724969862/sqrt(2); Ld[8]=-0.029783751299/sqrt(2);
  305. Ld[9]=0.063250562660/sqrt(2); Ld[10]=0.002499922093/sqrt(2); Ld[11]=-0.011031867509/sqrt(2);
  306. m_wavetransform.m_preoffset=5;
  307. m_wavetransform.m_aftoffset=6;
  308. break;
  309. case 14:
  310. Ld[0]=0.003792658534/sqrt(2); Ld[1]=-0.001481225915/sqrt(2); Ld[2]=-0.017870431651/sqrt(2);
  311. Ld[3]=0.043155452582/sqrt(2); Ld[4]=0.096014767936/sqrt(2); Ld[5]=-0.070078291222/sqrt(2);
  312. Ld[6]=0.024665659489/sqrt(2); Ld[7]=0.758162601964/sqrt(2); Ld[8]=1.085782709814/sqrt(2);
  313. Ld[9]=0.408183939725/sqrt(2); Ld[10]=-0.198056706807/sqrt(2); Ld[11]=-0.152463872896/sqrt(2);
  314. Ld[12]=0.005671342686/sqrt(2); Ld[13]=0.014521394762/sqrt(2);
  315. m_wavetransform.m_preoffset=7;
  316. m_wavetransform.m_aftoffset=6;
  317. break;
  318. case 16:
  319. Ld[0]=0.002672793393/sqrt(2); Ld[1]=-0.000428394300/sqrt(2); Ld[2]=-0.021145686528/sqrt(2);
  320. Ld[3]=0.005386388754/sqrt(2); Ld[4]=0.069490465911/sqrt(2); Ld[5]=-0.038493521263/sqrt(2);
  321. Ld[6]=-0.073462508761/sqrt(2); Ld[7]=0.515398670374/sqrt(2); Ld[8]=1.099106630537/sqrt(2);
  322. Ld[9]=0.680745347190/sqrt(2); Ld[10]=-0.086653615406/sqrt(2); Ld[11]=-0.202648655286/sqrt(2);
  323. Ld[12]=0.010758611751/sqrt(2); Ld[13]=0.044823623042/sqrt(2); Ld[14]=-0.000766690896/sqrt(2);
  324. Ld[15]=-0.004783458512/sqrt(2);
  325. m_wavetransform.m_preoffset=8;
  326. m_wavetransform.m_aftoffset=7;
  327. break;
  328. case 18:
  329. Ld[0]=0.001512487309/sqrt(2); Ld[1]=-0.000669141509/sqrt(2); Ld[2]=-0.014515578553/sqrt(2);
  330. Ld[3]=0.012528896242/sqrt(2); Ld[4]=0.087791251554/sqrt(2); Ld[5]=-0.025786445930/sqrt(2);
  331. Ld[6]=-0.270893783503/sqrt(2); Ld[7]=0.049882830959/sqrt(2); Ld[8]=0.873048407349/sqrt(2);
  332. Ld[9]=1.015259790832/sqrt(2); Ld[10]=0.337658923602/sqrt(2); Ld[11]=-0.077172161097/sqrt(2);
  333. Ld[12]=0.000825140929/sqrt(2); Ld[13]=0.042744433602/sqrt(2); Ld[14]=-0.016303351226/sqrt(2);
  334. Ld[15]=-0.018769396836/sqrt(2); Ld[16]=0.000876502539/sqrt(2); Ld[17]=0.001981193736/sqrt(2);
  335. m_wavetransform.m_preoffset=8;
  336. m_wavetransform.m_aftoffset=9;
  337. break;
  338. case 20:
  339. Ld[0]=0.001089170447/sqrt(2); Ld[1]=0.000135245020/sqrt(2); Ld[2]=-0.012220642630/sqrt(2);
  340. Ld[3]=-0.002072363923/sqrt(2); Ld[4]=0.064950924579/sqrt(2); Ld[5]=0.016418869426/sqrt(2);
  341. Ld[6]=-0.225558972234/sqrt(2); Ld[7]=-0.100240215031/sqrt(2); Ld[8]=0.667071338154/sqrt(2);
  342. Ld[9]=1.088251530500/sqrt(2); Ld[10]=0.542813011213/sqrt(2); Ld[11]=-0.050256540092/sqrt(2);
  343. Ld[12]=-0.045240772218/sqrt(2); Ld[13]=0.070703567550/sqrt(2); Ld[14]=0.008152816799/sqrt(2);
  344. Ld[15]=-0.028786231926/sqrt(2); Ld[16]=-0.001137535314/sqrt(2); Ld[17]=0.006495728375/sqrt(2);
  345. Ld[18]=0.000080661204/sqrt(2); Ld[19]=-0.000649589896/sqrt(2);
  346. m_wavetransform.m_preoffset=9;
  347. m_wavetransform.m_aftoffset=10;
  348. break;
  349. }
  350. }
  351. break;
  352. case 7://IDC_COIFLETS
  353. Ld[0]=-0.0802861503271*sqrt(2); Ld[1]=-0.0243085969067*sqrt(2); Ld[2]=0.362806341593*sqrt(2);
  354. Ld[3]=0.550576616156*sqrt(2); Ld[4]=0.229036357075*sqrt(2); Ld[5]=-0.0644368523121*sqrt(2);
  355. Ld[6]=-0.0115565483406*sqrt(2); Ld[7]=0.0381688330633*sqrt(2);
  356. m_wavetransform.m_preoffset=3;
  357. m_wavetransform.m_aftoffset=4;
  358. break;
  359. }
  360. for(int i=0;i<FWidth;i++)
  361. {
  362. Hd[i]=pow(-1,i+1)*Ld[-i-1+FWidth];
  363. }
  364. for(i=0;i<FWidth;i++)
  365. {
  366. Lr[i]=Ld[-i-1+FWidth];
  367. Hr[i]=Hd[-i-1+FWidth];
  368. }
  369. }
  370. pDoc->UpdateAllViews(NULL);
  371. EndWaitCursor();
  372. }
  374. void CMyView::OnCompression() 
  375. {
  376. // TODO: Add your command handler code here
  377. BeginWaitCursor();
  378. CMyDoc* pDoc = GetDocument();
  379. LPSTR lpDIB;
  380. // 指向DIB象素指针
  381. LPSTR    lpDIBBits;
  382. // 锁定DIB
  383. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL) pDoc->GetHDIB());
  384. // 找到DIB图像象素起始位置
  385. lpDIBBits = ::FindDIBBits(lpDIB);
  386. CWaveletLevel leveldlg;
  387. if(leveldlg.DoModal()==IDOK)
  388. {
  389. m_nLevel=leveldlg.m_nlevel;
  390. //分解
  391. m_wavetransform.DIBWavelet(lpDIB,lpDIBBits,Ld,Hd,FWidth,m_nLevel);
  392. pDoc->SetModifiedFlag(TRUE);
  393. pDoc->UpdateAllViews(NULL);
  394. }
  395. // 解除锁定
  396. ::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
  397. EndWaitCursor();
  398. }
  401. void CMyView::OnDiscompression() 
  402. {
  403. // TODO: Add your command handler code here
  404. CMyDoc* pDoc = GetDocument();
  405. LPSTR lpDIB;
  406. // 指向DIB象素指针
  407. LPSTR    lpDIBBits;
  408. // 锁定DIB
  409. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL) pDoc->GetHDIB());
  410. // 找到DIB图像象素起始位置
  411. lpDIBBits = ::FindDIBBits(lpDIB);
  412. //重建
  413. m_wavetransform.DIBDisWavelet(lpDIB,lpDIBBits,Lr,Hr,FWidth,m_nLevel);
  414. pDoc->SetModifiedFlag(TRUE);
  415. pDoc->UpdateAllViews(NULL);
  416. // 解除锁定
  417. ::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
  418. }
  421. void CMyView::OnDiprocFusion2()
  422. {   
  423. // TODO: Add your command handler code here
  424. //读取数字图像的文件头，获取图像的属性参数
  425.     CMyDoc* pDoc = GetDocument();
  426. LPSTR lpDIB1;
  427. // 指向DIB象素指针
  428. LPSTR    lpDIBBits1;
  429. LPSTR lpDIB2;
  430. // 指向DIB象素指针
  431. LPSTR    lpDIBBits2;
  432. // 锁定DIB
  433. lpDIB1 = (LPSTR) ::GlobalLock((HGLOBAL)pDoc->GetHDIB());
  434. // 找到DIB图像象素起始位置
  435. lpDIBBits1 = ::FindDIBBits(lpDIB1);
  436. unsigned long lWidth =::DIBWidth(lpDIB1);
  437. unsigned long lHeight = ::DIBHeight(lpDIB1);
  438.     Dwvlttrans(lpDIB1,lpDIBBits1);
  439. pDoc->SetModifiedFlag(TRUE);
  440. pDoc->UpdateAllViews(NULL);
  443. pDoc->Openanother();
  444. pDoc->SetModifiedFlag(TRUE);
  445.   //显示新打开的图像
  446.   pDoc->UpdateAllViews(NULL);
  447. lpDIB2 = (LPSTR) ::GlobalLock((HGLOBAL)pDoc->m_hMulDIB);
  448. // 找到DIB图像象素起始位置
  449. lpDIBBits2 = ::FindDIBBits(lpDIB2);
  450.     Dwvlttrans(lpDIB2,lpDIBBits2);
  452. pDoc->SetModifiedFlag(TRUE);
  453. pDoc->UpdateAllViews(NULL);
  455.    m_wavetransform.DIP_ImageFusion2(lpDIB1,lpDIBBits1,lpDIB2,lpDIBBits2,m_nLevel);
  456. pDoc->SetModifiedFlag(TRUE);
  457. pDoc->UpdateAllViews(NULL);
  459.     ::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
  460. ::GlobalUnlock((HGLOBAL) pDoc->m_hMulDIB);
  462.     lpDIB1 = (LPSTR) ::GlobalLock((HGLOBAL)pDoc->GetHDIB() );
  463. // 找到DIB图像象素起始位置
  464. lpDIBBits1 = ::FindDIBBits(lpDIB1);
  465.     m_wavetransform.DIBDisWavelet(lpDIB1,lpDIBBits1,Lr,Hr,FWidth,m_nLevel);
  466. pDoc->SetModifiedFlag(TRUE);
  467. pDoc->UpdateAllViews(NULL);
  468. // 解除锁定
  469. ::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
  470. ::GlobalUnlock((HGLOBAL) pDoc->m_hMulDIB);
  472. }
  473. void CMyView::OnDiprocFusion3()
  474. {
  475.   CMyDoc* pDoc = GetDocument();
  476. LPSTR lpDIB1;
  477. // 指向DIB象素指针
  478. LPSTR    lpDIBBits1;
  479. LPSTR lpDIB2;
  480. // 指向DIB象素指针
  481. LPSTR    lpDIBBits2;
  482. // 锁定DIB
  483. lpDIB1 = (LPSTR) ::GlobalLock((HGLOBAL)pDoc->GetHDIB());
  484. // 找到DIB图像象素起始位置
  485. lpDIBBits1 = ::FindDIBBits(lpDIB1);
  486. unsigned long lWidth =::DIBWidth(lpDIB1);
  487. unsigned long lHeight = ::DIBHeight(lpDIB1);
  488.     Swvlttrans(lpDIB1,lpDIBBits1);
  489. pDoc->SetModifiedFlag(TRUE);
  490. pDoc->UpdateAllViews(NULL);
  493. pDoc->Openanother();
  494. pDoc->SetModifiedFlag(TRUE);
  495.   //显示新打开的图像
  496.   pDoc->UpdateAllViews(NULL);
  497. lpDIB2 = (LPSTR) ::GlobalLock((HGLOBAL)pDoc->m_hMulDIB);
  498. // 找到DIB图像象素起始位置
  499. lpDIBBits2 = ::FindDIBBits(lpDIB2);
  500.     Swvlttrans(lpDIB2,lpDIBBits2);
  502. pDoc->SetModifiedFlag(TRUE);
  503. pDoc->UpdateAllViews(NULL);
  505.    m_wavetransform.DIP_ImageFusion2(lpDIB1,lpDIBBits1,lpDIB2,lpDIBBits2,m_nLevel);
  506. pDoc->SetModifiedFlag(TRUE);
  507. pDoc->UpdateAllViews(NULL);
  509.     ::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
  510. ::GlobalUnlock((HGLOBAL) pDoc->m_hMulDIB);
  512.     lpDIB1 = (LPSTR) ::GlobalLock((HGLOBAL)pDoc->GetHDIB() );
  513. // 找到DIB图像象素起始位置
  514. lpDIBBits1 = ::FindDIBBits(lpDIB1);
  515.     m_wavetransform.DIBDisWavelet(lpDIB1,lpDIBBits1,Lr,Hr,FWidth,m_nLevel);
  516. pDoc->SetModifiedFlag(TRUE);
  517. pDoc->UpdateAllViews(NULL);
  518. // 解除锁定
  519. ::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
  520. ::GlobalUnlock((HGLOBAL) pDoc->m_hMulDIB);
  522. }
  523. void CMyView::OnDiprocFusion4()
  524. {
  525. // TODO: Add your command handler code here
  526. //读取数字图像的文件头，获取图像的属性参数
  527.    CMyDoc* pDoc = GetDocument();
  528. LPSTR lpDIB1;
  529. // 指向DIB象素指针
  530. LPSTR    lpDIBBits1;
  531. LPSTR lpDIB2;
  532. // 指向DIB象素指针
  533. LPSTR    lpDIBBits2;
  534. // 锁定DIB
  535. lpDIB1 = (LPSTR) ::GlobalLock((HGLOBAL)pDoc->GetHDIB());
  536. // 找到DIB图像象素起始位置
  537. lpDIBBits1 = ::FindDIBBits(lpDIB1);
  538. unsigned long lWidth =::DIBWidth(lpDIB1);
  539. unsigned long lHeight = ::DIBHeight(lpDIB1);
  541.     Cwvlttrans(lpDIB1,lpDIBBits1);
  542. pDoc->SetModifiedFlag(TRUE);
  543. pDoc->UpdateAllViews(NULL);
  546. pDoc->Openanother();
  547. pDoc->SetModifiedFlag(TRUE);
  548.   //显示新打开的图像
  549.   pDoc->UpdateAllViews(NULL);
  550. lpDIB2 = (LPSTR) ::GlobalLock((HGLOBAL)pDoc->m_hMulDIB);
  551. // 找到DIB图像象素起始位置
  552. lpDIBBits2 = ::FindDIBBits(lpDIB2);
  553.     Cwvlttrans(lpDIB2,lpDIBBits2);
  555. pDoc->SetModifiedFlag(TRUE);
  556. pDoc->UpdateAllViews(NULL);
  558.     m_wavetransform.DIP_ImageFusion2(lpDIB1,lpDIBBits1,lpDIB2,lpDIBBits2,m_nLevel);
  559. pDoc->SetModifiedFlag(TRUE);
  560. pDoc->UpdateAllViews(NULL);
  562.     ::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
  563. ::GlobalUnlock((HGLOBAL) pDoc->m_hMulDIB);
  565.     lpDIB1 = (LPSTR) ::GlobalLock((HGLOBAL)pDoc->GetHDIB());
  566. // 找到DIB图像象素起始位置
  567. lpDIBBits1 = ::FindDIBBits(lpDIB1);
  568.     m_wavetransform.DIBDisWavelet(lpDIB1,lpDIBBits1,Lr,Hr,FWidth,m_nLevel);
  570. pDoc->SetModifiedFlag(TRUE);
  571. pDoc->UpdateAllViews(NULL);
  572. // 解除锁定
  573. ::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
  574. ::GlobalUnlock((HGLOBAL) pDoc->m_hMulDIB);
  576. }
  577. BOOL CMyView::Dwvlttrans(LPSTR lpDIB,LPSTR lpDIBBits)
  578. {
  579. CMyDoc* pDoc = GetDocument();
  580. Ld = new double[20];  //分解尺度函数
  581. Hd = new double[20];  //分解母函数
  582. Lr = new double[20];  //重建尺度函数
  583. Hr = new double[20];  //重建母函数
  584. CDaubechiesdlg daubechiesdlg;
  585.  if(daubechiesdlg.DoModal()==IDOK)
  586. {
  587. FWidth=2*daubechiesdlg.m_nlength;
  588. switch(FWidth)
  589. {
  590. case 2:
  591. Ld[0]=1/sqrt(2); Ld[1]=1/sqrt(2);
  592. m_wavetransform.m_preoffset=0;
  593. m_wavetransform.m_aftoffset=1;
  594. break;
  595. case 4:
  596. Ld[0]=0.4829629131445341; Ld[1]=0.8365163037378077; Ld[2]=0.2241438680420134;
  597. Ld[3]=-0.1294095225512603;
  598. m_wavetransform.m_preoffset=0;
  599. m_wavetransform.m_aftoffset=3;
  600. break;
  601. case 6:
  602. Ld[0]=0.3326705529500825; Ld[1]=0.8068915093110924; Ld[2]=0.4598775021184914;
  603. Ld[3]=-0.1350110200102546; Ld[4]=-0.0854412738820267; Ld[5]=0.035226218857095;
  604. m_wavetransform.m_preoffset=1;
  605. m_wavetransform.m_aftoffset=4;
  606. break;
  607. case 8:
  608. Ld[0]=0.230377813309; Ld[1]=0.714846570553; Ld[2]=0.630880767930;
  609. Ld[3]=-0.027983769417; Ld[4]=-0.187034811719; Ld[5]=0.030841381836;
  610. Ld[6]=0.032883011667; Ld[7]=-0.010597401785;
  611. m_wavetransform.m_preoffset=1;
  612. m_wavetransform.m_aftoffset=6;
  613. break;
  614. case 10:
  615. Ld[0]=0.160102397974; Ld[1]=0.603829269797; Ld[2]=0.724308528438; 
  616. Ld[3]=0.138428145901; Ld[4]=-0.242294887066; Ld[5]=-0.032244869585;
  617. Ld[6]=0.077571493840; Ld[7]=-0.006241490213; Ld[8]=-0.012580751999;
  618. Ld[9]=0.003335725285;
  619. m_wavetransform.m_preoffset=1;
  620. m_wavetransform.m_aftoffset=8;
  621. break;
  622. case 12:
  623. Ld[0]=0.111540743350;   Ld[1]=0.494623890398; Ld[2]=0.751133908021; 
  624. Ld[3]=0.315250351709; Ld[4]=-0.226264693965; Ld[5]=-0.129766867567;
  625. Ld[6]=0.097501605587; Ld[7]=0.027522865530; Ld[8]=-0.031582039318;
  626. Ld[9]=0.000553842201; Ld[10]=0.004777257511; Ld[11]=-0.001077301085;
  627. m_wavetransform.m_preoffset=1;
  628. m_wavetransform.m_aftoffset=10;
  629. break;
  630. case 14:
  631. Ld[0]=0.077852054085; Ld[1]=0.396539319482; Ld[2]=0.729132090846;
  632. Ld[3]=0.469782287405; Ld[4]=-0.143906003929; Ld[5]=-0.224036184994;
  633. Ld[6]=0.071309219267; Ld[7]=0.080612609151; Ld[8]=-0.038029936935;
  634. Ld[9]=-0.016574541631; Ld[10]=0.012550998556; Ld[11]=0.000429577973;
  635. Ld[12]=-0.001801640704; Ld[13]=0.000353713800;
  636. m_wavetransform.m_preoffset=1;
  637. m_wavetransform.m_aftoffset=12;
  638. break;
  639. case 16:
  640. Ld[0]=0.054415842243; Ld[1]=0.312871590914; Ld[2]=0.675630736297;
  641. Ld[3]=0.585354683654; Ld[4]=-0.015829105256; Ld[5]=-0.284015542962;
  642. Ld[6]=0.000472484574; Ld[7]=0.128747426620; Ld[8]=-0.017369301002;
  643. Ld[9]=-0.044088253931; Ld[10]=0.013981027917; Ld[11]=0.008746094047;
  644. Ld[12]=-0.004870352993; Ld[13]=-0.000391740373; Ld[14]=0.000675449406;
  645. Ld[15]=-0.000117476784;
  646. m_wavetransform.m_preoffset=1;
  647. m_wavetransform.m_aftoffset=14;
  648. break;
  649. case 18:
  650. Ld[0]=0.038077947364; Ld[1]=0.243834674613; Ld[2]=0.604823123690;
  651. Ld[3]=0.657288078051; Ld[4]=0.133197385825; Ld[5]=-0.293273783279;
  652. Ld[6]=-0.096840783223; Ld[7]=0.148540749338; Ld[8]=0.030725681479;
  653. Ld[9]=-0.067632829061; Ld[10]=0.000250947115; Ld[11]=0.022361662124;
  654. Ld[12]=-0.004723204758; Ld[13]=-0.004281503682; Ld[14]=0.001847646883;
  655. Ld[15]=0.000230385764; Ld[16]=-0.000251963189; Ld[17]=0.000039347320;
  656. m_wavetransform.m_preoffset=2;
  657. m_wavetransform.m_aftoffset=15;
  658. break;
  659. case 20:
  660. Ld[0]=0.026670057901; Ld[1]=0.188176800078; Ld[2]=0.527201188932;
  661. Ld[3]=0.688459039454; Ld[4]=0.281172343661; Ld[5]=-0.249846424327;
  662. Ld[6]=-0.195946274377; Ld[7]=0.127369340336; Ld[8]=0.093057364604;
  663. Ld[9]=-0.071394147166; Ld[10]=-0.029457536822; Ld[11]=0.033212674059;
  664. Ld[12]=0.003606553567; Ld[13]=-0.010733175483; Ld[14]=0.001395351747;
  665. Ld[15]=0.001992405295; Ld[16]=-0.000685856695; Ld[17]=-0.000116466855;
  666. Ld[18]=0.000093588670; Ld[19]=-0.000013264203;
  667. m_wavetransform.m_preoffset=2;
  668. m_wavetransform.m_aftoffset=17;
  669. break;
  670. }
  672. for(int i=0;i<FWidth;i++)
  673. {
  674. Hd[i]=pow(-1,i+1)*Ld[-i-1+FWidth];
  675. }
  676. for(i=0;i<FWidth;i++)
  677. {
  678. Lr[i]=Ld[-i-1+FWidth];
  679. Hr[i]=Hd[-i-1+FWidth];
  680. }
  681. }
  682. CWaveletLevel leveldlg;
  683. if(leveldlg.DoModal()==IDOK)
  684. {
  685. m_nLevel=leveldlg.m_nlevel;
  686. //分解
  687. m_wavetransform.DIBWavelet(lpDIB,lpDIBBits,Ld,Hd,FWidth,m_nLevel);
  688.     pDoc->SetModifiedFlag(TRUE);
  689. pDoc->UpdateAllViews(NULL);
  690. }
  691. return TRUE;// 返回
  693. }
  695. BOOL CMyView::Swvlttrans(LPSTR lpDIB,LPSTR lpDIBBits)
  696. {
  697. CMyDoc* pDoc = GetDocument();
  698. Ld = new double[20];  //分解尺度函数
  699. Hd = new double[20];  //分解母函数
  700. Lr = new double[20];  //重建尺度函数
  701. Hr = new double[20];  //重建母函数
  702. CSymletsdlg symdlg;
  703. if(symdlg.DoModal()==IDOK)
  704. { FWidth=2*symdlg.m_symlength;
  705.    switch(FWidth)
  706. {
  707. case 8:
  708. Ld[0]=-0.107148901418/sqrt(2); Ld[1]=-0.041910965126/sqrt(2); Ld[2]=0.703739068656/sqrt(2);
  709. Ld[3]=1.136658243409/sqrt(2); Ld[4]=0.421234534204/sqrt(2); Ld[5]=-0.140317624179/sqrt(2);
  710. Ld[6]=-0.017824701442/sqrt(2); Ld[7]=0.045570345896/sqrt(2);
  711. m_wavetransform.m_preoffset=3;
  712. m_wavetransform.m_aftoffset=4;
  713. break;
  714. case 10:
  715. Ld[0]=0.038654795955/sqrt(2); Ld[1]=0.041746864422/sqrt(2); Ld[2]=-0.055344186117/sqrt(2);
  716. Ld[3]=0.281990696854/sqrt(2); Ld[4]=1.023052966894/sqrt(2); Ld[5]=0.896581648380/sqrt(2);
  717. Ld[6]=0.023478923136/sqrt(2); Ld[7]=-0.247951362613/sqrt(2); Ld[8]=-0.029842499869/sqrt(2);
  718. Ld[9]=0.027632152958/sqrt(2);
  719. m_wavetransform.m_preoffset=4;
  720. m_wavetransform.m_aftoffset=5;
  721. break;
  722. case 12:
  723. Ld[0]=0.021784700327/sqrt(2); Ld[1]=0.004936612372/sqrt(2); Ld[2]=-0.166863215412/sqrt(2);
  724. Ld[3]=-0.068323121587/sqrt(2); Ld[4]=0.0694457972958/sqrt(2); Ld[5]=1.113892783926/sqrt(2);
  725. Ld[6]=0.477904371333/sqrt(2); Ld[7]=-0.102724969862/sqrt(2); Ld[8]=-0.029783751299/sqrt(2);
  726. Ld[9]=0.063250562660/sqrt(2); Ld[10]=0.002499922093/sqrt(2); Ld[11]=-0.011031867509/sqrt(2);
  727. m_wavetransform.m_preoffset=5;
  728. m_wavetransform.m_aftoffset=6;
  729. break;
  730. case 14:
  731. Ld[0]=0.003792658534/sqrt(2); Ld[1]=-0.001481225915/sqrt(2); Ld[2]=-0.017870431651/sqrt(2);
  732. Ld[3]=0.043155452582/sqrt(2); Ld[4]=0.096014767936/sqrt(2); Ld[5]=-0.070078291222/sqrt(2);
  733. Ld[6]=0.024665659489/sqrt(2); Ld[7]=0.758162601964/sqrt(2); Ld[8]=1.085782709814/sqrt(2);
  734. Ld[9]=0.408183939725/sqrt(2); Ld[10]=-0.198056706807/sqrt(2); Ld[11]=-0.152463872896/sqrt(2);
  735. Ld[12]=0.005671342686/sqrt(2); Ld[13]=0.014521394762/sqrt(2);
  736. m_wavetransform.m_preoffset=7;
  737. m_wavetransform.m_aftoffset=6;
  738. break;
  739. case 16:
  740. Ld[0]=0.002672793393/sqrt(2); Ld[1]=-0.000428394300/sqrt(2); Ld[2]=-0.021145686528/sqrt(2);
  741. Ld[3]=0.005386388754/sqrt(2); Ld[4]=0.069490465911/sqrt(2); Ld[5]=-0.038493521263/sqrt(2);
  742. Ld[6]=-0.073462508761/sqrt(2); Ld[7]=0.515398670374/sqrt(2); Ld[8]=1.099106630537/sqrt(2);
  743. Ld[9]=0.680745347190/sqrt(2); Ld[10]=-0.086653615406/sqrt(2); Ld[11]=-0.202648655286/sqrt(2);
  744. Ld[12]=0.010758611751/sqrt(2); Ld[13]=0.044823623042/sqrt(2); Ld[14]=-0.000766690896/sqrt(2);
  745. Ld[15]=-0.004783458512/sqrt(2);
  746. m_wavetransform.m_preoffset=8;
  747. m_wavetransform.m_aftoffset=7;
  748. break;
  749. case 18:
  750. Ld[0]=0.001512487309/sqrt(2); Ld[1]=-0.000669141509/sqrt(2); Ld[2]=-0.014515578553/sqrt(2);
  751. Ld[3]=0.012528896242/sqrt(2); Ld[4]=0.087791251554/sqrt(2); Ld[5]=-0.025786445930/sqrt(2);
  752. Ld[6]=-0.270893783503/sqrt(2); Ld[7]=0.049882830959/sqrt(2); Ld[8]=0.873048407349/sqrt(2);
  753. Ld[9]=1.015259790832/sqrt(2); Ld[10]=0.337658923602/sqrt(2); Ld[11]=-0.077172161097/sqrt(2);
  754. Ld[12]=0.000825140929/sqrt(2); Ld[13]=0.042744433602/sqrt(2); Ld[14]=-0.016303351226/sqrt(2);
  755. Ld[15]=-0.018769396836/sqrt(2); Ld[16]=0.000876502539/sqrt(2); Ld[17]=0.001981193736/sqrt(2);
  756. m_wavetransform.m_preoffset=8;
  757. m_wavetransform.m_aftoffset=9;
  758. break;
  759. case 20:
  760. Ld[0]=0.001089170447/sqrt(2); Ld[1]=0.000135245020/sqrt(2); Ld[2]=-0.012220642630/sqrt(2);
  761. Ld[3]=-0.002072363923/sqrt(2); Ld[4]=0.064950924579/sqrt(2); Ld[5]=0.016418869426/sqrt(2);
  762. Ld[6]=-0.225558972234/sqrt(2); Ld[7]=-0.100240215031/sqrt(2); Ld[8]=0.667071338154/sqrt(2);
  763. Ld[9]=1.088251530500/sqrt(2); Ld[10]=0.542813011213/sqrt(2); Ld[11]=-0.050256540092/sqrt(2);
  764. Ld[12]=-0.045240772218/sqrt(2); Ld[13]=0.070703567550/sqrt(2); Ld[14]=0.008152816799/sqrt(2);
  765. Ld[15]=-0.028786231926/sqrt(2); Ld[16]=-0.001137535314/sqrt(2); Ld[17]=0.006495728375/sqrt(2);
  766. Ld[18]=0.000080661204/sqrt(2); Ld[19]=-0.000649589896/sqrt(2);
  767. m_wavetransform.m_preoffset=9;
  768. m_wavetransform.m_aftoffset=10;
  769. break;
  770. }
  771.        for(int i=0;i<FWidth;i++)
  772. {
  773. Hd[i]=pow(-1,i+1)*Ld[-i-1+FWidth];
  774. }
  775. for(i=0;i<FWidth;i++)
  776. {
  777. Lr[i]=Ld[-i-1+FWidth];
  778. Hr[i]=Hd[-i-1+FWidth];
  779. }
  780. }
  781. CWaveletLevel leveldlg;
  782. if(leveldlg.DoModal()==IDOK)
  783. {
  784. m_nLevel=leveldlg.m_nlevel;
  785. //分解
  786. m_wavetransform.DIBWavelet(lpDIB,lpDIBBits,Ld,Hd,FWidth,m_nLevel);
  787.     pDoc->SetModifiedFlag(TRUE);
  788. pDoc->UpdateAllViews(NULL);
  789. }
  790. return TRUE;// 返回
  791. }
  793. BOOL CMyView::Cwvlttrans(LPSTR lpDIB,LPSTR lpDIBBits)
  794. {
  795. CMyDoc* pDoc = GetDocument();
  796. Ld = new double[20];  //分解尺度函数
  797. Hd = new double[20];  //分解母函数
  798. Lr = new double[20];  //重建尺度函数
  799. Hr = new double[20];  //重建母函数
  800. FWidth=8;
  801. Ld[0]=-0.0802861503271*sqrt(2); Ld[1]=-0.0243085969067*sqrt(2); Ld[2]=0.362806341593*sqrt(2);
  802. Ld[3]=0.550576616156*sqrt(2); Ld[4]=0.229036357075*sqrt(2); Ld[5]=-0.0644368523121*sqrt(2);
  803. Ld[6]=-0.0115565483406*sqrt(2); Ld[7]=0.0381688330633*sqrt(2);
  804. m_wavetransform.m_preoffset=3;
  805. m_wavetransform.m_aftoffset=4;
  806.        for(int i=0;i<FWidth;i++)
  807. {
  808. Hd[i]=pow(-1,i+1)*Ld[-i-1+FWidth];
  809. }
  810. for(i=0;i<FWidth;i++)
  811. {
  812. Lr[i]=Ld[-i-1+FWidth];
  813. Hr[i]=Hd[-i-1+FWidth];
  814. }
  815. CWaveletLevel leveldlg;
  816. if(leveldlg.DoModal()==IDOK)
  817. {
  818. m_nLevel=leveldlg.m_nlevel;
  819. //分解
  820. m_wavetransform.DIBWavelet(lpDIB,lpDIBBits,Ld,Hd,FWidth,m_nLevel);
  821.     pDoc->SetModifiedFlag(TRUE);
  822. pDoc->UpdateAllViews(NULL);
  823. }
  825. return TRUE;// 返回
  826. }
  828. void CMyView::OnImageHebing()
  829. {
  830. CHebingDlg dlg;
  831. dlg.DoModal();
  832. }
  833. void CMyView::OnDetectTemplate() 
  834. {
  835. // 获取文档
  836. CMyDoc* pDoc = GetDocument();
  838. // 指向DIB的指针
  839. LPSTR lpDIB,lpTemplateDIB;
  841. // 指向DIB象素指针
  842. LPSTR    lpDIBBits,lpTemplateDIBBits;
  844. //图像的宽度与高度
  845. long lWidth,lHeight;
  847. //模板的宽度与高度
  848. long lTemplateWidth,lTemplateHeight;
  850. HDIB hTemplateDIB;
  852. // 锁定DIB
  853. lpDIB = (LPSTR) ::GlobalLock((HGLOBAL) pDoc->GetHDIB());
  856. // 判断是否是8-bpp位图（这里为了方便，只处理8-bpp位图的水平镜像，其它的可以类推）
  857. if (::DIBNumColors(lpDIB) != 256)
  858. {
  859. // 提示用户
  860. MessageBox("目前只支持256色位图的平移！", "系统提示" , MB_ICONINFORMATION | MB_OK);
  862. // 解除锁定
  863. ::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
  865. // 返回
  866. return;
  867. }
  870. // 更改光标形状
  871. BeginWaitCursor();
  873. // 找到DIB图像象素起始位置
  874. lpDIBBits = ::FindDIBBits(lpDIB);
  876. lWidth = ::DIBWidth(lpDIB);
  877. lHeight = ::DIBHeight(lpDIB);
  878. CFileDialog dlg(TRUE,"bmp","*.bmp");
  879. if(dlg.DoModal() == IDOK)
  880. {
  881.  
  882.   CFile file;
  883.   CFileException fe;
  884.  
  885.   CString strPathName;
  886.  
  887. strPathName = dlg.GetPathName();
  888.  
  889. // 打开文件
  890. VERIFY(file.Open(strPathName, CFile::modeRead | CFile::shareDenyWrite, &fe));
  892. // 尝试调用ReadDIBFile()读取图像
  893. TRY
  894. {
  895.   hTemplateDIB = ::ReadDIBFile(file);
  896. }
  897. CATCH (CFileException, eLoad)
  898. {
  899. // 读取失败
  900.   file.Abort();
  901.  
  902. // 恢复光标形状
  903. EndWaitCursor();
  904.  
  905. // 报告失败
  906. //ReportSaveLoadException(strPathName, eLoad,
  907. // FALSE, AFX_IDP_FAILED_TO_OPEN_DOC);
  909. // 设置DIB为空
  910. hTemplateDIB = NULL;
  912. // 返回
  913. return;
  914. }
  915. END_CATCH
  916.  
  917. // 初始化DIB
  918. //InitDIBData();
  919.  
  920. // 判断读取文件是否成功
  921. if (hTemplateDIB == NULL)
  922. {
  923.   // 失败，可能非BMP格式
  924.   CString strMsg;
  925.   strMsg = "读取图像时出错！可能是不支持该类型的图像文件！";
  926.  
  927.   // 提示出错
  928.   MessageBox(strMsg, NULL, MB_ICONINFORMATION | MB_OK);
  929.   
  930.   // 恢复光标形状
  931.   EndWaitCursor();
  932.  
  933.   // 返回
  934.   return;
  935.   }
  936.   }
  937. else
  938. {
  939.   // 恢复光标形状
  940.   EndWaitCursor();
  942. return;
  943. }
  944. // 锁定DIB
  945. lpTemplateDIB = (LPSTR) ::GlobalLock((HGLOBAL) hTemplateDIB);
  947. // 判断是否是8-bpp位图（这里为了方便，只处理8-bpp位图的模板匹配，其它的可以类推）
  948. if (::DIBNumColors(lpTemplateDIB) != 256)
  949. {
  950. // 提示用户
  951. MessageBox("目前只支持256色位图！", "系统提示" , MB_ICONINFORMATION | MB_OK);
  953. // 解除锁定
  954. ::GlobalUnlock((HGLOBAL) hTemplateDIB);
  956. // 返回
  957. return;
  958. }
  961. // 更改光标形状
  962. BeginWaitCursor();
  964. // 找到DIB图像象素起始位置
  965. lpTemplateDIBBits = ::FindDIBBits(lpTemplateDIB);
  967. lTemplateWidth = ::DIBWidth(lpTemplateDIB);
  968. lTemplateHeight = ::DIBHeight(lpTemplateDIB);
  969. if(lTemplateHeight > lHeight || lTemplateWidth > lWidth )
  970. {
  971. // 提示用户
  972. MessageBox("模板尺寸大于源图像尺寸！", "系统提示" , MB_ICONINFORMATION | MB_OK);
  974. // 解除锁定
  975. ::GlobalUnlock((HGLOBAL) hTemplateDIB);
  977. // 返回
  978. return;
  980. }
  981. // 调用TemplateMatchDIB()函数进行模板匹配
  982. if (TemplateMatchDIB(lpDIBBits,lpTemplateDIBBits, lWidth,lHeight, lTemplateWidth,lTemplateHeight))
  983. {
  985. // 设置脏标记
  986. pDoc->SetModifiedFlag(TRUE);
  988. // 更新视图
  989. pDoc->UpdateAllViews(NULL);
  990. }
  991. else
  992. {
  993. // 提示用户
  994. MessageBox("分配内存失败！", "系统提示" , MB_ICONINFORMATION | MB_OK);
  995. }
  997. // 解除锁定
  998. ::GlobalUnlock((HGLOBAL) pDoc->GetHDIB());
  999. ::GlobalUnlock((HGLOBAL) hTemplateDIB);
  1001. // 恢复光标
  1002. EndWaitCursor();
  1004. }