File:Migrant caravan.pdf

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