gmane.comp.gis.gmt.user

Re: GMT Installation problem?

Paul Duncan — 2008-08-30T08:06:40

Hi Paul, On Sat, 2008-08-30 at 03:50 +0100, Paul Wessel wrote: Yup, that seems to have fixed that problem. Thanks very much for your help! Paul ~~~~~ To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: GMT Installation problem?

Paul Wessel — 2008-08-30T02:50:15

The key file is in the GMT_share archive which is one of the required installs (src + share + coast). -p On Aug 29, 2008, at 8:59 AM, Paul Duncan wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: GMT Installation problem?

Paul Duncan — 2008-08-29T18:59:10

Hi Paul, Thanks for getting back to me. On Fri, 2008-08-29 at 19:41 +0100, Paul Wessel wrote: I did try the solution proposed by Jay Theriot earlier (had to use find, as there's no locate on this box), but the file is not on the machine. The exact sequence of events for this install was 1) Untar GMT4.3.1_src.tar 2) Do the ./configure (with the destination directory set to /opt/GMT4.3.1. 3) make 4) make install as root 5) untar the various coastline files under /opt/GMT4.3.1. Thats it. Am I missing anything? I'm considering that since they're just text I can simply copy the files across from the working installation on my Kubuntu Linux system. Best Regards, Paul ~~~~~ To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: GMT Installation problem?

Paul Wessel — 2008-08-29T18:41:19

Strange, perhaps you are mixing old and new GMT components? PS_font_info.d lives in the share/pslib directory and that is where the programs expect to find it. It used to live in share a long time a go. -p On Aug 29, 2008, at 1:28 AM, Paul Duncan wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: grdimage can not plot all data under a polar coordinate system

Milton P. Plasencia Linares — 2008-08-29T17:30:06

Hi, In the psxy command (last) you must to cancel -K, because you not add more code below. Cheers, Milton ********************************** Milton P. Plasencia Linares Istituto Nazionale di Oceanografia e di Geofisica Sperimentale - OGS Borgo Grotta Gigante 42/C (34010) Sgonico - TRIESTE - ITALIA Tel: +39-040-2140305 Fax: +39-040-2140266 E-mail: mplasencia< at >ogs.trieste.it ASAIN (Antarctic Seismographic Argentinean Italian Network) ********************************* Quoting sqz2 SUBSEA7.COM>: Dear GMT Helper I'm trying to use the following script to plot satellite data collected by a station under a polar coordinate system. But the data whose elevations are larger than about 70 degree can not been displayed by "grdimage". The satellite track can be plot integrally by psxy. I can not solve the problem that the data over 70 degree elevation can not be plotted. Any suggestion and comments will be appreciated. Moreover, the test.dat is attached. http://www.nabble.com/file/p19224129/test.dat test.dat Best regards Qiang #!/bin/csh makecpt -Crainbow -T-1.6/1.6/0.05 -Z >! rain.cpt xyz2grd test.dat -Gtest.grd -R0/360/0/90 -I1/0.5 grdimage test.grd -R0/360/0/90 -JPa4ir -B45g45/10g10 -Crain.cpt -P -K >! test.ps psscale -Crain.cpt -D4.5i/2.0i/3.50i/0.2i -O -K -B0.4/:m: >> test.ps psxy test.dat -O -R -J -Sc.01 -L -Glightgrey -K >> test.ps gs test.ps & rm -f .gmt* *.cpt *.tmp -- View this message in context: http://www.nabble.com/grdimage-can-not-plot-all-data-under-a-polar-coordinate-system-tp19224129p19224129.html Sent from the GMT - User mailing list archive at Nabble.com. To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Email secured by **CeSIT** Check Point gateway ------------------------------------------------- This mail sent through IMP: http://horde.org/imp/ To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

grdimage can not plot all data under a polar coordinate system

sqz2 — 2008-08-29T17:18:58

Dear GMT Helper I'm trying to use the following script to plot satellite data collected by a station under a polar coordinate system. But the data whose elevations are larger than about 70 degree can not been displayed by "grdimage". The satellite track can be plot integrally by psxy. I can not solve the problem that the data over 70 degree elevation can not be plotted. Any suggestion and comments will be appreciated. Moreover, the test.dat is attached. http://www.nabble.com/file/p19224129/test.dat test.dat Best regards Qiang #!/bin/csh makecpt -Crainbow -T-1.6/1.6/0.05 -Z >! rain.cpt xyz2grd test.dat -Gtest.grd -R0/360/0/90 -I1/0.5 grdimage test.grd -R0/360/0/90 -JPa4ir -B45g45/10g10 -Crain.cpt -P -K >! test.ps psscale -Crain.cpt -D4.5i/2.0i/3.50i/0.2i -O -K -B0.4/:m: >> test.ps psxy test.dat -O -R -J -Sc.01 -L -Glightgrey -K >> test.ps gs test.ps & rm -f .gmt* *.cpt *.tmp

grdimage can not plot all data under a polar coordinate system

sqz2 — 2008-08-29T17:16:09

Dear GMT Helper I'm trying to use the following script to plot satellite data collected by a station under a polar coordinate system. But the data whose elevations are larger than about 70 degree can not been displayed by "grdimage". The satellite track can be plot integrally by psxy. I can not solve the problem that the data over 70 degree elevation can not be plotted. Any suggestion and comments will be appreciated. Moreover, the test.dat is attached. Best regards Qiang #!/bin/csh makecpt -Crainbow -T-1.6/1.6/0.05 -Z >! rain.cpt xyz2grd test.dat -Gtest.grd -R0/360/0/90 -I1/0.5 grdimage test.grd -R0/360/0/90 -JPa4ir -B45g45/10g10 -Crain.cpt -P -K >! test.ps psscale -Crain.cpt -D4.5i/2.0i/3.50i/0.2i -O -K -B0.4/:m: >> test.ps psxy test.dat -O -R -J -Sc.01 -L -Glightgrey -K >> test.ps gs test.ps & rm -f .gmt* *.cpt *.tmp http://www.nabble.com/file/p19224083/test.dat test.dat

grdimage can not display a whole image under polar coordinate system

sqz2 — 2008-08-29T17:05:07

Dear GMT Helper I'm trying to use the following script to plot satellite data collected by a station under a polar coordinate system. But the data whose elevations are larger than about 70 degree can not been displayed by "grdimage". The satellite track can be plot integrally by psxy. I can not solve the problem that the data over 70 degree elevation can not be plotted. Any suggestion and comments will be appreciated. Moreover, the test.dat is attached. http://www.nabble.com/file/p19223914/test.dat test.dat Best regards Qiang #!/bin/csh makecpt -Crainbow -T-1.6/1.6/0.05 -Z >! rain.cpt xyz2grd test.dat -Gtest.grd -R0/360/0/90 -I1/0.5 grdimage test.grd -R0/360/0/90 -JPa4ir -B45g45/10g10 -Crain.cpt -P -K >! test.ps psscale -Crain.cpt -D4.5i/2.0i/3.50i/0.2i -O -K -B0.4/:m: >> test.ps psxy test.dat -O -R -J -Sc.01 -L -Glightgrey -K >> test.ps gs test.ps & rm -f .gmt* *.cpt *.tmp

Re: GMT Installation problem?

Jay Theriot — 2008-08-29T15:35:22

I did a "locate PS_font_info.d" and then made a symbolic link to the file in the share directory. "Paul Duncan" wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Average grd file

Eric Fielding — 2008-08-29T15:31:27

Yacob, If you want coarser SRTM data, why not just get the SRTM30 product that already has 30 arcsecond resolution? That was produced as the full average of the 30x30 arcsecond area centered on each grid point. Your method of averaging only a small area will certainly not keep the total mass the same because of the aliasing that Paul mentioned. Regards, ++Eric On Aug 29, 2008, at 3:43 AM, yacob sen wrote: ----- Eric Fielding Eric.J.Fielding< at >jpl.nasa.gov http://science.jpl.nasa.gov/people/Fielding/ MS 300-233 direct phone: +1 (818) 354-9305 Jet Propulsion Lab fax: +1 (818) 354-9476 4800 Oak Grove Drive Pasadena, CA 91109 USA To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

GMT Installation problem?

Paul Duncan — 2008-08-29T11:28:10

Hello, I've recently had a bash at installing GMT4.3.1 on a Sun Blade 1500 running Solaris 10. The installation (including all three coastline fiels) appeared to have gone correctly, until I tried to run a command (pscoast), when I got the following error: GMT Fatal Error: /opt/GMT4.3.1/share/PS_font_info.d: No such file or directory I did a bit of searching around and it seems that others have had the same error message, but I couldn't find a definate solution. Is there a seperate package with fonts in that has to be installed? For info: I couldn't get the install_gmt script to work. I'm on a ship in the north Atlantic at the moment, so I did the individual downloads manually and that eventually worked. Thanks! Paul ~~~~~ To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Average grd file

yacob sen — 2008-08-29T10:43:14

Dear Paul, Thank you for your comments. My files are srtm 3 sec grd tiles (400 tiles in total) with grid registration mode. The idea is to average the nearest points in my case (9 points in each side of a given point with in the square, so the average would be for on single point = SUM(361_points)/361 , this I think leads me to a coarser resolution of 60 seconds). So in my case to work out the averaging of points in one single tiles, I need also data from 8 neighboring tiles to deal with the edges points. I hope this does not introduce aliasing and also this will keep the total mass the same both in the srtm 3 seconds and now the 60 seconds resolution product. Can this be done in GMT ? Kind regards yacob ----- Original Message ---- From: Paul Wessel hawaii.edu> To: GMT-HELP< at >HAWAII.EDU Sent: Friday, August 29, 2008 12:30:51 AM Subject: Re: [GMT-HELP] Average grd file Yacob- What you describe would introduce aliasing. That is never a good idea. If you want a coarse grid of 60 secs to be meaningful you need to average lots of points (e.g., a 60 sec full-width filter). -p On Aug 28, 2008, at 1:14 PM, yacob sen wrote: Dear All, I have a grd file with a 3 second resolution. I would like to make a simple average taking only the nearest few points to change to a coarser.grid to 60 seconds. I am not mean to use grdsample or grdfilter but a simple averaging. Can it be done using GMT. I am here intending to averaging a very big number of grd files. Regards Yacob. ----- Original Message ---- From: Eric Fielding AMPERSAND.JPL.NASA.GOV> To: GMT-HELP< at >HAWAII.EDU Sent: Thursday, August 28, 2008 6:59:27 PM Subject: Re: [GMT-HELP] Wrong Topogrphic Data of Tenerife Hello Michael, The name of the SRTM3 data tile is the location of the southwest corner of the tile, so you used the wrong region specification. It should have been: xyz2grd N28W017.hgt -R-17/-16/28/29 -I3c -N-32768 -ZTLhw -GN28W017.grd xyz2grd N27W017.hgt -R-17/-16/27/28 -I3c -N-32768 -ZTLhw -GN27W017.grd This explains why it was one degree off! ++Eric On Aug 28, 2008, at 9:39 AM, Michael Kremer wrote: ----- Eric Fielding Eric.J.Fielding< at >jpl.nasa.gov http://science.jpl.nasa.gov/people/Fielding/ MS 300-233 direct phone: +1 (818) 354-9305 Jet Propulsion Lab fax: +1 (818) 354-9476 4800 Oak Grove Drive Pasadena, CA 91109 USA 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 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

Re: Average grd file

Paul Wessel — 2008-08-28T23:30:51

Yacob- What you describe would introduce aliasing. That is never a good idea. If you want a coarse grid of 60 secs to be meaningful you need to average lots of points (e.g., a 60 sec full-width filter). -p On Aug 28, 2008, at 1:14 PM, yacob sen wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Average grd file

yacob sen — 2008-08-28T23:14:35

Dear All, I have a grd file with a 3 second resolution. I would like to make a simple average taking only the nearest few points to change to a coarser.grid to 60 seconds. I am not mean to use grdsample or grdfilter but a simple averaging. Can it be done using GMT. I am here intending to averaging a very big number of grd files. Regards Yacob. ----- Original Message ---- From: Eric Fielding AMPERSAND.JPL.NASA.GOV> To: GMT-HELP< at >HAWAII.EDU Sent: Thursday, August 28, 2008 6:59:27 PM Subject: Re: [GMT-HELP] Wrong Topogrphic Data of Tenerife Hello Michael, The name of the SRTM3 data tile is the location of the southwest corner of the tile, so you used the wrong region specification. It should have been: xyz2grd N28W017.hgt -R-17/-16/28/29 -I3c -N-32768 -ZTLhw -GN28W017.grd xyz2grd N27W017.hgt -R-17/-16/27/28 -I3c -N-32768 -ZTLhw -GN27W017.grd This explains why it was one degree off! ++Eric On Aug 28, 2008, at 9:39 AM, Michael Kremer wrote: ----- Eric Fielding Eric.J.Fielding< at >jpl.nasa.gov http://science.jpl.nasa.gov/people/Fielding/ MS 300-233 direct phone: +1 (818) 354-9305 Jet Propulsion Lab fax: +1 (818) 354-9476 4800 Oak Grove Drive Pasadena, CA 91109 USA 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

Re: old version

Paul Wessel — 2008-08-28T07:56:35

Check to see if one of the ftp sites has the old tarballs. Why would you want and older and more buggy version? -p On Aug 27, 2008, at 9:00 PM, Chandan Bora wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

pspolar - take-off angle

Daniela Kuehn — 2008-08-28T07:29:48

Hei, I'm trying to plot P-wave first-motion polarities on the inferior focal half-sphere by using "pspolar". The input to pspolar is meant to be (in columns): 1: station code 2: azimuth 3: take-off angle 4: polarity The azimuth seems to be measured clockwise from North, as expected. I'm not sure if I understand the meaning of the take-off angle correctly. So far I took it to be measured from the vertical upwards (as described in Aki&Richards), zero degrees meaning a ray leaving the source downwards, ninety degrees meaning leaving it horizontally, more than ninety degrees leaving it upwards. Onsets are plotted on the half-sphere by smaller and larger circles. I'd thought that larger points denote stations which directly lie on the inferior half-sphere whereas the smaller points are stations which are projected from the upper to the lower half-sphere. But instead a station related to a take-off angle of e.g. 100 degrees (which in my definition is situated on the upper half-sphere and is therefore projected down) is plotted with a larger circle...so I'm wrong with one of my conclusions...but with which one??? Thank you very much for your help!!!!! Daniela Kuehn To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

old version

Chandan Bora — 2008-08-28T07:00:39

How can I obtain GMT4.0? To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Wrong Topogrphic Data of Tenerife

Eric Fielding — 2008-08-28T17:59:27

Hello Michael, The name of the SRTM3 data tile is the location of the southwest corner of the tile, so you used the wrong region specification. It should have been: xyz2grd N28W017.hgt -R-17/-16/28/29 -I3c -N-32768 -ZTLhw -GN28W017.grd xyz2grd N27W017.hgt -R-17/-16/27/28 -I3c -N-32768 -ZTLhw -GN27W017.grd This explains why it was one degree off! ++Eric On Aug 28, 2008, at 9:39 AM, Michael Kremer wrote: ----- Eric Fielding Eric.J.Fielding< at >jpl.nasa.gov http://science.jpl.nasa.gov/people/Fielding/ MS 300-233 direct phone: +1 (818) 354-9305 Jet Propulsion Lab fax: +1 (818) 354-9476 4800 Oak Grove Drive Pasadena, CA 91109 USA To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Wrong Topogrphic Data of Tenerife

Michael Kremer — 2008-08-28T16:39:34

Re: Rotated: Geotiff to gmt

Eduardo A. Suárez — 2008-08-28T04:17:50

Hi Duncan, as you say, this isn't a gdal list, but from time to time we need it. Here are the commands: gdalwarp GMTtest.tiff GMTtest_northup.tiff gdalwarp -t_srs "+proj=latlong" GMTtest_northup.tiff GMTtest_latlon.tiff gdal_translate -b 1 -of GMT GMTtest_latlon.tiff GMTtest_latlon_red.grd gdal_translate -b 2 -of GMT GMTtest_latlon.tiff GMTtest_latlon_green.grd gdal_translate -b 3 -of GMT GMTtest_latlon.tiff GMTtest_latlon_blue.grd psbasemap -R136.3105940/136.3106153/-34.5029822/-34.5029652 -JM15 -Ba0.0001 -K -V > GMTtest.ps grdimage GMTtest_latlon_red.grd GMTtest_latlon_green.grd GMTtest_latlon_blue.grd -R -J -O -V >> GMTtest.ps The coordinates for -R were obtained with gdalinfo. Please download and use the binarias from http://fwtools.maptools.org (FWTools 2.0.6). Best regards, Edu.-

Re: Rotated: Geotiff to gmt

Paul Wessel — 2008-08-28T01:46:55

Hi Duncan- This is not GMT's strong suite. However, if you can reproject so that gdal can output a GMT grid then you need to get three grid (r, g, b) and you can then plot them with grdimage which can handle r, g, b grids (three at once). -p On Aug 27, 2008, at 3:39 PM, Duncan Mercer wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu