The formation of boron esters was investigated in peat-soluble humified materials with a detailed molecular-level description of boron-organic interactions. Thousands of individually baseline separated signals were obtained from the analysis of natural organic matter of peat samples, using Fourier transform ion cyclotron resonance mass spectrometry. This technique offers unsurpassed isotope-specific mass resolution that can lead to precise molecular formula assignments by means of mathematical data analysis and visualisation techniques, such as mass defect (Kendrick) or elemental ratio (van Krevelen) plots. The analysis of potential boron binding structures within the sample of natural organic matter was described based on prior results. Herein, we describe an algorithm that can be used to effectively distinguish and filter complexes through data obtained from boron-enriched systems with highly intricate mass spectra, such as natural organic matter.