java基础 编辑器IntelliJ IDEA(1).zip

package transform; public class GisCheckUtils { private static final Logger log = new Logger(); /** * 地球周长 */ private static final double L = 6381372 * Math.PI * 2; /** * 平面展开后，x轴等于周长 */ private static final double W = L; /** * y轴约等于周长一半 */ private static final double H = L / 2; /** * 米勒投影中的一个常数，范围大约在正负2.3之间 */ private static final double mill = 2.3; public static void main(String[] args) { /* 拾取坐标系统 可直接复制坐标测试 http://api.map.baidu.com/lbsapi/getpoint/index.html */ //线段一 SegmentLatLon segmentLatLonOne = new SegmentLatLon(); segmentLatLonOne.setStartLatLon(new LatLon(104.514, 30.7787)); segmentLatLonOne.setEndLatLon(new LatLon(104.53, 30.7776)); //线段二 SegmentLatLon segmentLatLonTwo = new SegmentLatLon(); segmentLatLonTwo.setStartLatLon(new LatLon(104.529, 30.7621)); segmentLatLonTwo.setEndLatLon(new LatLon(104.52412,30.78412)); //判断是否相交（延长线相交不算） boolean b = segIntersect(segmentLatLonOne, segmentLatLonTwo); log.info("判断是否相交: " + b); //交点坐标（延长线相交不算） LatLon intersectPoint = getIntersectPoint(segmentLatLonOne, segmentLatLonTwo); log.info("交点坐标：" + intersectPoint); } /** * 将经纬度转换成X和Y轴 * 米勒投影算法 * * @param lat 纬度 * @param lon 经度 * @return */ private static Point latLonToXy(double lat, double lon) { // 将经度从度数转换为弧度 double x = lon * Math.PI / 180; // 将纬度从度数转换为弧度 double y = lat * Math.PI / 180; // 米勒投影的转换 y = 1.25 * Math.log(Math.tan(0.25 * Math.PI + 0.4 * y)); // 弧度转为实际距离 x = (W / 2) + (W / (2 * Math.PI)) * x; y = (H / 2) - (H / (2 * mill)) * y; return new Point(x, y); } /** * xy轴转坐标 * * @param x * @param y * @return 坐标点 */ private static LatLon xyToLatLon(double x, double y) { //实际距离 转为弧度 x = (x - (W / 2)) / (W / (2 * Math.PI)); y = -1 * (y - (H / 2)) / (H / (2 * mill)); // 米勒投影的转换反转 y = (Math.atan(Math.pow(Math.E, y / 1.25)) - 0.25 * Math.PI) / 0.4; //将经度从弧度转换为度数 double lon = 180 / Math.PI * x; //将纬度从弧度转换为度数 double lat = 180 / Math.PI * y; return new LatLon(lon, lat); } /** * 获取两条直线相交的点 * * @param segmentLatLonOne 线段一 * @param segmentLatLonTwo 线段二 * @return 相交点坐标 为null 不相交 */ public static LatLon getIntersectPoint(SegmentLatLon segmentLatLonOne, SegmentLatLon segmentLatLonTwo) { //先判断有没有相交 (延长线不算) 不相交返回null 相交执行获取交点坐标 （如需获取延长线交点注释此方法）---- if (!segIntersect(segmentLatLonOne, segmentLatLonTwo)) { return null; } //---- //转换对象 Point[] points = segmentLatLonToPoint(segmentLatLonOne, segmentLatLonTwo); //获取两条直线相交的点 Point latLon = getIntersectPoint(points[0], points[1], points[2], points[3]); if (latLon != null) { return xyToLatLon(latLon.x, latLon.y); } return null; } /** * 获取两条直线相交的点 * * @param p1 线段一 开始点 * @param p2 线段一 结束点 * @param p3 线段二 开始点 * @param p4 线段二 结束点 * @return */ public static Point getIntersectPoint(Point p1, Point p2, Point p3, Point p4) { double A1 = p1.getY() - p2.getY(); double B1 = p2.getX() - p1.getX(); double C1 = A1 * p1.getX() + B1 * p1.getY(); double A2 = p3.getY() - p4.getY(); double B2 = p4.getX() - p3.getX(); double C2 = A2 * p3.getX() + B2 * p3.getY(); double det_k = A1 * B2 - A2 * B1; if (Math.abs(det_k) < 0.00001) { return null; } double a = B2 / det_k; double b = -1 * B1 / det_k; double c = -1 * A2 / det_k; double d = A1 / det_k; double x = a * C1 + b * C2; double y = c * C1 + d * C2; return new Point(x, y); } /** * 验证两条线有没有相交 * * @param segmentLatLonOne 线段1 * @param segmentLatLonTwo 线段2 * @return true 相交 */ public static boolean segIntersect(SegmentLatLon segmentLatLonOne, SegmentLatLon segmentLatLonTwo) { //转换对象 Point[] points = segmentLatLonToPoint(segmentLatLonOne, segmentLatLonTwo); //验证两条线有没有相交 return segIntersect(points[0], points[1], points[2], points[3]) > 0; } /** * 线段转换为点对象 * * @param segmentLatLonOne 线段一 * @param segmentLatLonTwo 线段二 * @return 点对象数组 */ private static Point[] segmentLatLonToPoint(SegmentLatLon segmentLatLonOne, SegmentLatLon segmentLatLonTwo) { //线段1 Double oneStartLat = segmentLatLonOne.getStartLatLon().getLat(); Double oneStartLon = segmentLatLonOne.getStartLatLon().getLon(); Double oneEndLat = segmentLatLonOne.getEndLatLon().getLat(); Double oneEndLon = segmentLatLonOne.getEndLatLon().getLon(); // 线段2 Double twoStartLat = segmentLatLonTwo.getStartLatLon().getLat(); Double twoStartLon = segmentLatLonTwo.getStartLatLon().getLon(); Double twoEndLat = segmentLatLonTwo.getEndLatLon().getLat(); Double twoEndLon = segmentLatLonTwo.getEndLatLon().getLon(); Point[] points = new Point[4]; //将经纬度转换成X和Y轴 points[0] = latLonToXy(oneStartLat, oneStartLon); points[1] = latLonToXy(oneEndLat, oneEndLon); points[2] = latLonToXy(twoStartLat, twoStartLon); points[3] = latLonToXy(twoEndLat, twoEndLon); return points; } /** * 验证两条线有没有相交 * * @param A 线段一 开始点 * @param B 线段一 结束点 * @param C 线段二 开始点 * @param D 线段二 结束点 * @return */ public static int segIntersect(Point A, Point B, Point C, Point D) { Point intersection = new Point(); if (Math.abs(B.getY() - A.getY()) + Math.abs(B.getX() - A.getX()) + Math.abs(D.getY() - C.getY()) + Math.abs(D.getX() - C.getX()) == 0) { if ((C.getX() - A.getX()) + (C.getY() - A.getY()) == 0) { log.info("ABCD是同一个点！"); } else { log.info("AB是一个点，CD是一个点，且AC不同！"); } return 0; } if (Math.abs(B.getY() - A.getY()) + Math.abs(B.getX() - A.getX()) == 0) { if ((A.getX() - D.getX()) * (C.getY() - D.getY()) - (A.getY() - D.getY()) * (C.getX() - D.getX()) == 0) { log.info("A、B是一个点，且在CD线段上！"); } else { log.info("A、B是一个点，且不在CD线段上！"); } return 0; } if (Math.abs(D.getY() - C.getY()) + Math.abs(D.getX() - C.getX()) == 0) { if ((D.getX() - B.getX()) * (A.getY() - B.getY()) - (D.getY() - B.getY()) * (A.getX() - B.getX()) == 0) { log.info("C、D是一个点，且在AB线段上！"); } else { log.info("C、D是一个点，且不在AB线段上！"); } return 0; } if ((B.getY() - A