import matplotlib.pyplot as plt
from astropy.table import Table
tbl = Table.read("../../../data/soil_spectra/exported_ICRAF_database/chemical_properties.TXT", format="csv")
plt.plot(tbl["PHH2O"], tbl["PHKCL"], ".")
plt.plot([2,11], [2,11], "k:")
plt.xlabel("pH (H2O)")
plt.ylabel("pH (KCl)")