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/10095 Rakia Meister 2008-10-07T13:45:02 http://comments.gmane.org/gmane.comp.gis.gmt.user/10087 I just came across a unified global topography/bathymetry gridded data set with 30 arc-sec resolution for topography and 2 arc-minute resolution for bathymetry: <http://www.gina.alaska.edu/page.xml?group=data&page=griddata> Evidently, one should refer to it as the GINA data set. This seems to be an improvement over the GLOBE and ETOPO2 data sets, which give equivalent resolution for topography and bathymetry, respectively. I intended to use GLOBE on land and ETOPO2 in the oceans, as backgrounds for seismicity plots, but mixed land/ocean regions are a problem. GINA is said to be easily used with GMT (already in .grd format). Has anyone tried it out in GMT? Is there any reason not to switch to this as the best general-purpose global gridded relief data set? Eric To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Eric Bergman 2008-10-03T23:13:10 http://comments.gmane.org/gmane.comp.gis.gmt.user/10086 Hello. I'm new to the GMT. I am reading the documentation and I have a little problem. I'm playing with the manual examples of psbasemap. I'm tring "Van der Grinten" projection. At first, I wanted to find a way to generate annotations on the two axes, not only the latitude. For this I modified the parameter-R: psbasemap -R0/360/-80/80 -JV180/7i -Bg30a30/g30a30 Sinusoidal projection work fine too: psbasemap -R0/360/-80/80 -Ji0/0.02i -Bg30a30/g30a30 Thanks for your attention. José María Michia To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu José María Michia 2008-10-03T19:49:20 http://comments.gmane.org/gmane.comp.gis.gmt.user/10083 Hi, Can GMT be used to create plot using netcdf files (*.cdf, *.nc) ? Regards Saji. P. K. Add more friends to your messenger and enjoy! Go to http://messenger.yahoo.com/invite/ To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Saji Kuttan 2008-10-03T08:23:57 http://comments.gmane.org/gmane.comp.gis.gmt.user/10081 Does anyone have the correct entry for the GLOBE DEM tiles for grdraster.info? Thanks, Eric To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Eric Bergman 2008-10-02T23:20:15 http://comments.gmane.org/gmane.comp.gis.gmt.user/10077 -----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi all, maybe I have missed something obvious. I would like to plot some text or symbols (pstext/psxy) and use -B to put a title above. But I don't want to have _any_ frame. As far as I see, e.g. ps[xy|test] -JX5i -R0/5/0/5 -B:."Title": /dev/null >! test.eps always draws at least the line of the frame. Even --FRAME_PEN=0p doesn't work. Is this wanted or would it be an improvement e.g. to be able to set -B-:."Title": to force the px*-commands not to draw the frame? Regards, Stephan - -- - -------------------------------------------------------------------- Stephan Eickschen Fon: +49 (0)228 933 999 70 Augustastrasse 89 Fax: +49 (0)228 933 999 71 53173 Bonn Cell: +49 (0)170 54 781 60 Germany Email: eickschen.haenel< at >t-online.de http://www.eiskultur.de - ------------------------------------------------------------- Stephan Eickschen 2008-10-02T18:43:18 http://comments.gmane.org/gmane.comp.gis.gmt.user/10076 -----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi all, I am wondering if I misunderstand something. I would like to put some text on a 3D map using pstext and have an outline around each text. I use the following commands: ==== snip ============================================= #!/bin/csh -f rm -f .gmt* set RJE="-JX5 -JZ2 -R0/10/0/10/0/1.1 -E60/30" set eps=test2.eps pstext $RJE -Z+ -W255o -B:"Band 1":/:"Band 2":/a0.5:"P":wsNEZ -V -P << EOF >! $eps 3 1 0.00 12 120 0 MC Class 1 4 7 0.25 12 120 0 MC Class 2 8 5 0.50 12 120 0 MC Class 3 EOF ps2raster -A -Tf $eps ==== snap ============================================= The test is set a the right position, but the outline is misplaced. It seems as if the z value is doubled: z=0 everything is fine ("Class 1") but for z=0.25 and z=0.5 the offset of the box from the text is obvious and is doubled as z itself is doubled... (s. attachment) Any comments? Regards, Stephan - -- - -------------------------------------------------------------------- Stephan Eicksc Stephan Eickschen 2008-10-02T16:35:35 http://comments.gmane.org/gmane.comp.gis.gmt.user/10037 Dear all I'm tring to install GMT for the first time in a laptop machine under linux ubuntu. I have used the install_gmt automatic script. The installation did not returns any error messages. Anyway if I try to run gmt programs, shell returns this path message error: bash: path: command not found I don't know why and I don't know what I have to do.. I think it is a "stupid" problem but I am just coming to approach unix. thanks to everybody help me Valerio Pasini To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu valerio pasini 2008-09-30T11:07:46 http://comments.gmane.org/gmane.comp.gis.gmt.user/10030 Hi all, I've used GMT 3.4.5-2 for years without trouble on a Mac (installed using Fink). I mainly draw visibility circles centered around points on Earth maps, see e.g. the attached "gmtcircle-Tiger.gif" file. When I upgraded my MacOS from 10.3 to 10.4, GMT continued to work fine. I have just upgraded to MacOS 10.5 (Leopard), and simply copied my /sw directory to my new Leopard volume. GMT (therefore still 3.4.5-2) appears to work fine, except the circles drawn with gmtcircle who look quite weird as seen on the attached "gmtcircle-Leopard.gif" file. I'm curious to know where this strange result comes from, and why gmtcircle behaves differently on Leopard than on Tiger. I searched the archive and found a 2006 message from Paul Wessel saying that gmtcircle is no longer required and can be replaced with psxy -SE, so I tried this option. It works the same way on Tiger and Leopard, but not exactly as I expected. A small part of the circle is missing, see the attached "psxySE.gif" file. Thanks in ad Herve Fagard 2008-09-29T15:23:48 http://comments.gmane.org/gmane.comp.gis.gmt.user/10023 Hi All, Does anyone know if it's possible to plot a degree symbol as part of an axis label - i.e. with the -B command, as one would with a units label? Thanks, Gavin ----------------------- Gavin P Hayes Post-doc Researcher USGS National Earthquake Information Center (Contracted by Synergetics, Inc.) ghayes< at >usgs.gov 303-273-8421 To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Gavin Hayes 2008-09-24T21:46:20 http://comments.gmane.org/gmane.comp.gis.gmt.user/10017 Hello, I'm trying to create a huge topographic map of europe using srtm3 data. Dowload and conversion seemed to work flawlessly; now i'm trying to paste 41 rows of 1°lat and 98°long together in one single grd grid, which will be used to create the map. Unfortunatly each row is 540 mb big, and after having pasted the 5 first rows, grdpaste dies badly, with the following error message: GMT fatal error: grdpaste could not allocate memory [3.2 Gb] I have an old amd athlon xp 1,6 ghz with 1gb RAM, and 2Gb SWAP. Do I need the full 41x540Mb (which will be the size of the final grid) amount of swap ? Is there no way to tell grdpaste to directly write the new row a the end of the final grid file, on the hdd, without mapping the whole thing into memory ? Thank you for your help, Felix To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu felix.schalck 2008-09-24T13:15:33 http://comments.gmane.org/gmane.comp.gis.gmt.user/10013 I'm trying to use default axis lengths in landscape view, but I keep getting a syntax error. The same problem doesn't occur with portrait view: psbasemap -JXv -R0/50/0/50 -B10/10 -P > test.ps runs cleanly and produces the expected plot while psbasemap -JXh -R0/50/0/50 -B10/10 > test.ps gives psbasemap: Warning: not a valid number and may not be decoded properly. psbasemap: GMT SYNTAX ERROR -J option. Correct syntax: -Jx<x-scale|width>[d|l|p<power>|t|T][/<y-scale|height>[d|l|p<power>|t|T]], scale in cm/units -Jz<z-scale>[l|p<power>], scale in cm/units Use / to specify separate x/y scaling (e.g., -Jx0.5/0.3.). Not allowed with 1:xxxxx Use -JX (and/or -JZ) to give axes lengths rather than scales I've tried it on a Sun with GMT4.3.1 and a Mac with GMT4.2.1 and get the same error on both. Can anyone see what I'm doing wrong? Thanks, Lara To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Lara Kalnins 2008-09-24T10:31:23 http://comments.gmane.org/gmane.comp.gis.gmt.user/10010 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 Mark Collier 2008-09-24T01:33:30 http://comments.gmane.org/gmane.comp.gis.gmt.user/10007 Several versions of GMT ago the following shell script worked fine. It used to draw an azimuthal equidistant map out to 180 degrees centered on user specified long and lat with an azimuthal grid on the center of the map. It also drew 2 points on the map. Now it complains about map limits from the pscoast and psxy call. Any ideas on what has changed and how to update the script? Thanks. Bob > stn_az_map -89.932 35.123 MEM pscoast: North is outside -90 to +90 degree range psxy: North is outside -90 to +90 degree range It should draw a map centered on Memphis and label it. #!/bin/sh ROOT=$HOME # this shell script does not have to be edited to run WORLDCOAST=0/360/0/180 WORLDGRID=0/360/0/90 RED=250/50/50 BLUE=50/50/255 GREEN=50/255/50 MOREPS=-K ADDPS=-O CONTINUEPS="-K -O" FILL=200 SCALE=1.75 XOFFSET=0.75 YOFFSET=1.5 GRIDCNTR=180/90/7/90 OUTPUTFILE=$0_$3.ps MAP="-R$WORLDCOAST -Je$1/$2/$SCALE/180" rm $OUTPUTFILE pscoast $MAP -S200/200/255 -G200/250/200 -W1 -Dc -P $MOREPS -X$XOFFSET -Y$YOFFSET > $OUTPUT Smalley 2008-09-22T19:02:40 http://comments.gmane.org/gmane.comp.gis.gmt.user/10005 To those who work with the CVS version of GMT: As of Saturday Sept 20 the GMT CVS server, pohaku.soest.hawaii.edu, has been fully restored. No information was lost. Cheers, Paul To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Paul Wessel 2008-09-21T20:45:42 http://comments.gmane.org/gmane.comp.gis.gmt.user/10002 Hello Paul , I am a newie using GMT, I just installed it and tried to create a map.ps , but I got this message "map.ps permission denied" MY operative system is ubuntu hardy. Do you know what have I done wrong? thanks for helping. Rengifo To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Rengifo Ortega 2008-09-18T20:23:14 http://comments.gmane.org/gmane.comp.gis.gmt.user/10001 Hi, I recently tried my hand at making a custom symbol for use with psxy. This was working fine under GMT 4.0, but, when I started using the same symbol definition with GMT 4.3.1 running under CentOS 5.2, I started to see spurious lines in my figure. These lines come and go and change shape depending on the size I specify for the symbol. I have attached my symbol definition file, aquarius.def , along with the images output by the following two command pipelines. (I use GhostScript to rasterize the PostScript from psxy to PNG format.) echo "2 2" \ | psxy -R0/4/0/4 -JX4 -P -X0 -Y0 -Skaquarius/2.15 \ | gs -sDEVICE=pngalpha -sOutputFile=- -q -r100 -g400x400 - \ > symbol_size.2.15.png echo "2 2" \ | psxy -R0/4/0/4 -JX4 -P -X0 -Y0 -Skaquarius/2.14 \ | gs -sDEVICE=pngalpha -sOutputFile=- -q -r100 -g400x400 - \ > symbol_size.2.14.png The only difference between the two commands is a hundredth of an inch in the symbol size specification; nevertheless, the first command produces the expected image while the s Norman Kuring 2008-09-19T17:38:38 http://comments.gmane.org/gmane.comp.gis.gmt.user/9994 Hi! I'm trying to generate a grid of a xyz data of 1800000 points and, either with the nearneighbor and surface technique, I get this error: *Fatal error: nearneighbor/surface could not allocate memory [661,71 Mb, n_items = 173464325 / 3,12 Gb, n_items = -944723638]* I'm guessing that with the triangulation technique I'll get the same error. I tried at another PC (with fewer RAM and hard drive memory) and it worked!, so I'm thinking maybe the problem lies in some missing .DLL file... or files of the GMT latest version. My OS is Windows XP Media Center Edition 2002 32-Bit with SP 3, the GMT version I'm using is 4.3.1. and the command lines that fail with error are: *nearneighbor datostodos_topografia.dat -GTopografiaGMT_nearneighbor.grd -I2500e/2500e -N4 -R-2567500/2000000/1105000/5900000r -S10000 -V* *surface datostodos_topografia.dat -GTopografiaGMT_surface.grd -I2500e/2500e -R-2567500/2000000/1105000/5900000r -V* I'd really appreciate your help. Thank you very much in advance. Angie Angie Bermùdez 2008-09-18T13:09:31 http://comments.gmane.org/gmane.comp.gis.gmt.user/9992 Hi all, The grdreformat command gives the following error message: grdreformat: GMT ERROR: Cannot read file /usr/local/GMT/share/gmt_formats.conf and in fact there is no gmt_formats.conf file in that location nor is it in any other place. I am using GMT 4.3.1 downloaded today. Is it a bug or have I missed something? Thanks in advance, Hernán Hernán De Angelis 2008-09-18T09:40:59 http://comments.gmane.org/gmane.comp.gis.gmt.user/9982 Hi All, I've been trying to merge 2 SRTM-4 tiles together using grdpaste. Given the information from grdinfo, the tiles have the same resolution and a common edge but I'm getting the error: grdpaste: Grids do not share a common edge! The files are apparently fine when used alone (function tested grdcut, grdimage) but they can't be merge together for some reason. Here are the file information: File1 Z_43_17.grd: Title: z Z_43_17.grd: Command: xyz2grd Z_43_17.ASC -GZ_43_17.grd -E -F Z_43_17.grd: Remark: Z_43_17.grd: Normal node registration used Z_43_17.grd: grdfile format: nf (# 18) Z_43_17.grd: x_min: 30 x_max: 34.9992 x_inc: 0.000833333 name: x nx: 6000 Z_43_17.grd: y_min: -25 y_max: -20.0008 y_inc: 0.000833333 name: y ny: 6000 Z_43_17.grd: z_min: -251 z_max: 2156 name: z Z_43_17.grd: scale_factor: 1 add_offset: 0 File2 Z_42_17.grd: Title: z Z_42_17.grd: Command: xyz2grd Z_42_17.ASC -GZ_42_17.grd -E -F Z_42_17.grd: Remark: Z_42_17.grd: Normal node registration used Z_42_17.grd: grdfile format: nf (# Anne-Sophie Meriaux 2008-09-17T21:08:15 http://comments.gmane.org/gmane.comp.gis.gmt.user/9978 Hi folks, I've about 800 swath bathymetry files that look good when we grid them using surface. Unfortunately the surface then extends well beyond the extents of the non-rectangular weirdly-shaped data. I am trying to think of a way to automatically clip the surface files to the data extent once each "chip" is finished. Any suggestions would be very welcome, Cheers, Mike Postdoctoral Researcher Byrd Polar Research Center Ohio State University www.polenet.org To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu Mike Willis 2008-09-17T20:01:01 Search the mailing list at Gmane query http://search.gmane.org/?group=$group=gmane.comp.gis.gmt.user