Lazy in-memory reprojection to a user-specified Coordinate Reference System (CRS)

reproject_when_opening.py

import xarray as xr
import rasterio
from rasterio.vrt import WarpedVRT

fname = '/data/peatlands_test_data/MOD11A1/LST_Day_1km/h17v03/MOD11A1_LST_Day_1km_h17v03_2019-08.tif'

# Lazy in-memory warping from Sinusoidal to WGS84 lat/lon
with rasterio.open(fname) as src:
    # print(src.profile)
    with WarpedVRT(src, crs='epsg:4326') as vrt:
        # print(vrt.profile)
        with xr.open_rasterio(vrt) as dataset_wgs84:
            print(dataset_wgs84)

# Lazy in-memory warping from Sinusoidal to BNG
# BNG definition
# from https://spatialreference.org/ref/epsg/osgb-1936-british-national-grid/
# User should pass those parameters?
bng_crs = (f' +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 '
           f' +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 '
           f' +units=m +no_defs ')

with rasterio.open(fname) as src:
    # print(src.profile)
    with WarpedVRT(src, crs=bng_crs) as vrt:
        # print(vrt.profile)
        with xr.open_rasterio(vrt) as dataset_bng:
            print(dataset_bng)