多边形中心坐标的计算

2024-06-20 16:18
文章标签 计算 多边形 坐标 中心

本文主要是介绍多边形中心坐标的计算,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

import java.util.ArrayList;
import java.util.List;
/**
 * 坐标以及电子围栏相关的工具类
 *
 */
public class PointUtil {
    
    public static void main(String[] args) {
        String str="114.316587,30.671626#114.325921,30.683437#114.342122,30.689342#114.359009,30.695246#114.377066,30.692884#114.393304,30.678417#114.412289,30.650953#114.426107,30.625253#114.433703,30.595853#114.444368,30.560458#114.466705,30.527415#114.457264,30.517327#114.454002,30.503097#114.458465,30.474631#114.430828,30.466313#114.397697,30.460362#114.343281,30.456778#114.305086,30.457354#114.269101,30.468743#114.23449,30.489598#114.182713,30.522872#114.17296,30.535364#114.167328,30.55022#114.165783,30.575754#114.167155,30.591715#114.161492,30.606935#114.160461,30.619049#114.170589,30.627469#114.206123,30.623481#114.233932,30.629832#114.184837,30.704361#114.187927,30.776472#114.233589,30.780122#114.238567,30.717035#114.248609,30.66835#114.284744,30.659813#";

    }
     /**
      * 传入字符串类型的坐标组;经纬度“,”分隔;坐标之间“#”分隔;以“#”好结束
      * @return
      */
     public static boolean is(String str,double x,double y){
         List<Point> list = getxyUtil(str);
         Point[] ps2 = new Point[list.size()];     
         for(int i=0;i<list.size();i++){     
             ps2[i]=list.get(i);
         }     
         Point n1 = new Point(x,y);  
         if(isPtInPoly(n1.getX() , n1.getY() , ps2)){     
             return true;
         }
         return false;
     }
     /**
      * 传入经度,维度;和需要判断的多边形坐标组
      * @param ALon
      * @param ALat
      * @param ps
      * @return
      */
        public static boolean isPtInPoly (double ALon , double ALat , Point[] ps) {  
            int iSum, iCount, iIndex;  
            double dLon1 = 0, dLon2 = 0, dLat1 = 0, dLat2 = 0, dLon;  
            if (ps.length < 3) {  
                return false;  
            }  
            iSum = 0;  
            iCount = ps.length;  
            for (iIndex = 0; iIndex<iCount;iIndex++) {  
                if (iIndex == iCount - 1) {  
                    dLon1 = ps[iIndex].getX();  
                    dLat1 = ps[iIndex].getY();  
                    dLon2 = ps[0].getX();  
                    dLat2 = ps[0].getY();  
                } else {  
                    dLon1 = ps[iIndex].getX();  
                    dLat1 = ps[iIndex].getY();  
                    dLon2 = ps[iIndex + 1].getX();  
                    dLat2 = ps[iIndex + 1].getY();  
                }  
                // 以下语句判断A点是否在边的两端点的水平平行线之间,在则可能有交点,开始判断交点是否在左射线上  
                if (((ALat >= dLat1) && (ALat < dLat2)) || ((ALat >= dLat2) && (ALat < dLat1))) {  
                    if (Math.abs(dLat1 - dLat2) > 0) {  
                        //得到 A点向左射线与边的交点的x坐标:  
                        dLon = dLon1 - ((dLon1 - dLon2) * (dLat1 - ALat) ) / (dLat1 - dLat2);  
                        // 如果交点在A点左侧(说明是做射线与 边的交点),则射线与边的全部交点数加一:  
                        if (dLon < ALon) {  
                            iSum++;  
                        }  
                    }  
                }  
            }  
            if ((iSum % 2) != 0) {  
                return true;  
            }  
            return false;  
        }  
        /**
         * 传入字符串坐标组;转为集合;114.288221,30.593268#114.321866,30.580559#114.311223,30.560015#
         * @param str
         * @return
         */
        public static List<Point> getxyUtil(String str){
            List<Point> list = new ArrayList<>();
            while (str.indexOf("#")!=-1) {
                String str1 = str.substring(0, str.indexOf("#"));//40.01,40.12
                Point pd = new Point(Double.valueOf(str1.substring(0, str1.indexOf(","))),Double.valueOf(str1.substring(str1.indexOf(",")+1)));
                list.add(pd);
                str=str.substring(str.indexOf("#")+1);
            }    
            while (str.indexOf(";")!=-1) {
                String str1 = str.substring(0, str.indexOf(";"));//40.01,40.12
                Point pd = new Point(Double.valueOf(str1.substring(0, str1.indexOf(","))),Double.valueOf(str1.substring(str1.indexOf(",")+1)));
                list.add(pd);
                str=str.substring(str.indexOf(";")+1);
            }    
            return list;
        } 
        /**
         * 获取多边形的中心
         * @param str
         * @return
         */
        public static PageData getCenterOfGravityPoint(String str ) { 
            PageData data = new PageData();
             List<Point> mPoints = getxyUtil(str);
            double area = 0.0;//多边形面积  
            double Gx = 0.0, Gy = 0.0;// 重心的x、y  
            for (int i = 1; i <= mPoints.size(); i++) {  
                double iLat = mPoints.get(i % mPoints.size()).getX();  
                double iLng = mPoints.get(i % mPoints.size()).getY();  
                double nextLat = mPoints.get(i - 1).getX();  
                double nextLng = mPoints.get(i - 1).getY();  
                double temp = (iLat * nextLng - iLng * nextLat) / 2.0;  
                area += temp;  
                Gx += temp * (iLat + nextLat) / 3.0;  
                Gy += temp * (iLng + nextLng) / 3.0;  
            }  
            Gx = Gx / area;  
            Gy = Gy / area; 
            data.put("x", String.valueOf(Gx).substring(0, String.valueOf(Gx).indexOf(".")+7));
            data.put("y", String.valueOf(Gy).substring(0, String.valueOf(Gy).indexOf(".")+7));
            return data;  
        } 
    
    
}

class  Point {
    private Double x;  
    private Double y;  
    public Point (Double x , Double y) {  
        this.x = x;  
        this.y = y;  
    }  
    public Double getX() {  
        return x;  
    }  
    public void setX(Double x) {  
        this.x = x;  
    }  
    public Double getY() {  
        return y;  
    }  
    public void setY(Double y) {  
        this.y = y;  
    }      
}

这篇关于多边形中心坐标的计算的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1078629

相关文章

Spring Cloud之注册中心Nacos的使用详解

《SpringCloud之注册中心Nacos的使用详解》本文介绍SpringCloudAlibaba中的Nacos组件,对比了Nacos与Eureka的区别,展示了如何在项目中引入SpringClo... 目录Naacos服务注册/服务发现引⼊Spring Cloud Alibaba依赖引入Naco编程s依

Python如何计算两个不同类型列表的相似度

《Python如何计算两个不同类型列表的相似度》在编程中,经常需要比较两个列表的相似度,尤其是当这两个列表包含不同类型的元素时,下面小编就来讲讲如何使用Python计算两个不同类型列表的相似度吧... 目录摘要引言数字类型相似度欧几里得距离曼哈顿距离字符串类型相似度Levenshtein距离Jaccard相

使用C#代码计算数学表达式实例

《使用C#代码计算数学表达式实例》这段文字主要讲述了如何使用C#语言来计算数学表达式,该程序通过使用Dictionary保存变量,定义了运算符优先级,并实现了EvaluateExpression方法来... 目录C#代码计算数学表达式该方法很长,因此我将分段描述下面的代码片段显示了下一步以下代码显示该方法如

如何用Java结合经纬度位置计算目标点的日出日落时间详解

《如何用Java结合经纬度位置计算目标点的日出日落时间详解》这篇文章主详细讲解了如何基于目标点的经纬度计算日出日落时间,提供了在线API和Java库两种计算方法,并通过实际案例展示了其应用,需要的朋友... 目录前言一、应用示例1、天安门升旗时间2、湖南省日出日落信息二、Java日出日落计算1、在线API2

跨国公司撤出在华研发中心的启示:中国IT产业的挑战与机遇

近日,IBM中国宣布撤出在华的两大研发中心,这一决定在IT行业引发了广泛的讨论和关注。跨国公司在华研发中心的撤出,不仅对众多IT从业者的职业发展带来了直接的冲击,也引发了人们对全球化背景下中国IT产业竞争力和未来发展方向的深思。面对这一突如其来的变化,我们应如何看待跨国公司的决策?中国IT人才又该如何应对?中国IT产业将何去何从?本文将围绕这些问题展开探讨。 跨国公司撤出的背景与

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

Codeforces Round #113 (Div. 2) B 判断多边形是否在凸包内

题目点击打开链接 凸多边形A, 多边形B, 判断B是否严格在A内。  注意AB有重点 。  将A,B上的点合在一起求凸包,如果凸包上的点是B的某个点,则B肯定不在A内。 或者说B上的某点在凸包的边上则也说明B不严格在A里面。 这个处理有个巧妙的方法,只需在求凸包的时候, <=  改成< 也就是说凸包一条边上的所有点都重复点都记录在凸包里面了。 另外不能去重点。 int