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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 图形/文字识别 > 查看源码
dipView.cpp
资源名称：vcc.rar　[点击查看]
上传用户：cdscwht
上传日期：2022-07-27
资源大小：264k
文件大小：79k
源码类别：

图形/文字识别

开发平台：

Visual Basic

dipView.cpp：源码内容
  1. // dipView.cpp : implementation of the CDipView class
  2. //
  4. #include "stdafx.h"
  5. #include "dip.h"
  7. #include "dipDoc.h"
  8. #include "dipView.h"
  9. #include "afxdlgs.h"
  10. #include "afxwin.h"
  11. #include "math.h"
  12. #include ".dipview.h"
  13. #include "stdlib.h"
  14. #include "malloc.h"
  15. #include "mainfrm.h"
  16. #include "charview.h"
  17. #include "recvchar.h"
  19. #include "training.h"
  22. #ifdef _DEBUG
  23. #define new DEBUG_NEW
  24. #undef THIS_FILE
  25. static char THIS_FILE[] = __FILE__;
  26. #endif
  28. /////////////////////////////////////////////////////////////////////////////
  29. // CDipView
  31. IMPLEMENT_DYNCREATE(CDipView, CView)
  33. BEGIN_MESSAGE_MAP(CDipView, CView)
  34. //{{AFX_MSG_MAP(CDipView)
  35. ON_COMMAND(ID_FILE_OPEN, OnFileOpen)
  36. ON_WM_PAINT()
  37. ON_COMMAND(ID_FILE_SAVE, OnFileSave)
  38. ON_COMMAND(ID_reopen, Onreopen)
  39. ON_WM_MOUSEMOVE()
  40. ON_WM_CREATE()
  41. ON_COMMAND(ID_LINE, OnLine)
  42. ON_COMMAND(ID_F3, Onfse)
  43. ON_COMMAND(ID_F4, Onfz)
  44. ON_COMMAND(ID_MoHu, OnMoHu)
  45. ON_COMMAND(ID_ForHandle, OnForHandle)
  46. ON_COMMAND(ID_HSearch, OnHSearch)
  47. ON_COMMAND(ID_Banlance, OnBanlance)
  48. ON_COMMAND(ID_SubVertical, OnSubVertical)
  49. ON_COMMAND(ID_SubHorizontal, OnSubHorizontal)
  50. ON_COMMAND(ID_MidValue, OnMidValue)
  51. ON_COMMAND(ID_Average, OnAverage)
  52. ON_COMMAND(ID_Bool, OnBool)
  53. ON_COMMAND(ID_Choose, OnChoose)
  54. ON_COMMAND(ID_TRAINFROMVEHICLE, OnTrainfromvehicle)
  55. ON_COMMAND(ID_SAVESA, OnSave)
  56. ON_COMMAND(ID_READSA, OnLoad)
  57. ON_COMMAND(ID_CHANGE, OnChange)
  58. ON_COMMAND(ID_TRAIN, OnTrain)
  59. ON_COMMAND(ID_SAVEBP, OnSavebp)
  60. ON_COMMAND(ID_VSearch, OnVSearch)
  61. ON_COMMAND(ID_Locate, OnLocate)
  62. ON_COMMAND(ID_Sobel, OnSobel)
  63. ON_COMMAND(ID_VSplit, OnVSplit)
  64. ON_COMMAND(ID_CharacterRec, OnCharacterRecg)
  65. ON_COMMAND(ID_READBP, OnReadbp)
  66. //}}AFX_MSG_MAP
  67. // Standard printing commands
  68. ON_COMMAND(ID_FILE_PRINT, CView::OnFilePrint)
  69. ON_COMMAND(ID_FILE_PRINT_DIRECT, CView::OnFilePrint)
  70. ON_COMMAND(ID_FILE_PRINT_PREVIEW, CView::OnFilePrintPreview)
  71. END_MESSAGE_MAP()
  73. /////////////////////////////////////////////////////////////////////////////
  74. // CDipView construction/destruction
  76. CDipView::CDipView()
  77. //: Height(NULL)
  78. //, TopRow(NULL)
  79. //: Widths(NULL)
  80. //, Left(NULL)
  81. //: Rect_X(NULL)
  82. //, Rect_Y(NULL)
  83. //, Rect_W(NULL)
  84. //, Rect_H(NULL)
  85. //, Rect_Areas(0)
  86. {
  87. // TODO: add construction code here
  88. //初始化
  89. m_pMemDC=new CDC; //分配内存DC
  90. m_pBitmap=new CBitmap; //分配背景图像的空间
  91. IsNewFile=true;
  92. IsNewBitmap=true;
  93. m_pBmp=NULL;
  94. Height=new int[20];
  95. TopRow=new int[20];
  96. Widths=new int[30];                     //Top' width
  97. Left=new int[30];                   //Top' top col
  98. Rect_X=new int[50];
  99.     Rect_Y=new int[50]; 
  100. Rect_W=new int[50];
  101.     Rect_H=new int[50];
  102. Rect_Areas=0;
  103. cRects=new crect[10];
  104. CharacterNum=0;
  105. cDatas=new cData[10];        //实际大小由CharacterNum确定
  106. char_ConnectPos=NULL;
  107. num_ConnectPos=NULL;
  108. han_ConnectPos=NULL;
  110. }
  112. CDipView::~CDipView()
  113. {
  114. //释放内存
  115. delete m_pBitmap;
  116. delete m_pMemDC;
  117. }
  119. BOOL CDipView::PreCreateWindow(CREATESTRUCT& cs)
  120. {
  121. // TODO: Modify the Window class or styles here by modifying
  122. //  the CREATESTRUCT cs
  123. return CView::PreCreateWindow(cs);
  124. }
  126. /////////////////////////////////////////////////////////////////////////////
  127. // CDipView drawing
  129. void CDipView::OnDraw(CDC* pDC)
  130. {
  131. CDipDoc* pDoc = GetDocument();
  132. ASSERT_VALID(pDoc);
  133. // TODO: add draw code for native data here
  135. //如果已经读入一幅新的图像，则向CDC中写入新的图像信息
  136. if ((pDoc->ImgData!=NULL)&&IsNewFile)
  137. {
  138. //
  139.    ChageStatuBar(&"正在打开文件 ");
  140. //
  141. int i,j,k;
  142. k=0;
  143. BYTE r,g,b;
  144. for(i=0;i<pDoc->ImgHeight;i++)
  145. for(j=0;j<pDoc->ImgWidth;j++)
  146. {
  147. //从调色板中得到r,g,b分量
  148. r=pDoc->palette[pDoc->ImgData[k]].rgbBlue;
  149. g=pDoc->palette[pDoc->ImgData[k]].rgbGreen;
  150. b=pDoc->palette[pDoc->ImgData[k]].rgbRed;
  151. pDC->SetPixel(j,i,(r<<16)|(g<<8)|b);
  152. k++;
  153. }
  154. IsNewFile=false; //写完后，将文件标志位置为假
  155.     //恢复鼠标
  157.  
  158.     char * path =pDoc->FileName.GetBuffer(256) ;
  159. ChageStatuBar(path);
  160. }
  162.        
  163. }
  165. /////////////////////////////////////////////////////////////////////////////
  166. // CDipView printing
  168. BOOL CDipView::OnPreparePrinting(CPrintInfo* pInfo)
  169. {
  170. // default preparation
  171. return DoPreparePrinting(pInfo);
  172. }
  174. void CDipView::OnBeginPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
  175. {
  176. // TODO: add extra initialization before printing
  177. }
  179. void CDipView::OnEndPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)
  180. {
  181. // TODO: add cleanup after printing
  182. }
  184. /////////////////////////////////////////////////////////////////////////////
  185. // CDipView diagnostics
  187. #ifdef _DEBUG
  188. void CDipView::AssertValid() const
  189. {
  190. CView::AssertValid();
  191. }
  192. //setpentext
  193. void CDipView::Dump(CDumpContext& dc) const
  194. {
  195. CView::Dump(dc);
  196. }
  198. CDipDoc* CDipView::GetDocument() // non-debug version is inline
  199. {
  200. ASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CDipDoc)));
  201. return (CDipDoc*)m_pDocument;
  202. }
  203. #endif //_DEBUG
  205. /////////////////////////////////////////////////////////////////////////////
  206. // CDipView message handlers
  208. //打开一个BMP文件
  209. void CDipView::OnFileOpen() 
  210. {
  211. // TODO: Add your command handler code here
  212. CFileDialog mydlg(true);
  213. CString filename;
  214. mydlg.m_ofn.lpstrFilter="BMP files(*.bmp)*.bmpJPG files(*.jpg)*.jpg";
  216. //显示文件对话框
  217. if (mydlg.DoModal()==IDOK)
  218. {
  219. filename=mydlg.GetPathName();
  221. //读入文件
  222. CDipDoc * myDoc;
  223. myDoc=GetDocument();            //enable you to call the number
  224.                                 //function of CDocment 
  225. myDoc->FileOpen(filename);
  227. //将文件标志位和图像标志位置为真
  228. IsNewFile=true;
  229. IsNewBitmap=true;
  231. //刷新屏幕
  232. myDoc->UpdateAllViews(NULL);
  233. }
  234. }
  236. void CDipView::OnPaint() 
  237. {
  238. CPaintDC dc(this); // device context for painting
  239. // TODO: Add your message handler code here
  240. this->GetActiveWindow()->SetWindowText("数字图像处理程序");
  241. //设定从内存DC到屏幕DC的拷贝区域
  242. CRect rc(0,0,::GetSystemMetrics(SM_CXSCREEN),::GetSystemMetrics(SM_CYSCREEN));
  244. CBitmap * pOldBitmap=(CBitmap *)(m_pMemDC->SelectObject(m_pBitmap));
  246. if (IsNewBitmap)
  247. m_pMemDC->Rectangle(-1,-1,2000,2000); //新图，清除背景
  248. else
  249. m_pMemDC->SelectObject(mymap); //旧图，直接调用内存中的图像
  251. m_pMemDC->SelectClipRgn(NULL);
  252. OnDraw(m_pMemDC); //向内存DC中画图像
  253. //拷贝到屏幕DC中
  254. dc.BitBlt(rc.left,rc.top,rc.Width(),rc.Height(),m_pMemDC,rc.left,rc.top,SRCCOPY);
  255. if (IsNewBitmap)
  256. {
  257. //保存历史图像，已备下次显示
  258. mymap=m_pMemDC->GetCurrentBitmap();
  259. IsNewBitmap=false;
  260. }
  262. m_pMemDC->SelectObject(pOldBitmap);
  264. // m_wait=theapp->LoadCursor(IDC_ARROW);
  265. // SetCursor(m_wait);
  266. // Do not call CView::OnPaint() for painting messages
  267. }
  269. void CDipView::OnInitialUpdate() 
  270. {
  271. CView::OnInitialUpdate();
  272. // ChageStatuBar(&"请打开一个BMP格式的位图文件");
  273. // TODO: Add your specialized code here and/or call the base class
  274. //初始化m_pMemDC和m_pBitmap
  275. if (m_pMemDC->GetSafeHdc()==NULL)
  276. {
  277. CClientDC dc(this);
  278. OnPrepareDC(&dc);
  279. //设定背景图像的大小
  280. CRect rc(0,0,::GetSystemMetrics(SM_CXSCREEN),::GetSystemMetrics(SM_CYSCREEN));
  281. m_pMemDC->CreateCompatibleDC(&dc);
  282. m_pBitmap->CreateCompatibleBitmap(&dc,rc.right,rc.bottom);
  283. }
  285. CDipDoc* pDoc=GetDocument();
  286. int i;
  287.    for( i=0;i<pDoc->ImgWidth*pDoc->ImgHeight;i++)
  288.    img[i]=pDoc->ImgData[i];
  289.    
  290. }
  294. void CDipView::OnFileSave() 
  295. {
  296. // TODO: Add your command handler code here
  297. CFileDialog mydlg(false);
  298. CString filename,promt;
  299. mydlg.m_ofn.lpstrFilter="BMP files(*.bmp)*.bmp";
  300. //打开存储对话框
  301. if (mydlg.DoModal()==IDOK)
  302. {
  303. filename=mydlg.GetPathName();
  304. if(filename.Find(".bmp")<=0)
  305. filename=filename+".bmp";
  307. //存储图像
  308. CDipDoc* pDoc = GetDocument();
  309. pDoc->FileSave(filename);
  310. }
  311. }
  313. //重新显示原始图像
  314. void CDipView::Onreopen() 
  315. {
  316. // TODO: Add your command handler code here
  317. CDipDoc* pDoc = GetDocument();
  318. pDoc->FileOpen();
  319. IsNewBitmap=true;
  320. IsNewFile=true;
  321. pDoc->UpdateAllViews(NULL);
  322. }
  324. void CDipView::OnMouseMove(UINT nFlags, CPoint point) 
  325. {
  326. if (AfxGetApp() != NULL && AfxGetApp()->m_pMainWnd != NULL) {
  327. char buffer[256];
  328. CStatusBar* pStatus = (CStatusBar*) 
  329. AfxGetApp()->m_pMainWnd->GetDescendantWindow(AFX_IDW_STATUS_BAR);
  331. va_list argptr;
  332. pStatus->SetPaneInfo(2,2,SBPS_STRETCH  ,7);       
  333. CString x;
  334. x.Format(_T("%d"),point.x); 
  336. CString y;
  337. y.Format(_T("%d"),point.y); 
  339.  CString zuobiao=" x= "+x + " , " + "y= " +y;
  340. char * fmt=zuobiao.GetBuffer(12);
  341. va_start(argptr, fmt);
  342. vsprintf(buffer, fmt, argptr);
  343. va_end(argptr);
  344.         if (pStatus != NULL)
  345. {
  346. pStatus->SetPaneText(2, buffer);
  347. pStatus->UpdateWindow();
  348. }
  349. }
  350. // TODO: Add your message handler code here and/or call default
  351. CView::OnMouseMove(nFlags, point);
  352. }
  354. void CDipView::ChageStatuBar(char * fmt)
  355. {
  357. if (AfxGetApp() != NULL && AfxGetApp()->m_pMainWnd != NULL) {
  358. char buffer[256];
  359. CStatusBar* pStatus = (CStatusBar*) 
  360. AfxGetApp()->m_pMainWnd->GetDescendantWindow(AFX_IDW_STATUS_BAR);
  362. va_list argptr;
  363. va_start(argptr, fmt);
  364. vsprintf(buffer, fmt, argptr);
  365. va_end(argptr);
  366. if (pStatus != NULL) {
  367. pStatus->SetPaneText(0, buffer);
  368. pStatus->UpdateWindow();
  369. }
  370. }
  373. }
  375. //DEL BOOL CDipView::OnEraseBkgnd(CDC* pDC) 
  376. //DEL {
  377. //DEL  // TODO: Add your message handler code here and/or call default
  378. //DEL  if (AfxGetApp() != NULL && AfxGetApp()->m_pMainWnd != NULL) {
  379. //DEL  char buffer[256];
  380. //DEL  CStatusBar* pStatus = (CStatusBar*) 
  381. //DEL  AfxGetApp()->m_pMainWnd->GetDescendantWindow(AFX_IDW_STATUS_BAR);
  382. //DEL 
  383. //DEL  va_list argptr;
  384. //DEL  // char * fmt =&" BMPBMPBMPBMP";
  385. //DEL //gai kuandu
  386. //DEL  pStatus->SetPaneInfo(2,2,SBPS_STRETCH  ,7);
  387. //DEL        //itoa(point.x,,10)
  388. //DEL 
  389. //DEL  char * fmt="请打开文件";
  390. //DEL  va_start(argptr, fmt);
  391. //DEL  vsprintf(buffer, fmt, argptr);
  392. //DEL  va_end(argptr);
  393. //DEL         if (pStatus != NULL)
  394. //DEL  {
  395. //DEL  // pStatus
  396. //DEL  pStatus->SetPaneText(2, buffer);
  397. //DEL  pStatus->UpdateWindow();
  398. //DEL  }
  399. //DEL  }
  400. //DEL  return CView::OnEraseBkgnd(pDC);
  401. //DEL }
  403. int CDipView::OnCreate(LPCREATESTRUCT lpCreateStruct) 
  404. {
  405. if (CView::OnCreate(lpCreateStruct) == -1)
  406. return -1;
  408. OnLoad();
  409. char_bp.Initialize(3434);
  410. char_bp.CreateBP(13,60,26);
  411. num_bp.Initialize(3434);
  412. num_bp.CreateBP(13,60,10);
  413. char_num_bp.Initialize(3434);
  414. char_num_bp.CreateBP(13,60,36);
  417. // TODO: Add your specialized creation code here
  418. // ChageStatuBar(&"请打开一个BMP格式的位图文件");
  419. return 0;
  420. }
  422. void CDipView::OnLine() 
  423. {
  424. // TODO: Add your command handler code here
  425. CDipDoc* pDoc = GetDocument();
  426. ASSERT_VALID(pDoc);
  427. Mdlg dlg;
  428. dlg.DoModal();
  429. a=dlg.m_edit1;
  430. b=dlg.m_edit2;
  431. for(int i=0;i<pDoc->ImgWidth*pDoc->ImgHeight;i++)
  432. {
  433. pDoc->ImgData[i]=a*pDoc->ImgData[i]+b;
  434. if(pDoc->ImgData[i]>255) pDoc->ImgData[i]=255;
  435. if(pDoc->ImgData[i]<0) pDoc->ImgData[i]=0;
  436. }
  438.      IsNewFile=true;
  439.      IsNewBitmap=true;
  440.  
  441. Invalidate();
  445. }
  447. void CDipView::Onfse() 
  448. {
  449. // TODO: Add your command handler code here
  450. CDipDoc* pDoc = GetDocument();
  451. ASSERT_VALID(pDoc);
  452. for(int i=0;i<pDoc->ImgWidth*pDoc->ImgHeight;i++)
  453. {
  454. //pDoc->ImgData[i]=a*pDoc->ImgData[i]+b;
  455.   pDoc->ImgData[i]=255-pDoc->ImgData[i];
  456. //if(pDoc->ImgData[i]<0) pDoc->ImgData[i]=0;
  457. }
  458.  IsNewFile=true;
  459.      IsNewBitmap=true;
  460.  
  461. Invalidate();
  463. }
  465. void CDipView::Onfz() 
  466. {
  467. // TODO: Add your command handler code here
  468. CDipDoc* pDoc = GetDocument();
  470. ASSERT_VALID(pDoc);
  471. BYTE * ImgData1=new BYTE[pDoc->ImgWidth*pDoc->ImgHeight];
  472. for(int i=0;i<pDoc->ImgWidth*pDoc->ImgHeight;i++)
  473. {
  474. //pDoc->ImgData[i]=a*pDoc->ImgData[i]+b;
  475.   ImgData1[pDoc->ImgWidth*pDoc->ImgHeight-i]=pDoc->ImgData[i];
  476. //if(pDoc->ImgData[i]<0) pDoc->ImgData[i]=0;
  477. }
  478. for( i=0;i<pDoc->ImgWidth*pDoc->ImgHeight;i++)
  479. {
  480. //pDoc->ImgData[i]=a*pDoc->ImgData[i]+b;
  481. pDoc->ImgData[i]=ImgData1[i];
  482. //if(pDoc->ImgData[i]<0) pDoc->ImgData[i]=0;
  483. }
  484.  IsNewFile=true;
  485.      IsNewBitmap=true;
  486.  
  487. Invalidate();
  488. }
  490. void CDipView::OnMoHu() 
  491. {
  492. // TODO: Add your command handler code here
  493. CDipDoc* pDoc = GetDocument();
  494. int h=pDoc->ImgHeight;
  495. int w=pDoc->ImgWidth;
  496. BYTE * g=new BYTE[h*w];
  497. float temp;
  498. //标示图象的坐标
  499. int x,y;
  500. //标示算子的坐标
  501. int i,j;
  503. BYTE mb[3][3]={{-1,0,0},{0,0,0},{0,0,1}};
  504.   
  505. for(x=1;x<h-1;x++)
  506. {
  507. for(y=1;y<w-1;y++)
  508. {
  509. temp=0;
  510.     for(i=0;i<3;i++)
  511.      {
  512.      for(j=0;j<3;j++)
  513. {
  514.      temp=temp+(float( mb[i][j])*( float (pDoc->ImgData[(x-1+i)*w+(y-1+j)])));
  515. }
  516. }
  517. g[x*w+y]=int(temp);
  518. }
  519. }
  520. for(x=1;x<h-1;x++)
  521. for(y=1;y<w-1;y++)
  522. pDoc->ImgData[x*w+y]=g[x*w+y];
  524. IsNewFile=true;
  525.     IsNewBitmap=true;  
  526. Invalidate();
  527. }
  529. void CDipView::OnForHandle() 
  530. {
  531. // TODO: Add your command handler code here
  532.  CDipDoc *p=GetDocument();
  533.  BYTE *temp=new BYTE[p->ImgHeight*p->ImgWidth];
  534.  int x,y;
  536. for(x=0;x<p->ImgHeight;x++)
  537. {
  538. for(y=0;y<p->ImgWidth;y++)
  539.  {
  540. temp[(x*p->ImgWidth)+y]=p->ImgData[(x*p->ImgWidth)+y];
  541. }
  542. }
  544.  long double total=0;
  545.  float aver=0;
  546.  for(x=0;x<p->ImgHeight;x++)
  547.  {
  548.  for(y=0;y<p->ImgWidth;y++)
  549.  {
  550.  total+=p->ImgData[x*p->ImgWidth+y];
  551.  }
  552.  }
  553.  aver=(BYTE)(total/(float)(p->ImgHeight*p->ImgWidth));
  555.  float delt; 
  556.  float sub=0; 
  557.  for(x=0;x<p->ImgHeight;x++)
  558.  {
  559.  for(y=0;y<p->ImgWidth;y++)
  560.  {
  561.  sub+=(p->ImgData[x*p->ImgWidth+y]-aver)*(p->ImgData[x*p->ImgWidth+y]-aver);  
  562.  }
  563.  }
  564.  delt=(float)(sqrt(sub/(float)(p->ImgHeight*p->ImgWidth)));
  566.  BYTE judge;
  567.  judge=(BYTE)(delt+aver);
  568.  for(x=0;x<p->ImgHeight;x++)
  569.  {
  570. for(y=0;y<p->ImgWidth;y++)
  571.  {
  572. if(p->ImgData[(x*p->ImgWidth)+y]>=judge)
  573. {
  574. p->ImgData[(x*p->ImgWidth)+y]=255;
  575. }
  576. else
  577. {
  578. p->ImgData[(x*p->ImgWidth)+y]=0;
  579. }
  580. }
  581. }
  582.  
  583.  free(temp);
  584.  IsNewBitmap=true;
  585.  IsNewFile=true;
  586.  p->UpdateAllViews(NULL);
  587. }
  589. void CDipView::OnHSearch() 
  590. {
  591. // TODO: Add your command handler code here
  592. CDipDoc *pDoc=GetDocument();
  594. int width,height;
  595. width=pDoc->ImgWidth;
  596. height=pDoc->ImgHeight;
  598. BYTE *g=new BYTE[width*height];
  600.     int *Count= new int[height];
  601. int *TempCount=new int[height];
  602. int *maxChanges=new int[10];
  603. int *Rows=new int[10];                                  
  605. int H[20];
  606. int Top[20];
  607.  
  608. for(int row=0;row<pDoc->ImgHeight;row++)
  609.  {
  610. for(int col=0;col<pDoc->ImgWidth;col++)
  611. {
  612.    g[row*width+col]=pDoc->ImgData[row*width+col];
  613. }
  614.  }
  615.  
  616. for(int c=0;c<height;c++)
  617. {
  618. Count[c]=0;
  619. TempCount[c]=0;
  620. }
  622. for(c=0;c<10;c++)
  623. {
  624.         maxChanges[c]=0;
  625.     Rows[c]=0; 
  626. H[c]=0;
  627. Top[c]=1;
  628. }
  630.     for(row=0;row<height;row++)
  631. {
  632. for(int col=0;col<width-1;col++)
  633. {
  634. if(pDoc->ImgData[row*width+col]!=pDoc->ImgData[row*width+col+1])
  635. Count[row]++;
  636. }
  637. TempCount[row]=Count[row];
  638. }
  640. /////////////////////////////////////Top 10 chage row
  641. for(int top=0;top<10;top++)
  642. {
  643. int Max=0;
  644. for(row=0;row<height;)
  645. {
  646. if(Max<TempCount[row])
  647. {
  648. Max=TempCount[row];
  649. }
  650. row=row+10 ;
  651. }
  652. maxChanges[top]=Max;
  654. int R=0;
  655.      for(row=0;row<height;)
  656. {
  657.      if(Max!=TempCount[row])
  658. {
  659.      row=row+10;
  660. }
  661. else
  662. {
  663. R=row;
  664. break;
  665. }
  666. }
  667. Rows[top]=R;
  668. TempCount[R]=0;
  669. }
  671. for(top=0;top<10;top++)
  672. {
  673. for(int range=-10;range<10;range++)
  674. {
  675. if((0.7*Count[Rows[top]]<=Count[Rows[top]+range])&&(Count[Rows[top]+range]<=1.3*Count[Rows[top]]))
  676. {
  677.  for(int col=0;col<width;col++)
  678.  {
  679.          g[(Rows[top]+range)*width+col]=255;
  680.  }
  681. }
  682. }
  683. }
  685.  for(row=0;row<height;row++)
  686.  {
  687. for(int col=0;col<pDoc->ImgWidth;col++)
  688. {
  689.     if(g[row*width+col]<255)
  690. {
  691.                g[row*width+col]=0;
  692. }
  693. }
  694.  }
  696. int Area=0;                   //标识区域的地点
  697. if(g[0*width+10]==255)
  698. {
  699. Area++;
  700. H[Area]++;
  701.     Top[Area]=0;
  702. }
  703. for(row=1;row<height;row++)
  704. {
  705. if((g[row*width+10]==255) && (g[(row-1)*width+10]==0))
  706. {
  707. Area++;
  708.         H[Area]++;
  709.      Top[Area]=row;
  710. }
  711. else 
  712. {
  713. if(g[row*width+10]==255&&g[(row-1)*width+10]==255)
  714.     {
  715.      H[Area]++;
  716.     }     
  717. }
  718.  }  
  720. for(int x=0;x<20;x++)
  721. {
  722. int h=H[x];
  723. if(h>20)                                //特定应用条件，拍摄所得车牌的区域高度、宽度均有一定的范围
  724. {                                       //分析图象可得
  725.     int r=0;
  726. int r0=0;
  727. r=Top[x];
  728. r0=Top[x]+h;
  730. if(r!=0&&r0!=height)
  731. {
  732. for(int col=0;col<width;col++)
  733.     {
  734. pDoc->ImgData[r*width+col]=255;
  735. pDoc->ImgData[r0*width+col]=255;
  736. }
  737. }
  738. else
  739. {
  740. if(r==0)
  741. {
  742. r=5;
  743. for(int col=0;col<width;col++)
  744.         {
  745. pDoc->ImgData[r*width+col]=255;
  746. pDoc->ImgData[r0*width+col]=255;
  747. }
  748. }
  749. else
  750. {
  751. if(r0==height)
  752. {
  753. r0=height-10;
  754. for(int col=0;col<width;col++)
  755. {
  756. pDoc->ImgData[r*width+col]=255;
  757. pDoc->ImgData[r0*width+col]=255;
  758. }
  759. }
  760. }
  761. }
  762. }
  763. }
  765. free(g);
  766. IsNewBitmap=true;
  767. IsNewFile=true;
  768. pDoc->UpdateAllViews(NULL);
  769. }
  771. void CDipView::OnBanlance() 
  772. {
  773. // TODO: Add your command handler code here
  774. CDipDoc *pDoc=GetDocument();
  776.  int width,height;
  777.  width=pDoc->ImgWidth;
  778.  height=pDoc->ImgHeight;
  779.  float num=0;
  781.  int row,col;
  782.  int count[256];
  783.  float s0[256];
  784.  int s[256];                        //256个灰度等级
  785.  int x=0;
  786.  int temp=0;
  788.  for( x=0;x<256;x++)
  789.  {
  790.  count[x]=0;
  791.  s0[x]=0.0;
  792.  s[x]=0;
  793.  }
  795.  for(row=0;row<pDoc->ImgHeight;row++)                          //每个灰度等级象素个数
  796.  {
  797. for(col=0;col<pDoc->ImgWidth;col++)
  798. {
  799. temp=pDoc->ImgData[row*width+col];
  800. count[temp]++;
  801. }
  802.  }
  804.  for(x=0;x<256;x++)                                       //对应灰度等级概率
  805.  {
  806.  num=(float)width*height;
  807.  s0[x]=(float)(count[x]/num);
  808.  }
  810.  for( x=1;x<256;x++)                                      //概率累计，合并灰度等级    
  811.  {
  812.  s0[x]+=s0[x-1];  
  813.  }
  815.      for( x=0;x<256;x++)                                      //映射
  816.  {
  817.   s[x]=(int)(s0[x]*255+0.5);
  818.  }
  820.   for(row=0;row<pDoc->ImgHeight;row++)                          //每个灰度等级象素个数
  821.  {
  822. for(col=0;col<pDoc->ImgWidth;col++)
  823. {
  824. int reflect=0;
  825. temp=pDoc->ImgData[row*width+col];
  826. reflect=s[temp];
  827. pDoc->ImgData[row*width+col]=reflect;
  828. }
  829.  }
  831.  IsNewBitmap=true;
  832.  IsNewFile=true;
  833.  pDoc->UpdateAllViews(NULL); 
  834. }
  836. void CDipView::OnSubVertical() 
  837. {
  838. // TODO: Add your command handler code here
  839. CDipDoc *pDoc=GetDocument();
  841. int width,height;
  842. width=pDoc->ImgWidth;
  843. height=pDoc->ImgHeight;
  845.     BYTE *g=new BYTE[width*height];
  847. int row,col,sub;
  849. for(row=0;row<height-1;row++)
  850. {
  851. for(col=0;col<width;col++)
  852. {
  853. sub=abs((pDoc->ImgData[row*width+col])-(pDoc->ImgData[(row+1)*width+col]));
  854. g[row*width+col]=sub;
  855. }
  856. }
  858. for(row=0;row<height;row++)
  859. {
  860. for(col=0;col<width-1;col++)
  861. {
  862. pDoc->ImgData[row*width+col]=g[row*width+col];
  863. }
  864. }
  865.  
  866. free(g);
  867. IsNewBitmap=true;
  868. IsNewFile=true;
  869. pDoc->UpdateAllViews(NULL);
  871. }
  873. void CDipView::OnSubHorizontal() 
  874. {
  875. // TODO: Add your command handler code here
  876. CDipDoc *pDoc=GetDocument();
  878. int width,height;
  879. width=pDoc->ImgWidth;
  880. height=pDoc->ImgHeight;
  882.     BYTE *g=new BYTE[width*height];
  884. int row,col,sub;
  886. for(row=0;row<height;row++)
  887. {
  888. for(col=0;col<width-1;col++)
  889. {
  890. sub=abs((pDoc->ImgData[row*width+col])-(pDoc->ImgData[row*width+col+1]));
  891. if(sub<15)
  892. {
  893. g[row*width+col]=0;
  894. }
  895. else
  896. {
  897. if(sub>15)
  898. {
  899. g[row*width+col]=sub+50;
  900. }
  901. // else
  902. // g[row*width+col]=sub;
  903. }
  904. }
  905. }
  907. for(row=0;row<height;row++)
  908. {
  909. for(col=0;col<width-1;col++)
  910. {
  911. pDoc->ImgData[row*width+col]=g[row*width+col];
  912. if(col==0)
  913. {
  914. pDoc->ImgData[row*width+col]=0;
  915. }
  916. }
  917. }
  918.  
  919. free(g);
  920. IsNewBitmap=true;
  921. IsNewFile=true;
  922. pDoc->UpdateAllViews(NULL);
  923. }
  926. void CDipView::OnMidValue() 
  927. {
  928. // TODO: Add your command handler code here
  929. CDipDoc *pDoc=GetDocument();
  931. int width,height;
  932. width=pDoc->ImgWidth;
  933. height=pDoc->ImgHeight;
  935.     BYTE *g=new BYTE[width*height];
  936. BYTE temp[25];
  938. int row,col;
  939. int u,v,a;
  941. for(row=2;row<height-2;row++)
  942. {
  943. for(col=2;col<width-2;col++)
  944. {
  945. for(u=0;u<5;u++)
  946. {
  947. for(v=0;v<5;v++)
  948. {
  949. temp[u*5+v]=pDoc->ImgData[(row-2+u)*width+(col-2)+v];
  950. }
  951. }
  953. for(u=0;u<25;u++)
  954. {
  955. for(v=u+1;v<25;v++)
  956. {
  957. if(temp[u]>temp[v])
  958. {
  959. a=temp[u];
  960. temp[u]=temp[v];
  961. temp[v]=a;
  962. }
  963. }
  964. }
  966. g[row*width+col]=temp[12];
  967. }
  968. }
  970. for(row=0;row<height;row++)
  971. {
  972. for(col=0;col<width-1;col++)
  973. {
  974. pDoc->ImgData[row*width+col]=g[row*width+col];
  975. }
  976. }
  978. IsNewBitmap=true;
  979. IsNewFile=true;
  980. pDoc->UpdateAllViews(NULL);
  981. }
  983. void CDipView::OnAverage() 
  984. {
  985. // TODO: Add your command handler code here
  986. CDipDoc *pDoc=GetDocument();
  988.  int width,height;
  989.  width=pDoc->ImgWidth;
  990.  height=pDoc->ImgHeight;
  992.  BYTE *g=new BYTE[width*height];
  994.  int row,col;
  995.  BYTE s[3][3]={{1,1,1},{1,1,1},{1,1,1}};
  997.  for(row=1;row<pDoc->ImgHeight-1;row++)
  998.  {
  999. for(col=1;col<pDoc->ImgWidth-1;col++)
  1000. {
  1001. long sum=0;
  1002. for(int i=0;i<3;i++)
  1003. {
  1004. for(int j=0;j<3;j++)
  1005. {
  1006. sum+=s[i][j]*pDoc->ImgData[((row-1+i)*pDoc->ImgWidth)+(col-1+j)];
  1007. }
  1008. }
  1009. g[(row*pDoc->ImgWidth)+col]=(BYTE)(sum/9);
  1010. }
  1011.  }
  1013.  for(row=0;row<height;row++)
  1014.  {
  1015. for(col=0;col<width;col++)
  1016. {
  1017. pDoc->ImgData[row*width+col]=g[row*width+col];
  1018. }
  1019.      }
  1021.  IsNewBitmap=true;
  1022.  IsNewFile=true;
  1023.  pDoc->UpdateAllViews(NULL); 
  1024. }
  1026. void CDipView::OnBool() 
  1027. {
  1028. // TODO: Add your command handler code here
  1029.  CDipDoc *pDoc=GetDocument();
  1031.  int width,height;
  1032.  width=pDoc->ImgWidth;
  1033.  height=pDoc->ImgHeight;
  1035.  int row,col;
  1037.  for(row=0;row<pDoc->ImgHeight;row++)
  1038.  {
  1039. for(col=0;col<pDoc->ImgWidth;col++)
  1040. {
  1041.     if(pDoc->ImgData[row*width+col]<255)
  1042. {
  1043.                pDoc->ImgData[row*width+col]=0;
  1044. }
  1045. }
  1046.  }
  1047.  
  1048.  IsNewBitmap=true;
  1049.  IsNewFile=true;
  1050.  pDoc->UpdateAllViews(NULL); 
  1051. }
  1053. void CDipView::OnVSearch()
  1054. {
  1055. // TODO: 在此添加命令处理程序代码
  1057. CDipDoc *pDoc=GetDocument();
  1059. int width,height;
  1060. width=pDoc->ImgWidth;
  1061. height=pDoc->ImgHeight;
  1063. BYTE *g=new BYTE[width*height];
  1065.     int *Count= new int[width];
  1066. int *TempCount=new int[width];
  1067. int *maxChanges=new int[30];
  1068. int *Cols=new int[30];                  //maxChanges' col                       
  1070. int *W=new int[30];                     //Top' width
  1071. int *left=new int[30];                   //Top' top col
  1072.  
  1073. for(int row=0;row<pDoc->ImgHeight;row++)
  1074.  {
  1075. for(int col=0;col<pDoc->ImgWidth;col++)
  1076. {
  1077.    g[row*width+col]=pDoc->ImgData[row*width+col];
  1078. }
  1079.  }
  1080.  
  1081. for(int c=0;c<width;c++)
  1082. {
  1083. Count[c]=0;
  1084. TempCount[c]=0;
  1085. }
  1087. for(c=0;c<30;c++)
  1088. {
  1089.         maxChanges[c]=0;
  1090.     Cols[c]=0; 
  1091. W[c]=0;
  1092. left[c]=-1;
  1093. }
  1095.     for(int col=0;col<width;col++)
  1096. {
  1097. for(int row=0;row<height-1;row++)
  1098. {
  1099. if(pDoc->ImgData[row*width+col]!=pDoc->ImgData[(row+1)*width+col])
  1100. Count[col]++;
  1101. }
  1102. TempCount[col]=Count[col];
  1103. }
  1105. /////////////////////////////////////Top 30 change col
  1106. for(int top=0;top<30;top++)
  1107. {
  1108. int Max=0;
  1109. for(int col=0;col<width;)
  1110. {
  1111. if(Max<TempCount[col])
  1112. {
  1113. Max=TempCount[col];
  1114. }
  1115. col=col+10 ;
  1116. }
  1117. maxChanges[top]=Max;
  1119. int C=0;
  1120.      for(col=0;col<width;)
  1121. {
  1122.      if(Max!=TempCount[col])
  1123. {
  1124.      col=col+10;
  1125. }
  1126. else
  1127. {
  1128. C=col;
  1129. break;
  1130. }
  1131. }
  1132. Cols[top]=C;
  1133. TempCount[C]=0;
  1134. }
  1136. for(top=0;top<30;top++)
  1137. {
  1138. for(int range=-20;range<20;range++)
  1139. {
  1140. if((Cols[top]+range)>=0&&(Cols[top]+range)<width)
  1141. {
  1142. if((0.8*Count[Cols[top]]<=Count[Cols[top]+range])&&(Count[Cols[top]+range]<=1.2*Count[Cols[top]]))
  1143. {
  1144. for(int row=0;row<height;row++)
  1145. {
  1146.          g[row*width+(Cols[top]+range)]=255;
  1147. }
  1148. }
  1149. }
  1150. else
  1151. {
  1152. continue;
  1153. }
  1155. }
  1156. }
  1158.  for(row=0;row<height;row++)
  1159.  {
  1160. for(int col=0;col<width;col++)
  1161. {
  1162.     if(g[row*width+col]<255)
  1163. {
  1164.                g[row*width+col]=0;
  1165. }
  1166. }
  1167.  }
  1169. int Area=-1;                   //标识区域的地点
  1170. if(g[10*width+0]==255)
  1171. {
  1172. Area++;
  1173. W[Area]++;
  1174.     left[Area]=0;
  1175. }
  1176. for(col=1;col<width;col++)
  1177. {
  1178. if(g[10*width+col]==255 && g[10*width+(col-1)]==0)
  1179. {
  1180. Area++;
  1181. W[Area]++;
  1182.      left[Area]=col;
  1183. }
  1184. else 
  1185. {
  1186. if((g[10*width+col]==255)&&(g[10*width+col-1]==255))
  1187. {
  1188.      W[Area]++;
  1189. }     
  1190. }
  1191.  }
  1193. for(int x=0;x<30;x++)
  1194. {
  1195. int w=W[x];
  1196. if(w>10)
  1197. {
  1198.     int c=0;
  1199. int c0=0;
  1200. c=left[x];
  1201. c0=left[x]+w;
  1203. if(c!=0&&c0!=width)
  1204. {
  1205. for(int row=0;row<height;row++)
  1206. {
  1207. pDoc->ImgData[row*width+c]=255;
  1208. pDoc->ImgData[row*width+c0]=255;
  1209. }
  1210. }
  1211. else
  1212. {
  1213. if(c==0)
  1214. {
  1215. c=10;
  1216. for(int row=0;row<height;row++)
  1217. {
  1218. pDoc->ImgData[row*width+c]=255;
  1219. pDoc->ImgData[row*width+c0]=255;
  1220. }
  1221. }
  1222. else
  1223. {
  1224. if(c0==width)
  1225. {
  1226. c0=width-10;
  1227. for(int row=0;row<height;row++)
  1228.     {
  1229. pDoc->ImgData[row*width+c]=255;
  1230. pDoc->ImgData[row*width+c0]=255;
  1231. }
  1232. }
  1233. }
  1234. }
  1235. }
  1236. }
  1238. free(g);
  1239. IsNewBitmap=true;
  1240. IsNewFile=true;
  1241. pDoc->UpdateAllViews(NULL);
  1242. }
  1244. /// <summary>
  1245. /// CDipView::OnLocate()的摘要说明。
  1246. /// 创建人  ：刘陈炜
  1247. /// 创建时间：2004.5.1
  1248. /// 函数功能：完成车牌定位的整个过程
  1249. /// 程序流程：1.产生副本
  1250. ///           2.水平差分提取边缘，寻找横向的车牌带状区域
  1251. ///           3.垂直差分提取边缘，寻找纵向的车牌带状区域
  1252. ///           4.利用先验知识标识车牌区域，进行车牌带状区域的选择，（横纵向）
  1253. ///           5.构造出车牌矩形域；
  1254. ///           6.再次利用利用跳变数选择车牌矩形区域，进行车牌区域的最终选择；
  1255. ///           7.精确定位车牌,即削弱车牌横向宽度；
  1256. ///           8.更新图象；
  1257. /// </summary> 
  1258. void CDipView::OnLocate()
  1259. {
  1260. // TODO: 在此添加命令处理程序代码
  1261. CDipDoc *pDoc=GetDocument();
  1263. int width,height;
  1264. width=pDoc->ImgWidth;
  1265. height=pDoc->ImgHeight;
  1267. BYTE *g=new BYTE[width*height];
  1269. CopyImg(pDoc,g,width,height);
  1270. //Banlance(g,width,height);
  1271. HSub(g,width,height);
  1272.     Search(g,width,height);
  1274. CopyImg(pDoc,g,width,height);
  1275. VSub(g,width,height);
  1276. VSearch(g,width,height);
  1278. Mark(pDoc,width,height);
  1279.     VMark(pDoc,width,height);
  1280. DrawRect(pDoc,width,height);
  1281. Choose(pDoc,width,height);
  1282. FurtherLocate(pDoc,width,height);
  1284. free(g);
  1285. IsNewBitmap=true;
  1286. IsNewFile=true;
  1287. pDoc->UpdateAllViews(NULL);    
  1288. }
  1290. // 均衡化图象
  1291. void CDipView::Banlance(BYTE* DisposeImg,int Width,int Height)
  1292. {
  1293.     float num=0;
  1294. int row,col;
  1295. int count[256];
  1296. float s0[256];
  1297. int s[256];                        //256个灰度等级
  1298. int x=0;
  1299. int temp=0;
  1301. for( x=0;x<256;x++)
  1302. {
  1303. count[x]=0;
  1304. s0[x]=0.0;
  1305. s[x]=0;
  1306. }
  1308. for(row=0;row<Height;row++)                          //每个灰度等级象素个数
  1309. {
  1310. for(col=0;col<Width;col++)
  1311. {
  1312. temp=DisposeImg[row*Width+col];
  1313. count[temp]++;
  1314. }
  1315. }
  1317. for(x=0;x<256;x++)                                       //对应灰度等级概率
  1318. {
  1319. num=(float)Width*Height;
  1320. s0[x]=(float)(count[x]/num);
  1321. }
  1323. for( x=1;x<256;x++)                                      //概率累计，合并灰度等级    
  1324. {
  1325. s0[x]+=s0[x-1];  
  1326. }
  1328.     for( x=0;x<256;x++)                                      //映射
  1329. {
  1330.  s[x]=(int)(s0[x]*255+0.5);
  1331. }
  1333. for(row=0;row<Height;row++)                          //每个灰度等级象素个数
  1334. {
  1335. for(col=0;col<Width;col++)
  1336. {
  1337. int reflect=0;
  1338. temp=DisposeImg[row*Width+col];
  1339. reflect=s[temp];
  1340. DisposeImg[row*Width+col]=reflect;
  1341. }
  1342. }
  1343. }
  1346. // 图象复制
  1347. void CDipView::CopyImg(CDipDoc * SourceImg, BYTE* DestImg, int Width, int Height)
  1348. {
  1349. for(int row=0;row<Height;row++)
  1350. {
  1351. for(int col=0;col<Width;col++)
  1352. {
  1353. DestImg[row*Width+col]=SourceImg->ImgData[row*Width+col];
  1354. }
  1355. }
  1356. }
  1358. // 图象复制
  1359. void CDipView::CopyImg(BYTE* SourceImg, BYTE* DestImg, int Width , int Height)
  1360. {
  1361.     for(int row=0;row<Height;row++)
  1362. {
  1363. for(int col=0;col<Width;col++)
  1364. {
  1365. DestImg[row*Width+col]=SourceImg[row*Width+col];
  1366. }
  1367. }
  1368. }
  1370. /// <summary>
  1371. /// CDipView::HSub(BYTE* DisposeImg, int Width , int Height)的摘要说明。
  1372. /// 创建人  ：刘陈炜
  1373. /// 创建时间：2004.5.1
  1374. /// 函数功能：水平差分提取图象边缘。
  1375. /// 程序流程：1.逐行扫描，sub=|f(x,y)-f(x,y+1)|；
  1376. ///           2.若差值小于20，则g(x,y)=0；
  1377. ///             若差值大于30，则增强原图象，令g(x,y)=g(x,y)+30；
  1378. ///           3.另最后一列全部为0，即黑色；
  1379. /// </summary> 
  1380. // 水平差分
  1381. void CDipView::HSub(BYTE* DisposeImg, int Width , int Height)
  1382. {
  1383. int row,col,sub;
  1385. for(row=0;row<Height;row++)
  1386. {
  1387. for(col=0;col<Width-1;col++)
  1388. {
  1389. sub=abs((DisposeImg[row*Width+col])-(DisposeImg[row*Width+col+1]));
  1390. if(sub<20)
  1391. {
  1392. DisposeImg[row*Width+col]=0;
  1393. }
  1394. else
  1395. {
  1396. if(sub>30)
  1397. {
  1398.                  DisposeImg[row*Width+col]=sub+30;
  1399. }
  1400. else
  1401. DisposeImg[row*Width+col]=sub;
  1402. }
  1403. }
  1404. }
  1406. for(row=0;row<Height;row++)
  1407. {
  1408. DisposeImg[row*Width+Width-1]=0;
  1409. }  
  1410. }
  1412. // 纵向差分
  1413. void CDipView::VSub(BYTE* DisposeImg, int Width , int Height)
  1414. {
  1415. int row,col,sub;
  1417. for(row=0;row<Height-1;row++)
  1418. {
  1419. for(col=0;col<Width;col++)
  1420. {
  1421. sub=abs((DisposeImg[row*Width+col])-(DisposeImg[(row+1)*Width+col]));
  1422. /* if(sub<5)
  1423. {
  1424. DisposeImg[row*Width+col]=0;
  1425. }
  1426. else
  1427. {
  1428. if(sub>10)
  1429. {
  1430.                  DisposeImg[row*Width+col]=sub+30;
  1431. }
  1432. else
  1433. DisposeImg[row*Width+col]=sub;
  1434. }*/
  1435. DisposeImg[row*Width+col]=sub;
  1436. }
  1437. }
  1439. for(col=0;col<Width-1;col++)
  1440. {
  1441. DisposeImg[(Height-1)*Width+col]=0;
  1442. }
  1443. }
  1445. /// <summary>
  1446. /// CDipView::Search(BYTE* DisposeImg, int width , int height)的摘要说明。
  1447. /// 创建人  ：刘陈炜
  1448. /// 创建时间：2004.5.1
  1449. /// 函数功能：完成图象车牌区域的初步定位。即寻找出差分扫描后水平跳变数最大的前10个区域。
  1450. /// 程序流程：1.利用间行扫描的方式找出跳变数最大的前10行；
  1451. ///           2.对 Top 10 行进行区域细粒度增长选择；
  1452. ///             其相似度的选取为：0.7~1.3 倍的该行跳变总数；
  1453. ///           3.二值化，threshold=255；
  1454. ///           4.任意选择一列，最好不要边缘列，可能有干扰（本程序选择为第十列），搜索出各个带状候选域，存入标志；
  1455. /// </summary> 
  1456. void CDipView::Search(BYTE* DisposeImg, int width , int height)
  1457. {
  1458. BYTE*TempImg=new BYTE[width*height];
  1460.     int *Count= new int[height];
  1461. int *TempCount=new int[height];
  1462. int *maxChanges=new int[10];
  1463. int *Rows=new int[10]; 
  1465. CopyImg(DisposeImg,TempImg,width,height);
  1467. for(int c=0;c<height;c++)
  1468. {
  1469. Count[c]=0;
  1470. TempCount[c]=0;
  1471. }
  1473. for(c=0;c<10;c++)
  1474. {
  1475.         maxChanges[c]=0;
  1476.     Rows[c]=0; 
  1477. Height[c]=0;
  1478. TopRow[c]=-1;
  1479. }
  1481.     for(int row=0;row<height;row++)
  1482. {
  1483. for(int col=0;col<width-1;col++)
  1484. {
  1485. if(DisposeImg[row*width+col]!=DisposeImg[row*width+col+1])
  1486. Count[row]++;
  1487. }
  1488. TempCount[row]=Count[row];
  1489. }
  1491. /////////////////////////////////////Top 10 chage row
  1492. for(int top=0;top<10;top++)
  1493. {
  1494. int Max=0;
  1495. for(row=0;row<height;)
  1496. {
  1497. if(Max<TempCount[row])
  1498. {
  1499. Max=TempCount[row];
  1500. }
  1501. row=row+10 ;
  1502. }
  1503. maxChanges[top]=Max;
  1505. int R=0;
  1506.      for(row=0;row<height;)
  1507. {
  1508.      if(Max!=TempCount[row])
  1509. {
  1510.      row=row+10;
  1511. }
  1512. else
  1513. {
  1514. R=row;
  1515. break;
  1516. }
  1517. }
  1518. Rows[top]=R;
  1519. TempCount[R]=0;
  1520. }
  1522. for(top=0;top<10;top++)
  1523. {
  1524. for(int range=-10;range<10;range++)
  1525. {
  1526. if((0.7*Count[Rows[top]]<=Count[Rows[top]+range])&&(Count[Rows[top]+range]<=1.3*Count[Rows[top]]))
  1527. {
  1528.  for(int col=0;col<width;col++)
  1529.  {
  1530.          DisposeImg[(Rows[top]+range)*width+col]=255;
  1531.  }
  1532. }
  1533. }
  1534. }
  1536.  for(row=0;row<height;row++)
  1537.  {
  1538. for(int col=0;col<width;col++)
  1539. {
  1540.     if(DisposeImg[row*width+col]<255)
  1541. {
  1542.                DisposeImg[row*width+col]=0;
  1543. }
  1544. }
  1545.  }
  1547. int Area=0;                   //标识区域的地点
  1548. if(DisposeImg[0*width+10]==255)
  1549. {
  1550. Area++;
  1551. Height[Area]++;
  1552.     TopRow[Area]=0;
  1553. }
  1554. for(row=1;row<height;row++)
  1555. {
  1556. if((DisposeImg[row*width+10]==255) && (DisposeImg[(row-1)*width+10]==0))
  1557. {
  1558. Area++;
  1559.         Height[Area]++;
  1560.      TopRow[Area]=row;
  1561. }
  1562. else 
  1563. {
  1564. if(DisposeImg[row*width+10]==255&&DisposeImg[(row-1)*width+10]==255)
  1565.     {
  1566.      Height[Area]++;
  1567.     }     
  1568. }
  1569. }
  1571. // Mark(TempImg,width,height);
  1573. // CopyImg(TempImg,DisposeImg,width,height);
  1575. free(TempImg);
  1576.     free(Count);
  1577. free(TempCount);
  1578. free(maxChanges);
  1579. free(Rows); 
  1580. }
  1582. /// <summary>
  1583. /// CDipView::Mark(CDipDoc * DisposeImg, int width , int height)的摘要说明。
  1584. /// 创建人  ：刘陈炜
  1585. /// 创建时间：2004.5.1
  1586. /// 函数功能：利用先验知识标识车牌区域，进行车牌区域的选择，（横向）
  1587. /// 函数说明：特定应用条件，拍摄所得车牌的区域高度、宽度均有一定的范围，其大小可以分析图象后可得
  1588. ///           根据统计得到：车牌的顶部极少出现在高度小于100个像素的区域；
  1589. ///                         车牌高度也一般大于20个像素。
  1590. ///           程序的最终结果为满足条件的经过扩展了的带状区域
  1591. /// 程序流程：1.舍弃区域顶部小于100的待选域
  1592. ///           2.选择区域高度大于20的待选域
  1593. ///           3.区域扩展，上下各增大5个像素。若顶部小于5，则顶部坐标不变，高度加5；
  1594. ///                                          若底部大于图象的底部，则底部坐标不变，高度加5；
  1595. /// </summary> 
  1596. void CDipView::Mark(CDipDoc * DisposeImg, int width , int height)
  1597. {
  1598. for(int x=0;x<20;x++)
  1599. {
  1600. int h=Height[x];
  1601. if(TopRow[x]<=5)
  1602. {
  1603. Height[x]=0;
  1604. TopRow[x]=-1;
  1605. }
  1606. }
  1607.  for(x=0;x<20;x++)
  1608. {
  1609. int h=Height[x];
  1610. if(h>20)                                
  1611. {                                       
  1612.     int r=0;
  1613. int r0=0;
  1614. r=TopRow[x]-5;
  1615. r0=TopRow[x]+h+5;
  1617. if((r>0)&&(r0<height))
  1618. {
  1619. // for(int col=0;col<width;col++)
  1620. //     {
  1621. // DisposeImg->ImgData[r*width+col]=255;
  1622. // DisposeImg->ImgData[r0*width+col]=255;
  1623. // }
  1625. TopRow[x]=r+5;
  1626. Height[x]=h;
  1627. }
  1628. else
  1629. {
  1630. if(r<=0)
  1631. {
  1632. r=r+5;
  1634. //for(int col=0;col<width;col++)
  1635.   //      {
  1636. // DisposeImg->ImgData[r*width+col]=255;
  1637. // DisposeImg->ImgData[r0*width+col]=255;
  1638. //}
  1639. TopRow[x]=r;
  1640. Height[x]=h;
  1641. }
  1642. else
  1643. {
  1644. if(r0>=height)
  1645. {
  1646. r0=r0-5;
  1647. //for(int col=0;col<width;col++)
  1648. //{
  1649. // DisposeImg->ImgData[r*width+col]=255;
  1650. // DisposeImg->ImgData[r0*width+col]=255;
  1651. //}
  1652. Height[x]=Height[x]-5;
  1653. }
  1654. }
  1655. }
  1656. }
  1657. }
  1659. for(x=0;x<20;x++)
  1660. {
  1661. int h=Height[x];
  1662. if(h<20||TopRow[x]<=100)
  1663. {
  1664. Height[x]=0;
  1665. TopRow[x]=-1;
  1666. }
  1667. }
  1668. }
  1670. /// <summary>
  1671. /// CDipView::VSearch(BYTE* DisposeImg, int width , int height)的摘要说明。
  1672. /// 创建人  ：刘陈炜
  1673. /// 创建时间：2004.5.1
  1674. /// 函数功能：水平查找后，纵向查找。完成图象车牌区域的初步定位。即寻找出差分扫描后垂直跳变数最大的前30个区域。
  1675. /// 程序流程：1.利用间行扫描的方式找出跳变数最大的前30行；
  1676. ///           2.对 Top 30 行进行区域增长；
  1677. ///             其相似度的选取为：0.8~1.2 倍的该行跳变总数；
  1678. ///           3.二值化，threshold=255；
  1679. ///           4.任意选择一列，（本程序选择为第十行），搜索出各个带状候选域，存入标志；
  1680. /// </summary> 
  1681. //
  1682. void CDipView::VSearch(BYTE* DisposeImg, int width, int height)
  1683. {
  1684. BYTE *g=new BYTE[width*height];
  1686.     int *Count= new int[width];
  1687. int *TempCount=new int[width];
  1688. int *maxChanges=new int[30];
  1689. int *Cols=new int[30];                  //maxChanges' col  
  1691. CopyImg(DisposeImg,g,width,height);
  1692.  
  1693. for(int c=0;c<width;c++)
  1694. {
  1695. Count[c]=0;
  1696. TempCount[c]=0;
  1697. }
  1699. for(c=0;c<30;c++)
  1700. {
  1701.         maxChanges[c]=0;
  1702.     Cols[c]=0; 
  1703. Widths[c]=0;
  1704. Left[c]=-1;
  1705. }
  1707.     for(int col=0;col<width;col++)
  1708. {
  1709. for(int row=0;row<height-1;row++)
  1710. {
  1711. if(DisposeImg[row*width+col]!=DisposeImg[(row+1)*width+col])
  1712. Count[col]++;
  1713. }
  1714. TempCount[col]=Count[col];
  1715. }
  1717. /////////////////////////////////////Top 30 change col
  1718. for(int top=0;top<30;top++)
  1719. {
  1720. int Max=0;
  1721. for(int col=0;col<width;)
  1722. {
  1723. if(Max<TempCount[col])
  1724. {
  1725. Max=TempCount[col];
  1726. }
  1727. col=col+10 ;
  1728. }
  1729. maxChanges[top]=Max;
  1731. int C=0;
  1732.      for(col=0;col<width;)
  1733. {
  1734.      if(Max!=TempCount[col])
  1735. {
  1736.      col=col+10;
  1737. }
  1738. else
  1739. {
  1740. C=col;
  1741. break;
  1742. }
  1743. }
  1744. Cols[top]=C;
  1745. TempCount[C]=0;
  1746. }
  1748. for(top=0;top<30;top++)
  1749. {
  1750. for(int range=-20;range<20;range++)
  1751. {
  1752. if((Cols[top]+range)>=0&&(Cols[top]+range)<width)
  1753. {
  1754. if((0.8*Count[Cols[top]]<=Count[Cols[top]+range])&&(Count[Cols[top]+range]<=1.2*Count[Cols[top]]))
  1755. {
  1756. for(int row=0;row<height;row++)
  1757. {
  1758.          DisposeImg[row*width+(Cols[top]+range)]=255;
  1759. }
  1760. }
  1761. }
  1762. else
  1763. {
  1764. continue;
  1765. }
  1767. }
  1768. }
  1770.  for(int row=0;row<height;row++)
  1771.  {
  1772. for(int col=0;col<width;col++)
  1773. {
  1774.     if(DisposeImg[row*width+col]<255)
  1775. {
  1776.                DisposeImg[row*width+col]=0;
  1777. }
  1778. }
  1779.  }
  1781. int Area=-1;                   //标识区域的地点
  1782. if(DisposeImg[10*width+0]==255)
  1783. {
  1784. Area++;
  1785. Widths[Area]++;
  1786.     Left[Area]=0;
  1787. }
  1788. for(col=1;col<width;col++)
  1789. {
  1790. if(DisposeImg[10*width+col]==255 && DisposeImg[10*width+(col-1)]==0)
  1791. {
  1792. Area++;
  1793. Widths[Area]++;
  1794.      Left[Area]=col;
  1795. }
  1796. else 
  1797. {
  1798. if((DisposeImg[10*width+col]==255)&&(DisposeImg[10*width+col-1]==255))
  1799. {
  1800.      Widths[Area]++;
  1801. }     
  1802. }
  1803.  }
  1804. }
  1806. // 纵向定位
  1807. void CDipView::VMark(CDipDoc* DisposeImg, int width, int height)
  1808. {
  1809. for(int x=0;x<30;x++)
  1810. {
  1811. int w=Widths[x];
  1812. if(w>100)
  1813. {
  1814.     int c=0;
  1815. int c0=0;
  1816. c=Left[x];
  1817. c0=Left[x]+w-10;
  1819. if(c!=0&&c0>0)
  1820. {
  1821. //for(int row=0;row<height;row++)
  1822. //{
  1823. // DisposeImg->ImgData[row*width+c]=255;
  1824. // DisposeImg->ImgData[row*width+c0]=255;
  1825. //}
  1826. Widths[x]=Widths[x]-10;
  1827. }
  1828. else
  1829. {
  1830. if(c<=0)
  1831. {
  1832. c=5;
  1833. //for(int row=0;row<height;row++)
  1834. //{
  1835. // DisposeImg->ImgData[row*width+c]=255;
  1836. // DisposeImg->ImgData[row*width+c0]=255;
  1837. //}
  1838. Left[x]=c;
  1839. }
  1840. else
  1841. {
  1842. if(c0<=0)
  1843. {
  1844. c0=c0+10;
  1845. //for(int row=0;row<height;row++)
  1846.  //   {
  1847. // DisposeImg->ImgData[row*width+c]=255;
  1848. // DisposeImg->ImgData[row*width+c0]=255;
  1849. //}
  1850. //Widths[x]=Widths[x]-5;
  1851. }
  1852. }
  1853. }
  1854. }
  1855. }
  1857. for(x=0;x<30;x++)
  1858. {
  1859. int w=Widths[x];
  1860. if(w<100)
  1861. {
  1862. Widths[x]=0;
  1863. Left[x]=-1;
  1864. }
  1865. }
  1866. }
  1868. /// <summary>
  1869. /// CDipView::DrawRect(CDipDoc* DisposeImg, int width, int height)的摘要说明。
  1870. /// 创建人  ：刘陈炜
  1871. /// 创建时间：2004.5.1
  1872. /// 函数功能：利用先验知识标识车牌区域，进行车牌区域的选择，（横纵向）
  1873. /// 函数说明：国家标准：车牌的宽高比近似为3：1；
  1874. /// 程序流程：1.初始化存储数组
  1875. ///           2.选择区域高宽比大于3的待选域
  1876. /// </summary> 
  1877. void CDipView::DrawRect(CDipDoc* DisposeImg, int width, int height)
  1878. {
  1879. for(int x=0;x<50;x++)
  1880. {
  1881. Rect_X[x]=-1;
  1882. Rect_Y[x]=-1;
  1883. Rect_W[x]=0;
  1884. Rect_H[x]=0;
  1885. }
  1887. Rect_Areas=0;
  1889. for(int m=0;m<20;m++)
  1890. {
  1891. for(int n=0;n<30;n++)
  1892. {
  1893. if((Widths[n]>3*Height[m])&&(Widths[n]*Height[m]!=0))
  1894. {
  1895. Rect_Areas++;
  1896. Rect_X[Rect_Areas-1]=Left[n];
  1897. Rect_Y[Rect_Areas-1]=TopRow[m];
  1898. Rect_W[Rect_Areas-1]=Widths[n];
  1899. Rect_H[Rect_Areas-1]=Height[m];
  1900. }
  1901. }
  1902. }
  1904. //for( x=0;x<Rect_Areas;x++)
  1905. //{
  1906. // for(int row=Rect_Y[x];row<(Rect_Y[x]+Rect_H[x]);row++)
  1907. // {
  1908. // DisposeImg->ImgData[row*width+Rect_X[x]]=255;
  1909. // DisposeImg->ImgData[row*width+Rect_X[x]+Rect_W[x]]=255;
  1910. // }
  1911. //}
  1913. //for( x=0;x<Rect_Areas;x++)
  1914. //{
  1915. // for(int col=Rect_X[x];col<(Rect_X[x]+Rect_W[x]);col++)
  1916. // {
  1917. // DisposeImg->ImgData[Rect_Y[x]*width+col]=255;
  1918. // DisposeImg->ImgData[(Rect_Y[x]+Rect_H[x])*width+col]=255;
  1919. // }
  1920. //}
  1921. }
  1924. void CDipView::OnChoose() 
  1925. {
  1926. // TODO: Add your command handler code here
  1927. }
  1929. /// <summary>
  1930. /// CDipView::Choose(CDipDoc *DisposeImg, int width, int height)的摘要说明。
  1931. /// 创建人  ：刘陈炜
  1932. /// 创建时间：2004.5.1
  1933. /// 函数功能：再次利用利用跳变数选择车牌区域，进行车牌区域的最终选择，（横纵向）
  1934. /// 程序流程：1.初始化存储数组
  1935. ///           2.计算每个伪车牌区域的跳变数
  1936. ///           3.选择具有最大跳变数的区域为车牌区
  1937. /// </summary> 
  1938. void CDipView::Choose(CDipDoc *DisposeImg, int width, int height)
  1939. {
  1940. int row=0;
  1941. int col=0;
  1943. int temp=0;
  1945. int *JumpNums=new int[Rect_Areas];
  1946. for(int area=0;area<Rect_Areas;area++)
  1947. {
  1948. JumpNums[area]=0;
  1949. }
  1951. for(area=0;area<Rect_Areas;area++)    //计算每个伪车牌区域的行跳变数
  1952. {
  1953. BYTE* TempImg=new BYTE[Rect_W[area]*Rect_H[area]];
  1955. for(row=0;row<Rect_H[area];row++)
  1956. {
  1957. for(col=0;col<Rect_W[area];col++)
  1958. {
  1959. TempImg[row*Rect_W[area]+col]=DisposeImg->ImgData[(Rect_Y[area]+row)*width+(Rect_X[area]+col)];
  1960. }
  1961. }
  1963. JumpNums[area]=JumpNum(TempImg,Rect_W[area],Rect_H[area]);
  1965. free(TempImg);
  1966. }
  1968. temp=MaxJump(JumpNums,Rect_Areas);
  1970. // 区域刻画
  1971. /* for(row=Rect_Y[temp];row<(Rect_Y[temp]+Rect_H[temp]);row++)
  1972. {
  1973. DisposeImg->ImgData[row*width+Rect_X[temp]]=255;
  1974. DisposeImg->ImgData[row*width+(Rect_X[temp]+Rect_W[temp])]=255;
  1975. }
  1977.     for(col=Rect_X[temp];col<(Rect_X[temp]+Rect_W[temp]);col++)
  1978. {
  1979. DisposeImg->ImgData[Rect_Y[temp]*width+col]=255;
  1980.      DisposeImg->ImgData[(Rect_Y[temp]+Rect_H[temp])*width+col]=255;
  1981. } */
  1983. //剔除伪车牌区
  1984. for(area=0;area<Rect_Areas;area++)    
  1985. {
  1986. if(area!=temp)
  1987. {
  1988. Rect_X[area]=-1;
  1989. Rect_Y[area]=-1;
  1990. Rect_W[area]=0;
  1991. Rect_H[area]=0;
  1992. }
  1993. }
  1995. //移至第一位
  1996. Rect_X[0]=Rect_X[temp];
  1997. Rect_Y[0]=Rect_Y[temp];
  1998. Rect_W[0]=Rect_W[temp];
  1999. Rect_H[0]=Rect_H[temp];
  2001. //清空
  2002. if(temp!=0)
  2003. {
  2004. Rect_X[temp]=-1;
  2005. Rect_Y[temp]=-1;
  2006. Rect_W[temp]=0;
  2007. Rect_H[temp]=0;
  2008. }
  2010. Rect_Areas=1;
  2011. }
  2013. /// <summary>
  2014. /// CDipView::JumpNum(BYTE *DisposeImg,int width,int height)的摘要说明。
  2015. /// 创建人  ：刘陈炜
  2016. /// 创建时间：2004.5.1
  2017. /// 函数功能：计算伪车牌区的跳变平均数
  2018. /// 程序流程：1.初始化；
  2019. ///           2.水平差分，边缘提取；
  2020. ///           3.选择车牌的中间1/3行来计算每个伪车牌区域的跳变平均数；
  2021. /// </summary>
  2022. /// <return>各个区域的平均跳变数</return>
  2023. int CDipView::JumpNum(BYTE *DisposeImg,int width,int height)
  2024. {
  2025. int startRow=0;
  2026. int endRow=0;
  2027. long JumpNum=0;
  2028. int averJump=0;
  2029. int r=0;
  2031. BYTE* g=new BYTE[width*height];
  2032.     CopyImg(DisposeImg,g,width,height);
  2034. HSub(g,width,height);
  2036. startRow=(int)(0+height/3);
  2037. endRow=(int)(0+2*height/3);
  2039. for(int row=startRow;row<endRow;row++)
  2040. {
  2041. for(int col=0;col<width-1;col++)
  2042. {
  2043. if(g[row*width+col]!=g[row*width+col+1])
  2044. JumpNum++;
  2045. }
  2046. }
  2048. r=endRow-startRow;
  2049. averJump=(int)(JumpNum/r);
  2051. free(g);
  2053. return averJump;
  2054. }
  2056. /// <summary>
  2057. /// CDipView::MaxJump(int *JumpNums,int num)的摘要说明。
  2058. /// 创建人  ：刘陈炜
  2059. /// 创建时间：2004.5.1
  2060. /// 函数功能：找出所有伪车牌区域中具有最大跳变平均数的区域号
  2061. /// 程序说明：选择排序
  2062. /// </summary>
  2063. /// <return>最大跳变平均数的区域号</return>
  2064. int CDipView::MaxJump(int *JumpNums,int num)
  2065. {
  2066. int area=0;
  2067. int temp=0;
  2068. int val=0;
  2070. val=JumpNums[0];
  2071. for(area=0;area<num;area++)
  2072. {
  2073. if(JumpNums[area]>val)
  2074. {
  2075. temp=area;
  2076. val=JumpNums[area];
  2077. }
  2078. }
  2080. return temp;
  2081. }
  2083. /// <summary>
  2084. /// CDipView::FurtherLocate(CDipDoc *DisposeImg,int width,int height)的摘要说明。
  2085. /// 创建人  ：刘陈炜
  2086. /// 创建时间：2004.5.1
  2087. /// 函数功能：精确定位车牌
  2088. /// 程序流程：
  2089. /// 程序说明：由于水平分割出的伪车牌区域已经比较准确，故仅作垂直方向的进一步定位（缩小宽度）
  2090. /// </summary>
  2091. void CDipView::FurtherLocate(CDipDoc *DisposeImg,int width,int height)
  2092. {
  2093. int row=0;
  2094. int col=0;
  2096. //copy 到一个暂存数组，以便处理操作且不破坏原图象
  2097. BYTE* TempImg=new BYTE[Rect_W[0]*Rect_H[0]];
  2098. for(row=0;row<Rect_H[0];row++)
  2099. {
  2100. for(col=0;col<Rect_W[0];col++)
  2101. {
  2102. TempImg[row*Rect_W[0]+col]=DisposeImg->ImgData[(Rect_Y[0]+row)*width+(Rect_X[0]+col)];
  2103. }
  2104. }
  2106. VFurthLocate(TempImg,Rect_W[0],Rect_H[0]);
  2108. /* int temp=0;
  2109. // 区域刻画
  2110. for(row=Rect_Y[temp];row<(Rect_Y[temp]+Rect_H[temp]);row++)
  2111. {
  2112. DisposeImg->ImgData[row*width+Rect_X[temp]]=255;
  2113. DisposeImg->ImgData[row*width+(Rect_X[temp]+Rect_W[temp])]=255;
  2114. }
  2116.     for(col=Rect_X[temp];col<(Rect_X[temp]+Rect_W[temp]);col++)
  2117. {
  2118. DisposeImg->ImgData[Rect_Y[temp]*width+col]=255;
  2119.      DisposeImg->ImgData[(Rect_Y[temp]+Rect_H[temp])*width+col]=255;
  2120. } */
  2122. free(TempImg);
  2123. }
  2125. /// <summary>
  2126. /// CDipView::VFurthLocate(BYTE *DisposeImg, int width, int height)的摘要说明。
  2127. /// 创建人  ：刘陈炜
  2128. /// 创建时间：2004.5.1
  2129. /// 函数功能：找出车牌的左右边缘
  2130. /// 程序流程：1.水平差分，二值化；
  2131. ///           2.垂直方向投影，统计各列的投影值；
  2132. ///           3.去除噪声点；
  2133. ///           4.找出左右边缘
  2134. ///           5.计算车牌宽度
  2135. ///           6.车牌精确定位赋值
  2136. /// </summary>
  2137. void CDipView::VFurthLocate(BYTE *DisposeImg, int width, int height)
  2138. {
  2139. int row=0;
  2140. int col=0;
  2141. int num=0;
  2143. int leftCol=0;
  2144. int rightCol=width-1;
  2146. HSub(DisposeImg,width,height); 
  2147. ForHandle(DisposeImg,width,height);                           //二值化
  2149.     //note the value of hProjection
  2150. int*HProj=new int[width];                        
  2151. for(col=0;col<width;col++)
  2152. {
  2153. HProj[col]=0;
  2154. }    
  2155. for(col=0;col<width;col++)
  2156. {
  2157. for(row=0;row<height;row++)
  2158. {
  2159. if(DisposeImg[row*width+col]==255)
  2160. HProj[col]++;
  2161. }
  2162. }
  2163. HProj[0]=0;
  2164. for(col=0;col<width;col++)                    //去除噪声点
  2165. {
  2166. if(HProj[col]<3)
  2167. {
  2168. HProj[col]=0;
  2169. }
  2170. }    
  2172. ///////////////////////////////////////////找出左边缘
  2173. int nextCol=0;
  2174. int w=0;
  2175. for(col=0;col<width;col++)
  2176. {
  2177. num=0;
  2178. if(HProj[col]==0&&HProj[col+1]!=0)
  2179. {
  2180. num++;
  2181. nextCol=FindNextCol_LR(HProj,width,col);
  2182. w=nextCol-col;
  2183. leftCol=col;
  2184. }
  2185. if(w<3)
  2186. {
  2187. num=num-1;
  2188. }
  2189. if(num==1)
  2190. {
  2191. break;
  2192. }
  2193. }
  2195. //////////////////////////////////////////找出右边缘
  2196. nextCol=0;
  2197. w=0;
  2198. for(col=width;col>=0;col--)
  2199. {
  2200. num=0;
  2201. if(HProj[col]==0&&HProj[col-1]!=0)
  2202. {
  2203. num++;
  2204. nextCol=FindNextCol_RL(HProj,width,col);
  2205. w=col-nextCol;
  2206. rightCol=col;
  2207. }
  2208. if(w<3)
  2209. {
  2210. num=num-1;
  2211. }
  2212. if(num==1)
  2213. {
  2214. break;
  2215. }
  2216. }
  2218. int blantCols=0;
  2219. blantCols=leftCol+(width-rightCol);
  2221. Rect_X[0]+=leftCol;
  2222. Rect_W[0]=width-blantCols;
  2224. free(HProj);
  2225. }
  2227. /// <summary>
  2228. /// CDipView::ForHandle(BYTE *DisposeImg, int width, int height)的摘要说明。
  2229. /// 创建人  ：刘陈炜
  2230. /// 创建时间：2004.5.1
  2231. /// 函数功能：二值化图象
  2232. /// 程序说明：threshold=average+delt；
  2233. /// 程序流程：1.计算图象均值；
  2234. ///           2.计算图象方差；
  2235. ///           3.二值化
  2236. /// </summary>
  2237. void CDipView::ForHandle(BYTE *DisposeImg, int width, int height)
  2238. {
  2239.  BYTE *temp=new BYTE[height*width];
  2240.  int x,y;
  2242.  CopyImg(DisposeImg,temp,width,height);
  2244.  long double total=0;
  2245.  float aver=0;
  2246.  for(x=0;x<height;x++)
  2247.  {
  2248.  for(y=0;y<width;y++)
  2249.  {
  2250.  total+=DisposeImg[x*width+y];
  2251.  }
  2252.  }
  2253.  aver=(BYTE)(total/(float)(height*width));
  2255.  float delt; 
  2256.  float sub=0; 
  2257.  for(x=0;x<height;x++)
  2258.  {
  2259.  for(y=0;y<width;y++)
  2260.  {
  2261.  sub+=(DisposeImg[x*width+y]-aver)*(DisposeImg[x*width+y]-aver);  
  2262.  }
  2263.  }
  2264.  delt=(float)(sqrt(sub/(float)(height*width)));
  2266.  BYTE judge;
  2267.  judge=(BYTE)(delt+aver);
  2268.  for(x=0;x<height;x++)
  2269.  {
  2270. for(y=0;y<width;y++)
  2271.  {
  2272. if(DisposeImg[(x*width)+y]>=judge)
  2273. {
  2274. DisposeImg[(x*width)+y]=255;
  2275. }
  2276. else
  2277. {
  2278. DisposeImg[(x*width)+y]=0;
  2279. }
  2280. }
  2281. }
  2282.  
  2283.  free(temp);
  2284. }
  2286. // 寻找left==255,right==0的列
  2287. int CDipView::FindNextCol_LR(int* HProj, int size, int objectCol)
  2288. {
  2289. int num=0;
  2290. int x=0;
  2291. for(x=objectCol;x<size;x++)
  2292. {
  2293. if(HProj[x]>=0&&HProj[x+1]==0)
  2294. {
  2295. num++;
  2296. }
  2298. if(num==1)
  2299. {
  2300. break;
  2301. }
  2302. }
  2304. if(num==0)
  2305. {
  2306. return objectCol; 
  2307. }
  2308. else
  2309. {
  2310. return x;
  2311. }
  2312. }
  2314. // 寻找right==255,left==0的列
  2315. int CDipView::FindNextCol_RL(int* HProj, int size, int objectCol)
  2316. {
  2317. int num=0;
  2318. int x=0;
  2319. for(x=objectCol;x>=0;x--)
  2320. {
  2321. if(HProj[x]>=0&&HProj[x-1]==0)
  2322. {
  2323. num++;
  2324. }
  2326. if(num==1)
  2327. {
  2328. break;
  2329. }
  2330. }
  2332. if(num==0)
  2333. {
  2334. return objectCol; 
  2335. }
  2336. else
  2337. {
  2338. return x;
  2339. }
  2340. }
  2342. void CDipView::OnSobel()
  2343. {
  2344. // TODO: 在此添加命令处理程序代码
  2345. }
  2348. /// <summary>
  2349. /// CDipView::OnVSplit()的摘要说明。
  2350. /// 创建人  ：刘陈炜
  2351. /// 创建时间：2004.5.1
  2352. /// 函数功能：车牌字符分割
  2353. /// 程序流程：1.区域扩展；
  2354. ///           2.投影分割；
  2355. ///           3.字符区域的获取；
  2356. /// </summary>
  2357. void CDipView::OnVSplit()
  2358. {
  2359. // TODO: 在此添加命令处理程序代码
  2360. CDipDoc *pDoc=GetDocument();
  2362. int width,height;
  2363. width=pDoc->ImgWidth;
  2364. height=pDoc->ImgHeight;
  2366. int row=0;
  2367. int col=0;
  2369. //extend area
  2370. int threshold=2;
  2372. Rect_X[0]=Rect_X[0];
  2373. Rect_Y[0]=Rect_Y[0];//-threshold;
  2374. Rect_W[0]=Rect_W[0]+7;
  2375. Rect_H[0]=Rect_H[0];//+2*threshold;
  2377. //车牌区域的高度、宽度
  2378. int w=Rect_W[0];
  2379. int h=Rect_H[0];
  2381. BYTE* TempImg=new BYTE[w*h];
  2382. for(row=0;row<h;row++)
  2383. {
  2384. for(col=0;col<w;col++)
  2385. {
  2386. TempImg[row*w+col]=pDoc->ImgData[(Rect_Y[0]+row)*width+(Rect_X[0]+col)];
  2387. }
  2388. }
  2390. //二值化车牌区域图象
  2391. ForHandle(TempImg,w,h);
  2393. //搜索各字符域
  2394. int *projArray=new int[w];
  2395. for(col=0;col<w;col++)
  2396. {
  2397. projArray[col]=0;
  2398. }
  2400. HProj(TempImg,projArray,w,h);
  2401. CharacterSplit(projArray,w);
  2402. // GetCharacter();
  2404. ///////////////////////////////标记出各个字符区域
  2405. int area;
  2406. area=CharacterNum;
  2407. for(;area>=0;area--)
  2408. {
  2409. if(cRects[area].w!=0&&cRects[area].x!=-1)
  2410. {
  2411. for(row=cRects[area].y;row<(cRects[area].y+cRects[area].h);row++)  //draw vertical line
  2412. {
  2413. pDoc->ImgData[row*width+cRects[area].x]=255;
  2414. pDoc->ImgData[row*width+(cRects[area].x+cRects[area].w)]=255;
  2415. }
  2417.      for(col=cRects[area].x;col<(cRects[area].x+cRects[area].w);col++)    //draw horizontal line
  2418. {
  2419. pDoc->ImgData[cRects[area].y*width+col]=255;
  2420. pDoc->ImgData[(cRects[area].y+cRects[area].h)*width+col]=255;
  2421. }
  2422. }
  2423. }
  2425. // 区域刻画
  2426. //int temp=0;
  2427. //for(row=Rect_Y[temp];row<(Rect_Y[temp]+Rect_H[temp]);row++)
  2428. //{
  2429. // pDoc->ImgData[row*width+Rect_X[temp]]=255;
  2430. // pDoc->ImgData[row*width+(Rect_X[temp]+Rect_W[temp])]=255;
  2431. //}
  2432. //
  2433.  //   for(col=Rect_X[temp];col<(Rect_X[temp]+Rect_W[temp]);col++)
  2434. //{
  2435. // pDoc->ImgData[Rect_Y[temp]*width+col]=255;
  2436.  //    pDoc->ImgData[(Rect_Y[temp]+Rect_H[temp])*width+col]=255;
  2437. //}
  2439. //for(row=Rect_Y[temp];row<(Rect_Y[temp]+Rect_H[temp]);row++)
  2440. //{
  2441. // for(col=0;col<w;col++)
  2442. // {
  2443. // if(projArray[col]==1)
  2444. // {
  2445. // pDoc->ImgData[row*width+Rect_X[temp]+col]=255;
  2446. // // pDoc->ImgData[row*width+(Rect_X[temp]+Rect_W[temp])]=255;
  2447. // }
  2448. // }
  2449. //}
  2451. free(TempImg);
  2452. free(projArray);
  2454. IsNewBitmap=true;
  2455. IsNewFile=true;
  2456. pDoc->UpdateAllViews(NULL); 
  2457. }
  2459. // 图像规定化0~100
  2460. void CDipView::GuiYiHua(BYTE* DisposeImg, int width , int height)
  2461. {
  2462. }
  2464. // 增强车牌区,削弱背景
  2465. void CDipView::Hance(BYTE* DisposeImg, int width , int height)
  2466. {
  2467. }
  2469. void CDipView::OnGuiYiHua()
  2470. {
  2471. // TODO: 在此添加命令处理程序代码
  2472. }
  2474. void CDipView::OnHance()
  2475. {
  2476. // TODO: 在此添加命令处理程序代码
  2477. }
  2479. /// <summary>
  2480. /// CDipView::HProj(BYTE* DisposeImg, int* projArray, int width , int height)的摘要说明。
  2481. /// 创建人  ：刘陈炜
  2482. /// 创建时间：2004.5.1
  2483. /// 函数功能：水平方向投影，分割出字符
  2484. /// 程序流程：1.统计各列白象素个数；
  2485. ///           2.平滑投影曲线；
  2486. ///           3.寻找波谷（trough）进行分割；
  2487. /// </summary>
  2488. void CDipView::HProj(BYTE* DisposeImg, int* projArray, int width , int height)
  2489. {
  2490. int row=0;
  2491. int col=0;
  2493. float alph=0.0;
  2494. int*E=new int[width];
  2495. long double*MSE=new long double[9];                         //针对每一个alph的deltl;
  2497. //统计各列白象素个数
  2498. for(col=0;col<width;col++)
  2499. {
  2500. for(row=0;row<height;row++)
  2501. {
  2502. if(DisposeImg[row*width+col]==255)
  2503. projArray[col]++;
  2504. }
  2505. }
  2507. //一次平滑
  2508. for(int x=1;x<10;x++)
  2509. {
  2510. alph=(float)(0.1*x);
  2511. E[0]=projArray[0];
  2512. for(col=1;col<width;col++)
  2513. {
  2514. E[col]=(int)(alph*projArray[col]+(1-alph)*E[col-1]);
  2515. }
  2517. for(col=0;col<width;col++)
  2518. {
  2519. MSE[x-1]+=(projArray[col]-E[col])*(projArray[col]-E[col]);
  2520. }
  2521. MSE[x-1]=MSE[x-1]/width;
  2522. }
  2524. long double min=0.0;
  2525. for(x=1;x<10;x++)
  2526. {
  2527. if(MSE[x-1]>min)
  2528. {
  2529. min=MSE[x-1];
  2530. }
  2531. }
  2532. for(x=1;x<10;x++)
  2533. {
  2534. if(MSE[x-1]==min)
  2535. {
  2536. alph=(float)(0.1*x);
  2537. }
  2538. }
  2540. E[0]=projArray[0];
  2541. for(col=1;col<width;col++)
  2542. {
  2543. E[col]=(int)(alph*projArray[col]+(1-alph)*E[col-1]);
  2544. }
  2546. for(col=1;col<width;col++)
  2547. {
  2548. projArray[col]=E[col];
  2549. E[col]=0;
  2550. }
  2552. //二次平滑
  2553. for(x=1;x<10;x++)
  2554. {
  2555. alph=(float)(0.1*x);
  2556. E[0]=projArray[0];
  2557. for(col=1;col<width;col++)
  2558. {
  2559. E[col]=(int)(alph*projArray[col]+(1-alph)*E[col-1]);
  2560. }
  2562. for(col=0;col<width;col++)
  2563. {
  2564. MSE[x-1]+=(projArray[col]-E[col])*(projArray[col]-E[col]);
  2565. }
  2566. MSE[x-1]=MSE[x-1]/width;
  2567. }
  2569. min=0.0;
  2570. for(x=1;x<10;x++)
  2571. {
  2572. if(MSE[x-1]<min)
  2573. {
  2574. min=MSE[x-1];
  2575. }
  2576. }
  2577. for(x=1;x<10;x++)
  2578. {
  2579. if(MSE[x-1]==min)
  2580. {
  2581. alph=(float)(0.1*x);
  2582. }
  2583. }
  2585. E[0]=projArray[0];
  2586. for(col=1;col<width;col++)
  2587. {
  2588. E[col]=(int)(alph*projArray[col]+(1-alph)*E[col-1]);
  2589. }
  2590. for(col=1;col<width;col++)
  2591. {
  2592. projArray[col]=E[col];
  2593. E[col]=0;
  2594. }
  2596.     //寻找波谷进行分割trough
  2597. FindTrough(projArray,width); //获得波谷位置
  2599. free(E);
  2600. free(MSE);
  2601. }
  2603. // 寻找波谷,划分字符区间
  2604. void CDipView::FindTrough(int* Array, int size)
  2605. {
  2606. int col=0;
  2607. bool*flage=new bool[size];
  2609. bool LFlage=false;                                //a wave start leftside
  2610. bool RFlage=false;                                //a wave end   rightside
  2611. int  LSide=0;                                     //left side
  2612. int  RSide=0;   //right side
  2613. int  length=0;                                    //area width
  2615. for(col=0;col<size;col++)
  2616. {
  2617. flage[col]=false;
  2618. }
  2620. for(col=1;col<size-1;col++)
  2621. {
  2622. if(Array[col]==0&&Array[col+1]!=0&&LFlage==false)
  2623. {
  2624. flage[col]=true;
  2625. LFlage=true;
  2626. }
  2627. if((Array[col-1]!=0)&&(Array[col]==0)&&(LFlage==true)&&(RFlage==false))
  2628. {
  2629. flage[col]=true;
  2630. RFlage=true;
  2631. }
  2632. if((LFlage==true)&&(RFlage==true))
  2633. {
  2634. LFlage=false;
  2635. RFlage=false;
  2636. }
  2637. }
  2639. for(col=0;col<size;col++)
  2640. {
  2641. Array[col]=0;
  2642. if(flage[col]==true)
  2643. {
  2644. Array[col]=1;
  2645. }
  2646. }
  2647. }
  2649. // 分割字符
  2650. void CDipView::CharacterSplit(int* Array, int size)
  2651. {
  2652. CDipDoc *pDoc=GetDocument();
  2654. int width,height;
  2655. width=pDoc->ImgWidth;
  2656. height=pDoc->ImgHeight;
  2658.     int x=Rect_X[0];
  2659. int y=Rect_Y[0];
  2660. int w=Rect_W[0];
  2661. int h=Rect_H[0];
  2663.     int num=0;
  2665. for(num=0;num<10;num++)
  2666. {
  2667. cRects[num].x=-1;
  2668. cRects[num].y=-1;
  2669. cRects[num].w=0;
  2670. cRects[num].h=0;
  2671. }
  2673. bool flage=false;
  2676. int row=0;
  2677. int col=0;
  2679. int area=0;
  2680. num=0;
  2681. for(col=0;col<size;col++)
  2682. {
  2683. if(Array[col]==1)
  2684. {
  2685. num++;
  2686. if(num%2!=0&&flage==false)
  2687. {
  2689. flage=true;
  2690. cRects[area].x=x+col;
  2691. }
  2692. else if(num%2==0&&flage==true)
  2693. {
  2694. flage=false;
  2695. cRects[area].w=col-(cRects[area].x-x);
  2696. cRects[area].y=y;
  2697. cRects[area].h=h;
  2698. area++;
  2699. }
  2700. }
  2701. }
  2703. CharacterNum=area;
  2705. // CharacterNum=area+1;
  2707. for(;area>=0;area--)
  2708. {
  2709. if(cRects[area].w!=0)
  2710. {
  2711. for(row=cRects[area].y;row<(cRects[area].y+cRects[area].h);row++)  //draw vertical line
  2712. {
  2713. pDoc->ImgData[row*width+cRects[area].x-1]=255;
  2714. pDoc->ImgData[row*width+(cRects[area].x+cRects[area].w)]=255;
  2715. }
  2717. for(col=cRects[area].x;col<(cRects[area].x+cRects[area].w);col++)    //draw horizontal line
  2718. {
  2719. pDoc->ImgData[(cRects[area].y-1)*width+col]=255;
  2720. pDoc->ImgData[(cRects[area].y+cRects[area].h)*width+col]=255;
  2721. }
  2722. }
  2723. }
  2725. //save character
  2726.     for(area=CharacterNum;area>=0;area--)
  2727. {
  2728. cDatas[area].w=cRects[area].w;
  2729. cDatas[area].h=cRects[area].h;
  2730. cDatas[area].Img=new BYTE[cRects[area].w*cRects[area].h];
  2732. for(row=cRects[area].y;row<(cRects[area].y+cRects[area].h);row++)  //draw vertical line
  2733. {
  2734. for(col=cRects[area].x;col<(cRects[area].x+cRects[area].w);col++)    //draw horizontal line
  2735. {
  2736. cDatas[area].Img[(row-cRects[area].y)*cRects[area].w+(col-cRects[area].x)]=pDoc->ImgData[row*width+col];
  2737. }
  2738. }
  2739. }
  2740. }
  2741. /*
  2742. // 读取样本库，传入一个sample类型的数组（从文件里读出的样本数据将放入其中），返回值为样本个数
  2743. int CDipView::LoadCharLib(sample* sa)
  2744. {
  2745. CFile cf;
  2747. if(cf.Open(".\mydata.dat",CFile::modeCreate|CFile::modeNoTruncate|CFile::modeReadWrite)==0)
  2748. {
  2749. AfxMessageBox("打开文件失败，n您最好退出程序");
  2750. }
  2752. int num=0;
  2753. try
  2754. {
  2755. sample sampleTemp;
  2757. cf.SeekToBegin();
  2758. DWORD dwBytesRemaining=(DWORD)(cf.GetLength());
  2759. while(dwBytesRemaining)
  2760. {
  2761. UINT nBytesRead=cf.Read(&sampleTemp,sizeof(sample));
  2762. sa[num]=sampleTemp;
  2763. dwBytesRemaining-=nBytesRead;
  2764. num++;
  2765. }
  2766. CString s;
  2767. // s.Format("共读取数据%d个,最后一个的序列号为%d",
  2768. // num,sampleTemp.serialnum);
  2769. // AfxMessageBox(s);
  2770. }
  2771. catch(CFileException *e)
  2772. {
  2773. e->ReportError();
  2774. e->Delete();
  2775. }
  2776. cf.Close();
  2777. return num;
  2778. }
  2780. // 缩放算法,一律缩放为20*36
  2781. void CDipView::Zoom(BYTE* DisposeImg, int width, int height)
  2782. {
  2783. if ((width==20)&&(height==36))
  2784. return;
  2786. // 源图像的宽度和高度
  2787. LONG lWidth=width;
  2788. LONG lHeight=height;
  2790. BYTE *temp_img=new BYTE[width*height];
  2792. // 循环变量
  2793. int i;
  2794. int j;
  2796. for(i = 0; i < height; i++)
  2797. {
  2798. for(j = 0; j < width; j++)
  2799. {
  2800. temp_img[j+i*width]=DisposeImg[j+i*width];
  2801. }
  2802. }
  2804. // 缩放后图像的宽度和高度
  2805. LONG lNewWidth;
  2806. LONG lNewHeight;
  2808. float XZRatio=(float)( 20.0/(float)(width));
  2809. float YZRatio=(float)(36.0/(float)(height));
  2811. // 象素在源坐标
  2812. LONG i0;
  2813. LONG j0;
  2815. //归一后的大小
  2816. lNewWidth = 20;
  2817. lNewHeight = 36;
  2819. height=lNewHeight;
  2820. width=lNewWidth;
  2823. for(i = 0;i< lNewHeight;i++)
  2824. {
  2826. for(j = 0;j< lNewWidth;j++)
  2827. {
  2828. i0 = (LONG) (i / YZRatio + 0.5);
  2829. j0 = (LONG) (j / XZRatio + 0.5);
  2831. // 判断是否在源图范围内
  2832. if( (j0 >= 0) && (j0 < lWidth) && (i0 >= 0) && (i0 < lHeight))
  2833. {
  2834. // 复制象素
  2835. DisposeImg[j+i*lWidth] = temp_img[j0 + i0*lWidth];
  2836. }
  2837. else
  2838. {
  2839. // 对于源图中没有的象素，直接赋值为255
  2840. // DisposeImg[j+i*lWidth] = 0;
  2841. }
  2842. }
  2843. }
  2844. delete [] temp_img;
  2845. }
  2847. // 二值化算法,二值化为0和1两种值
  2848. void CDipView::BinaryImg(BYTE* DisposeImg, int width , int height)
  2849. {
  2850.  BYTE *temp=new BYTE[height*width];
  2851.  int x,y;
  2853.  CopyImg(DisposeImg,temp,width,height);
  2855.  long double total=0;
  2856.  float aver=0;
  2857.  for(x=0;x<height;x++)
  2858.  {
  2859.  for(y=0;y<width;y++)
  2860.  {
  2861.  total+=DisposeImg[x*width+y];
  2862.  }
  2863.  }
  2864.  aver=(BYTE)(total/(float)(height*width));
  2866.  float delt; 
  2867.  float sub=0; 
  2868.  for(x=0;x<height;x++)
  2869.  {
  2870.  for(y=0;y<width;y++)
  2871.  {
  2872.  sub+=(DisposeImg[x*width+y]-aver)*(DisposeImg[x*width+y]-aver);  
  2873.  }
  2874.  }
  2875.  delt=(float)(sqrt(sub/(float)(height*width)));
  2877.  BYTE judge;
  2878.  judge=(BYTE)(delt+aver);
  2879.  for(x=0;x<height;x++)
  2880.  {
  2881. for(y=0;y<width;y++)
  2882.  {
  2883. if(DisposeImg[(x*width)+y]>=judge)
  2884. {
  2885. DisposeImg[(x*width)+y]=(BYTE)1;
  2886. }
  2887. else
  2888. {
  2889. DisposeImg[(x*width)+y]=(BYTE)0;
  2890. }
  2891. }
  2892.  }
  2893.  
  2894.  free(temp);
  2895. }
  2897. // Hilditch细化算法
  2898. void CDipView::ThinnerHilditch(void* image, unsigned long lx, unsigned long ly)
  2899. {
  2900. char *f, *g;
  2901. char n[10];
  2902. unsigned int counter;
  2903. short k, shori, xx, nrn;
  2904. unsigned long i, j;
  2905. long kk, kk11, kk12, kk13, kk21, kk22, kk23, kk31, kk32, kk33, size;
  2906. size = (long)lx * (long)ly;
  2907. g = (char *)malloc(size);
  2909. if(g == NULL)
  2910. {
  2911. printf("error in allocating memory!n");
  2912. return;
  2913. }
  2915. f = (char *)image;
  2916. for(i=0; i<lx; i++)
  2917. {
  2918. for(j=0; j<ly; j++)
  2919. {
  2920. kk=i*ly+j;
  2921. if(f[kk]!=0)
  2922. {
  2923. f[kk]=1;
  2924. g[kk]=f[kk];
  2925. }
  2926. }
  2927. }
  2929. counter = 1;
  2931. do
  2932. {
  2933. printf("%4d*",counter);
  2934. counter++;
  2935. shori = 0;
  2937. for(i=0; i<lx; i++)
  2938. {
  2939. for(j=0; j<ly; j++)
  2940. {
  2941. kk = i*ly+j;
  2942. if(f[kk]<0)
  2943. f[kk] = 0;
  2944. g[kk]= f[kk];
  2945. }
  2946. }
  2948. for(i=1; i<lx-1; i++)
  2949. {
  2950. for(j=1; j<ly-1; j++)
  2951. {
  2952. kk=i*ly+j;
  2954. if(f[kk]!=1)
  2955. continue;
  2957. kk11 = (i-1)*ly+j-1;
  2958. kk12 = kk11 + 1;
  2959. kk13 = kk12 + 1;
  2960. kk21 = i*ly+j-1;
  2961. kk22 = kk21 + 1;
  2962. kk23 = kk22 + 1;
  2963. kk31 = (i+1)*ly+j-1;
  2964. kk32 = kk31 + 1;
  2965. kk33 = kk32 + 1;
  2967. if((g[kk12]&&g[kk21]&&g[kk23]&&g[kk32])!=0)
  2968. continue;
  2970. nrn = g[kk11] + g[kk12] + g[kk13] + g[kk21] + g[kk23] + 
  2971. g[kk31] + g[kk32] + g[kk33];
  2973. if(nrn <= 1)
  2974. {
  2975. f[kk22] = 2;
  2976. continue;
  2977. }
  2979. n[4] = f[kk11];
  2980. n[3] = f[kk12];
  2981. n[2] = f[kk13];
  2982. n[5] = f[kk21];
  2983. n[1] = f[kk23];
  2984. n[6] = f[kk31];
  2985. n[7] = f[kk32];
  2986. n[8] = f[kk33];
  2987. n[9] = n[1];
  2988. xx = 0;
  2990. for(k=1; k<8; k=k+2)
  2991. {
  2992. if((!n[k])&&(n[k+1]||n[k+2]))
  2993. xx++;
  2994. }
  2996. if(xx!=1)
  2997. {
  2998. f[kk22] = 2;
  2999. continue;
  3000. }
  3002. if(f[kk12] == -1)
  3003. {
  3004. f[kk12] = 0;
  3005. n[3] = 0;
  3006. xx = 0;
  3008. for(k=1; k<8; k=k+2)
  3009. {
  3010. if((!n[k])&&(n[k+1]||n[k+2]))
  3011. xx++;
  3012. }
  3014. if(xx != 1)
  3015. {
  3016. f[kk12] = -1;
  3017. continue;
  3018. }
  3020. f[kk12] = -1;
  3021. n[3] = -1;
  3022. }
  3024. if(f[kk21]!=-1)
  3025. {
  3026. f[kk22] = -1;
  3027. shori = 1;
  3028. continue;
  3029. }
  3031. f[kk21] = 0;
  3032. n[5] = 0;
  3033. xx = 0;
  3035. for(k=1; k<8; k=k+2)
  3036. {
  3037. if((!n[k])&&(n[k+1]||n[k+2]))
  3038. {
  3039. xx++;
  3040. }
  3041. }
  3043. if(xx == 1)
  3044. {
  3045. f[kk21] = -1;
  3046. f[kk22] = -1;
  3047. shori =1;
  3048. }
  3049. else
  3050. f[kk21] = -1;
  3051. }
  3052. }
  3053. }while(shori);
  3055. free(g);
  3056. }
  3058. // 细化算法
  3059. void CDipView::ThinImage(BYTE* image, int width, int height)
  3060. {
  3061. LONG x,y,k;
  3062. BYTE image1[10000];
  3063. k=0;
  3065. LONG digitWidth = width;
  3066. LONG digitHeight = height;
  3068. for(x=0; x<digitWidth; x++)
  3069. {
  3070. image[x*width+0] = (BYTE)0;
  3071. image[x*width+digitHeight-1] = (BYTE)0;
  3072. }
  3074. for(y=0; y<digitHeight; y++)
  3075. {
  3076. image[0*width+y] = (BYTE)0;
  3077. image[(digitWidth-1)*width+y] = (BYTE)0;
  3078. }
  3080. for(x=0; x<digitWidth; x++)
  3081. {
  3082. for(y=0; y<digitHeight; y++)
  3083. {
  3084. image1[k] = image[x*width+y];
  3085. if(image1[k] != 0)
  3086. image1[k] = (BYTE)1;
  3087. k++;
  3088. }
  3089. }
  3091. ThinnerHilditch((void *)image1, digitWidth, digitHeight);
  3093. k=0;
  3094. for(x=0; x<digitWidth; x++)
  3095. {
  3096. for(y=0; y<digitHeight; y++)
  3097. {
  3098. image[x*width+y] = image1[k];
  3099. if(image[x*width+y]!=0)
  3100. image[x*width+y] = (BYTE)1;
  3101. k++;
  3102. }
  3103. }
  3104. }
  3106. // 获取特征值的函数
  3107. void CDipView::GetFeature(BYTE* image, long width, long height, double* feature, int size)
  3108. {
  3109. int i,j;
  3110. for(i=0; i<13; i++)
  3111. feature[i] = 0;
  3113. //图象是20×36大小的，分成9块
  3115. //第一块
  3116. for(i=0; i<7; i++)
  3117. {
  3118. for(j=0; j<12; j++)
  3119. {
  3120. if(image[i*width+j]==1)
  3121. feature[0]+=1.0;
  3122. }
  3123. }
  3125. //第二块
  3126. for(i=0; i<7; i++)
  3127. {
  3128. for(j=12; j<24; j++)
  3129. {
  3130. if(image[i*width+j]==1)
  3131. feature[1]+=1.0;
  3132. }
  3133. }
  3135. //第三块
  3136. for(i=0; i<7; i++)
  3137. {
  3138. for(j=24; j<36; j++)
  3139. {
  3140. if(image[i*width+j]==1)
  3141. feature[2]+=1.0;
  3142. }
  3143. }
  3145. //第四块
  3146. for(i=7; i<13; i++)
  3147. {
  3148. for(j=0; j<12; j++)
  3149. {
  3150. if(image[i*width+j]==1)
  3151. feature[3]+=1.0;
  3152. }
  3153. }
  3155. //第五块
  3156. for(i=7; i<13; i++)
  3157. {
  3158. for(j=12; j<24; j++)
  3159. {
  3160. if(image[i*width+j]==1)
  3161. feature[4]+=1.0;
  3162. }
  3163. }
  3165. //第六块
  3166. for(i=7; i<13; i++)
  3167. {
  3168. for(j=24; j<36; j++)
  3169. {
  3170. if(image[i*width+j]==1)
  3171. feature[5]+=1.0;
  3172. }
  3173. }
  3175. //第七块
  3176. for(i=13; i<20; i++)
  3177. {
  3178. for(j=0; j<12; j++)
  3179. {
  3180. if(image[i*width+j]==1)
  3181. feature[6]+=1.0;
  3182. }
  3183. }
  3185. //第八块
  3186. for(i=13; i<20; i++)
  3187. {
  3188. for(j=12; j<24; j++)
  3189. {
  3190. if(image[i*width+j]==1)
  3191. feature[7]+=1.0;
  3192. }
  3193. }
  3195. //第九块
  3196. for(i=13; i<20; i++)
  3197. {
  3198. for(j=24; j<36; j++)
  3199. {
  3200. if(image[i*width+j]==1)
  3201. feature[8]+=1.0;
  3202. }
  3203. }
  3205. //下面统计方向交点特征
  3207. for(j=0; j<36; j++)
  3208. {
  3209. if(image[7*width+j]==1)
  3210. feature[9]+=1.0;
  3211. }
  3213. for(j=0; j<36; j++)
  3214. {
  3215. if(image[13*width+j]==1)
  3216. feature[10]+=1.0;
  3217. }
  3219. for(i=0; i<20; i++)
  3220. {
  3221. if(image[i*width+12]==1)
  3222. feature[11]+=1.0;
  3223. }
  3225. for(i=0; i<20; i++)
  3226. {
  3227. if(image[i*width+24]==1)
  3228. feature[12]+=1.0;
  3229. }
  3230. }
  3232. // 识别函数，第一个参数为LoadCharLib读出的数组,
  3233. //第二个参数为LoadCharLib的返回值（即样本个数）,
  3234. //将上面获得的特征值传入第三个参数，返回为识别结果
  3235. CString CDipView::Recognize(sample* sa, int num, double* feature)
  3236. {
  3237. int j;
  3238. double min=100000.0;
  3239. CString Result;
  3241. for(int i=0;i<num;i++)
  3242. {
  3243. double diff = 0.0;
  3245. for(j=0; j<9; j++)
  3246. {
  3247. diff += fabs(feature[j] - sa[num].feature[j]);
  3248. }
  3250. for(j=9; j<13; j++)
  3251. {
  3252. diff += fabs(feature[j] - sa[num].feature[j])*9;
  3253. }
  3255. if(diff < min)
  3256. {
  3257. min = diff;
  3258. Result.Format("%s",sa[num].trueClass);
  3259. }
  3260. }
  3261. return Result;
  3262. }*/
  3264. // 字符识别
  3265. /*void CDipView::GetCharacter()
  3266. {
  3267. //CDipDoc *pDoc=GetDocument();
  3269. //int width,height;
  3270. //width=pDoc->ImgWidth;
  3271. //height=pDoc->ImgHeight;
  3273. //sample mysa[500];
  3274.  //   double myfeature[13];
  3275.  //   int num;
  3276.  //   CString result[10];
  3278. //BYTE*Img;
  3279. //int row=0;
  3280. //int col=0;
  3282. //int x=-1;
  3283. //int y=-1;
  3284. //int w=0;
  3285. //int h=0;
  3286. //int area=0;
  3287. //for(area=0;area<10;area++);
  3288. //{
  3289. // if(cRects[area].w>0)
  3290. // {
  3291. // x=cRects[area].x;
  3292. // y=cRects[area].y;
  3293. // w=cRects[area].w;
  3294. // h=cRects[area].h;
  3296.  //           Img=new BYTE[w*h];
  3297. // for(row=0;row<h;row++)
  3298. // {
  3299. // for(col=0;col<w;col++)
  3300. // {
  3301.  //        //           Img[row*w+col]=pDoc->ImgData[y*width+x];
  3302. // }
  3303. // }
  3305. //// //num=LoadCharLib(mysa);                         //获取样本库保存到mysa[500]中
  3306. //// //Zoom(Img,w,h);                                 //对切割出的图象进行大小归一化
  3307. //// //ThinImage(Img,w,h);                            //对大小归一化后的图象进行细化
  3308. //// //GetFeature(Img,w,h,myfeature,13);    //图象细化以后就可以提取特征了
  3309. //// //result[area]=Recognize(mysa,num,myfeature);
  3310. // }
  3311. //}
  3313. //delete []mysa;
  3314. //delete []myfeature;
  3315. //delete []result;
  3316. }*/
  3318. void CDipView::OnCharacterRecg()
  3319. {
  3320. CMainFrame* MainFrame=(CMainFrame*)this->GetParent()->GetParent(); 
  3321. CCharView* pchview=(CCharView*)MainFrame->m_wndSplitter.GetPane(0,1); 
  3323. CDigitClass digitTest;
  3324. CString classResult;
  3325. CString tempstr;
  3326. CString strInfo, strTemp;
  3327. int recgpos=1;
  3328. for(int k=0;k<CharacterNum;k++)
  3329. {
  3330. if ((float) cDatas[k].h/cDatas[k].w>4)
  3331. {
  3332. if (recgpos!=1&&recgpos<8)
  3333. {
  3334. classResult=classResult+"1";
  3335. recgpos++;
  3336. }
  3337. continue;
  3339. }
  3341. for(int x=0;x<cDatas[k].h;x++)
  3342. {
  3343. for(int y=0;y<cDatas[k].w;y++)
  3344. {
  3345. digitTest.digitarray[y][x]=cDatas[k].Img[x*cDatas[k].w+y];
  3346. }
  3347. }
  3348. digitTest.digitHeight=cDatas[k].h;
  3349. digitTest.digitWidth=cDatas[k].w;
  3350. digitTest.FixSize();
  3351. digitTest.BinaryDigit();
  3352. //digitTest.MarrBinary();
  3353. digitTest.ThinDigit_1();
  3355. pchview->m_result.SetSel(0,-1);
  3356. pchview->m_result.ReplaceSel("");
  3358. CString strInfo, strTemp;
  3359. strInfo.Format("");
  3360. int i,j;
  3361. for(i=0; i<digitTest.digitHeight; i++)
  3362. {
  3363. for(j=0; j<digitTest.digitWidth; j++)
  3364. {
  3365. if(digitTest.digitarray[j][i] == 1)
  3366. strInfo += "■";
  3367. else
  3368. strInfo += "  ";
  3369. }
  3370. strInfo += "rn";
  3371. }
  3372. pchview->m_result.SetSel(0,-1);
  3373. pchview->m_result.ReplaceSel(strInfo);
  3374. // AfxMessageBox("hello");
  3375. digitTest.GetFeature();
  3376. if (recgpos==1)
  3377. {
  3379. double min=100000.0;
  3380. POSITION pos=hanList.GetHeadPosition();
  3381. while(pos!=NULL)
  3382. {
  3383. Sample sa;
  3384. sa=hanList.GetNext(pos);
  3385. double diff = 0.0;
  3387. for(j=0; j<9; j++)
  3388. {
  3389. diff += fabs(digitTest.feature[j] - sa.feature[j]);
  3390. }
  3392. for(j=9; j<13; j++)
  3393. {
  3394. diff += fabs(digitTest.feature[j] - sa.feature[j])*9;
  3395. }
  3397. if(diff < min)
  3398. {
  3399. min = diff;
  3400. tempstr.Format("%s",sa.trueClass);
  3401. }
  3402. }
  3404. classResult=classResult+tempstr;
  3406. }
  3408. if (recgpos==2)
  3409. {
  3410. double temp[13];
  3411. double result1[26];
  3412. bool t=0;
  3413. for (i=0;i<13;i++)
  3414. {
  3415. temp[i]=digitTest.feature[i]/15;
  3416. }
  3417. char_bp.Identify(temp,13,result1,26);
  3418. for(i=0;i<26;i++)
  3419. {
  3420. if (ABS(result1[i]-1)<0.1)
  3421. {
  3422. tempstr.Format("%c",i+'A');
  3423. classResult=classResult+tempstr;
  3424. t=1;
  3425. }
  3427. }
  3428. if(t!=1)
  3429. classResult=classResult+'?';
  3430. }
  3431. // strTemp.Format("%d", classResult);
  3432. if (recgpos==3)
  3433. {
  3434. double temp[13];
  3435. double result1[36];
  3436. bool t=0;
  3437. for (i=0;i<13;i++)
  3438. {
  3439. temp[i]=digitTest.feature[i]/15;
  3440. }
  3441. char_num_bp.Identify(temp,13,result1,36);
  3442. for(i=0;i<36;i++)
  3443. {
  3444. if (ABS(result1[i]-1)<0.2)
  3445. {
  3446. if(i<=9)
  3447. {
  3448. tempstr.Format("%c",i+'0');
  3449. classResult=classResult+tempstr;
  3450. t=1;
  3451. }
  3452. else
  3453. {
  3454. tempstr.Format("%c",i+'A'-10);
  3455. classResult=classResult+tempstr;
  3456. t=1;
  3457. }
  3458. // tempstr.Format("%c",i);
  3459. }
  3460. }
  3461. if(t!=1)
  3462. classResult=classResult+'?';
  3463. }
  3465. if (recgpos>3&&recgpos<8)
  3466. {
  3467. double temp[13];
  3468. double result1[10];
  3469. bool t=0;
  3470. for (i=0;i<13;i++)
  3471. {
  3472. temp[i]=digitTest.feature[i]/15;
  3473. }
  3474. char_num_bp.Identify(temp,13,result1,10);
  3475. for(i=0;i<10;i++)
  3476. {
  3477. if (ABS(result1[i]-1)<0.1)
  3478. {
  3479. tempstr.Format("%c",i+'0');
  3480. classResult=classResult+tempstr;
  3481. t=1;
  3482. }
  3483. }
  3484. if(t!=1)
  3485. classResult=classResult+'?';
  3486. }
  3488. // strInfo += "nr识别结果为：  ";
  3489. // strInfo += classResult;
  3490. recgpos++;
  3491. // pchview->m_result.SetSel(0,-1);
  3492. //    pchview->m_result.ReplaceSel(strInfo);
  3495. }
  3496. strInfo += "nr识别结果为：  ";
  3497.     strInfo += classResult;
  3498. pchview->m_result.SetSel(0,-1);
  3499. pchview->m_result.ReplaceSel(strInfo);
  3503. }
  3505. void CDipView::OnTrainfromvehicle() 
  3506. {
  3507. OnLocate();
  3508. OnVSplit();
  3509. CMainFrame* MainFrame=(CMainFrame*)this->GetParent()->GetParent(); 
  3510. CCharView* pchview=(CCharView*)MainFrame->m_wndSplitter.GetPane(0,1); 
  3512. CDigitClass digitTrain;
  3513. CRecvChar dlg;
  3514. CString m_class;
  3515. for(int k=0;k<CharacterNum;k++)
  3516. {
  3517. for(int x=0;x<cDatas[k].h;x++)
  3518. {
  3519. for(int y=0;y<cDatas[k].w;y++)
  3520. {
  3521. digitTrain.digitarray[y][x]=cDatas[k].Img[x*cDatas[k].w+y];
  3522. }
  3523. }
  3524. digitTrain.digitHeight=cDatas[k].h;
  3525. digitTrain.digitWidth=cDatas[k].w;
  3526. digitTrain.FixSize();
  3527. digitTrain.BinaryDigit();
  3528. digitTrain.ThinDigit_1();
  3530. pchview->m_result.SetSel(0,-1);
  3531. pchview->m_result.ReplaceSel("");
  3533. CString strInfo, strTemp;
  3534. strInfo.Format("");
  3535. int i,j;
  3536. for(i=0; i<digitTrain.digitHeight; i++)
  3537. {
  3538. for(j=0; j<digitTrain.digitWidth; j++)
  3539. {
  3540. if(digitTrain.digitarray[j][i] == 1)
  3541. strInfo += "■";
  3542. else
  3543. strInfo += "  ";
  3544. }
  3545. strInfo += "rn";
  3546. }
  3547. pchview->m_result.SetSel(0,-1);
  3548. pchview->m_result.ReplaceSel(strInfo);
  3551. if(dlg.DoModal()==IDOK)
  3552. {
  3553. m_class=dlg.m_char;
  3554. m_class.TrimLeft();
  3555. m_class.TrimRight();
  3556. digitTrain.GetFeature();
  3557. Sample temp;
  3558. for(i=0; i<13; i++)
  3559. temp.feature[i] = digitTrain.feature[i];
  3560. memcpy(temp.trueClass,(const char*) m_class,3*sizeof(char));
  3561. if(temp.trueClass[1]=='')
  3562. {
  3563. if (temp.trueClass[0]>='0'&&temp.trueClass[0]<='9')
  3564. {
  3565. temp.serialnum=numList.GetCount();
  3566. numList.AddTail(temp);
  3567. }
  3569. if (temp.trueClass[0]>='A'&&temp.trueClass[0]<='Z')
  3570. {
  3571. temp.serialnum=charList.GetCount();
  3572. charList.AddTail(temp);
  3573. }
  3575. }
  3576. else if(temp.trueClass[2]=='')
  3577. {
  3578. temp.serialnum=hanList.GetCount();
  3579. hanList.AddTail(temp);
  3580. }
  3581. }
  3583. }
  3588. }
  3590. void CDipView::OnSave() 
  3591. {
  3592. CFile char_file,num_file,han_file;
  3593. if(char_file.Open(".\chdata.ch",CFile::modeCreate|CFile::modeNoTruncate|CFile::modeReadWrite)==0)
  3594. {
  3595. AfxMessageBox("打开字符文件失败");
  3596. }
  3599. try
  3600. {
  3601. Sample sampleTemp;
  3602. char_file.SeekToEnd();
  3603. if(char_ConnectPos!=NULL)
  3604. {
  3605. while(char_ConnectPos!=NULL)
  3606. {
  3607. sampleTemp=charList.GetNext(char_ConnectPos);
  3608. char_file.Write(&sampleTemp, sizeof(sample));
  3609. }
  3610. char_ConnectPos=charList.GetTailPosition();
  3612. }
  3613. else
  3614. {
  3615. // Sample sampleTemp;
  3616. POSITION pos=charList.GetHeadPosition();
  3617. while(pos!=NULL)
  3618. {
  3619. sampleTemp=charList.GetNext(pos);
  3620. char_file.Write(&sampleTemp, sizeof(sample));
  3621. }
  3622. char_ConnectPos=charList.GetTailPosition();
  3624. }
  3626. }
  3627. catch(CFileException *e){
  3628. //Alert user to error
  3629. e->Delete();
  3630. }//end try-catch
  3631. char_file.Close();
  3633. if(num_file.Open(".\numdata.num",CFile::modeCreate|CFile::modeNoTruncate|CFile::modeReadWrite)==0)
  3634. {
  3635. AfxMessageBox("打开数字文件失败");
  3636. }
  3638. try
  3639. {
  3640. Sample sampleTemp;
  3641. num_file.SeekToEnd();
  3642. if(num_ConnectPos!=NULL)
  3643. {
  3644. while(num_ConnectPos!=NULL)
  3645. {
  3646. sampleTemp=numList.GetNext(num_ConnectPos);
  3647. num_file.Write(&sampleTemp, sizeof(sample));
  3648. }
  3649. num_ConnectPos=numList.GetTailPosition();
  3651. }
  3652. else
  3653. {
  3654. // Sample sampleTemp;
  3655. POSITION pos=numList.GetHeadPosition();
  3656. while(pos!=NULL)
  3657. {
  3658. sampleTemp=numList.GetNext(pos);
  3659. num_file.Write(&sampleTemp, sizeof(sample));
  3660. }
  3661. num_ConnectPos=numList.GetTailPosition();
  3663. }
  3665. }
  3666. catch(CFileException *e){
  3667. //Alert user to error
  3668. e->Delete();
  3669. }//end try-catch
  3670. num_file.Close();
  3673. if(han_file.Open(".\handata.han",CFile::modeCreate|CFile::modeNoTruncate|CFile::modeReadWrite)==0)
  3674. {
  3675. AfxMessageBox("打开汉字文件失败");
  3676. }
  3678. try
  3679. {
  3680. Sample sampleTemp;
  3681. han_file.SeekToEnd();
  3682. if(han_ConnectPos!=NULL)
  3683. {
  3684. while(han_ConnectPos!=NULL)
  3685. {
  3686. sampleTemp=hanList.GetNext(han_ConnectPos);
  3687. han_file.Write(&sampleTemp, sizeof(sample));
  3688. }
  3689. han_ConnectPos=hanList.GetTailPosition();
  3691. }
  3692. else
  3693. {
  3694. // Sample sampleTemp;
  3695. POSITION pos=hanList.GetHeadPosition();
  3696. while(pos!=NULL)
  3697. {
  3698. sampleTemp=hanList.GetNext(pos);
  3699. han_file.Write(&sampleTemp, sizeof(sample));
  3700. }
  3701. han_ConnectPos=hanList.GetTailPosition();
  3703. }
  3705. }
  3706. catch(CFileException *e){
  3707. //Alert user to error
  3708. e->Delete();
  3709. }//end try-catch
  3710. han_file.Close();
  3712. }
  3714. void CDipView::OnLoad() 
  3715. {
  3716. if (charList.GetCount()!=0||numList.GetCount()!=0||hanList.GetCount()!=0)
  3717. {
  3718. int nReturn;
  3719. nReturn=
  3720. AfxMessageBox("训练样本不为空,现在读取会覆盖这些样本,确认覆盖吗？",MB_YESNO);
  3721. if(nReturn==IDNO)
  3722. return;
  3723. }
  3724. charList.RemoveAll();
  3725. numList.RemoveAll();
  3726. hanList.RemoveAll();
  3728. CFile char_file,num_file,han_file;
  3729. if(char_file.Open(".\chdata.ch",CFile::modeCreate|CFile::modeNoTruncate|CFile::modeReadWrite)==0)
  3730. {
  3731. AfxMessageBox("打开字符文件失败，您最好退出程序");
  3732. }
  3734. try
  3735. {
  3736. Sample sampleTemp;
  3737. char_file.SeekToBegin();
  3738. DWORD dwBytesRemaining=char_file.GetLength();
  3739. while(dwBytesRemaining)
  3740. {
  3741. UINT nBytesRead=char_file.Read(&sampleTemp,sizeof(sample));
  3742. charList.AddTail(sampleTemp);
  3743. dwBytesRemaining-=nBytesRead;
  3744. }
  3745. CString s;
  3746. s.Format("共读取字符数据%d个",charList.GetCount());
  3747. AfxMessageBox(s);
  3750. }
  3751. catch(CFileException *e){
  3753. e->ReportError();
  3754. e->Delete();
  3755. }
  3756. char_file.Close();
  3757. char_ConnectPos=charList.GetTailPosition();
  3759. if(num_file.Open(".\numdata.num",CFile::modeCreate|CFile::modeNoTruncate|CFile::modeReadWrite)==0)
  3760. {
  3761. AfxMessageBox("打开数字文件失败，您最好退出程序");
  3762. }
  3764. try
  3765. {
  3766. Sample sampleTemp;
  3767. num_file.SeekToBegin();
  3768. DWORD dwBytesRemaining=num_file.GetLength();
  3769. while(dwBytesRemaining)
  3770. {
  3771. UINT nBytesRead=num_file.Read(&sampleTemp,sizeof(sample));
  3772. numList.AddTail(sampleTemp);
  3773. dwBytesRemaining-=nBytesRead;
  3774. }
  3775. CString s;
  3776. s.Format("共读取数字数据%d个",numList.GetCount());
  3777. AfxMessageBox(s);
  3780. }
  3781. catch(CFileException *e){
  3783. e->ReportError();
  3784. e->Delete();
  3785. }
  3786. num_file.Close();
  3787. num_ConnectPos=numList.GetTailPosition();
  3789. if(han_file.Open(".\handata.han",CFile::modeCreate|CFile::modeNoTruncate|CFile::modeReadWrite)==0)
  3790. {
  3791. AfxMessageBox("打开汉字文件失败，您最好退出程序");
  3792. }
  3794. try
  3795. {
  3796. Sample sampleTemp;
  3797. han_file.SeekToBegin();
  3798. DWORD dwBytesRemaining=han_file.GetLength();
  3799. while(dwBytesRemaining)
  3800. {
  3801. UINT nBytesRead=han_file.Read(&sampleTemp,sizeof(sample));
  3802. hanList.AddTail(sampleTemp);
  3803. dwBytesRemaining-=nBytesRead;
  3804. }
  3805. CString s;
  3806. s.Format("共读取汉字数据%d个",hanList.GetCount());
  3807. AfxMessageBox(s);
  3810. }
  3811. catch(CFileException *e){
  3813. e->ReportError();
  3814. e->Delete();
  3815. }
  3816. han_file.Close();
  3817. han_ConnectPos=hanList.GetTailPosition();
  3819. }
  3821. void CDipView::OnChange() 
  3822. {
  3824. CFile cf;
  3825. if(cf.Open(".\mydata.dat",CFile::modeCreate|CFile::modeNoTruncate|CFile::modeReadWrite)==0)
  3826. {
  3827. AfxMessageBox("打开文件失败，n您最好退出程序");
  3828. }
  3831. try
  3832. {
  3833. Sample sampleTemp;
  3834. cf.SeekToBegin();
  3835. DWORD dwBytesRemaining=cf.GetLength();
  3836. while(dwBytesRemaining)
  3837. {
  3838. UINT nBytesRead=cf.Read(&sampleTemp,sizeof(sample));
  3839. if(sampleTemp.trueClass[1]=='')
  3840. {
  3841. if (sampleTemp.trueClass[0]>='0'&&sampleTemp.trueClass[0]<='9')
  3842. {
  3844. numList.AddTail(sampleTemp);
  3845. }
  3847. if (sampleTemp.trueClass[0]>='A'&&sampleTemp.trueClass[0]<='Z')
  3848. {
  3849. sampleTemp.serialnum=sampleTemp.serialnum-32;
  3850. charList.AddTail(sampleTemp);
  3851. }
  3853. }
  3854. else if(sampleTemp.trueClass[2]=='')
  3855. {
  3856. sampleTemp.serialnum=sampleTemp.serialnum-48;
  3857. hanList.AddTail(sampleTemp);
  3858. }
  3859. dwBytesRemaining-=nBytesRead;
  3860. }
  3861. // CString s;
  3862. // s.Format("共读取数据%d个,最后一个的序列号为%d",
  3863. // TemplateList.GetCount(),sampleTemp.serialnum);
  3864. // AfxMessageBox(s);
  3867. }
  3868. catch(CFileException *e){
  3870. e->ReportError();
  3871. e->Delete();
  3872. }
  3873. cf.Close();
  3875. }
  3877. void CDipView::OnTrain() 
  3878. {
  3880. Sample temp;
  3881. Sample sa;
  3882. double target[36];
  3883. POSITION pos;
  3884. m_pdlg=new CTraining;
  3885. m_pdlg->Create(IDD_TRAINING);
  3886. m_pdlg->m_progress.SetRange32(0,2000);
  3887. m_pdlg->ShowWindow(SW_SHOW);
  3889. for (int num=0;num<2000;num++)
  3890. {
  3891. pos=charList.GetHeadPosition();
  3893. while(pos!=NULL)
  3894. {
  3895. temp=charList.GetNext(pos);
  3896. for(int i=0; i<13; i++)
  3897. sa.feature[i] = temp.feature[i]/15;
  3898. double eo=1,eh=1;
  3899. if(temp.trueClass[1]=='')
  3900. {
  3901. memset(target,0,36*sizeof(double));
  3902. target[map(&temp.trueClass[0])-10]=1;
  3903. char_bp.Train(sa.feature,13,target,26,&eo,&eh);
  3904. }
  3906. }
  3907. m_pdlg->m_output.SetWindowText("正在训练字母网络，请稍候……");
  3908. m_pdlg->m_progress.SetPos(num);
  3910. }
  3911. // bp.Save("bp.dat");
  3912. for (num=0;num<2000;num++)
  3913. {
  3914. pos=numList.GetHeadPosition();
  3915. while(pos!=NULL)
  3916. {
  3917. temp=numList.GetNext(pos);
  3918. for(int i=0; i<13; i++)
  3919. sa.feature[i] = temp.feature[i]/15;
  3920. double eo=1,eh=1;
  3921. if(temp.trueClass[1]=='')
  3922. {
  3923. memset(target,0,36*sizeof(double));
  3924. target[map(&temp.trueClass[0])]=1;
  3925. num_bp.Train(sa.feature,13,target,10,&eo,&eh);
  3926. }
  3928. }
  3929. m_pdlg->m_output.SetWindowText("正在训练数字网络，请稍候……");
  3930. m_pdlg->m_progress.SetPos(num);
  3932. }
  3933. // bp.Save("bp.dat");
  3934. for (num=0;num<2000;num++)
  3935. {
  3936. pos=charList.GetHeadPosition();
  3938. while(pos!=NULL)
  3939. {
  3940. temp=charList.GetNext(pos);
  3941. for(int i=0; i<13; i++)
  3942. sa.feature[i] = temp.feature[i]/15;
  3943. double eo=1,eh=1;
  3944. if(temp.trueClass[1]=='')
  3945. {
  3946. memset(target,0,36*sizeof(double));
  3947. target[map(&temp.trueClass[0])]=1;
  3948. char_num_bp.Train(sa.feature,13,target,36,&eo,&eh);
  3949. }
  3951. }
  3953. pos=numList.GetHeadPosition();
  3954. while(pos!=NULL)
  3955. {
  3956. temp=numList.GetNext(pos);
  3957. for(int i=0; i<13; i++)
  3958. sa.feature[i] = temp.feature[i]/15;
  3959. double eo=1,eh=1;
  3960. if(temp.trueClass[1]=='')
  3961. {
  3962. memset(target,0,36*sizeof(double));
  3963. target[map(&temp.trueClass[0])]=1;
  3964. char_num_bp.Train(sa.feature,13,target,36,&eo,&eh);
  3965. }
  3967. }
  3968. m_pdlg->m_output.SetWindowText("正在训练字母、数字网络，请稍候……");
  3969. m_pdlg->m_progress.SetPos(num);
  3971. }
  3972. // bp.Save("bp.dat");
  3973. m_pdlg->DestroyWindow();
  3974. }
  3977. int CDipView::map(char *character)
  3978. {
  3979. if(*character>='0'&&*character<='9')
  3980. return *character-'0';
  3981. if(*character>='A'&&*character<='Z')
  3982. return *character-'A'+10;
  3983. }
  3985. void CDipView::OnSavebp() 
  3986. {
  3987. char_bp.Save(".\char.bp");
  3988. num_bp.Save(".\num.bp");
  3989. char_num_bp.Save(".\char_num.bp");
  3991. }
  3993. void CDipView::OnReadbp() 
  3994. {
  3995. char_bp.Read(".\char.bp");
  3996. num_bp.Read(".\num.bp");
  3997. char_num_bp.Read(".\char_num.bp");
  3999. }