双线性插值法实现图片的缩放以及旋转

最新推荐文章于 2024-03-05 11:09:48 发布
最新推荐文章于 2024-03-05 11:09:48 发布
阅读量452 收藏 10
点赞数 10
文章标签： 图像处理
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_64175477/article/details/134742324
版权

1.缩放

function result = myScale(f)
            R=f(:,:,1);
            G=f(:,:,2);
            B=f(:,:,3);
            [h,w]=size(R);
            kx=2.3;
            ky=1.6;
            H=ceil(ky*h);
            W=ceil(kx*w);
            for newx=1:W
                for newy=1:H
                    oldx=(newx-1)/kx;
                    oldy=(newy-1)/ky;
                    x=floor(oldx)+1;
                    y=floor(oldy)+1;
                    a=oldx-floor(oldx);
                    b=oldy-floor(oldy);
                    if x<w && y<h
                        r(newy,newx)=R(y,x)+b*(R(y+1,x)-R(y,x))+a*(R(y,x+1)-R(y,x))+a*b*(R(y+1,x+1)+R(y,x)-R(y,x+1)-R(y+1,x));
                    elseif x>=w && y<h
                        r(newy,newx)=R(y,x)+b*(R(y+1,x)-R(y,x));
                    elseif x<w && y>=h
                        r(newy,newx)=R(y,x)+a*(R(y,x+1)-R(y,x));
                    else
                        r(newy,newx)=R(y,x);
                    end
                end
            end

            for newx=1:W
                for newy=1:H
                    oldx=(newx-1)/kx;
                    oldy=(newy-1)/ky;
                    x=floor(oldx)+1;
                    y=floor(oldy)+1;
                    a=oldx-floor(oldx);
                    b=oldy-floor(oldy);
                    if x<w && y<h
                        g(newy,newx)=G(y,x)+b*(G(y+1,x)-G(y,x))+a*(G(y,x+1)-G(y,x))+a*b*(G(y+1,x+1)+G(y,x)-G(y,x+1)-G(y+1,x));
                    elseif x>=w && y<h
                        g(newy,newx)=G(y,x)+b*(G(y+1,x)-G(y,x));
                    elseif x<w && y>=h
                        g(newy,newx)=G(y,x)+a*(G(y,x+1)-G(y,x));
                    else
                        g(newy,newx)=G(y,x);
                    end
                end
            end

            for newx=1:W
                for newy=1:H
                    oldx=(newx-1)/kx;
                    oldy=(newy-1)/ky;
                    x=floor(oldx)+1;
                    y=floor(oldy)+1;
                    a=oldx-floor(oldx);
                    b=oldy-floor(oldy);
                    if x<w && y<h
                        blue(newy,newx)=B(y,x)+b*(B(y+1,x)-B(y,x))+a*(B(y,x+1)-B(y,x))+a*b*(B(y+1,x+1)+B(y,x)-B(y,x+1)-B(y+1,x));
                    elseif x>=w && y<h
                        blue(newy,newx)=B(y,x)+b*(B(y+1,x)-B(y,x));
                    elseif x<w && y>=h
                        blue(newy,newx)=B(y,x)+a*(B(y,x+1)-B(y,x));
                    else
                        blue(newy,newx)=B(y,x);
                    end
                end
            end
            result=cat(3,r,g,blue);
        end

2.旋转

function result = myRotate(f)
            R=f(:,:,1);
            G=f(:,:,2);
            B=f(:,:,3);
            [h,w]=size(R);
            ang=pi/2;
            cornerx=[0 w-1 w-1 0];
            cornery=[0 0 h-1 h-1];
            newcornerx=cornerx*cos(ang)+cornery*sin(ang);
            newcornery=-cornerx*sin(ang)+cornery*cos(ang);
            minx=min(newcornerx);
            miny=min(newcornery);
            H=ceil(max(newcornery)-miny+1);
            W=ceil(max(newcornerx)-minx+1);
            for newx=1:W
                for newy=1:H
                    oldx=(newx-1+minx)*cos(ang)-(newy-1+miny)*sin(ang);
                    oldy=(newx-1+minx)*sin(ang)-(newy-1+miny)*cos(ang);
                    if oldx<-0.5||oldy<-0.5||oldx>w-0.5||oldy>h-0.5
                        red(newy,newx)=0;
                    elseif oldx<0&&oldy<0
                        red(newy,newx)=R(1,1);
                    elseif oldx>w-1&&oldy>h-1
                        red(newy,newx)=R(h,w);
                    elseif oldx<0&&oldy>h-1
                        red(newy,newx)=R(h,1);
                    elseif oldy<0&&oldx>w-1
                        red(newy,newx)=R(1,w);
                    else
                        x=floor(oldx)+1;
                        y=floor(oldy)+1;
                        a=oldx-floor(oldx);
                        b=oldy-floor(oldy);
                        if x<=0&&y>0&&y<h
                            red(newy,newx)=R(y,1)+b*(R(y+1,1)-R(y,1));
                        end
                        if y<=0&&x<w&&x>0
                            red(newy,newx)=R(1,x)+a*(R(1,x+1)-R(1,x));
                        end
                        if x>=w&&y>0&&y<h
                            red(newy,newx)=R(y,w)+b*(R(y+1,w)-R(y,w));
                        end
                        if y>=h&&x<w&&x>0
                            red(newy,newx)=R(h,x)+a*(R(h,x+1)-R(h,x));
                        end
                        if x>0&&y>0&&y<h&&x<w
                            red(newy,newx)=R(y,x)+b*(R(y+1,x)-R(y,x))+a*(R(y,x+1)-R(y,x))+a*b*(R(y+1,x+1)+R(y,x)-R(y+1,x)-R(y,x+1));
                        end
                    end
                end
            end



            for newx=1:W
                for newy=1:H
                    oldx=(newx-1+minx)*cos(ang)-(newy-1+miny)*sin(ang);
                    oldy=(newx-1+minx)*sin(ang)-(newy-1+miny)*cos(ang);
                    if oldx<-0.5||oldy<-0.5||oldx>w-0.5||oldy>h-0.5
                        green(newy,newx)=0;
                    elseif oldx<0&&oldy<0
                        green(newy,newx)=G(1,1);
                    elseif oldx>w-1&&oldy>h-1
                        green(newy,newx)=G(h,w);
                    elseif oldx<0&&oldy>h-1
                        green(newy,newx)=G(h,1);
                    elseif oldy<0&&oldx>w-1
                        green(newy,newx)=G(1,w);
                    else
                        x=floor(oldx)+1;
                        y=floor(oldy)+1;
                        a=oldx-floor(oldx);
                        b=oldy-floor(oldy);
                        if x<=0&&y>0&&y<h
                            green(newy,newx)=G(y,1)+b*(G(y+1,1)-G(y,1));
                        end
                        if y<=0&&x<w&&x>0
                            green(newy,newx)=G(1,x)+a*(G(1,x+1)-G(1,x));
                        end
                        if x>=w&&y>0&&y<h
                            green(newy,newx)=G(y,w)+b*(G(y+1,w)-G(y,w));
                        end
                        if y>=h&&x<w&&x>0
                            green(newy,newx)=G(h,x)+a*(G(h,x+1)-G(h,x));
                        end
                        if x>0&&y>0&&y<h&&x<w
                            green(newy,newx)=G(y,x)+b*(G(y+1,x)-G(y,x))+a*(G(y,x+1)-G(y,x))+a*b*(G(y+1,x+1)+G(y,x)-G(y+1,x)-G(y,x+1));
                        end
                    end
                end
            end


            for newx=1:W
                for newy=1:H
                    oldx=(newx-1+minx)*cos(ang)-(newy-1+miny)*sin(ang);
                    oldy=(newx-1+minx)*sin(ang)-(newy-1+miny)*cos(ang);
                    if oldx<-0.5||oldy<-0.5||oldx>w-0.5||oldy>h-0.5
                        blue(newy,newx)=0;
                    elseif oldx<0&&oldy<0
                        blue(newy,newx)=B(1,1);
                    elseif oldx>w-1&&oldy>h-1
                        blue(newy,newx)=B(h,w);
                    elseif oldx<0&&oldy>h-1
                        blue(newy,newx)=B(h,1);
                    elseif oldy<0&&oldx>w-1
                        blue(newy,newx)=B(1,w);
                    else
                        x=floor(oldx)+1;
                        y=floor(oldy)+1;
                        a=oldx-floor(oldx);
                        b=oldy-floor(oldy);
                        if x<=0&&y>0&&y<h
                            blue(newy,newx)=B(y,1)+b*(B(y+1,1)-B(y,1));
                        end
                        if y<=0&&x<w&&x>0
                            blue(newy,newx)=B(1,x)+a*(B(1,x+1)-B(1,x));
                        end
                        if x>=w&&y>0&&y<h
                            blue(newy,newx)=B(y,w)+b*(B(y+1,w)-B(y,w));
                        end
                        if y>=h&&x<w&&x>0
                            blue(newy,newx)=B(h,x)+a*(B(h,x+1)-B(h,x));
                        end
                        if x>0&&y>0&&y<h&&x<w
                            blue(newy,newx)=B(y,x)+b*(B(y+1,x)-B(y,x))+a*(B(y,x+1)-B(y,x))+a*b*(B(y+1,x+1)+B(y,x)-B(y+1,x)-B(y,x+1));
                        end
                    end
                end
            end
            result=cat(3,red,green,blue);
        end

效果如下：

代码写的有点丑了，写完才意识到可以直接三通道在一个双重循环里完成的

如下：

result_image(newy, newx, 2) = G(y, x) + b*(G(y+1, x)-G(y, x)) + a*(G(y, x+1)-G(y, x)) + a*b*(G(y+1, x+1)+G(y, x)-G(y, x+1)-G(y+1, x));

Matlab用的不熟,写的时候头昏脑胀。

Jzhhhhh!
关注
  • 10
    点赞
  • 踩
  • 10
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
MATLAB中图像利用基本原理进行旋转(双线性插值法)
qq_42132231的博客
03-18 3578
本文采用基本原理对图像进行旋转 使用MATLAB语言， 基本推导过程可参考https://blog.csdn.net/qq_42655135/article/details/88601342 一、程序 clear;close all;clc; gray=rgb2gray(imread('test1.jpg'));%转化为灰度图像 M=size(gray,1); %计算原图像尺寸 行...
插值算法在图像缩放中的评估研究
07-27
图像缩放是数字图像处理的一个基本内容。介绍了图像处理中最常用的3种插值算法:最邻近插值法、双线性插值法、双立方插值法。通过应用不同的插值法对图像进行缩放,比较实验结果,从而验证双立方插值法效果较好。
FPGA图像处理5-基础功能-双线性插值（原理） fpga开发.pdf
04-22
FPGA图像处理5-基础功能-双线性插值（原理） fpga开发.pdf
matlab1_采用双线性插值的方法_一副图像逆时针旋转60°_
09-30
一副图像逆时针旋转60°，采用双线性插值的方法
Matlab学习笔记——双线性插值旋转图像
AndyNJU的博客
09-13 6230
用普通坐标变换的方法进行旋转图像时，由于取整后的坐标不能与原来一一对应，会出现空像素即黑点，例如I=imread('cameraman.tif'); x=pi/4; [m,n]=size(I); J=zeros(m,n); for n1=1:m for n2=1:n u=floor(n1*cos(x)-n2*sin(x))+m; v=floor(n1*sin(
第三周作业：matlab将一张图片进行顺时针旋转 20°，做水平镜像，做错切变换，缩小图像处理，并采用双线性插值方法
一名在读学生
03-19 4272
matlab将一张图片进行顺时针旋转 20°，做水平镜像，做错切变换，缩小图像处理，并采用双线性插值方法
基于双线性插值的图像旋转原理及MATLAB实现（非自带函数）
体会编程语言独到的美
04-29 5133
本文介绍了图像旋转的基本原理及MATLAB实现，在不借助MATLAB自带函数的情况下，自己书写了实现图像旋转步骤的几个函数，使用的插值方法为双线性插值。
图像缩放旋转 双线性和临域法
11-21
图像缩放旋转 双线性和临域法 本文主要讲解使用 MATLAB 对图像进行缩放和旋转的两种方法：双线性插值法和最近邻域法。 一、图像缩放 图像缩放是指将图像的尺寸缩放到指定的大小，缩放的方法有很多种，常用的有...
数字图像处理：最邻近、双线性、双三次插值的MATLAB实现
03-21
在MATLAB中，有几种不同的插值方法用于在缩放过程中填充新像素，包括最邻近插值、双线性插值和双三次插值。这三种方法各有优缺点，适用于不同的场景。 首先，我们来看最邻近插值（Nearest Neighbour Interpolation...
BMP旋转放缩(双线性插值，三次卷积插值)的c语言实现(含主程序可直接运行)
05-08
项目提供了C语言的实现，包括双线性插值和三次卷积插值两种方法，这两种都是提高图像缩放质量和保真度的重要技术。 双线性插值是一种在二维空间中进行插值的方法，它通过相邻四个像素的灰度值来估计新位置的像素值...
双线性插值进行图像缩放
08-14
双线性插值进行图像缩放
利用双线性插值实现的简单图像缩放
05-01
实现一个图像缩放函数，可以对输入图像进行任意倍数的缩放； 采用双线性插值进行重采样； X,Y方向的缩放倍数参函数参数的形式传入； 可以只考虑输入图像为3通道，8位深度的情况； 不能调用图像处理库的缩放函数来完成；
双线性插值图像放大缩小算法
10-17
具体讲解了双线性插值图像放大缩小算法，希望对大家有帮助
双线性插值算法实现图像旋转
05-17
用双线性插值算法估计出象素值，对于超出图像范围的象素返回255
MATLAB实现数图缩放：双线性内插法
陈小星的博客
09-16 3619
MATLAB实现数图缩放：双线性内插法
使用Matlab进行图像的读写、显示和缩放（最近临插值和双线性内插值法）
热门推荐
梦空间
10-18 2万+
上次我们开始进行数字图像处理这门课程的实验，直到现在才抽空出来写写文章，记录一下知识点。介绍一下，使用Matlab对数字图像的简单处理。 1、 读取与显示输入图像： %输入图像和显示图像 function []=readShow() ima=imread('test.jpg'); %输入图像 imshow(ima); %显示图像 end 显示输出图像 --> readShow() 输出
双线性插值缩放算法原理以及matlab与verilog的实现（一）
主要分享学习中的知识和遇到的问题
03-05 2632
视频图像缩放技术在数字图像处理领域中有着广泛的应用。现在各种液晶设备的分辨率不同，视频图像输入的分辨率也各不相同，想要在显示器上正确的显示出相应图像画面，就必须对输入的图像大小进行缩放调整到显示屏支持的分辨率。图像缩放算法有很多种，常见的包括：主流使用的缩放算法通常是双线性插值和双三次插值。双线性插值在实现简单的同时质量较高，适合一般性的图像缩放需求；双三次插值在质量上更为优秀，能够更好地保持图像细节，但计算量较大，适合对图像质量要求较高的场景。
图像缩放--双线性内插法及其python实现（图文代码）
Sail的博客
11-24 3886
双线性内插法 #数字图像处理 理论简介 双线性插值是图像内插缩放的一种方法。 简单来说，双线性插值即对于目标像素进行两个方向上的线性插值。 比如，首先在x上进行一次线性插值得到两个点，即上图得到的 R1R_1R1​，R2R_2R2​。再在y方向上进行一次线性插值，由 R1R_1R1​，R2R_2R2​ 得到目标点PPP。 故双线性插值需要四个最近邻像素点（原像素点，即已知灰度值），即上图中的 Q11,Q12,Q21,Q22Q_{11},Q_{12},Q_{21},Q_{22}Q11​,Q12​,Q21​,
07. 基于MFC实现双线性插值图像缩放算法
qq_27719791的博客
05-15 516
本博文内容是博文的一部分(博文实现了最简单的最近邻算法，本文在此基础上添加双线性插值缩放算法的实现。*双线性插值算法(Bilinear Interpolation)*原理如下图所示。所谓，即在两个方向分别进行一次线性插值，通过四个相邻像素插值得到待求像素，距离越近的像素贡献越大。：已知c00，c10，c01，c11为原图中的四邻像素，c点为待求像素。：双线性插值分两步实现。
matlab双线性插值法旋转图像
最新发布
06-27
在MATLAB中，双线性插值法（Bilinear Interpolation）是一种常用的图像旋转和缩放的方法，特别是在没有旋转矩阵或图像几何变换矩阵的情况下。旋转图像时，原图像的每个像素位置需要根据旋转后的坐标进行插值，以保持...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值