We propose a method to estimate temperature distribution non-destructively using combination of near-field acoustical holography (NAH) method and optimization method. NAH method is valid to calculate forward and back propagations of ultrasonic wave in a uniform medium. To estimate of temperature distribution, we proposed a modified method of NAH, called sectional NAH (SNAH), for using NAH in an inhomogeneous medium. In SNAH, a calculation space is discretized into a number of sections, so that the temperature in each section becomes nearly uniform, and NAH calculation is performed in each section. Calculation results using SNAH in the inhomogeneous medium well agreed with the calculation result using finite element method (FEM). By using SNAH, the temperature distribution is estimated as bellows. Firstly, sound fields, which are complex amplitude of harmonic oscillation, in three planes are measured. Next, sound field at one of measured plane is calculated from the other measured sound fields, which are calculated from other measurement planes using SNAH, with an initial sound velocity distribution. Then, difference between calculated and measured sound fields is minimized by optimizing the sound velocity distribution using multi-start downhill simplex method. In this method, a temperature distribution is obtained as the sound velocity distribution. In simulations, temperature distributions given as linear and Gaussian distributions were estimated. Validity of our proposal methods is confirmed by simulations.