영상처리 이미지변환의 총정리

// DP_서정회Doc.cpp : implementation of the CDP_Doc class // #include "stdafx.h" #include "DP_서정회.h" #include "math.h" #include "DP_서정회Doc.h" #ifdef _DEBUG #define new DEBUG_NEW #undef THIS_FILE static char THIS_FILE[] = __FILE__; #endif ///////////////////////////////////////////////////////////////////////////// // CDP_Doc IMPLEMENT_DYNCREATE(CDP_Doc, CDocument) BEGIN_MESSAGE_MAP(CDP_Doc, CDocument) //{{AFX_MSG_MAP(CDP_Doc) // NOTE - the ClassWizard will add and remove mapping macros here. // DO NOT EDIT what you see in these blocks of generated code! //}}AFX_MSG_MAP END_MESSAGE_MAP() ///////////////////////////////////////////////////////////////////////////// // CDP_Doc construction/destruction CDP_Doc::CDP_Doc() { // TODO: add one-time construction code here } CDP_Doc::~CDP_Doc() { } BOOL CDP_Doc::OnNewDocument() { if (!CDocument::OnNewDocument()) return FALSE; // TODO: add reinitialization code here // (SDI documents will reuse this document) return TRUE; } ///////////////////////////////////////////////////////////////////////////// // CDP_Doc serialization void CDP_Doc::Serialize(CArchive& ar) { if (ar.IsStoring()) { // TODO: add storing code here ar.Write(m_ResultImg, 256 * 256); } else { // TODO: add loading code here CFile* fp = ar.GetFile(); //파일 포인터를 얻는다 if(fp->GetLength() != 256 * 256) { AfxMessageBox("256 * 256 이미지가 아닙니다."); return; } ar.Read(m_OpenImg, fp->GetLength()); } } ///////////////////////////////////////////////////////////////////////////// // CDP_Doc diagnostics #ifdef _DEBUG void CDP_Doc::AssertValid() const { CDocument::AssertValid(); } void CDP_Doc::Dump(CDumpContext& dc) const { CDocument::Dump(dc); } #endif //_DEBUG ///////////////////////////////////////////////////////////////////////////// // CDP_Doc commands void CDP_Doc::DoubleLoad() { AfxMessageBox(("첫번째 이미지 선택하시오")); CFileDialog dlg(TRUE); if(dlg.DoModal()==IDOK) { CFile file; file.Open(dlg.GetFileName(),CFile::modeRead); file.Read(m_ImageBuffer1,sizeof(m_ImageBuffer1)); file.Close(); } AfxMessageBox(("두번째 이미지 선택하시오")); CFileDialog dlg2(TRUE); if(dlg2.DoModal()==IDOK) { CFile file; file.Open(dlg2.GetFileName(),CFile::modeRead); file.Read(m_ImageBuffer2,sizeof(m_ImageBuffer2)); file.Close(); } } void CDP_Doc::Binary() { for(int i=0;i<256;i++) { for(int j=0;j<256;j++) { if(m_OpenImg[i][j]>=128) { m_ResultImg[i][j]=255; } else m_ResultImg[i][j]=0; } } } void CDP_Doc::Histo_Equal() { int k=0,sum=0,total_pixels=0; int hist[256]; int sum_of_hist[256]; for(int z=0;z<256;z++) { hist[z]=0; sum_of_hist[z]=0; } for(int i=0;i<256;i++) { for(int j=0;j<256;j++) { k=m_OpenImg[i][j]; hist[k]=hist[k]+1; } } for(i=0;i<256;i++) { sum=sum+hist[i]; sum_of_hist[i]=sum; } total_pixels=256*256; for(i=0;i<256;i++) { for(int j=0;j<256;j++) { k=m_OpenImg[i][j]; m_ResultImg[i][j]=sum_of_hist[k]*(255.0/total_pixels); } } } void CDP_Doc::Histo_Stretch() { int histogram[256]; unsigned char LUT[256]; int lowthresh, highthresh; float scale_factor; lowthresh=0; highthresh=255; for(int i=0;i<256;i++) //히스토그램 명도 값의 빈도 수 조사 histogram [i]=0; for(i=0;i<256;i++) { for(int j=0;j<256;j++) { histogram[(int)m_OpenImg[i][j]]++; } } for(i=0;i<256;i++) //히스토그램 중 가장 낮은 명도 값 계싼 { if(histogram[i]) { lowthresh=i; break; } } for(i=255;i>=0;i--) //히스토그램 중 가장 높은 명도 값 계산 { if(histogram[i]) { highthresh=i; break; } } for(i=0;i<lowthresh;i++) //룩업 테이블 계산 LUT[i]=0; for(i=255;i>highthresh;i--) LUT[i]=255; scale_factor=(float)(255.0/(highthresh-lowthresh)); //스트레칭의 강도를 조절할 수 있는 factor계산 for(i=lowthresh;i<=highthresh;i++) LUT[i]=(unsigned char)((i-lowthresh)*scale_factor); for(i=0;i<256;i++) //룩업 테이블에서 계산된 출력 값 할당 { for(int j=0;j<256;j++) { m_ResultImg[i][j]=LUT[(unsigned char)m_OpenImg[i][j]]; } } } void CDP_Doc::ZoomIn() { int i=0,j=0,y=0,x=0,sum=0; for(y=0;y<256;y++) { for(x=0;x<256;x++) { i=y*scale_y; j=x*scale_x; m_scaleImg[i][j]=m_OpenImg[y][x]; } } for(y=0;y<256*scale_y;y++) { for(x=0;x<256*scale_x;x+=scale_x) { for(int m=1;m<scale_x;m++) { m_scaleImg[y][x-m]=m_scaleImg[y][x]; } } } for(x=0;x<256*scale_x;x++) { for(y=0;y<256*scale_y;y+=scale_y) { for(int n=1;n<scale_y;n++) { m_scaleImg[y-n][x]=m_scaleImg[y][x]; } } } y=0; x=0; for(i=1;i<512;i+=2) { for(j=1;j<512;j+=2) { sum+=m_scaleImg[y][(x)]; sum+=m_scaleImg[y][(x)+2]; sum+=m_scaleImg[y+2][(x)]; sum+=m_scaleImg[y+2][(x)+2]; sum=sum/4; m_scaleImg[i][(j)]=sum; sum=0; x+=2; } x=0; y+=2; } y=0; x=0; for(j=0;j<512;j+=2) { for(i=1;i<512;i+=2) { sum+=m_scaleImg[y][(x)]; sum+=m_scaleImg[y+2][(x)]; sum=sum/2; m_scaleImg[i][(j)]=sum; sum=0; y+=2; } y=0; x+=2; } } void CDP_Doc::ZoomOut() { int i=0,j=0,y=0,x=0; for(y=0;y<256;y++) { for(x=0;x<256;x++) { i=y*scale_y2; j=x*scale_x2; m_scaleImg2[i][j]=m_OpenImg[y][x]; } } } void CDP_Doc::Mirror() { int x2=0; for(int y=0;y<256;y++) { for(int x=0;x<256;x++) { x2=(255-x); m_ResultImg[y][x2]=m_OpenImg[y][x]; } } } void CDP_Doc::Flip() { int y2=0; for(int y=0;y<256;y++) { for(int x=0;x<256;x++) { y2=(255-y); m_ResultImg[y2][x]=m_OpenImg[y][x]; } } } void CDP_Doc::Rotate() { double i=0, j=0, m=0, n=0, x2=0, y2=0; int y=0, x=0, left=0, right = 0; double seta =0, pi=3.1415926535; seta=pi/4; for(y=0; y<256; y++) { for (x=0; x<256; x++) { i=-1*((y-128)*sin(seta))+128; j=(x-128)*cos(seta); x2=(i+j); if(x2>255)x2=255; if(x2<0)x2=0; m=((x-128)*sin(seta))+128; n=(y-128)*cos(seta); y2=(m+n); if(y2>255)y2=255; if(y2<0)y2=0; m_ResultImg[(int)y2][(int)x2]=m_OpenImg[y][x]; //m_ResultImg[(int)y2][(int)x2]=m_ResultImg[y][x]; } } for(y=1; y<255; y++) { for (x=1; x<255; x++) { left=m_ResultImg[y][x-1]; right=m_ResultImg[y][x+1]; if(left!=0 || right!=0 || m_ResultImg[y][x] == 0) { m_ResultImg[y][x] = (left+right) /2; } } } } void CDP_Doc::Sobel() { int centerValue1=0; int centerValue2=0; int sum=0; unsigned char m_tmpImg[258][258]; int mask1[3][3]={-1,0,1,-2,0,2,-1,0,1}; int mask2[3][3]={1,2,1,0,0,0,-1,-2,-1}; for(int k=1;k<257;k++) { for(int l=1;l<257;l++) { m_tmpImg[k][l]=m_OpenImg[k-1][l-1]; } } for(int row=0;row<256;row++) { for(int column=0;column<256;column++) { for(int i=0;i<3;i++) { for(int j=0;j<3;j++) { centerValue1+=m_tmpImg[i+row][j+column]*mask1[i][j]; centerValue2+=m_tmpImg[i+row][j+column]*mask2[i][j]; } } sum=abs(centerValue1)+abs(centerValue2); if(sum>255)sum=255; m_ResultImg[row][column]=(unsigned char)sum; centerValue1=0; centerValue2=0; sum=0; } } }

------------------------------------

 // DP_서정회View.cpp : implementation of the CDP_View class // #include "stdafx.h" #include "DP_서정회.h" #include "DP_서정회Doc.h" #include "DP_서정회View.h" #ifdef _DEBUG #define new DEBUG_NEW #undef THIS_FILE static char THIS_FILE[] = __FILE__; #endif ///////////////////////////////////////////////////////////////////////////// // CDP_View IMPLEMENT_DYNCREATE(CDP_View, CScrollView) BEGIN_MESSAGE_MAP(CDP_View, CScrollView) //{{AFX_MSG_MAP(CDP_View) ON_COMMAND(ID_C_SUM, OnCSum) ON_COMMAND(ID_Image_SUM, OnImageSUM) ON_COMMAND(ID_C_SUB, OnCSub) ON_COMMAND(ID_C_MUL, OnCMul) ON_COMMAND(ID_C_DIV, OnCDiv) ON_COMMAND(ID_Image_SUB, OnImageSUB) ON_COMMAND(ID_Image_MUL, OnImageMUL) ON_COMMAND(ID_Image_DIV, OnImageDIV) ON_COMMAND(ID_BINARY, OnBinary) ON_COMMAND(ID_HISTOEQUAL, OnHistoequal) ON_COMMAND(ID_HISTOSTRETCH, OnHistostretch) ON_COMMAND(ID_ZOOMIN, OnZoomin) ON_COMMAND(ID_ZOOMOUT, OnZoomout) ON_COMMAND(ID_MIRROR, OnMirror) ON_COMMAND(ID_FILP, OnFilp) ON_COMMAND(ID_ROTATE, OnRotate) ON_COMMAND(ID_SOBEL, OnSobel) //}}AFX_MSG_MAP // Standard printing commands ON_COMMAND(ID_FILE_PRINT, CScrollView::OnFilePrint) ON_COMMAND(ID_FILE_PRINT_DIRECT, CScrollView::OnFilePrint) ON_COMMAND(ID_FILE_PRINT_PREVIEW, CScrollView::OnFilePrintPreview) END_MESSAGE_MAP() ///////////////////////////////////////////////////////////////////////////// // CDP_View construction/destruction CDP_View::CDP_View() { // TODO: add construction code here } CDP_View::~CDP_View() { } BOOL CDP_View::PreCreateWindow(CREATESTRUCT& cs) { // TODO: Modify the Window class or styles here by modifying // the CREATESTRUCT cs return CScrollView::PreCreateWindow(cs); } ///////////////////////////////////////////////////////////////////////////// // CDP_View drawing void CDP_View::OnDraw(CDC* pDC) { CDP_Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); // TODO: add draw code for native data here for(int y=0;y<256;y++) { for(int x=0;x<256;x++) { // pDC->SetPixel(x,y,RGB(m_ViewImg[y][x], m_ViewImg[y][x], m_ViewImg[y][x])); pDC->SetPixel(x,y,RGB(pDoc->m_OpenImg[y][x], pDoc->m_OpenImg[y][x], pDoc->m_OpenImg[y][x])); } } if(zoom !=1 && zoom !=2) { for(y=0;y<256;y++) { for(int x=0;x<256;x++) { pDC->SetPixel(x+300,y,RGB(pDoc->m_ResultImg[y][x],pDoc->m_ResultImg[y][x], pDoc->m_ResultImg[y][x])); } } } if(zoom==1) { for(int y=0;y<scale_y*256;y++) { for(int x=0;x<scale_x*256;x++) { pDC->SetPixel(x+300,y,RGB(pDoc->m_scaleImg[y][x],pDoc->m_scaleImg[y][x], pDoc->m_scaleImg[y][x])); } } zoom=0; } if(zoom==2) { for(int y=0;y<scale_y2*256;y++) { for(int x=0;x<scale_x2*256;x++) { pDC->SetPixel(x+300,y,RGB(pDoc->m_scaleImg2[y][x],pDoc->m_scaleImg2[y][x],pDoc->m_scaleImg2[y][x])); } } zoom=0; } } void CDP_View::OnInitialUpdate() { CScrollView::OnInitialUpdate(); CSize sizeTotal; // TODO: calculate the total size of this view sizeTotal.cx = sizeTotal.cy = 100; SetScrollSizes(MM_TEXT, sizeTotal); CDP_Doc* pDoc=GetDocument(); for(int x=0;x<256;x++) { for(int y=0;y<256;y++) { m_ViewImg[x][y]=pDoc->m_OpenImg[x][y]; } } } ///////////////////////////////////////////////////////////////////////////// // CDP_View printing BOOL CDP_View::OnPreparePrinting(CPrintInfo* pInfo) { // default preparation return DoPreparePrinting(pInfo); } void CDP_View::OnBeginPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/) { // TODO: add extra initialization before printing } void CDP_View::OnEndPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/) { // TODO: add cleanup after printing } ///////////////////////////////////////////////////////////////////////////// // CDP_View diagnostics #ifdef _DEBUG void CDP_View::AssertValid() const { CScrollView::AssertValid(); } void CDP_View::Dump(CDumpContext& dc) const { CScrollView::Dump(dc); } CDP_Doc* CDP_View::GetDocument() // non-debug version is inline { ASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CDP_Doc))); return (CDP_Doc*)m_pDocument; } #endif //_DEBUG ///////////////////////////////////////////////////////////////////////////// // CDP_View message handlers void CDP_View::OnCSum() { // TODO: Add your command handler code here int data=0; CDP_Doc* pDoc=GetDocument(); ASSERT_VALID(pDoc); for(int x=0;x<256;x++) { for(int y=0;y<256;y++) { data=pDoc->m_OpenImg[x][y]+100; if(data>255) { pDoc->m_ResultImg[x][y]=255; } else { pDoc->m_ResultImg[x][y]; } } } Invalidate(FALSE); } void CDP_View::OnImageSUM() { // TODO: Add your command handler code here int data=0; CDP_Doc* pDoc=GetDocument(); ASSERT_VALID(pDoc); pDoc->DoubleLoad(); for(int x=0;x<256;x++) { for(int y=0;y<256;y++) { data=pDoc->m_ImageBuffer1[x][y]+pDoc->m_ImageBuffer2[x][y]; if(data>255) { pDoc->m_ResultImg[x][y]=255; } else { pDoc->m_ResultImg[x][y]=data; } pDoc->m_OpenImg[x][y]=data; } } Invalidate(FALSE); } void CDP_View::OnCSub() { // TODO: Add your command handler code here int data=0; CDP_Doc* pDoc=GetDocument(); ASSERT_VALID(pDoc); for(int x=0;x<256;x++) { for(int y=0;y<256;y++) { data=pDoc->m_OpenImg[x][y]-100; if(data<0) { pDoc->m_ResultImg[x][y]=0; } else { pDoc->m_ResultImg[x][y]; } } } Invalidate(FALSE); } void CDP_View::OnCMul() { // TODO: Add your command handler code here int data=0; CDP_Doc* pDoc=GetDocument(); ASSERT_VALID(pDoc); for(int x=0;x<256;x++) { for(int y=0;y<256;y++) { data=pDoc->m_OpenImg[x][y]*3; if(data>255) { pDoc->m_ResultImg[x][y]=255; } else { pDoc->m_ResultImg[x][y]; } } } Invalidate(FALSE); } void CDP_View::OnCDiv() { // TODO: Add your command handler code here int data=0; CDP_Doc* pDoc=GetDocument(); ASSERT_VALID(pDoc); for(int x=0;x<256;x++) { for(int y=0;y<256;y++) { data=pDoc->m_OpenImg[x][y]/2; pDoc->m_ResultImg[x][y]=data; } } Invalidate(FALSE); } void CDP_View::OnImageSUB() { // TODO: Add your command handler code here int data1=0, data2=0; CDP_Doc* pDoc=GetDocument(); ASSERT_VALID(pDoc); pDoc->DoubleLoad(); for(int x=0;x<256;x++) { for(int y=0;y<256;y++) { data1=pDoc->m_ImageBuffer1[x][y]-pDoc->m_ImageBuffer2[x][y]; data2=abs(pDoc->m_ImageBuffer1[x][y]-pDoc->m_ImageBuffer2[x][y]); if(data1<0) { pDoc->m_OpenImg[x][y]; } else { pDoc->m_OpenImg[x][y]=data1; } pDoc->m_ResultImg[x][y]=data2; } } Invalidate(FALSE); } void CDP_View::OnImageMUL() { // TODO: Add your command handler code here int data=0; CDP_Doc* pDoc=GetDocument(); ASSERT_VALID(pDoc); pDoc->DoubleLoad(); for(int x=0;x<256;x++) { for(int y=0;y<256;y++) { data=pDoc->m_ImageBuffer1[x][y]*pDoc->m_ImageBuffer2[x][y]; pDoc->m_ResultImg[x][y]=data; } } Invalidate(FALSE); } void CDP_View::OnImageDIV() { // TODO: Add your command handler code here int data=0; CDP_Doc* pDoc=GetDocument(); ASSERT_VALID(pDoc); pDoc->DoubleLoad(); for(int x=0;x<256;x++) { for(int y=0;y<256;y++) { data=pDoc->m_ImageBuffer1[x][y]/pDoc->m_ImageBuffer2[x][y]; if(data>255) { pDoc->m_OpenImg[x][y]=255; } else { pDoc->m_OpenImg[x][y]=data; } pDoc->m_ResultImg[x][y]=data; } } Invalidate(FALSE); } void CDP_View::OnBinary() { // TODO: Add your command handler code here CDP_Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); pDoc->Binary(); Invalidate(FALSE); } void CDP_View::OnHistoequal() { // TODO: Add your command handler code here CDP_Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); pDoc->Histo_Equal(); Invalidate(FALSE); } void CDP_View::OnHistostretch() { // TODO: Add your command handler code here CDP_Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); pDoc->Histo_Stretch(); Invalidate(FALSE); } void CDP_View::OnZoomin() { // TODO: Add your command handler code here CDP_Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); pDoc->ZoomIn(); CSize sizeTotal; sizeTotal.cx=300+256*scale_x; sizeTotal.cy=256*scale_y; SetScrollSizes(MM_TEXT, sizeTotal); zoom=1; Invalidate(FALSE); } void CDP_View::OnZoomout() { // TODO: Add your command handler code here CDP_Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); pDoc->ZoomOut(); zoom=2; Invalidate(TRUE); } void CDP_View::OnMirror() { // TODO: Add your command handler code here CDP_Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); pDoc->Mirror(); Invalidate(FALSE); } void CDP_View::OnFilp() { // TODO: Add your command handler code here CDP_Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); pDoc->Flip(); Invalidate(FALSE); } void CDP_View::OnRotate() { // TODO: Add your command handler code here CDP_Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); pDoc->Rotate(); Invalidate(TRUE); } void CDP_View::OnSobel() { // TODO: Add your command handler code here CDP_Doc* pDoc = GetDocument(); ASSERT_VALID(pDoc); pDoc->Sobel(); Invalidate(FALSE); }