%%% Inverse_Spherize
%%% 逆球面化 clc;
clear all;
close all; addpath('E:\PhotoShop Algortihm\Image Processing\PS Algorithm'); I=imread('4.jpg');
Image=double(I); [row, col,layer]=size(Image);
R=floor(min(row, col)/2);
K=pi/2;
Image_new=Image;
alpha=1;
for i=1:row
for j=1:col
y0=(row+1)/2-i;
x0=j-(col+1)/2;
Dis=x0.^2+y0.^2;
radius0=sqrt(Dis);
if(Dis<R*R)
if(x0~=0)
beta=atan(y0/x0);
if(x0<0)
beta=beta+pi;
end
else
beta=pi/2;
end
theta=radius0/R*K;
radius=R*sin(theta);
radius=(radius-radius0)*alpha+radius0;
x=radius*cos(beta);
y=radius*sin(beta);
x=col/2+x;
y=row/2-y;
if(x>1 && x<col && y>1 && y<row)
x1=floor(x);
y1=floor(y);
p=x-x1;
q=y-y1;
for k=1:layer
Image_new(i,j,k)=(1-p)*(1-q)*Image(y1,x1,k)+p*(1-q)*Image(y1,x1+1,k)...
+q*(1-p)*Image(y1+1,x1,k)+p*q*Image(y1+1,x1+1,k);
end
end
end
end
end imshow(Image_new/255);

原图

效果图

最新文章

  1. Sharing A Powerful Tool For Calculate Code Lines
  2. 有关docker的学习链接
  3. HDOJ(2056)&amp;HDOJ(1086)
  4. SpringMVC与Struts2区别与比较总结
  5. Atitit.如何选择技术职业方向
  6. Pythonxy下载链接
  7. POJ 2559 Largest Rectangle in a Histogram -- 动态规划
  8. Oracle OEM重建
  9. phpexcel的写出操作(生成excel表)
  10. HP的笔记本经常蓝屏崩溃 -------athr.sys
  11. 记一次使用修改字节码的方法解决java.lang.NoSuchMethodError
  12. github+hexo搭建自己的博客网站(五)进阶配置(畅言实现博客的评论)
  13. Unity3D各平台Application.xxxPath的路径
  14. UNIX网络编程——客户/服务器程序设计示范(四)
  15. pymongo 操作
  16. uva 202
  17. linux下expect环境安装以及简单脚本测试
  18. sourceInsight工具移除不掉项目 source Insight Add and Remove Project Files
  19. B/S开发介绍
  20. 分享:selenium(一) xpath

热门文章

  1. 【hiho一下】第一周 最长回文子串
  2. java.lang.IllegalStateException: ImageView no longer exists. You should not use this PhotoViewAttacher any more.
  3. eclipse 给jar包关联javadoc
  4. (比赛)B - Super Mobile Charger
  5. EasyNVR无插件摄像机直播之:摄像机网页低延时无插件直播实现
  6. python+NLTK 自然语言学习处理八:分类文本一
  7. maven编译问题-maven项目运行时找不到文件,解决方案之一
  8. Java并发之BlockingQueue
  9. 【Mybatis】Insert批量操作
  10. deviceToken的获取(二)