Mutlispec WCS for 2 dimensional

gistfile1.txt

class MultiSpecWCS(Model):


  def __init__(self, single_wcs=[list of 1d WCSs]):
    self.single_wcs = single_wcs
  
  def __call__(self, x, y):
    assert type(x) == int # we need to make sure that the user does not expect to interpolate between two adjacent spectra
    
    return self.single_wcs[x](y)