The Hubble Space Telescope (HST) has been operational for almost 30 years and throughout that time it has been bombarded by high energy charged particles colloquially referred to as cosmic rays. In this paper, we present a comprehensive study of more than 1.2 billion cosmic rays observed with HST using a custom written python package, \texttt{HSTcosmicrays}, that is available to the astronomical community. We analyzed $75,908$ dark calibration files taken as part of routine calibration programs for five different CCD imagers with operational coverage of Solar Cycle 23 and 24. We observe the expected modulation of galactic cosmic rays by solar activity. For the three imagers with the largest non-uniformity in thickness, we independently confirm the overall structure produced by fringing analyses by analyzing cosmic ray strikes across the detector field of view. We analyze STIS/CCD observations taken as HST crosses over the South Atlantic Anomaly and find a peak cosmic ray flux of $\sim1100$ $CR/s/cm^2$. We find strong evidence for two spatially confined regions over North America and Australia that exhibit increased cosmic ray fluxes at the $5\sigma$ level.