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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图象 > 查看源码
dibview.cpp
资源名称：lab5.rar　[点击查看]
上传用户：zbjinju
上传日期：2022-07-30
资源大小：11893k
文件大小：20k
源码类别：

图形图象

开发平台：

Visual C++

dibview.cpp：源码内容
  1. // dibview.cpp : implementation of the CDibView class
  2. //
  3. // This is a part of the Microsoft Foundation Classes C++ library.
  4. // Copyright (C) 1992-1998 Microsoft Corporation
  5. // All rights reserved.
  6. //
  7. // This source code is only intended as a supplement to the
  8. // Microsoft Foundation Classes Reference and related
  9. // electronic documentation provided with the library.
  10. // See these sources for detailed information regarding the
  11. // Microsoft Foundation Classes product.
  13. #include "stdafx.h"
  14. #include "diblook.h"
  16. #include "dibdoc.h"
  17. #include "dibview.h"
  18. #include "dibapi.h"
  19. #include "mainfrm.h"
  21. #include "BitmapInfoDlg.h"
  22. #include "ThresholdDialog.h"
  24. #include "DlgHistogram.h"
  25. #include <math.h>
  26. #define PI 3.14
  28. #include <afxtempl.h>
  30. #ifdef _DEBUG
  31. #undef THIS_FILE
  32. static char BASED_CODE THIS_FILE[] = __FILE__;
  33. #endif
  35. #define BEGIN_PROCESSING() INCEPUT_PRELUCRARI()
  37. #define END_PROCESSING(Title) SFARSIT_PRELUCRARI(Title)
  39. #define INCEPUT_PRELUCRARI() 
  40. CDibDoc* pDocSrc=GetDocument();
  41. CDocTemplate* pDocTemplate=pDocSrc->GetDocTemplate();
  42. CDibDoc* pDocDest=(CDibDoc*) pDocTemplate->CreateNewDocument();
  43. BeginWaitCursor();
  44. HDIB hBmpSrc=pDocSrc->GetHDIB();
  45. HDIB hBmpDest = (HDIB)::CopyHandle((HGLOBAL)hBmpSrc);
  46. if ( hBmpDest==0 ) {
  47. pDocTemplate->RemoveDocument(pDocDest);
  48. return;
  49. }
  50. BYTE* lpD = (BYTE*)::GlobalLock((HGLOBAL)hBmpDest);
  51. BYTE* lpS = (BYTE*)::GlobalLock((HGLOBAL)hBmpSrc);
  52. int iColors = DIBNumColors((char *)&(((LPBITMAPINFO)lpD)->bmiHeader)); 
  53. RGBQUAD *bmiColorsDst = ((LPBITMAPINFO)lpD)->bmiColors;
  54. RGBQUAD *bmiColorsSrc = ((LPBITMAPINFO)lpS)->bmiColors;
  55. BYTE * lpDst = (BYTE*)::FindDIBBits((LPSTR)lpD);
  56. BYTE * lpSrc = (BYTE*)::FindDIBBits((LPSTR)lpS);
  57. int dwWidth  = ::DIBWidth((LPSTR)lpS);
  58. int dwHeight = ::DIBHeight((LPSTR)lpS);
  59. int w=WIDTHBYTES(dwWidth*((LPBITMAPINFOHEADER)lpS)->biBitCount);
  61. #define BEGIN_SOURCE_PROCESSING 
  62. CDibDoc* pDocSrc=GetDocument();
  63. BeginWaitCursor();
  64. HDIB hBmpSrc=pDocSrc->GetHDIB();
  65. BYTE* lpS = (BYTE*)::GlobalLock((HGLOBAL)hBmpSrc);
  66. int iColors = DIBNumColors((char *)&(((LPBITMAPINFO)lpS)->bmiHeader)); 
  67. RGBQUAD *bmiColorsSrc = ((LPBITMAPINFO)lpS)->bmiColors;
  68. BYTE * lpSrc = (BYTE*)::FindDIBBits((LPSTR)lpS);
  69. int dwWidth  = ::DIBWidth((LPSTR)lpS);
  70. int dwHeight = ::DIBHeight((LPSTR)lpS);
  71. int w=WIDTHBYTES(dwWidth*((LPBITMAPINFOHEADER)lpS)->biBitCount);
  75. #define END_SOURCE_PROCESSING
  76. ::GlobalUnlock((HGLOBAL)hBmpSrc);
  77.     EndWaitCursor();
  78. /////////////////////////////////////////////////////////////////////////////
  81. //---------------------------------------------------------------
  82. #define SFARSIT_PRELUCRARI(Titlu)
  83. ::GlobalUnlock((HGLOBAL)hBmpDest);
  84. ::GlobalUnlock((HGLOBAL)hBmpSrc);
  85.     EndWaitCursor();
  86. pDocDest->SetHDIB(hBmpDest);
  87. pDocDest->InitDIBData();
  88. pDocDest->SetTitle((LPCSTR)Titlu);
  89. CFrameWnd* pFrame=pDocTemplate->CreateNewFrame(pDocDest,NULL);
  90. pDocTemplate->InitialUpdateFrame(pFrame,pDocDest);
  91. /////////////////////////////////////////////////////////////////////////////
  92. // CDibView
  94. IMPLEMENT_DYNCREATE(CDibView, CScrollView)
  96. BEGIN_MESSAGE_MAP(CDibView, CScrollView)
  97. //{{AFX_MSG_MAP(CDibView)
  98. ON_COMMAND(ID_EDIT_COPY, OnEditCopy)
  99. ON_UPDATE_COMMAND_UI(ID_EDIT_COPY, OnUpdateEditCopy)
  100. ON_COMMAND(ID_EDIT_PASTE, OnEditPaste)
  101. ON_UPDATE_COMMAND_UI(ID_EDIT_PASTE, OnUpdateEditPaste)
  102. ON_MESSAGE(WM_DOREALIZE, OnDoRealize)
  103. ON_COMMAND(ID_PROCESSING_PARCURGERESIMPLA, OnProcessingParcurgereSimpla)
  104. //}}AFX_MSG_MAP
  106. // Standard printing commands
  107. ON_COMMAND(ID_FILE_PRINT, CScrollView::OnFilePrint)
  108. ON_COMMAND(ID_FILE_PRINT_PREVIEW, CScrollView::OnFilePrintPreview)
  109. ON_COMMAND(ID_PROCESSING_LAB1, &CDibView::OnProcessingLab1)
  110. ON_COMMAND(ID_LAB2_BMPINFO, &CDibView::OnLab2Bmpinfo)
  111. ON_COMMAND(ID_LAB2_GRAYSCALE, &CDibView::OnLab2Grayscale)
  112. ON_COMMAND(ID_LAB2_SORTPALETTE, &CDibView::OnLab2Sortpalette)
  113. ON_COMMAND(ID_LAB2_BLACKANDWHITE, &CDibView::OnLab2Blackandwhite)
  114. ON_COMMAND(ID_LAB3_HISTOGRAM, &CDibView::OnLab3Histogram)
  115. ON_COMMAND(ID_LAB3_PRAGURIMULTIPLE, &CDibView::OnLab3Pragurimultiple)
  116. ON_COMMAND(ID_LAB3_PRAGURIMULTIPLEFS, &CDibView::OnLab3Pragurimultiplefs)
  117. ON_COMMAND(ID_PROCESSING_LAB4, &CDibView::OnProcessingLab4)
  118. ON_COMMAND(ID_PROCESSING_LAB5, &CDibView::OnProcessingLab5)
  119. END_MESSAGE_MAP()
  121. /////////////////////////////////////////////////////////////////////////////
  122. // CDibView construction/destruction
  124. CDibView::CDibView()
  125. {
  126. }
  128. CDibView::~CDibView()
  129. {
  130. }
  132. /////////////////////////////////////////////////////////////////////////////
  133. // CDibView drawing
  135. void CDibView::OnDraw(CDC* pDC)
  136. {
  137. CDibDoc* pDoc = GetDocument();
  139. HDIB hDIB = pDoc->GetHDIB();
  140. if (hDIB != NULL)
  141. {
  142. LPSTR lpDIB = (LPSTR) ::GlobalLock((HGLOBAL) hDIB);
  143. int cxDIB = (int) ::DIBWidth(lpDIB);         // Size of DIB - x
  144. int cyDIB = (int) ::DIBHeight(lpDIB);        // Size of DIB - y
  145. ::GlobalUnlock((HGLOBAL) hDIB);
  146. CRect rcDIB;
  147. rcDIB.top = rcDIB.left = 0;
  148. rcDIB.right = cxDIB;
  149. rcDIB.bottom = cyDIB;
  150. CRect rcDest;
  151. if (pDC->IsPrinting())   // printer DC
  152. {
  153. // get size of printer page (in pixels)
  154. int cxPage = pDC->GetDeviceCaps(HORZRES);
  155. int cyPage = pDC->GetDeviceCaps(VERTRES);
  156. // get printer pixels per inch
  157. int cxInch = pDC->GetDeviceCaps(LOGPIXELSX);
  158. int cyInch = pDC->GetDeviceCaps(LOGPIXELSY);
  160. //
  161. // Best Fit case -- create a rectangle which preserves
  162. // the DIB's aspect ratio, and fills the page horizontally.
  163. //
  164. // The formula in the "->bottom" field below calculates the Y
  165. // position of the printed bitmap, based on the size of the
  166. // bitmap, the width of the page, and the relative size of
  167. // a printed pixel (cyInch / cxInch).
  168. //
  169. rcDest.top = rcDest.left = 0;
  170. rcDest.bottom = (int)(((double)cyDIB * cxPage * cyInch)
  171. / ((double)cxDIB * cxInch));
  172. rcDest.right = cxPage;
  173. }
  174. else   // not printer DC
  175. {
  176. rcDest = rcDIB;
  177. }
  178. ::PaintDIB(pDC->m_hDC, &rcDest, pDoc->GetHDIB(),
  179. &rcDIB, pDoc->GetDocPalette());
  180. }
  181. }
  183. /////////////////////////////////////////////////////////////////////////////
  184. // CDibView printing
  186. BOOL CDibView::OnPreparePrinting(CPrintInfo* pInfo)
  187. {
  188. // default preparation
  189. return DoPreparePrinting(pInfo);
  190. }
  192. /////////////////////////////////////////////////////////////////////////////
  193. // CDibView commands
  196. LRESULT CDibView::OnDoRealize(WPARAM wParam, LPARAM)
  197. {
  198. ASSERT(wParam != NULL);
  199. CDibDoc* pDoc = GetDocument();
  200. if (pDoc->GetHDIB() == NULL)
  201. return 0L;  // must be a new document
  203. CPalette* pPal = pDoc->GetDocPalette();
  204. if (pPal != NULL)
  205. {
  206. CMainFrame* pAppFrame = (CMainFrame*) AfxGetApp()->m_pMainWnd;
  207. ASSERT_KINDOF(CMainFrame, pAppFrame);
  209. CClientDC appDC(pAppFrame);
  210. // All views but one should be a background palette.
  211. // wParam contains a handle to the active view, so the SelectPalette
  212. // bForceBackground flag is FALSE only if wParam == m_hWnd (this view)
  213. CPalette* oldPalette = appDC.SelectPalette(pPal, ((HWND)wParam) != m_hWnd);
  215. if (oldPalette != NULL)
  216. {
  217. UINT nColorsChanged = appDC.RealizePalette();
  218. if (nColorsChanged > 0)
  219. pDoc->UpdateAllViews(NULL);
  220. appDC.SelectPalette(oldPalette, TRUE);
  221. }
  222. else
  223. {
  224. TRACE0("tSelectPalette failed in CDibView::OnPaletteChangedn");
  225. }
  226. }
  228. return 0L;
  229. }
  231. void CDibView::OnInitialUpdate()
  232. {
  233. CScrollView::OnInitialUpdate();
  234. ASSERT(GetDocument() != NULL);
  236. SetScrollSizes(MM_TEXT, GetDocument()->GetDocSize());
  237. }
  240. void CDibView::OnActivateView(BOOL bActivate, CView* pActivateView,
  241. CView* pDeactiveView)
  242. {
  243. CScrollView::OnActivateView(bActivate, pActivateView, pDeactiveView);
  245. if (bActivate)
  246. {
  247. ASSERT(pActivateView == this);
  248. OnDoRealize((WPARAM)m_hWnd, 0);   // same as SendMessage(WM_DOREALIZE);
  249. }
  250. }
  252. void CDibView::OnEditCopy()
  253. {
  254. CDibDoc* pDoc = GetDocument();
  255. // Clean clipboard of contents, and copy the DIB.
  257. if (OpenClipboard())
  258. {
  259. BeginWaitCursor();
  260. EmptyClipboard();
  261. SetClipboardData (CF_DIB, CopyHandle((HANDLE) pDoc->GetHDIB()) );
  262. CloseClipboard();
  263. EndWaitCursor();
  264. }
  265. }
  269. void CDibView::OnUpdateEditCopy(CCmdUI* pCmdUI)
  270. {
  271. pCmdUI->Enable(GetDocument()->GetHDIB() != NULL);
  272. }
  275. void CDibView::OnEditPaste()
  276. {
  277. HDIB hNewDIB = NULL;
  279. if (OpenClipboard())
  280. {
  281. BeginWaitCursor();
  283. hNewDIB = (HDIB) CopyHandle(::GetClipboardData(CF_DIB));
  285. CloseClipboard();
  287. if (hNewDIB != NULL)
  288. {
  289. CDibDoc* pDoc = GetDocument();
  290. pDoc->ReplaceHDIB(hNewDIB); // and free the old DIB
  291. pDoc->InitDIBData();    // set up new size & palette
  292. pDoc->SetModifiedFlag(TRUE);
  294. SetScrollSizes(MM_TEXT, pDoc->GetDocSize());
  295. OnDoRealize((WPARAM)m_hWnd,0);  // realize the new palette
  296. pDoc->UpdateAllViews(NULL);
  297. }
  298. EndWaitCursor();
  299. }
  300. }
  303. void CDibView::OnUpdateEditPaste(CCmdUI* pCmdUI)
  304. {
  305. pCmdUI->Enable(::IsClipboardFormatAvailable(CF_DIB));
  306. }
  308. void CDibView::OnProcessingParcurgereSimpla() 
  309. {
  310. // TODO: Add your command handler code here
  311. BEGIN_PROCESSING();
  314. // Makes a grayscale image by equalizing the R, G, B components from the LUT
  315. for (int k=0;  k < iColors ; k++)
  316. bmiColorsDst[k].rgbRed=bmiColorsDst[k].rgbGreen=bmiColorsDst[k].rgbBlue=k;
  319. //  Goes through the bitmap pixels and performs their negative
  320. for (int i=0;i<dwHeight;i++)
  321. for (int j=0;j<dwWidth;j++)
  322.   {
  323. lpDst[i*w+j]= 255 - lpSrc[i*w+j]; //makes image negative
  324.   }
  327. END_PROCESSING("Operation name");
  328. }
  332. void CDibView::OnProcessingLab1()
  333. {
  334. BEGIN_PROCESSING();
  337. // Makes a grayscale image by equalizing the R, G, B components from the LUT
  338. //for (int k=0;  k < iColors ; k++)
  339. // bmiColorsDst[k].rgbRed=bmiColorsDst[k].rgbGreen=bmiColorsDst[k].rgbBlue=k;
  342. //  Goes through the bitmap pixels and performs their negative
  343. for (int i=0;i<dwHeight;i++)
  344. for (int j=0;j<dwWidth;j++)
  345.   {
  346. lpDst[i*w+j]= 5*lpSrc[i*w+j]; //makes image negative
  347.   }
  350. END_PROCESSING("lab1");
  351. }
  353. void CDibView::OnLab2Bmpinfo()
  354. {
  355. BEGIN_SOURCE_PROCESSING;
  356. CBitmapInfoDlg dlgBmpHeader;
  358. LPBITMAPINFO pBitmapInfoSrc = (LPBITMAPINFO) lpS;
  359. CString buffer;
  360. for(int i=0;i<iColors;++i){
  361. buffer.Format("%3d. t%3dt%3dt%3drn",i,
  362. bmiColorsSrc[i].rgbRed,
  363. bmiColorsSrc[i].rgbGreen,
  364. bmiColorsSrc[i].rgbBlue
  365. );
  366. dlgBmpHeader.m_TEXT += buffer;
  367. }
  369. dlgBmpHeader.DoModal();
  370. END_SOURCE_PROCESSING;
  371. }
  373. void CDibView::OnLab2Grayscale()
  374. {
  375. BEGIN_PROCESSING();
  377. for (int k=0;  k < iColors ; k++){
  378. int mean = (bmiColorsDst[k].rgbRed + bmiColorsDst[k].rgbGreen + bmiColorsDst[k].rgbBlue)/3;
  379. bmiColorsDst[k].rgbRed=bmiColorsDst[k].rgbGreen=bmiColorsDst[k].rgbBlue=mean;
  380. }
  382. END_PROCESSING("GrayScale");
  383. }
  385. void CDibView::OnLab2Sortpalette()
  386. {
  387. BEGIN_PROCESSING();
  389. int g[256];
  390. for (int k=0;  k < iColors ; k++){
  391. g[k] = bmiColorsDst[k].rgbRed;
  392. bmiColorsDst[k].rgbRed=bmiColorsDst[k].rgbGreen=bmiColorsDst[k].rgbBlue=k;
  393. }
  395. for (int i=0;i<dwHeight;i++)
  396. for (int j=0;j<dwWidth;j++){
  397. int k = lpDst[i*w+j];
  398. lpDst[i*w+j]= g[k];
  399. }
  401. END_PROCESSING("SortPalette");
  403. }
  405. void CDibView::OnLab2Blackandwhite()
  406. {
  407. BYTE threshold;
  408. CThresholdDialog dlgThresh;
  409. if(dlgThresh.DoModal() == IDOK){
  410. threshold = dlgThresh.m_thresh;
  411. BEGIN_PROCESSING();
  413. for (int i=0;i<dwHeight;i++)
  414. for (int j=0;j<dwWidth;j++){
  415. if(lpDst[i*w+j]<threshold)
  416. lpDst[i*w+j]=0;
  417. else
  418. lpDst[i*w+j]=255;
  419. }
  421. END_PROCESSING("Black and White");
  422. }
  423. }
  425. void CDibView::OnLab3Histogram()
  426. {
  427. BEGIN_SOURCE_PROCESSING; 
  429. int histValues[256]; 
  430. float FDPValues[256]; 
  432. int i,j;
  433. for(i=0;i<256;++i){
  434. histValues[i] = 0;
  435. }
  437. for (i=0;i<dwHeight;i++)
  438. for (j=0;j<dwWidth;j++){
  439. histValues[ lpSrc[i*w+j] ]++;
  440. }
  442. int m = dwHeight * dwWidth;
  443. for(i=0; i<256; ++i)
  444. FDPValues[i] = (float)histValues[i]/m;
  446. CDlgHistogram dlg;
  447. memcpy(dlg.m_Histogram.values,histValues,sizeof(histValues)); 
  448. dlg.DoModal();
  450. END_SOURCE_PROCESSING;
  451. }
  453. //cel mai apropiat maxim din histograma
  454. int closestMax(int *h, int k){
  455. if(h[k] == 1)
  456. return k;
  457. int i=1;
  458. while(1){
  459. if(k-i < 0)
  460. return 0;
  461. if(h[k-i] == 1)
  462. return k-i;
  463. if(k+i > 255)
  464. return 255;
  465. if(h[k+i] == 1)
  466. return k+i;
  468. ++i;
  469. }
  470. return 0;
  471. }
  473. void CDibView::OnLab3Pragurimultiple()
  474. {
  475. int WH = 5;
  476. double TH = 0.0003;
  478. BEGIN_PROCESSING();
  480. int histValues[256]; 
  481. double FDPValues[256]; 
  482. int locMax[256];
  485. int i,j,k;
  486. for(i=0;i<256;++i){
  487. histValues[i] = 0;
  488. locMax[i] = 0;
  489. }
  491. for (i=0;i<dwHeight;i++)
  492. for (j=0;j<dwWidth;j++){
  493. histValues[ lpSrc[i*w+j] ]++;
  494. }
  496. int m = dwHeight * dwWidth;
  497. for(i=0; i<256; ++i)
  498. FDPValues[i] = (double)histValues[i]/m;
  500. //calculam maximii locali
  501. //locMax va contine 1 dak e maxim local, 0 dak nu
  503. locMax[0] = 1;
  504. locMax[255] = 1;
  506. double v,max;
  507. for(k = WH; k<255-WH; ++k){
  508. v = 0.0;
  509. max = FDPValues[k-WH];
  510. for(i=k-WH; i<=k+WH; ++i){
  511. v += FDPValues[i];
  512. if(FDPValues[i] > max)
  513. max = FDPValues[i];
  514. }
  515. v /= (2*WH+1);
  516. if(FDPValues[k] > v+TH && FDPValues[k] >= max)
  517. locMax[k] = 1;
  518. }
  520. for (int i=0;i<dwHeight;i++)
  521. for (int j=0;j<dwWidth;j++){
  522. k = lpDst[i*w+j];
  523. lpDst[i*w+j]= closestMax(locMax, k);
  524. }
  528. END_PROCESSING("PraguriMultiple");
  529. }
  531. void CDibView::OnLab3Pragurimultiplefs()
  532. {
  533. int WH = 5;
  534. double TH = 0.0003;
  536. BEGIN_PROCESSING();
  538. int histValues[256]; 
  539. double FDPValues[256]; 
  540. int locMax[256];
  543. int i,j,k;
  544. for(i=0;i<256;++i){
  545. histValues[i] = 0;
  546. locMax[i] = 0;
  547. }
  549. for (i=0;i<dwHeight;i++)
  550. for (j=0;j<dwWidth;j++){
  551. histValues[ lpSrc[i*w+j] ]++;
  552. }
  554. int m = dwHeight * dwWidth;
  555. for(i=0; i<256; ++i)
  556. FDPValues[i] = (double)histValues[i]/m;
  558. //calculam maximii locali
  559. //locMax va contine 1 dak e maxim local, 0 dak nu
  561. locMax[0] = 1;
  562. locMax[255] = 1;
  564. double v,max;
  565. for(k = WH; k<255-WH; ++k){
  566. v = 0.0;
  567. max = FDPValues[k-WH];
  568. for(i=k-WH; i<=k+WH; ++i){
  569. v += FDPValues[i];
  570. if(FDPValues[i] > max)
  571. max = FDPValues[i];
  572. }
  573. v /= (2*WH+1);
  574. if(FDPValues[k] > v+TH && FDPValues[k] >= max)
  575. locMax[k] = 1;
  576. }
  578. int pixel_vechi, pixel_nou, eroare;
  579. BYTE *source = (BYTE*)malloc(dwHeight * w * sizeof BYTE);
  580. memcpy(source, lpSrc, dwHeight * w * sizeof BYTE);
  581. for (int i=0;i<dwHeight;i++)
  582. for (int j=0;j<dwWidth;j++){
  583. pixel_vechi = source[i*w+j];
  584. pixel_nou = closestMax(locMax, pixel_vechi);
  585. lpDst[i*w+j] = pixel_nou;
  586. eroare = pixel_vechi - pixel_nou;
  587. if(j < (dwWidth-1))
  588. source[i*w+(j+1)] = source[i*w+(j+1)] + 7*eroare/16;
  589. if(i < (dwHeight-1) && j > 0)
  590. source[(i+1)*w+(j-1)] = source[(i+1)*w+(j-1)] + 3*eroare/16;
  591. if(i < (dwHeight-1))
  592. source[(i+1)*w+j] = source[(i+1)*w+j] + 5*eroare/16;
  593. if(i < (dwHeight-1) && j < (dwWidth-1))
  594. source[(i+1)*w+(j+1)] = source[(i+1)*w+(j+1)] + eroare/16;
  595. }
  597. free(source);
  599. END_PROCESSING("PraguriMultiple Floyd-Steinberg");
  600. }
  602. void CDibView::OnProcessingLab4()
  603. {
  604.   BEGIN_PROCESSING();
  605.   CString msg;
  606.   int A=0;
  607.   int perim=0;
  608.   int ri,ci;
  609.   float arct;
  610.   float x,y,x1,x2;
  611.   int E=0;
  612.   int cmax=0,cmin=1000,rmax=0,rmin=1000;
  613.   //aria
  614.   for (int i=1;i<dwHeight;i++)
  615.   for (int j=1;j<dwWidth; j++)
  616.   if (lpSrc[i*w+j]==0){
  617.                  A++;
  618.                 
  619.   }
  620.  //centrul de masa 
  621.   int xcm=0;
  622.   for(int i=0;i<dwHeight;i++)
  623.   for(int j=0; j<dwWidth;j++)
  624.   if (lpSrc[i*w+j]==0)
  625.   xcm=xcm+i;
  626.   xcm=xcm/A;
  627.   int ycm=0;
  628.   for(int i=0;i<dwHeight;i++)
  629.   for(int j=0; j<dwWidth;j++)
  630.   if (lpSrc[i*w+j]==0)
  631.   ycm=ycm+j;
  632.   ycm=ycm/A;
  634.   lpDst[xcm*w+ycm]=2;
  636.   bmiColorsDst[2].rgbRed=255;
  637.   bmiColorsDst[2].rgbGreen=0;
  638.   bmiColorsDst[2].rgbBlue=0;
  640.   bmiColorsDst[3].rgbRed=0;
  641.   bmiColorsDst[3].rgbGreen=255;
  642.   bmiColorsDst[3].rgbBlue=0;
  644.   bmiColorsDst[4].rgbRed=0;
  645.   bmiColorsDst[4].rgbGreen=0;
  646.   bmiColorsDst[4].rgbBlue=255;
  648.   //ungiul axei
  649.   for(int i=0;i<dwHeight;i++)
  650.   for(int j=0;j<dwWidth;j++)
  651.           if (lpSrc[i*w+j]==0)
  652.   {
  653.  y+=(i-xcm)*(j-ycm);
  654.  x1+=(j-ycm)*(j-ycm);
  655.  x2+=(i-xcm)*(i-xcm);
  656.   }
  658.    x=(x1-x2);
  659.    arct =(180.0*atan2(y,x))/(2.0*PI);
  661.    //Perimetrul
  662.     for(int i=0;i<dwHeight;i++)
  663.   for(int j=0;j<dwWidth;j++)
  664.   {
  665.   if (lpSrc[i*w+j]==0) 
  666.   {   
  667.   if  ((lpSrc[(i+1)*w+j]>1) || (lpSrc[(i-1)*w+j]>1) || (lpSrc[i*w+(j+1)]>1) || (lpSrc[i*w+(j-1)]>1))
  668.   {
  669.   perim++;
  670.   lpDst[i*w+j]=3;
  671.   }
  672.   }
  674.   }
  675.    // Elongatia
  676.    for(int i=0;i<dwHeight;i++)
  677.   for(int j=0;j<dwWidth;j++)
  678.          if (lpSrc[i*w+j]==0) 
  679.   { 
  680.              if (i>rmax)
  681.     rmax=i;
  682.  if (i<rmin)
  683.     rmin=i;
  684.  if (j>cmax)
  685.     cmax=j;
  686.  if (j<cmin)
  687.     cmin=j;
  689.   }
  690.    E=(cmax-cmin+1)/(rmax-rmin+1);
  691.  
  692.   // proiectiile
  694.    int *HX=new int[dwHeight];
  695.    memset(HX,0,sizeof(int)*dwHeight);
  696.     
  697.    int *WX=new int[dwWidth];
  698.    memset(WX,0,sizeof(int)*dwWidth);
  700.   for(int i=0;i<dwHeight;i++)
  701.   for(int j=0;j<dwWidth;j++)
  702.   {
  703.   if (lpSrc[i*w+j]==0) 
  704.            HX[i]++;
  705.   }
  707. for(int i=0;i<dwHeight;i++)
  708.   for(int j=0;j<dwWidth;j++)
  709.   {
  710.   if (lpSrc[i*w+j]==0) 
  711.            WX[j]++;
  712.   }
  714.    for (int i=0;i<dwHeight;i++)
  715.    for(int j=0;j<HX[i];j++)
  716.           lpDst[i*w+j]=4;
  718.  for (int i=0;i<dwWidth;i++)
  719.    for(int j=0;j<WX[i];j++)
  720.           lpDst[i*w+j]=4;
  722.   msg.Format("Aria %dnCm(%d,%d)nO%fnP %dnE %dn",A,xcm,ycm,arct,perim,E);
  723.   AfxMessageBox(msg);
  724.   END_PROCESSING("Calcule");
  725.  
  726. }
  728. void CDibView::OnProcessingLab5()
  729. {
  730. BEGIN_PROCESSING();
  731.    int  *Labels;
  732.    Labels=new int[w*dwHeight];
  733.    memset(Labels,0,(w*dwHeight)*sizeof(int));
  734.    int crtLabel=1;
  735.    CArray<CPoint,CPoint&>Pair;
  737.       for(int i=dwHeight-2;i>=0;i--)
  738.     for(int j=1;j<dwWidth-1;j++)
  739.  {
  740. int B=lpSrc[(i+1)*w+(j-1)];
  741.  int C=lpSrc[(i+1)*w+j];
  742.  int D=lpSrc[(i+1)*w+(j+1)];
  743.  int A=lpSrc[i*w+(j-1)];
  744.  int X=lpSrc[i*w+j];
  745.  if(X==0)
  746.  {
  747.               if(B!=0)
  748.   if(A!=0)
  749.   if(C!=0)
  750.   if(D!=0)
  751.   Labels[i*w+j]=crtLabel++;  
  752.   else
  753.   Labels[i*w+j]=Labels[(i+1)*w+(j+1)];  
  754.   else
  755.   Labels[i*w+j]=Labels[(i+1)*w+j];
  756.   else 
  757.   if(C!=0)
  758.   if(D!=0)
  759.   Labels[i*w+j]=Labels[i*w+(j-1)]; 
  760.   else
  761.   {
  762.   Labels[i*w+j]=Labels[i*w+(j-1)];
  763.   Pair.Add(CPoint(Labels[i*w + (j-1)],Labels[(i+1)*w+ (j+1)]));
  764.   }
  765.   else
  766.   {
  767.   Labels[i*w+j]=Labels[i*w+(j-1)];
  768.   Pair.Add(CPoint(Labels[i*w + (j-1)],Labels[(i+1)*w+ j])); 
  769.   }
  770.   else if ((C!=0)&&(D==0)) 
  771.          {
  772. Labels[i*w+j]=Labels[(i+1)*w+(j-1)];
  773. Pair.Add(CPoint(Labels[(i+1)*w+(j-1)],Labels[(i+1)*w+ (j+1)])); 
  774.      }
  775.   else
  776.     Labels[i*w+j]=Labels[(i+1)*w+(j-1)]; 
  777.            }
  778. }
  780. int *Union;
  781. Union =new int[crtLabel];
  782.         for (int i=0; i<crtLabel; i++)
  783. Union[i]=i;
  784. int la,lb;
  785.  
  786. for (int i=0; i<Pair.GetSize(); i++)
  787. {
  788. la=Union[Pair[i].x];
  789. lb=Union[Pair[i].y];
  790. for (int i=0;i<crtLabel;i++)
  791. {
  792. if (Union[i]==la)
  793. Union[i]=lb;
  794. }
  795. }
  798. for(int i=dwHeight-2;i>=0;i--)
  799. for(int j=1;j<dwWidth-1;j++){
  800. if(Labels[i*w+j]>0)
  801. lpDst[i*w+j]=Union[Labels[i*w+j]]%254+1;
  802. }
  803.              
  805. for(int k=1;k<254;k++){
  806. bmiColorsDst[k].rgbRed=rand()%255;
  807. bmiColorsDst[k].rgbBlue=rand()%255;
  808. bmiColorsDst[k].rgbGreen=rand()%255;
  809. }
  811. END_PROCESSING ("Lab 5");
  812. }