Sorry, Candice, I misunderstood your question.
geomorph works with landmark shape data. One possibility would be to do a subsampling of your outline data to produce semilandmarks, using the digit.curves() function of geomorph, and possibly complement these with some fixed landmarks, which could then be used in a landmark-based shape analysis. In that case, you would need to superimpose your data using gpagen() with the curves argument, and then do the PCA using plotTangentSpace. But this is an analytically completely different solution compared to outline analysis using EFA.
Alternatively, the PCA and scatter-plotting part of plotTangentSpace are both based on the prcomp() base function of R, so producing plots with your EFA-treated data should be straightforward by setting uniform plot() arguments in R.
Hope this helps,
Antigoni