In many cases, 1D inversion is still an important step in transient electromagnetic data processing. Potential issues may arise in the calculation of apparent resistivity using induced electromotive force (EMF) due to overshoot and the presence of multi-valued functions. Obtaining reliable and consistent inversion results using a uniform half-space as the initial model is challenging, especially when aiming for efficient inversion. Focusing on these problems, we use the land-based transient electromagnetic (TEM) sounding data, which was acquired by using a large fixed-loop transmitter, and adopt a quasi-2D inversion scheme to generate improved images of the subsurface resistivity structure. First, we have considered directly using magnetic field data or converting induced EMF into magnetic field, and then calculating the apparent resistivity over the whole zone. Next, a resistivity profile that varies with depth is obtained through fast smoke ring imaging. This profile serves as the initial model for the subsequent optimal inversion. The inversion scheme uses a nonlinear least-squares method, incorporating lateral and vertical constraints, to produce a quasi-2D subsurface image. The potentiality of the proposed methodology has been exemplified through the interpretation of synthetic data derived from a 3D intricate resistivity model, as well as field data obtained from a TEM survey conducted in a coalmine field. In both cases, the inversion process yields quasi-2D subsurface images that exhibit a reasonable level of accuracy. These images appear to be less moulded by 3D effects and demonstrate a satisfactory level of agreement with the known target area.