Right now only res_high = [1/240 1/240] is possible.
For this to work correctly, we have to rewrite many parts of the code.
A_lu will not be a 2-D matrix any more, but a 3-D matrix, where each pixel is replaced by a 1x17 vector that counts the occurrence of certain land use types. In the trivial case of res_high = [1/240 1/240], those vectors will be filled with 16 zeros and 1 one. Or it will be replaced with the most common land use type (easy).
For float values (elevation, bathymetry), we would calculate the average.
Right now only res_high = [1/240 1/240] is possible. For this to work correctly, we have to rewrite many parts of the code. A_lu will not be a 2-D matrix any more, but a 3-D matrix, where each pixel is replaced by a 1x17 vector that counts the occurrence of certain land use types. In the trivial case of res_high = [1/240 1/240], those vectors will be filled with 16 zeros and 1 one. Or it will be replaced with the most common land use type (easy). For float values (elevation, bathymetry), we would calculate the average.