A 3D-model approach has been developed to describe the electron density of the topside ionosphere and plasmasphere based on Global Navigation Satellite System (GNSS) measurements onboard low Earth orbit satellites. Electron density profiles derived from ionospheric Radio Occultation (RO) data are extrapolated to the upper ionosphere and plasmasphere based on a linear Vary-Chap function and Total Electron Content (TEC) measurements. A final update is then obtained by applying tomographic algorithms to the slant TEC measurements. Since the background specification is created with RO data, the proposed approach does not require using any external ionospheric/plasmaspheric model to adapt to the most recent data distributions. We assessed the model accuracy in 2013 and 2018 using independent TEC data, in situ electron density measurements, and ionosondes. A systematic better specification was obtained in comparison to NeQuick, with improvements around 15% in terms of electron density at 800 km, 26% at the top-most region (above 10,000 km) and 26% to 55% in terms of TEC, depending on the solar activity level. Our investigation shows that the developed model follows a known variation of electron density with respect to geographic/geomagnetic latitude, altitude, solar activity level, season, and local time, revealing the approach as a practical and useful tool for describing topside ionosphere and plasmasphere using satellite-based GNSS data.