Solid Immersion (SI) microscopy is a modern imaging modality that overcomes the Abbe diffraction limit and offers novel applications in various branches of visible, infrared, terahertz, and millimeter-wave optics. Despite the widespread use, SI microscopy usually results in qualitative imaging. Indeed, it presents only the raw distributions (in the image plane) of the backscattered field intensity, while unlocking the information about the physical properties of an imaged object, such as its complex refractive index (RI) distribution, requires resolving the inverse problem and remains a daunting task. In this paper, a method for resolving the SI microscopy inverse problem is developed, capable of reconstructing the RI distribution at the object imaging plane with subwavelength spatial resolution, while performing only intensity measurements. The sample RI is retrieved via minimization of the error function that characterizes discrepancy between the experimental data and the predictions of analytical model. This model incorporates all the key features of the electromagnetic-wave interaction with the SI lens and an imaged object, including contributions of the evanescent and ordinary-reflected waves, as well as effects of light polarization and wide beam aperture. The model is verified numerically, using the finite-element frequency-domain method, and experimentally, using the in-house reflection-mode continuous-wave terahertz SI microscope. Spatial distributions of the terahertz RIs of different low-absorbing optical materials and highly absorbing biological objects were studied and compared to a priori known data to demonstrate the potential of the novel SI microscopy modality. Given the linear nature of the Maxwell’s equations, the developed method can be applied for subwavelength-resolution SI microscopy at other spectral ranges.