Processing downloaded MODIS hdf4 data without using Gdal or rasterio

I tested a new workflow purely based on pyhdf and xarray, see example fileexample.nc

@daniel-eggert Is frontend able to handle the example dataset?