Round maps with R - r

Question

I need to draw a round map in an R object in PowerBI.
Now I use two objects, a transparent circle image and a R visual object (map).
This is the code:
dataset <- data.frame(03. Ciudad origen, lon, lat, 02. Ciudad destino, lon.1, lat.1, zoomMapa)
dataset <- unique(dataset)
par(fg="#f0ecf8")
colnames(dataset) <- c("ciudad_base","lon_base","lat_base","apto_dest","lon_dest","lat_dest","zoomMapa")
library(stringr)
library(maps)
library(geosphere)
library(magrittr)
offset <- unique(dataset$zoomMapa)
plot_my_connection=function( dep_lon, dep_lat, arr_lon, arr_lat, ...){
inter <- gcIntermediate(c(dep_lon, dep_lat), c(arr_lon, arr_lat), n=50, addStartEnd=TRUE, > breakAtDateLine=F)
inter=data.frame(inter)
diff_of_lon=abs(dep_lon) + abs(arr_lon)
if(diff_of_lon > 180){
lines(subset(inter, lon>=0), ...)
lines(subset(inter, lon<0), ...)
}else{
lines(inter, ...)
}
}
dif_lon=function(dep_lon, arr_lon){
ret = 0
if (dep_lon< 0 & arr_lon<0 ){
ret = abs(min(dep_lon, arr_lon)-max(dep_lon, arr_lon))
}
if (dep_lon> 0 & arr_lon>0 ){
ret = abs(max(dep_lon, arr_lon)-min(dep_lon, arr_lon))
}
if (dep_lon> 0 & arr_lon<0 ){
ret = dep_lon - arr_lon
}
if (dep_lon<0 & arr_lon>0 ){
ret = - dep_lon + arr_lon
}
return (ret)
}
dif_lat=function(dep_lat, arr_lat){
ret = 0
if (dep_lat< 0 & arr_lat<0 ){
ret = abs(min(dep_lat, arr_lat)-max(dep_lat, arr_lat))
}
if (dep_lat> 0 & arr_lat>0 ){
ret = abs(max(dep_lat, arr_lat)-min(dep_lat, arr_lat))
}
if (dep_lat> 0 & arr_lat<0 ){
ret = dep_lat - arr_lat
}
if (dep_lat<0 & arr_lat>0 ){
ret = - dep_lat + arr_lat
}
return (ret)
}
lonmin= min(min(dataset$lon_base), min(dataset$lon_dest)) - dataset$zoomMapa
lonmax = max(max(dataset$lon_base), max(dataset$lon_dest)) + dataset$zoomMapa
latmin= min(min(dataset$lat_base), min(dataset$lat_dest)) - dataset$zoomMapa
latmax=max(max(dataset$lat_base), max(dataset$lat_dest)) + dataset$zoomMapa
diflon = dif_lon(lonmin,lonmax)
diflat = dif_lat(latmin, latmax)
off_lon = 0
off_lat = 0
if (diflon > diflat)
off_lat = (diflon-diflat)/2
if (diflon < diflat)
off_lon = (diflat-diflon)/2
lonmin= lonmin - off_lon
lonmax = lonmax + off_lon
latmin= latmin - off_lat
latmax = latmax + off_lat
par(mar=c(0,0,0,0))
map('world', boundary = TRUE, regions = "", interior = TRUE, bg="white", col="#f0ecf8", fill=TRUE, mar = c(0.1, 0.1, 0, 0.1), xlim=c(lonmin,lonmax), ylim=c(latmin,latmax) )
points(x=as.double(dataset$lon_base1), y=as.double(dataset$lat_base1), col="#1a2732", cex=1, pch=1)
text(dataset$ciudad_base1, x=as.double(dataset$lon_base1), y=as.double(dataset$lat_base1), col="#1a2732", cex=1, pos=1)
for(i in 1:nrow(dataset)){
plot_my_connection(as.double(dataset$lon_base[i]), as.double(dataset$lat_base[i]), as.double(dataset$lon_dest[i]),as.double( dataset$lat_dest[i]), col="#1a2732", lwd=1)
points(x=as.double(dataset$lon_dest[i]), y=as.double(dataset$lat_dest[i]), col="#1a2732", cex=1, pch=1)
text(dataset$apto_dest[i], x=as.double(dataset$lon_dest[i]), y=as.double(dataset$lat_dest[i]), col="#1a2732", cex=1, pos=1)
}
The main problem is that for some routes I don't get the map filled, only the rectangle part.
Any tip there?
Thank you very much.