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

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 = imread('2013.08.01_13.54_.bmp');
 
X = linspace(min(D), max(D), 10240);
Ai = zeros(size(A,1), size(X,2)); 
d = ones(size(D));
for i = 2:length(D)
    [u d] = min( abs(X-D(i)) ); 
    Ai(d) = A(i); 
end
 
imshow(uint8(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('ampl.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)) );     
    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); 
A = zeros(size(I,1), size(I,2)*10);
for i = 1:size(I,2)
    [u d] = min( abs(X-Dx(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    % j = d(i-1):d(i); 
    A(:,d) = I(:,i); % записываем столбец
end
for i = 1:size(I,1)
    [u d] = min( abs(Y-Dy(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    % j = d(i-1):d(i); 
    B(d,:) = A(i,:); % записываем повторяющиеся строки
end
imshow(uint8(B))

clc; clear all;
format short g
 
I = imread('_google_maps4.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)); 
% расширяем по столбцам
d = 1;
for i = 1:size(I,2)
    [u d] = min( abs(X-Dx(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    
    Ai(:,d) = I(:,i); % записываем столбец
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))
imwrite(uint8(B),'D:\matlab\17.png')

clc; clear all;
format short g
 
I = imread('1527.bmp');
I=rgb2gray(I);
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 = 1:size(I,2)
    [u d] = min( abs(X-Dx(i)) ); 
    % d - номер столбца с единицами, номер элемента в массиве m
    % j = d(i-1):d(i); 
    Ai(:,d) = I(:,i); % записываем столбец
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))

fname = '1527.bmp'; % путь и имя файла
I = imread(fname);
if ndims(I)==3 % если изображение цветное
   I = rgb2gray(I);
end