Solving large tomographic linear systems: size reduction and error estimation