Over 85% of all chemical industry products are made using catalysts (1), with the overwhelming majority of these employing heterogeneous catalysts (2) functioning at the gas-solid interface (3). Consequently, optimizing catalytic reactor design attracts much effort. Such optimization relies on heat transfer and fluid dynamics modeling coupled to surface reaction kinetics (4). The complexity of these systems demands many approximations, which can only be tested with experimental observations (5,6) of quantities such as temperature, pressure, concentrations, flow rates, etc. One essential measurement is a map of the spatial variation in temperature throughout the catalyst bed. We present here the first non-invasive maps of gas temperatures in catalyst-filled reactors, including high spatial resolution maps in microreactors enabled by parahydrogen. The thermal maps reveal energy flux patterns whose length scale correlates with the catalyst packing.