as I couldn't find any solution for this, I ask for some ideas for my problem on this platform:
I have two raster datasets (a World Soil Map and Global Horizontal Irradiance), which I merged together in QGIS in order to get values of both channels for every pixel.
The Soil Map has 30 soil type classes and the GHI raster has been restricted to 16 classes.
I would like to get an overview of the frequency of each combination (30x16).
The next example shows the values of a random point. As the upper table is simplified, here you get an actual point of the raster: Channel 1 shows the soil type (57 represents Gleysols) and channel 2 gives the value of irradiance (Long-term average of daily sum <1752 kWh/m^2):
Do you have any recommendations how to process a frequency matrix of this? I also want to do this using vector data and compare the results. (My R and Python-Skills are rather limited, but I'm willing to learn)