Focused ultrasound (FUS) is a non-invasive and non-ionizing technique which deploys ultrasound waves to induce bio-effects. When paired with acoustically active particles such as microbubbles (MBs), it can open the blood brain barrier (BBB) to facilitate drug delivery otherwise inhibited due to the presence of BBB. One of the parameters that affects the FUS beam propagation is the beam incidence angle on the skull. Prior work by our group has shown that, as incidence angles deviate from 90°, FUS focal pressures attenuate and result in a smaller BBB opening volume. The incidence angles calculated in our prior studies were in 2D and used skull information from CT. The study presented herein develops methods to calculate incidence angle in 3D in non-human primate (NHP) skull fragments using harmonic ultrasound imaging without using ionizing radiation. Our results show that ultrasound harmonic imaging is capable of accurately depicting features such as sutures and eye-sockets of the skull. Furthermore, we were able to reproduce previously reported relationships between the incidence angle and FUS beam attenuation. We also show feasibility of performing ultrasound harmonic imaging in in-vivo non-human primates. The all-ultrasound method presented herein combined with our neuronavigation system stands to increase more widespread adoption of FUS and render it accessible by eliminating the need for CT cranial mapping.