векторы на изображении в матлабе

Пожалуйста, взгляните на документацию getpts, она возвращает [X,Y], с Point_i=[X(i), Y(i)].

Поскольку в вашем случае две линии параллельны, они не пересекаются, и поэтому d - это разница x двух выбранных точек. Или это был просто частный случай?

Первый черновик будет выглядеть так:

% Read and display image
i = imread('Sample.jpg');
im = i;
f = figure;
imshow(im);
% Get user clicked points
[X, Y] = getpts(f);
% Draw parallel vertical lines at given pts. parallel to the y-axis/post
hold on;
line([X(1) X(1)], get(gca, 'ylim'));
line([X(2) X(2)], get(gca, 'ylim'));

d=X(2)-X(1);

Для остальной части кода, пожалуйста, будьте более конкретными. Поскольку вы хотите показать вектор нормали изображения, хотите ли вы построить все это в трехмерном графике?

РЕДАКТИРОВАТЬ: что касается комментариев, окончательный код был изменен для отображения 2D-изображения на 3D-графике. Дальнейшее улучшение можно сделать, но это должно быть довольно просто (например, взгляните на это или это для стрелок)

% Read and display image
myImg = imread('Sample.jpg');
f1 = figure;
imshow(myImg);
imSize=size(myImg);
% Get user clicked points
[X, Y] = getpts(f1);
xLim=get(gca, 'xlim');
yLim=get(gca, 'ylim');
close(f1);

% Draw parallel vertical lines at given pts. parallel to the y-axis/post in the image
lWidth = 5; %width of line in Pixel, should be odd
lColor = [255,0,0]; %Linecolor in RGB
pColor = [0,255,0]; %Pointcolor in RGB
vColor = [0,0,255]; %Vectorcolor in RGB

% image positions of points in pixel
xP1=int16((X(1)-xLim(1))/(xLim(2)-xLim(1))*imSize(2));
xP2=int16((X(2)-xLim(1))/(xLim(2)-xLim(1))*imSize(2));
yP1=int16((Y(1)-yLim(1))/(yLim(2)-yLim(1))*imSize(1));
yP2=int16((Y(2)-yLim(1))/(yLim(2)-yLim(1))*imSize(1));

% Draw lines in image
for xLine = [xP1,xP2]
    for i = max([1,xLine-(lWidth-1)/2]):min([imSize(2),xLine+(lWidth-1)/2])
        for j=1:imSize(1)
            myImg(j,i,1:3)=lColor;
        end
    end
end

%mark Points
Points=[[xP1,yP1];[xP2,yP2]];
for k = 1:2
    Pos=Points(k,:);
    for i = max([1,Pos(1)-(lWidth-1)/2]):min([imSize(2),Pos(1)+(lWidth-1)/2])
        for j = max([1,Pos(2)-(lWidth-1)/2]):min([imSize(1),Pos(2)+(lWidth-1)/2])
            myImg(j,i,1:3)=pColor;
        end
    end
end

%3D plot - Some rotations required to get a nice axisvalue alignment
g = hgtransform('Matrix',makehgtform('translate',[0,imSize(1),0])*makehgtform('axisrotate',[0,1,0],pi)*makehgtform('zrotate',pi));
f2 = image(g,myImg);
axis([0,imSize(2),0,imSize(1),-0.1,1.1]) %axis range
view(40,30) %view point

%add normal vectors
hold on
r=(lWidth-1)/2;
[X,Y,Z] = cylinder(r);
xC1=X+double(xP1);
xC2=X+double(xP2);
yC1=imSize(1)-Y-double(yP1);
yC2=imSize(1)-Y-double(yP2);
mesh(xC1,yC1,Z,'facecolor',vColor/255,'edgecolor',vColor/255);
mesh(xC2,yC2,Z,'facecolor',vColor/255,'edgecolor',vColor/255);