def predict(model, X, width, height):
    """Prediction with a fitted model.

    Parameters
    ----------
    model : classifier
        Sklearn classifier (fitted).
    X : array
        Input data.
    width : int
        Raster width.
    height : int
        Raster height.

    Returns
    -------
    proba : 2d array
        RF probabilistic output of shape (height, width).
    """
    proba = model.predict_proba(X)[:, 0]
    proba = proba.reshape((height, width))
    return proba