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

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > GDI/图象编程 > 查看源码
图像拼接技术Doc.cpp
资源名称：txpj.rar　[点击查看]
上传用户：pureled
上传日期：2013-05-27
资源大小：1078k
文件大小：21k
源码类别：

GDI/图象编程

开发平台：

Visual C++

图像拼接技术Doc.cpp：源码内容
  1. // 图像拼接技术Doc.cpp : implementation of the CMyDoc class
  2. //
  4. #include "stdafx.h"
  5. #include "图像拼接技术.h"
  6. #include "图像拼接技术Doc.h "
  7. #include "windowsx.h"
  8. #include "Start.h"
  9. #include "Other.h"
  10. #include "ZuoBiao.h"
  11. #include "Weizhi.h"
  12. #include "Wuizhi2.h"
  13. #include "Weizhi3.h"
  14. #include "RotateAngle.h"
  15. #include <math.h>
  16. #ifdef _DEBUG
  17. #define new DEBUG_NEW
  18. #define PI 3.1415926535
  19. #define RADIAN(angle) ((angle)*PI/180.0)
  20. #undef THIS_FILE
  21. static char THIS_FILE[] = __FILE__;
  22. #endif
  24. /////////////////////////////////////////////////////////////////////////////
  25. // CMyDoc
  27. IMPLEMENT_DYNCREATE(CMyDoc, CDocument)
  29. BEGIN_MESSAGE_MAP(CMyDoc, CDocument)
  30. //{{AFX_MSG_MAP(CMyDoc)
  31. ON_COMMAND(ID_FILE_OTHEROPEN, OnFileOtheropen)
  32. ON_COMMAND(ID_PINJIE, OnPinjie)
  33. ON_COMMAND(ID_SUANFA1, OnSuanfa1)
  34. ON_COMMAND(ID_SUANFA2, OnSuanfa2)
  35. ON_COMMAND(ID_ROTATEPINJIE, OnRotatepinjie)
  36. ON_COMMAND(ID_SUANFA3, OnSuanfa3)
  37. ON_COMMAND(ID_CHUIZHIPINJIE, OnChuizhipinjie)
  38. ON_COMMAND(ID_SAVEBITMAP, OnSavebitmap)
  39. //}}AFX_MSG_MAP
  40. END_MESSAGE_MAP()
  42. /////////////////////////////////////////////////////////////////////////////
  43. // CMyDoc construction/destruction
  45. CMyDoc::CMyDoc()
  46. {
  47. // TODO: add one-time construction code here
  49. }
  51. CMyDoc::~CMyDoc()
  52. {
  54. }
  56. BOOL CMyDoc::OnNewDocument()
  57. {
  58. if (!CDocument::OnNewDocument())
  59. return FALSE;
  61. // TODO: add reinitialization code here
  62. // (SDI documents will reuse this document)
  64. return TRUE;
  65. }
  69. /////////////////////////////////////////////////////////////////////////////
  70. // CMyDoc serialization
  72. void CMyDoc::Serialize(CArchive& ar)
  73. {
  74. if (ar.IsStoring())
  75. {
  76. // TODO: add storing code here
  77. }
  78. else
  79. {
  80. // TODO: add loading code here
  81. }
  82. }
  84. /////////////////////////////////////////////////////////////////////////////
  85. // CMyDoc diagnostics
  87. #ifdef _DEBUG
  88. void CMyDoc::AssertValid() const
  89. {
  90. CDocument::AssertValid();
  91. }
  93. void CMyDoc::Dump(CDumpContext& dc) const
  94. {
  95. CDocument::Dump(dc);
  96. }
  97. #endif //_DEBUG
  99. /////////////////////////////////////////////////////////////////////////////
  100. // CMyDoc commands
  107. BOOL CMyDoc::OnOpenDocument(LPCTSTR lpszPathName) 
  108. {
  109. if (!CDocument::OnOpenDocument(lpszPathName))
  110. return FALSE;
  112. // TODO: Add your specialized creation code here
  114. CFile file;
  115. LONG lBitsSize;
  116. HLOCAL hNew;
  118. if(file.Open(lpszPathName,CFile::modeRead)==NULL)
  119. {
  120. MessageBox(NULL,"打开文件失败!","提示",MB_OK);
  121. return FALSE;
  122. }
  124. file.Read(&lbmpHeader,sizeof(BITMAPFILEHEADER));
  125. if(lbmpHeader.bfType!=(WORD)(('M'<<8)|'B'))
  126. {
  127. MessageBox(NULL,"不是位图文件！","提示",MB_OK);
  128. file.Close();
  129. return FALSE;
  130. }
  132. lBitsSize=lbmpHeader.bfSize;
  134. lbmpInfo=(BITMAPINFO*)GlobalAllocPtr(GHND,lBitsSize-
  135. sizeof(BITMAPFILEHEADER));//开辟空间
  136. file.Read(lbmpInfo,lBitsSize-sizeof(BITMAPFILEHEADER));
  137. if((lbmpInfo->bmiHeader.biBitCount!=8)||(lbmpInfo->bmiHeader.biCompression!=BI_RGB))
  138. {
  139. MessageBox(NULL,"不是256色位图或是压缩的图像","提示",MB_OK);
  140. lbmpInfo=NULL;
  141. file.Close();
  142. return FALSE;
  143. }
  144. size.cx=lbmpInfo->bmiHeader.biWidth;
  145. size.cy=lbmpInfo->bmiHeader.biHeight;
  146. rgbPal=(LPBYTE)lbmpInfo+sizeof(BITMAPINFOHEADER);
  147. lpBits=(LPBYTE)lbmpInfo+sizeof(BITMAPINFOHEADER)+256*sizeof(RGBQUAD);
  148. TrueWidth=(size.cx+3)/4*4;
  149. file.Close();
  151. //获取第一图的位置
  152. Start XY;
  153. XY.DoModal();
  154. m_x=XY.m_x;
  155. m_y=XY.m_y;
  157. //分配内存
  158. hNew=LocalAlloc(LHND,size.cy*TrueWidth);
  159. Lp1=(BYTE*)LocalLock(hNew);
  161. //存放第一幅图像
  162. memcpy(Lp1,lpBits,size.cy*TrueWidth);
  164. return TRUE;
  165. }
  167. void CMyDoc::ShowBitmap(CDC *pDC)
  168. {
  169. ::StretchDIBits(pDC->m_hDC,m_x,m_y,size.cx,size.cy,0,0,
  170. size.cx,size.cy,lpBits,lbmpInfo,DIB_RGB_COLORS,SRCCOPY);
  171. }
  173. void CMyDoc::OnFileOtheropen() 
  174. {
  175. // TODO: Add your command handler code here
  177. CString FilePath;
  178. CFile file;
  179. BITMAPFILEHEADER lbmpHeader;
  180. LONG lBitsSize;
  181. HLOCAL hNew;
  182. CFileDialog dlgOpen(TRUE,".bmp",FilePath,OFN_HIDEREADONLY | OFN_OVERWRITEPROMPT,
  183. "位图文件(*.bmp)|*.bmp||");
  185. if(dlgOpen.DoModal()==IDOK)
  186. FilePath=dlgOpen.GetPathName();
  188. if(file.Open(FilePath,CFile::modeRead)==NULL)
  189. {
  190. MessageBox(NULL,"打开文件失败!","提示",MB_OK);
  191. return;
  192. }
  194. file.Read(&lbmpHeader,sizeof(BITMAPFILEHEADER));
  195. if(lbmpHeader.bfType!=(WORD)(('M'<<8)|'B'))
  196. {
  197. MessageBox(NULL,"不是位图文件！","提示",MB_OK);
  198. file.Close();
  199. return;
  200. }
  202. lBitsSize=lbmpHeader.bfSize;
  204. lbmpInfo2=(BITMAPINFO*)GlobalAllocPtr(GHND,lBitsSize-
  205. sizeof(BITMAPFILEHEADER));
  206. file.Read(lbmpInfo2,lBitsSize-sizeof(BITMAPFILEHEADER));
  207. if((lbmpInfo2->bmiHeader.biBitCount!=8)||(lbmpInfo2->bmiHeader.biCompression!=BI_RGB))
  208. {
  209. MessageBox(NULL,"不是256色位图或是压缩的图像","提示",MB_OK);
  210. lbmpInfo=NULL;
  211. file.Close();
  212. return;
  213. }
  214. size2.cx=lbmpInfo2->bmiHeader.biWidth;
  215. size2.cy=lbmpInfo2->bmiHeader.biHeight;
  216. lpBits2=(LPBYTE)lbmpInfo2+sizeof(BITMAPINFOHEADER)+256*sizeof(RGBQUAD);
  217. TrueWidth2=(size2.cx+3)/4*4;
  218. file.Close();
  220. //获取第二幅图的位置
  221. Other XY;
  222. XY.DoModal();
  223. m2_x=XY.m_x;
  224. m2_y=XY.m_y;
  225.     
  226.     // 分配内存
  227. hNew=LocalAlloc(LHND,size2.cy*TrueWidth2);
  228. Lp2=(BYTE*)LocalLock(hNew);
  230.     //存放第二幅图像
  231. memcpy(Lp2,lpBits2,size2.cy*TrueWidth2);
  232.         
  233. UpdateAllViews(NULL);
  234. }
  236. void CMyDoc::ShowOther(CDC *pDC)
  237. {
  238. ::StretchDIBits(pDC->m_hDC,m2_x,m2_y,size2.cx,size2.cy,0,0,
  239. size2.cx,size2.cy,lpBits2,lbmpInfo2,DIB_RGB_COLORS,SRCCOPY);
  240. }
  244. void CMyDoc::Moban(int a,int b,int c,int d)
  245. {
  246. HLOCAL hNew;
  247. int i;
  248. int j;
  249. m3_x=a;
  250. m3_y=b;
  251.     Height=c;
  252. Width=(d+3)/4*4;
  254. // 分配内存
  255. hNew=LocalAlloc(LHND,Height*Width);
  256. Lp3=(BYTE*)LocalLock(hNew);
  257.         
  258. for (i=0;i<Height;i++)
  259. for (j=0;j<Width;j++)
  260. {
  261. Lp3[i*Width+j]=lpBits2[(size2.cy-m3_y-Height+i)*TrueWidth2+m3_x+j];
  263. }
  264. }
  269. /*void CMyDoc::TemplateMatch()
  270. {
  271. int i;
  272. int j;
  273. int m;
  274. int n;
  275.     //中间结果 
  276.     double dsst;
  277. double dss;
  278. double dst;
  279. //相似性测度
  280. double R;
  281. //最大相似性测度
  282. double maxR;
  283. //最大相似性出现的位置
  284.     int maxWidth;
  285. int maxHeight;
  288. dst=0;
  290. for (n=0;n<Height;n++)
  291. {
  292. for (m=0;m<Width;m++)
  293. {
  294. dst+=(double)DataTemp[n][m];
  295. }
  296. }
  297. maxR=0.0;
  298. for(j=0;j<m_y-Height+1;j++)
  299. {
  300. for(i=0;i<m_x-Width+1;i++)
  301. {
  302. dsst=0;
  303. dss=0;
  304. for(n=0;n<Height;n++)
  305. {
  306. for(m=0;m<Width;m++)
  307. {
  308. dss+=DataMap1[j+n][i+m];
  309. dsst+=DataTemp[n][m];
  310. }
  311. }
  312. //计算相似性
  313. R=dsst/(sqrt(dss)*sqrt(dst));
  314.             //与最大相似性比较
  315. if(R>maxR)
  316. {
  317. maxR=R;
  318. maxWidth=i;
  319. maxHeight=j;
  320. }
  321. }
  322. }
  323.     Match_x=maxWidth;
  324. Match_y=maxHeight;
  326.    
  327. }*/
  329. void CMyDoc::OnPinjie() 
  330. {
  331. // TODO: Add your command handler code here
  332.     //预处理？
  333. HLOCAL hNew;
  335.     int L;
  336. L=size.cx-Match_x+m3_x;
  338. //分配空间
  339.     hNew=LocalAlloc(LHND,size.cy*((2*size.cx-L+3)/4*4));
  340.     NewP=(BYTE*)LocalLock(hNew);
  341.     
  342. //进行拼接  
  343.     for(int j=0;j<size.cy;j++)
  344. for(int i=0;i<TrueWidth;i++)
  345. {
  346.         NewP[j*((2*size.cx-L+3)/4*4)+i]=Lp1[j*TrueWidth+i];
  347. }
  349.     for(int n=0;n<size2.cy;n++)
  350. for(int m=L;m<TrueWidth2;m++)
  351. {
  352. NewP[n*((2*size.cx-L+3)/4*4)+TrueWidth2+m-L]=Lp2[n*TrueWidth2+m];
  353. }
  356. m_bmpHead->bfSize=(DWORD)(sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFO)+
  357. 256*sizeof(RGBQUAD)+(DWORD)(size.cy*((2*size.cx-L+3)/4*4)));
  359. m_bmpInfo=(BITMAPINFO*)GlobalAllocPtr(GHND,m_bmpHead->bfSize-sizeof(BITMAPFILEHEADER));
  360. memcpy(m_bmpInfo,lbmpInfo,m_bmpHead->bfSize-sizeof(BITMAPFILEHEADER));
  361. m_bmpInfo->bmiHeader.biWidth=(2*size.cx-L+3)/4*4;
  362. m_bmpInfo->bmiHeader.biHeight=size.cy;
  363. m_bmpInfo->bmiHeader.biSizeImage=size.cy*((2*size.cx-L+3)/4*4);
  365.         //释放内存
  366.    
  367. UpdateAllViews(NULL);
  370. }
  372. void CMyDoc::ShowNew(CDC *pDC)
  373. {   
  374. ::StretchDIBits(pDC->m_hDC,0,500,m_bmpInfo->bmiHeader.biWidth,m_bmpInfo->bmiHeader.biHeight,0,0,
  375. m_bmpInfo->bmiHeader.biWidth,m_bmpInfo->bmiHeader.biHeight,NewP,m_bmpInfo,DIB_RGB_COLORS,SRCCOPY);
  376. }
  379. void CMyDoc::OnSuanfa1() 
  380. {
  381. // TODO: Add your command handler code here
  382.     int a1,a2,a3,a4;
  383.     CZuoBiao XY;
  384. XY.DoModal();
  385. a1=XY.mx;
  386. a2=XY.my;
  387. x1=XY.mx;
  388. y1=XY.my;
  389. a3=XY.h;
  390. a4=XY.w;
  392. // 更改光标形状
  393. BeginWaitCursor();
  395. //获取模板图像
  396.     Moban(a1,a2,a3,a4);
  397. int i;
  398. int j;
  399. int m;
  400. int n;
  402.     //中间结果 
  403.     double dsst;
  404. double dss;
  405. double dst;
  407. //相似性测度
  408. double R;
  410. //最大相似性测度
  411. double maxR;
  413. //最大相似性出现的位置
  414.     int maxWidth;
  415. int maxHeight;
  418. dst=0;
  420. for (n=0;n<Height;n++)
  421. {
  422. for (m=0;m<Width;m++)
  423. {
  424. dst+=(double)Lp3[n*Width+m]*Lp3[n*Width+m];
  425. }
  426. }
  428. maxR=0.0;
  430. for(j=0;j<size.cy-Height+1;j++)
  431. {
  432. for(i=0;i<TrueWidth-Width+1;i++)
  433. {
  434. dsst=0;
  435. dss=0;
  436. for(n=0;n<Height;n++)
  437. {
  438. for(m=0;m<Width;m++)
  439. {
  440. dss+=Lp1[(j+n)*TrueWidth+i+m]*Lp1[(j+n)*TrueWidth+i+m];
  441. dsst+=Lp3[n*Width+m]*Lp1[(j+n)*TrueWidth+i+m];
  442. }
  443. }
  445. //计算相似性
  446. R=dsst/(sqrt(dss)*sqrt(dst));
  448.             //与最大相似性比较
  449. if(R>maxR)
  450. {
  451. maxR=R;
  452. maxWidth=i;
  453. maxHeight=size.cy-Height-j;
  454. }
  455. }
  456. }
  458.     Match_x=maxWidth;
  459. Match_y=maxHeight;
  461. Weizhi d;
  462. d.match_x=Match_x;
  463. d.match_y=Match_y;
  464. d.DoModal();
  465. // 恢复光标形状
  466. EndWaitCursor();
  468. UpdateAllViews(NULL);
  469. }
  471. void CMyDoc::OnSuanfa2() 
  472. {
  473. // TODO: Add your command handler code here
  475. //定义旋转变量
  476.     //float srcX1,srcY1,srcX2,srcY2,srcX3,srcY3,srcX4,srcY4;
  477. //float dstX1,dstY1,dstX2,dstY2,dstX3,dstY3,dstX4,dstY4;
  478. float sina,cosa;
  479. //float num1,num2;
  480. float RotateAngle;
  481. //int Wnew,Hnew;
  482. int k;
  483. int x0,y0,x2,y2;
  484.     float x1,y1,x,y;
  485. double zj;
  486.    
  487.     int a1,a2,a3,a4;
  488.     CZuoBiao XY;
  489. XY.DoModal();
  490. a1=XY.mx;
  491. a2=XY.my;
  492. a3=XY.h;
  493. a4=XY.w;
  495. // 更改光标形状
  496. BeginWaitCursor();
  498. //获取模板图像
  499.     Moban(a1,a2,a3,a4);
  500. int i;
  501. int j;
  502. int n;
  503. int m;
  505.     //中间结果 
  506.     double dsst;
  507. double dss;
  508. double dst;
  510. //相似性测度
  511. double R;
  513. //最大相似性测度
  514. double maxR;
  515. double Max;
  517. //最大相似性出现的位置
  518.     int maxWidth;
  519. int maxHeight;
  520. int matchWidth;
  521. int matchHeight;
  524. dst=0;
  526. for (n=0;n<Height;n++)
  527. {
  528. for (m=0;m<Width;m++)
  529. {
  530. dst+=(double)Lp3[n*Width+m]*Lp3[n*Width+m];
  531. }
  532. }
  533. maxR=0.0;
  534. Max=0.0;
  536.     //以0。1度步进来旋转
  537. for (k=1;k<=30;k++)
  538. {
  539. RotateAngle=(float)RADIAN(k);
  540.         cosa=(float)cos((double)RotateAngle);
  541. sina=(float)sin((double)RotateAngle);
  545. for(j=50;j<size.cy-Height+1;j++)
  546. {
  547.      for(i=50;i<TrueWidth-Width+1;i++)
  548. {
  549.      dsst=0;
  550.      dss=0;
  551.       for(n=0;n<Height;n++)
  552. {
  553.      for(m=0;m<Width;m++)
  554. {
  555. x0=j+n;
  556. y0=i+m;
  557. //x1=x0*cosa-y0*sina+num1;
  558. //y1=x0*sina-y0*cosa+num2;
  559. x1=(x0-i)*cosa-(y0-j)*sina+i;
  560. y1=(x0-i)*sina+(y0-j)*cosa+j;
  561.                         x2=(int)(x1*10/10);
  562. y2=(int)(y1*18/10);
  563. x=x1-x2;
  564. y=y1-y2;
  565. /*zj=(Lp1[x2*TrueWidth+y2+1]-Lp1[x2*TrueWidth+y2])*y+
  566. (Lp1[(x2+1)*TrueWidth+y2]-Lp1[x2*TrueWidth+y2])*x+
  567. (Lp1[(x2+1)*TrueWidth+y2+1]+Lp1[x2*TrueWidth+y2]-
  568.                              Lp1[x2*TrueWidth+y2+1]-Lp1[(x2+1)*TrueWidth+y2])*x*y+
  569.                              Lp1[x2*TrueWidth+y2];*/
  570. zj=Lp1[x2*TrueWidth+y2];
  571. dss+=zj*zj;
  572. dsst+=Lp3[n*Width+m]*zj;
  573. }
  574. }
  576. //计算相似性
  577.      R=dsst/(sqrt(dss)*sqrt(dst));
  579.     //与最大相似性比较
  580.      if(R>maxR)
  581. {
  582.      maxR=R;
  583.      maxWidth=i;
  584.       maxHeight=size.cy-Height-j;
  585. }
  586. }
  587. }
  589. if(maxR>Max)
  590. {
  591. Max=maxR;
  592. matchWidth=maxWidth;
  593. matchHeight=maxHeight;
  594. matchAngle=(float)(k);
  595. }
  596. }
  598. //获得最佳匹配位置
  599.         Match_x=maxWidth;
  600.      Match_y=maxHeight;
  601.         MatchRotateAngle=matchAngle; 
  603.       Wuizhi2 e;
  604.      e.match2_x=Match_x;
  605.      e.match2_y=Match_y;
  606. e.angle=MatchRotateAngle;
  608.      e.DoModal();
  610. // 恢复光标形状
  611. EndWaitCursor();
  613.      UpdateAllViews(NULL);
  615. }
  618.                               
  621. void CMyDoc::OnRotatepinjie() 
  622. {
  623. // TODO: Add your command handler code here
  624.     float srcX1,srcY1,srcX2,srcY2,srcX3,srcY3,srcX4,srcY4;
  625. float dstX1,dstY1,dstX2,dstY2,dstX3,dstY3,dstX4,dstY4;
  626. float sina,cosa;
  627. float num1,num2;
  628. float Rotateangle;
  629. int Wnew,Hnew,TrueWnew;
  630. int i,j,i0,j0;
  631. int k,l;
  632. int L;
  633. int value;
  634. int a;
  635. int bb;
  636. // float a,b;
  638.     HLOCAL htemp;
  639. HLOCAL yuantu;
  640. HLOCAL hNew;
  641. LPBYTE lpPtr;
  642. LPBYTE lpsrc; 
  643.                   
  644. RotateAngle A;
  645. A.DoModal();
  646.     bb=A.Angle;
  647. Rotateangle=float(360-bb); //360-MatchRotateAngle;
  649.     switch (bb)
  650. {   case  1: value=5; a=0;break;
  651. case  2: value=7; a=0;break;
  652. case  4: value=15;a=0;break;
  653. case  3: value=12;a=0;break;
  654.         case  5: value=20;a=0;break;
  655. case  6: value=26;a=0;break;
  656. case  7: value=28;a=0;break;
  657. case  8: value=32;a=0;break;
  658. case  9: value=37;a=0;break;
  659. case 10: value=41;a=0;break;
  660. case 11: value=44;a=0;break;
  661. case 12: value=49;a=0;break;
  662. case 13: value=52;a=0;break;
  663. case 14: value=54;a=0;break;
  664. case 15: value=56;a=0;break;
  665. case 16: value=60;a=0;break;
  666. case 17: value=65;a=0;break;
  667. case 18: value=70;a=0;break;
  668. case 19: value=72;a=0;break;
  669. case 20: value=74;a=0;break;
  670. case 21: value=77;a=0;break;
  671. case 22: value=79;a=0;break;
  672. case 23: value=82;a=0;break;
  673. case 24: value=85;a=0;break;
  674. case 25: value=89;a=0;break;
  675. case 26: value=92;a=0;break;
  676. case 27: value=95;a=0;break;
  677. case 28: value=98;a=0;break;
  678. case 29: value=102;a=0;break;
  679. case 30: value=105;a=3;break;
  680. }
  681.     cosa=(float)cos((double)RADIAN(Rotateangle));
  682. sina=(float)sin((double)RADIAN(Rotateangle));
  683.   /*  cosa=(float)(((Match_x-Match3_x)*(x1-x2)+(Match_y-Match3_x)*(y1-y2))/((x1-x2)*(x1-x2)+
  684.      (y1-y2)*(y1-y2)));
  685. sina=(float)(((Match_x-Match3_x)*(y1-y2)-(x1-x2)*(Match_y-Match3_x))/((x1-x2)*(x1-x2)+
  686.          (y1-y2)*(y1-y2)));*/
  687. // a=acos(cosa);
  688. // b=asin(sina);//旋转2度是7，5度是20，30度是104，10度是41,20度是79，30度是104
  691. //原图的4个角点的坐标
  692.     srcX1=(float)(-(TrueWidth-1)/2);
  693.     srcY1=(float)( (size.cy-1)/2);
  694.     srcX2=(float)( (TrueWidth-1)/2);
  695.     srcY2=(float)( (size.cy-1)/2);
  696.     srcX3=(float)(-(TrueWidth-1)/2);
  697.     srcY3=(float)(-(size.cy-1)/2);
  698.     srcX4=(float)( (TrueWidth-1)/2);
  699.     srcY4=(float)(-(size.cy-1)/2);
  700.    
  702.     //新图4个角点的坐标
  703. dstX1= cosa*srcX1+sina*srcY1;
  704. dstY1=-sina*srcX1+cosa*srcY1;
  705.     dstX2= cosa*srcX2+sina*srcY2;
  706.     dstY2=-sina*srcX2+cosa*srcY2;
  707.     dstX3= cosa*srcX3+sina*srcY3;
  708.     dstY3=-sina*srcX3+cosa*srcY3;
  709. dstX4= cosa*srcX4+sina*srcY4;
  710. dstY4=-sina*srcX4+cosa*srcY4;
  712.     //新图的高和宽
  713.     Wnew=(int)(max(fabs(dstX4-dstX1),fabs(dstX3-dstX2))+0.5);
  714. Hnew=(int)(max(fabs(dstY4-dstY1),fabs(dstY3-dstY2))+0.5);
  716. //真实宽度
  717.     TrueWnew=(Wnew+3)/4*4;
  719.     num1=(float)(-0.5*(Wnew-1)*cosa-0.5*(Hnew-1)*sina+0.5*(TrueWidth-1));
  720.     num2=(float)( 0.5*(Wnew-1)*sina-0.5*(Hnew-1)*cosa+0.5*(size.cy-1));
  722.     //分配空间
  723.     htemp=LocalAlloc(LHND,Hnew*TrueWnew);
  724.     lpPtr=(BYTE*)LocalLock(htemp);
  725.     
  726. //内存中的象素都是白色
  727. //memcpy(lpPtr,BYTE(255),Hnew*TrueWnew);
  729. for(i=0;i<Hnew;i++)
  730. {
  731. for(j=0;j<Wnew;j++)
  732. {
  733.             //计算在原图中的坐标 
  734. i0=(int)(-((float)j)*sina+((float)i)*cosa+num2+0.5);
  735. j0=(int)( ((float)j)*cosa+((float)i)*sina+num1+0.5);
  736.             
  737.     //判断是否在原图范围内
  738. if((j0>=0)&&(j0<TrueWidth)&&(i0>=0)&&(i0<size.cy))
  739. {
  740. //复制图象
  741. lpPtr[TrueWnew*(Hnew-1-i)+j]=Lp1[TrueWidth*(size.cy-1-i0)+j0];
  742. }
  743. else
  744. {
  745. //对于图中没有的图象，直接赋为255
  746. lpPtr[TrueWnew*(Hnew-1-i)+j]=255;
  747. }
  748. }
  749. }
  751.       // 分配内存
  752.     yuantu=LocalAlloc(LHND,size2.cy*TrueWidth2);
  753.     lpsrc=(BYTE*)LocalLock(yuantu);
  755.         for (k=0;k<size2.cy;k++)
  756.      for (l=0;l<TrueWidth2;l++)
  757. {
  758.      lpsrc[k*TrueWidth2+l]=lpPtr[(Hnew-value-size2.cy+k)*TrueWnew+value+l];
  760. }
  762.        //旋转2度是7，5度是20，30度是104，10度是41,20度是79，30度是104
  764.      //L=size.cx-Match_x+m3_x;
  765.         L=128;
  766.      //分配空间
  767.         hNew=LocalAlloc(LHND,size2.cy*((2*size2.cx-L+3)/4*4));
  768.         NewP=(BYTE*)LocalLock(hNew);
  769.     
  770.     //进行拼接  
  771.         for(int p=0;p<size2.cy;p++)
  772.    for(int q=0;q<TrueWidth2;q++)
  773.    {
  774.           NewP[p*((2*size2.cx-L+3)/4*4)+q]=lpsrc[p*TrueWidth2+q];
  775.    }
  777.         for(int n=0;n<size2.cy;n++)
  778.      for(int m=L;m<TrueWidth2;m++)
  779. {
  780.     NewP[n*((2*size2.cx-L+3)/4*4)+TrueWidth2+m-L]=Lp2[(n+a)*TrueWidth2+m];
  781. }
  783. //开辟空间
  784. m_bmpHead=(BITMAPFILEHEADER*)GlobalAllocPtr(GHND,sizeof(BITMAPFILEHEADER));
  785. memcpy(m_bmpHead,&lbmpHeader,sizeof(BITMAPFILEHEADER));
  787. m_bmpHead->bfOffBits=(DWORD)(sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFO)+256*sizeof(RGBQUAD));
  788. m_bmpHead->bfSize=(DWORD)(sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFO)+
  789. 256*sizeof(RGBQUAD)+(DWORD)(size2.cy*((2*size2.cx-L+3)/4*4)));
  791. m_bmpInfo=(BITMAPINFO*)GlobalAllocPtr(GHND,m_bmpHead->bfSize-sizeof(BITMAPFILEHEADER));
  792.       memcpy(m_bmpInfo,lbmpInfo,m_bmpHead->bfSize-sizeof(BITMAPFILEHEADER));
  793. m_bmpInfo->bmiHeader.biWidth=(2*size2.cx-L+3)/4*4;
  794. m_bmpInfo->bmiHeader.biHeight=size2.cy;
  795. m_bmpInfo->bmiHeader.biSizeImage=size2.cy*((2*size2.cx-L+3)/4*4);
  796.                                              
  797. // lpBits=(LPBYTE)lbmpInfo+sizeof(BITMAPINFOHEADER)+256*sizeof(RGBQUAD);
  798. // memcpy(lpBits,NewP,lbmpInfo->bmiHeader.biSizeImage);
  800.      //   size.cx=lbmpInfo->bmiHeader.biWidth;
  801. // size.cy=lbmpInfo->bmiHeader.biHeight;
  803.         //释放内存
  804.    
  805. UpdateAllViews(NULL);
  807. }
  812.      
  817. void CMyDoc::OnSuanfa3() 
  818. {
  819. // TODO: Add your command handler code here
  820.     int a1,a2,a3,a4;
  821.     CZuoBiao XY;
  822. XY.DoModal();
  823. a1=XY.mx;
  824. a2=XY.my;
  825. x2=XY.mx;
  826. y2=XY.my;
  827. a3=XY.h;
  828. a4=XY.w;
  830. // 更改光标形状
  831. BeginWaitCursor();
  833. //获取模板图像
  834.     Moban(a1,a2,a3,a4);
  835. int i;
  836. int j;
  837. int m;
  838. int n;
  840.     //中间结果 
  841.     double dsst;
  842. double dss;
  843. double dst;
  845. //相似性测度
  846. double R;
  848. //最大相似性测度
  849. double maxR;
  851. //最大相似性出现的位置
  852.     int maxWidth;
  853. int maxHeight;
  856. dst=0;
  858. for (n=0;n<Height;n++)
  859. {
  860. for (m=0;m<Width;m++)
  861. {
  862. dst+=(double)Lp3[n*Width+m]*Lp3[n*Width+m];
  863. }
  864. }
  866. maxR=0.0;
  868. for(j=0;j<size.cy-Height+1;j++)
  869. {
  870. for(i=0;i<TrueWidth-Width+1;i++)
  871. {
  872. dsst=0;
  873. dss=0;
  874. for(n=0;n<Height;n++)
  875. {
  876. for(m=0;m<Width;m++)
  877. {
  878. dss+=Lp1[(j+n)*TrueWidth+i+m]*Lp1[(j+n)*TrueWidth+i+m];
  879. dsst+=Lp3[n*Width+m]*Lp1[(j+n)*TrueWidth+i+m];
  880. }
  881. }
  883. //计算相似性
  884. R=dsst/(sqrt(dss)*sqrt(dst));
  886.             //与最大相似性比较
  887. if(R>maxR)
  888. {
  889. maxR=R;
  890. maxWidth=i;
  891. maxHeight=size.cy-Height-j;
  892. }
  893. }
  894. }
  896.     Match3_x=maxWidth;
  897. Match3_y=maxHeight;
  899. Weizhi3 f;
  900. f.match3_x=Match3_x;
  901. f.match3_y=Match3_y;
  902. f.DoModal();
  904. // 恢复光标形状
  905. EndWaitCursor();
  907. UpdateAllViews(NULL);
  909. }
  911. void CMyDoc::OnChuizhipinjie() 
  912. {
  913. // TODO: Add your command handler code here
  914. //预处理？
  915. HLOCAL hNew;
  917.     int L;
  918. L=size.cy-Match_y+m3_y;
  920. //分配空间
  921.     hNew=LocalAlloc(LHND,TrueWidth*(2*size.cy-L));
  922.     NewP=(BYTE*)LocalLock(hNew);
  923.     
  924. //进行拼接  
  925.     for(int j=0;j<size.cy;j++)
  926. for(int i=0;i<TrueWidth;i++)
  927. {
  928.         NewP[j*TrueWidth+i]=Lp2[j*TrueWidth+i];
  929. }
  931.     for(int n=L;n<size.cy;n++)
  932. for(int m=0;m<TrueWidth;m++)
  933. {
  934. NewP[(n+size.cy-L)*TrueWidth+m]=Lp1[n*TrueWidth+m];
  935. }
  938. m_bmpHead->bfSize=(DWORD)(sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFO)+
  939. 256*sizeof(RGBQUAD)+(DWORD)(TrueWidth*(2*size.cy-L)));
  941. m_bmpInfo=(BITMAPINFO*)GlobalAllocPtr(GHND,m_bmpHead->bfSize-sizeof(BITMAPFILEHEADER));
  942. memcpy(m_bmpInfo,lbmpInfo,m_bmpHead->bfSize-sizeof(BITMAPFILEHEADER));
  943. m_bmpInfo->bmiHeader.biWidth=TrueWidth;
  944. m_bmpInfo->bmiHeader.biHeight=2*size.cy-L;
  945. m_bmpInfo->bmiHeader.biSizeImage=TrueWidth*(2*size.cy-L);
  947.         //释放内存
  948.    
  949. UpdateAllViews(NULL);
  952. }
  954. /*void CMyDoc::DeleteContents() 
  955. {
  956. // TODO: Add your specialized code here and/or call the base class
  957. int liCount=m_oalines.GetSize;
  958. int liPos;
  960. if(liCount)
  961. {
  962. for (liPos=0;liPos<liCount;liPos++)
  963. delete m_oaLines[liPos];
  964. m_oaLines.RemoveAll();
  965. }
  966.         
  968. CDocument::DeleteContents();
  969. }*/
  971. //DEL BOOL CMyDoc::OnSaveDocument(LPCTSTR lpszPathName) 
  972. //DEL {
  973. //DEL  // TODO: Add your specialized code here and/or call the base class
  974. //DEL   CFile fp;
  975. //DEL  fp.Open(lpszPathName,CFile::modeWrite|CFile::modeCreate);
  976. //DEL  fp.Write((LPSTR)&lbmpHeader,sizeof(BITMAPFILEHEADER));
  977. //DEL  fp.Write(lpBits,lbmpHeader.bfSize-sizeof(BITMAPFILEHEADER));
  978. //DEL  fp.Close();
  979. //DEL 
  980. //DEL 
  981. //DEL  SetModifiedFlag(FALSE);
  982. //DEL  return TRUE;
  983. //DEL 
  984. //DEL  //return CDocument::OnSaveDocument(lpszPathName);
  985. //DEL }
  989. void CMyDoc::SaveBitmap()
  990. { 
  991.   CFile fp;
  992.   CFileDialog dlg(NULL,"BMP","*.bmp");
  993.   if(dlg.DoModal()==IDOK)
  994.   {
  995.       fp.Open(dlg.GetPathName(),CFile::modeWrite|CFile::modeCreate);
  996.     fp.Write(m_bmpHead,sizeof(BITMAPFILEHEADER));
  997.       fp.Write(m_bmpInfo,sizeof(BITMAPINFO));
  998.   fp.Write(rgbPal,256*sizeof(RGBQUAD));
  999.   fp.Write(NewP,m_bmpHead->bfSize-sizeof(BITMAPFILEHEADER));
  1000.   }
  1001.   else
  1002.   return;   
  1003.   fp.Close();
  1004. }
  1006. void CMyDoc::OnSavebitmap() 
  1007. {
  1008. // TODO: Add your command handler code here
  1009. SaveBitmap();
  1010. }
  1012. BOOL CMyDoc::OnSaveDocument(LPCTSTR lpszPathName) 
  1013. {
  1014. // TODO: Add your specialized code here and/or call the base class
  1016. return CDocument::OnSaveDocument(lpszPathName);
  1017. }