[GRASS-user] Importing GLOBE DEM data

gotnidea_Business_De · October 28, 2006, 5:44am

Hi,

I have been trying to import GLOBE DEM file from:

http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912

I don’t know what to do next. I have downloaded the file and unzipped it to get three files .bil, .hdr and one color file. How should I go about using it?

Can anyone please guide me? I am using WinGRASS6.3 CVS Native on Windows XP,

thanks in advance,

Abhinav

Markus_Neteler1 · October 28, 2006, 10:04am

Hi,

r.in.gdal should import it, see

http://grass.itc.it/grass63/manuals/html63_user/r.in.gdal.html

Markus

On Sat, Oct 28, 2006 at 11:14:34AM +0530, gotnidea! Business Development Centre wrote:

   Hi,

   I have been trying to import GLOBE DEM file from:

   [1]http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912

   I don't know what to do next. I have downloaded the file and unzipped it
   to get three files .bil, .hdr and one color file. How should I go about
   using it?

   Can anyone please guide me? I am using WinGRASS6.3 CVS Native on Windows
   XP,

   thanks in advance,

   Abhinav

References

   Visible links
   1. http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912

_______________________________________________
grassuser mailing list
grassuser@grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

--
Markus Neteler <neteler itc it> http://mpa.itc.it/markus/
ITC-irst - Centro per la Ricerca Scientifica e Tecnologica
MPBA - Predictive Models for Biol. & Environ. Data Analysis
Via Sommarive, 18 - 38050 Povo (Trento), Italy

gotnidea_Business_De · October 28, 2006, 10:27am

Yes, it should. But I fail to find any documentation as to how it is done. Even on the page you have mentioned, in example it is given like this:

"To import GLOBE DEM tiles (approx 1km resolution, better than GTOPO30 DEM data), the user has to download additionally the related HDR file(s). Finally, e.g. the ‘terrain’ color table can be assigned to the imported map with r.colors. "

How am I supposed to do it?

Abhinav

On 10/28/06, Markus Neteler <neteler@itc.it> wrote:

Hi,

r.in.gdal should import it, see

http://grass.itc.it/grass63/manuals/html63_user/r.in.gdal.html

Markus

On Sat, Oct 28, 2006 at 11:14:34AM +0530, gotnidea! Business Development Centre wrote:

Hi,

I have been trying to import GLOBE DEM file from:

[1]http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912

I don’t know what to do next. I have downloaded the file and unzipped it
to get three files .bil, .hdr and one color file. How should I go about
using it?

Can anyone please guide me? I am using WinGRASS6.3 CVS Native on Windows
XP,

thanks in advance,

Abhinav

References

Visible links

http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912

grassuser mailing list
grassuser@grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

–
Markus Neteler http://mpa.itc.it/markus/
ITC-irst - Centro per la Ricerca Scientifica e Tecnologica
MPBA - Predictive Models for Biol. & Environ. Data Analysis
Via Sommarive, 18 - 38050 Povo (Trento), Italy

grassuser mailing list
grassuser@grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

mlennert · October 28, 2006, 12:16pm

gotnidea! Business Development Centre wrote:

Yes, it should. But I fail to find any documentation as to how it is done. Even on the page you have mentioned, in example it is given like this:
"To import GLOBE DEM tiles <http://www.ngdc.noaa.gov/seg/topo/gltiles.shtml> (approx 1km resolution, better than GTOPO30 DEM data), the user has to download additionally the related HDR file(s) <ftp://ftp.ngdc.noaa.gov/GLOBE_DEM/data/elev/esri/hdr/>\. Finally, e.g. the 'terrain' color table can be assigned to the imported map with /r.colors/. "
How am I supposed to do it?

In general, please don't ask questions like "How am I supposed to it", but rather questions like: "After having read the manual I tried the following command, but I got the following error message" or (if this is your only problem) "How do I run r.in.gdal from the GUI ?"

If the man page is not clear, please tell us what information you are missing exactly.

In order to access r.in.gdal from the GUI, go to File->Import->Raster Map->"Multiple formats using GDAL".

Your data already contains the hdr and clr files.
But if I understand correctly, the data file actually does not contain the projection information, so if you use r.in.gdal with the location= option, it will just create an xy (unprojected) location.

So you have to create a location with the relevant projection first (see http://www.ngdc.noaa.gov/mgg/topo/report/s11/s11H.html).

However, I agree that this is currently a bit difficult in winGRASS...

A solution at this stage might be to use g.proj to create the new location (but I don't have access to wingrass right now to test).

In an existing location (for example spearfish) go to Config->Manage Projections->"Show projection information..." and enter the following parameters:

- check "Create new projection file" flag
- PROJ.4 projection description: +proj=longlat +datum=wgs84
- Name of new location to create: GUJARAT (or whatever you like)

Click on Run.

Now open GRASS in the new location and run r.in.gdal from there. You will have to check the "Override projection" flag.

Moritz

Abhinav

On 10/28/06, *Markus Neteler* <neteler@itc.it <mailto:neteler@itc.it>> wrote:

    Hi,

    r.in.gdal should import it, see

    http://grass.itc.it/grass63/manuals/html63_user/r.in.gdal.html

    Markus

    On Sat, Oct 28, 2006 at 11:14:34AM +0530, gotnidea! Business
    Development Centre wrote:
     > Hi,
     >
     > I have been trying to import GLOBE DEM file from:
     >
     > [1]http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912
     >
     > I don't know what to do next. I have downloaded the file and
    unzipped it
     > to get three files .bil, .hdr and one color file. How should I
    go about
     > using it?
     >
     > Can anyone please guide me? I am using WinGRASS6.3 CVS Native
    on Windows
     > XP,
     >
     > thanks in advance,
     >
     > Abhinav
     >
     > References
     >
     > Visible links
     > 1.
    http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912
    <http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912>

     > _______________________________________________
     > grassuser mailing list
     > grassuser@grass.itc.it <mailto:grassuser@grass.itc.it>
     > http://grass.itc.it/mailman/listinfo/grassuser

    --
    Markus Neteler <neteler itc it> http://mpa.itc.it/markus/
    ITC-irst - Centro per la Ricerca Scientifica e Tecnologica
    MPBA - Predictive Models for Biol. & Environ. Data Analysis
    Via Sommarive, 18 - 38050 Povo (Trento), Italy

    _______________________________________________
    grassuser mailing list
    grassuser@grass.itc.it <mailto:grassuser@grass.itc.it>
    http://grass.itc.it/mailman/listinfo/grassuser
    <http://grass.itc.it/mailman/listinfo/grassuser>

------------------------------------------------------------------------

_______________________________________________
winGRASS mailing list
winGRASS@grass.itc.it
http://grass.itc.it/mailman/listinfo/wingrass

mlennert · October 28, 2006, 2:20pm

gotnidea! Business Development Centre wrote:

I would be more precise with exactly what I am not able to get this time around. But I hope you consider that I am still a newbie with GRASS:

I understand, no harm in not knowing how to do things. But formulating the right question makes it much easier for us to help you.

Your data already contains the hdr and clr files.

To start with, I am not finding a way as to how these two files could be used to define the data I have downloaded. Where exactly they are to be placed or how GRASS would locate them. I am unable to find any documentation which could guide me in this regard. The documentation at NGDC is primarily for UNIX.

I think that r.in.gdal uses them automatically, but I'm not 100% sure.

    In an existing location (for example spearfish) go to Config->Manage
    Projections->"Show projection information..." and enter the following
    parameters:

    - check "Create new projection file" flag
    - PROJ.4 projection description: +proj=longlat +datum=wgs84
    - Name of new location to create: GUJARAT (or whatever you like)

    Click on Run.

You can actually open GRASS as many times as you want as long as you open it with a different location. So technically you are not obliged to close the spearfish session in order to be able to open the gujarat session. But unless you have good reasons to keep both open, I would recommend closing any GRASS session you are not actively using, just to reduce the amount of open windows and the risk of working in the wrong GRASS instance.

Now open GRASS in the new location and run r.in.gdal from there. You
will have to check the "Override projection" flag.

When I restarted GRASS with new location. I am able to generate the map, but it is all white (Blank). I believe I need to use some command for adding colours and may be I have to use the header file too.

I think it is probably due to the region settings, not to color issues. Try "zoom to the selected map" in the Map Display window (10th button from the left).

Moritz

mlennert · October 28, 2006, 3:37pm

(Please continue to send your mails to the users list.)

gotnidea! Business Development Centre wrote:

Yes by zooming it to selected map I could generate the outline of
what I find to be a map of Gujarat. But it is just in two colours:
Black and White.

I don't get the same. Did you have the .clr file in the same directory
as the .bil file when you imported that with r.in.gdal ?

How can I get colours?

As Markus suggested, you can use r.colors.

The documentation refers to r.colors. Which I am able to get from
Raster>>Manage Map Colors>>... Now which option should I undertake
and how can I link it to the color file I have downloaded.

I am sorry, but this again sounds like you did not read the manual (I am not saying you didn't, but the way you formulate your question just makes it sound that way). Did you try to run the command ? With which settings ? Please try to work out these things on your own, before you ask others.

The manual:
http://grass.itc.it/grass63/manuals/html63_user/r.colors.html (or just push 'help' button in the module window).

In brief, to do what Markus suggested:

set

map=YourMap
rules=terrain

To see the new colors applied, you have to push the "Redraw all layers" button (2nd from the left) in the Map Display.

Moritz

Glynn_Clements1 · October 28, 2006, 5:46pm

gotnidea! Business Development Centre wrote:

I have been trying to import GLOBE DEM file from:

http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912

I don't know what to do next. I have downloaded the file and unzipped it to
get three files .bil, .hdr and one color file. How should I go about using
it?

Can anyone please guide me? I am using WinGRASS6.3 CVS Native on Windows XP,

If you are new to GRASS, I would suggest using the Cygwin version
instead of the native version.

The native Windows version is labelled "experimental" for a reason. It
is essentially an "alpha-test" version, i.e. not even at a "beta-test"
level of stability.

--
Glynn Clements <glynn@gclements.plus.com>

Maciej_Sieczka2 · October 28, 2006, 9:09pm

Moritz Lennert wrote:

You can actually open GRASS as many times as you want as long as you
open it with a different location.

Just for the record - only the *mapset* has to be different. The
location might be the same.

Cheers,

P.S.

Moritz, congrats on your patience :).

H_B · October 29, 2006, 4:53am

I have been trying to import GLOBE DEM file from:
http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912

I don't know what to do next. I have downloaded the file and unzipped
it to get three files .bil, .hdr and one color file. How should I go
about using it?

Import with r.in.gdal:
http://gdal.org/frmt_various.html#EHdr

You should also be able to import it with r.in.bin, but you will
probably have to set some parameters by hand that way.

Hamish

Luigi_Ponti · October 29, 2006, 5:32am

Hi Abhinav,

Importing the GLOBE DEM was the very first thing I did with GRASS (cygwin). I used Huidae Cho's instructions:
http://geni.ath.cx/grass.html#toc21

And it wasn't that hard. Did you try that?
Regards,

Luigi

gotnidea! Business Development Centre wrote:

Hi,
I have been trying to import GLOBE DEM file from:
http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912
I don't know what to do next. I have downloaded the file and unzipped it to get three files .bil, .hdr and one color file. How should I go about using it?
Can anyone please guide me? I am using WinGRASS6.3 CVS Native on Windows XP,
thanks in advance,
Abhinav
------------------------------------------------------------------------

_______________________________________________
winGRASS mailing list
winGRASS@grass.itc.it
http://grass.itc.it/mailman/listinfo/wingrass

gotnidea_Business_De · October 29, 2006, 6:19am

Thanks to all for the guidance.

I have been able to import and view the data using r.in.gdal. The only problem is that it is in Black and White. I tried to use this command:

r.colors map=Gujarat color=rules rast=Gujarat rules=terrain

I hope this is what Moritz meant:

In brief, to do what Markus suggested:

set

map=YourMap
rules=terrain

To see the new colors applied, you have to push the “Redraw all layers”
button (2nd from the left) in the Map Display.

The map is now in shades of grey.

The other Data I downloaded is : freegis_worlddata-0.1_simpl.tar

It is also in Black and White with only the outlines and crosses in black. I have made a separate location for this data. I am unable to have a coloured map in it also. It’s GRASS how to file (which comes along with the data mentions:

Create location

Location:

Mapset:

Database: <a path to a writable directory

Coordinate system: B, Latitude-Longitude

description:

geodatic datum: no

ellipsoid: grs80

region:

north: 90

south: -90

west: -180

east: 180

grid resolution: not important

When I tried to do this, I am getting this error : ‘clear’ is not recognized as an internal or external command, operable program or batch file.

However I did try to import the data as is given in the second point:

import data

v.in.ogr dsn=/usr/share/geo-data/country_geodata/countries_simpl.shp output=countries

This was comfortably done.

you can now load /usr/share/geo-data/country_geodata/pop.gm with the GRASS

data manager to have a first look on the data and to print them.

How can I undertake this step. I am a bit clueless.

Print colorcoded maps

#start monitor fo visualization

d.mon x0

print china in red

d.vect map=countries type=area where=“fibs=‘CH’” fcolor=255:0:0

print countries with total population < 5000000 in green

d.vect map=countries type=area where=“total<5000000” fcolor=0:255:0

erase all objects from monitor

d.erase

stop monitor

d.mon stop=x0

This last step mentioned in the How to… is again a bit confusing. I hope this is what is supposed to be written in the Output - GIS.m and run. Is it?

Thanks everyone for your patience and efforts towards helping me out. As for going back to Cygwin, I would only take it up as the last option available.

Thanks again for the help

regards

Abhinav

On 10/29/06, Hamish <hamish_nospam@yahoo.com> wrote:

I have been trying to import GLOBE DEM file from:
http://www.ngdc.noaa.gov/cgi-bin/seg/ff/nph-newform.pl/seg/topo/customdatacd?disp=one&sid=peak.140.172.184.251.0610272224463912

I don’t know what to do next. I have downloaded the file and unzipped
it to get three files .bil, .hdr and one color file. How should I go
about using it?

Import with r.in.gdal:
http://gdal.org/frmt_various.html#EHdr

You should also be able to import it with r.in.bin, but you will
probably have to set some parameters by hand that way.

Hamish

mlennert · October 29, 2006, 12:21pm

gotnidea! Business Development Centre wrote:

Thanks to all for the guidance.
I have been able to import and view the data using r.in.gdal. The only problem is that it is in Black and White. I tried to use this command:
r.colors map=Gujarat color=rules rast=Gujarat rules=terrain
I hope this is what Moritz meant:

no:

r.colors map=Gujarat rules=terrain

the color= and rast= options are two different ways of setting color tables, using all of them at the same time does not make sense (maybe the module should actually complain about this ?).

The other Data I downloaded is : freegis_worlddata-0.1_simpl.tar

In case you did not know: this is vector data.

It is also in Black and White with only the outlines and crosses in black. I have made a separate location for this data. I am unable to have a coloured map in it also. It's GRASS how to file (which comes along with the data mentions:

    1. Create location

    Location: <some name>

    Mapset: <a name>

    Database: <a path to a writable directory

    Coordinate system: B, Latitude-Longitude

    description: <some description>

    geodatic datum: no

    ellipsoid: grs80

    region:

    north: 90

    south: -90

    west: -180

    east: 180

    grid resolution: not important

When I tried to do this, I am getting this error : 'clear' is not recognized as an internal or external command, operable program or batch file.

Yes, as mentioned earlier on the list, creation of locations at startup currently does not work.

    3. you can now load /usr/share/geo-data/country_geodata/pop.gm with
    the GRASS

    data manager to have a first look on the data and to print them.

How can I undertake this step. I am a bit clueless.

I don't think you can as it is probably for the old gui (d.m).

    4. Print colorcoded maps

    #start monitor fo visualization

    d.mon x0

    # print china in red

    d.vect map=countries type=area where="fibs='CH'" fcolor=255:0:0

    # print countries with total population < 5000000 in green

    d.vect map=countries type=area where="total<5000000" fcolor=0:255:0

    # erase all objects from monitor

    d.erase

    # stop monitor

    d.mon stop=x0

This last step mentioned in the How to... is again a bit confusing. I hope this is what is supposed to be written in the Output - GIS.m and run. Is it?

No, these are commands for the command line. But you cannot use them in wingrass as they depend on the old GRASS monitor display system which is not available in the wingrass gui (gis.m - GIS Manager).

Currently you cannot use any display commands from the command line if you use the GIS Manager GUI (this is true for *nix as well).

All of these commands are now integrated into the gui.

To reproduce this example:

- click on the "Add vector layer" button
- click on the vector1 label

Beneath the layer tree, you will see appear the vector panel.

Fill out the relevant fields (e.g. "Vector name", "SQL query", etc)

If I have the time one of these days, I might make some detailed movies, but not now.

Moritz

gotnidea_Business_De · October 30, 2006, 4:21am

Thanks.

By using the vector panel I could add colors to the map freegis_worlddata. But I couldn’t do the same with my Gujarat map I downloaded from the NGDC.

Now, I believe, I should ask one question which really forced me to work on GRASS in the first place. I want to have a display like:

http://www.eicinformation.org/data/sampledata.asp?id=2

Is it possible on GRASS?? If at all it is not possible on a web page, can I reproduce a similar map on my desktop, linking the data with SQL or even by using the existing postgreSQL?

I wish to achieve a level of competency in GRASS (or any other software which is freely available)which could help me produce such a work.

regards
Abhinav

On 10/29/06, Moritz Lennert <mlennert@club.worldonline.be> wrote:

gotnidea! Business Development Centre wrote:

Thanks to all for the guidance.

I have been able to import and view the data using r.in.gdal. The only
problem is that it is in Black and White. I tried to use this command:

r.colors map=Gujarat color=rules rast=Gujarat rules=terrain

I hope this is what Moritz meant:

no:

r.colors map=Gujarat rules=terrain

the color= and rast= options are two different ways of setting color
tables, using all of them at the same time does not make sense (maybe
the module should actually complain about this ?).

The other Data I downloaded is : freegis_worlddata-0.1_simpl.tar

In case you did not know: this is vector data.

It is also in Black and White with only the outlines and crosses in
black. I have made a separate location for this data. I am unable to
have a coloured map in it also. It’s GRASS how to file (which comes
along with the data mentions:

Create location

Location:

Mapset:

Database: <a path to a writable directory

Coordinate system: B, Latitude-Longitude

description:

geodatic datum: no

ellipsoid: grs80

region:

north: 90

south: -90

west: -180

east: 180

grid resolution: not important

When I tried to do this, I am getting this error : ‘clear’ is not
recognized as an internal or external command, operable program or batch
file.

Yes, as mentioned earlier on the list, creation of locations at startup
currently does not work.

you can now load /usr/share/geo-data/country_geodata/pop.gm with
the GRASS

data manager to have a first look on the data and to print them.

How can I undertake this step. I am a bit clueless.

I don’t think you can as it is probably for the old gui ( d.m).

Print colorcoded maps

#start monitor fo visualization

d.mon x0

print china in red

d.vect map=countries type=area where=“fibs=‘CH’” fcolor=255:0:0

print countries with total population < 5000000 in green

d.vect map=countries type=area where=“total<5000000” fcolor=0:255:0

erase all objects from monitor

d.erase

stop monitor

d.mon stop=x0

This last step mentioned in the How to… is again a bit confusing. I
hope this is what is supposed to be written in the Output - GIS.m and
run. Is it?

No, these are commands for the command line. But you cannot use them in
wingrass as they depend on the old GRASS monitor display system which is
not available in the wingrass gui ( gis.m - GIS Manager).

Currently you cannot use any display commands from the command line if
you use the GIS Manager GUI (this is true for *nix as well).

All of these commands are now integrated into the gui.

To reproduce this example:

click on the “Add vector layer” button

click on the vector1 label

Beneath the layer tree, you will see appear the vector panel.

Fill out the relevant fields ( e.g. “Vector name”, “SQL query”, etc)

If I have the time one of these days, I might make some detailed movies,
but not now.

Moritz

Luca_Penasa · October 30, 2006, 10:50am

im trying to export a map using the ps.map command...
im using grass 6.3 from cvs and the output is;

WARNING: Data exists after final 'end' instruction!
PS-PAINT: scale set to 1 : -168897.
PS-PAINT: reading raster file <mappa in licheni> ...
PS-PAINT: reading labels file <ibl in licheni> ...
PS-PAINT: reading vector points file <sites_203 in licheni> ...
PS-PAINT: reading text file ...
ERROR: PSmap: do_scale(): shouldn't happen

Any idea?

luca

mlennert · October 30, 2006, 10:09am

Luca Penasa wrote:

im trying to export a map using the ps.map command...
im using grass 6.3 from cvs and the output is;

WARNING: Data exists after final 'end' instruction!

here's the problem: you have an end somewhere where you probably shouldn't and so it does not take into account any subsequent commands.

If you need more help you should post your ps.map instructions file.

Moritz

Luca_Penasa · October 30, 2006, 12:14pm

Moritz Lennert wrote:

Luca Penasa wrote:

im trying to export a map using the ps.map command...
im using grass 6.3 from cvs and the output is;

WARNING: Data exists after final 'end' instruction!

here's the problem: you have an end somewhere where you probably shouldn't and so it does not take into account any subsequent commands.

If you need more help you should post your ps.map instructions file.

Moritz

i solve this deleting an 'end' at eof. but the problem is this:

PS-PAINT: scale set to 1 : -168897.
PS-PAINT: reading raster file <mappa in licheni> ...
PS-PAINT: reading labels file <ibl in licheni> ...
PS-PAINT: reading vector points file <sites_203 in licheni> ...
PS-PAINT: reading text file ...
*ERROR: PSmap: do_scale(): shouldn't happen *

H_B · October 30, 2006, 11:23am

Luca Penasa wrote:

PS-PAINT: scale set to 1 : -168897.
PS-PAINT: reading raster file <mappa in licheni> ...
PS-PAINT: reading labels file <ibl in licheni> ...
PS-PAINT: reading vector points file <sites_203 in licheni> ...
PS-PAINT: reading text file ...
*ERROR: PSmap: do_scale(): shouldn't happen *

... is this a lat/lon location? ps.map scale is broken (meaningless?)
for those. (see bug #3096)

otherwise we would need to see your ps.map instructions file,
'g.region -p' settings, and 'g.proj -p' settings to help.

Hamish

mlennert · October 30, 2006, 12:24pm

Luca Penasa wrote:

Moritz Lennert wrote:

Luca Penasa wrote:

im trying to export a map using the ps.map command...
im using grass 6.3 from cvs and the output is;

WARNING: Data exists after final 'end' instruction!

here's the problem: you have an end somewhere where you probably shouldn't and so it does not take into account any subsequent commands.

If you need more help you should post your ps.map instructions file.

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

Moritz

Moritz

i solve this deleting an 'end' at eof. but the problem is this:

PS-PAINT: scale set to 1 : -168897.
PS-PAINT: reading raster file <mappa in licheni> ...
PS-PAINT: reading labels file <ibl in licheni> ...
PS-PAINT: reading vector points file <sites_203 in licheni> ...
PS-PAINT: reading text file ...
*ERROR: PSmap: do_scale(): shouldn't happen *

_______________________________________________
grassuser mailing list
grassuser@grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

Luca_Penasa · October 30, 2006, 3:36pm

Hamish wrote:

Luca Penasa wrote:


PS-PAINT: scale set to 1 : -168897.
PS-PAINT: reading raster file <mappa in licheni> ...
PS-PAINT: reading labels file <ibl in licheni> ...
PS-PAINT: reading vector points file <sites_203 in licheni> ...
PS-PAINT: reading text file ...
*ERROR: PSmap: do_scale(): shouldn't happen *

... is this a lat/lon location? ps.map scale is broken (meaningless?)
for those. (see bug #3096)

Its a gauss-boaga location...
now iv installed a old version of grass (6.0.2)... with this version the same instructions file make ps.map works... cold not to be a bug?

otherwise we would need to see your ps.map instructions file,
'g.region -p' settings, and 'g.proj -p' settings to help.

Hamish

Luca_Penasa · October 30, 2006, 7:08pm

Markus Neteler wrote:

Luca,
please post
  g.region -p
and
  g.proj -w
to me.

thanks
markus

Luca Penasa wrote on 10/30/2006 04:36 PM:


Hamish wrote:


Luca Penasa wrote:

PS-PAINT: scale set to 1 : -168897.
PS-PAINT: reading raster file <mappa in licheni> ...
PS-PAINT: reading labels file <ibl in licheni> ...
PS-PAINT: reading vector points file <sites_203 in licheni> ...
PS-PAINT: reading text file ...
*ERROR: PSmap: do_scale(): shouldn't happen *


... is this a lat/lon location? ps.map scale is broken (meaningless?)
for those. (see bug #3096)


Its a gauss-boaga location...
now iv installed a old version of grass (6.0.2)... with this version
the same instructions file make ps.map works... cold not to be a bug?


otherwise we would need to see your ps.map instructions file,
'g.region -p' settings, and 'g.proj -p' settings to help.

Hamish

_______________________________________________
grassuser mailing list
grassuser@grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

GRASS 6.0.2 (lichTest_OK):~ > g.region -p
projection: 99 (Transverse Mercator)
zone: 0
datum: rome40
ellipsoid: international
north: 5108040
south: 5099460
west: 1662320
east: 1666520
nsres: 0.5
ewres: 0.5
rows: 17160
cols: 8400

GRASS 6.0.2 (lichTest_OK):~ > g.proj -w
PROJCS["Transverse Mercator",
    GEOGCS["international",
        DATUM["Monte_Mario",
            SPHEROID["international",6378388,297],
            TOWGS84[-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",9.00000000333333],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",1500000],
    PARAMETER["false_northing",0],
    UNIT["meter",1]]