|
»
grdrotater problem
|
View:
New views
7 Messages
—
Rating Filter:
Alert me
|
 |
grdrotater problem
by Andy_Alvey
() ::
Rate this Message:
I am new to using GMT and I was told that I could preform plate tectonic reconstructions of gridded data using the grdrotater function. I have tried using this and the program appears to rotate my map (using a clipping polygon) into the correct position. my problem is that the data inside the rotated grid is being destroyed to give a series of horizontal lines. (see attached .ps file). is this a problem with the input data or something that i am doing.
the script I am using is: grdrotater ct.grd -T-34.96/53.06/49.45 -V -Fsa_cob.xy -Gctrot.grd >rotct_path.xy grdimage ct.grd -R-90/40/-65/60 -Cct.cpt -JM6i -B20g10 -P -K -Q > Plotrot.ps grdimage ctrot.grd -R-90/40/-65/60 -Cct.cpt -JM6i -Q -B20g10 -P -K -O >> Plotrot.ps pscoast -R-90/40/-65/60 -JM6i -P -O -K -Di -W >> plotrot.ps any info is greatly appreciated cheers Andy Plotrot.ps |
|
|
Minimum region of interest?
by Peter Schmidt-3
::
Rate this Message:
Hmmm
Is there any neat way to find the minimum region of interest (i.e. to append to the -R option) that includes a known data set for a known projection? Or do I need to fiddle with this myself for every new data set I wish to plot. If not existing yet maybe this is something that can be added to the wishlist, i.e. given a data set data.xy, something along the lines: >>> region data.xy -Jparameters xmin xmax ymin ymax At least I assume that such a command could be written since the -N option exist (and -S option for mapproject) , which tricks me to believe that GMT somehow can decide if a point falls outside the map-boundaries for a given projection, hence given this it should be possible to find the minimum map-boundaries that encloses the data-set. Has anyone looked into this, if not and there will be no attempts, any pointers on where to start looking if I decide to try it myself would be highly appreciated as I'm not fluent in C (yet ;-) regards Peter (and yes I know that such a little tool is not really along the lines of GMT's philosophy, but I think it would be very useful, especially for newbies) -- ******************************************************************************** Peter Schmidt Tel: +46-18-4717104 PhD Student Mobile: +46-73-3190975 Dept. of Earth Sciences e-mail: peter.schmidt@... Geophysics http://www.geofys.uu.se/?q=view/home/75 Uppsala University Villavagen 16 SE-75236 Uppsala ******************************************************************************** To unsubscribe, send the message "signoff gmt-help" to listserv@... |
|
|
Re: Minimum region of interest?
by Walter H. F. Smith
::
Rate this Message:
Have you looked at minmax and its options?
On Nov 10, 2009, at 12:29 PM, Peter Schmidt wrote: > Hmmm > > Is there any neat way to find the minimum region of interest (i.e. to > append to the -R option) that includes a known data set for a known > projection? Or do I need to fiddle with this myself for every new data > set I wish to plot. If not existing yet maybe this is something > that can > be added to the wishlist, i.e. given a data set data.xy, something > along > the lines: > >>>> region data.xy -Jparameters > xmin xmax ymin ymax > > At least I assume that such a command could be written since the -N > option exist (and -S option for mapproject) , which tricks me to > believe > that GMT somehow can decide if a point falls outside the map- > boundaries > for a given projection, hence given this it should be possible to find > the minimum map-boundaries that encloses the data-set. Has anyone > looked > into this, if not and there will be no attempts, any pointers on where > to start looking if I decide to try it myself would be highly > appreciated as I'm not fluent in C (yet ;-) > > > regards Peter > > > (and yes I know that such a little tool is not really along the > lines of > GMT's philosophy, but I think it would be very useful, especially for > newbies) > > -- > ********************************************************************** > ********** > Peter Schmidt Tel: +46-18-4717104 > PhD Student Mobile: +46-73-3190975 > Dept. of Earth Sciences e-mail: peter.schmidt@... > Geophysics http://www.geofys.uu.se/?q=view/home/75 > Uppsala University > Villavagen 16 > SE-75236 Uppsala > ********************************************************************** > ********** > > To unsubscribe, send the message "signoff gmt-help" to > listserv@... To unsubscribe, send the message "signoff gmt-help" to listserv@... |
|
|
Re: Minimum region of interest?
by Peter Schmidt-3
::
Rate this Message:
Yes, that will certainly work for Cartesian coordinates (Jx), not really
sure if this is the case for geographic projections. Hmmm , maybe I should try it to see what happens, I'll give it a go... BTW. thanks for the ninja fast reply, must have been less then 5 secs, new record :-D regP Walter H. F. Smith wrote: > Have you looked at minmax and its options? > > > > > On Nov 10, 2009, at 12:29 PM, Peter Schmidt wrote: > >> Hmmm >> >> Is there any neat way to find the minimum region of interest (i.e. to >> append to the -R option) that includes a known data set for a known >> projection? Or do I need to fiddle with this myself for every new data >> set I wish to plot. If not existing yet maybe this is something that can >> be added to the wishlist, i.e. given a data set data.xy, something along >> the lines: >> >>>>> region data.xy -Jparameters >> xmin xmax ymin ymax >> >> At least I assume that such a command could be written since the -N >> option exist (and -S option for mapproject) , which tricks me to believe >> that GMT somehow can decide if a point falls outside the map-boundaries >> for a given projection, hence given this it should be possible to find >> the minimum map-boundaries that encloses the data-set. Has anyone looked >> into this, if not and there will be no attempts, any pointers on where >> to start looking if I decide to try it myself would be highly >> appreciated as I'm not fluent in C (yet ;-) >> >> >> regards Peter >> >> >> (and yes I know that such a little tool is not really along the lines of >> GMT's philosophy, but I think it would be very useful, especially for >> newbies) >> >> -- >> ******************************************************************************** >> >> Peter Schmidt Tel: +46-18-4717104 >> PhD Student Mobile: +46-73-3190975 >> Dept. of Earth Sciences e-mail: peter.schmidt@... >> Geophysics http://www.geofys.uu.se/?q=view/home/75 >> Uppsala University >> Villavagen 16 >> SE-75236 Uppsala >> ******************************************************************************** >> >> >> To unsubscribe, send the message "signoff gmt-help" to >> listserv@... > > To unsubscribe, send the message "signoff gmt-help" to > listserv@... > -- ******************************************************************************** Peter Schmidt Tel: +46-18-4717104 PhD Student Mobile: +46-73-3190975 Dept. of Earth Sciences e-mail: peter.schmidt@... Geophysics http://www.geofys.uu.se/?q=view/home/75 Uppsala University Villavagen 16 SE-75236 Uppsala ******************************************************************************** To unsubscribe, send the message "signoff gmt-help" to listserv@... |
|
|
Re: Minimum region of interest?
by Walter H. F. Smith
::
Rate this Message:
If you want something in projected coordinates, it may work to choose
a -R that is sure to be large enough, use mapproject with that -R and your desired -J, run the result through minmax, and then put the output of that through the inverse projection. In this way you can try to narrow down to a better choice of -R. w On Nov 10, 2009, at 12:43 PM, Peter Schmidt wrote: > Yes, that will certainly work for Cartesian coordinates (Jx), not > really > sure if this is the case for geographic projections. Hmmm , maybe I > should try it to see what happens, I'll give it a go... > > BTW. thanks for the ninja fast reply, must have been less then 5 secs, > new record :-D > > regP > > Walter H. F. Smith wrote: >> Have you looked at minmax and its options? >> >> >> >> >> On Nov 10, 2009, at 12:29 PM, Peter Schmidt wrote: >> >>> Hmmm >>> >>> Is there any neat way to find the minimum region of interest >>> (i.e. to >>> append to the -R option) that includes a known data set for a known >>> projection? Or do I need to fiddle with this myself for every new >>> data >>> set I wish to plot. If not existing yet maybe this is something >>> that can >>> be added to the wishlist, i.e. given a data set data.xy, >>> something along >>> the lines: >>> >>>>>> region data.xy -Jparameters >>> xmin xmax ymin ymax >>> >>> At least I assume that such a command could be written since the -N >>> option exist (and -S option for mapproject) , which tricks me to >>> believe >>> that GMT somehow can decide if a point falls outside the map- >>> boundaries >>> for a given projection, hence given this it should be possible to >>> find >>> the minimum map-boundaries that encloses the data-set. Has anyone >>> looked >>> into this, if not and there will be no attempts, any pointers on >>> where >>> to start looking if I decide to try it myself would be highly >>> appreciated as I'm not fluent in C (yet ;-) >>> >>> >>> regards Peter >>> >>> >>> (and yes I know that such a little tool is not really along the >>> lines of >>> GMT's philosophy, but I think it would be very useful, especially >>> for >>> newbies) >>> >>> -- >>> ******************************************************************** >>> ************ >>> >>> Peter Schmidt Tel: +46-18-4717104 >>> PhD Student Mobile: +46-73-3190975 >>> Dept. of Earth Sciences e-mail: peter.schmidt@... >>> Geophysics http://www.geofys.uu.se/?q=view/ >>> home/75 >>> Uppsala University >>> Villavagen 16 >>> SE-75236 Uppsala >>> ******************************************************************** >>> ************ >>> >>> >>> To unsubscribe, send the message "signoff gmt-help" to >>> listserv@... >> >> To unsubscribe, send the message "signoff gmt-help" to >> listserv@... >> > > > -- > ********************************************************************** > ********** > Peter Schmidt Tel: +46-18-4717104 > PhD Student Mobile: +46-73-3190975 > Dept. of Earth Sciences e-mail: peter.schmidt@... > Geophysics http://www.geofys.uu.se/?q=view/home/75 > Uppsala University > Villavagen 16 > SE-75236 Uppsala > ********************************************************************** > ********** > > To unsubscribe, send the message "signoff gmt-help" to > listserv@... To unsubscribe, send the message "signoff gmt-help" to listserv@... |
|
|
Re: Minimum region of interest?
by Mike Willis
::
Rate this Message:
In a similar vein, a tool that gives just the R values for where two
(or more) grids intersect would be useful. On Tue, Nov 10, 2009 at 12:45 PM, Walter H. F. Smith <Walter.HF.Smith@...> wrote: > If you want something in projected coordinates, it may work to choose a -R > that is sure to be large enough, use mapproject with that -R and your > desired -J, run the result through minmax, and then put the output of that > through the inverse projection. In this way you can try to narrow down to a > better choice of -R. > > w > > > > On Nov 10, 2009, at 12:43 PM, Peter Schmidt wrote: > >> Yes, that will certainly work for Cartesian coordinates (Jx), not really >> sure if this is the case for geographic projections. Hmmm , maybe I >> should try it to see what happens, I'll give it a go... >> >> BTW. thanks for the ninja fast reply, must have been less then 5 secs, >> new record :-D >> >> regP >> >> Walter H. F. Smith wrote: >>> >>> Have you looked at minmax and its options? >>> >>> >>> >>> >>> On Nov 10, 2009, at 12:29 PM, Peter Schmidt wrote: >>> >>>> Hmmm >>>> >>>> Is there any neat way to find the minimum region of interest (i.e. to >>>> append to the -R option) that includes a known data set for a known >>>> projection? Or do I need to fiddle with this myself for every new data >>>> set I wish to plot. If not existing yet maybe this is something that can >>>> be added to the wishlist, i.e. given a data set data.xy, something along >>>> the lines: >>>> >>>>>>> region data.xy -Jparameters >>>> >>>> xmin xmax ymin ymax >>>> >>>> At least I assume that such a command could be written since the -N >>>> option exist (and -S option for mapproject) , which tricks me to believe >>>> that GMT somehow can decide if a point falls outside the map-boundaries >>>> for a given projection, hence given this it should be possible to find >>>> the minimum map-boundaries that encloses the data-set. Has anyone looked >>>> into this, if not and there will be no attempts, any pointers on where >>>> to start looking if I decide to try it myself would be highly >>>> appreciated as I'm not fluent in C (yet ;-) >>>> >>>> >>>> regards Peter >>>> >>>> >>>> (and yes I know that such a little tool is not really along the lines of >>>> GMT's philosophy, but I think it would be very useful, especially for >>>> newbies) >>>> >>>> -- >>>> >>>> ******************************************************************************** >>>> >>>> Peter Schmidt Tel: +46-18-4717104 >>>> PhD Student Mobile: +46-73-3190975 >>>> Dept. of Earth Sciences e-mail: peter.schmidt@... >>>> Geophysics http://www.geofys.uu.se/?q=view/home/75 >>>> Uppsala University >>>> Villavagen 16 >>>> SE-75236 Uppsala >>>> >>>> ******************************************************************************** >>>> >>>> >>>> To unsubscribe, send the message "signoff gmt-help" to >>>> listserv@... >>> >>> To unsubscribe, send the message "signoff gmt-help" to >>> listserv@... >>> >> >> >> -- >> >> ******************************************************************************** >> Peter Schmidt Tel: +46-18-4717104 >> PhD Student Mobile: +46-73-3190975 >> Dept. of Earth Sciences e-mail: peter.schmidt@... >> Geophysics http://www.geofys.uu.se/?q=view/home/75 >> Uppsala University >> Villavagen 16 >> SE-75236 Uppsala >> >> ******************************************************************************** >> >> To unsubscribe, send the message "signoff gmt-help" to listserv@... > > To unsubscribe, send the message "signoff gmt-help" to listserv@... > To unsubscribe, send the message "signoff gmt-help" to listserv@... |
|
|
Re: Minimum region of interest?
by Peter Schmidt-3
::
Rate this Message:
That might actually work, thanks again
\p Walter H. F. Smith wrote: > If you want something in projected coordinates, it may work to choose > a -R that is sure to be large enough, use mapproject with that -R and > your desired -J, run the result through minmax, and then put the > output of that through the inverse projection. In this way you can > try to narrow down to a better choice of -R. > > w > > > > On Nov 10, 2009, at 12:43 PM, Peter Schmidt wrote: > >> Yes, that will certainly work for Cartesian coordinates (Jx), not really >> sure if this is the case for geographic projections. Hmmm , maybe I >> should try it to see what happens, I'll give it a go... >> >> BTW. thanks for the ninja fast reply, must have been less then 5 secs, >> new record :-D >> >> regP >> >> Walter H. F. Smith wrote: >>> Have you looked at minmax and its options? >>> >>> >>> >>> >>> On Nov 10, 2009, at 12:29 PM, Peter Schmidt wrote: >>> >>>> Hmmm >>>> >>>> Is there any neat way to find the minimum region of interest (i.e. to >>>> append to the -R option) that includes a known data set for a known >>>> projection? Or do I need to fiddle with this myself for every new data >>>> set I wish to plot. If not existing yet maybe this is something >>>> that can >>>> be added to the wishlist, i.e. given a data set data.xy, something >>>> along >>>> the lines: >>>> >>>>>>> region data.xy -Jparameters >>>> xmin xmax ymin ymax >>>> >>>> At least I assume that such a command could be written since the -N >>>> option exist (and -S option for mapproject) , which tricks me to >>>> believe >>>> that GMT somehow can decide if a point falls outside the >>>> map-boundaries >>>> for a given projection, hence given this it should be possible to find >>>> the minimum map-boundaries that encloses the data-set. Has anyone >>>> looked >>>> into this, if not and there will be no attempts, any pointers on where >>>> to start looking if I decide to try it myself would be highly >>>> appreciated as I'm not fluent in C (yet ;-) >>>> >>>> >>>> regards Peter >>>> >>>> >>>> (and yes I know that such a little tool is not really along the >>>> lines of >>>> GMT's philosophy, but I think it would be very useful, especially for >>>> newbies) >>>> >>>> -- >>>> ******************************************************************************** >>>> >>>> >>>> Peter Schmidt Tel: +46-18-4717104 >>>> PhD Student Mobile: +46-73-3190975 >>>> Dept. of Earth Sciences e-mail: peter.schmidt@... >>>> Geophysics http://www.geofys.uu.se/?q=view/home/75 >>>> Uppsala University >>>> Villavagen 16 >>>> SE-75236 Uppsala >>>> ******************************************************************************** >>>> >>>> >>>> >>>> To unsubscribe, send the message "signoff gmt-help" to >>>> listserv@... >>> >>> To unsubscribe, send the message "signoff gmt-help" to >>> listserv@... >>> >> >> >> -- >> ******************************************************************************** >> >> Peter Schmidt Tel: +46-18-4717104 >> PhD Student Mobile: +46-73-3190975 >> Dept. of Earth Sciences e-mail: peter.schmidt@... >> Geophysics http://www.geofys.uu.se/?q=view/home/75 >> Uppsala University >> Villavagen 16 >> SE-75236 Uppsala >> ******************************************************************************** >> >> >> To unsubscribe, send the message "signoff gmt-help" to >> listserv@... > > To unsubscribe, send the message "signoff gmt-help" to > listserv@... > -- ******************************************************************************** Peter Schmidt Tel: +46-18-4717104 PhD Student Mobile: +46-73-3190975 Dept. of Earth Sciences e-mail: peter.schmidt@... Geophysics http://www.geofys.uu.se/?q=view/home/75 Uppsala University Villavagen 16 SE-75236 Uppsala ******************************************************************************** To unsubscribe, send the message "signoff gmt-help" to listserv@... |
| Free embeddable forum powered by Nabble | Forum Help |
