Расширение изображения с преобразованием наклонных дальностей в земные координаты

clear, clc
A = dlmread('А.txt'); % амплитуды
R = 750 + 1.4*(0:1023); % расстояния
h = 600; % высота
D = sqrt( R.^2 - h^2 ); % исходные координаты
Di = linspase(min(D), max(D), 10240 ); % расширенный массив 
Ai = spline(D, A, Di); % интерполируем
plot(Di, Ai)

clc; clear all;
format short g
n = 0:1024-1;
R = 933 + 1.4*n;
n1 = 0:10240-1;
c=cos(600./(933+n1*0.14));
D = 714.5+n1.*(0.14./c);
R(n)=A(n);
A=load('A.txt');

clc; clear all;
format short g
n = 0:1024-1;
R = 933 + 1.4*n;
c = cos(600./(933+n*0.14));
D = 714.5 + n.*(0.14./c);
A = dlmread('А.txt');
 
X = linspace(min(D), max(D), 10240);
Ai = zeros(size(X));
for i = 1:length(D)
    [u d] = min( abs(X-D(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    Ai(d) = A(i); % записываем столбец
end
plot(X,Ai)

clc; clear all;
format short g
n = 0:1024-1;
R = 933 + 1.4*n;
c = cos(600./(933+n*0.14));
D = 714.5 + n.*(0.14./c);
A = dlmread('А.txt');
 
X = linspace(min(D), max(D), 10240); % новый расширенный массив
 
Ai = zeros(size(X)); 
d = ones(size(D));
for i = 2:length(D)
    [u d(i)] = min( abs(X-D(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    Ai( d(i-1):d(i) ) = A(i); % записываем столбец
end
plot(X,Ai)

clc; clear all;
format short g
n = 0:1024-1;
R = 933 + 1.4*n;
c = cos(600./(933+n*0.14));
D = 714.5 + n.*(0.14./c);
A = dlmread('А.txt');
Img = imread('2013.08.01_13.54_.bmp');
 
X = linspace(min(D), max(D), 10240); % новый расширенный массив
 
Ai = zeros(size(Img,1), size(X,2)); 
d = ones(size(D));
for i = 2:length(D)
    [u d(i)] = min( abs(X-D(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    j = d(i-1):d(i); 
    Ai(:,j) = Img(:,i*ones(size(j))); % записываем столбец
end
imshow(uint8(Ai))

clc; clear all;
format short g
 
I = imread('2013.08.01_13.54_.bmp');
 
nx = 0:size(I,2)-1;
ny = 0:size(I,1)-1;
Rx = 933 + 1.4*nx;
Ry = 933 + 1.4*ny;
cx = cos(600./(933+nx*0.14));
cy = cos(600./(933+ny*0.14));
Dx = 714.5 + nx.*(0.14./cx);
Dy = 714.5 + ny.*(0.14./cy);
 
X = linspace(min(Dx), max(Dx), 10*size(I,2)); % новый расширенный массив
Y = linspace(min(Dy), max(Dy), 10*size(I,1)); % новый расширенный массив
 
B = zeros(size(I)*10); 
% расширяем по столбцам
d = 1;
for i = 2:size(I,2)
    [u d(i)] = min( abs(X-Dx(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    j = d(i-1):d(i); 
    Ai(:,j) = I(:,i*ones(size(j))); % записываем столбец
end
 
% расширяем по строкам
d = 1;
for i = 2:size(I,1)
    [u d(i)] = min( abs(Y-Dy(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    j = d(i-1):d(i); 
    B(j,:) = repmat(Ai(i,:),[length(j),1]); % записываем строку
end
imshow(uint8(B))

clc; clear all;
format short g
 
I = imread('2013.08.01_13.54_.bmp');
 
nx = 0:size(I,2)-1;
ny = 0:size(I,1)-1;
Rx = 933 + 1.4*nx;
Ry = 933 + 1.4*ny;
cx = cos(600./(933+nx*0.14));
cy = cos(600./(933+ny*0.14));
Dx = 714.5 + nx.*(0.14./cx);
Dy = 714.5 + ny.*(0.14./cy);
 
X = linspace(min(Dx), max(Dx), 10*size(I,2)); % новый расширенный массив
Y = linspace(min(Dy), max(Dy), 10*size(I,1)); % новый расширенный массив
 
B = zeros(size(I)*10); 
d = 1;
A = zeros(size(I,1), size(I,2)*10);
for i = 2:size(I,2)
    [u d(i)] = min( abs(X-Dx(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    j = d(i-1):d(i); 
    A(:,j) = I(:,(i-1)*ones(size(j))); % записываем столбец
end
d = 1;
for i = 2:size(I,1)
    [u d(i)] = min( abs(Y-Dy(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    j = d(i-1):d(i); 
    B(j,:) = repmat(A(i-1,:),[length(j),1]); % записываем повторяющиеся строки
end
imshow(uint8(B))