The Hubble Space Telescope (HST) has been operational for over 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, 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. We model the observed energy-loss distributions to derive an estimate of 534 ± 117 MeV for the kinetic energy of the typical cosmic ray impacting HST. For the three imagers with the largest nonuniformity 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 particle flux of ∼1100 particle s−1 cm−2. We find strong evidence for two spatially confined regions over North America and Australia that exhibit increased cosmic ray particle fluxes at the 5σ level.