data(world)
plot(world)
# use different projections
azimuth_equaldist <- "+proj=aeqd +lat_0=90 +lon_0=45"
mollweide_atlantic <- "+proj=moll +lon_0=11.5"
mollweide_pacific <- "+proj=moll +lon_0=151"
plot(sf::st_transform(world, crs = azimuth_equaldist))
Run the code above in your browser using DataLab