@classmethod
def from_gdal_dataset(cls, gdal_dataset, **kwargs):
    """Create VRT from GDAL Dataset with the same size/georeference but wihout bands.

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

        Returns
        -------
        vrt : VRT

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