def intensity_localization_fv(data3dr, voxelsize_mm, seeds=None, unique_cls
    =None):
    """
    Use organ_localizator features plus intensity features

    :param data3dr:
    :param voxelsize_mm:
    :param seeds:
    :param unique_cls:
    :return:
    """
    import scipy
    import numpy as np
    import os.path as op
    try:
        from lisa import organ_localizator
    except:
        import organ_localizator
    f0 = scipy.ndimage.filters.gaussian_filter(data3dr, sigma=0.5).reshape(
        -1, 1)
    f1 = scipy.ndimage.filters.gaussian_filter(data3dr, sigma=3).reshape(-1, 1)
    import organ_localizator
    fvall = organ_localizator.localization_fv(data3dr, voxelsize_mm)
    fv = np.concatenate([f0, f1, fvall], 1)
    if seeds is not None:
        sd = seeds.reshape(-1, 1)
        selection = np.in1d(sd, unique_cls)
        fv = fv[selection]
        sd = sd[selection]
        return fv, sd
    return fv
