polygon, _ := geoencoding.Decode([]byte(`MULTIPOLYGON(((121.308225482526 31.3977223795235,121.320481 31.399024,121.322129 31.386769,121.304536 31.385516,121.300412 31.393475,121.308225482526 31.3977223795235)),((121.308225482526 31.3977223795235,121.308156 31.397715,121.308337 31.397783,121.308225482526 31.3977223795235)))`), geoencoding.WKT)
line, _ := geoencoding.Decode([]byte(`linestring(121.318481 31.40068, 121.317523 31.393429, 121.324714 31.392735, 121.325045 31.390363)`), geoencoding.WKT)
a := planar.NewMegrezAlgorithm()
i, e := a.Intersects(polygon, line)
fmt.Println(i, e)
