Passing scpos to GbmHealPix should generate geo_loc the same way it generates geo_radius

[eburnsastro]

Example code below. This generates a GbmHealPix from a Gaussian, adds a systematic, then attempts to remove Earth. After running the below script, 'example_map_stat_only.geo_radius' gives 67.5 deg, demonstrating it is using the x,y,z information in scpos to generate the radius. I believe this can similarly be used to generate geo_loc, with the (ra_dec) for the geo center being the spherical representation calculated by converting the negative of the scpos

from gdt.missions.fermi.gbm.localization import GbmHealPix
example_RA = 178.017
example_Dec = -18.117
example_Err = 7.88


# generation of Gaussian skymap using the data tools
example_map_stat_only = GbmHealPix.from_gaussian(example_RA, example_Dec, 
                            example_Err, scpos=array([1656.69361922, 6634.3140278 ,  875.01112923]))
skyplot = EquatorialPlot()
skyplot.add_localization(example_map_stat_only, detectors=None)
plt.show()

# as the data tools have convolution built in to add a systematic
example_map_full_errors = example_map_stat_only.convolve(GBM_GND_Model)
skyplot = EquatorialPlot()
skyplot.add_localization(example_map_full_errors, detectors=None)
plt.show()

# example_map_full_errors_no_earth = example_map_full_errors.remove_earth()
example_map_full_errors_no_earth = GbmHealPix.remove_earth(example_map_full_errors)
skyplot = EquatorialPlot()
skyplot.add_localization(example_map_full_errors_no_earth, detectors=None)
plt.show()

[eburnsastro]

It looks like this may already be implemented in the GBM data tools, but not gdt-fermi

[eburnsastro]

From talking with Adam, it is not calculating georadius, it is just using a default value