Psxy error when specifying "landscape" orientation with gmtset

[Malley1]

Description of the problem

Hi,

I am running gmt6.0.0. I found an error in psxy while drawing lines or polygons - when PS_PAGE_ORIENTATION is set to landscape, strange vertical lines appear on the plot. See attached examples which only differ by gmt gmtset PS_PAGE_ORIENTATION = portrait vs. PS_PAGE_ORIENTATION = landscape. I am aware that the "convert" line also contains a "-rotate 90" flag which is why the images aren't actually appearing in portrait or landscape respectively but that's not the issue.

#!/bin/bash
source /home/co323/.bashrc

gmt gmtset PS_PAGE_ORIENTATION = landscape

testno="35"
outputstring="output_8"

psfile="./timeslice_plots/all_africa_drainage_blt${testno}_${outputstring}_landscape.ps"

drainagedir="/space/co323/landscape_evolution/badlands/inversions/blt${testno}/${outputstring}/data/raw_data"
basins=`ls $drainagedir`

rgn="-R-24/64/-39/44"
proj="-JM12c"

coast1="/space/co323/afr_arab_mad_riv_inv/plotting/coast_full_afr_mad.dat"
d=">"

gmt psbasemap $rgn $proj -Ba0 -K > $psfile

datadir="riv_data_blt${testno}_${outputstring}"
riverlist=`ls $datadir`

coast1="/space/co323/afr_arab_mad_riv_inv/plotting/coast_full_afr_mad.dat"
d=">"
awk -v d=$d '{if (NF==2) print ($1*1000)-5e6, ($2*1000)-5e6 ; else print d}' $coast1 | invproj +proj=aea +lon_0=25e +lat_0=0 +lat_1=23s +lat_2=20n +datum=WGS84 +units=m +no_defs +to +proj=latlon -f '%.8f' | awk -v d=$d '{if ($1 != "25.00000000" && $2 != "-0.00000000>") print $1, $2 ; else print d}' > coastfile_ll.temp

gmt psclip coastfile_ll.temp $rgn $proj -O -K >> $psfile

echo "plotting all_ll_blt${testno}_${outputstring}.txt..."
gmt psxy all_ll_blt${testno}_${outputstring}.txt $rgn $proj -W0.5p,blue -O -K >> $psfile

gmt psclip -C $rgn $proj -O -K >> $psfile

gmt pscoast -Dh -W1p,black -A6000/0/0 -O -K $rgn $proj >> $psfile 


ps2jpg

eog $jpgfile &

ps2jpg is:

ps2jpg () {
	if [ $# -eq 1 ] ; then
		psfile="$1"
	fi
	jpgfile=$(echo $psfile | sed s/"\.ps"/".jpg"/)

	convert -density 300 -quality 100 -rotate 90 -trim -matte -bordercolor white -border 5x5 $psfile $jpgfile
}

The following plot differs only in that the gmtset line now reads:

gmt gmtset PS_PAGE_ORIENTATION = portrait

(No error message)

System information

• Operating system: ubuntu 16.04 LTS
• Version of GMT: 6.0.0

[welcome]

👋 Thanks for opening your first issue here! Please make sure you filled out the template with as much detail as possible. We appreciate that you took the time to contribute!

Please make sure you read our Contributing Guide and abide by our Code of Conduct.

[joa-quim]

@Malley1 thank you for the bug report. But we need something that we can reproduce. Can you make as a minimum working example?

And BTW, no need for a ps2jpg we prefer our psconvert program.

[PaulWessel]

This is clearly a problem with your PS conversion. The PS_ORIENTATION, whatever its setting, only adds a translate/rotate component - it cannot possibly effect how lines are drawn. Yet, if you can make a small subset of your data and a script available that fails for you then we can at least test it.

[Malley1]

Paul and Joaquim thanks for the input, just checked the .ps file and it is indeed fine so the problem is with ImageMagick "convert" - it's just strange that the bug doesn't appear in portrait. I have included the entire gmt script and the "convert" command, and two example images - aren't they visible? The data being plotted is just longitudes and latitudes in an ascii file which I've attached, although it is large.
all_ll_blt35_output_8.txt

[PaulWessel]

Thanks, but since there is not anything we can do about convert etc we will close this issue. Recommend using psconvert instead of converters not designed for high-quality graphics.