http://blog.gmane.org/gmane.comp.gis.gmt.user hourly 1 1901-01-01T00:00+00:00 http://gmane.org/img/gmane-25t.png http://gmane.org http://comments.gmane.org/gmane.comp.gis.gmt.user/9891 Hi, I guess this is a bug: Changing the color of the frame, for example with gmtset BASEMAP_FRAME_RGB 0/100/0 has for effect to change also three other variables as following GRID_PEN_PRIMARY = 0.25p,0/100/0 GRID_PEN_SECONDARY = 0.5p,0/100/0 TICK_PEN = 0.5p,0/100/0 However, if after drawing my basemap I reset the frame color to black (gmtset BASEMAP_FRAME_RGB 0/0/0), the other three are not reset and remain as above. Is this wanted? Best regards --- Bernard Giroux, ing., Ph.D. Chercheur http://geo.polymtl.ca/~giroux To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Bernard Giroux 2008-09-07T03:08:56 http://comments.gmane.org/gmane.comp.gis.gmt.user/9888 John Robbins 2008-09-05T15:10:14 http://comments.gmane.org/gmane.comp.gis.gmt.user/9887 Hi all- In our quest for world domination I am pleased to announce that we now have ftp mirror sites on all continents (well, Antarctica's gotta wait a while longer...). Earlier this year we added a new site in Sydney, Australia operated by the School of Geosciences at Sydney University, and now we are pleased to announce a new mirror in Africa hosted by TENET, the Tertiary Education & Research Networks of South Africa. The GMT web pages and install_gmt script have been updated to reflect the new choices. Sincerely, The GMT Team To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Paul Wessel 2008-09-05T04:17:56 http://comments.gmane.org/gmane.comp.gis.gmt.user/9884 I am using the most recent release of gmt available via fink, I believe this is release 4.2.0-2 on an Intel MacBook Pro. I have downloaded grav.16.1 1 minute gravity grid from Sandwell's web site. I want to extract a subset of the data from this file using img2mercgrd. For example, I want to extract a square of data that extends from 0 to 10 East, and -10 to 0 North. If I use the command /sw/bin/img2mercgrd grav.img.16.1 -R0/10/-10/0 -D -Gtoss.grd -T1 -S0.1 -V -m1.0 I get the following img2mercgrd expects grav.img.16.1 to be 21600 by 17280 pixels spanning 0/360.0/-80.738009/80.738009 img2mercgrd: To fit [averaged] input, your -R0/10/-10/0 is adjusted to -R0/10/-10.0152710647/0 img2mercgrd: The output will be 600 by 604 pixels. Created 600 by 604 Mercatorized .grd file. Min, Max values are -69 227.2 This appears to be exactly what I am trying to extract- so I feel good. If I then use grdinfo- rowley% /sw/bin/grdinfo toss.grd toss.grd: Title: Data from Altimetry toss.grd: Command: /sw/bin/img2mercgrd -V -R0/10/-10/0 grav.img.16.1 - D -Gtoss.grd -T1 -S0.1 -m1.0 toss.grd: Remark: Spherical Mercator Projected with -Jm1 - R0/10/-10.0152710647/0 toss.grd: Pixel node registration used toss.grd: grdfile format: nf (# 18) toss.grd: x_min: 0 x_max: 10 x_inc: 0.0166667 name: Spherical Mercator projected Longitude, -Jm1, length from 0 nx: 600 toss.grd: y_min: 0 y_max: 10.0667 y_inc: 0.0166667 name: Spherical Mercator projected Latitude, -Jm1, length from -10.0152710647 ny: 604 toss.grd: z_min: -69 z_max: 227.2 name: meters, mGal, Eotvos, or micro- radians, depending on img file and -S. toss.grd: scale_factor: 1 add_offset: 0 Again, this looks fine, but if I now export this grid as an xyz file using grd2xyz with the command- rowley% /sw/bin/grd2xyz toss.grd -V > toss.xyz grd2xyz: Working on file toss.grd grd2xyz: 362400 values extracted But - if I more toss.xyz I get the following at the top of the file 0.00833333 10.0583 -4 0.025 10.0583 -4 0.0416667 10.0583 -3.8 0.0583333 10.0583 -3.6 0.075 10.0583 -3.4 0.0916667 10.0583 -3.4 0.108333 10.0583 -3.4 0.125 10.0583 -3.4 0.141667 10.0583 -3.8 0.158333 10.0583 -4 0.175 10.0583 -4.2 0.191667 10.0583 -4.2 0.208333 10.0583 -4 0.225 10.0583 -3.6 Thus longitude is fine, but latitude is completely off. The northernmost latitude should be 0, with the southernmost latitude of -10.0152710647 The end of the file toss.xyz is 9.8250.008333335.4 9.841670.008333335.4 9.858330.008333335.2 9.8750.008333335 9.891670.008333334.8 9.908330.008333335 9.9250.008333335 9.941670.008333335.2 9.958330.008333335.4 9.9750.008333335.8 9.991670.008333336 Again longitude is fine but latitude is not where I want my extracted data set to be from. What am I doing wrong. I have read and reread the img2mercgrd and can't find what I am missing. Thanks, David David Rowley Professor Department of the Geophysical Sciences The University of Chicago 5734 S. Ellis Avenue Chicago, IL 60637 phone: 1-773-702-8146 fax:1-773-702-9505 David Rowley Professor Department of the Geophysical Sciences The University of Chicago 5734 S. Ellis Avenue Chicago, IL 60637 phone: 1-773-702-8146 fax:1-773-702-9505 To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu David Rowley 2008-09-04T20:56:25 http://comments.gmane.org/gmane.comp.gis.gmt.user/9881 I will be out of the office starting 09/03/2008 and will not return until 10/06/2008. I am presently at sea on the Coast Guard Cutter Healy. I will be away from the office from September 3 thru October 6,, 2008. During this time my email address will be: jonathan.childs< at >healy.polarscience.net My USGS email will not be forwarded to the boat. While on board the Healy, incoming email messages are restricted to 100 kb, including attachments. Please limit email sent to the Healy to the minimum necessary. Thank you. To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Jonathan R Childs 2008-09-03T17:00:42 http://comments.gmane.org/gmane.comp.gis.gmt.user/9868 Hi, when I run pscoast -R-82/-57/-53/-20r -JA-67.5/-38.5/14 -Df -G255/242/223 -S208/229/242 -W1/0/84/223 -N1/1/192/80/80 -V -O I get this *** glibc detected *** ./pscoast: munmap_chunk(): invalid pointer: 0x081b44d0 *** ======= Backtrace: ========= /lib/i686/cmov/libc.so.6[0xb7da94f4] ./pscoast[0x8075beb] ./pscoast[0x8075d7b] ./pscoast[0x808f39c] ./pscoast[0x804baa2] /lib/i686/cmov/libc.so.6(__libc_start_main+0xe5)[0xb7d51455] ./pscoast[0x8049f21] ======= Memory map: ======== 08048000-08104000 r-xp 00000000 08:01 1367430 /root/gmt/GMT4.3.1/src/pscoast 08104000-08127000 rw-p 000bc000 08:01 1367430 /root/gmt/GMT4.3.1/src/pscoast 08127000-08258000 rw-p 08127000 00:00 0 [heap] b7c86000-b7cd3000 rw-p b7d13000 00:00 0 b7cd3000-b7cd5000 rw-p b7cd3000 00:00 0 b7cda000-b7ce4000 r-xp 00000000 08:01 205311 /lib/i686/cmov/libnss_files-2.7.so b7ce4000-b7ce6000 rw-p 00009000 08:01 205311 /lib/i686/cmov/libnss_files-2.7.so b7ce6000-b7cef000 r-xp 00000000 08:01 205313 /lib/i686/cmov/libnss_nis-2.7.so b7cef000-b7cf1000 rw-p 00008000 08:01 205313 /lib/i686/cmov/libnss_nis-2.7.so b7cf1000-b7d06000 r-xp 00000000 08:01 205306 /lib/i686/cmov/libnsl-2.7.so b7d06000-b7d08000 rw-p 00014000 08:01 205306 /lib/i686/cmov/libnsl-2.7.so b7d08000-b7d0a000 rw-p b7d08000 00:00 0 b7d0a000-b7d11000 r-xp 00000000 08:01 205307 /lib/i686/cmov/libnss_compat-2.7.so b7d11000-b7d13000 rw-p 00006000 08:01 205307 /lib/i686/cmov/libnss_compat-2.7.so b7d3a000-b7d3b000 rw-p b7d3a000 00:00 0 b7d3b000-b7e90000 r-xp 00000000 08:01 205298 /lib/i686/cmov/libc-2.7.so b7e90000-b7e91000 r--p 00155000 08:01 205298 /lib/i686/cmov/libc-2.7.so b7e91000-b7e93000 rw-p 00156000 08:01 205298 /lib/i686/cmov/libc-2.7.so b7e93000-b7e97000 rw-p b7e93000 00:00 0 b7e97000-b7ebb000 r-xp 00000000 08:01 205303 /lib/i686/cmov/libm-2.7.so b7ebb000-b7ebd000 rw-p 00023000 08:01 205303 /lib/i686/cmov/libm-2.7.so b7ebd000-b7eec000 r-xp 00000000 08:01 224192 /usr/lib/libnetcdf.so.4.0.0 b7eec000-b7eed000 rw-p 0002f000 08:01 224192 /usr/lib/libnetcdf.so.4.0.0 b7eed000-b7eee000 rw-p b7eed000 00:00 0 b7f05000-b7f11000 r-xp 00000000 08:01 100546 /lib/libgcc_s.so.1 b7f11000-b7f12000 rw-p 0000b000 08:01 100546 /lib/libgcc_s.so.1 b7f12000-b7f15000 rw-p b7f12000 00:00 0 b7f15000-b7f2f000 r-xp 00000000 08:01 192949 /lib/ld-2.7.so b7f2f000-b7f31000 rw-p 0001a000 08:01 192949 /lib/ld-2.7.so bf916000-bf930000 rw-p bffe6000 00:00 0 [stack] ffffe000-fffff000 r-xp 00000000 00:00 0 [vdso] or segment violation. The same problem arises in other machine, both are Debian Lenny. GMT: 4.3.1 OS: Linux 2.6.24-1-686 #1 SMP Thu May 8 02:16:39 UTC 2008 i686 GNU/Linux CPU: Intel(R) Core(TM)2 Duo CPU E6750 < at > 2.66GHz Memory: 2GB gcc-4.3 4.3.1-2 libc6-i686 2.7-13 I've tried other projections: Azimutal Lambert: error Azimutal Equidistant: error Gnomonic: error Orthographic: error If I change slightly the region boundaries or lon0 or lat0 (0.01 arc second) pscoast runs smoothly. I recompiled pscoast with -g (compiler) and without -s (linker) when the same results (crash), when I ran it thru valgrind I got: Azimutal Lambert: nothing, pscoast doesn't crash Azimutal Equidistant: nothing, pscoast doesn't crash Gnomonic: ==26714== Invalid read of size 87 ==26714== at 0x80943AF: GMT_polygon_is_open (gmt_support.c:8343) ==26714== by 0x8075B63: GMT_rect_clip (gmt_map.c:4957) ==26714== by 0x8075D7A: GMT_clip_to_map (gmt_map.c:4760) ==26714== by 0x808F39B: GMT_prep_polygons (gmt_shore.c:695) ==26714== by 0x804BAA1: main (pscoast.c:741) ==26714== Address 0x4689600 is 8 bytes before a block of size 928 alloc'd ==26714== at 0x4021E22: calloc (vg_replace_malloc.c:397) ==26714== by 0x80A1DD1: GMT_memory_func (gmt_support.c:2440) ==26714== by 0x807578E: GMT_rect_clip (gmt_map.c:4922) ==26714== by 0x8075D7A: GMT_clip_to_map (gmt_map.c:4760) ==26714== by 0x808F39B: GMT_prep_polygons (gmt_shore.c:695) ==26714== by 0x804BAA1: main (pscoast.c:741) Orthograpic: nothing, pscoast doesn't crash Cheers, Edu.- Eduardo A. Suárez 2008-09-01T20:55:28 http://comments.gmane.org/gmane.comp.gis.gmt.user/9854 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 sqz2 2008-08-29T17:18:58 http://comments.gmane.org/gmane.comp.gis.gmt.user/9853 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 sqz2 2008-08-29T17:16:09 http://comments.gmane.org/gmane.comp.gis.gmt.user/9852 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 sqz2 2008-08-29T17:05:07 http://comments.gmane.org/gmane.comp.gis.gmt.user/9849 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 Paul Duncan 2008-08-29T11:28:10 http://comments.gmane.org/gmane.comp.gis.gmt.user/9846 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 <ericf< at >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 yacob sen 2008-08-28T23:14:35 http://comments.gmane.org/gmane.comp.gis.gmt.user/9844 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 Daniela Kuehn 2008-08-28T07:29:48 http://comments.gmane.org/gmane.comp.gis.gmt.user/9843 How can I obtain GMT4.0? To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Chandan Bora 2008-08-28T07:00:39 http://comments.gmane.org/gmane.comp.gis.gmt.user/9838 Hi I too am trying to see if I can import some geotiff images into GMT. Unfortunately my images are not axis aligned. So if I try to use gdal_translate, then I get the following error... I assume that I can probably work around the problem by reprojecting into a correctly aligned window, but I haven't managed to figure out the best way to do this yet. Does anyone have examples they can send me, or links to similar? Also gdal_transform tells me that it's only able to deal with a single band from the 3 colour bands in the geotiff. Is there an easy way to work around this and get the correctly located and coloured image into GMT? An example of the images I want to try and import (taken by autonomous robotic submarine) is available at. http://www-personal.acfr.usyd.edu.au/d.mercer/geotiff/GMTtest.tif Am I asking too much of GMT, would I be better sticking to one of the GIS packages? (GRASS / arcGIS etc.) Any helpful suggestions on any of the above gratefully received. Thanks Duncan 2008/8/27 Eduardo A. Suárez <esuarez< at >fcaglp.fcaglp.unlp.edu.ar>: Duncan Mercer 2008-08-28T01:39:28 http://comments.gmane.org/gmane.comp.gis.gmt.user/9836 Michael Kremer 2008-08-27T16:01:30 http://comments.gmane.org/gmane.comp.gis.gmt.user/9835 Dear Joaquim Luis I use GMT on Fedora8, WinXP and WinVista. I can compile GMT on linux but no success on windows. Would you please explain the steps required to compile GMT using MS VC++ 6? Thank you. To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Chandan Bora 2008-08-27T11:55:09 http://comments.gmane.org/gmane.comp.gis.gmt.user/9833 Hello, I have been looking at an image I have created as a Geotiff and need to transfer it to a grdfile in gmt to create a psimage, the geotiff has the header .twf with the header information, 28.50000000000000 0.00000000000000 0.00000000000000 -28.50000000000000 730084.50000000000000 -130530.00000000000000 is there anyway to convert this to a grdfile, Thanks Richard To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Richard Wall 2008-08-27T08:54:06 http://comments.gmane.org/gmane.comp.gis.gmt.user/9831 Ian Lilly* 2008-08-27T04:13:12 http://comments.gmane.org/gmane.comp.gis.gmt.user/9828 I was thinking of ways to produce more decorative maps... I like GMT's approach of overlaying PostScript files. I thought one approach would be to create a decorative overlay in a program like Inkscape, and superimpose it over GMT's output. After some googling, though, that doesn't seem entirely straightforward. Has anyone tried to overlay arbitrary PostScript files below? (I would like to avoid early rasterization if I can). Thanks, Adam To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Adam Baker 2008-08-26T17:57:49 http://comments.gmane.org/gmane.comp.gis.gmt.user/9826 Hello, I am currently trying to convert a Geotiff file from the program ENVI into a suitable grdfile for GMT. The file can be transferred to an ASCII format but the header is not read in GMT. Does anyone know of a way to convert these files into a grdfile? I have also added my correspondence with Paul below, Thankyou, Richard r.j.wall< at >dur.ac.uk You may want to repost your original question, modified by our correspondence, to the gmt-help< at >hawaii.edu list. Your posting seems to have been sent to the request address which is me only. You need to subscribe to post for spam reasons. If you end up hacking to make it work, you can see what xyz2grd expects for an ESRI grid by creating a dummy one: grdmath -R0/10/0/10 -I1 0 = dummy.nc grd2xyz -E dummy.nc > dummy.txt Cheers, Paul On Aug 26, 2008, at 6:52 AM, r.j.wall< at >durham.ac.uk wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Richard Wall 2008-08-26T17:27:18 http://comments.gmane.org/gmane.comp.gis.gmt.user/9823 Hello, sorry I am writing here but I did not succeed in the registration mailing list. What I want to know is How to convert tha bathymetry grids obtained from this site: ftp://topex.ucsd.edu/pub/srtm30_plus/data/ I couldn't find the way to get a GMT grid format. Could you help me on this? I am a very new GMT user. Thanks. Eugenia To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Maria Eugenia 2008-08-24T16:30:08 Search the mailing list at Gmane query http://search.gmane.org/?group=$group=gmane.comp.gis.gmt.user