R/sl.plot.text.R

In helgegoessling/spheRlab: Spherical Geometry, Analysis, and Plotting of Geoscientific Data on Arbitrary Grids

sl.plot.text <-
function (plot.init.res,lon,lat,labels=seq(1,length(lon)),col="black",ignore.visibility=FALSE,adj=NULL,pos=NULL,offset=0.5,vfont=NULL,cex=1,font=NULL,srt=0) {
	
	projection = plot.init.res$projection
	
	if (projection == "platon" || projection == "3D") {
	  
	  if (!is.null(plot.init.res$transform.function)) {
	    lonlat.trans = plot.init.res$transform.function(lon,lat)
	    lon = lonlat.trans$lon
	    lat = lonlat.trans$lat
	  }
		
		npir = 1
		repeat {
			pir = plot.init.res[[npir]]
			if (!is.list(pir)) {return()}
			if (is.null(pir$projection)) {return()}
			sl.plot.text(pir,lon,lat,labels,col,ignore.visibility,adj,pos,offset,vfont,cex,font,srt)
			npir = npir + 1
		}
		
	}
	
	if (!is.null(plot.init.res$transform.function)) {
	  lonlat.trans = plot.init.res$transform.function(lon,lat)
	  lon = lonlat.trans$lon
	  lat = lonlat.trans$lat
	}
	
	vsr.res = sl.vis.shift.rot(plot.init.res,lon,lat)
	x = vsr.res$x
	y = vsr.res$y
	rot.lon = vsr.res$rot.lon
	rot.lat = vsr.res$rot.lat
	visible = vsr.res$visible
	
	if (ignore.visibility) {
		visible[] = TRUE
	}
	
	xshift = plot.init.res$xshift
	yshift = plot.init.res$yshift
	
	if (any(visible)) {
	  x = x[visible]
	  y = y[visible]
	  labels = labels[visible]
	  if (plot.init.res$projection == "mollweide") {
	    xy = sl.proj.mollweide(lon = x, lat = y)
	    x = xy$x
	    y = xy$y
	  }
		text(x=x+xshift,y=y+yshift,labels=labels,adj=adj,pos=pos,offset=offset,vfont=vfont,cex=cex,col=col,font=font,srt=srt)
	}
	
}