and with some algebra (that you are free to check), I get
rate = POINT.EST/DOF
shape = DOF
which you can then plug into base R functions like dgamma().
The standard deviation would be given by sqrt() of the variance above.
If I now understand your previous question, throughout ctmm, I try to use CI formulas that are closest to the true sampling distribution. So for area estimates, I use gamma CIs, which are exact for the IID isotropic case and respect the lower bound of zero area.