pyvo

import pyvo import os # 1. Search the VO registry for services # Example: Find services related to galaxies and spectroscopy print("Searching the VO registry...") results = pyvo.registry.search(keywords=['galaxy', 'spectroscopy']) print(f"Found {len(results)} services matching 'galaxy' and 'spectroscopy'.") # 2. Access a TAP service (Table Access Protocol) # For a stable quickstart, we'll use a widely known public TAP service. # Note: Service availability can change. This is an example. tap_service_url = os.environ.get('PYVO_TAP_SERVICE_URL', 'https://tap.voservices.net/tap') tap_service = pyvo.dal.TAPService(tap_service_url) print(f"\nAttempting to connect to TAP service: {tap_service.url}") # 3. Execute an ADQL (Astronomical Data Query Language) query # Example: Query the ivoa.ObsCore table for Gaia data query = "SELECT TOP 5 ra, dec, dataproduct_type FROM ivoa.ObsCore WHERE obs_collection LIKE '%Gaia%'" try: tap_results = tap_service.search(query) print(f"\nTAP query returned {len(tap_results)} results.") # Results are returned as an Astropy Table object (pyvo.dal.TableSet) print("\nFirst 5 results (Astropy Table):") print(tap_results.to_table()) except Exception as e: print(f"\nError during TAP query: {e}") print("The TAP service might be temporarily unavailable or the query invalid.")