Using the remote sensing technology and the visual interpretation, this paper extracted the glacier boundaries of four phases in the border area between China and Bhutan in eastern Himalaya based on the Landsat image and the DEM image from 1990 to 2015 to investigate the relationship between the glacier boundary change and climate change. The results show that from 1990 to 2015, the shrunk area of the glacier is 191. 87 km~2,and the glacier region retreat rate on average reaches to 0. 43% per year. According to the time division, the average annual retreat rate of glaciers was increased gradually, indicating that the glacier ablation rate was gradually accelerated. The total area of glaciers in the region showed a decreasing trend,and the glacier area of different grades showed a trend of retreat. The distribution of glaciers is mainly north,north-west and north-east oriented while the main direction of retreat is east, southeast and northeast oriented. based on meteorological grid data, the change rates of air temperature and precipitation in winter and summer were retrieved by using the fitting operation on the air temperature and precipitation in different time periods. From 1990 to 2015, there was not obvious variation in precipitation while the temperature was significantly raised up. Therefore, it is believed that the main factor of glacier retreat is the rising temperature. Two specific adjacent glaciers were chosen to estimate the surface velocity using the CIAS software and it is found that the surface velocity of glacier is the fastest near the equilibrium line. The glacier flow rate was gradually slowed down to the highest altitude,and the glacier flow rate was the minimum. The glacier flow rate was gradually decreased along both sides of the glacier main stream. The surface velocity estimation shows that the chosen ice lake in the Bhutan region may speed up the melting of the glacier. There is a negative correlation between the change of glacier area and the change of temperature. There is a positive correlation between the change of ice front area and temperature change.