docs: def join(rasters):
      """
      This method takes a list of rasters and returns a raster that is constructed of all of them
      """
  
      raster = rasters[0]  # using the first raster to understand what is the type of data we have
      mask_band = None
      nodata = None
      with raster._raster_opener(raster.source_file) as r:
          nodata = r.nodata
          mask_flags = r.mask_flag_enums
      per_dataset_mask = all([rasterio.enums.MaskFlags.per_dataset in flags for flags in mask_flags])
      if per_dataset_mask and nodata is None:
          mask_band = 0
      return GeoRaster2.from_rasters(rasters, relative_to_vrt=False, nodata=nodata, mask_band=mask_band)