matlab练习程序（图像球面化）

心已赠人 2021-12-16 10:39 280阅读 0赞

十一没什么事干，接着看图像算法。

这个球面化算法最初是在ps上的球面化滤镜中看到的，感觉挺有意思，就研究了一下。

算法的详细推导可以在[这篇博客][Link 1]中找到，我比较懒，只在纸上推了一遍，就不在博客上编辑了。

不过这里还是要把逆变换公式写一下。

公式如下：

![eq1][]

![eq2][]

其中R为球的半径，x,y为目标图像像素坐标，xx,yy为源图像像素坐标。

原图：

![011935171448086.jpg][]

球面化后：

![011936167374359.jpg][]

matlab代码如下：

clear all;close all;clc;
    
    img=imread('lena.jpg');
    [h w]=size(img);
    imshow(img);
    
    imgn=zeros(h,w);
    R=h/2;
    cenX=w/2;
    cenY=h/2;
    theta=pi;
    for y=1-cenY:h-cenY
        for x=1-cenX:w-cenX
            
            disX=1.3*x;     %系数为1则半径为h/2
            disY=1.3*y;
          
            dis=disX^2+disY^2;
            r=sqrt(disX^2+disY^2);
            if r<=R        
                xx=2*R*disX*acos(sqrt(R^2-dis)/R)/(theta*r)+cenX;
                yy=2*R*disY*acos(sqrt(R^2-dis)/R)/(theta*r)+cenY;
                          
                xx=round(xx);
                yy=round(yy);
                if xx>=1 && xx<=w && yy>=1 && yy<=h
                    imgn(y+cenY,x+cenX)=img(yy,xx);
                end
            else
                imgn(y+cenY,x+cenX)=img(y+cenY,x+cenX);
            end
            
        end
    end
    
    figure;
    imshow(imgn,[])

转载于:https://www.cnblogs.com/tiandsp/p/4003524.html

[Link 1]: http://blog.csdn.net/xeral/article/details/6593416
[eq1]: http://latex.codecogs.com/gif.latex?xx%3D2*R*x*%5Carccos%20%28%5Csqrt%7BR%5E2-x%5E2-y%5E2%7D/R%7D%29/%28%5Cpi%20*%5Csqrt%7Bx%5E2+y%5E2%7D%29
[eq2]: http://latex.codecogs.com/gif.latex?yy%3D2*R*y*%5Carccos%20%28%5Csqrt%7BR%5E2-x%5E2-y%5E2%7D/R%7D%29/%28%5Cpi%20*%5Csqrt%7Bx%5E2+y%5E2%7D%29
[011935171448086.jpg]: /images/20211214/3af2f3165a6a42f69f7601b0f7904ff8.png
[011936167374359.jpg]: /images/20211214/6f47339f944b4246b1066a7c0a34b13c.png