We calculate full 3D pressure fields for inhomogeneous nanoscale systems using molecular dynamics simulation data. The fields represent systems with increasing level of complexity, ranging from semivesicles and vesicles to membranes characterized by coexistence of two phases, including also a protein-membrane complex. We show that the 3D pressure field is distinctly different for curved and planar bilayers, the pressure field depends strongly on the phase of the membrane, and that an integral protein modulates the tension and elastic properties of the membrane.