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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 行业应用 > 交通/航空行业 > 查看源码
Bmptestview.cpp
资源名称：vehiclehsv.rar　[点击查看]
上传用户：cjd055
上传日期：2013-04-01
资源大小：608k
文件大小：65k
源码类别：

交通/航空行业

开发平台：

Visual C++

Bmptestview.cpp：源码内容
  1. // BmpTestView.cpp : implementation of the CBmpTestView class
  2. //
  3. #include "stdafx.h"
  4. #include "BmpTest.h"
  5. #include "stack.h"
  6. #include "BmpTestDoc.h"
  7. #include "BmpTestView.h"
  8. #include <math.h>
  10. #ifdef _DEBUG
  11. #define new DEBUG_NEW
  12. #undef THIS_FILE
  13. static char THIS_FILE[] = __FILE__;
  14. #endif
  15. #define MIN(a,b) (a<b? a:b)
  16. #define MAX(a,b) (a>b? a:b)
  17. #define NO_HUE -1;
  18. int m_Grade[256],m_Gradmax=0,m_Grade_Blue[256],m_Grade_Green[256],m_Grade_Red[256],b[1024][1024];
  19. extern void str(int quotient,char addr[15]);
  20. extern BOOL  Contour(CBmpTestDoc* pDoc,double *pAnew_mon,BOOL bmoment,double n,double *protate);
  21. extern BOOL  GCC(CBmpTestDoc* pDoc,double *n,POINT *gcc);
  22. extern void Sh_Rotate(CBmpTestDoc* pDoc,double ro_angle,double size,POINT *pav_cent);
  23. DISTANCE distbmp[50];
  25. /////////////////////////////////////////////////////////////////////////////
  26. // CBmpTestView
  28. IMPLEMENT_DYNCREATE(CBmpTestView, CScrollView)
  30. BEGIN_MESSAGE_MAP(CBmpTestView, CScrollView)
  31. //{{AFX_MSG_MAP(CBmpTestView)
  32. ON_COMMAND(ID_BT_BLACKNESS, OnBtBlackness)
  33. ON_COMMAND(ID_BT_DSTINVERT, OnBtDstinvert)
  34. ON_COMMAND(ID_BT_MERGECOPY, OnBtMergecopy)
  35. ON_COMMAND(ID_BT_MERGEPAINT, OnBtMergepaint)
  36. ON_COMMAND(ID_BT_NOTSRCCOPY, OnBtNotsrccopy)
  37. ON_COMMAND(ID_BT_NOTSRCERASE, OnBtNotsrcerase)
  38. ON_COMMAND(ID_BT_PATCOPY, OnBtPatcopy)
  39. ON_COMMAND(ID_BT_PATINVERT, OnBtPatinvert)
  40. ON_COMMAND(ID_BT_PATPAINT, OnBtPatpaint)
  41. ON_COMMAND(ID_BT_SRCAND, OnBtSrcand)
  42. ON_COMMAND(ID_BT_SRCCOPY, OnBtSrccopy)
  43. ON_COMMAND(ID_BT_SRCERASE, OnBtSrcerase)
  44. ON_COMMAND(ID_BT_SRCINVERT, OnBtSrcinvert)
  45. ON_COMMAND(ID_BT_SRCPAINT, OnBtSrcpaint)
  46. ON_COMMAND(ID_BT_WHITENESS, OnBtWhiteness)
  47. ON_COMMAND(ID_ZFT_SHOW, OnZftShow)
  48. ON_COMMAND(ID_ZFT_BALANCE, OnZftBalance)
  49. ON_COMMAND(ID_IMAGE_SHOW, OnImageShow)
  50. ON_COMMAND(ID_DTSM_IN, OnDtsmIn)
  51. ON_COMMAND(ID_DTSM_SHOW, OnDtsmShow)
  52. ON_COMMAND(ID_DTSM_OUT, OnDtsmOut)
  53. ON_COMMAND(ID_REC_FILT, OnRecFilt)
  54. ON_COMMAND(ID_MID_FILT, OnMidFilt)
  55. ON_COMMAND(ID_EDGE_ROBERT, OnEdgeRobert)
  56. ON_COMMAND(ID_EDGE_LAPLACE, OnEdgeLaplace)
  57. ON_COMMAND(ID_GAUSE_FILT, OnGauseFilt)
  58. ON_COMMAND(ID_EDGE_MARR, OnEdgeMarr)
  59. ON_COMMAND(ID_STRU_ADD, OnStruAdd)
  60. ON_COMMAND(ID_STRU_SUB, OnStruSub)
  61. ON_COMMAND(ID_STRU_OPEN, OnStruOpen)
  62. ON_COMMAND(ID_STRU_CLOSE, OnStruClose)
  63. ON_COMMAND(ID_EDIT_BACKUP, OnEditBackup)
  64. ON_COMMAND(ID_EDIT_BLUE, OnEditBlue)
  65. ON_COMMAND(ID_EDIT_GREEN, OnEditGreen)
  66. ON_COMMAND(ID_EDIT_RED, OnEditRed)
  67. ON_COMMAND(ID_EDGE_COLORPREWITT, OnEdgeColorprewitt)
  68. ON_COMMAND(ID_EDGE_COLORROBERT, OnEdgeColorrobert)
  69. ON_COMMAND(ID_MID_FILT_LINE, OnMidFiltLine)
  70. ON_COMMAND(ID_DENSITY_SCAN, OnDensityScan)
  71. ON_COMMAND(ID_ROW_DILATION, OnRowDilation)
  72. ON_COMMAND(ID_EDGE_COLORPREWITT1, OnEdgeColorprewitt1)
  73. ON_COMMAND(ID_EDGE_COLORROBERT1, OnEdgeColorrobert1)
  74. ON_COMMAND(ID_ROW_EROSION, OnRowErosion)
  75. ON_COMMAND(ID_CLOSE, OnClose)
  76. ON_COMMAND(ID_GET, OnGet)
  77. ON_COMMAND(ID_RELA_SCAN, OnRelaScan)
  78. ON_COMMAND(ID_DENSITY_SCAN_MID, OnDensityScanMid)
  79. ON_COMMAND(ID_AREA_MARK, OnAreaMark)
  80. ON_COMMAND(ID_LOCATING, OnLocating)
  81. ON_COMMAND(ID_VEHICLE_BITMAP, OnVehicleBitmap)
  82. ON_COMMAND(ID_BLUE_CUT, OnBlueCut)
  83. //}}AFX_MSG_MAP
  84. // Standard printing commands
  85. ON_COMMAND(ID_FILE_PRINT, CScrollView::OnFilePrint)
  86. ON_COMMAND(ID_FILE_PRINT_DIRECT, CScrollView::OnFilePrint)
  87. ON_COMMAND(ID_FILE_PRINT_PREVIEW, CScrollView::OnFilePrintPreview)
  88. END_MESSAGE_MAP()
  90. /////////////////////////////////////////////////////////////////////////////
  91. // CBmpTestView construction/destruction
  93. CBmpTestView::CBmpTestView()
  94. {
  95. // TODO: add construction code here
  96. m_Rop=SRCCOPY;
  97. }
  99. CBmpTestView::~CBmpTestView()
  100. {
  101. }
  103. BOOL CBmpTestView::PreCreateWindow(CREATESTRUCT& cs)
  104. {
  105. // TODO: Modify the Window class or styles here by modifying
  106. //  the CREATESTRUCT cs
  108. return CScrollView::PreCreateWindow(cs);
  109. }
  111. /////////////////////////////////////////////////////////////////////////////
  112. // CBmpTestView drawing
  114. void CBmpTestView::OnDraw(CDC* pDC)
  115. {
  116. CPen pen,*open;
  117. CBmpTestDoc* pDoc = GetDocument();
  118. ASSERT_VALID(pDoc);
  119. // TODO: add draw code for native data here
  120. //设置滚动视图大小，使之与Dib大小一致：
  121. CSize sizeTotal;
  122. sizeTotal.cx = pDoc->m_Dib.Width;
  123. sizeTotal.cy = pDoc->m_Dib.Height+100;
  124. SetScrollSizes(MM_TEXT, sizeTotal);
  125. //绘制Dib：
  126. if (pDoc->m_BM)
  127. {
  128. //在滚动视图中绘制背景：
  129.      CBitmap pat;
  130.     CBrush br,* obr;
  131. CRect rc;
  132. rc.top=0;
  133. rc.left=0;
  134. rc.right=sizeTotal.cx;
  135. rc.bottom=sizeTotal.cy;
  136. GetClientRect(&rc);  //获得矩形区域大小
  137. pat.LoadBitmap(IDB_BITMAP2); //IDB_BITMAP2->pat
  138. br.CreatePatternBrush(&pat); //创建模式刷, pat->br
  139. obr=pDC->SelectObject(&br);
  140. pDC->FillRect(&rc,&br); //用画刷填充整个矩形区域
  141. pDC->SelectObject(obr);
  142. br.DeleteObject(); 
  143. CRect dibrect;
  144. dibrect.top=0;
  145. dibrect.left=0;
  146. dibrect.right=sizeTotal.cx;
  147. dibrect.bottom=sizeTotal.cy;
  148. pDoc->m_Dib.Draw(pDC->GetSafeHdc(),dibrect,dibrect);
  149.     }
  150. /*绘制DDB：
  151. CBitmap ddb;
  152. ddb.LoadBitmap( IDB_BITMAP1 );
  153.     // Calculate bitmap size using a BITMAP structure.
  154.     BITMAP bm;
  155.     ddb.GetObject( sizeof(BITMAP), &bm );
  156.     // Create a memory DC, select the bitmap into the
  157.     // memory DC, and BitBlt it into the view.
  158.     CDC dcMem;
  159.     dcMem.CreateCompatibleDC( pDC );
  160.     CBitmap* pbmpOld = dcMem.SelectObject( &ddb );
  161.     pDC->BitBlt( dibrect.Width()/2,dibrect.Height()/2,
  162.  dibrect.Width()/2+bm.bmWidth,dibrect.Height()/2+bm.bmHeight,
  163.                  &dcMem,
  164.  0,0,
  165.  m_Rop 
  166.  ); 
  167.     // Reselect the original bitmap into the memory DC.
  168.     dcMem.SelectObject( pbmpOld ); */
  169. else
  170. {
  171.      CBrush br,* obr;
  172.     CRect rc;
  173. CString s;
  174. rc.top=0;
  175. rc.left=0;
  176. rc.right=sizeTotal.cx;
  177. rc.bottom=sizeTotal.cy;
  178.      GetClientRect(&rc);
  179.     br.CreateSolidBrush(RGB(255,255,255));
  180.     obr=pDC->SelectObject(&br);
  181.      pDC->FillRect(&rc,&br);
  182.     pDC->SelectObject(obr);
  183.     br.DeleteObject(); 
  184. if(pDoc->m_Dib.BPP>=16)
  185. {
  186. pen.CreatePen(0,2,RGB(0,0,255));
  187.      open=pDC->SelectObject(&pen);
  188. pDC->MoveTo(10,150);
  189. pDC->LineTo(520,150);
  190.      int i,j=10,k=0;
  191. for ( i=0;i<256;i=i+20)
  192. {
  193. s.Format("%d",i);
  194.                pDC->TextOut(10+i*2, 150,  s);
  195. }
  196. for ( i=0;i<256;i++)
  197. {
  198.      pDC->MoveTo(j,150);
  199.     if (m_Gradmax!=0) k=150-(int)m_Grade_Blue[i]*150/m_Gradmax;
  200.          pDC->LineTo(j,k);
  201.          j=j+2;
  202. }
  203. pen.DeleteObject();
  204. pen.CreatePen(0,2,RGB(0,255,0));
  205.      open=pDC->SelectObject(&pen);
  206. pDC->MoveTo(10,310);
  207. pDC->LineTo(520,310);
  208. j=10;
  209. for ( i=0;i<256;i=i+20)
  210. {
  211. s.Format("%d",i);
  212.                pDC->TextOut(10+i*2, 310, s);
  213. }
  214. for ( i=0;i<256;i++)
  215. {
  216.      pDC->MoveTo(j,310);
  217.     if (m_Gradmax!=0) k=310-(int)m_Grade_Green[i]*150/m_Gradmax;
  218.          pDC->LineTo(j,k);
  219.          j=j+2;
  220. }
  221. pen.DeleteObject();
  222. pen.CreatePen(0,2,RGB(255,0,0));
  223.      open=pDC->SelectObject(&pen);
  224. pDC->MoveTo(10,470);
  225. pDC->LineTo(520,470);
  226. for ( i=0;i<256;i=i+20)
  227. {
  228. s.Format("%d",i);
  229.                pDC->TextOut(10+i*2, 470, s);
  230. }
  231. j=10;
  232. for ( i=0;i<256;i++)
  233. {
  234.      pDC->MoveTo(j,470);
  235.     if (m_Gradmax!=0) k=470-(int)m_Grade_Red[i]*150/m_Gradmax;
  236.          pDC->LineTo(j,k);
  237.          j=j+2;
  238. }
  239. j=10;
  240. pDC->SelectObject(open);
  241. pen.DeleteObject();
  242. pDC->MoveTo(j,10);
  243. }
  244. else 
  245. {
  246.      pen.CreatePen(0,2,RGB(0,0,0));
  247.      open=pDC->SelectObject(&pen);
  248. pDC->MoveTo(10,260);
  249. pDC->LineTo(520,260);
  250. pDC->TextOut(200, 300, " 的直方图");
  251.      int i,j=10,k=0;
  252. for ( i=0;i<256;i=i+20)
  253. {
  254. s.Format("%d",i);
  255. pDC->TextOut(10+i*2, 260, s);
  256. }
  257. for ( i=0;i<256;i++)
  258. {
  259.      pDC->MoveTo(j,260);
  260.     if (m_Gradmax!=0) k=260-(int)m_Grade[i]*250/m_Gradmax;
  261.          pDC->LineTo(j,k);
  262.          j=j+2;
  263. }
  264. pDC->SelectObject(open);
  265. pen.DeleteObject();
  266. pDC->MoveTo(j,10);
  267. }
  268. }
  269. }
  271. void CBmpTestView::OnInitialUpdate()
  272. {
  273. CScrollView::OnInitialUpdate();
  275. CSize sizeTotal;
  276. // TODO: calculate the total size of this view
  277. sizeTotal.cx = sizeTotal.cy = 100;
  278. SetScrollSizes(MM_TEXT, sizeTotal);
  279. }
  281. /////////////////////////////////////////////////////////////////////////////
  282. // CBmpTestView printing
  284. BOOL CBmpTestView::OnPreparePrinting(CPrintInfo* pInfo)
  285. {
  286. // default preparation
  287. return DoPreparePrinting(pInfo);
  288. }
  290. void CBmpTestView::OnBeginPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
  291. {
  292. // TODO: add extra initialization before printing
  293. }
  295. void CBmpTestView::OnEndPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
  296. {
  297. // TODO: add cleanup after printing
  298. }
  300. /////////////////////////////////////////////////////////////////////////////
  301. // CBmpTestView diagnostics
  303. #ifdef _DEBUG
  304. void CBmpTestView::AssertValid() const
  305. {
  306. CScrollView::AssertValid();
  307. }
  309. void CBmpTestView::Dump(CDumpContext& dc) const
  310. {
  311. CScrollView::Dump(dc);
  312. }
  314. CBmpTestDoc* CBmpTestView::GetDocument() // non-debug version is inline
  315. {
  316. ASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CBmpTestDoc)));
  317. return (CBmpTestDoc*)m_pDocument;
  318. }
  319. #endif //_DEBUG
  321. /////////////////////////////////////////////////////////////////////////////
  322. // CBmpTestView message handlers
  324. void CBmpTestView::OnBtBlackness() 
  325. {
  326. // TODO: Add your command handler code here
  327. m_Rop=BLACKNESS;
  328. Invalidate();
  329. }
  331. void CBmpTestView::OnBtDstinvert() 
  332. {
  333. // TODO: Add your command handler code here
  334. m_Rop=DSTINVERT;
  335. Invalidate();
  336. }
  338. void CBmpTestView::OnBtMergecopy() 
  339. {
  340. // TODO: Add your command handler code here
  341. m_Rop=MERGECOPY;
  342. Invalidate();
  343. }
  345. void CBmpTestView::OnBtMergepaint() 
  346. {
  347. // TODO: Add your command handler code here
  348. m_Rop=MERGEPAINT;
  349. Invalidate();
  350. }
  352. void CBmpTestView::OnBtNotsrccopy() 
  353. {
  354. // TODO: Add your command handler code here
  355. m_Rop=NOTSRCCOPY;
  356. Invalidate();
  357. }
  359. void CBmpTestView::OnBtNotsrcerase() 
  360. {
  361. // TODO: Add your command handler code here
  362. m_Rop=NOTSRCERASE;
  363. Invalidate();
  364. }
  366. void CBmpTestView::OnBtPatcopy() 
  367. {
  368. // TODO: Add your command handler code here
  369. m_Rop=PATCOPY;
  370. Invalidate();
  371. }
  373. void CBmpTestView::OnBtPatinvert() 
  374. {
  375. // TODO: Add your command handler code here
  376. m_Rop=PATINVERT;
  377. Invalidate();
  378. }
  380. void CBmpTestView::OnBtPatpaint() 
  381. {
  382. // TODO: Add your command handler code here
  383. m_Rop=PATPAINT;
  384. Invalidate();
  385. }
  387. void CBmpTestView::OnBtSrcand() 
  388. {
  389. // TODO: Add your command handler code here
  390. m_Rop=SRCAND;
  391. Invalidate();
  392. }
  394. void CBmpTestView::OnBtSrccopy() 
  395. {
  396. // TODO: Add your command handler code here
  397. m_Rop=SRCCOPY;
  398. Invalidate();
  399. }
  401. void CBmpTestView::OnBtSrcerase() 
  402. {
  403. // TODO: Add your command handler code here
  404. m_Rop=SRCERASE;
  405. Invalidate();
  406. }
  408. void CBmpTestView::OnBtSrcinvert() 
  409. {
  410. // TODO: Add your command handler code here
  411. m_Rop=SRCINVERT;
  412. Invalidate();
  413. }
  415. void CBmpTestView::OnBtSrcpaint() 
  416. {
  417. // TODO: Add your command handler code here
  418. m_Rop=SRCPAINT;
  419. Invalidate();
  420. }
  422. void CBmpTestView::OnBtWhiteness() 
  423. {
  424. // TODO: Add your command handler code here
  425. m_Rop=WHITENESS;
  426. Invalidate();
  427. }
  429. void CBmpTestView::OnZftShow(CDC *pDC) 
  430. {
  431. // TODO: Add your command handler code here
  432. /* CRect rc;
  433. GetClientRect(&rc); */
  434. CBmpTestDoc* pDoc = GetDocument();
  435. ASSERT_VALID(pDoc);  
  436. if(!pDoc->m_Dib.m_BmpLoaded)
  437. {
  438. AfxMessageBox("还没有打开图像文件");
  439. return;
  440. }
  441.     int i;
  442.     long j;
  443. if (pDoc->m_Dib.BPP>=16)
  444.     {   
  445.    m_Gradmax=0;
  446.    for ( i=0;i<256;i++)
  447.    {
  448.       m_Grade_Blue[i]=0;
  449.   m_Grade_Green[i]=0;
  450.   m_Grade_Red[i]=0;
  451.    }
  452.        for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j=j+3)
  453.    {
  454.       m_Grade_Blue[int(*(pDoc->m_Dib.m_Buffer+j))]++;
  455.   m_Grade_Green[int(*(pDoc->m_Dib.m_Buffer+j+1))]++;
  456.   m_Grade_Red[int(*(pDoc->m_Dib.m_Buffer+j+2))]++;
  457.    }
  458.        for ( i=0;i<256;i++)
  459.    {
  460.      if (m_Grade_Blue[i]>m_Gradmax)
  461. m_Gradmax=m_Grade_Blue[i];
  462.  if (m_Grade_Green[i]>m_Gradmax)
  463. m_Gradmax=m_Grade_Green[i];
  464.  if (m_Grade_Red[i]>m_Gradmax)
  465. m_Gradmax=m_Grade_Red[i];
  466.    }
  467. }
  468.     else if (pDoc->m_Dib.BPP==8)
  469. {
  470.    m_Gradmax=0;
  471.        for ( i=0;i<256;i++)
  472.    {
  473.       m_Grade[i]=0;
  474.    }
  475.        for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  476.    {
  477.       m_Grade[int(*(pDoc->m_Dib.m_Buffer+j))]++;
  478.    }
  479.    m_Gradmax=0;
  480.        for ( i=0;i<256;i++)
  481.    {
  482.      if (m_Grade[i]>m_Gradmax)
  483. m_Gradmax=m_Grade[i];
  484.    }
  485. }
  486. else
  487. {
  488.    long int black=0, white=0;
  489.    int a;
  490.    for ( i=0;i<256;i++)
  491.    {
  492.       m_Grade[i]=0;
  493.    }
  494.        for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  495.    {
  496.    a=128;
  497.        for ( i=0;i<8;i++)
  498.    {
  499.     if( ((int)*(pDoc->m_Dib.m_Buffer+j) & a)==0)
  500. black++;
  501.     else
  502. white++;
  503.     a=a>>1;
  504.    }
  505.    }
  506.    if (black>white)
  507. m_Gradmax=black;
  508.    else
  509.    m_Gradmax=white;
  510.        m_Grade[0]=black;
  511.    m_Grade[255]=white;
  512. }
  513. pDoc->m_BM=false;
  514. Invalidate();
  515. }
  517. void CBmpTestView::OnZftBalance() 
  518. {
  519. // TODO: Add your command handler code here
  521. CBmpTestDoc* pDoc = GetDocument();
  522. ASSERT_VALID(pDoc);  
  523.     int i;
  524. double m_GradAdd[256];
  525. long j;
  526.     m_GradAdd[0]=m_Grade[0];
  527. // calculate integral
  528. for ( i=1;i<256;i++)
  529. {
  530.     m_GradAdd[i]=0;
  531. }
  532. for ( i=1;i<256;i++)
  533. {
  534.     m_GradAdd[i]=m_GradAdd[i-1]+m_Grade[i];
  535. }
  536. for ( i=1;i<256;i++)
  537. {
  538.     m_GradAdd[i]=m_GradAdd[i]*255/m_GradAdd[255];
  539. }
  540. // 归一化阈值面积function
  541.     for ( j=0;j<(long int) pDoc->m_Dib.BufferSize;j++)
  542. {
  543.     *(pDoc->m_Dib.m_Buffer+j)=(int) (m_GradAdd[(*(pDoc->m_Dib.m_Buffer+j))]*255/m_GradAdd[255]);
  544. }
  545. for ( i=0;i<256;i++)
  546. {
  547.     m_Grade[i]=0;
  548. }
  549.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  550. {
  551.     m_Grade[int(*(pDoc->m_Dib.m_Buffer+j))]++;
  552. }
  553. for ( i=0;i<256;i++)
  554. {
  555.     if (m_Grade[i]>m_Gradmax)
  556. m_Gradmax=m_Grade[i];
  557. }
  558. pDoc->m_BM=false;
  559. Invalidate();
  560. }
  562. void CBmpTestView::OnImageShow() 
  563. {
  564. // TODO: Add your command handler code here
  565. CBmpTestDoc* pDoc = GetDocument();
  566. ASSERT_VALID(pDoc);  
  567. pDoc->m_BM=true;
  568. Invalidate();
  569. }
  571. void CBmpTestView::OnDtsmIn() 
  572. {
  573. // TODO: Add your command handler code here
  574. CBmpTestDoc* pDoc = GetDocument();
  575. ASSERT_VALID(pDoc);  
  576. if(pDoc->m_Dib.BPP!=8)
  577. {
  578.        AfxMessageBox( "此图像为非256色图像，目前该算子还不支持这种图像格式");
  579.    return;
  580. }
  581. float T0,u1,u2,T1,w1,w2,u;
  582. int i,decimal, sign;
  583. long j;
  584. char string[80]="二值化最佳阈值为",*string1,*string2;
  585.   for ( i=0;i<256;i++)
  586. {
  587.     m_Grade[i]=0;
  588. }
  589.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  590. {
  591.     m_Grade[int(*(pDoc->m_Dib.m_Buffer+j))]++;
  592. }
  593. u=0;
  594. for(i=0;i<256;i++)
  595. u=u+(float) i*m_Grade[i]/(int) pDoc->m_Dib.BufferSize;
  596. T0=u;T1=0;
  597.     while (fabs(T0-T1)>0.01)
  598. {
  599. w1=0;w2=0;u1=0;u2=0;
  600. for(i=0;i<=int(T0);i++)
  601. {
  602. w1=w1+m_Grade[i];
  603. u1=u1+i*m_Grade[i];
  604. }
  605. for(i=int(T0)+1;i<256;i++)
  606. {
  607. w2=w2+m_Grade[i];
  608. u2=u2+i*m_Grade[i];
  609. }
  610. T1=T0;u1=u1/w1;u2=u2/w2;T0=0.5*(u1+u2);
  611. }
  612. string2=_fcvt(T0, 0, &decimal, &sign );
  613. string1=strcat(string,string2);
  614.     AfxMessageBox(string1);
  615.     for ( j=0;j<(long int) pDoc->m_Dib.BufferSize;j++)
  616. {
  617. *(pDoc->m_Dib.m_Buffer+j)=(*(pDoc->m_Dib.m_Buffer+j)>=T0)?255:0;
  618. }
  619. }
  622. void CBmpTestView::OnDtsmShow() 
  623. {
  624. // TODO: Add your command handler code here
  625. CBmpTestView::OnImageShow();
  626. }
  628. void CBmpTestView::OnDtsmOut() 
  629. {
  630. // TODO: Add your command handler code here
  631. CBmpTestDoc* pDoc = GetDocument();
  632. ASSERT_VALID(pDoc);  
  633. float u,uT[255],w[255],sigma[255],sigmamax;
  634. int i,T,Topt;
  635. long j;
  636. char string[80]="类间方差法二值化最佳阈值为",*string1,*string2,buffer[10];
  637.   for ( i=0;i<256;i++)
  638. {
  639.     m_Grade[i]=0;
  640. }
  641.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  642. {
  643.     m_Grade[int(*(pDoc->m_Dib.m_Buffer+j))]++;
  644. }
  645. u=0;
  646. for(i=0;i<256;i++)
  647. {
  648. u=u+i*m_Grade[i];
  649. }
  650. u=u/pDoc->m_Dib.BufferSize;
  651. sigmamax=0;Topt=0;
  652. for(T=0;T<256;T++)
  653. {
  654. w[T]=0;uT[T]=0;
  655. for(i=0;i<=T;i++)
  656. {
  657. w[T]=w[T]+m_Grade[i];
  658. uT[T]=uT[T]+i*m_Grade[i];
  659. }
  660. uT[T]=uT[T]/pDoc->m_Dib.BufferSize;
  661. w[T]=w[T]/pDoc->m_Dib.BufferSize;
  662. sigma[T]=(u*w[T]-uT[T])*(u*w[T]-uT[T])/w[T]*(1-w[T]);
  663. if(sigma[T]>sigmamax)
  664. {
  665. sigmamax=sigma[T];Topt=T;
  666. }
  667. }
  668. string2=_itoa(Topt, buffer,10);
  669. string1=strcat(string,string2);
  670.     AfxMessageBox(string1);
  671. }
  673. void CBmpTestView::OnRecFilt() 
  674. {
  675. // TODO: Add your command handler code here
  676. CBmpTestDoc* pDoc = GetDocument();
  677. ASSERT_VALID(pDoc);  
  678. int i,j,image[4][4],Height,Width;
  679. //filtering
  680. Height=(int) pDoc->m_Dib.Height;
  681. Width=(int) pDoc->m_Dib.Bpl;
  682.   for (i=1;i<Height-1;i++)
  683. {
  684.     for (j=1;j<Width +1;j++)
  685. {
  686. image[2][2]=(int) *(pDoc->m_Dib.m_Buffer1+i*Width+j);
  687. image[2][1]=(int) *(pDoc->m_Dib.m_Buffer1+i*Width+j-1);
  688. image[2][3]=(int) *(pDoc->m_Dib.m_Buffer1+i*Width+j+1);
  689. image[1][2]=(int) *(pDoc->m_Dib.m_Buffer1+(i-1)*Width+j);
  690. image[1][1]=(int) *(pDoc->m_Dib.m_Buffer1+(i-1)*Width+j-1);
  691. image[1][3]=(int) *(pDoc->m_Dib.m_Buffer1+(i-1)*Width+j+1);
  692. image[3][2]=(int) *(pDoc->m_Dib.m_Buffer1+(i+1)*Width+j);
  693. image[3][1]=(int) *(pDoc->m_Dib.m_Buffer1+(i+1)*Width+j-1);
  694. image[3][3]=(int) *(pDoc->m_Dib.m_Buffer1+(i+1)*Width+j+1);
  695. *(pDoc->m_Dib.m_Buffer+i*Width+j)=(image[1][1]+image[1][2]+image[1][3]+image[2][1]+image[2][2]+image[2][3]+image[3][1]+image[3][2]+image[3][3])/9;
  696. }
  697. }
  698. pDoc->m_BM=true;
  699. Invalidate();
  700. }
  702. /*void queue(int *im,int n)
  703. //int im[],n;
  704. {
  705. int i,j,medium;
  706. for(i=1;i<n;i++)
  707. {
  708. for(j=i+1;j<n+1;j++)
  709. {
  710. if(im(j)>im(i))
  711. {
  712. medium=im(i);
  713. im(i)=im(j);
  714. im(j)=medium;
  715. }
  716. }
  717. }
  718. }*/
  720. void CBmpTestView::OnMidFilt() 
  721. {
  722. // TODO: Add your command handler code here
  723. CBmpTestDoc* pDoc = GetDocument();
  724. ASSERT_VALID(pDoc);  
  725. BYTE * m_Buffer1;
  726. int i,Height,Width,i1,j1,j;
  727. //filtering
  728. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  729.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  730. {
  731.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  732. }
  733. if(pDoc->m_Dib.BPP==8)
  734. {
  735. int image[10],medium;
  736. Height=(int) pDoc->m_Dib.Height;
  737. Width=(int) pDoc->m_Dib.Bpl;
  738.   for (i=1;i<Height-1;i++)
  739. {
  740. for (j=1;j<Width +1;j++)
  741. {
  742. image[5]=(int) *(m_Buffer1+i*Width+j);
  743. image[4]=(int) *(m_Buffer1+i*Width+j-1);
  744. image[6]=(int) *(m_Buffer1+i*Width+j+1);
  745. image[2]=(int) *(m_Buffer1+(i-1)*Width+j);
  746. image[1]=(int) *(m_Buffer1+(i-1)*Width+j-1);
  747. image[3]=(int) *(m_Buffer1+(i-1)*Width+j+1);
  748. image[8]=(int) *(m_Buffer1+(i+1)*Width+j);
  749. image[7]=(int) *(m_Buffer1+(i+1)*Width+j-1);
  750. image[9]=(int) *(m_Buffer1+(i+1)*Width+j+1);
  751.          for(i1=1;i1<9;i1++)
  752. {
  753.          for(j1=i1+1;j1<10;j1++)
  754. {
  755.      if(image[j1]>image[i1])
  756. {
  757.      medium=image[i1];
  758.          image[i1]=image[j1];
  759.          image[j1]=medium;
  760. }
  761. }
  762. }
  763. *(pDoc->m_Dib.m_Buffer+i*Width+j)=image[5];
  764. }
  765. }
  766. }
  767. else
  768. {
  769. // 24位彩色中值滤波
  770. int image[10][4],medium[4],k;
  771. Height=(int) pDoc->m_Dib.Height;
  772. Width=(int) pDoc->m_Dib.Bpl;
  773.   for (i=1;i<Height-1;i++)
  774. {
  775. for (j=1;j<(int)Width/3 ;j++)
  776. {
  777. for(k=0;k<3;k++)
  778. {
  779. image[5][k+1]=(int) *(m_Buffer1+i*Width+j*3+k);
  780. image[4][k+1]=(int) *(m_Buffer1+i*Width+(j-1)*3+k);
  781. image[6][k+1]=(int) *(m_Buffer1+i*Width+(j+1)*3+k);
  782. image[2][k+1]=(int) *(m_Buffer1+(i-1)*Width+j*3+k);
  783. image[1][k+1]=(int) *(m_Buffer1+(i-1)*Width+(j-1)*3+k);
  784. image[3][k+1]=(int) *(m_Buffer1+(i-1)*Width+(j+1)*3+k);
  785. image[8][k+1]=(int) *(m_Buffer1+(i+1)*Width+j*3+k);
  786. image[7][k+1]=(int) *(m_Buffer1+(i+1)*Width+(j-1)*3+k);
  787. image[9][k+1]=(int) *(m_Buffer1+(i+1)*Width+(j+1)*3+k);
  788. }
  789. //冒泡排序
  790. for(i1=1;i1<9;i1++)
  791. {
  792.          for(j1=i1+1;j1<10;j1++)
  793. {
  794.      if((image[j1][1]+image[j1][2]+image[j1][3])/3>(image[i1][1]+image[i1][2]+image[i1][3])/3)
  795. {
  796. for(k=1;k<4;k++)
  797. {
  798.      medium[k]=image[i1][k];
  799. image[i1][k]=image[j1][k];
  800. image[j1][k]=medium[k];
  801. }
  802. }
  803. }
  804. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=image[5][1];
  805. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=image[5][2];
  806. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=image[5][3];
  807. }
  808. }
  809. }
  810. }
  811. if(m_Buffer1)
  812. {
  813. delete [] m_Buffer1;
  814. m_Buffer1=0;
  815. }
  816. pDoc->m_BM=true;
  817. Invalidate();
  819. }
  821. void CBmpTestView::OnEdgeRobert() 
  822. {
  823. // TODO: Add your command handler code here
  824. #include  <math.h>
  825. CBmpTestDoc* pDoc = GetDocument();
  826. ASSERT_VALID(pDoc);
  827. if(pDoc->m_Dib.BPP>8)
  828. {
  829. AfxMessageBox("该算子不能处理8位图像");
  830. return;
  831. }
  832. BYTE * m_Buffer1;
  833. int i,j,image[10],Height,Width,a,b;
  834. //filtering
  835. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  836.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  837. {
  838.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  839. }
  840. Height=(int) pDoc->m_Dib.Height;
  841. Width=(int) pDoc->m_Dib.Bpl;
  842.   for (i=1;i<Height-1;i++)
  843. {
  844.     for (j=1;j<Width -1;j++)
  845. {
  846. image[1]=(unsigned int) *(m_Buffer1+(i-1)*Width+j-1);
  847. image[3]=(unsigned int) *(m_Buffer1+(i-1)*Width+j+1);
  848. image[7]=(unsigned int) *(m_Buffer1+(i+1)*Width+j-1);
  849. image[9]=(unsigned int) *(m_Buffer1+(i+1)*Width+j+1);
  850.             a=abs(image[1]-image[9]);
  851. b=abs(image[3]-image[7]);
  852. *(pDoc->m_Dib.m_Buffer+i*Width+j)=((a>b)?a:b);
  853. }
  854. }
  855. if(m_Buffer1)
  856. {
  857. delete [] m_Buffer1;
  858. m_Buffer1=0;
  859. }
  860. pDoc->m_BM=true;
  861. Invalidate();
  862. }
  864. void CBmpTestView::OnEdgeLaplace() 
  865. {
  866. // TODO: Add your command handler code here
  867. CBmpTestDoc* pDoc = GetDocument();
  868. ASSERT_VALID(pDoc);  
  869. if(pDoc->m_Dib.BPP>8)
  870. {
  871. AfxMessageBox("该算子只能处理8位图像");
  872. return;
  873. }
  874. int i,j,image[10],Height,Width,a;
  875. //filtering
  876. Height=(int) pDoc->m_Dib.Height;
  877. Width=(int) pDoc->m_Dib.Bpl;
  878. // Total=(int) pDoc->m_Dib.BufferSize;
  879.   for (i=1;i<Height-1;i++)
  880. {
  881.     for (j=1;j<Width +1;j++)
  882. {
  883. image[5]=(unsigned int) *(pDoc->m_Dib.m_Buffer1+i*Width+j);
  884. image[4]=(unsigned int) *(pDoc->m_Dib.m_Buffer1+i*Width+j-1);
  885. image[6]=(unsigned int) *(pDoc->m_Dib.m_Buffer1+i*Width+j+1);
  886. image[2]=(unsigned int) *(pDoc->m_Dib.m_Buffer1+(i-1)*Width+j);
  887. image[1]=(unsigned int) *(pDoc->m_Dib.m_Buffer1+(i-1)*Width+j-1);
  888. image[3]=(unsigned int) *(pDoc->m_Dib.m_Buffer1+(i-1)*Width+j+1);
  889. image[8]=(unsigned int) *(pDoc->m_Dib.m_Buffer1+(i+1)*Width+j);
  890. image[7]=(unsigned int) *(pDoc->m_Dib.m_Buffer1+(i+1)*Width+j-1);
  891. image[9]=(unsigned int) *(pDoc->m_Dib.m_Buffer1+(i+1)*Width+j+1);
  892. a=image[2]+image[4]+image[6]+image[8]-4*image[5];
  893. *(pDoc->m_Dib.m_Buffer+i*Width+j)=((a>0)?a:0);
  894. }
  895. }
  896. pDoc->m_BM=true;
  897. Invalidate();
  898. }
  900. void CBmpTestView::OnGauseFilt() 
  901. {
  902. // TODO: Add your command handler code here
  904. CBmpTestDoc* pDoc = GetDocument();
  905. ASSERT_VALID(pDoc);  
  906. if(pDoc->m_Dib.BPP>8)
  907. {
  908. AfxMessageBox("该算子只能处理8位图像");
  909. return;
  910. }
  911. int i,j,image[4][4],Height,Width;
  912. double a0,a1,a2;
  913. a0=1/sqrt(2*3.14);
  914. a1=0.6065/sqrt(2*3.14);
  915. a2=exp(-1)/sqrt(2*3.14);
  916. //filtering
  917. Height=(int) pDoc->m_Dib.Height;
  918. Width=(int) pDoc->m_Dib.Bpl;
  919.   for (i=1;i<Height-1;i++)
  920. {
  921.     for (j=1;j<Width +1;j++)
  922. {
  923. image[2][2]=(int) *(pDoc->m_Dib.m_Buffer1+i*Width+j);
  924. image[2][1]=(int) *(pDoc->m_Dib.m_Buffer1+i*Width+j-1);
  925. image[2][3]=(int) *(pDoc->m_Dib.m_Buffer1+i*Width+j+1);
  926. image[1][2]=(int) *(pDoc->m_Dib.m_Buffer1+(i-1)*Width+j);
  927. image[1][1]=(int) *(pDoc->m_Dib.m_Buffer1+(i-1)*Width+j-1);
  928. image[1][3]=(int) *(pDoc->m_Dib.m_Buffer1+(i-1)*Width+j+1);
  929. image[3][2]=(int) *(pDoc->m_Dib.m_Buffer1+(i+1)*Width+j);
  930. image[3][1]=(int) *(pDoc->m_Dib.m_Buffer1+(i+1)*Width+j-1);
  931. image[3][3]=(int) *(pDoc->m_Dib.m_Buffer1+(i+1)*Width+j+1);
  932. *(pDoc->m_Dib.m_Buffer+i*Width+j)=(int)(image[1][1]*a2+image[1][2]*a1+image[1][3]*a2+image[2][1]*a1+image[2][2]*a0+image[2][3]*a1+image[3][1]*a2+image[3][2]*a1+image[3][3]*a2);
  933. }
  934. }
  935. pDoc->m_BM=true;
  936. Invalidate();
  938. }
  939. void CBmpTestView::OnEdgeMarr() 
  940. {
  941. // TODO: Add your command handler code here
  942. OnGauseFilt();
  943. OnEdgeLaplace();
  944. }
  947. void CBmpTestView::OnStruAdd() 
  948. {
  949. // 二值形态学：膨胀，结构和
  950. CBmpTestDoc* pDoc = GetDocument();
  951. ASSERT_VALID(pDoc);  
  952. int i,j,image[4][4],Height,Width,a,k;
  953. //filtering
  954. Height=(int) pDoc->m_Dib.Height;
  955. Width=(int) pDoc->m_Dib.Bpl;
  956. BYTE * m_Buffer1;
  957. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  958.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  959. {
  960.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  961. }
  962. m_Gradmax=0;
  963.     for ( i=0;i<256;i++)
  964. {
  965.     m_Grade[i]=0;
  966. }
  967.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  968. {
  969.     m_Grade[int(*(pDoc->m_Dib.m_Buffer+j))]++;
  970. }
  971. m_Gradmax=0;
  972.     for ( i=1;i<256;i++)
  973. {
  974.    if (m_Grade[i]>m_Gradmax)
  975.    {
  976.   m_Gradmax=m_Grade[i];
  977.   k=i;
  978.    }
  979. }
  980.   for (i=1;i<Height-1;i++)
  981. {
  982.     for (j=1;j<Width +1;j++)
  983. {
  984. image[2][2]=(unsigned int) *(m_Buffer1+i*Width+j);
  985. image[2][1]=(unsigned int) *(m_Buffer1+i*Width+j-1);
  986. image[2][3]=(unsigned int) *(m_Buffer1+i*Width+j+1);
  987. image[1][2]=(unsigned int) *(m_Buffer1+(i-1)*Width+j);
  988. image[1][1]=(unsigned int) *(m_Buffer1+(i-1)*Width+j-1);
  989. image[1][3]=(unsigned int) *(m_Buffer1+(i-1)*Width+j+1);
  990. image[3][2]=(unsigned int) *(m_Buffer1+(i+1)*Width+j);
  991. image[3][1]=(unsigned int) *(m_Buffer1+(i+1)*Width+j-1);
  992. image[3][3]=(unsigned int) *(m_Buffer1+(i+1)*Width+j+1);
  993. a=image[1][1]||image[1][2]||image[1][3]||image[2][1]||image[2][2]||image[2][3]||image[3][1]||image[3][2]||image[3][3];
  994. *(pDoc->m_Dib.m_Buffer+i*Width+j)=(a?k:0);
  995. }
  996. }
  997. if(m_Buffer1)
  998. {
  999. delete [] m_Buffer1;
  1000. m_Buffer1=0;
  1001. }
  1002. pDoc->m_BM=true;
  1003. Invalidate();
  1004. }
  1006. void CBmpTestView::OnStruSub() 
  1007. {
  1008. // TODO: Add your command handler code here
  1009. CBmpTestDoc* pDoc = GetDocument();
  1010. ASSERT_VALID(pDoc);  
  1011. int i,j,image[4][4],Height,Width,a,k;
  1012. //filtering
  1013. Height=(int) pDoc->m_Dib.Height;
  1014. Width=(int) pDoc->m_Dib.Bpl;
  1015. //如下代码可防止不同图像调色板的不同
  1016. BYTE * m_Buffer1;
  1017. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  1018.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  1019. {
  1020.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  1021. }
  1022. m_Gradmax=0;
  1023.     for ( i=0;i<256;i++)
  1024. {
  1025.     m_Grade[i]=0;
  1026. }
  1027.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  1028. {
  1029.     m_Grade[int(*(pDoc->m_Dib.m_Buffer+j))]++;
  1030. }
  1031. m_Gradmax=0;
  1032.     for ( i=1;i<256;i++)
  1033. {
  1034.    if (m_Grade[i]>m_Gradmax)
  1035.    {
  1036.   m_Gradmax=m_Grade[i];
  1037.   k=i;
  1038.    }
  1039. }
  1040.   for (i=1;i<Height-1;i++)
  1041. {
  1042.     for (j=1;j<Width +1;j++)
  1043. {
  1044. image[2][2]=(int) *(m_Buffer1+i*Width+j);
  1045. image[2][1]=(int) *(m_Buffer1+i*Width+j-1);
  1046. image[2][3]=(int) *(m_Buffer1+i*Width+j+1);
  1047. image[1][2]=(int) *(m_Buffer1+(i-1)*Width+j);
  1048. image[1][1]=(int) *(m_Buffer1+(i-1)*Width+j-1);
  1049. image[1][3]=(int) *(m_Buffer1+(i-1)*Width+j+1);
  1050. image[3][2]=(int) *(m_Buffer1+(i+1)*Width+j);
  1051. image[3][1]=(int) *(m_Buffer1+(i+1)*Width+j-1);
  1052. image[3][3]=(int) *(m_Buffer1+(i+1)*Width+j+1);
  1053. a=image[1][1]&&image[1][2]&&image[1][3]&&image[2][1]&&image[2][2]&&image[2][3]&&image[3][1]&&image[3][2]&&image[3][3];
  1054. *(pDoc->m_Dib.m_Buffer+i*Width+j)=(a?k:0);
  1055. }
  1056. }
  1057. if(m_Buffer1)
  1058. {
  1059. delete [] m_Buffer1;
  1060. m_Buffer1=0;
  1061. }
  1062. pDoc->m_BM=true;
  1063. Invalidate();
  1064. }
  1066. void CBmpTestView::OnStruOpen() 
  1067. {
  1068. // TODO: Add your command handler code here
  1069. OnStruSub();
  1070. OnStruAdd();
  1071. }
  1073. void CBmpTestView::OnStruClose() 
  1074. {
  1075. // TODO: Add your command handler code here
  1076. OnStruAdd();
  1077. OnStruSub();
  1078. }
  1080. void CBmpTestView::OnEditBackup() 
  1081. {
  1082. // TODO: Add your command handler code here
  1083. CBmpTestDoc* pDoc = GetDocument();
  1084. ASSERT_VALID(pDoc);  
  1085.     pDoc->m_Dib.restoreimage();
  1086. pDoc->m_BM=true;
  1087. Invalidate();
  1088. }
  1090. void CBmpTestView::OnEditBlue() 
  1091. {
  1092. // TODO: Add your command handler code here
  1093. CBmpTestDoc* pDoc = GetDocument();
  1094. ASSERT_VALID(pDoc);
  1095. if(pDoc->m_Dib.BPP<=8)
  1096. {
  1097.        AfxMessageBox( "此图像为256色图像，目前暂时不能分出独立通道");
  1098.    return;
  1099. }
  1100. if(pDoc->m_Dib.GetBlueImg())
  1101. Invalidate();
  1102. }
  1104. void CBmpTestView::OnEditGreen() 
  1105. {
  1106. // TODO: Add your command handler code here
  1107. CBmpTestDoc* pDoc = GetDocument();
  1108. ASSERT_VALID(pDoc);  
  1109. if(pDoc->m_Dib.BPP<=8)
  1110. {
  1111.        AfxMessageBox( "此图像为256色图像，目前暂时不能分出独立通道");
  1112.    return;
  1113. }
  1114. if(pDoc->m_Dib.GetGreenImg())
  1115. Invalidate();
  1116. }
  1118. void CBmpTestView::OnEditRed() 
  1119. {
  1120. // TODO: Add your command handler code here
  1121. CBmpTestDoc* pDoc = GetDocument();
  1122. ASSERT_VALID(pDoc);  
  1123. if(pDoc->m_Dib.BPP<=8)
  1124. {
  1125.        AfxMessageBox( "此图像为256色图像，目前暂时不能分出独立通道");
  1126.    return;
  1127. }
  1128. if(pDoc->m_Dib.GetRedImg())
  1129. Invalidate();
  1130. }
  1132. void CBmpTestView::OnEdgeColorprewitt() 
  1133. {
  1134. // TODO: 一种强有力的非线性算子，加强了纵向，但也造成横向断裂
  1135. # include <math.h>
  1136. CBmpTestDoc* pDoc = GetDocument();
  1137. ASSERT_VALID(pDoc);  
  1138. if(pDoc->m_Dib.BPP<=8)
  1139. {
  1140.        AfxMessageBox( "此图像为256色图像，目前该算子还不支持这种图像格式");
  1141.    return;
  1142. }
  1143. int i,j,k,image[10][4],Height,Width;
  1144. unsigned int a;
  1145. BYTE * m_Buffer1;
  1146. double e1,e2,e3;
  1147. //filtering
  1148. Height=(int) pDoc->m_Dib.Height;
  1149. Width=(int) pDoc->m_Dib.Bpl;
  1150. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  1151.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  1152. {
  1153.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  1154. }
  1155.   for (i=1;i<Height-1;i++)
  1156. {
  1157.     for (j=1;j<(int)Width/3 -1;j++)
  1158. {
  1159. for (k=0;k<3;k++)
  1160. {
  1161. image[1][k+1]=(unsigned int) *(m_Buffer1+(i-1)*Width+(j-1)*3+k);
  1162. image[3][k+1]=(unsigned int) *(m_Buffer1+(i-1)*Width+(j+1)*3+k);
  1163. image[4][k+1]=(unsigned int) *(m_Buffer1+i*Width+(j-1)*3+k);
  1164. image[6][k+1]=(unsigned int) *(m_Buffer1+i*Width+(j+1)*3+k);
  1165. image[7][k+1]=(unsigned int) *(m_Buffer1+(i+1)*Width+(j-1)*3+k);
  1166. image[9][k+1]=(unsigned int) *(m_Buffer1+(i+1)*Width+(j+1)*3+k);
  1167. }
  1168.             e1=sqrt(pow((image[1][1]-image[3][1]),2)+pow((image[1][2]-image[3][2]),2)+pow((image[1][3]-image[3][3]),2));
  1169. e2=sqrt(pow((image[4][1]-image[6][1]),2)+pow((image[4][2]-image[6][2]),2)+pow((image[4][3]-image[6][3]),2));
  1170. e3=sqrt(pow((image[7][1]-image[9][1]),2)+pow((image[7][2]-image[9][2]),2)+pow((image[7][3]-image[9][3]),2));
  1171. a=(unsigned int)(e1+e2+e3);
  1172. a=((a>255)?255:a);
  1173. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=a;
  1174. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=a;
  1175. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=a;
  1176. }
  1177. }
  1178. if(m_Buffer1)
  1179. {
  1180. delete [] m_Buffer1;
  1181. m_Buffer1=0;
  1182. }
  1183. pDoc->m_BM=true;
  1184. Invalidate();
  1185. }
  1187. void CBmpTestView::OnEdgeColorrobert() 
  1188. {
  1189. // TODO: 一种强烈的非线性算子，加强了双45度方向，保持了字符笔划的完整性
  1190. # include <math.h>
  1191. CBmpTestDoc* pDoc = GetDocument();
  1192. ASSERT_VALID(pDoc);  
  1193. if(pDoc->m_Dib.BPP<=8)
  1194. {
  1195.        AfxMessageBox( "此图像为256色图像，目前该算子还不支持这种图像格式");
  1196.    return;
  1197. }
  1198. int i,j,k,image[10][4],Height,Width;
  1199. unsigned int a;
  1200. BYTE * m_Buffer1;
  1201. double e2,e1;
  1202. //filtering
  1203. Height=(int) pDoc->m_Dib.Height;
  1204. Width=(int) pDoc->m_Dib.Bpl;
  1205. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  1206.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  1207. {
  1208.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  1209. }
  1210.   for (i=1;i<Height-1;i++)
  1211. {
  1212.     for (j=1;j<(int)Width/3 -1;j++)
  1213. {
  1214. for (k=0;k<3;k++)
  1215. {
  1216. image[1][k+1]=(unsigned int) *(m_Buffer1+(i-1)*Width+(j-1)*3+k);
  1217. image[3][k+1]=(unsigned int) *(m_Buffer1+(i-1)*Width+(j+1)*3+k);
  1218. image[7][k+1]=(unsigned int) *(m_Buffer1+(i+1)*Width+(j-1)*3+k);
  1219. image[9][k+1]=(unsigned int) *(m_Buffer1+(i+1)*Width+(j+1)*3+k);
  1220. }
  1221.             e1=pow((image[1][1]-image[9][1]),2)+pow((image[1][2]-image[9][2]),2)+pow((image[1][3]-image[9][3]),2);
  1222. e2=pow((image[3][1]-image[7][1]),2)+pow((image[3][2]-image[7][2]),2)+pow((image[3][3]-image[7][3]),2);
  1223. a=(unsigned int)sqrt(e1+e2);
  1224. a=((a>255)?255:a);
  1225. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=a;
  1226. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=a;
  1227. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=a;
  1228. }
  1229. }
  1230. if(m_Buffer1)
  1231. {
  1232. delete [] m_Buffer1;
  1233. m_Buffer1=0;
  1234. }
  1235. pDoc->m_BM=true;
  1236. Invalidate();
  1238. }
  1240. void CBmpTestView::OnMidFiltLine() 
  1241. {
  1242. // TODO: Add your command handler code here
  1243. CBmpTestDoc* pDoc = GetDocument();
  1244. ASSERT_VALID(pDoc);  
  1245. BYTE * m_Buffer1;
  1246. int i,Height,Width,i1,j1,j;
  1247. //filtering
  1248. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  1249.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  1250. {
  1251.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  1252. }
  1253. if(pDoc->m_Dib.BPP==8)
  1254. {
  1255. int image[6],medium;
  1256. Height=(int) pDoc->m_Dib.Height;
  1257. Width=(int) pDoc->m_Dib.Bpl;
  1258. for (j=0;j<Width +1;j++)
  1259. {
  1260.   for (i=2;i<Height-2;i++)
  1261. {
  1262. image[3]=(int) *(m_Buffer1+i*Width+j);
  1263. image[2]=(int) *(m_Buffer1+(i-1)*Width+j);
  1264. image[4]=(int) *(m_Buffer1+(i+1)*Width+j);
  1265. image[1]=(int) *(m_Buffer1+(i-2)*Width+j);
  1266. image[5]=(int) *(m_Buffer1+(i+2)*Width+j);
  1267. //冒泡排序
  1268. for(i1=1;i1<5;i1++)
  1269. {
  1270.          for(j1=i1+1;j1<6;j1++)
  1271. {
  1272.      if(image[j1]>image[i1])
  1273. {
  1274.      medium=image[i1];
  1275.          image[i1]=image[j1];
  1276.          image[j1]=medium;
  1277. }
  1278. }
  1279. }
  1280. *(pDoc->m_Dib.m_Buffer+i*Width+j)=image[3];
  1281. }
  1282. }
  1283. }
  1284. else
  1285. {
  1286. // 24位彩色中值滤波
  1287. int image[6][4],medium[4],k;
  1288. Height=(int) pDoc->m_Dib.Height;
  1289. Width=(int) pDoc->m_Dib.Bpl;
  1290. for (j=0;j<(int)Width/3;j++)
  1291. {
  1292.   for (i=2;i<Height-2;i++)
  1293. {
  1294. for(k=0;k<3;k++)
  1295. {
  1296. image[3][k+1]=(int) *(m_Buffer1+i*Width+j*3+k);
  1297. image[2][k+1]=(int) *(m_Buffer1+(i-1)*Width+j*3+k);
  1298. image[4][k+1]=(int) *(m_Buffer1+(i+1)*Width+j*3+k);
  1299. image[1][k+1]=(int) *(m_Buffer1+(i-2)*Width+j*3+k);
  1300. image[5][k+1]=(int) *(m_Buffer1+(i+2)*Width+j*3+k);
  1301. }
  1302. //冒泡排序
  1303. for(i1=1;i1<5;i1++)
  1304. {
  1305.          for(j1=i1+1;j1<6;j1++)
  1306. {
  1307.      if((image[j1][1]+image[j1][2]+image[j1][3])/3>(image[i1][1]+image[i1][2]+image[i1][3])/3)
  1308. {
  1309. for(k=1;k<4;k++)
  1310. {
  1311.      medium[k]=image[i1][k];
  1312. image[i1][k]=image[j1][k];
  1313. image[j1][k]=medium[k];
  1314. }
  1315. }
  1316. }
  1317. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=image[3][1];
  1318. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=image[3][2];
  1319. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=image[3][3];
  1320. }
  1321. }
  1322. }
  1323. }
  1324. if(m_Buffer1)
  1325. {
  1326. delete [] m_Buffer1;
  1327. m_Buffer1=0;
  1328. }
  1329. pDoc->m_BM=true;
  1330. Invalidate();
  1331. }
  1333. void CBmpTestView::OnDensityScan() 
  1334. {
  1335. // TODO: Add your command handler code here
  1336. CBmpTestDoc* pDoc = GetDocument();
  1337. ASSERT_VALID(pDoc);  
  1338. # include <math.h>
  1339. int i,j,Height,Width;
  1340. if(pDoc->m_Dib.BPP==8)
  1341. return;
  1342. Height=(int) pDoc->m_Dib.Height;
  1343. Width=(int) pDoc->m_Dib.Bpl;
  1344. int image[2048][4],b[2048],count;
  1345. double diff;
  1346. //扫描从图像底部开始
  1347. for (i=0;i<Height;i++)
  1348. {
  1349. for(j=0;j<(int)Width/3;j++)
  1350. {
  1351. image[j][1]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3);
  1352. image[j][2]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1);
  1353. image[j][3]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2);
  1354. if(image[j][1]<60)
  1355. {
  1356. image[j][1]=0;
  1357. image[j][2]=0;
  1358. image[j][3]=0;
  1359. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1360. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1361. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1362. }
  1363. } 
  1364. j=0;
  1365. count=0;
  1366. while (j<(int)Width/3)
  1367. {
  1368. //因每点三分量均相等，故可简化 diff=sqrt(pow((image[j][1]-image[j+3][1]),2)+pow((image[j][2]-image[j+3][2]),2)+pow((image[j][3]-image[j+3][3]),2));
  1369. diff=abs(image[j][1]-image[j+1][1])*1.7;
  1370. b[j]=(((int)diff>150)?1:0);
  1371. count+=b[j];
  1372. j++;
  1373. }
  1374. //变化15次已是最高限了，再高就会使车牌出现空线
  1375. if (count<15) 
  1376. {
  1377.      for(j=0;j<(int)Width/3;j++)
  1378. {
  1379. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1380. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1381. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1382. }
  1383. }
  1384. else
  1385. {
  1386. //分区域扫描法
  1387. count=0;
  1388. //第一节
  1389. for(j=0;j<(int)Width/45;j++)
  1390. count+=b[j];
  1391. //剩余部分的简化算法
  1392. do
  1393. {
  1394. //变化6次是最佳经验值
  1395. if(count<4)
  1396. {
  1397. *(pDoc->m_Dib.m_Buffer+i*Width+(j-(int)Width/45)*3)=0;
  1398. *(pDoc->m_Dib.m_Buffer+i*Width+(j-(int)Width/45)*3+1)=0;
  1399. *(pDoc->m_Dib.m_Buffer+i*Width+(j-(int)Width/45)*3+2)=0;
  1400. }
  1401. count=count+b[j]-b[j-(int)Width/45];
  1402. j++;
  1403. }while(j<(int)Width/3);
  1404. for(j=j-(int)Width/45;j<(int)Width/3;j++)
  1405. {
  1406. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1407. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1408. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1409. } 
  1410. //进一步扫描，消除图像中的残留点
  1411. for(j=0;j<(int)Width/3;j++)
  1412. {
  1413. image[j][1]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3);
  1414. } 
  1415. j=0;
  1416. count=0;
  1417. while (j<(int)Width/3)
  1418. {
  1419. //因每点三分量均相等，故可简化 diff=sqrt(pow((image[j][1]-image[j+3][1]),2)+pow((image[j][2]-image[j+3][2]),2)+pow((image[j][3]-image[j+3][3]),2));
  1420. diff=abs(image[j][1]-image[j+1][1])*1.7;
  1421. b[j]=(((int)diff>220)?1:0);
  1422. count+=b[j];
  1423. j++;
  1424. }
  1425. //变化3次已是最高限了，再高就会使车牌出现空线
  1426. if (count<8) 
  1427. {
  1428.      for(j=0;j<(int)Width/3;j++)
  1429. {
  1430. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1431. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1432. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1433. } 
  1434. } 
  1435. } 
  1436. }
  1437.   for (i=Height-5;i<Height;i++)
  1438. {
  1439.     for (j=0;j<(int)Width/3 ;j++)
  1440. {
  1441. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1442. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1443. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1444. }
  1445. } 
  1446. //反向扫描
  1447. for (i=0;i<Height;i++)
  1448. {
  1449.     for (j=0;j<(int)Width/3 ;j++)
  1450. {
  1451. image[j][1]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3);
  1452. image[j][2]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1);
  1453. image[j][3]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2);
  1454. } 
  1455. j=(int)Width/3;
  1456. while (j>=1)
  1457. {
  1458. //因每点三分量均相等，故可简化 diff=sqrt(pow((image[j][1]-image[j+3][1]),2)+pow((image[j][2]-image[j+3][2]),2)+pow((image[j][3]-image[j+3][3]),2));
  1459. diff=abs(image[j][1]-image[j-1][1])*1.7;
  1460. b[j]=(((int)diff>60)?1:0);
  1461. j--;
  1462. }
  1463. //变化15次已是最高限了，再高就会使车牌出现空线
  1464. //分区域扫描法
  1465. count=0;
  1466. //第一节
  1467. for(j=(int)Width/3-1;j>=(int)Width/3-(int)Width/45;j--)
  1468. count+=b[j];
  1469. //剩余部分的简化算法
  1470. do
  1471. {
  1472. //变化6次是最佳经验值
  1473. if(count<4)
  1474. {
  1475. *(pDoc->m_Dib.m_Buffer+i*Width+(j+(int)Width/45)*3)=0;
  1476. *(pDoc->m_Dib.m_Buffer+i*Width+(j+(int)Width/45)*3+1)=0;
  1477. *(pDoc->m_Dib.m_Buffer+i*Width+(j+(int)Width/45)*3+2)=0;
  1478. }
  1479. count=count+b[j]-b[j+(int)Width/45];
  1480. j--;
  1481. }while(j>=1); 
  1482. } 
  1483. pDoc->m_BM=true;
  1484. Invalidate();
  1485. }
  1487. void CBmpTestView::OnRowDilation() 
  1488. {
  1489. // TODO: Add your command handler code here
  1490. // 形态学：横向膨胀，结构和,两边各扩展6个像素
  1491. CBmpTestDoc* pDoc = GetDocument();
  1492. ASSERT_VALID(pDoc);  
  1493. if(pDoc->m_Dib.BPP==8)
  1494. return;
  1495. int i,j,i1,j1,Height,Width,a;
  1496. Height=(int) pDoc->m_Dib.Height;
  1497. Width=(int) pDoc->m_Dib.Bpl;
  1498. BYTE * m_Buffer1;
  1499. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  1500.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  1501. {
  1502.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  1503. }
  1504.   for (i=2;i<Height-3;i++)
  1505. {
  1506.     for (j=4;j<(int)Width/3-5;j++)
  1507. {
  1508. a=0;
  1509. for (i1=-2;i1<3;i1++)
  1510. {
  1511. for(j1=-4;j1<5;j1++)
  1512. {
  1513. if ((unsigned int) *(m_Buffer1+(i+i1)*Width+(j+j1)*3)>50)
  1514. a=1;
  1515. }
  1516. }
  1517. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=(a?255:0);
  1518. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=(a?255:0);
  1519. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=(a?255:0);
  1520. }
  1521. }
  1522. if(m_Buffer1)
  1523. {
  1524. delete [] m_Buffer1;
  1525. m_Buffer1=0;
  1526. }
  1527. pDoc->m_BM=true;
  1528. Invalidate();
  1529. }
  1531. void CBmpTestView::OnEdgeColorprewitt1() 
  1532. {
  1533. // TODO: Add your command handler code here
  1534. # include <math.h>
  1535. CBmpTestDoc* pDoc = GetDocument();
  1536. ASSERT_VALID(pDoc);  
  1537. if(pDoc->m_Dib.BPP<=8)
  1538. {
  1539.        AfxMessageBox( "此图像为256色图像，目前该算子还不支持这种图像格式");
  1540.    return;
  1541. }
  1542. int i,j,k,image[5][4],Height,Width;
  1543. unsigned int a;
  1544. BYTE * m_Buffer1;
  1545. double e1,e2;
  1546. //filtering
  1547. Height=(int) pDoc->m_Dib.Height;
  1548. Width=(int) pDoc->m_Dib.Bpl;
  1549. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  1550.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  1551. {
  1552.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  1553. }
  1554.   for (i=0;i<Height-1;i++)
  1555. {
  1556.     for (j=0;j<(int)Width/3 -1;j++)
  1557. {
  1558. for (k=0;k<3;k++)
  1559. {
  1560. image[1][k+1]=(unsigned int) *(m_Buffer1+i*Width+j*3+k);
  1561. image[2][k+1]=(unsigned int) *(m_Buffer1+i*Width+(j+1)*3+k);
  1562. image[3][k+1]=(unsigned int) *(m_Buffer1+(i+1)*Width+j*3+k);
  1563. image[4][k+1]=(unsigned int) *(m_Buffer1+(i+1)*Width+(j+1)*3+k);
  1564. }
  1565.             e1=sqrt(pow((image[1][1]-image[2][1]),2)+pow((image[1][2]-image[2][2]),2)+pow((image[1][3]-image[2][3]),2));
  1566. e2=sqrt(pow((image[3][1]-image[4][1]),2)+pow((image[3][2]-image[4][2]),2)+pow((image[3][3]-image[4][3]),2));
  1567. a=(unsigned int)(e1+e2);
  1568. a=((a>255)?255:a);
  1569. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=a;
  1570. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=a;
  1571. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=a;
  1572. }
  1573. }
  1574. if(m_Buffer1)
  1575. {
  1576. delete [] m_Buffer1;
  1577. m_Buffer1=0;
  1578. }
  1579. pDoc->m_BM=true;
  1580. Invalidate();
  1581. }
  1583. void CBmpTestView::OnEdgeColorrobert1() 
  1584. {
  1585. // TODO: Add your command handler code here
  1587. #include <math.h>
  1588. CBmpTestDoc* pDoc = GetDocument();
  1589. ASSERT_VALID(pDoc);  
  1590. if(pDoc->m_Dib.BPP<=8)
  1591. {
  1592.        AfxMessageBox( "此图像为256色图像，目前该算子还不支持这种图像格式");
  1593.    return;
  1594. }
  1595. int i,j,k,image[5][4],Height,Width;
  1596. unsigned int a;
  1597. BYTE * m_Buffer1;
  1598. double e2,e1;
  1599. //filtering
  1600. Height=(int) pDoc->m_Dib.Height;
  1601. Width=(int) pDoc->m_Dib.Bpl;
  1602. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  1603.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  1604. {
  1605.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  1606. }
  1607.   for (i=0;i<Height-1;i++)
  1608. {
  1609.     for (j=0;j<(int)Width/3 -1;j++)
  1610. {
  1611. for (k=0;k<3;k++)
  1612. {
  1613. image[1][k+1]=(unsigned int) *(m_Buffer1+i*Width+j*3+k);
  1614. image[2][k+1]=(unsigned int) *(m_Buffer1+i*Width+(j+1)*3+k);
  1615. image[3][k+1]=(unsigned int) *(m_Buffer1+(i+1)*Width+j*3+k);
  1616. image[4][k+1]=(unsigned int) *(m_Buffer1+(i+1)*Width+(j+1)*3+k);
  1617. }
  1618.             e1=sqrt(pow((image[1][1]-image[4][1]),2)+pow((image[1][2]-image[4][2]),2)+pow((image[1][3]-image[4][3]),2));
  1619. e2=sqrt(pow((image[2][1]-image[3][1]),2)+pow((image[2][2]-image[3][2]),2)+pow((image[2][3]-image[3][3]),2));
  1620. a=(unsigned int)(e1+e2);
  1621. a=((a>255)?255:a);
  1622. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=a;
  1623. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=a;
  1624. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=a;
  1625. }
  1626. }
  1627. if(m_Buffer1)
  1628. {
  1629. delete [] m_Buffer1;
  1630. m_Buffer1=0;
  1631. }
  1632. pDoc->m_BM=true;
  1633. Invalidate();
  1634. }
  1636. void CBmpTestView::OnRowErosion() 
  1637. {
  1638. // TODO: Add your command handler code here
  1639. // 形态学：横向收缩，结构差,两边各收缩4个像素
  1640. CBmpTestDoc* pDoc = GetDocument();
  1641. ASSERT_VALID(pDoc);  
  1642. if(pDoc->m_Dib.BPP==8)
  1643. return;
  1644. int i,j,i1,j1,Height,Width,a;
  1645. Height=(int) pDoc->m_Dib.Height;
  1646. Width=(int) pDoc->m_Dib.Bpl;
  1647. BYTE * m_Buffer1;
  1648. m_Buffer1=new BYTE[pDoc->m_Dib.BufferSize];
  1649.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  1650. {
  1651.     *(m_Buffer1+j)=*(pDoc->m_Dib.m_Buffer+j);
  1652. }
  1653.   for (i=2;i<Height-3;i++)
  1654. {
  1655.     for (j=4;j<(int)Width/3-5;j++)
  1656. {
  1657. a=0;
  1658. for (i1=-2;i1<3;i1++)
  1659. {
  1660. for(j1=-4;j1<5;j1++)
  1661. {
  1662. if ((unsigned int) *(m_Buffer1+(i+i1)*Width+(j+j1)*3)<50)
  1663. a=1;
  1664. }
  1665. }
  1666. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=(a?0:255);
  1667. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=(a?0:255);
  1668. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=(a?0:255);
  1669. }
  1670. }
  1671. if(m_Buffer1)
  1672. {
  1673. delete [] m_Buffer1;
  1674. m_Buffer1=0;
  1675. }
  1676. pDoc->m_BM=true;
  1677. Invalidate();
  1678. }
  1680. void CBmpTestView::OnClose() 
  1681. {
  1682. // TODO: Add your command handler code here
  1683. OnRowDilation();
  1684. OnRowDilation();
  1685. OnRowErosion();
  1686. // OnRowDilation();
  1687. }
  1689. void CBmpTestView::OnGet() 
  1690. {
  1691. // TODO: Add your command handler code here
  1692. CBmpTestDoc* pDoc = GetDocument();
  1693. ASSERT_VALID(pDoc);  
  1694. if(pDoc->m_Dib.BPP==8)
  1695. return;
  1696. int j;
  1697.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  1698. {
  1699.     *(pDoc->m_Dib.m_Buffer+j)=(*(pDoc->m_Dib.m_Buffer+j))&(*(pDoc->m_Dib.m_Buffer1+j));
  1700. }
  1701. pDoc->m_BM=true;
  1702. Invalidate();
  1703. }
  1705. void CBmpTestView::OnRelaScan() 
  1706. {
  1707. // TODO: Add your command handler code here
  1708. CBmpTestDoc* pDoc = GetDocument();
  1709. ASSERT_VALID(pDoc);  
  1710. # include <math.h>
  1711. int i,j,i1,j1,Height,Width;
  1712. //    char s[15]="";
  1713. if(pDoc->m_Dib.BPP==8)
  1714. return;
  1715. Height=(int) pDoc->m_Dib.Height;
  1716. Width=(int) pDoc->m_Dib.Bpl;
  1717. int image[1024],count;
  1718. double diff;
  1719. CString strr,s;
  1720. LPCTSTR fn;
  1721. CFileDialog fd(FALSE,
  1722.    NULL,
  1723.    NULL,
  1724.    OFN_HIDEREADONLY | OFN_OVERWRITEPROMPT,
  1725.    "data files(*.dat)|*.dat|All files(*.*)|*.*||"
  1726.    );
  1727. if(fd.DoModal()==IDOK)
  1728. {
  1729. strr=fd.GetPathName();
  1730. if(strr.Find(".dat")==-1)
  1731. {
  1732. strr=strr+".dat";
  1733. }
  1734. }
  1735. fn=strr.GetBuffer(100);
  1736. if(fn=="")
  1737. {
  1738. return;
  1739. }
  1740. FILE * fp;
  1741. if((fp=fopen(fn,"wb"))==NULL)
  1742. {
  1743. return;
  1744. }
  1746. //扫描从图像底部开始
  1747. for (i=0;i<Height;i++)
  1748. {
  1749. for(j=0;j<(int)Width/3;j++)
  1750. {
  1751. image[j]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3);
  1752. if(image[j]<60)
  1753. {
  1754. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1755. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1756. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1757. }
  1758. } 
  1759. j=0;
  1760. count=0;
  1761. while (j<(int)Width/3)
  1762. {
  1763. //因每点三分量均相等，故可简化 diff=sqrt(pow((image[j][1]-image[j+3][1]),2)+pow((image[j][2]-image[j+3][2]),2)+pow((image[j][3]-image[j+3][3]),2));
  1764. diff=abs(image[j]-image[j+1])*1.7;
  1765. b[i][j]=(((int)diff>150)?1:0);
  1766. count+=b[i][j];
  1767. j++;
  1768. }
  1769. //变化15次已是最高限了，再高就会使车牌出现空线
  1770. if (count<15) 
  1771. {
  1772.      for(j=0;j<(int)Width/3;j++)
  1773. {
  1774. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1775. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1776. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1777. }
  1778. }
  1779. } 
  1780. //将模板内计数并存入文件中以观查
  1781. for (i=(int)Width/144;i<Height-(int)Width/144;i++)  //车牌纵横比约为3:1
  1782. {
  1783. for(j=(int)Width/48;j<(int)5*Width/16;j++)   //模板大小取车牌最小值，图像宽度的1/8
  1784. {
  1785. //   str(b[i][j],s);   我的数字——>字符串函数,优点是速度快
  1786. count=0;
  1787. for(i1=i-(int)Width/144;i1<i+(int)Width/144;i1++)
  1788. {
  1789. for(j1=j-(int)Width/48;j1<j+(int)Width/48;j1++)
  1790. count+=b[i1][j1];
  1791. }
  1792. s.Format("%d",count);
  1793.     fprintf(fp,"%s ",s); 
  1794. }
  1795.         fprintf(fp,"n"); 
  1796. }
  1797. pDoc->m_BM=true;
  1798. Invalidate();
  1799. fclose(fp);
  1800. }
  1804. void CBmpTestView::OnDensityScanMid() 
  1805. {
  1806. // TODO: Add your command handler code here
  1807. CBmpTestDoc* pDoc = GetDocument();
  1808. ASSERT_VALID(pDoc);  
  1809. # include <math.h>
  1810. int i,j,Height,Width;
  1811. if(pDoc->m_Dib.BPP==8)
  1812. return;
  1813. Height=(int) pDoc->m_Dib.Height;
  1814. Width=(int) pDoc->m_Dib.Bpl;
  1815. int image[2048][4],b[2048],count;
  1816. double diff;
  1817. //扫描从图像底部开始
  1818. for (i=0;i<Height;i++)
  1819. {
  1820. for(j=0;j<(int)Width/3;j++)
  1821. {
  1822. image[j][1]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3);
  1823. image[j][2]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1);
  1824. image[j][3]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2);
  1825. if(image[j][1]<60)
  1826. {
  1827. image[j][1]=0;
  1828. image[j][2]=0;
  1829. image[j][3]=0;
  1830. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1831. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1832. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1833. }
  1834. } 
  1835. j=0;
  1836. count=0;
  1837. while (j<(int)Width/3)
  1838. {
  1839. //因每点三分量均相等，故可简化 diff=sqrt(pow((image[j][1]-image[j+3][1]),2)+pow((image[j][2]-image[j+3][2]),2)+pow((image[j][3]-image[j+3][3]),2));
  1840. diff=abs(image[j][1]-image[j+1][1])*1.7;
  1841. b[j]=(((int)diff>150)?1:0);
  1842. count+=b[j];
  1843. j++;
  1844. }
  1845. //变化15次已是最高限了，再高就会使车牌出现空线
  1846. if (count<15) 
  1847. {
  1848.      for(j=0;j<(int)Width/3;j++)
  1849. {
  1850. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1851. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1852. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1853. }
  1854. }
  1855. else
  1856. {
  1857. //分区域扫描法
  1858. count=0;
  1859. //第一节
  1860. for(j=0;j<(int)Width/45;j++)
  1861. {
  1862. count+=b[j];
  1863. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1864. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1865. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1866. }
  1867. //剩余部分的简化算法
  1868. do
  1869. {
  1870. //变化6次是最佳经验值
  1871. if(count<5)
  1872. {
  1873. *(pDoc->m_Dib.m_Buffer+i*Width+(j-(int)Width/90)*3)=0;
  1874. *(pDoc->m_Dib.m_Buffer+i*Width+(j-(int)Width/90)*3+1)=0;
  1875. *(pDoc->m_Dib.m_Buffer+i*Width+(j-(int)Width/90)*3+2)=0;
  1876. }
  1877. count=count+b[j]-b[j-(int)Width/45];
  1878. j++;
  1879. }while(j<(int)Width/3);
  1880. for(j=j-(int)Width/45;j<(int)Width/3;j++)
  1881. {
  1882. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1883. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1884. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1885. } 
  1886. //进一步扫描，消除图像中的残留点
  1887. for(j=0;j<(int)Width/3;j++)
  1888. {
  1889. image[j][1]=(unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3);
  1890. } 
  1891. j=0;
  1892. count=0;
  1893. while (j<(int)Width/3)
  1894. {
  1895. //因每点三分量均相等，故可简化 diff=sqrt(pow((image[j][1]-image[j+3][1]),2)+pow((image[j][2]-image[j+3][2]),2)+pow((image[j][3]-image[j+3][3]),2));
  1896. diff=abs(image[j][1]-image[j+1][1])*1.7;
  1897. b[j]=(((int)diff>150)?1:0);
  1898. count+=b[j];
  1899. j++;
  1900. }
  1901. //变化3次已是最高限了，再高就会使车牌出现空线
  1902. if (count<6) 
  1903. {
  1904.      for(j=0;j<(int)Width/3;j++)
  1905. {
  1906. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1907. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1908. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1909. } 
  1910. } 
  1911. } 
  1912. }
  1913.   for (i=Height-5;i<Height;i++)
  1914. {
  1915.     for (j=0;j<(int)Width/3 ;j++)
  1916. {
  1917. *(pDoc->m_Dib.m_Buffer+i*Width+j*3)=0;
  1918. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+1)=0;
  1919. *(pDoc->m_Dib.m_Buffer+i*Width+j*3+2)=0;
  1920. }
  1921. } 
  1922. pDoc->m_BM=true;
  1923. Invalidate();
  1925. }
  1928. void rgbToHSV(float r,float g,float b,float *h,float *s,float *v)
  1929. {
  1930. float max=MAX(r,MAX(g,b)),min=MIN(r,MIN(g,b));
  1931. float delta=max-min;
  1932. *v=max;
  1933. if(max!=0.0) *s=delta/max;
  1934. else *s=0.0;
  1935. if(*s==0.0) *h=NO_HUE
  1936. else
  1937. {
  1938. if(r==max) *h=(g-b)/delta;
  1939. else if(g==max) *h=2+(b-r)/delta;
  1940. else  *h=4+(r-g)/delta;
  1941. *h=*h*60;
  1942. if(*h<0) *h=*h+360;
  1943. }
  1944. }
  1947. void CBmpTestView::OnAreaMark() 
  1948. {
  1949. // TODO: Add your command handler code here
  1950. CBmpTestDoc* pDoc = GetDocument();
  1951. ASSERT_VALID(pDoc);  
  1952. HWND      hWnd;
  1953. int left[100],right[100],top[100],bottom[100],*pIntensity,*pSaturation,*pHue,SaturationI,HueI,IntensityI,i,j,k,premark,min,mark,Blue,Red,Green,markmax,Height,Width,Wid,Hei,*pHeight,*pWidth,*Widmax,*Heimax,CountofBlue,CountofGreen,CountofRed,maxofIntensity,maxofHue,maxofSat,max;
  1954. bool begin;
  1955. float Saturation,Hue,Intensity, ExpofIntensity,SigmaofIntensity,ExpofHue,SigmaofHue,ExpofSat,SigmaofSat;
  1956. float T0ofInt,T0,u1,u2,T1,w1,w2,u,sig1,sig2,u1ofInt,sig1ofInt,u2ofInt,sig2ofInt,distance,simulate;
  1957. CString s;
  1958. POINT  CurP,NeighborP,SeedP;
  1959. BYTE *lpTempPtr,*lpTempPtr1;
  1960. MYSTACK SeedFillStack;
  1961. if(pDoc->m_Dib.BPP<16)
  1962. return;
  1963. Height=(int) pDoc->m_Dib.Height;
  1964. Width=(int) pDoc->m_Dib.Bpl;
  1965. //建立堆栈，大小与图像相同
  1966. SeedFillStack.InitStack(hWnd,Height*Width);
  1967. mark=0;
  1968. //初始化与标号
  1969. for (i=0;i<Height;i++)
  1970. {
  1971. for(j=0;j<(int)Width/3;j++)
  1972. {
  1973. b[i][j]=0;
  1974. }
  1975. }
  1976. for (i=0;i<Height;i++)
  1977. {
  1978. for(j=0;j<(int)Width/3;j++)
  1979. {
  1980. //种子法和迷宫法并求出每个区域的最左，最右，最上，最下
  1981. if((unsigned int) *(pDoc->m_Dib.m_Buffer+i*Width+j*3)!=0 && b[i][j]==0)
  1982. {
  1983. mark++;
  1984. left[mark]=(int)Width/3-1;right[mark]=0;bottom[mark]=Height-1;top[mark]=0;
  1985. SeedP.x=j;
  1986. SeedP.y=i;
  1987. SeedFillStack.MyPush(SeedP);
  1988. while(!SeedFillStack.IsStackEmpty())
  1989. {
  1990. CurP=SeedFillStack.MyPop();
  1991. b[CurP.y][CurP.x]=mark;
  1992. if(CurP.y>top[mark])
  1993. top[mark]=CurP.y;
  1994. if(CurP.y<bottom[mark])
  1995. bottom[mark]=CurP.y;
  1996. if(CurP.x>right[mark])
  1997. right[mark]=CurP.x;
  1998. if(CurP.x<left[mark])
  1999. left[mark]=CurP.x;
  2000. lpTempPtr=pDoc->m_Dib.m_Buffer+CurP.y*Width+CurP.x*3;
  2001. *lpTempPtr=255;
  2002. *(lpTempPtr+1)=255;
  2003. *(lpTempPtr+2)=255;
  2004. //标记颜色
  2005. /* switch(mark%4)
  2006. {
  2007. case 0:
  2008. *lpTempPtr=255;
  2009. *(lpTempPtr+1)=0;
  2010. *(lpTempPtr+2)=0;
  2011. break;
  2012. case 1:
  2013. *lpTempPtr=0;
  2014. *(lpTempPtr+1)=255;
  2015. *(lpTempPtr+2)=0;
  2016. break;
  2017. case 2:
  2018. *lpTempPtr=0;
  2019. *(lpTempPtr+1)=0;
  2020. *(lpTempPtr+2)=255;
  2021. break;
  2022. case 3:
  2023. *lpTempPtr=127;
  2024. *(lpTempPtr+1)=127;
  2025. *(lpTempPtr+2)=127;
  2026. break;
  2027. } */
  2028. //left neighbour
  2029. if(CurP.x>0 )
  2030. {
  2031. NeighborP.x=CurP.x-1;
  2032. NeighborP.y=CurP.y;
  2033. lpTempPtr1=lpTempPtr-3;
  2034. if(*lpTempPtr1!=0 && b[NeighborP.y][NeighborP.x]==0)
  2035. SeedFillStack.MyPush(NeighborP);
  2036. }
  2037. //down neighbour
  2038. if(CurP.y>0)
  2039. {
  2040. NeighborP.x=CurP.x;
  2041. NeighborP.y=CurP.y-1;
  2042. lpTempPtr1=lpTempPtr-Width;
  2043. if(*lpTempPtr1!=0 && b[NeighborP.y][NeighborP.x]==0)
  2044. SeedFillStack.MyPush(NeighborP);
  2045. }
  2046. //right neighbour
  2047. if(CurP.x<(int)Width/3-1)
  2048. {
  2049. NeighborP.x=CurP.x+1;
  2050. NeighborP.y=CurP.y;
  2051. lpTempPtr1=lpTempPtr+3;
  2052. if(*lpTempPtr1!=0 && b[NeighborP.y][NeighborP.x]==0)
  2053. SeedFillStack.MyPush(NeighborP);
  2054. }
  2055. //up neighbour
  2056. if(CurP.y<Height-1)
  2057. {
  2058. NeighborP.x=CurP.x;
  2059. NeighborP.y=CurP.y+1;
  2060. lpTempPtr1=lpTempPtr+Width;
  2061. if(*lpTempPtr1!=0 && b[NeighborP.y][NeighborP.x]==0)
  2062. SeedFillStack.MyPush(NeighborP);
  2063. }
  2064. }
  2065. }
  2066. }  //列结束
  2067. } //行结束
  2068. SeedFillStack.DeInitStack();
  2069. //计算每一区域的平均高度和宽度及区域面积
  2070. markmax=mark+1;
  2071. pWidth=new int[markmax];
  2072. pHeight=new int[markmax];
  2073. Heimax=new int[markmax];
  2074. Widmax=new int[markmax];
  2075. pIntensity=new int[256*markmax];
  2076. pSaturation=new int[256*markmax];
  2077. pHue=new int[256*markmax];
  2078. //先计算各个区域的平均宽度
  2079. for(i=0;i<markmax;i++)
  2080. {
  2081. pWidth[i]=0;
  2082. pHeight[i]=0;
  2083. Heimax[i]=0;
  2084. Widmax[i]=0;
  2085. for (j=0;j<256;j++)
  2086. {
  2087. pIntensity[j+i*256]=0;
  2088. pSaturation[j+i*256]=0;
  2089. pHue[j+i*256]=0;
  2090. }
  2091. }
  2092. for (i=0;i<Height;i++)
  2093. {
  2094. Wid=0;
  2095. premark=0;
  2096. for(j=0;j<(int)Width/3;j++)
  2097. {
  2098. mark=b[i][j];
  2099. if(mark!=0)
  2100. {
  2101. if(mark!=premark)
  2102. {
  2103. pWidth[mark]++;
  2104. Wid++;
  2105. premark=mark;
  2106. }
  2107. else
  2108. {
  2109. pWidth[mark]++;
  2110. Wid++;
  2111. }
  2112. }
  2113. else
  2114. {
  2115. if((Wid>Widmax[premark])&&(premark!=0))
  2116. Widmax[premark]=Wid;
  2117. Wid=0;
  2118. premark=0;
  2119. } 
  2120. }
  2121. } 
  2122. //计算各区域的平均高度
  2123. for(j=0;j<(int)Width/3;j++)
  2124. {
  2125. Hei=0;
  2126. premark=0;
  2127. for (i=0;i<Height;i++)
  2128. {
  2129. mark=b[i][j];
  2130. if(mark!=0)
  2131. {
  2132. if(mark!=premark)
  2133. {
  2134. pHeight[mark]++;
  2135. Hei++;
  2136. premark=mark;
  2137. }
  2138. else
  2139. {
  2140. pHeight[mark]++;
  2141. Hei++;
  2142. }
  2143. }
  2144. else
  2145. {
  2146. if((Hei>Heimax[premark])&&(premark!=0))
  2147. Heimax[premark]=Hei;
  2148. Hei=0;
  2149. premark=0;
  2150. } 
  2151. }
  2152. } 
  2153. //删除无用的区域
  2154. for(i=1;i<markmax;i++)
  2155. {
  2156. //求各区域平均宽度及平均高度
  2157. pWidth[i]/=Heimax[i];
  2158. pHeight[i]/=Widmax[i];
  2159. if(pWidth[i]>Width/9||pWidth[i]<Width/30)
  2160. Widmax[i]=0;
  2161. else
  2162. {
  2163. if(abs(pWidth[i]/pHeight[i]-2.5)>2)
  2164. Widmax[i]=0;
  2165. }
  2166. }
  2167. for (i=0;i<Height;i++)
  2168. {
  2169. for(j=0;j<(int)Width/3;j++)
  2170. {
  2171. mark=b[i][j];
  2172. if(mark!=0)
  2173. {
  2174. if(Widmax[mark]==0)
  2175. {
  2176. lpTempPtr=pDoc->m_Dib.m_Buffer+i*Width+j*3;
  2177. *lpTempPtr=0;
  2178. *(lpTempPtr+1)=0;
  2179. *(lpTempPtr+2)=0;
  2180. b[i][j]=0;
  2181. }
  2182. else
  2183. {
  2184. lpTempPtr=pDoc->m_Dib.m_Buffer+i*Width+j*3;
  2185. *lpTempPtr=255;
  2186. *(lpTempPtr+1)=255;
  2187. *(lpTempPtr+2)=255;
  2188. }
  2189. }
  2190. }
  2191. }
  2192.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  2193. {
  2194.     *(pDoc->m_Dib.m_Buffer+j)=(*(pDoc->m_Dib.m_Buffer+j))&(*(pDoc->m_Dib.m_Buffer1+j));
  2195. }
  2196. pDoc->m_BM=true;
  2197. Invalidate();
  2198. AfxMessageBox("下面是光强、色度、色饱和度直方图显示！");
  2199. m_Gradmax=0;
  2200.   for (i=0;i<Height;i++)
  2201. {
  2202. for(j=0;j<(int)Width/3;j++)
  2203. {
  2204. mark=b[i][j];
  2205. if(mark!=0)
  2206. {
  2207. lpTempPtr=pDoc->m_Dib.m_Buffer+i*Width+j*3;
  2208.     Blue=(int) *lpTempPtr;
  2209. Green=(int) *(lpTempPtr+1);
  2210. Red=(int) *(lpTempPtr+2);
  2211. rgbToHSV(Red,Green,Blue,&Hue,&Saturation,&Intensity);
  2212. //如下是三个量的量化过程
  2213. HueI=(int) (Hue*255/360);
  2214. SaturationI=(int) (Saturation*255);
  2215. IntensityI= (int) Intensity;
  2216. // 直方图统计过程
  2217. pIntensity[IntensityI+mark*256]++;
  2218. pHue[HueI+mark*256]++;
  2219. pSaturation[SaturationI+mark*256]++; 
  2220. }
  2221. }
  2222. }
  2223. for(j=0;j<markmax;j++)
  2224. {
  2225. m_Gradmax=0;
  2226. CountofBlue=0;
  2227. CountofGreen=0;
  2228. CountofRed=0;
  2229.     for ( i=0;i<256;i++)
  2230. {
  2231. m_Grade_Blue[i]=pIntensity[i+j*256];
  2232. m_Grade_Green[i]=pHue[i+j*256];
  2233. m_Grade_Red[i]=pSaturation[i+j*256];
  2234. if (m_Grade_Blue[i]>m_Gradmax)
  2235. m_Gradmax=m_Grade_Blue[i];
  2236. if (m_Grade_Green[i]>m_Gradmax)
  2237. m_Gradmax=m_Grade_Green[i];
  2238. if (m_Grade_Red[i]>m_Gradmax)
  2239. m_Gradmax=m_Grade_Red[i];
  2240. CountofBlue=CountofBlue+m_Grade_Blue[i];
  2241. CountofGreen=CountofGreen+m_Grade_Green[i];
  2242. CountofRed=CountofRed+m_Grade_Red[i];
  2243. }
  2245. if (m_Gradmax>0)
  2246. {
  2247. pDoc->m_BM=false;
  2248. Invalidate();
  2249. ExpofIntensity=0;
  2250. ExpofHue=0;
  2251. SigmaofIntensity=0;
  2252. SigmaofHue=0;
  2253. ExpofSat=0;SigmaofSat=0;
  2254. // Calculate estimate and maxpoint of Intensity and Hue
  2255. max=0;min=0;
  2256.     for ( i=0;i<256;i++)
  2257. {
  2258. if(max<m_Grade_Blue[i])
  2259. {
  2260. max=m_Grade_Blue[i];
  2261. maxofIntensity=i;
  2262. }
  2263. if(min<m_Grade_Red[i])
  2264. {
  2265. min=m_Grade_Red[i];
  2266. maxofSat=i;
  2267. }
  2268. ExpofIntensity=ExpofIntensity+(float) i*m_Grade_Blue[i]/CountofBlue;
  2269. ExpofHue=ExpofHue+(float) i*m_Grade_Green[i]/CountofGreen;
  2270. ExpofSat=ExpofSat+(float) i*m_Grade_Red[i]/CountofRed;
  2271. }
  2272. max=0;
  2273. for ( i=0;i<256;i++)
  2274. {
  2275. if(max<m_Grade_Green[i])
  2276. {
  2277. max=m_Grade_Green[i];
  2278. maxofHue=i;
  2279. }
  2280. SigmaofIntensity=SigmaofIntensity+(float) (i-ExpofIntensity)*(i-ExpofIntensity)*m_Grade_Blue[i]/CountofBlue;
  2281. SigmaofHue=SigmaofHue+(float) (i-ExpofHue)*(i-ExpofHue)*m_Grade_Green[i]/CountofGreen;
  2282. SigmaofSat=SigmaofSat+(float) (i-ExpofSat)*(i-ExpofSat)*m_Grade_Red[i]/CountofRed;
  2283. }
  2284. SigmaofIntensity=sqrt(SigmaofIntensity);
  2285. SigmaofHue=sqrt(SigmaofHue);
  2286. SigmaofSat=sqrt(SigmaofSat);
  2287. //estimate if there is dual peak of Intensity
  2288. T0=ExpofIntensity;T1=0;
  2289. while (fabs(T0-T1)>0.01)
  2290. {
  2291. w1=0;w2=0;u1=0;u2=0;sig1=0;sig2=0;
  2292. for(i=0;i<=int(T0);i++)
  2293. {
  2294. w1=w1+(float)m_Grade_Blue[i]/CountofBlue;
  2295. u1=u1+(float)i*m_Grade_Blue[i]/CountofBlue;
  2296. }
  2297. u1=u1/w1;
  2298. for(i=0;i<=int(T0);i++)
  2299. {
  2300. sig1=sig1+(float)(i-u1)*(i-u1)*m_Grade_Blue[i]/CountofBlue;
  2301. }
  2302. for(i=int(T0)+1;i<256;i++)
  2303. {
  2304. w2=w2+(float)m_Grade_Blue[i]/CountofBlue;
  2305. u2=u2+(float)i*m_Grade_Blue[i]/CountofBlue;
  2306. }
  2307. u2=u2/w2;
  2308. for(i=int(T0)+1;i<256;i++)
  2309. {
  2310. sig2=sig2+(float)(i-u2)*(i-u2)*m_Grade_Blue[i]/CountofBlue;
  2311. }
  2312. T1=T0;T0=0.5*(u1+u2);sig1=sqrt(sig1);sig2=sqrt(sig2);
  2313. }
  2314. T0ofInt=T0;u1ofInt=u1;sig1ofInt=sig1;u2ofInt=u2;sig2ofInt=sig2;
  2315. //estimate if there is dual peak of Hue
  2316. T0=ExpofHue;T1=0;
  2317. while (fabs(T0-T1)>0.01)
  2318. {
  2319. w1=0;w2=0;u1=0;u2=0;sig1=0;sig2=0;
  2320. for(i=0;i<=int(T0);i++)
  2321. {
  2322. w1=w1+(float)m_Grade_Green[i]/CountofGreen;
  2323. u1=u1+(float)i*m_Grade_Green[i]/CountofGreen;
  2324. }
  2325. u1=u1/w1;
  2326. for(i=0;i<=int(T0);i++)
  2327. {
  2328. sig1=sig1+(float)(i-u1)*(i-u1)*m_Grade_Green[i]/CountofGreen;
  2329. }
  2330. for(i=int(T0)+1;i<256;i++)
  2331. {
  2332. w2=w2+(float)m_Grade_Green[i]/CountofGreen;
  2333. u2=u2+(float)i*m_Grade_Green[i]/CountofGreen;
  2334. }
  2335. u2=u2/w2;
  2336. for(i=int(T0)+1;i<256;i++)
  2337. {
  2338. sig2=sig2+(float)(i-u2)*(i-u2)*m_Grade_Green[i]/CountofGreen;
  2339. }
  2340. T1=T0;T0=0.5*(u1+u2);sig1=sqrt(sig1);sig2=sqrt(sig2);
  2341. }
  2342. s.Format("ExpofIntensity=%f,SigmaofIntensity=%f,T0ofIntensity=%f,u1ofInt=%f,sig1ofInt=%f,u2ofInt=%f,sig2ofInt=%f,!!!!!!,ExpofHue=%f,SigmaofHue=%f,T0=%f,u1=%f,sig1=%f,u2=%f,sig2=%f,%d,%d",ExpofIntensity,SigmaofIntensity,T0ofInt,u1ofInt,sig1ofInt,u2ofInt,sig2ofInt,ExpofHue,SigmaofHue,T0,u1,sig1,u2,sig2,maxofIntensity,maxofHue);
  2343. AfxMessageBox(s);
  2344. // judge if Intensity is dual peak or single peak
  2345. if (abs(ExpofIntensity-maxofIntensity)>15)
  2346. // dual peak
  2347.             {
  2348. if (abs(u1ofInt-maxofIntensity)<15)
  2349. {
  2350. ExpofIntensity=u1ofInt;
  2351. SigmaofIntensity=sig1ofInt;
  2352. }
  2353. else
  2354. {
  2355. if (abs(u2ofInt-maxofIntensity)<15)
  2356. {
  2357. ExpofIntensity=u2ofInt;
  2358. SigmaofIntensity=sig2ofInt;
  2359. }
  2360. // else
  2361. // ExpofIntensity=-1;
  2362. }
  2363. }
  2364. // judge if Hue is dual peak or single peak
  2365. if (abs(ExpofHue-maxofHue)>15)
  2366. // dual peak
  2367.             {
  2368. if (abs(u1-maxofHue)<15)
  2369. {
  2370. ExpofHue=u1;
  2371. SigmaofHue=sig1;
  2372. }
  2373. else
  2374. {
  2375. if (abs(u2-maxofHue)<15)
  2376. {
  2377. ExpofHue=u2;
  2378. SigmaofHue=sig2;
  2379. }
  2380. else
  2381. ExpofHue=-1;
  2382. }
  2383. }
  2384. s.Format("ExpofIntensity=%f,SigmaofIntensity=%f,ExpofHue=%f,SigmaofHue=%f,ExpofSat=%f,SigmaofSat=%f",ExpofIntensity,SigmaofIntensity,ExpofHue,SigmaofHue,ExpofSat,SigmaofSat);
  2385. AfxMessageBox(s);
  2386. // estimate the color of vehicle plate
  2387. //deal with blue
  2388. begin=false;
  2389. if(abs(ExpofHue-150)<19||(ExpofHue==-1&&abs(maxofHue-140)<9))
  2390. {
  2391. if((ExpofIntensity>60)&&(ExpofIntensity<242))
  2392. {
  2393. s.Format("车牌为兰色");
  2394. begin=true;
  2395. }
  2396. else
  2397. {
  2398. if(ExpofIntensity<=60)
  2399. {
  2400. s.Format("车牌为黑色");
  2401. begin=true;
  2402. }
  2403. else
  2404. {
  2405. s.Format("无用区域，删除");
  2406. for (i=0;i<Height;i++)
  2407. {
  2408. for(k=0;k<(int)Width/3;k++)
  2409. {
  2410. mark=b[i][k];
  2411. if(mark==j)
  2412. {
  2413. lpTempPtr=pDoc->m_Dib.m_Buffer+i*Width+k*3;
  2414. *lpTempPtr=0;
  2415. *(lpTempPtr+1)=0;
  2416. *(lpTempPtr+2)=0;
  2417. b[i][k]=0;
  2418. }
  2419. }
  2420. }
  2421. }
  2422. }
  2423. }
  2424. else
  2425. { //deal with yellow
  2426. if(abs(ExpofHue-33)<5||(ExpofHue==-1&&abs(maxofHue-33)<5))
  2427. {
  2428. if ((ExpofIntensity>60)&&(ExpofIntensity<242))
  2429. {
  2430. s.Format("车牌为黄色");
  2431. begin=true;
  2432. }
  2433. }
  2434. else
  2435. {
  2436. if(ExpofIntensity<=60&&ExpofIntensity>=0)
  2437. {
  2438. s.Format("车牌为黑色");
  2439. begin=true;
  2440. }
  2441. else
  2442. {
  2443. if(ExpofIntensity>242||(ExpofIntensity==-1&&maxofIntensity>242))
  2444. {
  2445. s.Format("车牌为白色");
  2446. begin=true;
  2447. }
  2448. else
  2449. {
  2450. s.Format("无用区域，删除");
  2451. for (i=0;i<Height;i++)
  2452. {
  2453. for(k=0;k<(int)Width/3;k++)
  2454. {
  2455. mark=b[i][k];
  2456. if(mark==j)
  2457. {
  2458. lpTempPtr=pDoc->m_Dib.m_Buffer+i*Width+k*3;
  2459. *lpTempPtr=0;
  2460. *(lpTempPtr+1)=0;
  2461. *(lpTempPtr+2)=0;
  2462. b[i][k]=0;
  2463. }
  2464. }
  2465. }
  2466. }
  2467. }
  2468. }
  2469. }
  2470. AfxMessageBox(s);
  2471. if(begin)
  2472. {
  2473.   for (i=bottom[j]-3;i<=top[j]+3;i++)
  2474. {
  2475. for(k=left[j];k<=right[j];k++)
  2476. {
  2477. lpTempPtr=pDoc->m_Dib.m_Buffer+i*Width+k*3;
  2478. *lpTempPtr=255;
  2479. *(lpTempPtr+1)=255;
  2480. *(lpTempPtr+2)=255;
  2481. }
  2482. }
  2483. // s.Format("left=%d,right=%d,bottom=%d,top=%d",left[j],right[j],bottom[j],top[j]);
  2484. // AfxMessageBox(s);
  2485.   for (i=bottom[j]-3;i<=top[j]+3;i++)
  2486. {
  2487. for(k=left[j];k<=right[j];k++)
  2488. {
  2489. lpTempPtr=pDoc->m_Dib.m_Buffer1+i*Width+k*3;
  2490.     Blue=(int) *lpTempPtr;
  2491. Green=(int) *(lpTempPtr+1);
  2492. Red=(int) *(lpTempPtr+2);
  2493. rgbToHSV(Red,Green,Blue,&Hue,&Saturation,&Intensity);
  2494. //如下是三个量的量化过程
  2495. HueI=(int) (Hue*255/360);
  2496. SaturationI=(int) (Saturation*255);
  2497. IntensityI= (int) Intensity;
  2498. lpTempPtr=pDoc->m_Dib.m_Buffer+i*Width+k*3;
  2499. distance=sqrt(pow(SaturationI*cos(HueI*2*pi/255)/255-ExpofSat*cos(ExpofHue*2*pi/255)/255,2)+pow(SaturationI*sin(HueI*2*pi/255)/255-ExpofSat*sin(ExpofHue*2*pi/255)/255,2)+pow((IntensityI-ExpofIntensity)/255,2));
  2500. simulate=1-distance/sqrt(5);
  2501. // s.Format("simulate=%f",simulate);
  2502. // AfxMessageBox(s);
  2503. if(simulate>0.85)
  2504. {
  2505. *lpTempPtr=255;
  2506. *(lpTempPtr+1)=255;
  2507. *(lpTempPtr+2)=255;
  2508. }
  2509. else
  2510. {
  2511. *lpTempPtr=0;
  2512. *(lpTempPtr+1)=0;
  2513. *(lpTempPtr+2)=0;
  2514. }
  2515. }
  2516. }
  2517. }
  2518. }
  2519. }
  2520.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j++)
  2521. {
  2522.     *(pDoc->m_Dib.m_Buffer+j)=(*(pDoc->m_Dib.m_Buffer+j))&(*(pDoc->m_Dib.m_Buffer1+j));
  2523. }
  2524. delete []pWidth;
  2525. delete []pHeight;
  2526. delete []Widmax;
  2527. delete []Heimax;
  2528. delete []pIntensity;
  2529. delete []pSaturation;
  2530. delete []pHue;
  2531. pDoc->m_BM=true;
  2532. Invalidate();
  2533. }
  2535. void CBmpTestView::OnLocating() 
  2536. {
  2537. // TODO: Add your command handler code here
  2538. OnEdgeColorprewitt1();
  2539. OnDensityScanMid();
  2540. OnMidFiltLine();
  2541. OnClose();
  2542. OnAreaMark();
  2543. }
  2545. void CBmpTestView::OnVehicleBitmap() 
  2546. {
  2547. // TODO: Add your command handler code here
  2548. CBmpTestDoc* pDoc = GetDocument();
  2549. ASSERT_VALID(pDoc);  
  2550. int i,j,Height,Width,th;
  2551. long s;
  2552. if(!pDoc->m_Dib.m_BmpLoaded)
  2553. {
  2554. AfxMessageBox("还没有打开图像文件");
  2555. return;
  2556. }
  2557. if(pDoc->m_Dib.BPP==8)
  2558. return;
  2559. Height=(int) pDoc->m_Dib.Height;
  2560. Width=(int) pDoc->m_Dib.Bpl;
  2561. s=0;
  2562.     for ( i=0;i<256;i++)
  2563. {
  2564. m_Grade_Blue[i]=0;
  2565. }
  2566.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j=j+3)
  2567. {
  2568.     m_Grade_Blue[int(*(pDoc->m_Dib.m_Buffer+j))]++;
  2569. }
  2570. i=255;
  2571. do {
  2572. s+=m_Grade_Blue[i];
  2573. i--;
  2574. }while(s<(int) pDoc->m_Dib.BufferSize/120);
  2575. th=i;
  2576.     for ( j=0;j<(int) pDoc->m_Dib.BufferSize;j=j+3)
  2577. {
  2578.     if(*(pDoc->m_Dib.m_Buffer+j)>th)
  2579. {
  2580. // *(pDoc->m_Dib.m_Buffer+j)=255;
  2581. // *(pDoc->m_Dib.m_Buffer+j+1)=255;
  2582. // *(pDoc->m_Dib.m_Buffer+j+2)=255;
  2583. }
  2584. else
  2585. {
  2586. *(pDoc->m_Dib.m_Buffer+j)=0;
  2587. *(pDoc->m_Dib.m_Buffer+j+1)=0;
  2588. *(pDoc->m_Dib.m_Buffer+j+2)=0;
  2589. }
  2590. }
  2591. pDoc->m_BM=true;
  2592. Invalidate();
  2593. }
  2595. void CBmpTestView::OnBlueCut() 
  2596. {
  2597. // TODO: Add your command handler code here
  2598. CBmpTestDoc* pDoc = GetDocument();
  2599. ASSERT_VALID(pDoc);  
  2600. int i,j,k,Blue,Red,Green,Height,Width;
  2601. float Saturation,Hue,Intensity;
  2602. CString s;
  2603. BYTE *lpTempPtr,*lpTempPtr1;
  2604. if(pDoc->m_Dib.BPP<16)
  2605. return;
  2606. Height=(int) pDoc->m_Dib.Height;
  2607. Width=(int) pDoc->m_Dib.Bpl;
  2608.   for (i=0;i<Height;i++)
  2609. {
  2610. for(j=0;j<(int)Width/3;j++)
  2611. {
  2612. lpTempPtr=pDoc->m_Dib.m_Buffer+i*Width+j*3;
  2613.     Blue=(int) *lpTempPtr;
  2614. Green=(int) *(lpTempPtr+1);
  2615. Red=(int) *(lpTempPtr+2);
  2616. rgbToHSV(Red,Green,Blue,&Hue,&Saturation,&Intensity);
  2617. if((Hue>240)||(Hue<200))
  2618. {
  2619. *lpTempPtr=0;
  2620. *(lpTempPtr+1)=0;
  2621. *(lpTempPtr+2)=0;
  2622. }
  2623. }
  2624. }
  2625. pDoc->m_BM=true;
  2626. Invalidate();
  2627. }