Discovering hidden geothermal signatures using non-negative matrix factorization with customized k-means clustering