map2lonlat(xusr, yusr, tolerance=1e-4)longitude and latitdue.mapproject function.  The inversion is done by using
    optim to find a good match between a point in geographical
    space and plot space, and this match is not always good, being dependent on
    a starting guess.mapPlot for general information on plotting maps,
    including other functions.