Как преобразовать декартово изображение в Полярное изображение

Question

Я использую приведенный ниже код для преобразования изображения в Полярный к Декартовому.

---

for ii=1:500; 
    for jj=1:500 
     phase(ii,jj)=((ii-250)^2/1000+(jj-250)^2/1000); 
    end 
end 
Cartesian=cos(phase); 
[imrow,imcol]=size(Cartesian); 

% choose the center of the image 
rcent=250; 
ccent=250; 

rmax=sqrt((imrow-rcent)^2+(imcol-ccent)^2); 
%prepare the gridspace in the transformed coordinate 
[r,theta]=meshgrid(linspace(0,rmax,imrow),linspace(0,2*pi,imcol)); 
%corresponding locations in the original coordinate 
x=r.*cos(theta)+ccent; 
y=r.*sin(theta)+rcent; 

%sample original fringe pattern using interpolation 
Polar=interp2(Cartesian,x,y); 

subplot(1,2,1); imagesc(Cartesian) ; axis square 
subplot(1,2,2); imagesc(Polar) ; axis square 

---

Теперь я хотел бы преобразовать декартовы в Заполярном. Я бы очень признателен за любую помощь.

С уважением,

J. Cooper

Answer 1

Спасибо за ответы. С уважением

Я чрезвычайно новичок в MATLAB. Я много пытался в Google. Однако я не нашел ответа. Они противоположны моей цели. Вышеприведенный код относится к одной из тем, которые я имел в Google. На самом деле, я пробовал это:

r1 = sqrt(x.^2+y.^2); 
theta1 = atan(y./x); 
Cartesian2 = interp2(Polar,r1,theta1); 

subplot(1,3,1); imagesc(Cartesian) ; axis square 
subplot(1,3,2); imagesc(Polar) ; axis square 
subplot(1,3,3); imagesc(Cartesian2) ; axis square 

Но это не работает.