Multispectral images obtained from Unmanned Aerial Vehicles (UAVs) have become strategic tools in precision agriculture, particularly for analyzing spatial variability in soil attributes. This study aimed to evaluate the spatial distribution of soil electrical (EC) and total organic carbon (TOC) in irrigated forage cactus areas in the Brazilian semiarid region, using field measurements and UAV-based multispectral imagery. The study was conducted in a communal agricultural settlement located in the Mimoso Alluvial Valley (MAV), where EC and TOC were measured at 96 points, and seven biophysical indices were derived from UAV multispectral imagery. Geostatistical models, including cokriging with spectral indices (NDVI, EVI, GDVI, SAVI, and NDSI), were applied to map soil attributes at different spatial scales. Cokriging improved the spatial prediction of EC and TOC by reducing uncertainty and increasing mapping accuracy. The standard deviation of EC decreased from 1.39 (kriging) to 0.67 (cokriging with EVI), and for TOC from 15.55 to 8.78 (cokriging with NDVI and NDSI), reflecting a 43.5% reduction in uncertainty. The indices, EVI, NDVI, and NDSI, showed strong potential in representing and enhancing the spatial variability in soil attributes. NDVI and NDSI were particularly effective at finer grid resolutions, supporting more efficient irrigation strategies and sustainable agricultural practices.