How to crop a map in eps

Hello:

I want to crop a square from a circle map. In this case, it's a
gnomonic projection (images attached). By changing manually the
bounding box I managed to crop the map as I want it, but when I upload
it in AI or another drawing program, it doesn't respect the bounding
box, so the safest way is to crop the map itself.

How can I do it?

Thanks in advance for any hint.

-- A.

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

GMT Warning: Character encoding error

Hi,

I am a beginner in Linux and GMT. I am learning GMT from the GMT tutorial
document available along with the GMT software. I am trying to generate a
map using pscoast command, but I get the following error

santosh <at> santosh-laptop:~/tutorial$ pscoast -R-90/-70/0/20 -JM6i -B5g5 -P
-Gchocolate > map3.ps
GMT Warning: Selected character encoding does not have suitable degree
symbol - will use space instead
GMT Warning: Selected character encoding does not have minute symbol (single
quote) - will use space instead
GMT Warning: Selected character encoding does not have second symbol (double
quote) - will use space instead
Segmentation fault

It creates a ps document, but it doesn't open. Can anyone suggest me how to
fix this problem?
I would appreciate any help.

Thanks,
Santosh

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

grdfilter gaussian

Hello,

I've a grid file with 3 columns, long, lat and value. The area is
0/360/-70/70 and the resolution is 2 by 2 degrees. As it doesn't look very
pretty (see attachement) I'd like to smooth it by 4 or possibly 5 degrees,
using a Gaussian filter.

The command grdfilter diff_jan96_jul96.grd -Gfiltered_diff_96.grd -D4 (or 5)
-Fg4 -V produces the error message

grdfilter: Input nx,ny = (181 71), output nx,ny = (181 71), filter nx,ny =
(182 1)
grdfilter: Filter type is Gaussian.
grdfilter: Processing output line 70
grdfilter: Unable to estimate value at 4466 nodes, set to NaN.

Any help would be greatly appreciated.

Many thanks,

Rakia

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

GINA

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

map border problem (projection "Van der Grinten")

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

Netcdf with GMT

grdraster.info for GLOBE

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

ps* -B:."Title" without frame?

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
--------------------------------------------------------------------

Question about pstext -Z+ -W

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 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
--------------------------------------------------------------------

Re: Large PS files conversion ?

Walter,

please, correct me if I am wrong.
Walter H. F. Smith schrieb am 30.09.2008 16:41:
> Stephan,
> 
> Are we misunderstanding each other?
I think so - but I fear the problem is between my ears.

> ps2raster 4.3.2 (CVS Jun 25 2008 11:06:24) correctly computes
> the size of the necessary bounding box, even if it was wrong
> in the original postscript file.  The png and pdf it creates
> are also OK.
> 
> What did you need to fix in CVS?

As far as I see:
* your system (GMT, gs, ...) works correctly - the boundingbox is
  tightly adjusted to te large plot;
* my system crops the plot to a4 size;
* I mentioned that I made an CVS update (of my system) to be sure
  to have the most recent version;
* the bounding box reported by ps2raster on my system is wrong
  although it has the right values in the original eps file;
All this leads me to the assumption, that gs on my system makes
something strange...

ps2raster: Processing testa0.eps: Find HiResBoundingBox [45.9 51 595
842]... Convert to PNG... Done.

ps2raster: Processing testa0.eps: Find HiResBoundingBox [45.9 51 595
842]... Convert to PDF... Done.

grep -i bound testa0.eps

%%BoundingBox: (atend)
%%HiResBoundingBox: (atend)
%%BoundingBox: 0 0 3262 2128
%%HiResBoundingBox: 0 0 3262.11 2128.25

Confused,

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
--------------------------------------------------------------------

path installation problem