This project demonstrates the creation of a choropleth map to visualize population density across Indonesian provinces using R. The analysis includes importing population density data, preprocessing province names, integrating administrative boundary data from GADM, merging spatial and tabular datasets, and producing a thematic map with sf, dplyr, ggplot2, and viridis. This project showcases fundamental geospatial data processing and visualization techniques in R.