I would like to generate a spatially resolved solar seasonality factor based on raster datasets from the Global Solar Atlas (https://globalsolaratlas.info/map) in QGIS. The data basis are average global horizontal irradiation GHI per month as twelve raster data sets for a specific country. Now I want to divide the irradiance of the highest month by the irradiance of the lowest month for each raster cell to get a seasonality factor.
However, this does not seem to be possible within the raster calculator. What is the best way to perform such an operation with multiple raster datasets in QGIS? Do I have to convert them into vector datasets and perform spatial joins or is there an easier direct way for raster files?