<reponame>jorisvandenbossche/DS-python-geospatial
suitable_locations = land_use_residential * (dem_gent > 10) * (1 - roads_buffer_arr)
suitable_locations = land_use_residential * (dem_gent > 10) * (1 - roads_buffer_arr)