An identification method of urban functional areas based on mixing degree of functions and integrated learning includes the following steps: 1) performing data acquisition and preprocessing; 2) constructing 10 indicator features of an urban functional area identification system; 3) structuring the indicator features: acquiring, by a spatial statistical tool, the 10 indicator features corresponding to each parcel; 4) constructing an independent variable dataset; 5) labeling response variables; 6) dividing a training dataset into a plurality of training subsets according to the mixing degree of functions; 7) training a Stacking-based integrated learning model; and 8) joining an attribute in one table to another table, so as to complete the identification of the urban functional areas on each parcel. The identification method divides the training dataset by grading the mixing degree of functions, and makes predictions based on the prediction dataset with corresponding mixing degree of functions.