Re: Extracting a line from a grid
Paul Wessel <[email protected]>
| Newsgroups | gmane.comp.gis.gmt.user |
|---|---|
| Message-ID | <[email protected]> |
Hm, I must be running a more recent version than you (I only get the 3 cols but if you are running the Nov 5 release it may have a bug we have since fixed); so you are also getting distances along the line. Until the next release (or update yours) you could use -o0,1,3 to get what you want. -p On Feb 25, 2014, at 10:00 AM, Shahar Barak <[email protected]> wrote: > Thanks, Paul! > I used: > gmt grdtrack -G${grd_file} -E${LON1}/${LAT1}/${LON2}/${LAT2}+i0.005d > but I expected an output of 3 columns: lon / lat / interpolated-data-point. > But I got 4 columns in the output. The first two seems like the lons and lats, not sure about the other two. > -117 32.25 0 0.0597692318261 > -116.995762712 32.2526483051 0.00499681204029 0.0604514658204 > -116.991525424 32.2552966102 0.00999362408057 0.0611412584514 > -116.987288136 32.2579449153 0.0149904361209 0.0618330594925 > -116.983050847 32.2605932203 0.0199872481611 0.0625284501071 > -116.978813559 32.2632415254 0.0249840602014 0.0632297094842 > -116.974576271 32.2658898305 0.0299808722417 0.0639385271851 > -116.970338983 32.2685381356 0.034977684282 0.0646567940722 > > gmt grdinfo ${grd_file} -C > tomo_slice.grd -120 -114 32 36 -0.333035200834 0.176299631596 0.005 0.005 1201 801 > > > > On Tue, Feb 25, 2014 at 11:18 AM, Paul Wessel <[email protected]> wrote: > Yes, see grdtrack -E > > On Feb 25, 2014, at 8:58 AM, Shahar Barak <[email protected]> wrote: > > > Hi, > > > > Is there a way to use GMT to extract data points from a grid along a line, given the 2 end points? > > > > Thanks, > > > > Shahar > > To unsubscribe, send the message "signoff gmt-help" to [email protected] > > To unsubscribe, send the message "signoff gmt-help" to [email protected] > > To unsubscribe, send the message "signoff gmt-help" to [email protected] To unsubscribe, send the message "signoff gmt-help" to [email protected]