The dose data produced by treatment plan system (TPS) in intensity-modulated radiation therapy (IMRT) has many gradient edge points. Considering this feature we proposed a new interpolation algorithm called treatment plan dose interpolation algorithm based on gradient feature in intensity-modulated radiation therapy (TDAGI), which improves the Canny algorithm to detect the gradient edge points and non-edge points by using the gradient information in the dose data plane. For each gradient edge point, the corresponding gradient profile was traced and the profile's sharpness was calculated, and for each non-edge point, the dispersion was calculated. With the sharpness or dispersion, the kernel coefficients of bi-cubic interpolation can be obtained and can be used as the central point to complete the bi-cubic interpolation calculation. Compared with bi-cubic interpolation and bilinear interpolation, the TDAGI algorithm is more accurate. Furthermore, the TDAGI algorithm has the advantage of gradient keeping. Therefore, TDAGI can be used as an alternative method in the dose interpolation of TPS in IMRT.