Merge pull request #4 from TUW-GEO/develop

Develop

requirements.txt

@@ -3,7 +3,7 @@ astropy>=2.0.1
 numpy>=1.13.1
 toolz>=0.8.2
 pyproj>=1.9.5.1
-pytileproj>=0.0.3
+pytileproj>=0.0.4
 
 # gdal - install via conda
 # numba - install via conda

setup.cfg

@@ -40,7 +40,7 @@ packages =
 # py.test options when running `python setup.py test`
 addopts = tests
 
-[pytest]
+[tool:pytest]
 # Options for py.test:
 # Specify command line options as you would do when invoking py.test directly.
 # e.g. --cov-report html (or xml) for html/xml output or --junitxml junit.xml

tests/test_geom_intersect.py

@@ -0,0 +1,60 @@
+from osgeo import ogr, osr
+import pytest
+
+from equi7grid.equi7grid import Equi7Grid
+
+def test_geom_intersect():
+    geom_roi = setup_geom_roi()
+
+    # reference file list
+    ref_tiles = ['EU500M_E012N012T6', 'EU500M_E012N018T6', 'EU500M_E018N006T6', 'EU500M_E024N006T6',
+                 'EU500M_E024N012T6', 'EU500M_E030N006T6', 'EU500M_E030N012T6', 'EU500M_E030N024T6',
+                 'EU500M_E036N006T6', 'EU500M_E036N012T6', 'EU500M_E036N018T6', 'EU500M_E036N024T6',
+                 'EU500M_E036N030T6', 'EU500M_E036N036T6', 'EU500M_E036N042T6', 'EU500M_E042N000T6',
+                 'EU500M_E042N006T6', 'EU500M_E042N012T6', 'EU500M_E042N018T6', 'EU500M_E042N024T6',
+                 'EU500M_E042N030T6', 'EU500M_E042N042T6', 'EU500M_E048N000T6', 'EU500M_E048N006T6',
+                 'EU500M_E048N012T6', 'EU500M_E048N018T6', 'EU500M_E048N024T6', 'EU500M_E048N030T6',
+                 'EU500M_E048N036T6', 'EU500M_E054N000T6', 'EU500M_E054N006T6', 'EU500M_E054N012T6',
+                 'EU500M_E054N018T6', 'EU500M_E054N024T6', 'EU500M_E054N030T6', 'EU500M_E054N036T6',
+                 'EU500M_E054N042T6', 'EU500M_E054N048T6', 'EU500M_E060N000T6', 'EU500M_E060N006T6',
+                 'EU500M_E060N012T6', 'EU500M_E060N018T6', 'EU500M_E060N024T6', 'EU500M_E060N030T6',
+                 'EU500M_E060N036T6', 'EU500M_E060N048T6', 'EU500M_E066N000T6', 'EU500M_E066N006T6',
+                 'EU500M_E066N012T6', 'EU500M_E066N018T6', 'EU500M_E066N024T6', 'EU500M_E066N030T6',
+                 'EU500M_E066N036T6', 'EU500M_E072N000T6', 'EU500M_E072N006T6', 'EU500M_E072N012T6',
+                 'EU500M_E072N018T6', 'EU500M_E072N024T6', 'EU500M_E072N030T6', 'EU500M_E078N006T6',
+                 'EU500M_E078N012T6']
+
+    grid = Equi7Grid(500)
+    res_tiles = grid._search_sgrid_tiles(geom_roi, 'EU', True)
+
+    assert sorted(ref_tiles) == sorted(res_tiles)
+
+def setup_geom_roi():
+    ring_global = ogr.Geometry(ogr.wkbLinearRing)
+    ring_global.AddPoint(-180.0, -89.999928)
+    ring_global.AddPoint(-0.000072, -89.999928)
+    ring_global.AddPoint(179.999856, -89.999928)
+    ring_global.AddPoint(179.999856, 0.000036)
+    ring_global.AddPoint(179.999856, 90.0)
+    ring_global.AddPoint(-0.000072, 90.0)
+    ring_global.AddPoint(-180, 90)
+    ring_global.AddPoint(-180.0, 0.000036)
+    ring_global.AddPoint(-180.0, -89.999928)
+
+    poly_global = ogr.Geometry(ogr.wkbPolygon)
+    poly_global.AddGeometry(ring_global)
+
+    geom_global_wkt = '''GEOGCS[\"WGS 84\",
+                           DATUM[\"WGS_1984\",
+                                 SPHEROID[\"WGS 84\", 6378137, 298.257223563,
+                                          AUTHORITY[\"EPSG\", \"7030\"]],
+                                 AUTHORITY[\"EPSG\", \"6326\"]],
+                           PRIMEM[\"Greenwich\", 0],
+                           UNIT[\"degree\", 0.0174532925199433],
+                           AUTHORITY[\"EPSG\", \"4326\"]]'''
+    geom_global_sr = osr.SpatialReference()
+    geom_global_sr.ImportFromWkt(geom_global_wkt)
+    poly_global.AssignSpatialReference(geom_global_sr)
+
+    return poly_global
+