Apparent resistivity imaging is a typical rapid imaging method in the ground-airborne frequency-domain electromagnetic method. At present, the apparent resistivity is typically calculated by the measured magnetic field, however, this imaging method exhibits limited capability in recognizing the center of three-dimensional anomalies. Therefore, this paper proposed the calculation of apparent resistivity using magnetic field gradients. To solve the problem of random artificial anomalies that existed during the calculated process, this paper presents a hybrid least square support vector machine (LSSVM) and Northern Goshawk optimization (NGO) to establish the mapping relationship between the magnetic field gradient and apparent resistivity variation. This approach enables accurate prediction of apparent resistivity variations and effectively resolves the challenge of correcting background resistivity. Furthermore, three typical theoretical models and field examples are used to predict the apparent resistivity variation, the imaging results demonstrate that the proposed NGO-LSSVM algorithm is a feasible and efficient tool for predicting the apparent resistivity variation with high accuracy. This study provides a novel and efficient imaging method, which facilitates the application of ground-airborne frequency-domain electromagnetics for high-resolution detection requirements, such as mineral exploration.