@classmethod
def copy_dataset(cls, gdal_dataset, **kwargs):
    """Create VRT with bands and georefernce as a full copy of input GDAL Dataset

        Parameters
        ----------
        gdal_dataset : GDAL.Dataset
            input dataset
        **kwargs : dict
            arguments for VRT()

        Returns
        -------
        vrt : VRT

        """
    vrt = cls.__new__(cls)
    vrt._copy_from_dataset(gdal_dataset, **kwargs)
    return vrt
