def METHOD_NAME(OBJ_RA, OBJ_DEC, binVect_RA, binVect_DEC, N_OBJ_th=2):
    RA_DEC_meshgrid = np.meshgrid(binVect_RA, binVect_DEC)
    OBJ_dNdRAdDEC = np.histogram2d(x=OBJ_DEC, y=OBJ_RA,
                                   bins=[binVect_DEC, binVect_RA])[0]
    med_OBJ_dNdRAdDEC = np.median(OBJ_dNdRAdDEC)
    skyArea_meshgrid = AreaFunc(RA_DEC_meshgrid[1][0:-1, 0:-1] * np.pi / 180.,
                                RA_DEC_meshgrid[1][1:, 0: -1] * np.pi / 180.,
                                RA_DEC_meshgrid[0][0:-1, 0: -1] * np.pi / 180.,
                                RA_DEC_meshgrid[0][0:-1, 1:] * np.pi / 180.)
    area_OK = (OBJ_dNdRAdDEC) >= N_OBJ_th
    res_area = np.sum(skyArea_meshgrid[area_OK])
    return res_area, med_OBJ_dNdRAdDEC