# Problem with large buffers in Geopandas

I'm trying to create a large buffer (10,000km) around a specific point using Geopandas but the resulting shape is incorrect. I've changed the `crs` of the geodataframe to a projected crs (`epsg:4087`, World Equidistant Cylindrical, meters) before calculating the buffer, but the result is still not right. Here's the code:

``````import pandas as pd
import numpy as np
import geopandas as gpd

# This is the point around which I want to generate the buffer
pyongyang = pd.DataFrame({'city': ['Pyongyang'], 'country': ["Democratic People's Republic of Korea"], 'longitude': [125.7552106631059], 'latitude': [39.038430756002995]})

# Create the GeoDataFrame
pyongyangGeo = gpd.GeoDataFrame(pyongyang, geometry=gpd.points_from_xy(pyongyang.longitude, pyongyang.latitude), crs='epsg:4326')

# Change the crs to a projected crs in meters
pyongyangGeo = pyongyangGeo.to_crs('epsg:4087')

# Calculate the buffer (10,000km)
maxReach = pyongyangGeo.buffer(10000000)

# Create GeoDataFrame with the buffer
maxReachGeo = gpd.GeoDataFrame(pyongyang, geometry=maxReach, crs='epsg:4087')

# Change the crs for the buffer to WGS 84 to visualize it
maxReachGeo = maxReachGeo.to_crs('epsg:4326')

# Plot it
maxReachGeo.plot()
``````

And this is the resulting plot:

As you can see there, the buffer seems to wrap around itself and remove the overlapping portions.

• You probably want to use a projection centred on or near N. Korea, other wise you run in to edge effects when your buffer falls of the edge of the world Commented Jun 21, 2021 at 16:02

As answered above in the comments, the problem was in the projection. After choosing a projection centered near North Korea (`epsg:5179`) the buffer works well.