This dataset contains three R scripts used for spatial analysis and modeling in the study on greenhouse gas emissions from drained peatlands in the European Union. Mainscript_Peatcover_maps.R: Processes spatial data to generate high-resolution peatland cover maps across EU+ countries. It integrates multiple geospatial datasets and applies masking and classification procedures to identify peatland areas. Emission_script.R: Calculates greenhouse gas emissions from drained peatlands using peatland cover data and emission factor allocations. It employs parallel processing techniques to efficiently estimate emissions across large spatial domains. Hotspot_map.R: Generates hotspot maps of greenhouse gas emissions by aggregating spatial emission data, applying statistical downscaling, and visualizing emissions using "bi-variate" mapping techniques. The script integrates multiple geospatial sources to refine emission estimates. These scripts form the computational backbone for mapping and quantifying peatland-related emissions in the EU+ region.