Hi,
You can obtain the Vth directly from de I-V curves.
PLot IDS vs VGS, at fixed small VDS(assumig that the MOSFET is operating in the triode region). The threshold voltage is estimated by linearly extrapolating back to the x-axis. For a MOSFET operating in the saturation region you can stimated the threshold voltage in the same way, also by calcuating gm=B(VGS-VTh).