gmane.comp.gis.gmt.user

Mark Collier | 24 Sep 03:29

problem drawing line with psxy

Hi,

I am using an SGI Altix 4700 system running SLES10.

We have just installed GMT 4.3.1, however, the problem existed with an earlier version 4.2.0. The problem is that a strange extra leg is drawn rather than just the line between New York and London. I guess it is something to do with crossing the main meridian. I have checked the coordinates with Google Earth and tried alternative ways of specifying the lat/lon (signed coordinates), as well as other projections.

pscoast -P -Y10c -K -R-180/180/-60/90 -G205 -Dc -W0.25p -C -B30/30 -Jcyl_stere/140/0/0.015i > globe.ps

psxy test.in -A -O -M -R-180/180/-60/90 -Jcyl_stere/140/0/0.015i -W2p,red -Gred >> globe.ps

The input file test.in has the contents:

73.9666666667w,40.7833333333n,0.2,New York

0.0833333333333w,51.5333333333n,0.2,London

Either a bug or I’m doing something silly. Regards,

Mark.

To unsubscribe, send the message "signoff gmt-help" to listserv <at> hawaii.edu

headers

Paul Wessel | 24 Sep 10:42

Re: problem drawing line with psxy

Hi Mark-

Mostly silly: You are specifying -Gred to fill a polygon red. I am not sure what that means with 2 points but removing the -Goption seems to fix the problem.

-p

On Sep 23, 2008, at 3:33 PM, Mark Collier wrote:

Hi,

I am using an SGI Altix 4700 system running SLES10.

We have just installed GMT 4.3.1, however, the problem existed with an earlier version 4.2.0. The problem is that a strange extra leg is drawn rather than just the line between New York and London. I guess it is something to do with crossing the main meridian. I have checked the coordinates with Google Earth and tried alternative ways of specifying the lat/lon (signed coordinates), as well as other projections.

pscoast -P -Y10c -K -R-180/180/-60/90 -G205 -Dc -W0.25p -C -B30/30 -Jcyl_stere/140/0/0.015i > globe.ps

psxy test.in -A -O -M -R-180/180/-60/90 -Jcyl_stere/140/0/0.015i -W2p,red -Gred >> globe.ps

The input file test.in has the contents:

73.9666666667w,40.7833333333n,0.2,New York

0.0833333333333w,51.5333333333n,0.2,London

Either a bug or I’m doing something silly. Regards,

Mark.

To unsubscribe, send the message "signoff gmt-help" to listserv <at> hawaii.edu

To unsubscribe, send the message "signoff gmt-help" to listserv <at> hawaii.edu

headers

Patrick Sibanda | 24 Sep 11:05

Re: Labeling the latitude and longitude

Dear Paul,

I asked this question a few days ago but I have not found any help yet. 
Is there a way to doing this?

I would like to label my axes as latitude and longitude but this is not 
working. I am trying to use the following piece of code  and I expect is 
to label my axes but this does seem to work, am I doing something wrong?

Here is the piece of code

psbasemap -R-180/180/-90/90 -JX25.8d/15d 
-Ba30:"Latitude":/30:"Longitude":WESN -K -X1.5 -Y1.5 > test.eps
pscoast -J -R -Dc -N1 -G230-0.1-0.45 -A1000 -S0/0/100 -W0.25p -K -O >> 
test.eps

Please help  me!

Kind regards
Patrick

-- 
========================================
Patrick Sibanda
Dept of Physics and Electronics
Rhodes University

Hermanus Magnetic Observatory
P.O. Box 32
Hermanus 7200
RSA
Email : psibanda <at> hmo.ac.za

Phone : +27 (0)28 312 1196  Web: http://www.star.ac.za
Fax   : +27 (0)28 312 2039     : http://www.hmo.ac.za
--------------------------------------------------------------------------
"Arguments are to be avoided, they are always vulgar and often convincing"

To unsubscribe, send the message "signoff gmt-help" to listserv <at> hawaii.edu

headers

José Antonio Álvarez Gómez | 24 Sep 14:26

Re: Labeling the latitude and longitude

Hi Patrick

I have been trying your code and find that when you don't use the "d" flag (meaning the data is in geographical
coordinates) in the JX the code works properly, although you can't draw the coastline. I tried too making a
Mercator map (JM) and a equidistant cylindrical (JQ), but the axes titles don't appear neither.

It seems that the labels of the axes don't work for geographical projections. I am not a developer of the
program, just a user, and I'm not sure if this is a bug or if this is the way it works.

Maybe you can duplicate the psbasemap in order to plot the labels first and then setting the geographical
"d" flag:

psbasemap -R-180/180/-90/90 -JX25.8/15 -Ba30:"Latitude":/30:"Longitude":WESN -K -X1.5 -Y1.5 > test.eps
psbasemap -R -JX25.8d/15d -Bwesn -K -O >> test.eps
pscoast -J -R -Dc -N1 -G230-0.1-0.45 -A1000 -S0/0/100 -W0.25p -K -O >> test.eps

Regards

Jose

-----Mensaje original-----
De: World-Wide GMT Usage and Help Mailing List [mailto:GMT-HELP <at> HAWAII.EDU] En nombre de Patrick Sibanda
Enviado el: miércoles, 24 de septiembre de 2008 11:06
Para: GMT-HELP <at> HAWAII.EDU
Asunto: Re: [GMT-HELP] Labeling the latitude and longitude

Dear Paul,

I asked this question a few days ago but I have not found any help yet.
Is there a way to doing this?

I would like to label my axes as latitude and longitude but this is not
working. I am trying to use the following piece of code  and I expect is
to label my axes but this does seem to work, am I doing something wrong?

Here is the piece of code

psbasemap -R-180/180/-90/90 -JX25.8d/15d
-Ba30:"Latitude":/30:"Longitude":WESN -K -X1.5 -Y1.5 > test.eps
pscoast -J -R -Dc -N1 -G230-0.1-0.45 -A1000 -S0/0/100 -W0.25p -K -O >>
test.eps

Please help  me!

Kind regards
Patrick

--
========================================
Patrick Sibanda
Dept of Physics and Electronics
Rhodes University

Hermanus Magnetic Observatory
P.O. Box 32
Hermanus 7200
RSA
Email : psibanda <at> hmo.ac.za

Phone : +27 (0)28 312 1196  Web: http://www.star.ac.za
Fax   : +27 (0)28 312 2039     : http://www.hmo.ac.za
--------------------------------------------------------------------------
"Arguments are to be avoided, they are always vulgar and often convincing"

To unsubscribe, send the message "signoff gmt-help" to listserv <at> hawaii.edu

To unsubscribe, send the message "signoff gmt-help" to listserv <at> hawaii.edu