nb2. Query a catalogue

import astropy
import astropy.coordinates
import pandas as pd
pd.set_option('display.max_columns', None)
import numpy as np
import pyvo
import matplotlib.pyplot as plt

service = pyvo.dal.TAPService("https://almascience.eso.org/tap")      # for the EU ALMA TAP service

# service = pyvo.dal.TAPService("https://almascience.nao.ac.jp/tap")  # for the EA ALMA TAP service
# service = pyvo.dal.TAPService("https://almascience.nrao.edu/tap")   # for the NA ALMA TAP service

def query_coordinate_list(service, ra_list, dec_list, radius):
    """Runs a single query for an entire list of coordinates. Query only 100 or less coordinates in one go.
    
       service        pyvo TAPService instance           
       ra_list        List of Right Ascension values in decimal degrees
       dec            List of Declination values in decimal degrees
       radius         Search radius in decimal degrees
           
       returns        pandas table
    """
    
    conditionlist = [f"INTERSECTS(CIRCLE('ICRS',{c[0]},{c[1]},{radius}),s_region)=1" for c in zip(ra_list, dec_list)]
    
    query = f"""
            SELECT *
            FROM ivoa.obscore
            WHERE {" OR ".join(conditionlist)}"""
    
    return service.search(query).to_table().to_pandas()

def query_resolved_source_name_list(service, source_name_list, radius):
    """Resolves the list of source names with Sesame and runs a coordinate search on them
    
       service           pyvo TAPService instance           
       source_name_list  list with source names
           
       returns           pandas table
    """
    
    coordinates = []
    for source in source_name_list:
        try:
            coordinates.append(astropy.coordinates.SkyCoord.from_name(source))
        except astropy.coordinates.name_resolve.NameResolveError:
            print(f"Warning: source '{source}' could not be resolved by Sesame.")
            
    return query_coordinate_list(service, [c.ra.degree for c in coordinates], [c.dec.degree for c in coordinates], radius)
    

def query_alma_source_name_list(service, source_name_list):
    """queries for list of source names as given by the PI in the ALMA proposals
    
       service            pyvo TAPService instance           
       source_name_lis    complete source name 
           
       returns            pandas table
    """
    
    query = f"""  
            SELECT *  
            FROM ivoa.obscore  
            WHERE target_name in ('{"','".join(source_name_list)}') """
    
    return service.search(query).to_table().to_pandas()

ra_list = highzqso_catalog['RAJ2000']
dec_list = highzqso_catalog['DEJ2000']

result = pd.concat([query_coordinate_list(service, ra_list[i:i+100], dec_list[i:i+100], 0.01) for i in range(0, len(ra_list), 100)])

query = """SELECT DISTINCT(target_name) from ivoa.obscore where target_name like 'NGC____'"""
source_name_list = list(service.search(query).to_table().to_pandas()['(target_name)'].values)

query = """SELECT DISTINCT(target_name) from ivoa.obscore where target_name like 'NGC____'"""
source_name_list = list(service.search(query).to_table().to_pandas()['(target_name)'].values)