Implement "Tuned Geodetic Information" procedure

[marco-buttu]

The EarthLocation class of Astropy library gives information about observatories. That's for example the information about SRT:

>>> from astropy.coordinates import EarthLocation
>>> srt = EarthLocation.of_site("SRT")
>>> srt.info.name
'Sardinia Radio Telescope'
>>> srt.height
<Quantity 671.6665 m>

The height of SRT according to Astropy is 671.7 meters. We realized that this height is not the same saved on the FITS files generated from DISCOS, as you can see here:

>>> from astropy.io import fits
>>> hdul = fits.open('tests/data/xyz.fits')
>>> hdul[0].header["ANTENNA"]
'SRT'
>>> hdul[0].header["SiteHeight"]
650.0

I would like to write a procedure that, given an observatory name, returns the average geodetic information of that observatory, computed taking in account the values from Astropy and from a FITS file of that observatory. The returned information contains: observatory name, latitude, longitude, and height.

[marco-buttu]

Fixed in PR #4