Mapping large wildfires (LW) is essential for environmental applications and enhances the understanding of the dynamics of affected areas. Remote sensing techniques supported by machine learning and time series have been increasingly used in studies addressing this issue and have shown potential for this type of analysis. The main aim of this article is to develop a methodology for mapping LW in northwestern Portugal using a machine learning algorithm and time series from Landsat images. For the burnt area classification, we initially used the Fourier harmonic model to define outliers in the time series that represented pixels of possible burnt areas and, then, we applied the random forest classifier for the LW classification. The results indicate that the harmonic analysis provided estimates with the actual observed values of the NBR index; thus, the pixels classified by random forest were only those that were masked, collaborated in the processing, and reduced possible spectral confusion between targets with similar behaviour. The burnt area maps revealed that ~23.5% of the territory was burnt at least once from 2001 to 2020. The temporal variability of the burnt area indicated that, on average, 6.504 hectares were affected by LW within the 20 years. The annual burnt area varied over the years, with the minimum annual area detected in 2014 (679.5 hectares) and the maximum mapped area detected in 2005 (73,025.1 hectares). We concluded that the process of defining the mask with the outliers considerably reduced the universe of pixels to be classified within each image, which leaves the training of the classifier focused on separating the set of pixels into two groups with very similar spectral characteristics, thus contributing so that the separation of groups with similar spectral behaviour was performed automatically and without great sampling effort. The method showed satisfactory accuracy results with little omission for burnt areas.