Feel free to improve the map!

(File karavane.r)

#R CMD BATCH karavane.r;pdfcrop migrant_caravan.pdf migrant_caravan.pdf library(rnaturalearth) library(plotrix) liste<-read.table("karavane.tsv",header=TRUE,sep="\t",dec=",") liste[,1]<-as.Date(liste[,1]) zeck<-.75#Schriftgröße vek<-c(1,10,17,18,21,22)#Genannte Städte band<-.95*max(max(liste$Longitude)-min(liste$Longitude),max(liste$Latitude)-min(liste$Latitude))/2#edge length map mittex<-(max(liste$Longitude)+min(liste$Longitude))/2 mittey<-(max(liste$Latitude)+min(liste$Latitude))/2 map<-ne_countries(continent = "north america",scale=10)#Precise map for publication #map<-ne_countries(continent = "north america",scale=110)#Small file for testing farbe<-function(farbe.n){ farbe.f<-"red" if(farbe.n%%2==0){farbe.f<-"green"} farbe.f }#Allow for differen colors per stage #Plotfunktion: bild<-function(){if(require(sp)){plot(map,ylim=c(mittey-band,mittey+band),xlim=c(mittex-band,mittex+band))}}#produces the map marke<-function(tag){ text(liste$Longitude[tag],liste$Latitude[tag],labels=liste$Plak[tag],cex=1,pos=3)#Stadt text(liste$Longitude[tag],liste$Latitude[tag]-.6,labels=liste$KD[tag],cex=zeck,pos=3)#Datum text(liste$Longitude[tag],liste$Latitude[tag]-1.2,labels=liste$People[tag],cex=zeck,pos=3)#Stärke }  dicke<-.00003#stream thickness  punkt<-function(punkt.n,punkt.r){ punkt.p<-c(liste$Longitude[punkt.n],liste$Latitude[punkt.n]) punkt.weg<-c(cos(punkt.r),sin(punkt.r))*dicke*liste$People[punkt.n] punkt.p+punkt.weg }  #flow strich<-function(strich.n){ grund<-c(liste$Longitude[strich.n]-liste$Longitude[strich.n-1],liste$Latitude[strich.n]-liste$Latitude[strich.n-1]) gradient<-atan2(grund[2],grund[1]) länge<-norm(grund,type="2") tüte<-asin((liste$People[strich.n]-liste$People[strich.n-1])*dicke/länge) eck<-matrix(c(1:8),4) eck[1,]<-punkt(strich.n,gradient+tüte+pi/2) eck[2,]<-punkt(strich.n,gradient-tüte-pi/2) eck[3,]<-punkt(strich.n-1,gradient-tüte-pi/2) eck[4,]<-punkt(strich.n-1,gradient+tüte+pi/2) polygon(eck[,1],eck[,2],col=farbe(strich.n),border=farbe(strich.n)) draw.circle(liste$Longitude[strich.n],liste$Latitude[strich.n],liste$People[strich.n]*dicke,col=farbe(strich.n),border=farbe(strich.n)) draw.circle(liste$Longitude[strich.n-1],liste$Latitude[strich.n-1],liste$People[strich.n-1]*dicke,col=farbe(strich.n),border=farbe(strich.n)) }  #lines() #strich<-function(strich.n){lines(c(liste$Longitude[strich.n-1],liste$Longitude[strich.n]),c(liste$Latitude[strich.n-1],liste$Latitude[strich.n]),lwd=liste$People[strich.n-1]/1000,col=farbe(strich.n))}  #makes the path ganz<-function(){ bild()#Plot erstellen lapply(c(2:dim(liste)[1]),strich)#insert path lapply(vek,marke)#Marken einfügen } #ganz() pdf("migrant_caravan.pdf")#Outputdatei öffnen ganz()#Jetzt wird gedruckt dev.off()#Outputdatei schließen 

(File karavane.tsv, as derived from w:en:Central American migrant caravans#Table)

Date	KD	Place	Plak	Latitude	Longitude	People 2018-10-12	12 Oct.	San Pedro Sula	Sula	15,052018	-88,033333	160 2018-10-13	13 Oct.	Santa Rosa de Copan	Copan	14,766660	-88,783333	1000 2018-10-14	14 Oct.	Ocotepeque	Ocotepeque	14,433333	-89,183333	1700 2018-10-15	15 Oct.	Aguas Calientes (Honduras)	Calientes	14,531281	-89,283139	1600 2018-10-16	16 Oct.	Chiquimula	Chiquimula	14,783333	-89,533333	2000 2018-10-18	18 Oct.	Ciudad de Guatemala	Guatemala	14,613333	-90,535278	3000 2018-10-19	19 Oct.	Ciudad Tecun Uman	Tecun Uman	14,666667	-92,150000	3000 2018-10-21	21 Oct.	Tapachula	Tapachula	14,900000	-92,266667	7000 2018-10-24	24 Oct.	Huixtla	Huixtla	15,133333	-92,466667	5800 2018-10-25	25 Oct.	Mapastepec	Mapastepec	15,427778	-92,898333	4500 2018-10-26	26 Oct.	Tonala	Tonala	16,100000	-93,750000	4500 2018-10-27	27 Oct.	Arriaga	Arriaga	16,235000	-93,895556	3500 2018-10-28	28 Oct.	Tapanatepec	Tapanatepec	16,366667	-94,200000	4000 2018-10-30	30 Oct.	Santiago Niltepec	Niltepec	16,566667	-94,616667	4000 2018-10-31	31 Oct.	Juchitán de Zaragoza	Juchitán	16,433333	-95,016667	5000 2018-11-02	2 Nov.	Matías Romero	Matías Romero	16,866667	-95,033333	4000 2018-11-03	3 Nov.	Isla, Veracruz	Isla	18,033333	-95,533333	4000 2018-11-05	5 Nov.	Ciudad de México	México	19,433333	-99,133333	6000 2018-11-10	10 Nov.	Santiago de Querétaro	Querétaro	20,587500	-100,392778	6531 2018-11-11	11 Nov.	Irapuato	Irapuato	20,666667	-101,350000	1 2018-11-13	13 Nov.	Guadalajara	Guadalajara	20,676667	-103,347500	2000 2018-11-15	15 Nov.	Tijuana	Tijuana	32,525000	-117,033333	1500