polygon.go 2.1 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970
  1. // Copyright 2015 by caixw, All rights reserved.
  2. // Use of this source code is governed by a MIT
  3. // license that can be found in the LICENSE file.
  4. package identicon
  5. var (
  6. // 4个元素分别表示cos(0),cos(90),cos(180),cos(270)
  7. cos = []float64{1, 0, -1, 0}
  8. // 4个元素分别表示sin(0),sin(90),sin(180),sin(270)
  9. sin = []float64{0, 1, 0, -1}
  10. )
  11. // 将points中的所有点,以x,y为原点旋转angle个角度。
  12. // angle取值只能是[0,1,2,3],分别表示[0,90,180,270]
  13. func rotate(points []float64, x, y float64, angle int) {
  14. if angle < 0 || angle > 3 {
  15. panic("rotate:参数angle必须0,1,2,3三值之一")
  16. }
  17. for i := 0; i < len(points); i += 2 {
  18. px := points[i] - x
  19. py := points[i+1] - y
  20. points[i] = px*cos[angle] - py*sin[angle] + x
  21. points[i+1] = px*sin[angle] + py*cos[angle] + y
  22. }
  23. }
  24. // 判断某个点是否在多边形之内,不包含构成多边形的线和点
  25. // x,y 需要判断的点坐标
  26. // points 组成多边形的所顶点,每两个元素表示一点顶点,其中最后一个顶点必须与第一个顶点相同。
  27. func pointInPolygon(x float64, y float64, points []float64) bool {
  28. if len(points) < 8 { // 只有2个以上的点,才能组成闭合多边形
  29. return false
  30. }
  31. // 大致算法如下:
  32. // 把整个平面以给定的测试点为原点分两部分:
  33. // - y>0,包含(x>0 && y==0)
  34. // - y<0,包含(x<0 && y==0)
  35. // 依次扫描每一个点,当该点与前一个点处于不同部分时(即一个在y>0区,一个在y<0区),
  36. // 则判断从前一点到当前点是顺时针还是逆时针(以给定的测试点为原点),如果是顺时针r++,否则r--。
  37. // 结果为:2==abs(r)。
  38. r := 0
  39. x1, y1 := points[0], points[1]
  40. prev := (y1 > y) || ((x1 > x) && (y1 == y))
  41. for i := 2; i < len(points); i += 2 {
  42. x2, y2 := points[i], points[i+1]
  43. curr := (y2 > y) || ((x2 > x) && (y2 == y))
  44. if curr == prev {
  45. x1, y1 = x2, y2
  46. continue
  47. }
  48. mul := (x1-x)*(y2-y) - (x2-x)*(y1-y)
  49. if mul > 0 {
  50. r++
  51. } else if mul < 0 {
  52. r--
  53. }
  54. x1, y1 = x2, y2
  55. prev = curr
  56. }
  57. return r == 2 || r == -2
  58. }