High-frequency skywave propagation can be heavily impacted by D- and E-region dynamics requiring accurate global measurements to optimize performance. A standard measurement technique is to use ionosondes, but they are unable to measure below 1 MHz and are only available at a limited number of land-locked sites around the globe. In contrast, the Global Navigation Satellite System radio occultation (GNSS-RO) bottom-up method is a new approach specifically designed to generate electron density profiles in the D- and E- region ionosphere. It takes advantage of satellite constellations that currently provide over 20,000 daily measurements and global coverage. In this paper, GNSS-RO profiles were compared against ionosonde profiles at four sites covering a wide latitudinal range, and FIRI modeled profiles corresponding to the same latitude and local solar time. This comparison was completed using daytime profiles when sporadic-E (Es ) was not present. The average GNSS-RO profile is found to be a few kilometers higher in altitude than the ionosonde profiles at the minimum frequency, f min. When the ionosonde profiles are shifted so that the altitudes match at f min, they are in good agreement up to the E-region peak altitude, hmE. Below f min, the GNSS-RO profile is in good agreement with the FIRI profile, indicating that the profiles can measure the D- to E- transition region. The frequency of the E-region peak, f oE, showed general agreement between the GNSS-RO and ionosonde measurements; however, the hmE agreement was weaker and the GNSS-RO profiles tend to have an hmE in a narrow altitude range for all profiles. Virtual heights were simulated for the GNSS-RO profiles using a numerical ray tracer for direct comparison with ionosonde observations, which showed agreement for many of the virtual heights near f min, but also indicated a positive bias in the GNSS-RO virtual heights that may be due to low f oE or elevated hmE estimates. For a quiet ionosphere, the shifted GNSS-RO electron density profiles were a good match for both measured ionosonde profiles and modeled FIRI profiles and the method is capable of providing global coverage of the D- and E-regions. Future work will require more data for seasonal and morning–afternoon comparisons as well as comparisons for the disturbed ionosphere when the sporadic-E layer is present.