UTM to GMT. First usage

[jgarcia-2]

Hello all,
This is my first attempt to use GMT. I'm trying to pass some UTM data to
GMT and plot them, and I cannot find my way. It's a region of about 80km x
60km (lon,lat) in the Southeast of Spain, belonging to the UTM Zone and
Latitude Band 30S. My data combine non available data with height data
from a DEM, e.g.;

643225 4196375 NaN
643275 4196375 NaN
660625 4194925 522.3460083008
660675 4194925 521.3060302734
...

I've tried the commands:

makecpt -T0/1100/100 -Crelief > col.cpt
mapproject dem_gmt.txt -JU30S/10 -R-1.52/-0.61/37.5/38 -F  -I  > test
psxy test  -JU30S/10 -R-1.52/-0.61/37.5/38 -B -Ccol.cpt -Sc0.1  -K > test.ps
pscoast -R -JM -Dh -W2/100/100/250 -O >> test.ps

The coast line is perfect! Just what I needed, but the grid does not
appear. It seems my mapproject command is wrong. Could you help me in
these first steps?

In reality my data are not WGS84, but European Datum 1950; anyway perhaps
for my initial purposes of plotting the small spatial lag would not be
noticed. If the difference is noticeable, how would I explicit my xy data
are ED1950?


Thanks,

Javier
------------------

To unsubscribe, send the message "signoff gmt-help" to listserv@...

[Paul Wessel]

You need to use the -C option in mapproject so that your x/ys are  
considered relative to projection origin and not the lower left corner.

To deal with ED1950 etc, see mapproject -T (and use -Qd to see which  
datums are available).
-p

On Nov 6, 2008, at 12:54 AM, Javier Garcia-Pintado wrote:

> Hello all,
> This is my first attempt to use GMT. I'm trying to pass some UTM  
> data to
> GMT and plot them, and I cannot find my way. It's a region of about  
> 80km x
> 60km (lon,lat) in the Southeast of Spain, belonging to the UTM Zone  
> and
> Latitude Band 30S. My data combine non available data with height data
> from a DEM, e.g.;
>
> 643225 4196375 NaN
> 643275 4196375 NaN
> 660625 4194925 522.3460083008
> 660675 4194925 521.3060302734
> ...
>
> I've tried the commands:
>
> makecpt -T0/1100/100 -Crelief > col.cpt
> mapproject dem_gmt.txt -JU30S/10 -R-1.52/-0.61/37.5/38 -F  -I  > test
> psxy test  -JU30S/10 -R-1.52/-0.61/37.5/38 -B -Ccol.cpt -Sc0.1  -K >  
> test.ps
> pscoast -R -JM -Dh -W2/100/100/250 -O >> test.ps
>
> The coast line is perfect! Just what I needed, but the grid does not
> appear. It seems my mapproject command is wrong. Could you help me in
> these first steps?
>
> In reality my data are not WGS84, but European Datum 1950; anyway  
> perhaps
> for my initial purposes of plotting the small spatial lag would not be
> noticed. If the difference is noticeable, how would I explicit my xy  
> data
> are ED1950?
>
>
> Thanks,
>
> Javier
> ------------------
>
> To unsubscribe, send the message "signoff gmt-help" to listserv@...

To unsubscribe, send the message "signoff gmt-help" to listserv@...

[jgarcia-2]

Thank you, Paul!
Perfect!!!
Javier
-------------------

> You need to use the -C option in mapproject so that your x/ys are
> considered relative to projection origin and not the lower left corner.
>
> To deal with ED1950 etc, see mapproject -T (and use -Qd to see which
> datums are available).
> -p
>
> On Nov 6, 2008, at 12:54 AM, Javier Garcia-Pintado wrote:
>
>> Hello all,
>> This is my first attempt to use GMT. I'm trying to pass some UTM
>> data to
>> GMT and plot them, and I cannot find my way. It's a region of about
>> 80km x
>> 60km (lon,lat) in the Southeast of Spain, belonging to the UTM Zone
>> and
>> Latitude Band 30S. My data combine non available data with height data
>> from a DEM, e.g.;
>>
>> 643225 4196375 NaN
>> 643275 4196375 NaN
>> 660625 4194925 522.3460083008
>> 660675 4194925 521.3060302734
>> ...
>>
>> I've tried the commands:
>>
>> makecpt -T0/1100/100 -Crelief > col.cpt
>> mapproject dem_gmt.txt -JU30S/10 -R-1.52/-0.61/37.5/38 -F  -I  > test
>> psxy test  -JU30S/10 -R-1.52/-0.61/37.5/38 -B -Ccol.cpt -Sc0.1  -K >
>> test.ps
>> pscoast -R -JM -Dh -W2/100/100/250 -O >> test.ps
>>
>> The coast line is perfect! Just what I needed, but the grid does not
>> appear. It seems my mapproject command is wrong. Could you help me in
>> these first steps?
>>
>> In reality my data are not WGS84, but European Datum 1950; anyway
>> perhaps
>> for my initial purposes of plotting the small spatial lag would not be
>> noticed. If the difference is noticeable, how would I explicit my xy
>> data
>> are ED1950?
>>
>>
>> Thanks,
>>
>> Javier
>> ------------------
>>
>> 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@...