The invention discloses a joint inversion method based on magnetotelluric and DC resistivity data, and the method comprises the steps: determining an inversion region according to observation data, carrying out the subdivision to obtain an initial grid, and setting an initial value of an inversion parameter vector; performing Gaussian-Newton inversion iteration based on the current value of the inversion parameter vector to calculate a model change amount and a linear search step length, and then calculating an iterative update value of the inversion parameter vector; judging whether an inversion termination condition is met or not, and if not, carrying out next iteration; if so, judging whether a progressive grid inversion termination condition is reached or not, and if so, ending the inversion; and if not, carrying out refined subdivision on the grid cell, updating the inversion parameter vector value, and returning to carry out Gauss-Newton inversion iteration until inversion is finished. According to the invention, a joint inversion technology of inversion grid adaptive adjustment is researched and developed, the problem of multiplicity of solutions of joint inversion interpretation of a magnetotelluric method and direct current resistivity data is effectively reduced, and the joint interpretation accuracy is improved.