Infinite values in geostationary bounding box crash intersection function

[mraspaud]

Infinite values in geostationary bounding box crash intersection function

[gerritholl]

Related: #393

[gerritholl]

NB: this gets triggered from Satpy if you use pad_data=False for geostationary data:

import hdf5plugin
import pathlib
import os
from satpy import Scene
from glob import glob
bd = pathlib.Path("/media/nas/x21308/MTG_test_data/2022_05_MTG_Testdata/RC0040/")
sc = Scene(filenames=[os.fspath(f) for f in bd.glob("*BODY*003[789].nc")],reader="fci_l1c_nc")
sc.load(["vis_04"], pad_data=False)
ls = sc.resample("nqeuro1km")
ls.save_datasets()