Can't plot data points out of the [-180, 180, -90, 90] range on geographic maps for external wrappers

[seisman]

See https://forum.generic-mapping-tools.org/t/plotting-points-out-of-90-degree-range/5479 for the initial issue report.

To reproduce the issue:

import pygmt
fig = pygmt.Figure()
fig.basemap(region="d", projection="Q10i", frame="afg")
fig.coast(land="#666666", water="skyblue")
fig.plot(x=60, y=110, style="c0.3c", fill="red", pen="black", no_clip=True)
fig.show()

The error plot [WARNING]: Latitude (110) at line # 1 exceeds -|+ 90! - set to NaN comes from the gmtlib_process_binary_input function.

gmt/src/gmt_io.c

Line 4839 in 65e0917

GMT_Report (GMT->parent, GMT_MSG_WARNING, "Latitude (%g) at line # %" PRIu64 " exceeds -|+ 90! - set to NaN\n", GMT->current.io.curr_rec[col_no], GMT->current.io.rec_no);

The equivalent CLI version works as expected. The reason is that, for GMT CLI, gmtio_ascii_input is called to read the data points, but for external wrappers, data are passed in binary buffers, so the gmtio_bin_input is called instead.

I guess the same issue exists for GMT.jl, so it should be declared as a GMT bug.

[joa-quim]

Yep. This works

C:\v>echo 60 110 | gmt plot -Rd -JQ10i -Baf -BWSen -N -Gred -Sc7p -png lixo

but this, not

Julia> plot([60 110], region="d", projection="Q10i", marker="c0.3c", fill="red", no_clip=true, Vd=1)
        psxy  -Rd -JQ10i -Baf -BWSen -N -Gred -Sc7p -P -K > c:\TMP/GMTjl_j.ps
psxy [WARNING]: Latitude (110) at line # 1 exceeds -|+ 90! - set to NaN

[joa-quim]

This is rather binary vs ascii issue than an externals.

echo 60 110 | gmtconvert -bo | gmt plot -bi -Rd -JQ10i -Baf -BWSen -N -Gred -Sc7p -png lixo
plot [WARNING]: Latitude (110) at line # 1 exceeds -|+ 90! - set to NaN

[joa-quim]

The problem is that when gmtio_bin_input() gets called from GMT_Get_Record() in psxy, it knows nothing on who called it nor -N so the lat limits checking pops up. The ascii branch does not have those checks and holds on till psxy checks if the point is to be skipped or not in function of -N.

Not easy to fix this.

[joa-quim]

And is restricted to geogs. Cartesians work

gmt psxy lixo.bin -bi -R-180/180/-90/90 -JX10 -Baf -BWSen -N -Gred -Sc7p -png lixo

[joa-quim]

Fixed by #8663