A diffusion kurtosis imaging method, comprising: acquiring MR signals of a scanned object in multiple directions with a non-zero b value representing diffusion intensity, each direction of the multiple directions being represented with a direction vector; fitting the MR signals using an unconstrained optimization algorithm to obtain elements of a first diffusion tensor and elements of a first kurtosis tensor; determining diffusion tensor imaging parameters and/or kurtosis tensor imaging parameters based on the elements of the first diffusion tensor and the elements of the first kurtosis tensor, wherein the diffusion tensor imaging parameters include at least one of: the diffusion coefficients in each direction, a mean diffusion coefficient, an axial diffusion coefficient, a radial diffusion coefficient, or a diffusion anisotropy coefficient, and the kurtosis tensor imaging parameters include at least one of: the kurtosis coefficients in each direction, a mean kurtosis coefficient, an axial kurtosis coefficient, a radial kurtosis coefficient, or a kurtosis anisotropy coefficient; and generating a diffusion parameter image based on the diffusion tensor imaging parameters and/or generating a kurtosis parameter image based on the kurtosis tensor imaging parameters, characterized in that determining the diffusion tensor imaging parameters and/or the kurtosis tensor imaging parameters based on the elements of the first diffusion tensor and the elements of the first kurtosis tensor comprises: calculating (S201, S401, S601, S701, S1001) the diffusion coefficients in each direction and/or the kurtosis coefficients in each direction based on the elements of the first diffusion tensor and the elements of the first kurtosis tensor, correcting (S202, S302, S402, S503, S603, S702, S802, S1002, S1004) the diffusion coefficients and/or the kurtosis coefficients outside preset ranges to obtain elements of a second diffusion tensor and/or elements of a second kurtosis tensor, and determining (S203, S303, S403, S604, S703, S1005) the diffusion tensor imaging parameters and/or the kurtosis tensor imaging parameters based on the elements of the second diffusion tensor and/or the elements of the second kurtosis tensor. acquiring MR signals of a scanned object in multiple directions with a non-zero b value representing diffusion intensity, each direction of the multiple directions being represented with a direction vector; fitting the MR signals using an unconstrained optimization algorithm to obtain elements of a first diffusion tensor and elements of a first kurtosis tensor; determining diffusion tensor imaging parameters and/or kurtosis tensor imaging parameters based on the elements of the first diffusion tensor and the elements of the first kurtosis tensor, wherein the diffusion tensor imaging parameters include at least one of: the diffusion coefficients in each direction, a mean diffusion coefficient, an axial diffusion coefficient, a radial diffusion coefficient, or a diffusion anisotropy coefficient, and the kurtosis tensor imaging parameters include at least one of: the kurtosis coefficients in each direction, a mean kurtosis coefficient, an axial kurtosis coefficient, a radial kurtosis coefficient, or a kurtosis anisotropy coefficient; and generating a diffusion parameter image based on the diffusion tensor imaging parameters and/or generating a kurtosis parameter image based on the kurtosis tensor imaging parameters, characterized in that determining the diffusion tensor imaging parameters and/or the kurtosis tensor imaging parameters based on the elements of the first diffusion tensor and the elements of the first kurtosis tensor comprises: calculating (S201, S401, S601, S701, S1001) the diffusion coefficients in each direction and/or the kurtosis coefficients in each direction based on the elements of the first diffusion tensor and the elements of the first kurtosis tensor, correcting (S202, S302, S402, S503, S603, S702, S802, S1002, S1004) the diffusion coefficients and/or the kurtosis coefficients outside preset ranges to obtain elements of a second diffusion tensor and/or elements of a second kurtosis tensor, and determining (S203, S303, S403, S604, S703, S1005) the diffusion tensor imaging parameters and/or the kurtosis tensor imaging parameters based on the elements of the second diffusion tensor and/or the elements of the second kurtosis tensor.