qt图像算法—图像的旋转之c++实现(不调包)

1.基本原理

本章所讲的图像旋转算法就是让图像按照其中心点旋转指定角度。图像在旋转之后,其内容不会变形,只是角度改变了,但图像的大小会变化(图像一般会变大,能包含住原先的内容)。图像的旋转变化,比较复杂,涉及到多次矩阵变换。

1.旋转公式

在下图中将坐标为(x0, y0)的像素点顺时针旋转角度后,其坐标变成(x, y),其中r表示像素坐标距离原点的距离,表示旋转前像素点与原点连线的度数。

可推导出公式:

2.坐标变换

图像坐标系与数学坐标系是不相同,所以在旋转过程中需要进行两次坐标变换。

1)旋转操作前

将图像坐标系转到数学坐标系,公式中,W和H为原始图像的宽高。

2)旋转操作后

将旋转后图像的数学坐标系转到图像坐标系(相当于转回来),公式中,Wnew和Hnew分别表示输出图像的宽高。

3.旋转公式

有了前面的铺垫,我们知道图像的旋转要经过如下三步才能完成:

1.输入图像的坐标系转数学坐标系

2.通过旋转关系计算出旋转后的坐标

3.旋转坐标系(数学坐标系)转图像坐标系

总的公式可见下面(输出坐标映射到原坐标):

2.代码实现(代码是我以前自学图像处理时写的,代码很粗糙没做任何优化,但很好理解)

cpp 复制代码
/*图像的旋转函数(最临近插值法) angle为旋转度数,以弧度表示*/
QImage* MainWindow::RotateNormal(QImage* image,double angle)
{
    int srcX1, srcX2, srcX3, srcX4;
    int srcY1, srcY2, srcY3, srcY4;

    srcX1 = 0;
    srcY1 = 0;
    srcX2 = image->width() - 1;
    srcY2 = 0;
    srcX3 = 0;
    srcY3 = image->height() - 1;
    srcX4 = image->width() - 1;
    srcY4 = image->height() - 1;

    double fSin = sin(angle);
    double fCos = cos(angle);

    double tranX1, tranX2, tranX3, tranX4;
    double tranY1, tranY2, tranY3, tranY4;

    tranX1 = fCos * srcX1 + fSin * srcY1;
    tranY1 = -fSin * srcX1 + fCos * srcY1;
    tranX2 = fCos * srcX2 + fSin * srcY2;
    tranY2 = -fSin * srcX2 + fCos * srcY2;
    tranX3 = fCos * srcX3 + fSin * srcY3;
    tranY3 = -fSin * srcX3 + fCos * srcY3;
    tranX4 = fCos * srcX4 + fSin * srcY4;
    tranY4 = -fSin * srcX4 + fCos * srcY4;

    unsigned int outWidth = (unsigned int)(max(fabs(tranX4 - tranX1), fabs(tranX3 - tranX2)) + 1.5);
    unsigned int outHeight = (unsigned int)(max(fabs(tranY4 - tranY1), fabs(tranY3 - tranY2)) + 1.5);

    QImage* newImage = new QImage(outWidth,outHeight,QImage::Format_ARGB32);

    double num1 = -0.5 * outWidth * fCos - 0.5 * outHeight * fSin + 0.5 * image->width();
    double num2 =  0.5 * outWidth * fSin - 0.5 * outHeight * fCos + 0.5 * image->height();

    unsigned char* copyPixel = NULL;
    unsigned char* objPixel = NULL;
    int x = 0;
    int y = 0;

    for(long j = 0; j < (long)outHeight; j++)
    {
        for(long i = 0; i <(long)outWidth; i++)
        {
            x = (int)(i * fCos + j * fSin + num1 + 0.5);
            y = (int)(-i * fSin + j * fCos + num2 + 0.5);

            if(x == image->width())
            {
                x--;
            }
            if(y == image->height())
            {
                y--;
            }
            copyPixel = image->bits() + y * image->width() * 4  + x * 4;
            objPixel = newImage->bits() + j * outWidth * 4 + i * 4;

            if(x >= 0 && x < image->width() && y >=0 && y < image->height())
            {
                memcpy(objPixel, copyPixel, 4);
            }
        }
    }
    return newImage;
}

/*图像的旋转函数(双线性插值法) angle为旋转度数,以弧度表示*/
QImage* MainWindow::RotateInterpolation(QImage* image,double angle)
{
    int srcX1, srcX2, srcX3, srcX4;
    int srcY1, srcY2, srcY3, srcY4;

    srcX1 = 0;
    srcY1 = 0;
    srcX2 = image->width() - 1;
    srcY2 = 0;
    srcX3 = 0;
    srcY3 = image->height() - 1;
    srcX4 = image->width() - 1;
    srcY4 = image->height() - 1;

    double fSin = sin(angle);
    double fCos = cos(angle);

    double tranX1, tranX2, tranX3, tranX4;
    double tranY1, tranY2, tranY3, tranY4;

    tranX1 = fCos * srcX1 + fSin * srcY1;
    tranY1 = -fSin * srcX1 + fCos * srcY1;
    tranX2 = fCos * srcX2 + fSin * srcY2;
    tranY2 = -fSin * srcX2 + fCos * srcY2;
    tranX3 = fCos * srcX3 + fSin * srcY3;
    tranY3 = -fSin * srcX3 + fCos * srcY3;
    tranX4 = fCos * srcX4 + fSin * srcY4;
    tranY4 = -fSin * srcX4 + fCos * srcY4;

    long outWidth = (unsigned int)(max(fabs(tranX4 - tranX1), fabs(tranX3 - tranX2)) + 1.5);
    long outHeight = (unsigned int)(max(fabs(tranY4 - tranY1), fabs(tranY3 - tranY2)) + 1.5);

    QImage* newImage = new QImage(outWidth,outHeight,QImage::Format_ARGB32);

    double num1 = -0.5 * outWidth * fCos - 0.5 * outHeight * fSin + 0.5 * image->width();
    double num2 =  0.5 * outWidth * fSin - 0.5 * outHeight * fCos + 0.5 * image->height();

    double x = 0.0;
    double  y = 0.0;

    int r,g,b;
    for (long  j = 0; j < outHeight; j++)
    {

        for(long i =0; i < outWidth; i++)
        {
            x = (i * fCos + j * fSin + num1 + 0.5);
            y = (-i * fSin + j * fCos + num2 + 0.5);


            if (x > image->width() || x < 0 || y > image->height() || y < 0)
                continue;
            int x1, x2, y1, y2;
            x1= ( int)x;
            x2 = x1 + 1;
            y1 = ( int)y;
            y2 = y1 + 1;

            QColor oldcolor1;
            QColor oldcolor2;
            QColor oldcolor3;
            QColor oldcolor4;
            double u, v;
            u = x - x1;
            v = y - y1;
            if ((x >= image->width() - 1 ) && (y >= image->height() - 1 ))
            {
            oldcolor1 = QColor(image->pixel(x1,y1));
            r = oldcolor1.red();
            g = oldcolor1.green();
            b = oldcolor1.blue();
            }
            else if (x >= image->width() - 1)
            {
                oldcolor1 = QColor(image->pixel(x1,y1));
                oldcolor3 = QColor(image->pixel(x1,y2));
                r = oldcolor1.red() * (1 - v) + oldcolor3.red() * v;
                g = oldcolor1.green() * (1 - v) + oldcolor3.green() * v;
                b = oldcolor1.blue() * (1 - v) + oldcolor3.blue() * v;
            }
            else if (x > image->height() - 1)
            {
                oldcolor1 = QColor(image->pixel(x1,y1));
                oldcolor2 = QColor(image->pixel(x2,y1));
                r = oldcolor1.red() * (1 - u) + oldcolor2.red() * u;
                g = oldcolor1.green() * (1 - u) + oldcolor2.green() * u;
                b = oldcolor1.blue() * (1 - u) + oldcolor2.blue() * u;
            }
            else
            {
                oldcolor1 = QColor(image->pixel(x1,y1));
                oldcolor2 = QColor(image->pixel(x2,y1));
                oldcolor3 = QColor(image->pixel(x1,y2));
                oldcolor4 = QColor(image->pixel(x2,y2));
                int r1,g1,b1;
                r = oldcolor1.red() * (1 - u) + oldcolor2.red() * u;
                g = oldcolor1.green() * (1 - u) + oldcolor2.green() * u;
                b = oldcolor1.blue() * (1 - u) + oldcolor2.blue() * u;

                r1 = oldcolor3.red() * (1 - u) + oldcolor4.red() * u;
                g1 = oldcolor3.green() * (1 - u) + oldcolor4.green() * u;
                b1 = oldcolor3.blue() * (1 - u) + oldcolor4.blue() * u;

                r = r * (1 - v) + r1 * v;
                g = g * (1 - v) + g1 * v;
                b = b * (1 - v) + b1 * v;
            }

          newImage->setPixel(i, j, qRgb(r, g, b));
        }
    }
    return newImage;
}

3.项目源码下载:

整套算法系列:
https://blog.csdn.net/u013289254/category_12811658.html?spm=1001.2014.3001.5482https://blog.csdn.net/u013289254/category_12811658.html?spm=1001.2014.3001.5482项目源码下载地址:关注WX【AI街潜水的八角】,回复【qt图像算法】即可下载

整套项目源码内容包含

1.根据算法原理,编写纯c++源码,不调用外源库opencv 等;

2.包括各种图像处理的基本算法,包含腐蚀膨胀,缩放,转置,镜像,平移,均衡变化,灰度拉升,灰度阈值,灰度非线性,转灰度,灰度线性,旋转,简单平滑,高斯平滑,轮廓跟踪,种子算法,hough直线检测,拉普拉斯,带方向边缘检测,常规边缘检测(梯度算子、Roberts算子和Sobel算子),中值滤波,反色操作等;

3.程序中有完整的注释,便于大家很好理解代码。

相关推荐
zoujiahui_20185 分钟前
Python 包与环境管理工具 uv
开发语言·python·uv
小小龙学IT5 分钟前
Qt 模型/视图框架(Model/View)深度解析:从 MVC 演化到数据展示架构的王者
qt·架构·mvc
吴声子夜歌18 分钟前
Java面试题——JVM(二)
java·开发语言·jvm
CRMEB20 分钟前
商城大促活动策划全流程:从定目标到复盘四阶段
java·开发语言·人工智能·ai·开源·php
denglinqings21 分钟前
Java开发从入门到精通所有课程
java·开发语言
程序猿阿森1 小时前
Python 闭包与装饰器:从入门到精通
开发语言·python·面试
余额瞒着我当琳1 小时前
C++--vector底层,leetcode118杨辉三角实战,对于传统reserve的坑点,迭代器失效
开发语言·c++
曹牧1 小时前
C#:字符串做20位截断
开发语言·c#
Quz1 小时前
QML ComboBox 完全自定义样式:背景、内容、下拉项与箭头
qt
醉城夜风~1 小时前
(超详细)C++多态(Polymorphism)
开发语言·c++