Background: Malaria is a major public health burden in Southeastern Bangladesh, particularly in the Chittagong Hill Tracts region. Malaria is endemic in 13 districts of Bangladesh and the highest prevalence occurs in Khagrachari (15.47%). Methods. A risk map was developed and geographic risk factors identified using a Bayesian approach. The Bayesian geostatistical model was developed from previously identified individual and environmental covariates (p < 0.2; age, different forest types, elevation and economic status) for malaria prevalence using WinBUGS 1.4. Spatial correlation was estimated within a Bayesian framework based on a geostatistical model. The infection status (positives and negatives) was modeled using a Bernoulli distribution. Maps of the posterior distributions of predicted prevalence were developed in geographic information system (GIS). Results: Predicted high prevalence areas were located along the north-eastern areas, and central part of the study area. Low to moderate prevalence areas were predicted in the southwestern, southeastern and central regions. Individual age and nearness to fragmented forest were associated with malaria prevalence after adjusting the spatial auto-correlation. Conclusion: A Bayesian analytical approach using multiple enabling technologies (geographic information systems, global positioning systems, and remote sensing) provide a strategy to characterize spatial heterogeneity in malaria risk at a fine scale. Even in the most hyper endemic region of Bangladesh there is substantial spatial heterogeneity in risk. Areas that are predicted to be at high risk, based on the environment but that have not been reached by surveys are identified.