A metrology and data analysis protocol is described for high throughput determination of thermochromic metal-insulator phase diagrams for lightly substituted VO2 thin films. The technique exploits the abrupt change in near infrared optical properties, measured in reflection, as an indicator of the temperature- or impurity-driven metal-insulator transition. Transition metal impurities were introduced in a complementary combinatorial synthesis process for producing thin film libraries with the general composition space V(1-x-y)M(x)M'(y)O2, with M and M' being transition metals and x and y varying continuously across the library. The measurement apparatus acquires reflectance spectra in the visible or near infrared at arbitrarily many library locations, each with a unique film composition, at temperatures of 1 °C-85 °C. Data collection is rapid and automated; the measurement protocol is computer controlled to automate the collection of thousands of reflectance spectra, representing hundreds of film compositions at tens of different temperatures. A straightforward analysis algorithm is implemented to extract key information from the thousands of spectra such as near infrared thermochromic transition temperatures and regions of no thermochromic transition; similarly, reflectance to the visible spectrum generates key information for materials selection of smart window materials. The thermochromic transition for 160 unique compositions on a thin film library with the general formula V(1-x-y)M(x)M'(y)O2 can be measured and described in a single 20 h experiment. The resulting impurity composition-temperature phase diagrams will contribute to the understanding of metal-insulator transitions in doped VO2 systems and to the development of thermochromic smart windows.