Re: Extracting a line from a grid
Paul Wessel <[email protected]>
| Newsgroups | gmane.comp.gis.gmt.user |
|---|---|
| Message-ID | <[email protected]> |
You are running GMT4. There is no such grdtrack option in in GMT 4. The email train explains a workaround for how that sort of thing was done in GMT4. Get GMT5 if you wish to use GMT5 options. -p On Feb 25, 2014, at 7:08 PM, Sanjay <[email protected]> wrote: > Now i used > grdtrack -bou.grd -p89/20/98/30+0.004d > but still error > -p is unrecognised option > > > Regards > sanjay > > > On Wednesday, February 26, 2014 10:32 AM, Paul Wessel <[email protected]> wrote: > And if you read the email train you would see grdtrack -E is a GMT 5 feature. > GMT 4 scheme is also discussed in the previous emails. > -p > > On Feb 25, 2014, at 6:22 PM, Sanjay <[email protected]> wrote: > >> When i tried to extract the profile from grdtrack >> error message came with E is unrecognized option. >> i am using GMT 4.5.7 >> Regards >> sanjay >> >> >> >> >> On Wednesday, February 26, 2014 1:42 AM, Paul Wessel <[email protected]> wrote: >> 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] >> >> >> 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]