Estimating the size of key populations, including female sex workers (FSW) and men who have sex with men (MSM), can inform planning and resource allocation for HIV programs at local and national levels. In geographic areas where direct population size estimates (PSEs) for key populations have not been collected, small area estimation (SAE) can help fill in gaps using supplemental data sources known as auxiliary data. However, routinely collected program data have not historically been used as auxiliary data to generate subnational estimates for key populations, including in Namibia. To systematically generate regional size estimates for FSW and MSM in Namibia, we used a consensus-informed estimation approach with local stakeholders that included the integration of routinely collected HIV program data provided by key populations' HIV service providers. We used quarterly program data reported by key population implementing partners, including counts of the number of individuals accessing HIV services over time, to weight existing PSEs collected through bio-behavioral surveys using a Bayesian triangulation approach. SAEs were generated through simple imputation, stratified imputation, and multivariable Poisson regression models. We selected final estimates using an iterative qualitative ranking process with local key population implementing partners. Extrapolated national estimates for FSW ranged from 4777 to 13,148 across Namibia, comprising 1.5% to 3.6% of female individuals aged between 15 and 49 years. For MSM, estimates ranged from 4611 to 10,171, comprising 0.7% to 1.5% of male individuals aged between 15 and 49 years. After the inclusion of program data as priors, the estimated proportion of FSW derived from simple imputation increased from 1.9% to 2.8%, and the proportion of MSM decreased from 1.5% to 0.75%. When stratified imputation was implemented using HIV prevalence to inform strata, the inclusion of program data increased the proportion of FSW from 2.6% to 4.0% in regions with high prevalence and decreased the proportion from 1.4% to 1.2% in regions with low prevalence. When population density was used to inform strata, the inclusion of program data also increased the proportion of FSW in high-density regions (from 1.1% to 3.4%) and decreased the proportion of MSM in all regions. Using SAE approaches, we combined epidemiologic and program data to generate subnational size estimates for key populations in Namibia. Overall, estimates were highly sensitive to the inclusion of program data. Program data represent a supplemental source of information that can be used to align PSEs with real-world HIV programs, particularly in regions where population-based data collection methods are challenging to implement. Future work is needed to determine how best to include and validate program data in target settings and in key population size estimation studies, ultimately bridging research with practice to support a more comprehensive HIV response.