@classmethod
def from_lonlat(cls, lon, lat, add_gcps=True, **kwargs):
    """Create VRT from longitude, latitude arrays

        Create VRT with dataset without bands but with GEOLOCATION metadata and Geolocation
        object. Geolocation contains 2 2D arrays with lon/lat values given at regular pixel/line
        steps. GCPs can bee created from lon/lat arrays and added to the dataset

        Parameters
        ----------
        lon : numpy.ndarray
            array with longitudes
        lat : numpy.ndarray
            array with latitudes
        add_gcps : bool
            Create GCPs from lon/lat arrays and add to dataset
        **kwargs : dict
            arguments for VRT()

        Returns
        -------
        vrt : VRT

        """
    vrt = cls.__new__(cls)
    vrt._init_from_lonlat(lon, lat, add_gcps, **kwargs)
    return vrt
