Based on the DESI Legacy Imaging Surveys released data and available spectroscopic redshifts, we identify 1.58 million clusters of galaxies by searching for the overdensity of the stellar mass distribution of galaxies within redshift slices around preselected massive galaxies, with 877,806 clusters being found for the first time. The identified clusters have an equivalent mass of M_500_>=0.47x10^14^M_{sun}_ with an uncertainty of 0.2dex. The redshift distribution of clusters extends to z~1.5, and 338,841 clusters have spectroscopic redshifts. Our cluster sample includes most of the rich optical clusters in previous catalogs, more than 95% of the massive Sunyaev-Zeldovich clusters, and 90% of the ROSAT and eROSITA X-ray clusters. From the light distributions of the member galaxies, we derive the dynamical state parameters for 27,685 rich clusters and find no significant evolution of the dynamical state with redshift. We find that the stellar mass of the brightest cluster galaxies grows by a factor of 2 since z=1.