scientific_documentation.md

@@ -50,7 +50,9 @@ The computation of flood damages is simple for road areas and green areas. In th
 ## Finding flooded buildings
 When identifying flooded buildings, we often have the challenge that the simulated flood areas contain "holdes" in building locations, because in the model water cannot flow through the buildings. A direct overlay of building and flood depth raster will thus have the result that no buildings are considered flooded. 
 
-To avoid this issue, the flood depth in buildings is here determined using a modified flood depth raster. First, a $`k x k`$ max neighborhood filter is applied to the flood depth raster. For each pixel, the filter finds the maximal flood depth from all pixels in a $`k x k`$ neighborhood around the considered pixel and applies it to the pixel under consideration. In the figure (topright subfigure), the result of this process is shown for a $`3 x 3`$ max filter. Subsequently, a modified flood raster is created where the values from the filtered raster are applied in building locations and the values from the original raster in all other locations. This process is illustrated in the bottomright subfigure. 
+To avoid this issue, the flood depth in buildings is here determined using a modified flood depth raster. First, a $`k x k`$ max neighborhood filter is applied to the flood depth raster. For each pixel, the filter finds the maximal flood depth from all pixels in a $`k x k`$ neighborhood around the considered pixel and applies it to the pixel under consideration. In the figure (topright subfigure), the result of this process is shown for a $`3 x 3`$ max filter. In the program, a filter size of `(floor(float(Dx)/float(resolution)/2)+1)*2` is applied, where `Dx` is the pixel size of the original flood depth raster and `resolution` is the pixel length of the rasters used for the computation of flood damages. 
+
+Subsequently, a modified flood raster is created where the values from the filtered raster are applied in building locations and the values from the original raster in all other locations. This process is illustrated in the bottomright subfigure. 
 
 ![Replacing the values of the flood depth raster in building locations by the maximal flood depths around the building.](img/3x3_Filter.png "Replacing the values of the flood depth raster in building locations by the maximal flood depths around the building.")