from matplotlib import pyplot as plt
from m4opt.synphot import bandpass_from_svo
import numpy as np
from astropy import units as u
from astropy.visualization import quantity_support

filter_id = "Palomar/ZTF.r"
bandpass = bandpass_from_svo(filter_id)
wavelength = np.arange(5000, 8000) * u.angstrom
transmission = bandpass(wavelength)

with quantity_support():
    ax = plt.axes()
    ax.plot(wavelength, transmission)
    ax.set_title(filter_id)