am using the Invest Urban Cooling module and noticed something strange when checking the intermediate results. The dc_park_ini.tif
file seems to be calculated using green_area
with an exponential distance weighting, rather than using cci
. According to the documentation, the formula should be:
CCparki=∑j∈radius from igj⋅CCj⋅e(−d(i,j)dcool)CC_{park_i} = \sum_{j \in \text{radius from } i} g_j \cdot CC_j \cdot e^{\left( \frac{-d(i,j)}{d_{cool}} \right)}CCparki=j∈radius from i∑gj⋅CCj⋅e(dcool−d(i,j))
However, after further inspection, I found the following code snippet in the source on GitHub:
python
Copy code
green_area_decay_kernel_distance = int(numpy.round(
float(args['green_area_cooling_distance']) / cell_size))
cc_park_raster_path = os.path.join(
intermediate_dir, f'cc_park{file_suffix}.tif')
cc_park_task = task_graph.add_task(
func=convolve_2d_by_exponential,
args=(
green_area_decay_kernel_distance,
task_path_prop_map['green_area'][1],
cc_park_raster_path),
target_path_list=[cc_park_raster_path],
dependent_task_list=[
task_path_prop_map['green_area'][0]],
task_name='calculate T air')
This seems to indicate that green_area
is being used in the calculation, not cci
. Could I be misunderstanding the code? I would appreciate some clarification.