Exp: Add RTree index for polygon (#4)

Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "geometry-rs"
-version = "0.1.2"
+version = "0.2.0"
 edition = "2021"
 license-file = "LICENSE"
 description = "geometry utils"
@@ -15,6 +15,7 @@ keywords = ["geometry"]
 [dependencies]
 
 float_next_after = "1.0.0"
+rtree_rs = "0.1.4"
 
 [dev-dependencies]
 serde = { version = "1.0", features = ["derive"] }

README.md

@@ -3,54 +3,10 @@
 Rewrite parts of [tidwall/geometry](https://github.com/tidwall/geometry) to Rust
 for [ringsaturn/tzf-rs](https://github.com/ringsaturn/tzf-rs).
 
-```toml
-[dependencies]
-geometry-rs = "0.1.2"
+```bash
+cargo add geometry-rs
 ```
 
-```rust
-use std::vec;
+## TODO
 
-use geometry_rs;
-
-fn main() {
-    let poly = geometry_rs::Polygon::new(
-        vec![
-            geometry_rs::Point {
-                x: 90.48826291293898,
-                y: 45.951129815858565,
-            },
-            geometry_rs::Point {
-                x: 90.48826291293898,
-                y: 27.99437617512571,
-            },
-            geometry_rs::Point {
-                x: 122.83201291294,
-                y: 27.99437617512571,
-            },
-            geometry_rs::Point {
-                x: 122.83201291294,
-                y: 45.951129815858565,
-            },
-            geometry_rs::Point {
-                x: 90.48826291293898,
-                y: 45.951129815858565,
-            },
-        ],
-        vec![],
-    );
-
-    let p_out = geometry_rs::Point {
-        x: 130.74216916294148,
-        y: 37.649011392900306,
-    };
-
-    print!("{:?}\n", poly.contains_point(p_out));
-
-    let p_in = geometry_rs::Point {
-        x: 99.9804504129416,
-        y: 39.70716466970461,
-    };
-    print!("{:?}\n", poly.contains_point(p_in));
-}
-```
+- [ ] Use compressed data of RTree index, since too many memory costs.

benches/bench.rs

@@ -31,6 +31,7 @@ mod benches_polygon {
                 },
             ],
             vec![],
+            false,
         );
         return poly;
     }

benches/bench_az.rs

@@ -5,7 +5,28 @@ mod benches_az_polygon {
     extern crate test;
     use test::Bencher;
 
-    // test data copy from https://github.com/unitedstates/districts/blob/gh-pages/states/AZ/shape.geojson
+    /* Test data copy from https://github.com/unitedstates/districts/blob/gh-pages/states/AZ/shape.geojson
+
+    Python code to generate the data:
+        ```python
+        import json
+
+        TPL = """geometry_rs::Point {x: {lng}, y: {lat}},"""
+        
+        with open("./az.geojson") as f:
+            data = json.loads(f.read())
+        
+        coordinates = data["coordinates"][0][0]
+        # print(len(coordinates))
+        gens = []
+        for coord in coordinates:
+            # gens.append(TPL.format(lng=coord[0], lat=coord[1]))
+            gens.append(TPL.replace("{lng}", str(coord[0])).replace("{lat}", str(coord[1])))
+        gen_text = "\n".join(gens)
+        print(gen_text)
+        ```
+    */
+
     fn load_poly() -> geometry_rs::Polygon{
         let poly = geometry_rs::Polygon::new(
             vec![
@@ -1594,6 +1615,7 @@ mod benches_az_polygon {
 
             ],
             vec![],
+            true,
         );
         return poly;
     }