2011年09月29日

来源：互联网发布：新疆大学网络管理系统编辑：程序博客网时间：2024/05/20 12:22

转载自（http://www.cnblogs.com/HappyXie/archive/2011/03/02/1969434.html）

//OpenCV 4 下的图像任意角度的旋转

//需要inter公司的OpenCV的支持.

//OpenCV 4下的图像任意角度的旋转
//待旋转的图像IplImage* Img_old
//返回的旋转后图像 IplImage* Img_tmp.
//旋转的角度,单位度.
//三种不同的方法.其中方法二没有完全测试,方法一可以满足大部分需要
IplImage * FitRotate (IplImage* Img_old, double angle,intmethod)

{

IplImage* Img_tmp = NULL;

double anglerad = (CV_PI* (angle/180)) ;
int newheight =int (fabs((sin(anglerad)*Img_old->width )) + fabs((cos(anglerad)*Img_old->height )) );
int newwidth =int (fabs((sin(anglerad)*Img_old->height)) + fabs((cos(anglerad)*Img_old->width)) );

    Img_tmp =cvCreateImage(cvSize(newwidth,newheight), IPL_DEPTH_8U, 3);
cvFillImage(Img_tmp,0);//目的图像使用扩展的大小

IplImage* dst = cvCloneImage( Img_old );//目的图像与原图像等大


    floatm[6];
       CvMat M = cvMat( 2, 3, CV_32F, m );

   if(1==method)
   {
     //方法一提取象素四边形，使用子象素精度

    int w =Img_old->width;
         int h = Img_old->height;

           m[0] = (float)(cos(angle*CV_PI/180.));
           m[1] = (float)(sin(angle*CV_PI/180.));
           m[2] = w*0.5f;
           m[3] = -m[1];
           m[4] = m[0];
           m[5] = h*0.5f;

          cvGetQuadrangleSubPix( Img_old, dst, &M, 1,cvScalarAll(0));

          cvGetQuadrangleSubPix( Img_old, Img_tmp, &M,CV_INTER_LINEAR, cvScalarAll(0));//+CV_WARP_FILL_OUTLIERS

    //方法一提取象素四边形，使用子象素精度
   }

        if(2==method)
   {
   //方法二使用二维旋转的仿射变换矩阵存在问题要求输入和输出图像一样大旋转中心不对

     CvPoint2D32f center;
   center.x=float(Img_old->width/2.0+0.5);//float(Img_tmp->width/2.0+0.5);
   center.y=float(Img_old->height/2.0+0.5);//float(Img_tmp->height/2.0+0.5);
   cv2DRotationMatrix( center,angle,1, &M);

   cvWarpAffine( Img_old, dst,&M,CV_INTER_LINEAR,cvScalarAll(0) );//小图
      //小目标图像

   //对图像进行扩展
  //   只能一定角度以内不同象限的不同对待
          int dx=int((newwidth -Img_old->width )/2+0.5);
        int dy=int((newheight-Img_old->height)/2+0.5);

        uchar* old_ptr,*temp_ptr;

           for( int y=0 ; y<Img_old->height;y++) //为了不越界
   {
          for (int x=0 ; x< Img_old->width;x++)
    {
       old_ptr =&((uchar*)(Img_old->imageData +Img_old->widthStep*y))[(x)*3];
    temp_ptr =&((uchar*)(Img_tmp->imageData +Img_tmp->widthStep*(y+dy)))[(x+dx)*3];
    temp_ptr[0]=old_ptr[0]; //green
    temp_ptr[1]=old_ptr[1]; //blue
    temp_ptr[2]=old_ptr[2]; //Red
    }
   }


   center.x=float(Img_tmp->width/2.0+0.5);
   center.y=float(Img_tmp->height/2.0+0.5);
   cv2DRotationMatrix( center,angle,1, &M);

   IplImage* temp = cvCloneImage(Img_tmp );//生成输出图像
   cvWarpAffine( Img_tmp,temp    ,&M,CV_INTER_LINEAR+CV_WARP_FILL_OUTLIERS,cvScalarAll(0));//大图
     Img_tmp=cvCloneImage( temp );

   //问题
   //cvWarpAffine( Img_tmp,Img_tmp,&M,CV_INTER_LINEAR+CV_WARP_FILL_OUTLIERS,cvScalarAll(0));//大图

   //方法二使用二维旋转的仿射变换矩阵
   }

         if(3==method)
    {
    //方法三透视变换
     CvPoint2D32f src_point[4];
     CvPoint2D32f dst_point[4];

     src_point[0].x=0.0;                   src_point[0].y=0.0;
     src_point[1].x=0.0;                   src_point[1].y=(float) Img_old->height;
     src_point[2].x=(float) Img_old->width;src_point[2].y=(float) Img_old->height;
     src_point[3].x=(float) Img_old->width;src_point[3].y=0.0;

           dst_point[0].x=0;
     dst_point[0].y=(float) fabs((sin(anglerad)*Img_old->width ));

     dst_point[1].x=(float) fabs((sin(anglerad)*Img_old->height));
     dst_point[1].y=(float) fabs((sin(anglerad)*Img_old->width ))+(float) fabs((cos(anglerad)*Img_old->height));

     dst_point[2].x=(float) fabs((sin(anglerad)*Img_old->height))+(float) fabs((cos(anglerad)*Img_old->width));
     dst_point[2].y=(float) fabs((cos(anglerad)*Img_old->height));

     dst_point[3].x=(float) fabs((cos(anglerad)*Img_old->width));
     dst_point[3].y=0;


     floatnewm[9];
           CvMat newM = cvMat( 3, 3, CV_32F, newm );
     cvWarpPerspectiveQMatrix(src_point,dst_point,&newM);

     cvWarpPerspective(Img_old,dst,&newM,CV_INTER_LINEAR+CV_WARP_FILL_OUTLIERS,cvScalarAll(0) );
     cvWarpPerspective(Img_old,Img_tmp,&newM,CV_INTER_LINEAR+CV_WARP_FILL_OUTLIERS,cvScalarAll(0) );

   //方法三透视变换
    }


   // cvNamedWindow( "dst_litter", 1 );
   // cvShowImage( "dst_litter", dst );

// cvNamedWindow( "dst_big", 1 );
// cvShowImage( "dst_big", Img_tmp );

    return Img_tmp;

   }

0 0