shakasom

countries.crs

shakasom

# Look at the geometry column: decimal degrees
countries.geometry[:5]

shakasom

# Project the data into Mercator Projection epsg=3395
countries_projected = countries.to_crs({'init': 'epsg:3395'})

# See the geometry column of the projected countries
countries_projected['geometry'][:5] 

shakasom

fig, ax = plt.subplots(figsize=(14,12))
countries_projected.plot(ax=ax)
# Cities are still in WGS84
cities.plot(ax=ax, color='red'); 

shakasom

# Save projected countries 
shp = '1-introData/countries_epsg3395.shp'
countries_projected.to_file(shp)

shakasom

# Read the data

# 1. Read Countries --> Geopackage Format
file = '1-introData/countries.gpkg'
countries = gpd.read_file(file)

# Cities: Point data
cities = gpd.read_file('1-introData/cities.geojson')

shakasom

# GeoDataFrame
geoDataFrame = countries[['geometry', 'NAME']]
geoDataFrame.head()
type(geoDataFrame)

shakasom

# GeoDataFrame
geoDataFrame = countries[['geometry', 'NAME']]
geoDataFrame.head()
type(geoDataFrame)

shakasom

# Series
series = countries['NAME']
series[:5]
print(type(series))

# GeoSeries
geoSeries = countries['geometry']
geoSeries[:5]
print(type(geoSeries))

shakasom

fig, ax = plt.subplots(figsize=(12,12))
geoDataFrame.centroid.plot(ax=ax)
plt.tight_layout()
plt.show()