I came across Johnson's system of distributions. Here, I found the following figure:
The author says:
There is a unique Johnson distribution corresponding to each feasible combination of $\beta_1$ and $\beta_2$[, i.e. skewness and kurtosis].
My question is as follows. How can I draw a sample from a distribution that corresponds to prespecified values of $\beta_1$ and $\beta_2$? Is there any python or R package that can do that?