Since you have the structures the best option, in my opinion, is Pymol.
Open Pymol. Load up your protein of choice. Download color_h.py which is a script from the University of Osaka that colours the residues according to Eisenberg's scale of hydrophobicity. Load this into Pymol by
File->Run->PATH/TO/color_h.py. Then in Pymol run the following commands
Note that the script can be edited to use whatever scales you like. @mimat points you toward protscale and I would suggest the same. The Eisenberg scale is a rather old consensus scale. More up to date scales might be more appropriate in your case.
To start quantifying this, Pymol has tools for calculating solvent accessibility for example.
Additionally PDBSum and PDBe-PISA can estimate solvent accesible surface area.
The quickest and easiest way to generate electrostatics is with the inbuilt vacuum electrostatics tools. In Pymol, select the
action button (A)->generate->vacuum electrostatics->protein contact potential.
Since the Pymol inbuilt electrostatics make a lot of assumptions, for publication it is best to get a much more accurately quantified surface map. However, ABPS is a bit more involved. Read the wiki if this sounds more like what you need.