Matlab 图像转极坐标系
实心圆环
imgVP1=flip(imgVP1,1); % 水平翻转
polarVP1=polarVolinPlot(imgVP1); % 调用函数
空心圆环
[m,n,~]=size(imgVP2);
% 填充空心部分
for i=m:m+300
for j=1:n
imgVP2(i,j,:)=[255,255,255];
end
end
imgVP2=flip(imgVP2,1);
polarVP2=polarVolinPlot(imgVP2);
function [polarImg] = polarImg(img0)
% polarImg
Rimg0=double(img0(:,:,1));
Gimg0=double(img0(:,:,2));
Bimg0=double(img0(:,:,3));
[m,n]=size(Rimg0);
%原图坐标
[t,r]=meshgrid(linspace(-pi,pi,n),1:m);
M=2*m;N=2*n;
[NA,MA]=meshgrid((1:N)-n-0.5,(1:M)-m-0.5);
T=atan2(NA,MA);
R=sqrt(MA.^2+NA.^2);%变换坐标
Rpolar=interp2(t,r,Rimg0,T,R,'linear',0);
Gpolar=interp2(t,r,Gimg0,T,R,'linear',0);
Bpolar=interp2(t,r,Bimg0,T,R,'linear',0);
polarImg=uint8(cat(3,Rpolar,Gpolar,Bpolar));
polarImg=imrotate(polarImg, -90);
% 图片裁剪
if M<N
polarImg = imcrop(polarImg,[0,(N-M)/2,M,M]);
else
polarImg = imcrop(polarImg,[0,(M-N)/2,N,N]);
end
end
---------------------
最新文章
- Android开发学习之路-Git的极简教程？
- Inventory Costing in AX 2009
- tomcat，zookeeper，activeMQ，Kafka设置jvm参数
- SSH Tunneling
- hdu4578 Transformation
- Unsupported configuration attributes: [FILE_UPLOAD]
- 转：已知2个整形数据a,b.不使用if,?:以及其他任何条件判断的语法，找出a跟b中数据的大者。
- Xmanager连接CentOS的远程桌面
- JAVA简单的SWING及AWT
- Lotto
- Android studio一些常见技巧(不断更新)
- 抛弃JQ，回归原生js……
- 第一次使用github、git工具，本地仓库、远程仓库使用
- LightGBM算法（转载）
- spring的历史和设计科学
- spawn
- PS切图采坑
- python拉格朗日插值
- Latex使用的注意事项
- linux Ctrl+z和Ctrl+c的区别