def get_ndv_fn(fn): ds = gdal.Open(fn, gdal.GA_ReadOnly) return get_ndv_ds(ds) #Want to modify to handle multi-band images and return list of ndv