We aim at identifying the clusters members by deriving membership probabilities for the sources within 1 degree of the clusters center, going further away than equivalent previous studies. We measure accurate proper motions and multi-wavelength (optical and near-infrared) photometry using ground based archival images of the cluster. We use these measurements to compute membership probabilities. The list of candidate members from Barrado+2001 is used as training set to identify the clusters locus in a multi-dimensional space made of proper motions, luminosities and colors. The final catalog includes 338892 sources with multi-wavelength photometry. Approximately half (194452) were detected at more than two epochs and we measured their proper motion and used it to derive membership probability. A total of 4349 candidate members with membership probabilities greater than 50% are found in this sample in the luminosity range between 10 and 22mag. The slow proper motion of the cluster and the overlap of its sequence with the field and background sequences in almost all color-magnitude and color-color diagrams complicate the analysis and the contamination level is expected to be significant. Our study nevertheless provides a coherent and quantitative membership analysis of Messier 35 based on a large fraction of the best ground-based data sets obtained over the past 18 years. As such, it represents a valuable input for follow-up studies using in particular the Kepler K2 photometric time series.