gmane.comp.gis.gmt.user

Extract minimum depth from bathymetry in vicinity of track

Maarten Vanneste — 2008-12-04T23:35:57

Dear all, I am currently looking into iceberg ploughmarks and the variation of their depression across (a small part of) the shelf. I have some very-high-resolution bathymetric grids (down to 10 cm by 10 cm bin size). I also have coordinates of approximate positions (x,y,index) of the ploughmark. I have written a shell script which extracts bathymetric profiles (grdtrack) across the ploughmark (perpendicular) at regular intervals/distances, filters them (filter1d, remove outliers and smoothens), and analyses the characteristics of the ploughmark (defines rims, calculates area, etc.). The script works very satisfactorily, and the code (fortran) to define the ploughmark rims works reasonably well. However, from inspecting the data, I see that there are misfits. These errors are primarily caused by the fact that the profile in places falls outside the ploughmark, and thus one of the basic assumptions of the algorithm (1. a ploughmark is a depression and 2. the position falls within the ploughmark) becomes inval

Re: Bug when drawing ellipse

Paul Wessel — 2008-12-04T23:21:01

My educated guess would be Jan 1, 2009 Paul Wessel On Dec 4, 2008, at 1:10 PM, Roman Klaehne wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Bug when drawing ellipse

Roman Klaehne — 2008-12-04T23:10:32

Hi Remko, thats great! It works! Do you know when the next GMT version is built? Best regards /roman Remko Scharroo schrieb: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Bug when drawing ellipse

Remko Scharroo — 2008-12-04T22:38:28

Hi again, Indeed there was an error in the PostScript prologue in the definition of the circle (which is also used for the ellipse). Copy the attached file PSL_prologue.ps into your share/pslib directory, or update from CVS. If you do the latter, you need to update both pslib.c and PSL_prologue.ps. Regards, Remko To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu On 4 Dec 2008, at 16:48, Remko Scharroo wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Bug when drawing ellipse

Remko Scharroo — 2008-12-04T21:48:55

Hi Roman, Your input file is fine, the software is fine too. I discovered there is something wrong with the PostScript code that generates either the star or the ellipse symbol. Somewhere a "stroke" needs to be added. I'm working on it. Remko On 4 Dec 2008, at 15:28, Roman Klaehne wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Bug when drawing ellipse

Roman Klaehne — 2008-12-04T21:33:27

Hi Milton, the strange thing is, that with all other symbols except ellipse there are no problems with my input files and all is drawn as I want. Best regards /roman Milton P. Plasencia Linares schrieb: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Bug when drawing ellipse

Milton P. Plasencia Linares — 2008-12-04T21:32:10

Hi Roman, The -M flag plot multiple segment files and i think no multi symbols how you want to plot, for this reason I have divided in two commands psxy for plotting each symbol. Sorry I don't how to plot the symbols without the line using your nodes.data Saluti ********************************** 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 Roman Klaehne ZIB.DE>: Hi Milton, I think the -M flag says GMT that the input files have the format like my nodes.data file: > -G255/0/0 -W0.40p,0/0/0,solid 7.13 51.29 12.0 ap > -G255/0/0 -W0.40p,0/0/0,solid 10.53 48.22 0 4.0 8.0 ep I have to generate the GMT data automatically supporting different colors, sizes, shapes, etc. for nodes and links. S

Re: Bug when drawing ellipse

Roman Klaehne — 2008-12-04T20:28:49

Hi Milton, I think the -M flag says GMT that the input files have the format like my nodes.data file: 7.13 51.29 12.0 ap 10.53 48.22 0 4.0 8.0 ep I have to generate the GMT data automatically supporting different colors, sizes, shapes, etc. for nodes and links. So I have to use such input files. The -S option is not allowed to be added to the segment headers; the node shape information has to write behind the coordinates. Do you have any idea how to write your script so that it accepts input files with equal information link in my nodes.data file? I think the -S option will be a problem. Best regards /roman Milton P. Plasencia Linares schrieb: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Bug when drawing ellipse

Milton P. Plasencia Linares — 2008-12-04T20:06:58

Hi Roman, please, see the use the -M flag, perhaps here is the problem in connection with the file nodes.data, Try this script, plot the star and the ellipse without problems Cheers, #!/bin/csh # psbasemap -R5.55/13.72/47.59/52.83 -JM16.00c+ -Ba1f1:Longitude:/a1f1:Latitude:WeSn \ -X1.5 -Y6 -P -V -K >! ellipse.ps pscoast -R -JM -W -Di -Ia -Na -A20.00 -K -O -P -V >> ellipse.ps psxy -R -JM -W2 -Sa.3 -Gred -M -O -K -V -P <> ellipse.ps 7.13 51.29 END psxy -R -JM -W4 -Se.05 -Gred -M -O -V -P <> ellipse.ps 10.53 48.22 0 0.2 0.4 END \rm -f .gmt* ********************************** 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 Roman Klaehne ZIB.DE>: Hi, in the attachment you will find a sim

Re: Forcing surface to output flat surfaces?

Walter H. F. Smith — 2008-12-04T16:30:00

Gmt-team: In this instance, should 'surface' abort, or simply generate a trivial flat grid, bypassing all the spline iteration steps? Seems like the latter to me. w Walter H F Smith Chairman, GEBCO SCDB Geophysicist, Laboratory for Satellite Altimetry NOAA NESDIS code E/RA-31 1335 East West Hwy, room 5408 Silver Spring MD 20910-3226 tel 301-713-2857, extension 126 fax 301-713-3136 Walter.HF.Smith< at >noaa.gov On Dec 4, 2008, at 10:50 AM, Peter Schmidt wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Forcing surface to output flat surfaces?

Peter Schmidt — 2008-12-04T15:50:13

Yepp, I've tried it, got an error telling me something like: "surface: All data are defined on a flat surface, aborting" and an empty grd file I'll try out Joaquims last suggestion. Thanks for your quick response regP Walter H. F. Smith wrote:

Re: ps2raster EPS BoundingBox and PAPER_MEDIA

Remko Scharroo — 2008-12-04T15:35:52

Yes, that would explain it, since it would prevent GMT from writing out the BoundingBox at the end of the file. Remko To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

AW: [GMT-HELP] ps2raster EPS BoundingBox and PAPER_MEDIA

Andrea Heiker — 2008-12-04T15:28:03

Hi, I worked it out. I suppose the reason was, that I had the -K option in the last plot command. Andrea -----Ursprüngliche Nachricht----- Von: World-Wide GMT Usage and Help Mailing List [mailto:GMT-HELP< at >HAWAII.EDU] Im Auftrag von Joaquim Luis Gesendet: Donnerstag, 4. Dezember 2008 15:02 An: GMT-HELP< at >HAWAII.EDU Betreff: Re: [GMT-HELP] ps2raster EPS BoundingBox and PAPER_MEDIA Andrea, It works fine here, but I'm using the latest CVS and there were quite some improvements done to ps2raster. So, updating is probably what you need. Joaquim work. 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

Bug when drawing ellipse

Roman Klaehne — 2008-12-04T15:10:59

Hi, in the attachment you will find a simple script, which draws a star and a ellipse on a pscoast-map. The problem is, that, additionally to the symbols, a line connecting these symbols is drawn. I'am using a nodes.data file of the following content: 7.13 51.29 12.0 ap 10.53 48.22 0 4.0 8.0 ep I think the problem is the ellipse symbol. If I change it to a circle the "strange line" is not drawn. Also if I change the pen color of the star a little, the line is not drawn. I would say, that this behavior seems not to be correct. In the attachment you will find: gmtScript.sh (script for creating the ps-code) nodes.data (data file with the content above) out.ps (the created ps file) out.pdf (pdf created with ps2pdf) The command for creating the PS/PDF is: ./gmtScript.sh > out.ps; ps2pdf out.ps out.pdf ; acroread out.pdf Other information: GMT-Version: 4.3.1 OS: Linux OpenSuse Thanks in advance for a little help /Roman Klaehne Konrad Zuse Zentrum Berlin (ZIB) To unsubscribe, send the message "

Re: Forcing surface to output flat surfaces?

Joaquim Luis — 2008-12-04T14:44:45

Yes, thew grdmath solution implies that you must know in advance the cte value OK, if you data is on a flat surface, horizontal or not, than "triangulate" should be the right tool. Since internally it builds triangles, you will get a flat surface. Joaquim To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Forcing surface to output flat surfaces?

Walter H. F. Smith — 2008-12-04T14:39:15

Have you tried it? What happens if you feed surface such a data set? I think surface is supposed to remove a plane by least squares in the first place, and then interpolate the residual with [tensioned] splines, and then add the plane back in at the end. So maybe it works fine? w On Dec 4, 2008, at 9:31 AM, Peter Schmidt wrote: To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: Forcing surface to output flat surfaces?

Peter Schmidt — 2008-12-04T14:31:40

Hmmm, I guess by that you mean, writ e some extra code to check the data first, if on a flat surface use grdmath to construct a gridfile with that value. But what if the surface is flat but not at a constant value (i.e. inclined) regards Peter Joaquim Luis wrote:

Re: Forcing surface to output flat surfaces?

Joaquim Luis — 2008-12-04T14:04:02

Peter, Use grdmath grdmath -Rg -I1 10 = lixo.grd creates a global grid with a constant value of 10 Joaquim To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Re: ps2raster EPS BoundingBox and PAPER_MEDIA

Joaquim Luis — 2008-12-04T14:02:07

Andrea, It works fine here, but I'm using the latest CVS and there were quite some improvements done to ps2raster. So, updating is probably what you need. Joaquim To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu

Forcing surface to output flat surfaces?

Peter Schmidt — 2008-12-04T13:40:20

Hi Is there some option (or other way) to force surface to output a grd file if all defined coordinates are located on a flat surface? I want to use it to interpolate an uneven spaced set of data onto an even spaced set, and in some cases the input data all lie on a plane. (so since I'm lazy and GMT usually offers a quick solution it would be nice to be able to use surface in all cases) Regards Peter

ps2raster EPS BoundingBox and PAPER_MEDIA

Andrea Heiker — 2008-12-04T09:48:04

Hello, I want to convert my PS to EPS with e tight bounding box an it doesn't work. If I', not setting gmtset PAPER_MEDIA to A4+ the bounding box is the whole papersize A4 (but should be only my map+ legend) even if I'm using the -A -Te options for ps2raster. Code: gmtset ANNOT_FONT_SIZE_PRIMARY 10 LABEL_FONT_SIZE 10 PLOT_DEGREE_FORMAT ddd:mm:ssF . some code . ps2raster -A -Te my_map.ps But If I set the PAPER_MEDIA to A4+ no bounding box is created and my ps2ratser failed. Code gmtset PAPER_MEDIA A4+ ANNOT_FONT_SIZE_PRIMARY 10 LABEL_FONT_SIZE 10 PLOT_DEGREE_FORMAT ddd:mm:ssF . some code . ps2raster -A -Te my_map.ps The error message is: GMT FATAL ERROR: The file my_map.ps has no BoundingBox in the first or last 20 lines. Use the -A option. What am I doing wrong? Andrea To unsubscribe, send the message "signoff gmt-help" to listserv< at >hawaii.edu