Rasterio with NetCDF, VSIS3, and earthaccess

import numpy as np
import rasterio

def open_data():
    src = "NETCDF:earthaccess_data/20020601090000-JPL-L4_GHRSST-SSTfnd-MUR-GLOB-v02.0-fv04.1.nc:analysed_sst"
    with rasterio.open(src) as src_da:
        ma = src_da.read(1, masked=True)
        ma = ma.astype("float32", casting="unsafe")
        np.multiply(ma, src_da.scales[0], out=ma, casting="unsafe")
        np.add(ma, src_da.offsets[0], out=ma, casting="unsafe")
    return ma.filled(fill_value=np.nan)

if __name__ == "__main__":
    da = open_data()