The data is from:

https://simplemaps.com/data/world-cities

We're proud to offer a simple, accurate and up-to-date database of the world's cities and towns. We've built it from the ground up using authoritative sources such as the NGIA, US Geological Survey, US Census Bureau, and NASA.

Our database is:

• Up-to-date: It was last refreshed on May 11, 2025.
• Comprehensive: Over 4 million unique cities and towns from every country in the world (about 48 thousand in basic database).
• Accurate: Cleaned and aggregated from official sources. Includes latitude and longitude coordinates.
• Simple: A single CSV file, concise field names, only one entry per city.

It is estimated that more than 8 billion people live on Earth and the population is likely to hit more than 9 billion by 2050. Approximately 55 percent of Earth’s human population currently live in areas classified as urban. That number is expected to grow by 2050 to 68 percent, according to the United Nations (UN).The largest cities in the world include Tōkyō, Japan; New Delhi, India; Shanghai, China; México City, Mexico; and São Paulo, Brazil. Each of these cities classifies as a megacity, a city with more than 10 million people. The UN estimates the world will have 43 megacities by 2030.Most cities' populations are growing as people move in for greater economic, educational, and healthcare opportunities. But not all cities are expanding. Those cities whose populations are declining may be experiencing declining fertility rates (the number of births is lower than the number of deaths), shrinking economies, emigration, or have experienced a natural disaster that resulted in fatalities or forced people to leave the region.This Global Cities map layer contains data published in 2018 by the Population Division of the United Nations Department of Economic and Social Affairs (UN DESA). It shows urban agglomerations. The UN DESA defines an urban agglomeration as a continuous area where population is classified at urban levels (by the country in which the city resides) regardless of what local government systems manage the area. Since not all places record data the same way, some populations may be calculated using the city population as defined by its boundary and the metropolitan area. If a reliable estimate for the urban agglomeration was unable to be determined, the population of the city or metropolitan area is used.Data Citation: United Nations Department of Economic and Social Affairs. World Urbanization Prospects: The 2018 Revision. Statistical Papers - United Nations (ser. A), Population and Vital Statistics Report, 2019, https://doi.org/10.18356/b9e995fe-en.

This world cities layer presents the locations of many cities of the world, both major cities and many provincial capitals.Population estimates are provided for those cities listed in open source data from the United Nations and US Census.

This city boundary shapefile was extracted from Esri Data and Maps for ArcGIS 2014 - U.S. Populated Place Areas. This shapefile can be joined to 500 Cities city-level Data (GIS Friendly Format) in a geographic information system (GIS) to make city-level maps.

This map features boundaries and places for the World, including countries, 1st order administrative areas, and cities. The map layers are delivered as features, which you can click on for attribute information or re-symbolize as you choose.

This global accessibility map enumerates land-based travel time to the nearest densely-populated area for all areas between 85 degrees north and 60 degrees south for a nominal year 2015. Densely-populated areas are defined as contiguous areas with 1,500 or more inhabitants per square kilometre or a majority of built-up land cover types coincident with a population centre of at least 50,000 inhabitants. This map was produced through a collaboration between MAP (University of Oxford), Google, the European Union Joint Research Centre (JRC), and the University of Twente, Netherlands.The underlying datasets used to produce the map include roads (comprising the first ever global-scale use of Open Street Map and Google roads datasets), railways, rivers, lakes, oceans, topographic conditions (slope and elevation), landcover types, and national borders. These datasets were each allocated a speed or speeds of travel in terms of time to cross each pixel of that type. The datasets were then combined to produce a "friction surface"; a map where every pixel is allocated a nominal overall speed of travel based on the types occurring within that pixel. Least-cost-path algorithms (running in Google Earth Engine and, for high-latitude areas, in R) were used in conjunction with this friction surface to calculate the time of travel from all locations to the nearest (in time) city. The cities dataset used is the high-density-cover product created by the Global Human Settlement Project. Each pixel in the resultant accessibility map thus represents the modelled shortest time from that location to a city. Authors: D.J. Weiss, A. Nelson, H.S. Gibson, W. Temperley, S. Peedell, A. Lieber, M. Hancher, E. Poyart, S. Belchior, N. Fullman, B. Mappin, U. Dalrymple, J. Rozier, T.C.D. Lucas, R.E. Howes, L.S. Tusting, S.Y. Kang, E. Cameron, D. Bisanzio, K.E. Battle, S. Bhatt, and P.W. Gething. A global map of travel time to cities to assess inequalities in accessibility in 2015. (2018). Nature. doi:10.1038/nature25181

Processing notes: Data were processed from numerous sources including OpenStreetMap, Google Maps, Land Cover mapping, and others, to generate a global friction surface of average land-based travel speed. This accessibility surface was then derived from that friction surface via a least-cost-path algorithm finding at each location the closest point from global databases of population centres and densely-populated areas. Please see the associated publication for full details of the processing.

Source: https://map.ox.ac.uk/research-project/accessibility_to_cities/

This dataset includes a total of 635 annotated image patches from historical city maps. It is designed for the semantic segmentation of the maps into 5 semantic classes (building blocks, non-built, water, road network, background frame). 330 patches are taken from maps of the city of Paris, while the 305 others are taken from a balanced corpus of city maps from 90 countries all around the world.

Please read the detailed informations about data collection methodology, associated metadata and annotation ontology in README.md hereunder :

The "Major Cities" layer is derived from the "World Cities" dataset provided by ArcGIS Data and Maps group as part of the global data layers made available for public use.

"Major cities" layer specifically contains National and Provincial capitals that have the highest population within their respective country. Cities were filtered based on the STATUS (“National capital”, “National and provincial capital”, “Provincial capital”, “National capital and provincial capital enclave”, and “Other”). Majority of these cities within larger countries have been filtered at the highest levels of POP_CLASS (“5,000,000 and greater” and “1,000,000 to 4,999,999”). However, China for example, was filtered with cities over 11 million people due to many highly populated cities. Population approximations are sourced from US Census and UN Data.

Disclaimer: The designations employed and the presentation of material at this site do not imply the expression of any opinion whatsoever on the part of the Secretariat of the United Nations concerning the legal status of any country, territory, city or area or of its authorities, or concerning the delimitation of its frontiers or boundaries.

Data publication: 2021-03-12

Citation:

Credits: ESRI, CIA World Factbook, GMI, NIMA, UN Data, UN Habitat, US Census Bureau

Contact points:

Resource Contact: ESRI - ArcGIS Data and Maps

Metadata Contact: Justeen De Ocampo

Resource constraints:

Creative Commons Attribution-NonCommercial-ShareAlike 3.0 IGO (CC BY-NC- SA 3.0 IGO)

Online resources:

World Cities layer from ArcGIS Data & Maps

ArcGIS Data and Maps group background and available datasets.

This National Geographic Style Map (World Edition) web map provides a reference map for the world that includes administrative boundaries, cities, protected areas, highways, roads, railways, water features, buildings, and landmarks, overlaid on shaded relief and a colorized physical ecosystems base for added context to conservation and biodiversity topics. Alignment of boundaries is a presentation of the feature provided by our data vendors and does not imply endorsement by Esri, National Geographic or any governing authority.This basemap, included in the ArcGIS Living Atlas of the World, uses the National Geographic Style vector tile layer and the National Geographic Style Base and World Hillshade raster tile layers.The vector tile layer in this web map is built using the same data sources used for other Esri Vector Basemaps. For details on data sources contributed by the GIS community, view the map of Community Maps Basemap Contributors. Esri Vector Basemaps are updated monthly.Use this MapThis map is designed to be used as a basemap for overlaying other layers of information or as a stand-alone reference map. You can add layers to this web map and save as your own map. If you like, you can add this web map to a custom basemap gallery for others in your organization to use in creating web maps. If you would like to add this map as a layer in other maps you are creating, you may use the tile layers referenced in this map.

The Human Geography Map (World Edition) web map provides a detailed vector basemap with a monochromatic style and content adjusted to support Human Geography information. Where possible, the map content has been adjusted so that it observes WCAG contrast criteria.This basemap, included in the ArcGIS Living Atlas of the World, uses 3 vector tile layers:Human Geography Label, a label reference layer including cities and communities, countries, administrative units, and at larger scales street names.Human Geography Detail, a detail reference layer including administrative boundaries, roads and highways, and larger bodies of water. This layer is designed to be used with a high degree of transparency so that the detail does not compete with your information. It is set at approximately 50% in this web map, but can be adjusted.Human Geography Base, a simple basemap consisting of land areas in a very light gray only.The vector tile layers in this web map are built using the same data sources used for other Esri Vector Basemaps. For details on data sources contributed by the GIS community, view the map of Community Maps Basemap Contributors. Esri Vector Basemaps are updated monthly.Learn more about this basemap from the cartographic designer in Introducing a Human Geography Basemap.Use this MapThis map is designed to be used as a basemap for overlaying other layers of information or as a stand-alone reference map. You can add layers to this web map and save as your own map. If you like, you can add this web map to a custom basemap gallery for others in your organization to use in creating web maps. If you would like to add this map as a layer in other maps you are creating, you may use the tile layer item referenced in this map.

Reference values from [18] are provided in the top row, mean and standard deviation (St. Dev.) across all Urban Atlas cities are provided on top of the individual city listings.

A feature layer of world cities with labels, for illustrative purposes only for use as a reference layer. This feature layer is pointing to the Political_Map_World_Cities_Features layer provided by Maps.com. The symbology and labels were modified slightly in this version.This layer is used as a reference layer in NOAA NCEI's VIIRS Nighttime Imagery map viewer, displayed in the 3D global view.

World Cities represents a base map layer of the locations of cities for the world. The cities include national capitals, provincial capitals, major population centers, and landmark cities.

Planisphere with the location of the largest cities

The dataset and the validation are fully described in a Nature Scientific Data Descriptor https://www.nature.com/articles/s41597-019-0265-5

If you want to use this dataset in an interactive environment, then use this link https://mybinder.org/v2/gh/GeographerAtLarge/TravelTime/HEAD

The following text is a summary of the information in the above Data Descriptor.

The dataset is a suite of global travel-time accessibility indicators for the year 2015, at approximately one-kilometre spatial resolution for the entire globe. The indicators show an estimated (and validated), land-based travel time to the nearest city and nearest port for a range of city and port sizes.

The datasets are in GeoTIFF format and are suitable for use in Geographic Information Systems and statistical packages for mapping access to cities and ports and for spatial and statistical analysis of the inequalities in access by different segments of the population.

These maps represent a unique global representation of physical access to essential services offered by cities and ports.

The datasets travel_time_to_cities_x.tif (where x has values from 1 to 12) The value of each pixel is the estimated travel time in minutes to the nearest urban area in 2015. There are 12 data layers based on different sets of urban areas, defined by their population in year 2015 (see PDF report).

travel_time_to_ports_x (x ranges from 1 to 5)

The value of each pixel is the estimated travel time to the nearest port in 2015. There are 5 data layers based on different port sizes.

Format Raster Dataset, GeoTIFF, LZW compressed Unit Minutes

Data type Byte (16 bit Unsigned Integer)

No data value 65535

Flags None

Spatial resolution 30 arc seconds

Spatial extent

Upper left -180, 85

Lower left -180, -60 Upper right 180, 85 Lower right 180, -60 Spatial Reference System (SRS) EPSG:4326 - WGS84 - Geographic Coordinate System (lat/long)

Temporal resolution 2015

Temporal extent Updates may follow for future years, but these are dependent on the availability of updated inputs on travel times and city locations and populations.

Methodology Travel time to the nearest city or port was estimated using an accumulated cost function (accCost) in the gdistance R package (van Etten, 2018). This function requires two input datasets: (i) a set of locations to estimate travel time to and (ii) a transition matrix that represents the cost or time to travel across a surface.

The set of locations were based on populated urban areas in the 2016 version of the Joint Research Centre’s Global Human Settlement Layers (GHSL) datasets (Pesaresi and Freire, 2016) that represent low density (LDC) urban clusters and high density (HDC) urban areas (https://ghsl.jrc.ec.europa.eu/datasets.php). These urban areas were represented by points, spaced at 1km distance around the perimeter of each urban area.

Marine ports were extracted from the 26th edition of the World Port Index (NGA, 2017) which contains the location and physical characteristics of approximately 3,700 major ports and terminals. Ports are represented as single points

The transition matrix was based on the friction surface (https://map.ox.ac.uk/research-project/accessibility_to_cities) from the 2015 global accessibility map (Weiss et al, 2018).

Code The R code used to generate the 12 travel time maps is included in the zip file that can be downloaded with these data layers. The processing zones are also available.

Validation The underlying friction surface was validated by comparing travel times between 47,893 pairs of locations against journey times from a Google API. Our estimated journey times were generally shorter than those from the Google API. Across the tiles, the median journey time from our estimates was 88 minutes within an interquartile range of 48 to 143 minutes while the median journey time estimated by the Google API was 106 minutes within an interquartile range of 61 to 167 minutes. Across all tiles, the differences were skewed to the left and our travel time estimates were shorter than those reported by the Google API in 72% of the tiles. The median difference was −13.7 minutes within an interquartile range of −35.5 to 2.0 minutes while the absolute difference was 30 minutes or less for 60% of the tiles and 60 minutes or less for 80% of the tiles. The median percentage difference was −16.9% within an interquartile range of −30.6% to 2.7% while the absolute percentage difference was 20% or less in 43% of the tiles and 40% or less in 80% of the tiles.

This process and results are included in the validation zip file.

Usage Notes The accessibility layers can be visualised and analysed in many Geographic Information Systems or remote sensing software such as QGIS, GRASS, ENVI, ERDAS or ArcMap, and also by statistical and modelling packages such as R or MATLAB. They can also be used in cloud-based tools for geospatial analysis such as Google Earth Engine.

The nine layers represent travel times to human settlements of different population ranges. Two or more layers can be combined into one layer by recording the minimum pixel value across the layers. For example, a map of travel time to the nearest settlement of 5,000 to 50,000 people could be generated by taking the minimum of the three layers that represent the travel time to settlements with populations between 5,000 and 10,000, 10,000 and 20,000 and, 20,000 and 50,000 people.

The accessibility layers also permit user-defined hierarchies that go beyond computing the minimum pixel value across layers. A user-defined complete hierarchy can be generated when the union of all categories adds up to the global population, and the intersection of any two categories is empty. Everything else is up to the user in terms of logical consistency with the problem at hand.

The accessibility layers are relative measures of the ease of access from a given location to the nearest target. While the validation demonstrates that they do correspond to typical journey times, they cannot be taken to represent actual travel times. Errors in the friction surface will be accumulated as part of the accumulative cost function and it is likely that locations that are further away from targets will have greater a divergence from a plausible travel time than those that are closer to the targets. Care should be taken when referring to travel time to the larger cities when the locations of interest are extremely remote, although they will still be plausible representations of relative accessibility. Furthermore, a key assumption of the model is that all journeys will use the fastest mode of transport and take the shortest path.

This web map references the live tiled map service from the OpenStreetMap project. OpenStreetMap (OSM) is an open collaborative project to create a free editable map of the world. Volunteers gather location data using GPS, local knowledge, and other free sources of information such as free satellite imagery, and upload it. The resulting free map can be viewed and downloaded from the OpenStreetMap server: http://www.OpenStreetMap.org. See that website for additional information about OpenStreetMap. It is made available as a basemap for GIS work in Esri products under a Creative Commons Attribution-ShareAlike license.Tip: This service is one of the basemaps used in the ArcGIS.com map viewer and ArcGIS Explorer Online. Simply click one of those links to launch the interactive application of your choice, and then choose Open Street Map from the Basemap control to start using this service. You'll also find this service in the Basemap gallery in ArcGIS Explorer Desktop and ArcGIS Desktop 10.

The Nova Map (World Edition) web map provides a detailed world basemap featuring a dark background with glowing blue symbology and colors that are reminiscent of science-fiction shows, where one is looking at a map of the world on a 'head's up' device or a map that would be projected from a transparent glass wall. The map is designed with a grid pattern across the ocean and stripes or square stippled patterns for land use features visible at larger scales. Additional graphics in the oceans presents a futuristic user interface. The futuristic and less terrestrial feel theme continues with the geometric patterns, starburst city dot symbols, and cool color scheme. The fonts displayed are clean and squarish (san serif) with a futuristic, science-fiction, or high technology appearance.This basemap, included in the ArcGIS Living Atlas of the World, uses the Nova vector tile layer.The vector tile layer in this web map is built using the same data sources used for other Esri Vector Basemaps. For details on data sources contributed by the GIS community, view the map of Community Maps Basemap Contributors. Esri Vector Basemaps are updated monthly.Use this MapThis map is designed to be used as a basemap for overlaying other layers of information or as a stand-alone reference map. You can add layers to this web map and save as your own map. If you like, you can add this web map to a custom basemap gallery for others in your organization to use in creating web maps. If you would like to add this map as a layer in other maps you are creating, you may use the tile layer referenced in this map.

This vector web map features outline maps of the World. The maps can be used for coloring and other fun activities by budding cartographers. These outline maps are great for teaching children about our World. Have them color and label countries, regions and bodies of water. Limited labels appear on the map at large scales. After coloring the city maps, children can do further research to learn more about these places. These maps are also available in a printable PDF format. See this blog with more details on how to work with the vector maps in ArcGIS Pro.For other creatively designed Esri vector basemaps, see the ArcGIS Living Atlas of the World gallery.

This 24-page booklet features outline maps of the World. The maps can be used for coloring and other fun activities by budding cartographers. These outline maps are great for teaching children about our World. Have them color and label countries, regions and bodies of water. After coloring the city maps, children can do further research to learn more about these places. Map pages created from this web map available on ArcGIS.com.See this blog with more details on how to work with the vector maps in ArcGIS Pro.

The Human Geography Dark Map (World Edition) web map provides a detailed world basemap with a dark monochromatic style and content adjusted to support human geography information. Where possible, the map content has been adjusted so that it observes WCAG contrast criteria.This basemap, included in the ArcGIS Living Atlas of the World, uses 3 vector tile layers:Human Geography Dark Label, a label reference layer including cities and communities, countries, administrative units, and at larger scales street names.Human Geography Dark Detail, a detail reference layer including administrative boundaries, roads and highways, and larger bodies of water. This layer is designed to be used with a high degree of transparency so that the detail does not compete with your information. It is set at approximately 50% in this web map, but can be adjusted.Human Geography Dark Base, a simple basemap consisting of land areas in a very dark gray only.The vector tile layers in this web map are built using the same data sources used for other Esri Vector Basemaps. For details on data sources contributed by the GIS community, view the map of Community Maps Basemap Contributors. Esri Vector Basemaps are updated monthly.Learn more about this basemap from the cartographic designer in A Dark Version of the Human Geography Basemap.Use this MapThis map is designed to be used as a basemap for overlaying other layers of information or as a stand-alone reference map. You can add layers to this web map and save as your own map. If you like, you can add this web map to a custom basemap gallery for others in your organization to use in creating web maps. If you would like to add this map as a layer in other maps you are creating, you may use the tile layers referenced in this map.