fix(hairy-tiles): fix wrong recognition of no-data values (#44)

Author: vwvbrand

readme.md

@@ -177,7 +177,7 @@ Error: /src/data/cat_aggr_buf_30m/input/vector/osm_merged_2012.gpkg_transformed.
 **Solution**: run the command again with the same parameters.
 - **Warning**: if you delete stressors, you would have to calculate them again using the 3d command to proceed with the 4th command.
 
-Not each issue is known so far, so if you have faced something different, please let us know by opening a new issue.
+Not each issue is known so far, so if you have faced something different, please let us know by opening a new issue. **All relevant issues are mentioned** [**here**](https://github.com/AD4GD/pilot-2/issues). 
 
 #### Further development
 This tool is mostly complete, but a few improvements could be made if there is demand from users.

src/data/shared/input/osm_data/case_study_albera/gpkg_temp/osm_merged_2017.gpkg

@@ -447,6 +447,9 @@ def mask_raster_with_raster(self, input_raster, mask_raster, nodata_value, outpu
         # get raster properties
         in_band = in_ds.GetRasterBand(1)
         mask_band = mask_ds.GetRasterBand(1)
+        mask_nodata_val = mask_band.GetNoDataValue()
+        print(f"Nodata value of the masking raster dataset: {mask_nodata_val}")
+
         out_driver = gdal.GetDriverByName('GTiff')
         if output_raster is None:
             out_ds = gdal.Open(input_raster, gdal.GA_Update)
@@ -458,7 +461,7 @@ def mask_raster_with_raster(self, input_raster, mask_raster, nodata_value, outpu
 
         # create a mask array from the mask raster
         mask_data = mask_band.ReadAsArray()
-        mask_data[mask_data == 0] = nodata_value
+        mask_data[mask_data == mask_nodata_val] = nodata_value
 
         # apply the mask to the input raster
         in_data = in_band.ReadAsArray()