补充数学知识:
直线方程一般式 Ax + By + C = 0
点斜式 y-y1 = k(x - x1)
斜切式 y = kx +b
两点式 (y-y1)/(y2-y1) =(x-x1)/(x2-x1) x1不等于x2 y1不等于要y2
截距式 x/a + y/b = 1 a,b都不等于0
斜率,亦称“角系数”,表示一条直线相对于横坐标轴的倾斜程度。一条直线与某平面直角坐标系横坐标轴正半轴方向的夹角的正切值即该直线相对于该坐标系的斜率。
表达式
k=tanα=(y2-y1)/(x2-x1)
对于一次函数y=kx+b(斜截式)k即该函数图像的斜率。
相关公式
当直线L的斜率存在时,斜截式y=kx+b 当k=0时 y=b
当直线L的斜率存在时,点斜式y2—y1=k(X2—X1),
当直线L在两坐标轴上存在非零截距时,有截距式ax + by +c =0
对于任意函数上任意一点,其斜率等于其切线与x轴正方向的夹角,即tanα 斜率计算:ax+by+c=0中,k=-a/b.
直线斜率公式:k=(y2-y1)/(x2-x1)
两条垂直相交直线的斜率相乘积为-1:k1*k2=-1
证明:已知一点(x0,y0) 到直线l Ax + By +c =0 的直线距离
PointToLineDis.java
packagecom.mapbar.algorithm;
importcom.vividsolutions.jts.geom.Coordinate;
importcom.vividsolutions.jts.geom.Geometry;
importcom.vividsolutions.jts.geom.GeometryFactory;
importcom.vividsolutions.jts.geom.Point;
importcom.vividsolutions.jts.io.ParseException;
importcom.vividsolutions.jts.io.WKTReader;
/**
* 点到直线,线段的距离
*
* @author chenlly
*
*/
publicclassPointToLineDis {
privatestaticGeometryFactory factory =newGeometryFactory();
privatestaticWKTReader reader =newWKTReader();
// 直线方程一般式 Ax + By + C = 0;
privatedoubleA;
privatedoubleB;
privatedoubleC;
/**
* 求直线方程的一般式
*
* @param point1
* @param point2
* 直线l经过的两个点
*/
publicvoidlineExp(Point point1, Point point2) {
/**
* 由起始点和终止点构成的直线方程一般式的系数A
*/
A = (point1.getY() - point2.getY())
/ Math.sqrt(Math.pow((point1.getY() - point2.getY()), 2)
+ Math.pow((point1.getX() - point2.getX()), 2));
/**
* 由起始点和终止点构成的直线方程一般式的系数
*/
B = (point2.getX() - point1.getX())
/ Math.sqrt(Math.pow((point1.getY() - point2.getY()), 2)
+ Math.pow((point1.getX() - point2.getX()), 2));
/**
* 由起始点和终止点构成的直线方程一般式的系数
*/
C = (point1.getX() * point2.getY() - point2.getX() * point1.getY())
/ Math.sqrt(Math.pow((point1.getY() - point2.getY()), 2)
+ Math.pow((point1.getX() - point2.getX()), 2));
}
/**
* 点到直线方程的距离 此公式需要证明
*
* @param x
* @param y
* @return
*/
publicdoublealLine(doublex,doubley) {
doubled = Math.abs(A * (x) + B * (y) + C)/Math.sqrt(Math.pow(A,2)+Math.pow(B,2));
returnd;
}
/**
* 点到线段的距离 计算线段和计算点到直线的距离类似,不同Y点在包含l线段的投影可能不是线段上的点,投影可能在起点之前或者终点之后
*
* @param x
* @param y
* @return
*/
publicdoublealSegmentDis(doublex,doubley, Point point1, Point point2) {
doubledis =0;
doublea, b, c;
a = lineDis(point1.getX(), point1.getY(), point2.getX(), point2.getY());// 线段的长度
b = lineDis(point1.getX(), point1.getY(), x, y); // point1到点的距离
c = lineDis(point2.getX(), point2.getY(), x, y);//point2到点的距离
if(c + b == a) {// 点在线段上
dis = 0;
returndis;
}
if(c * c >= a * a + b * b) {// 组成直角三角形或钝角三角形,投影在point1延长线上,
dis = b;
returndis;
}
if(b * b >= a * a + c * c) {// 组成直角三角形或钝角三角形,投影在point2延长线上,
dis = c;
returndis;
}
// 组成锐角三角形,则求三角形的高
doublep = (a + b + c) /2;// 半周长
doubles = Math.sqrt(p * (p - a) * (p - b) * (p - c));// 海伦公式求面积
dis = 2* s / a;// 返回点到线的距离(利用三角形面积公式求高)
returndis;
}
/**
* 计算两点之间的距离
* @param x1
* @param y1
* @param x2
* @param y2
* @return
*/
privatedoublelineDis(doublex1,doubley1,doublex2,doubley2) {
doublelineLength =0;
lineLength = Math.sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2));
returnlineLength;
}
publicstaticvoidmain(String[] args)throwsParseException {
PointToLineDis pl = newPointToLineDis();
// 直线经过的两个点
String line = "LINESTRING (-5 6, -4 8)";
Geometry geometry = reader.read(line);
Coordinate[] coords = geometry.getCoordinates();
Point point1 = factory.createPoint(coords[0]);
Point point2 = factory.createPoint(coords[1]);
// 求直线方程一般式
pl.lineExp(point1, point2);
// 点(x,y)
doublex =2;
doubley = -1;
doubled1 = pl.alLine(x, y);
System.out.println("result:"+ d1);// result:9.391485505499116
doublex1 = -3;
doubley1 =3;
doubled2 = pl.alSegmentDis(x1, y1, point1, point2);
System.out.println("result:"+ d2);// result:3.605551275463989
}
}