  We present an analysis of late-time Hubble Space Telescope Wide Field Camera
3 and Wide Field Planetary Camera 2 observations of the site of the Type Ic SN
2007gr in NGC 1058. The SN is barely recovered in the late-time WFPC2
observations, while a possible detection in the later WFC3 data is debatable.
These observations were used to conduct a multiwavelength study of the
surrounding stellar population. We fit spatial profiles to a nearby bright
source that was previously proposed to be a host cluster. We find that, rather
than being an extended cluster, it is consistent with a single point-like
object. Fitting stellar models to the observed spectral energy distribution of
this source, we conclude it is A1-A3 Yellow Supergiant, possibly corresponding
to a star with $M_{ZAMS} = 40M_{\odot}$. SN 2007gr is situated in a massive
star association, with diameter of $\approx 300\,\mathrm{pc}$. We present a
Bayesian scheme to determine the properties of the surrounding massive star
population, in conjunction with the Padova isochrones. We find that the stellar
population, as observed in either the WFC3 and WFPC2 observations, can be well
fit by two age distributions with mean ages: ~6.3 Myr and ~50 Myr. The stellar
population is clearly dominated by the younger age solution (by factors of 3.5
and 5.7 from the WFPC2 and WFC3 observations, respectively), which corresponds
to the lifetime of a star with $M_{ZAMS} \sim 30M_{\odot}$. This is strong
evidence in favour of the hypothesis that SN 2007gr arose from a massive
progenitor star, possibly capable of becoming a Wolf-Rayet star.
