panlei
2019-11-16 0d95e05cf47aa7b24a57a673db46d5607016016c
ruleserver/geoPolygon.go
@@ -2,6 +2,7 @@
import (
   "math"
   "ruleprocess/structure"
)
func min(num1 float64, num2 float64) float64 {
@@ -17,12 +18,6 @@
      return num2
   }
   return num1
}
//Point 坐标点
type Point struct {
   X float64 `json:"x"`
   Y float64 `json:"y"`
}
//Rect 检测目标
@@ -42,7 +37,7 @@
//PintIsInPolygon 判断点是否在多边形内部
//point为要判断的坐标点
//polygon是多边形各点数组
func pintIsInPolygon(point Pointfloat, polygon []Point, widthScale float64, heightScale float64) bool {
func PintIsInPolygon(point structure.Pointfloat, polygon []structure.Point, widthScale float64, heightScale float64) bool {
   var nCross int = 0
   for i := 0; i < len(polygon); i++ {
@@ -72,17 +67,17 @@
}
//GetLocation 将一个给定起始坐标,宽度长度的矩形区域均分为n方份并返回中心坐标(n为单边平分数值)和面积
func getLocation(rect Rect, n int) ([]Pointfloat, float64) {
func GetLocation(rect structure.Rect, n int) ([]structure.Pointfloat, float64) {
   xArr := make([]float64, n) // 用切片不用数组,数组不能用变量定义长度
   yArr := make([]float64, n)
   pointArr := make([]Pointfloat, 0, n*n)
   pointArr := make([]structure.Pointfloat, 0, n*n)
   for i := 0; i < n; i++ {
      xArr[i] = rect.X + (rect.Width/float64(2*n))*float64(2*i+1)
      yArr[i] = rect.Y + (rect.Height/float64(2*n))*float64(2*i+1)
   }
   for i := 0; i < n; i++ {
      for j := 0; j < n; j++ {
         point := Pointfloat{X: xArr[i], Y: yArr[j]}
         point := structure.Pointfloat{X: xArr[i], Y: yArr[j]}
         pointArr = append(pointArr, point)
      }
   }
@@ -91,7 +86,7 @@
}
//ComputePolygonArea 计算任意多边形面积
func computePolygonArea(polygon []Point) float64 {
func ComputePolygonArea(polygon []structure.Point) float64 {
   pointNum := len(polygon)
   var s float64 = 0
   if pointNum < 3 {
@@ -105,18 +100,18 @@
}
//PgsInterPercent calculate percent of two polygon intersection  计算两个多边形的重叠占比
func PgsInterPercent(pgpts []Point, box Rect, widthScale float64, heightScale float64) (percent float64) {
func PgsInterPercent(pgpts []structure.Point, box structure.Rect, widthScale float64, heightScale float64) (percent float64) {
   areapts, areaBox := getLocation(box, 10)
   areapts, areaBox := GetLocation(box, 10)
   var count = 0
   for _, pts := range areapts {
      if pintIsInPolygon(pts, pgpts, widthScale, heightScale) {
      if PintIsInPolygon(pts, pgpts, widthScale, heightScale) {
         count++
      }
   }
   perInterBox := float64(count) / float64(len(areapts)) // 重合面积占矩形的比例
   areaInter := perInterBox * areaBox
   areaPg := computePolygonArea(pgpts)
   areaPg := ComputePolygonArea(pgpts)
   perInterPg := areaInter / areaPg // 重合面积占多边形区域的比例
   // 哪个占的比例大按哪个计算
   if perInterBox > perInterPg {