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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > 生物技术 > 查看源码
FaceDetectDlg.cpp
资源名称：face.rar　[点击查看]
上传用户：joy_gong
上传日期：2020-11-30
资源大小：459k
文件大小：44k
源码类别：

生物技术

开发平台：

Visual C++

FaceDetectDlg.cpp：源码内容
  1. // FaceDetectDlg.cpp : implementation file
  2. //
  4. #include "stdafx.h"
  5. #include "FaceDetect.h"
  6. #include "FaceDetectDlg.h"
  7. #include "ReplaceDlg.h"
  9. #ifdef _DEBUG
  10. #define new DEBUG_NEW
  11. #undef THIS_FILE
  12. static char THIS_FILE[] = __FILE__;
  13. #endif
  15. #include "AddSampleDlg.h"
  17. /////////////////////////////////////////////////////////////////////////////
  18. // CAboutDlg dialog used for App About
  20. class CAboutDlg : public CDialog
  21. {
  22. public:
  23. CAboutDlg();
  25. // Dialog Data
  26. //{{AFX_DATA(CAboutDlg)
  27. enum { IDD = IDD_ABOUTBOX };
  28. //}}AFX_DATA
  30. // ClassWizard generated virtual function overrides
  31. //{{AFX_VIRTUAL(CAboutDlg)
  32. protected:
  33. virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
  34. //}}AFX_VIRTUAL
  36. // Implementation
  37. protected:
  38. //{{AFX_MSG(CAboutDlg)
  39. //}}AFX_MSG
  40. DECLARE_MESSAGE_MAP()
  41. };
  43. CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
  44. {
  45. //{{AFX_DATA_INIT(CAboutDlg)
  46. //}}AFX_DATA_INIT
  47. }
  49. void CAboutDlg::DoDataExchange(CDataExchange* pDX)
  50. {
  51. CDialog::DoDataExchange(pDX);
  52. //{{AFX_DATA_MAP(CAboutDlg)
  53. //}}AFX_DATA_MAP
  54. }
  56. BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
  57. //{{AFX_MSG_MAP(CAboutDlg)
  58. // No message handlers
  59. //}}AFX_MSG_MAP
  60. END_MESSAGE_MAP()
  62. /////////////////////////////////////////////////////////////////////////////
  63. // CFaceDetectDlg dialog
  65. CFaceDetectDlg::CFaceDetectDlg(CWnd* pParent /*=NULL*/)
  66. : CDialog(CFaceDetectDlg::IDD, pParent)
  67. {
  68. //{{AFX_DATA_INIT(CFaceDetectDlg)
  69. // NOTE: the ClassWizard will add member initialization here
  70. //}}AFX_DATA_INIT
  71. // Note that LoadIcon does not require a subsequent DestroyIcon in Win32
  72. m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
  73. }
  75. void CFaceDetectDlg::DoDataExchange(CDataExchange* pDX)
  76. {
  77. CDialog::DoDataExchange(pDX);
  78. //{{AFX_DATA_MAP(CFaceDetectDlg)
  79. // NOTE: the ClassWizard will add DDX and DDV calls here
  80. //}}AFX_DATA_MAP
  81. }
  83. BEGIN_MESSAGE_MAP(CFaceDetectDlg, CDialog)
  84. //{{AFX_MSG_MAP(CFaceDetectDlg)
  85. ON_WM_SYSCOMMAND()
  86. ON_WM_PAINT()
  87. ON_WM_QUERYDRAGICON()
  88. ON_BN_CLICKED(IDC_BTN_BINARY, OnBtnBinary)
  89. ON_BN_CLICKED(IDC_BTN_EDGE, OnBtnEdge)
  90. ON_BN_CLICKED(IDC_BTN_FACEHAIR, OnBtnFacehair)
  91. ON_BN_CLICKED(IDC_BTN_HISTOGRAM_FACE, OnBtnHistogramFace)
  92. ON_BN_CLICKED(IDC_BTN_HISTOGRAM_H, OnBtnHistogramH)
  93. ON_BN_CLICKED(IDC_BTN_HISTOGRAM_HAIR, OnBtnHistogramHair)
  94. ON_BN_CLICKED(IDC_BTN_HISTOGRAM_V, OnBtnHistogramV)
  95. ON_BN_CLICKED(IDC_BTN_LIKEHOOD, OnBtnLikehood)
  96. ON_BN_CLICKED(IDC_BTN_MARK_EYE, OnBtnMarkEye)
  97. ON_BN_CLICKED(IDC_BTN_MARK_FACE_1, OnBtnMarkFace1)
  98. ON_BN_CLICKED(IDC_BTN_MARK_FACE_2, OnBtnMarkFace2)
  99. ON_BN_CLICKED(IDC_BTN_MARK_MOUSE, OnBtnMarkMouse)
  100. ON_BN_CLICKED(IDC_BTN_MARK_NOSE, OnBtnMarkNose)
  101. ON_BN_CLICKED(IDC_BTN_OPENFILE, OnBtnOpenfile)
  102. //}}AFX_MSG_MAP
  103. END_MESSAGE_MAP()
  105. /////////////////////////////////////////////////////////////////////////////
  106. // CFaceDetectDlg message handlers
  108. BOOL CFaceDetectDlg::OnInitDialog()
  109. {
  110. CDialog::OnInitDialog();
  112. // Add "About..." menu item to system menu.
  114. // IDM_ABOUTBOX must be in the system command range.
  115. ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
  116. ASSERT(IDM_ABOUTBOX < 0xF000);
  118. CMenu* pSysMenu = GetSystemMenu(FALSE);
  119. if (pSysMenu != NULL)
  120. {
  121. CString strAboutMenu;
  122. strAboutMenu.LoadString(IDS_ABOUTBOX);
  123. if (!strAboutMenu.IsEmpty())
  124. {
  125. pSysMenu->AppendMenu(MF_SEPARATOR);
  126. pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
  127. }
  128. }
  130. // Set the icon for this dialog.  The framework does this automatically
  131. //  when the application's main window is not a dialog
  132. SetIcon(m_hIcon, TRUE); // Set big icon
  133. SetIcon(m_hIcon, FALSE); // Set small icon
  135. // TODO: Add extra initialization here
  136. CWnd *pWnd0= GetDlgItem(IDC_BMPSHOW);
  137. pDCShow = pWnd0->GetDC();
  139. m_pMainDib = new CDib();
  141. m_tOriPixelArray = NULL;
  142. m_tResPixelArray = NULL;
  144. m_pResMap = NULL;
  146. m_nWndWidth = 0;
  147. m_nWndHeight= 0;
  148. m_sFileName = "";
  150. m_bSelectByMan = false;
  151. m_bLBottonDown = false;
  153. m_ManLeft = -1;
  154. m_ManRight = -1;
  155. m_ManTop = -1;
  156. m_ManBottom = -1;
  158. m_bFaceOK = false;
  159. m_bShowFace = false;
  160. m_rFaceRegion.left = m_rFaceRegion.right = m_rFaceRegion.top = m_rFaceRegion.bottom = 0;
  162. m_bManualMarkFacial = false;
  163. m_bLeftEyeOK = m_bRightEyeOK = m_bLeftNostrilOK = m_bRightNostrilOK =
  164. m_bLeftEyeLeftCornerOK = m_bLeftEyeRightCornerOK = m_bRightEyeLeftCornerOK = 
  165. m_bRightEyeRightCornerOK = m_bLeftMouthCornerOK = m_bRightMouthCornerOK = false;
  167. m_bMidMouthOK = m_bMidNoseOK = false;
  169. m_LeftEye = m_RightEye = m_LeftEyeLeftCorner = m_LeftEyeRightCorner = 
  170. m_LeftNostril = m_RightNostril = m_RightEyeLeftCorner = m_RightEyeRightCorner =
  171. m_LeftMouthCorner = m_RightMouthCorner = m_MidMouth = m_MidNose = CPoint(-1,-1);
  173. return TRUE;  // return TRUE  unless you set the focus to a control
  174. }
  176. void CFaceDetectDlg::OnSysCommand(UINT nID, LPARAM lParam)
  177. {
  178. if ((nID & 0xFFF0) == IDM_ABOUTBOX)
  179. {
  180. CAboutDlg dlgAbout;
  181. dlgAbout.DoModal();
  182. }
  183. else
  184. {
  185. CDialog::OnSysCommand(nID, lParam);
  186. }
  187. }
  189. // If you add a minimize button to your dialog, you will need the code below
  190. //  to draw the icon.  For MFC applications using the document/view model,
  191. //  this is automatically done for you by the framework.
  193. void CFaceDetectDlg::OnPaint() 
  194. {
  195. if (IsIconic())
  196. {
  197. CPaintDC dc(this); // device context for painting
  199. SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);
  201. // Center icon in client rectangle
  202. int cxIcon = GetSystemMetrics(SM_CXICON);
  203. int cyIcon = GetSystemMetrics(SM_CYICON);
  204. CRect rect;
  205. GetClientRect(&rect);
  206. int x = (rect.Width() - cxIcon + 1) / 2;
  207. int y = (rect.Height() - cyIcon + 1) / 2;
  209. // Draw the icon
  210. dc.DrawIcon(x, y, m_hIcon);
  212. if(m_tResPixelArray==NULL) return;
  213. }
  214. else
  215. {
  216. CDialog::OnPaint();
  217. }
  218. }
  220. // The system calls this to obtain the cursor to display while the user drags
  221. //  the minimized window.
  222. HCURSOR CFaceDetectDlg::OnQueryDragIcon()
  223. {
  224. return (HCURSOR) m_hIcon;
  225. }
  227. ////////////////////////////////////////////////////////////////////////////////
  228. // 画十字形标记
  229. // 参数：  pDC－CDC指针
  230. //         point－要画的点的坐标
  231. //         crColor－标记得颜色
  232. ////////////////////////////////////////////////////////////////////////////////
  233. void CFaceDetectDlg::DrawCross(CDC *pDC, CPoint point, COLORREF crColor)
  234. {
  235. CPen pen,*oldPen;
  236. pen.CreatePen(PS_SOLID,1,crColor);
  237. oldPen = (CPen*)pDC->SelectObject(&pen);
  238. pDC->MoveTo(point.x-7,point.y);
  239. pDC->LineTo(point.x+7,point.y);
  240. pDC->MoveTo(point.x,point.y-7);
  241. pDC->LineTo(point.x,point.y+7);
  242. pDC->SelectObject(oldPen);
  243. pen.DeleteObject();
  244. }
  246. ////////////////////////////////////////////////////////////////////////////////
  247. // 拷贝位图
  248. // 参数：  dest－目标位图指针
  249. //         source－源位图指针
  250. ////////////////////////////////////////////////////////////////////////////////
  251. bool CFaceDetectDlg::CopyBitMap(RGBQUAD **dest, RGBQUAD **source)
  252. {
  253. if(source==NULL || dest==NULL) 
  254. return false;
  255. for(int i=0; i<m_nWndHeight; i++)
  256. for(int j=0; j<m_nWndWidth; j++)
  257. dest[i][j]=source[i][j];
  259. return true;
  260. }
  262. ////////////////////////////////////////////////////////////////////////////////
  263. // 生成新的位图
  264. ////////////////////////////////////////////////////////////////////////////////
  265. void CFaceDetectDlg::MakeBitMap()
  266. {
  267. CClientDC ClientDC(pDCShow->GetWindow());
  268. if(m_pResMap!=NULL) delete m_pResMap;
  270. m_pResMap=new CBitmap();
  271. m_pResMap->CreateCompatibleBitmap(&ClientDC,m_nWndWidth,m_nWndHeight);
  273. CDC  dc;
  274. dc.CreateCompatibleDC(&ClientDC);
  275. dc.SelectObject(m_pResMap);
  277. for(int i=0; i<m_nWndHeight; i++)
  278. for(int j=0; j<m_nWndWidth; j++)
  279. dc.SetPixelV(j,i,RGB(m_tResPixelArray[i][j].rgbRed,m_tResPixelArray[i][j].rgbGreen,m_tResPixelArray[i][j].rgbBlue));
  281. if(m_bFaceOK && m_bShowFace)
  282. {
  283. CBrush Pen;
  284. Pen.CreateSolidBrush(RGB(255,0,0));
  285. dc.FrameRect(m_rFaceRegion,&Pen);
  286.   Pen.DeleteObject();
  287. }
  289. if(m_bLeftEyeOK) DrawCross(&dc,m_LeftEye,RGB(255,0,0));
  290. if(m_bRightEyeOK) DrawCross(&dc,m_RightEye,RGB(255,0,0));
  291. if(m_bLeftEyeLeftCornerOK) DrawCross(&dc,m_LeftEyeLeftCorner,RGB(255,0,255));
  292. if(m_bLeftEyeRightCornerOK) DrawCross(&dc,m_LeftEyeRightCorner,RGB(255,255,0));
  293. if(m_bRightEyeLeftCornerOK) DrawCross(&dc,m_RightEyeLeftCorner,RGB(255,0,255));
  294. if(m_bRightEyeRightCornerOK) DrawCross(&dc,m_RightEyeRightCorner,RGB(255,255,0));
  295. if(m_bLeftNostrilOK) DrawCross(&dc,m_LeftNostril,RGB(0,255,0));
  296. if(m_bRightNostrilOK) DrawCross(&dc,m_RightNostril,RGB(0,255,0));
  297. if(m_bMidNoseOK) DrawCross(&dc,m_MidNose,RGB(0,255,0));
  298. if(m_bLeftMouthCornerOK) DrawCross(&dc,m_LeftMouthCorner,RGB(0,0,255));
  299. if(m_bRightMouthCornerOK) DrawCross(&dc,m_RightMouthCorner,RGB(0,0,255));
  300. if(m_bMidMouthOK) DrawCross(&dc,m_MidMouth,RGB(0,0,255));
  302. dc.DeleteDC();
  304. MyDraw();
  305. }
  307. ////////////////////////////////////////////////////////////////////////////////////
  308. //读原图的数据
  309. ////////////////////////////////////////////////////////////////////////////////////
  310. void CFaceDetectDlg::LoadOriPixel(CDib *pDib)
  311. {
  312. BYTE *colorTable;
  313. colorTable = (BYTE *)pDib->m_pDibBits;
  314. int byteBitCount  = pDib->GetBiBitCount()/8;
  316. m_tOriPixelArray  = new RGBQUAD*[m_nWndHeight];
  317. m_tResPixelArray  = new RGBQUAD*[m_nWndHeight];
  318. for(int l=0 ; l<m_nWndHeight; l++)
  319. {
  320. m_tOriPixelArray[l] = new RGBQUAD[m_nWndWidth];
  321. m_tResPixelArray[l] = new RGBQUAD[m_nWndWidth];
  322. }
  324. int count = 0;
  325. for(int i=m_nWndHeight-1; i>=0; i--)
  326. {
  327. for(int j=0; j<m_nWndWidth; j++)
  328. {
  329. m_tOriPixelArray[i][j].rgbBlue =colorTable[count++];
  330. m_tOriPixelArray[i][j].rgbGreen=colorTable[count++];
  331. m_tOriPixelArray[i][j].rgbRed  =colorTable[count++];
  332. m_tOriPixelArray[i][j].rgbReserved = 0;
  333. m_tResPixelArray[i][j]=m_tOriPixelArray[i][j];
  334. count += byteBitCount-3;
  335. }
  336. count += (4-(m_nWndWidth*byteBitCount)%4)%4;
  337. }
  339. method1 = new CLikelyHood(m_tOriPixelArray,m_nWndWidth,m_nWndHeight);
  340. method2 = new CHairFace(m_tOriPixelArray,m_nWndWidth,m_nWndHeight);
  341. }
  343. ////////////////////////////////////////////////////////////////////////////////
  344. // 给位图赋值
  345. // 参数：  target－目标位图指针
  346. //         Val－要赋予的值
  347. ////////////////////////////////////////////////////////////////////////////////
  348. void CFaceDetectDlg::SetPixelArray(RGBQUAD **target, int Val)
  349. {
  350. for(int i=0; i<m_nWndHeight; i++)
  351. for(int j=0; j<m_nWndWidth; j++)
  352. {
  353. target[i][j].rgbRed =  target[i][j].rgbBlue = target[i][j].rgbGreen = Val;
  354. }
  355. }
  357. ////////////////////////////////////////////////////////////////////////////////
  358. // 边界检测
  359. ////////////////////////////////////////////////////////////////////////////////
  360. void CFaceDetectDlg::DoLOG(int left, int right, int top, int bottom, RGBQUAD **source, RGBQUAD **target)
  361. {
  362. int i,j;
  364. double **result;  
  365. result = new double*[m_nWndHeight];
  366. for(int l=0 ; l<m_nWndHeight; l++)
  367. {
  368. result[l] = new double[m_nWndWidth];
  369. for(j=0; j<m_nWndWidth; j++)
  370. result[l][j] = source[l][j].rgbRed;
  371. }
  373. for(i=0; i<m_nWndHeight; i++)
  374. for(j=0; j<m_nWndWidth; j++)
  375. {
  376. double r,g,temp;
  377. temp = source[i][j].rgbGreen+source[i][j].rgbRed+source[i][j].rgbBlue;
  378. r = (double)source[i][j].rgbRed/temp;
  379. g = (double)source[i][j].rgbGreen/temp;
  380. if(g<0.398 && g > 0.246 && r<0.664 && r>0.233 && r>g && g>=0.5*(1-r))
  381. {
  382. target[i][j].rgbRed = 255;  //face
  383. }
  384. else target[i][j].rgbRed = 0;
  385. }
  387. for(i=top+2; i<bottom-2; i++)
  388. for(j=left+2; j<right-2; j++)
  389. {
  390. result[i][j] =  
  391.      (0-2.0/24.0)*((unsigned char)target[i-2][j-2].rgbRed) + 
  392.                  (0-4.0/24.0)*((unsigned char)target[i-2][j-1].rgbRed) + 
  393.      (0-4.0/24.0)*((unsigned char)target[i-2][j].rgbRed)   +
  394.      (0-4.0/24.0)*((unsigned char)target[i-2][j+1].rgbRed) +
  395.      (0-2.0/24.0)*((unsigned char)target[i-2][j+2].rgbRed) +
  396.  (0-4.0/24.0)*((unsigned char)target[i-1][j-2].rgbRed) + 
  397.  (8.0/24.0)  *((unsigned char)target[i-1][j].rgbRed)   +
  398.  (0-4.0/24.0)*((unsigned char)target[i-1][j+2].rgbRed) +
  399.  (0-4.0/24.0)*((unsigned char)target[i][j-2].rgbRed)   + 
  400.  (8.0/24.0)  *((unsigned char)target[i][j-1].rgbRed)   + 
  401.  (1.0)       *((unsigned char)target[i][j].rgbRed)     +
  402.  (8.0/24.0)  *((unsigned char)target[i][j+1].rgbRed)   +
  403.  (0-4.0/24.0)*((unsigned char)target[i][j+2].rgbRed)   +
  404.  (0-4.0/24.0)*((unsigned char)target[i+1][j-2].rgbRed) + 
  405.  (8.0/24.0)  *((unsigned char)target[i+1][j].rgbRed)   +
  406.  (0-4.0/24.0)*((unsigned char)target[i+1][j+2].rgbRed) +
  407.  (0-2.0/24.0)*((unsigned char)target[i+2][j-2].rgbRed) + 
  408.  (0-4.0/24.0)*((unsigned char)target[i+2][j-1].rgbRed) + 
  409.  (0-4.0/24.0)*((unsigned char)target[i+2][j].rgbRed)   +
  410.  (0-4.0/24.0)*((unsigned char)target[i+2][j+1].rgbRed) +
  411.  (0-2.0/24.0)*((unsigned char)target[i+2][j+2].rgbRed);
  412. }
  414. SetPixelArray(target,255);
  416. for(i=top+1; i<bottom-1; i++)
  417. for(j=left+1; j<right-1; j++)
  418. {
  419. int positive = 0;   
  420. int negtive  = 0;
  421. for(int m=-1;m<=1;m++)
  422. for(int n=-1;n<=1;n++)
  423. if(m!=0 || n!=0)
  424. {
  425. if(result[i+m][j+n]<-5)negtive++;
  426. if(result[i+m][j+n]>=5)positive++;
  427. }
  428. if(positive>2 && negtive>2) 
  429. {
  430. target[i][j].rgbBlue = target[i][j].rgbGreen = target[i][j].rgbRed = 0;
  431. }
  433. }
  435. if(result!=NULL)
  436. {
  437. for (int i=0 ;i<m_nWndHeight;i++)
  438. if(result[i]!=NULL) delete result[i];
  439. delete result;
  440. }
  441. }
  443. ////////////////////////////////////////////////////////////////////////////////
  444. // 二值化
  445. ////////////////////////////////////////////////////////////////////////////////
  446. void CFaceDetectDlg::OnBtnBinary() 
  447. {
  448. SetCursor(LoadCursor(NULL,IDC_WAIT));
  450. if(!method1->CalBinary())
  451. {
  452. AfxMessageBox("请先计算相似度！");
  453. SetCursor(LoadCursor(NULL,IDC_ARROW));
  454. return;
  455. }
  457. m_bShowFace = false;
  458. for(int i=0; i<m_nWndHeight; i++)
  459. for(int j=0; j<m_nWndWidth;  j++)
  460. {
  461. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  462. m_tResPixelArray[i][j].rgbRed  = (int)(method1->m_pBinaryArray[i][j]*255);
  463. }
  464. MakeBitMap();
  465. SetCursor(LoadCursor(NULL,IDC_ARROW));
  466. MyDraw();
  467. }
  469. ////////////////////////////////////////////////////////////////////////////////
  470. // 边缘提取
  471. ////////////////////////////////////////////////////////////////////////////////
  472. void CFaceDetectDlg::OnBtnEdge() 
  473. {
  474. if(!m_bFaceOK)
  475. {
  476. AfxMessageBox("请先确定脸部区域");
  477. return;
  478. }
  479. //左右眼的水平区域
  480. int nLeft,nRight,nTop,nBottom;
  482. nLeft = m_rFaceRegion.left-5  > 0 ? m_rFaceRegion.left-5:0;
  483. nRight = m_rFaceRegion.right+5 < m_nWndWidth? m_rFaceRegion.right+5:m_nWndWidth-1;
  484. nTop = m_rFaceRegion.top-5   > 0 ? m_rFaceRegion.top-5:0;
  485. nBottom = m_rFaceRegion.bottom+5< m_nWndHeight?m_rFaceRegion.bottom+5:m_nWndHeight-1;
  486. //边缘检查
  487. DoLOG(nLeft,nRight,nTop,nBottom,m_tOriPixelArray,m_tResPixelArray);
  488. MakeBitMap();
  489. }
  491. ////////////////////////////////////////////////////////////////////////////////
  492. // 求取头发和脸部区域
  493. ////////////////////////////////////////////////////////////////////////////////
  494. void CFaceDetectDlg::OnBtnFacehair() 
  495. {
  496. m_bShowFace = false;
  497. SetCursor(LoadCursor(NULL,IDC_WAIT));
  498. method2->MarkHairFace();
  499. for(int i=0; i<m_nWndHeight; i++)
  500. for(int j=0; j<m_nWndWidth;  j++)
  501. {
  502. switch(method2->m_pBinaryArray[i][j])
  503. {
  504. case 0:
  505. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen = 0;
  506. m_tResPixelArray[i][j].rgbRed  = 255;
  507. break;
  508. case 1:
  509. m_tResPixelArray[i][j].rgbBlue = 255;
  510. m_tResPixelArray[i][j].rgbGreen=m_tResPixelArray[i][j].rgbRed=0;
  511. break;
  512. case 2:
  513. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  514. m_tResPixelArray[i][j].rgbRed  = 0;
  515. break;
  516. }
  517. }
  518. MakeBitMap();
  519. SetCursor(LoadCursor(NULL,IDC_ARROW));
  520. }
  522. ////////////////////////////////////////////////////////////////////////////////
  523. // 脸部区域直方图
  524. ////////////////////////////////////////////////////////////////////////////////
  525. void CFaceDetectDlg::OnBtnHistogramFace() 
  526. {
  527. if(!method2->m_bBinaryOK)
  528. {
  529. AfxMessageBox("请先计算二值化图！");
  530. return;
  531. }
  533. m_bShowFace = false;
  534. SetCursor(LoadCursor(NULL,IDC_WAIT));
  535. for(int j=0; j<m_nWndWidth;  j++)
  536. {
  537. int count = 0;
  538. for(int i=0; i<m_nWndHeight; i++)
  539. {
  540. if(method2->m_pBinaryArray[i][j] == 0) count++;
  541. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  542. m_tResPixelArray[i][j].rgbRed  = 255;
  543. }
  544. for(i=m_nWndHeight-1; i>=m_nWndHeight-count;i--)
  545. {
  546. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  547. m_tResPixelArray[i][j].rgbRed  = 0;
  548. }
  549. }
  551. MakeBitMap();
  552. SetCursor(LoadCursor(NULL,IDC_ARROW));
  553. }
  555. ////////////////////////////////////////////////////////////////////////////////
  556. // 水平方向直方图
  557. ////////////////////////////////////////////////////////////////////////////////
  558. void CFaceDetectDlg::OnBtnHistogramH() 
  559. {
  560. if(!method1->m_bBinaryReady)
  561. {
  562. AfxMessageBox("请先计算二值图");
  563. return;
  564. }
  565. m_bShowFace = false;
  566. SetCursor(LoadCursor(NULL,IDC_WAIT));
  567. for(int j=0; j<m_nWndWidth;  j++)
  568. {
  569. int count = 0;
  570. for(int i=0; i<m_nWndHeight; i++)
  571. {
  572. if(method1->m_pBinaryArray[i][j] == 1) count++;
  573. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  574. m_tResPixelArray[i][j].rgbRed  = 255;
  575. }
  576. for(i=m_nWndHeight-1; i>=m_nWndHeight-count;i--)
  577. {
  578. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  579. m_tResPixelArray[i][j].rgbRed  = 0;
  580. }
  581. }
  582. MakeBitMap();
  583. SetCursor(LoadCursor(NULL,IDC_ARROW));
  584. }
  586. ////////////////////////////////////////////////////////////////////////////////
  587. // 头发的直方图
  588. ////////////////////////////////////////////////////////////////////////////////
  589. void CFaceDetectDlg::OnBtnHistogramHair() 
  590. {
  591. if(!method2->m_bBinaryOK)
  592. {
  593. AfxMessageBox("请先计算二值图！");
  594. return;
  595. }
  596. m_bShowFace = false;
  597. SetCursor(LoadCursor(NULL,IDC_WAIT));
  598. for(int j=0; j<m_nWndWidth;  j++)
  599. {
  600. int count = 0;
  601. for(int i=0; i<m_nWndHeight; i++)
  602. {
  603. if(method2->m_pBinaryArray[i][j] == 1) count++;
  604. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  605. m_tResPixelArray[i][j].rgbRed  = 255;
  606. }
  607. for(i=m_nWndHeight-1; i>=m_nWndHeight-count;i--)
  608. {
  609. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  610. m_tResPixelArray[i][j].rgbRed  = 0;
  611. }
  612. }
  613. MakeBitMap();
  614. SetCursor(LoadCursor(NULL,IDC_ARROW));
  615. }
  617. ////////////////////////////////////////////////////////////////////////////////
  618. // 垂直方向的直方图
  619. ////////////////////////////////////////////////////////////////////////////////
  620. void CFaceDetectDlg::OnBtnHistogramV() 
  621. {
  622. if(!method1->m_bBinaryReady)
  623. {
  624. AfxMessageBox("请先计算二值图");
  625. return;
  626. }
  628. m_bShowFace = false;
  629. SetCursor(LoadCursor(NULL,IDC_WAIT));
  630. for(int i=0; i<m_nWndHeight; i++)
  631. {
  632. int count = 0;
  633. for(int j=0; j<m_nWndWidth;  j++)
  634. {
  635. if(method1->m_pBinaryArray[i][j] == 1) count++;
  636. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  637. m_tResPixelArray[i][j].rgbRed  = 255;
  638. }
  639. for(j=0; j<count;  j++)
  640. {
  641. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  642. m_tResPixelArray[i][j].rgbRed  = 0;
  643. }
  644. }
  645. MakeBitMap();
  646. SetCursor(LoadCursor(NULL,IDC_ARROW));
  647. }
  649. ////////////////////////////////////////////////////////////////////////////////
  650. // 计算相似度
  651. ////////////////////////////////////////////////////////////////////////////////
  652. void CFaceDetectDlg::OnBtnLikehood() 
  653. {
  654. m_bShowFace = false;
  655. SetCursor(LoadCursor(NULL,IDC_WAIT));
  657. method1->CalLikeHood();
  658. for(int i=0; i<m_nWndHeight; i++)
  659. for(int j=0; j<m_nWndWidth;  j++)
  660. {
  661. m_tResPixelArray[i][j].rgbBlue = m_tResPixelArray[i][j].rgbGreen =
  662. m_tResPixelArray[i][j].rgbRed  = (int)(method1->m_pLikeliHoodArray[i][j]*255);
  663. }
  664. MakeBitMap();
  666. SetCursor(LoadCursor(NULL,IDC_ARROW));
  667. }
  669. ////////////////////////////////////////////////////////////////////////////////
  670. // 标记眼睛区域
  671. ////////////////////////////////////////////////////////////////////////////////
  672. void CFaceDetectDlg::OnBtnMarkEye() 
  673. {
  674. int i,j;
  675. if(!m_bFaceOK)
  676. {
  677. AfxMessageBox("请先确定脸部区域");
  678. return;
  679. }
  680. //左右眼的水平区域
  681. CPoint LeftEyeAreaH(-1,-1),RightEyeAreaH(-1,-1);
  682. CPoint LeftEyeAreaV(-1,-1),RightEyeAreaV(-1,-1);
  684. int nLeft,nRight,nTop,nBottom;
  686. nLeft = m_rFaceRegion.left-5  > 0 ? m_rFaceRegion.left-5:0;
  687. nRight = m_rFaceRegion.right+5 < m_nWndWidth? m_rFaceRegion.right+5:m_nWndWidth-1;
  688. nTop = m_rFaceRegion.top-5   > 0 ? m_rFaceRegion.top-5:0;
  689. nBottom = m_rFaceRegion.bottom+5< m_nWndHeight?m_rFaceRegion.bottom+5:m_nWndHeight-1;
  690. //边缘检查
  691. DoLOG(nLeft,nRight,nTop,nBottom,m_tOriPixelArray,m_tResPixelArray);
  693. ///////////////////////////////////
  694. //确认两个眼睛的水平区域
  695.     //////////////////////////////////
  696. int nSlidWinWidth  = (m_rFaceRegion.right - m_rFaceRegion.left)/6/2;
  697. int nSlidWinHeight = (m_rFaceRegion.bottom - m_rFaceRegion.top)/15/2;
  698. int nMidFaceH = (m_rFaceRegion.right+m_rFaceRegion.left)/2;
  699. int nMidFaceV = (m_rFaceRegion.bottom+m_rFaceRegion.top)/2;
  701. int *tempArray = new int[m_nWndWidth]; 
  702. for(i = 0; i<m_nWndWidth; i++) tempArray[i] = 0;
  704. for(i=nMidFaceV-nSlidWinHeight; i > m_rFaceRegion.top+6*nSlidWinHeight; i--)
  705. for(j=m_rFaceRegion.left+nSlidWinWidth; j<m_rFaceRegion.right-nSlidWinWidth; j++)
  706. {
  707. int count = 0;
  708. for(int p= -nSlidWinHeight ;p<nSlidWinHeight;p++)
  709. for(int q= -nSlidWinWidth ;q<nSlidWinWidth;q++)
  710. {
  711. if(m_tResPixelArray[i+p][j+q].rgbRed == 0) count++;
  712. }
  713. if(count >= nSlidWinWidth*nSlidWinHeight/3)
  714. {
  715. m_tResPixelArray[i][j].rgbRed = 255;
  716. tempArray[j] ++;
  717. }
  718. }
  720. MakeBitMap();
  721. AfxMessageBox("眼睛的区域鉴别");
  723. CList<CPoint,CPoint&> myList1(sizeof(CPoint));
  724. CList<CPoint,CPoint&> myList2(sizeof(CPoint));
  725. int flag = 0;
  726. CPoint tPoint(-1,-1);
  727. for(i = 0; i<m_nWndWidth; i++)
  728. {
  729. if(tempArray[i] > 0 && flag ==0)
  730. {
  731. tPoint.x = i;
  732. flag = 1;
  733. }
  734. if(tempArray[i] == 0 && flag ==1)
  735. {
  736. tPoint.y = i;
  737. myList1.AddTail(tPoint);
  738. flag = 0;
  739. }
  740. }
  741. delete tempArray;
  742. //去掉长度太小的候选者
  743. for(i=0; i<myList1.GetCount();i++)
  744. {
  745. CPoint temp(-1,-1);
  746. temp = myList1.GetAt(myList1.FindIndex(i));
  747. int minVal = (m_rFaceRegion.right - m_rFaceRegion.left)/20;
  748. if((temp.y-temp.x)>=minVal)
  749. myList2.AddTail(temp);
  750. }
  751. myList1.RemoveAll();
  752.     //合并相邻很紧的区域
  753. bool quit = 1;
  754. while(quit)
  755. {
  756. bool doJoin = false;
  757. for(int i=0; i<myList2.GetCount()-1;i++)
  758. {
  759. CPoint temp1(-1,-1),temp2(-1,-1);
  760. temp1 = myList2.GetAt(myList2.FindIndex(i));
  761. temp2 = myList2.GetAt(myList2.FindIndex(i+1));
  762. if((temp2.x-temp1.y)<=(m_rFaceRegion.right - m_rFaceRegion.left)/40)
  763. {
  764. temp1.y = temp2.y;
  765. myList2.RemoveAt(myList2.FindIndex(i));
  766. myList2.RemoveAt(myList2.FindIndex(i));
  767. if(i == 0) myList2.AddHead(temp1);
  768. else     myList2.InsertAfter(myList2.FindIndex(i-1),temp1);
  769. doJoin = true;
  770. break;
  771. }
  772. }
  773. if(!doJoin) quit = 0;
  774. }
  776. //没有找到眼睛区域
  777. if(myList2.GetCount()<2) 
  778. {
  779. CPoint t=myList2.GetHead();
  780. if((t.y-t.x)>(m_rFaceRegion.right - m_rFaceRegion.left)/2)
  781. {
  782. LeftEyeAreaH.x = t.x; 
  783. LeftEyeAreaH.y = t.x+(t.y-t.x)/3; 
  784. RightEyeAreaH.x = t.y-(t.y-t.x)/3;
  785. RightEyeAreaH.y = t.y; 
  786. }
  787. else
  788. {
  789. AfxMessageBox("确认眼睛位置失败，请手动标定");
  790. return;
  791. }
  792. }
  793. //仅有两个区域
  794. else if(myList2.GetCount()==2)
  795. {
  796. LeftEyeAreaH = myList2.GetHead();
  797. RightEyeAreaH = myList2.GetTail();
  798. }
  799. else  //多于两个区域
  800. {
  801. int ldis = -100000;
  802. int rdis = 100000;
  803. for(i=0; i<myList2.GetCount();i++)
  804. {
  805. CPoint temp(-1,-1);
  806. temp = myList2.GetAt(myList2.FindIndex(i));
  807. //右眼
  808. if((temp.x+temp.y)/2 > nMidFaceH)
  809. {
  810. if(((temp.x+temp.y)/2-nMidFaceH)<rdis)
  811. {
  812. rdis = (temp.x+temp.y)/2-nMidFaceH;
  813. RightEyeAreaH = temp;
  814. }
  815. }
  816. //左眼
  817. else
  818. {
  819. if(((temp.x+temp.y)/2-nMidFaceH)>ldis)
  820. {
  821. ldis = (temp.x+temp.y)/2-nMidFaceH;
  822. LeftEyeAreaH = temp;
  823. }
  824. }
  825. }
  826. }
  827. myList2.RemoveAll();
  828. ///////////////////////////////////
  829. //确认两个眼睛的垂直区域
  830.     //////////////////////////////////
  831. //左眼
  832. if(LeftEyeAreaH != CPoint(-1,-1))
  833. {
  834. int *tArray = new int[m_nWndHeight]; 
  835. int i,j;
  836. for(i = 0; i<m_nWndHeight; i++) tArray[i] = 0;
  838. for(i=nMidFaceV-nSlidWinHeight; i > m_rFaceRegion.top+6*nSlidWinHeight; i--)
  839. for(j=LeftEyeAreaH.x; j<=LeftEyeAreaH.y;j++)
  840. if(m_tResPixelArray[i][j].rgbRed == 255 && m_tResPixelArray[i][j].rgbGreen == 0)
  841. tArray[i] ++;
  843. CList<CPoint,CPoint&> myListA(sizeof(CPoint));
  844. CList<CPoint,CPoint&> myListB(sizeof(CPoint));
  845. int flag = 0;
  846. CPoint tPoint(-1,-1);
  847. for(i = nMidFaceV-nSlidWinHeight; i > m_rFaceRegion.top+6*nSlidWinHeight; i--)
  848. {
  849. if(tArray[i] > 0 && flag ==0)
  850. {
  851. tPoint.x = i;
  852. flag = 1;
  853. }
  854. if(tArray[i] == 0 && flag ==1)
  855. {
  856. tPoint.y = i;
  857. myListA.AddTail(tPoint);
  858. flag = 0;
  859. }
  860. }
  861. delete tArray;
  862. //去掉长度太小的候选者
  863. for(i=0; i<myListA.GetCount();i++)
  864. {
  865. CPoint temp(-1,-1);
  866. temp = myListA.GetAt(myListA.FindIndex(i));
  867. int minVal = (m_rFaceRegion.bottom - m_rFaceRegion.top)/100;
  868. if((temp.x-temp.y)>=minVal)
  869. myListB.AddTail(temp);
  870. }
  871. myListA.RemoveAll();
  872. //合并相邻很紧的区域
  873. bool quit = 1;
  874. while(quit)
  875. {
  876. bool doJoin = false;
  877. for(int i=0; i<myListB.GetCount()-1;i++)
  878. {
  879. CPoint temp1(-1,-1),temp2(-1,-1);
  880. temp1 = myListB.GetAt(myListB.FindIndex(i));
  881. temp2 = myListB.GetAt(myListB.FindIndex(i+1));
  882. if((temp1.y-temp2.x)<=(m_rFaceRegion.bottom - m_rFaceRegion.top)/100)
  883. {
  884. temp1.y = temp2.y;
  885. myListB.RemoveAt(myListB.FindIndex(i));
  886. myListB.RemoveAt(myListB.FindIndex(i));
  887. if(i == 0) myListB.AddHead(temp1);
  888. else     myListB.InsertAfter(myListB.FindIndex(i-1),temp1);
  889. doJoin = true;
  890. break;
  891. }
  892. }
  893. if(!doJoin) quit = 0;
  894. }
  895. if(myListB.GetCount()==0)
  896. {
  897. AfxMessageBox("无法确定左眼的位置");
  898. }
  899. else
  900. {
  901. LeftEyeAreaV = myListB.GetHead();
  903. double sumX = 0.0;
  904. double sumY = 0.0;
  905. int sum = 0;
  906. m_LeftEyeLeftCorner.x = 100000;
  907. m_LeftEyeRightCorner.x = -1;
  909. for(i=LeftEyeAreaV.x; i>= LeftEyeAreaV.y;i--)
  910. for(j=LeftEyeAreaH.x; j<=LeftEyeAreaH.y;j++)
  911. if(m_tResPixelArray[i][j].rgbGreen == 0)
  912. {
  913. if(j<m_LeftEyeLeftCorner.x)
  914. {
  915. m_LeftEyeLeftCorner.x = j;
  916. m_LeftEyeLeftCorner.y = i;
  917. }
  918. if(j>m_LeftEyeRightCorner.x)
  919. {
  920. m_LeftEyeRightCorner.x = j;
  921. m_LeftEyeRightCorner.y = i;
  922. }
  923. sumX += j;
  924. sumY += i;
  925. sum++;
  926. }
  928. m_LeftEye.x = (int)(sumX/sum);
  929. m_LeftEye.y = (int)(sumY/sum);
  931. m_bLeftEyeOK = TRUE;
  932. m_bLeftEyeLeftCornerOK = TRUE;
  933. m_bLeftEyeRightCornerOK =TRUE;
  934. }
  935. myListB.RemoveAll();
  936. }
  937. //右眼
  938. if(RightEyeAreaH != CPoint(-1,-1))
  939. {
  940. int *tArray = new int[m_nWndHeight]; 
  941. int i,j;
  942. for(i = 0; i<m_nWndHeight; i++) tArray[i] = 0;
  944. for(i=nMidFaceV-nSlidWinHeight; i > m_rFaceRegion.top+6*nSlidWinHeight; i--)
  945. for(j=RightEyeAreaH.x; j<=RightEyeAreaH.y;j++)
  946. if(m_tResPixelArray[i][j].rgbRed == 255 && m_tResPixelArray[i][j].rgbGreen == 0)
  947. tArray[i] ++;
  949. CList<CPoint,CPoint&> myListA(sizeof(CPoint));
  950. CList<CPoint,CPoint&> myListB(sizeof(CPoint));
  951. int flag = 0;
  952. CPoint tPoint(-1,-1);
  953. for(i = nMidFaceV-nSlidWinHeight; i > m_rFaceRegion.top+6*nSlidWinHeight; i--)
  954. {
  955. if(tArray[i] > 0 && flag ==0)
  956. {
  957. tPoint.x = i;
  958. flag = 1;
  959. }
  960. if(tArray[i] == 0 && flag ==1)
  961. {
  962. tPoint.y = i;
  963. myListA.AddTail(tPoint);
  964. flag = 0;
  965. }
  966. }
  967. delete tArray;
  968. //去掉长度太小的候选者
  969. for(i=0; i<myListA.GetCount();i++)
  970. {
  971. CPoint temp(-1,-1);
  972. temp = myListA.GetAt(myListA.FindIndex(i));
  973. int minVal = (m_rFaceRegion.bottom - m_rFaceRegion.top)/100;
  974. if((temp.x-temp.y)>=minVal)
  975. myListB.AddTail(temp);
  976. }
  977. myListA.RemoveAll();
  978. //合并相邻很紧的区域
  979. bool quit = 1;
  980. while(quit)
  981. {
  982. bool doJoin = false;
  983. for(int i=0; i<myListB.GetCount()-1;i++)
  984. {
  985. CPoint temp1(-1,-1),temp2(-1,-1);
  986. temp1 = myListB.GetAt(myListB.FindIndex(i));
  987. temp2 = myListB.GetAt(myListB.FindIndex(i+1));
  988. if((temp1.y-temp2.x)<=(m_rFaceRegion.bottom - m_rFaceRegion.top)/50)
  989. {
  990. temp1.y = temp2.y;
  991. myListB.RemoveAt(myListB.FindIndex(i));
  992. myListB.RemoveAt(myListB.FindIndex(i));
  993. if(i == 0) myListB.AddHead(temp1);
  994. else     myListB.InsertAfter(myListB.FindIndex(i-1),temp1);
  995. doJoin = true;
  996. break;
  997. }
  998. }
  999. if(!doJoin) quit = 0;
  1000. }
  1001. if(myListB.GetCount()==0)
  1002. {
  1003. AfxMessageBox("无法确定右眼的位置");
  1004. }
  1005. else
  1006. {
  1007. if(myListB.GetCount()==1)
  1008. RightEyeAreaV = myListB.GetHead();
  1009. else
  1010. {
  1011. CPoint tt =  myListB.GetHead();
  1012. int index = myListB.GetCount();
  1013. while(tt.y > LeftEyeAreaV.x && index > 0)
  1014. {
  1015. index --;
  1016. tt = myListB.GetAt(myListB.FindIndex(myListB.GetCount()-index)); 
  1017. }
  1018. RightEyeAreaV = tt;
  1019. }
  1021. double sumX = 0.0;
  1022. double sumY = 0.0;
  1023. int sum = 0;
  1024. m_RightEyeLeftCorner.x = 100000;
  1025. m_RightEyeRightCorner.x = -1;
  1027. for(i=RightEyeAreaV.x; i>=RightEyeAreaV.y;i--)
  1028. for(j=RightEyeAreaH.x; j<=RightEyeAreaH.y;j++)
  1029. if(m_tResPixelArray[i][j].rgbGreen == 0)
  1030. {
  1031. if(j<m_RightEyeLeftCorner.x)
  1032. {
  1033. m_RightEyeLeftCorner.x = j;
  1034. m_RightEyeLeftCorner.y = i;
  1035. }
  1036. if(j>m_RightEyeRightCorner.x)
  1037. {
  1038. m_RightEyeRightCorner.x = j;
  1039. m_RightEyeRightCorner.y = i;
  1040. }
  1041. sumX += j;
  1042. sumY += i;
  1043. sum++;
  1044. }
  1045. m_RightEye.x = (int)(sumX/sum);
  1046. m_RightEye.y = (int)(sumY/sum);
  1048. m_bRightEyeOK = TRUE;
  1049. m_bRightEyeLeftCornerOK = TRUE;
  1050. m_bRightEyeRightCornerOK =TRUE;
  1052. }
  1053. myListB.RemoveAll();
  1054. }
  1055. CopyBitMap(m_tResPixelArray,m_tOriPixelArray);
  1056. MakeBitMap();
  1057. }
  1059. ////////////////////////////////////////////////////////////////////////////////
  1060. // 第一种方法标记脸部区域
  1061. ////////////////////////////////////////////////////////////////////////////////
  1062. void CFaceDetectDlg::OnBtnMarkFace1() 
  1063. {
  1064. if(!method1->m_bBinaryReady)
  1065. {
  1066. AfxMessageBox("请先计算二值化图！");
  1067. return;
  1068. }
  1069. m_bShowFace = true;
  1070. SetCursor(LoadCursor(NULL,IDC_WAIT));
  1071. int *temp = new int[m_nWndWidth];
  1072. int max = 0;
  1073. int pos = -1;
  1074. for(int j=0; j<m_nWndWidth;  j++)
  1075. {
  1076. int count = 0;
  1077. for(int i=0; i<m_nWndHeight; i++)
  1078. {
  1079. if(method1->m_pBinaryArray[i][j] == 1) count++;
  1080. }
  1081. temp[j] = count;
  1082. if(count > max)
  1083. {
  1084. max = count;
  1085. pos = j;
  1086. }
  1087. }
  1088. int left,right,l,top,bottom;
  1089. for(l=pos; l>=0; l--)
  1090. {
  1091. if(temp[l]<max*0.2||l==0)
  1092. {
  1093. left = l;
  1094. break;
  1095. }
  1096. }
  1097. for(l=pos; l<m_nWndWidth; l++)
  1098. {
  1099. if(temp[l]<max*0.3||l==m_nWndWidth-1)
  1100. {
  1101. right = l;
  1102. break;
  1103. }
  1104. }
  1105. for(int i=0; i<m_nWndHeight; i++)
  1106. {
  1107. int count = 0;
  1108. for(l = left;l<=right;l++)
  1109. {
  1110. if(method1->m_pBinaryArray[i][l] == 1) count++;
  1111. }
  1112. if(count>=(right-left)*0.5)
  1113. {
  1114. top = i;
  1115. break;
  1116. }
  1117. }
  1118. bottom = (int)(top+(right-left)*1.5)>=m_nWndHeight? m_nWndHeight-1:(int)(top+(right-left)*1.5);
  1120. CopyBitMap(m_tResPixelArray,m_tOriPixelArray);
  1122. for(i=top;i<=bottom;i++)
  1123. {
  1124. m_tResPixelArray[i][left].rgbBlue=255;
  1125. m_tResPixelArray[i][left].rgbGreen = m_tResPixelArray[i][left].rgbRed = 0;
  1126. m_tResPixelArray[i][right].rgbBlue=255;
  1127. m_tResPixelArray[i][right].rgbGreen = m_tResPixelArray[i][right].rgbRed = 0;
  1128. }
  1129. for(j=left;j<=right;j++)
  1130. {
  1131. m_tResPixelArray[top][j].rgbBlue=255;
  1132. m_tResPixelArray[top][j].rgbGreen = m_tResPixelArray[top][j].rgbRed = 0;
  1133. m_tResPixelArray[bottom][j].rgbBlue=255;
  1134. m_tResPixelArray[bottom][j].rgbGreen = m_tResPixelArray[bottom][j].rgbRed = 0;
  1135. }
  1136. MakeBitMap();
  1137. SetCursor(LoadCursor(NULL,IDC_ARROW));
  1139. if(m_bFaceOK)
  1140. {
  1141. ReplaceDlg dlg;
  1142. if(dlg.DoModal()==IDOK)
  1143. {
  1144. CopyBitMap(m_tResPixelArray,m_tOriPixelArray);
  1145. CRect rect(left,top,right,bottom);
  1146. m_rFaceRegion = rect;
  1147. MakeBitMap();
  1148. }
  1149. }
  1150. else
  1151. {
  1152. m_bFaceOK = true;
  1153. CopyBitMap(m_tResPixelArray,m_tOriPixelArray);
  1154. CRect rect(left,top,right,bottom);
  1155. m_rFaceRegion = rect;
  1156. MakeBitMap();
  1157. }
  1158. }
  1160. ////////////////////////////////////////////////////////////////////////////////
  1161. // 第二种方法标记脸部区域
  1162. ////////////////////////////////////////////////////////////////////////////////
  1163. void CFaceDetectDlg::OnBtnMarkFace2() 
  1164. {
  1165. if(!method2->m_bBinaryOK)
  1166. {
  1167. AfxMessageBox("请先计算二值化图");
  1168. return;
  1169. }
  1170. m_bShowFace = true;
  1171. SetCursor(LoadCursor(NULL,IDC_WAIT));
  1172. int *numR,*numB, i, j, left,right,top,bottom;
  1173. int maxnumR = 0, maxnumB = 0;
  1175. numR = new int[m_nWndWidth];
  1176. numB = new int[m_nWndWidth];
  1177. for(j=0;j<m_nWndWidth ;j++)
  1178. {
  1179. int countR = 0, countB = 0;
  1180. for(i=0;i<m_nWndHeight;i++)
  1181. {
  1182. if(method2->m_pBinaryArray[i][j] == 0)
  1183. countR++;
  1184. if(method2->m_pBinaryArray[i][j] == 1)
  1185. countB++;
  1186. }
  1187. if(countR > maxnumR) maxnumR = countR;
  1188. numR[j] = countR;
  1189. if(countB > maxnumB) maxnumB = countB;
  1190. numB[j] = countB;
  1191. }
  1192. CList<CPoint,CPoint> myListR(sizeof(CPoint));
  1193. CList<CPoint,CPoint> myListB(sizeof(CPoint));
  1194. CPoint tempR,tempB;
  1195. int flagR = 0,flagB = 0;
  1196. for(j=0;j<m_nWndWidth ;j++)
  1197. {
  1198. if(flagR == 0)
  1199. {
  1200. if(numR[j]>maxnumR/2)
  1201. {
  1202. flagR = 1;
  1203. tempR.x = j;
  1204. }
  1205. }
  1206. else
  1207. {
  1208. if(numR[j]<=maxnumR/2 || j==m_nWndWidth-1)
  1209. {
  1210. flagR = 0;
  1211. tempR.y = j;
  1212. myListR.AddTail(tempR);
  1213. }
  1214. }
  1215. if(flagB == 0)
  1216. {
  1217. if(numB[j]>maxnumB/5)
  1218. {
  1219. flagB = 1;
  1220. tempB.x = j;
  1221. }
  1222. }
  1223. else
  1224. {
  1225. if(numB[j]<=maxnumB/5 || j==m_nWndWidth-1)
  1226. {
  1227. flagB = 0;
  1228. tempB.y = j;
  1229. if(myListB.GetCount() > 1 && (tempB.x-myListB.GetTail().y)<20)
  1230. myListB.SetAt(myListB.GetTailPosition(),CPoint(myListB.GetTail().x,j));
  1231. else
  1232. myListB.AddTail(tempB);
  1233. }
  1234. }
  1236. }
  1237. if(numR!=NULL)delete numR;
  1238. if(numB!=NULL)delete numB;
  1240. int *hairmark, k;
  1241. hairmark = new int[m_nWndWidth];
  1242. for(j=0;j<m_nWndWidth ;j++) hairmark[j]=0;
  1244. for(k=0;k<myListB.GetCount();k++)
  1245. {
  1246. CPoint temp = myListB.GetAt(myListB.FindIndex(k));
  1247. if((temp.y-temp.x)>m_nWndWidth/10)
  1248. {
  1249. for(int t = temp.x;t<=temp.y;t++)
  1250. hairmark[t] = 1;
  1251. }
  1252. }
  1254. for(k=0;k<myListR.GetCount();k++)
  1255. {
  1256. CPoint temp = myListR.GetAt(myListR.FindIndex(k));
  1257. int templeft=-1;
  1258. int tempright=-1;
  1260. if((temp.y-temp.x)>m_nWndWidth/10)
  1261. {
  1262. for(int t=temp.x;t<=temp.y;t++)
  1263. {
  1264. if(hairmark[t]==1)
  1265. {
  1266. int endpos=t+(temp.y-temp.x)/5;
  1267. if(endpos > temp.y)endpos = temp.y;
  1268. int yes = 1;
  1269. for(int q=t;q<=endpos;q++)
  1270. if(hairmark[q]==0) yes = 0;
  1271. if(yes == 1)
  1272. {
  1273. templeft = t;
  1274. break;
  1275. }
  1276. }
  1277. }
  1278. for(int p=temp.y;p>=temp.x;p--)
  1279. {
  1280. if(hairmark[p]==1)
  1281. {
  1282. int beginpos=p-(temp.y-temp.x)/5;
  1283. if(beginpos < temp.x)beginpos = temp.x;
  1284. int yes = 1;
  1285. for(int q=p;q>=beginpos;q--)
  1286. if(hairmark[q]==0) yes = 0;
  1287. if(yes == 1)
  1288. {
  1289. tempright = p;
  1290. break;
  1291. }
  1292. }
  1293. }
  1294. }
  1295. if(templeft!=-1 && tempright!=-1)
  1296. {
  1297. left = templeft;
  1298. right = tempright;
  1299. break;
  1300. }
  1301. }
  1302. if(hairmark !=NULL) delete hairmark;
  1303. myListR.RemoveAll();
  1304. myListB.RemoveAll();
  1306. if(left-m_nWndWidth/50>0) left-=(int)m_nWndWidth/50;
  1307. else left = 0;
  1308. if(right+m_nWndWidth/40>m_nWndWidth) right=m_nWndWidth-1;
  1309. else right += m_nWndWidth/40;
  1311. for(i=0; i<m_nWndHeight; i++)
  1312. {
  1313. int count = 0;
  1314. for(int l = left;l<=right;l++)
  1315. {
  1316. if(method2->m_pBinaryArray[i][l] == 0) count++;
  1317. }
  1318. if(count>=(right-left)*0.5)
  1319. {
  1320. top = i;
  1321. break;
  1322. }
  1323. }
  1324. bottom = (int)(top+(right-left)*1.5)>=m_nWndHeight? m_nWndHeight-1:(int)(top+(right-left)*1.5);
  1326. CopyBitMap(m_tResPixelArray,m_tOriPixelArray);
  1328. for(i=top;i<=bottom;i++)
  1329. {
  1330. m_tResPixelArray[i][left].rgbBlue=255;
  1331. m_tResPixelArray[i][left].rgbGreen = m_tResPixelArray[i][left].rgbRed = 0;
  1332. m_tResPixelArray[i][right].rgbBlue=255;
  1333. m_tResPixelArray[i][right].rgbGreen = m_tResPixelArray[i][right].rgbRed = 0;
  1334. }
  1335. for(j=left;j<=right;j++)
  1336. {
  1337. m_tResPixelArray[top][j].rgbBlue=255;
  1338. m_tResPixelArray[top][j].rgbGreen = m_tResPixelArray[top][j].rgbRed = 0;
  1339. m_tResPixelArray[bottom][j].rgbBlue=255;
  1340. m_tResPixelArray[bottom][j].rgbGreen = m_tResPixelArray[bottom][j].rgbRed = 0;
  1341. }
  1342. MakeBitMap();
  1343. SetCursor(LoadCursor(NULL,IDC_ARROW));
  1345. if(m_bFaceOK)
  1346. {
  1347. ReplaceDlg dlg;
  1348. if(dlg.DoModal()==IDOK)
  1349. {
  1350. CopyBitMap(m_tResPixelArray,m_tOriPixelArray);
  1351. CRect rect(left,top,right,bottom);
  1352. m_rFaceRegion = rect;
  1353. MakeBitMap();
  1354. }
  1355. }
  1356. else
  1357. {
  1358. m_bFaceOK = true;
  1359. CopyBitMap(m_tResPixelArray,m_tOriPixelArray);
  1360. CRect rect(left,top,right,bottom);
  1361. m_rFaceRegion = rect;
  1362. MakeBitMap();
  1363. }
  1364. }
  1366. ////////////////////////////////////////////////////////////////////////////////
  1367. // 标记嘴巴
  1368. ////////////////////////////////////////////////////////////////////////////////
  1369. void CFaceDetectDlg::OnBtnMarkMouse() 
  1370. {
  1371. int i,j;
  1372. if(!(m_bLeftEyeOK&&m_bRightEyeOK))
  1373. {
  1374. AfxMessageBox("请先确定眼睛");
  1375. return;
  1376. }
  1377. //左右眼的水平区域
  1378. int nLeft,nRight,nTop,nBottom;
  1380. nLeft = m_rFaceRegion.left-5  > 0 ? m_rFaceRegion.left-5:0;
  1381. nRight = m_rFaceRegion.right+5 < m_nWndWidth? m_rFaceRegion.right+5:m_nWndWidth-1;
  1382. nTop = m_rFaceRegion.top-5   > 0 ? m_rFaceRegion.top-5:0;
  1383. nBottom = m_rFaceRegion.bottom+5< m_nWndHeight?m_rFaceRegion.bottom+5:m_nWndHeight-1;
  1385. SetPixelArray(m_tResPixelArray,0);
  1386. for(i=nTop; i<=nBottom; i++)
  1387. for(j=nLeft; j<=nRight; j++)
  1388. {
  1389. BYTE R,G,B;
  1390. double temp,dlta;
  1391. R = m_tOriPixelArray[i][j].rgbRed;
  1392. G = m_tOriPixelArray[i][j].rgbGreen;
  1393. B = m_tOriPixelArray[i][j].rgbBlue;
  1394. if((R==G) && (G==B)) temp = 0;
  1395. else temp = 0.5*(2*R-G-B)/sqrt((R-G)*(R-G)+(R-B)*(G-B));
  1396. dlta = acos(temp);
  1397. if(dlta < 0.2)
  1398. {
  1399. m_tResPixelArray[i][j].rgbRed = 255;  
  1400. }
  1401. else m_tResPixelArray[i][j].rgbRed = 0;
  1402. }
  1404. MakeBitMap();
  1405. AfxMessageBox("嘴的肤色鉴定");
  1406. //双目斜角
  1407. double tanThta;
  1408. if(m_RightEye.y == m_LeftEye.y) tanThta = 0;
  1409. else tanThta = (m_RightEye.y - m_LeftEye.y)/(m_RightEye.x - m_LeftEye.x);
  1410. //双目距离
  1411. int EyesDis = (m_RightEye.x-m_LeftEye.x)*(m_RightEye.x-m_LeftEye.x);
  1412. EyesDis += (m_RightEye.y-m_LeftEye.y)*(m_RightEye.y-m_LeftEye.y);
  1413. EyesDis = (int)sqrt(EyesDis);
  1414. //双目平均高度
  1415. int EyeV    = (m_RightEye.y + m_LeftEye.y)/2;
  1416. //可能的嘴的区域
  1417. int MouthUp   = (EyeV+1.0*EyesDis) > nBottom ? nBottom:(int)(EyeV+1.0*EyesDis);
  1418. int MouthDown = (EyeV+1.5*EyesDis) > nBottom ? nBottom:(int)(EyeV+1.5*EyesDis);
  1421. int* Y_Arry = new int[MouthDown-MouthUp];
  1422. for(i =0 ;i < MouthDown-MouthUp ;i++) Y_Arry[i] = 0;
  1423. int* X_Arry = new int[EyesDis];
  1424. for(i =0 ;i < EyesDis ;i++) X_Arry[i] = 0;
  1425. for(i = MouthUp ; i < MouthDown; i++)
  1426. for(j = m_LeftEye.x; j< m_RightEye.x; j++)
  1427. {
  1428. if(m_tResPixelArray[i][j].rgbRed == 255)
  1429. {
  1430. Y_Arry[i-MouthUp] ++;
  1431. X_Arry[j-m_LeftEye.x] ++;
  1432. }
  1433. }
  1435. int maxY = 0;
  1436. for(i =0 ;i < MouthDown-MouthUp ;i++)
  1437. {
  1438. if(Y_Arry[i]>maxY)
  1439. {
  1440. maxY = Y_Arry[i];
  1441. m_MidMouth.y =  i+MouthUp - (MouthDown-MouthUp)/10;
  1442. }
  1443. }
  1444. m_LeftMouthCorner.y  =(int)(m_MidMouth.y - tanThta*EyesDis/2);
  1445. m_RightMouthCorner.y =(int)(m_MidMouth.y + tanThta*EyesDis/2);
  1446. for(i =0 ;i < EyesDis ;i++)
  1447. {
  1448. if(X_Arry[i]>0)
  1449. {
  1450. m_LeftMouthCorner.x = i+m_LeftEye.x;
  1451. break;
  1452. }
  1453. }
  1454. for(i = EyesDis -1; i >=0 ;i--)
  1455. {
  1456. if(X_Arry[i]>0)
  1457. {
  1458. m_RightMouthCorner.x = m_LeftEye.x+i;
  1459. break;
  1460. }
  1461. }
  1463. //唇中点较薄
  1464. int min = 1000000;
  1465. for(i = (int)(3*EyesDis/7+0.5) ; i <= (int)(4*EyesDis/7+0.5);i++)
  1466. {
  1467. if(X_Arry[i]<min)
  1468. {
  1469. min = X_Arry[i];
  1470. m_MidMouth.x = m_LeftEye.x+i;
  1471. }
  1472. }
  1473. m_MidMouth.x = (int)(m_MidMouth.x+(m_LeftEye.x+EyesDis/2))/2;
  1475. m_bMidMouthOK = TRUE;
  1476. m_bLeftMouthCornerOK = TRUE;
  1477. m_bRightMouthCornerOK = TRUE;
  1479. CopyBitMap(m_tResPixelArray,m_tOriPixelArray);
  1480. MakeBitMap();
  1481. }
  1483. ////////////////////////////////////////////////////////////////////////////////
  1484. // 标记鼻子
  1485. ////////////////////////////////////////////////////////////////////////////////
  1486. void CFaceDetectDlg::OnBtnMarkNose() 
  1487. {
  1488. int i,j;
  1489. if(!(m_bLeftEyeOK&&m_bRightEyeOK))
  1490. {
  1491. AfxMessageBox("请先确定眼睛");
  1492. return;
  1493. }
  1494. //左右眼的水平区域
  1495. int nLeft,nRight,nTop,nBottom;
  1497. nLeft = m_rFaceRegion.left-5  > 0 ? m_rFaceRegion.left-5:0;
  1498. nRight = m_rFaceRegion.right+5 < m_nWndWidth? m_rFaceRegion.right+5:m_nWndWidth-1;
  1499. nTop = m_rFaceRegion.top-5   > 0 ? m_rFaceRegion.top-5:0;
  1500. nBottom = m_rFaceRegion.bottom+5< m_nWndHeight?m_rFaceRegion.bottom+5:m_nWndHeight-1;
  1502. SetPixelArray(m_tResPixelArray,0);
  1503. for(i=nTop; i<=nBottom; i++)
  1504. for(j=nLeft; j<=nRight; j++)
  1505. {
  1506. double Y;
  1507. Y = 0.30*m_tOriPixelArray[i][j].rgbRed+0.59*m_tOriPixelArray[i][j].rgbGreen
  1508. +0.11*m_tOriPixelArray[i][j].rgbBlue;
  1509. if(Y<100)
  1510. {
  1511. m_tResPixelArray[i][j].rgbRed = 255;  
  1512. }
  1513. else m_tResPixelArray[i][j].rgbRed = 0;
  1514. }
  1516. MakeBitMap();
  1517. AfxMessageBox("鼻子的肤色鉴定");
  1518. //双目斜角
  1519. double tanThta;
  1520. if(m_RightEye.y == m_LeftEye.y) tanThta = 0;
  1521. else tanThta = (m_RightEye.y - m_LeftEye.y)/(m_RightEye.x - m_LeftEye.x);
  1522. //双目距离
  1523. int EyesDis = (m_RightEye.x-m_LeftEye.x)*(m_RightEye.x-m_LeftEye.x);
  1524. EyesDis += (m_RightEye.y-m_LeftEye.y)*(m_RightEye.y-m_LeftEye.y);
  1525. EyesDis = (int)sqrt(EyesDis);
  1526. //双目平均高度
  1527. int EyeV    = (m_RightEye.y + m_LeftEye.y)/2;
  1528. //可能的鼻子的区域
  1529. int NoseUp   = (EyeV+0.5*EyesDis) > nBottom ? nBottom:(int)(EyeV+0.5*EyesDis);
  1530. int NoseDown = (EyeV+0.8*EyesDis) > nBottom ? nBottom:(int)(EyeV+0.8*EyesDis);
  1532. int* Y_Arry = new int[NoseDown-NoseUp];
  1533. for(i =0 ;i < NoseDown-NoseUp ;i++) Y_Arry[i] = 0;
  1535. int* X_Arry = new int[EyesDis];
  1536. for(i =0 ;i < EyesDis ;i++) X_Arry[i] = 0;
  1538. for(i = NoseUp ; i < NoseDown; i++)
  1539. for(j = m_LeftEye.x+EyesDis/5; j< m_RightEye.x-EyesDis/5; j++)
  1540. {
  1541. if(m_tResPixelArray[i][j].rgbRed == 255)
  1542. {
  1543. Y_Arry[i-NoseUp] ++;
  1544. X_Arry[j-m_LeftEye.x] ++;
  1545. }
  1546. }
  1548. int maxY = 0;
  1549. for(i =0 ;i < NoseDown-NoseUp ;i++)
  1550. {
  1551. if(Y_Arry[i]>maxY)
  1552. {
  1553. maxY = Y_Arry[i];
  1554. m_MidNose.y =  i+NoseUp;
  1555. }
  1556. }
  1557. m_LeftNostril.y  =(int)(m_MidNose.y - tanThta*EyesDis/2);
  1558. m_RightNostril.y =(int)(m_MidNose.y + tanThta*EyesDis/2);
  1559. for(i =0 ;i < EyesDis ;i++)
  1560. {
  1561. if(X_Arry[i]>0)
  1562. {
  1563. m_LeftNostril.x = i+m_LeftEye.x;
  1564. break;
  1565. }
  1566. }
  1567. for(i = EyesDis-1; i >=0 ;i--)
  1568. {
  1569. if(X_Arry[i]>0)
  1570. {
  1571. m_RightNostril.x = i+m_LeftEye.x;
  1572. break;
  1573. }
  1574. }
  1575. //唇中点较薄
  1576. int min = 1000000;
  1577. for(i = (int)(EyesDis/3+0.5) ; i <= (int)(2*EyesDis/3+0.5);i++)
  1578. {
  1579. if(X_Arry[i]<min)
  1580. {
  1581. min = X_Arry[i];
  1582. m_MidNose.x = m_LeftEye.x+i;
  1583. }
  1584. }
  1585. m_MidNose.x = (m_MidNose.x+(m_LeftEye.x+EyesDis/2))/2;
  1587. m_bMidNoseOK = TRUE;
  1588. m_bLeftNostrilOK = TRUE;
  1589. m_bRightNostrilOK = TRUE;
  1591. CopyBitMap(m_tResPixelArray,m_tOriPixelArray);
  1592. MakeBitMap();
  1593. }
  1595. ////////////////////////////////////////////////////////////////////////////////
  1596. // 打开文件
  1597. ////////////////////////////////////////////////////////////////////////////////
  1598. void CFaceDetectDlg::OnBtnOpenfile() 
  1599. {
  1600. CAddSampleDlg FileDlg(TRUE, "", NULL, 
  1601. OFN_HIDEREADONLY | OFN_OVERWRITEPROMPT | OFN_ALLOWMULTISELECT, 
  1602. "BMP人脸图像(*.bmp)|*.bmp|所有文件(*.*)|*.*||",
  1603. AfxGetMainWnd());
  1605. CString strFile;
  1606. if (FileDlg.DoModal () != IDOK)
  1607. return;
  1609. POSITION pos = FileDlg.GetStartPosition();
  1610. strFile = FileDlg.GetNextPathName(pos);
  1611. m_pMainDib->Open(strFile);
  1613. m_nWndWidth = m_pMainDib->GetWidth();
  1614. m_nWndHeight= m_pMainDib->GetHeight();
  1615. m_sFileName = strFile;
  1617. m_rFaceRegion.left = m_rFaceRegion.right = m_rFaceRegion.top = m_rFaceRegion.bottom = 0;
  1619. m_bLeftEyeOK = m_bRightEyeOK = m_bLeftNostrilOK = m_bRightNostrilOK =
  1620. m_bLeftEyeLeftCornerOK = m_bLeftEyeRightCornerOK = m_bRightEyeLeftCornerOK = 
  1621. m_bRightEyeRightCornerOK = m_bLeftMouthCornerOK = m_bRightMouthCornerOK = false;
  1623. m_bMidMouthOK = m_bMidNoseOK = false;
  1625. m_LeftEye = m_RightEye = m_LeftEyeLeftCorner = m_LeftEyeRightCorner = 
  1626. m_LeftNostril = m_RightNostril = m_RightEyeLeftCorner = m_RightEyeRightCorner =
  1627. m_LeftMouthCorner = m_RightMouthCorner = m_MidMouth = m_MidNose = CPoint(-1,-1);
  1629. SetCursor(LoadCursor(NULL,IDC_WAIT));
  1630. //获取像素的值
  1631. LoadOriPixel(m_pMainDib);
  1632. MakeBitMap();
  1633. SetCursor(LoadCursor(NULL,IDC_ARROW));
  1634. }
  1636. ////////////////////////////////////////////////////////////////////////////////
  1637. // 画图
  1638. ////////////////////////////////////////////////////////////////////////////////
  1639. void CFaceDetectDlg::MyDraw()
  1640. {
  1641. CRect rc;
  1642. pDCShow->GetWindow()->GetClientRect(&rc);
  1643. pDCShow->Rectangle(&rc);
  1644. CDC dc;
  1645. CBitmap *pOldBitmap;
  1646. dc.CreateCompatibleDC(pDCShow);
  1647. pOldBitmap=dc.SelectObject(m_pResMap);
  1648. pDCShow->StretchBlt(0,0,m_nWndWidth,m_nWndHeight,&dc,0,0,m_nWndWidth,m_nWndHeight,SRCCOPY);
  1649. dc.SelectObject(pOldBitmap);
  1650. dc.DeleteDC();
  1651. }