use same resampling for mask (#106)

* use same resampling for mask

* update version [skip ci]

CHANGES.txt

@@ -1,4 +1,10 @@
-1.2.5 (2019-05-16)
+1.2.6 (2019-05-10)
+------------------
+
+- Use same resampling method for mask and for data (#105)
+
+
+1.2.5 (2019-05-08)
 ------------------
 
 - add tile_edge_padding option to be passed to rio_tiler.utils._tile_read to reduce sharp edges that occur due to resampling (#104)

rio_tiler/__init__.py

@@ -1,3 +1,3 @@
 """rio-tiler."""
 
-__version__ = "1.2.5"
+__version__ = "1.2.6"

rio_tiler/utils.py

@@ -350,7 +350,11 @@ def _tile_read(
             resampling=Resampling[resampling_method],
         )
 
-        mask = vrt.dataset_mask(out_shape=(tilesize, tilesize), window=out_window)
+        mask = vrt.dataset_mask(
+            out_shape=(tilesize, tilesize),
+            window=out_window,
+            resampling=Resampling[resampling_method],
+        )
 
         return data, mask
 

tests/test_utils.py

@@ -738,3 +738,19 @@ def test_aligned_with_internaltile():
 
     with rasterio.open(COG_WEB_TILED) as src_dst:
         assert utils._requested_tile_aligned_with_internal_tile(src_dst, bounds, 256)
+
+
+def test_tile_read_validMask():
+    """Dataset mask should be the same as the actual mask."""
+    address = "{}_B2.TIF".format(LANDSAT_PATH)
+
+    bounds = (
+        -8844681.416934313,
+        3757032.814272982,
+        -8766409.899970293,
+        3835304.331237001,
+    )
+    tilesize = 128
+    arr, mask = utils.tile_read(address, bounds, tilesize, nodata=0)
+    masknodata = (arr[0] != 0).astype(np.uint8) * 255
+    np.testing.assert_array_equal(mask, masknodata)