In this paper, a ghost fluid thermal lattice Boltzmann method is developed to simulate Dirichlet and Neumann thermal boundary conditions at curved boundaries. As such, a new formulation for both thermal boundary conditions is developed using a bilinear interpolation method. The presented method is also formulated to address the special cases that arise when the values of the macroscopic variables are interpolated at the image points surrounded by many solid nodes as well as the fluid nodes. The results of the presented method are compared to those available in the literature from conventional numerical methods, and excellent agreement is observed.