Using GMT 6.4.0
I have some table data I would like to plot using psxy. The data are already projected in a map projection with units of meters, so I have tried using the -fp[e] flag to indicate this to GMT. But my data wasn’t appearing on my plot. Running in verbose mode, I noticed that GMT uses mapproject quietly in the background to transform my data to lon-lat before plotting. The default behavior of mapproject is to assume that the projected coordinates are relative to the lower-left corner of the map domain. But in my case, the projected coordinates are relative to the projection center. It does not appear to be possible to specify the origin of the coordinate space in the -fp flag. Am I correct?
Of course there are workarounds for this problem – I can preprocess my data with mapproject -C -I and pipe to psxy, or I could add an offset to the x and y coordinates, e.g., using awk. But neither of these options are really satisfying. Is there any other way to specify that my coordinate origin is the projection center?
Thanks, Jed