How to draw orthogonal lines ?

16 次查看（过去 30 天）
SUSHMA MB
SUSHMA MB 2015-8-14
评论： Image Analyst 2015-8-15
I have two points, one is start and the other is end point. I have connected these two points by a straight line. Now i want to draw orthogonal lines over the line. How can i draw it? I have attached a picture for reference. As shown in the picture, i want eight orthogonal planes placed in equidistant. Please help me with this problem.

请先登录，再回答此问题。

回答（3 个）

blaat
blaat 2015-8-14
The line between your two points can be described by
y = a (x - x1) + b,
where
a = (y2 - y1)/(x2 - x1)
b = y1,
if we call your two points (x1, y1) and (x2, y2). Lines perpendicular to the original line will have a slope of -1/a and can be expressed as:
y = -1/a (x - x0) + y0,
where (x0, y0) is the point on the original line where it intersects the orthogonal line.
Equidistant points on the line can be easily computed using linspace():
num_orth = 8;
x_orth = linspace(x1, x2, num_orth);
y_orth = linspace(y1, y2, num_orth);
Mike Garrity
Mike Garrity 2015-8-14
编辑：Mike Garrity 2015-8-14
Let's say we have the following line:
pt1 = 10*randn(1,2);
pt2 = 10*randn(1,2);
line([pt1(1), pt2(1)],[pt1(2),pt2(2)])
The points where four equally spaced orthogonal lines cross it are:
n = 4;
t = linspace(0,1,n+2); % evenly spaced parameters
t = t(2:(end-1)); % we don't need the start and end points
v = pt2 - pt1;
x = pt1(1) + t*v(1); % p(t) = p1 + t*(p2-p1)
y = pt1(2) + t*v(2);
h = line(x,y);
h.LineStyle = 'none';
h.Marker = 'o';
Next we need to normalize that vector:
delete(h)
v = v / norm(v);
And then we rotate it by 90 degrees. That's just swapping the X & Y components of v, and changing the sign of one:
for i=1:n
line([x(i)+v(2), x(i)-v(2)],[y(i)-v(1), y(i)+v(1)]);
end
The one catch at this point is that the axes might be using different scale factors for the X & Y. This will make the lines look like they're not orthogonal, even if they are mathematically. We can fix this by calling:
axis equal
  2 个评论
Mike Garrity
Mike Garrity 2015-8-14
Things are a bit more interesting in 3D. I discussed that case on the MATLAB Graphics blog in these two posts ( post1 , post2 ).
SUSHMA MB
SUSHMA MB 2015-8-14
编辑：SUSHMA MB 2015-8-14
Hi...Thank you for the answer. It works. But my points are positive, but why am i getting negative axis. Can i adjust the axis, like i want only positive axis. Like i have points: start_point = [5.95 37.55]; goal_point = [35.62 5.73]; and i want to define the length of the orthogonal lines, i.e., 5 units towards both the side from the point of intersection. Is it possible?

请先登录，再进行评论。

Mike Garrity
Mike Garrity 2015-8-14
I'm afraid I don't understand the question "why am i getting negative axis".
The length of the orthogonal lines goes at the point where I divided v by norm(v). That was to get them to be unit length. Just multiply by the length you want:
start_point = [5.95 37.55]
goal_point = [35.62 5.73]
line([start_point(1), goal_point(1)],[start_point(2), goal_point(2)],'Marker','o')
n = 8;
t = linspace(0,1,n+2);
t = t(2:(end-1));
v = goal_point - start_point;
x = start_point(1) + t*v(1);
y = start_point(2) + t*v(2);
v = 5*v / norm(v);
for i=1:n
line([x(i)+v(2), x(i)-v(2)],[y(i)-v(1), y(i)+v(1)]);
end
axis equal
  3 个评论
显示 1更早的评论
SUSHMA MB
SUSHMA MB 2015-8-15
And how to further divide each of theses orthogonal lines in an equal spacing. I mean i want to further divide the line into an equal spacing of 8, but along the orthogonal lines. And this time i only want the points, not necessary to visualize the points on the orthogonal lines. I hope i made myself clear. Thank you in advance.
Image Analyst
Image Analyst 2015-8-15
It's just simple 10th grade algebra. But if you want the MATLAB way, just use linspace():
xEquallySpaced = linspace(x(1), x(2), numPoints);
yEquallySpaced = linspace(y(1), y(2), numPoints);

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by