code: def load_tiff(filename, crs=None, apply_transform=False, nan_nodata=False, **kwargs):
      
      try:
          import xarray as xr
      except:
          raise ImportError('Loading tiffs requires xarray to be installed')
  
      with warnings.catch_warnings():
          warnings.filterwarnings('ignore')
          da = xr.open_rasterio(filename)
      return from_xarray(da, crs, apply_transform, nan_nodata, **kwargs)
docs: Returns an RGB or Image element loaded from a geotiff file.
  
      The data is loaded using xarray and rasterio. If a crs attribute
      is present on the loaded data it will attempt to decode it into
      a cartopy projection otherwise it will default to a non-geographic
      HoloViews element.
  
      Parameters
      ----------
      filename: string
        Filename pointing to geotiff file to load
      crs: Cartopy CRS or EPSG string (optional)
        Overrides CRS inferred from the data
      apply_transform: boolean
        Whether to apply affine transform if defined on the data
      nan_nodata: boolean
        If data contains nodata values convert them to NaNs
      **kwargs:
        Keyword arguments passed to the HoloViews/GeoViews element
  
      Returns
      -------
      element: Image/RGB/QuadMesh element