Brain temperature, a crucial yet under-researched neurophysiological parameter, is governed by the equilibrium between cerebral oxygen metabolism and hemodynamics. Therapeutic hypothermia has been demonstrated as an effective intervention for acute brain injuries, enhancing survival rates and prognosis. The success of this treatment hinges on the precise regulation of brain temperature. However, the absence of comprehensive brain temperature monitoring methods during therapy, combined with a limited understanding of human brain heat transmission mechanisms, significantly hampers the advancement of hypothermia-based neuroprotective therapies. Leveraging the principles of bioheat transfer and MRI technology, this study conducted quantitative analyses of brain heat transfer during mild hypothermia therapy. Utilizing MRI, we reconstructed brain structures, estimated cerebral blood flow and oxygen consumption parameters, and developed a brain temperature calculation model founded on bioheat transfer theory. Employing computational cerebral hemodynamic simulation analysis, we established an intracranial arterial fluid dynamics model to predict brain temperature variations across different therapeutic hypothermia modalities. We introduce a noninvasive, spatially resolved, and optimized mathematical bio-heat model that synergizes model-predicted and MRI-derived data for brain temperature prediction and imaging. Our findings reveal that the brain temperature images generated by our model reflect distinct spatial variations across individual participants, aligning with experimentally observed temperatures.