The propagation of large amplitude ion-acoustic solitary waves (IASWs) in a fully relativistic plasma consisting of cold ions and ultra-relativistic hot electrons and positrons is investigated using the Sagdeev pseudopotential method in a relativistic hydrodynamics model. The effects of streaming speed of the plasma fluid, thermal energy, positron density, and positron temperature on large amplitude IASWs are studied by analysis of the pseudopotential structure. It is found that in regions in which the streaming speed of the plasma fluid is larger than that of the solitary wave, by increasing the streaming speed of the plasma fluid, the depth and width of the potential well increase, resulting in narrower solitons with larger amplitude. This behavior is opposite to the case where the streaming speed of the plasma fluid is less than that of the solitary wave. On the other hand, an increase in the thermal energy results in wider solitons with smaller amplitude, because the depth and width of the potential well decrease in that case. Additionally, the maximum soliton amplitude increases and the width becomes narrower as a result of an increase in positron density. It is shown that varying the positron temperature does not have a considerable effect on the width and amplitude of IASWs. The existence of stationary soliton-like arbitary amplitude waves is also predicted in fully relativistic electron-positron-ion (EPI) plasmas. The effects of streaming speed of the plasma fluid, thermal energy, positron density, and positron temperature on these kinds of solitons are the same for large amplitude IASWs.