[GRASSLIST:8825] GRASS 6 in a nutshell tutorial

Luigi_Ponti · October 31, 2005, 7:32pm

Hi there,

while going through the Neteler's GRASS 6 in a nutshell tutorial I get the message:

"Cannot open database /home/Administrator/grassdata/spearfish60/user1/dbf/ by driver dbf"

when at section 3.1.2. I try to import vector ESRI SHAPE files form the package tiger2000_UTM13_nad27.tar.gz, via the "import various formats using ogr" option in the GIS manager.

Does anybody know what the problem could be?

Thanks,

Luigi

============================================================================
System Information

Platform Unix
----------------------------------------------------------------------------
Operating System CYGWIN_NT-5.1
----------------------------------------------------------------------------
OS Version 1.5.18(0.132/4/2)
----------------------------------------------------------------------------
Processor Type (arch) i686
----------------------------------------------------------------------------
Full OS ID CYGWIN_NT-5.1-1.5.18(0.132/4/2)-i686
----------------------------------------------------------------------------
Tcl Version/Patchlevel 8.4.7
----------------------------------------------------------------------------
Name of Executable /usr/local/bin/wish
----------------------------------------------------------------------------
Name of current Script /usr/local/grass6.0.0/etc/dm/tksys.tcl
----------------------------------------------------------------------------

Markus_Neteler1 · November 2, 2005, 10:55am

Luigi Ponti wrote:

Hi there,

while going through the Neteler's GRASS 6 in a nutshell tutorial I get
the message:

"Cannot open database
/home/Administrator/grassdata/spearfish60/user1/dbf/ by driver dbf"

when at section 3.1.2. I try to import vector ESRI SHAPE files form
the package tiger2000_UTM13_nad27.tar.gz, via the "import various
formats using ogr" option in the GIS manager.

Does anybody know what the problem could be?

Thanks,

Luigi,

please check if the directory
/home/Administrator/grassdata/spearfish60/user1/dbf/
exists. If not, please create it. From current GRASS 6.1-CVS onwards the
dbf/ directory in the current
mapset is created automatically.

Hope this helps

Markus

Luigi_Ponti · November 3, 2005, 11:58pm

Markus Neteler wrote:

Luigi Ponti wrote:

Hi there,

while going through the Neteler's GRASS 6 in a nutshell tutorial I get
the message:

"Cannot open database
/home/Administrator/grassdata/spearfish60/user1/dbf/ by driver dbf"

when at section 3.1.2. I try to import vector ESRI SHAPE files form
the package tiger2000_UTM13_nad27.tar.gz, via the "import various
formats using ogr" option in the GIS manager.

Does anybody know what the problem could be?

Thanks,

Luigi,

please check if the directory
/home/Administrator/grassdata/spearfish60/user1/dbf/
exists. If not, please create it. From current GRASS 6.1-CVS onwards the
dbf/ directory in the current
mapset is created automatically.

Hope this helps

Markus

Thanks Markus,

It worked out perfectly fine.

Now I would need to create my own location with a state boundary and an elevation model for California (this would provide the basis to map the output of my ecological model --reminder). Does anybody know if there is online files available for download that would allow me to accomplish such task? I have been trying to do this by OGR importing ESRI .shp files from http://gis.ca.gov/ into a location with projection based on .prj file content, but the new layer does not show up in the monitor. Any hints?

Also when I try to access GRASS help, the following message shows up:
"help /usr/local/grass6.0.0/scripts/g.manual: line 71: iexplore: command not found".
Is it that GRASS cannot find the browser?

Finally, when trying to access the man version via "g.manual -m d.rast", the shell says:
"fopen: No such file or directory
Cannot open man page /home/Administrator/C:\cygwin/usr/local/grass6.0.0/man/man1/d.rast.1.gz
No manual entry for C:\cygwin/usr/local/grass6.0.0/man/man1/d.rast.1".
Is there any fix to this?

Ciao,

Luigi

Ian_MacMillan1 · November 4, 2005, 8:13am

Luigi, I am not sure about your manual problems, but you can grab data of California from several location for free. Try http://data.geocomm.com/ for a vector outline of the state. You have to create an account on that site, but it is free, and they have a lot of useful data.

The best DEM of the state can be had from the USGS seamless server here: http://seamless.usgs.gov. DEMs of a variety of resolutions are available there.

If you have already imported data and can't view it, have you made sure to set your region?
g.region rast=your_rast
or
g.region vect=your_vect
then
d.rast rast=your_rast

Hope this helps,
Ian

On Nov 3, 2005, at 3:58 PM, Luigi Ponti wrote:

Markus Neteler wrote:

Luigi Ponti wrote:

Hi there,

while going through the Neteler's GRASS 6 in a nutshell tutorial I get
the message:

"Cannot open database
/home/Administrator/grassdata/spearfish60/user1/dbf/ by driver dbf"

when at section 3.1.2. I try to import vector ESRI SHAPE files form
the package tiger2000_UTM13_nad27.tar.gz, via the "import various
formats using ogr" option in the GIS manager.

Does anybody know what the problem could be?

Thanks,

Luigi,

please check if the directory
/home/Administrator/grassdata/spearfish60/user1/dbf/
exists. If not, please create it. From current GRASS 6.1-CVS onwards the
dbf/ directory in the current
mapset is created automatically.

Hope this helps

Markus

Thanks Markus,

It worked out perfectly fine.

Now I would need to create my own location with a state boundary and an elevation model for California (this would provide the basis to map the output of my ecological model --reminder). Does anybody know if there is online files available for download that would allow me to accomplish such task? I have been trying to do this by OGR importing ESRI .shp files from http://gis.ca.gov/ into a location with projection based on .prj file content, but the new layer does not show up in the monitor. Any hints?

Also when I try to access GRASS help, the following message shows up:
"help /usr/local/grass6.0.0/scripts/g.manual: line 71: iexplore: command not found".
Is it that GRASS cannot find the browser?

Finally, when trying to access the man version via "g.manual -m d.rast", the shell says:
"fopen: No such file or directory
Cannot open man page /home/Administrator/C:\cygwin/usr/local/grass6.0.0/man/man1/d.rast.1.gz
No manual entry for C:\cygwin/usr/local/grass6.0.0/man/man1/d.rast.1".
Is there any fix to this?

Ciao,

Luigi

>
What happens if a big asteroid hits Earth? Judging from realistic simulations involving a sledge hammer and a common laboratory frog, we can assume it will be pretty bad.
- Dave Barry

-------------------------------------------------------------
This message has been scanned by Postini anti-virus software.

Luigi_Ponti · November 9, 2005, 7:10pm

Thanks Ian,

I got the vector .shp outline of the state from the California Spatial Information Library http://gis.ca.gov/index.epl since the one you recemmended me on http://data.geocomm.com/ was in .e00 format and GRASS would not import it directly.

The DEM I got it from NOAA/GLOBE http://www.ngdc.noaa.gov/mgg/topo/globe.html where you can customize your download area as opposed to the DEM on http://seamless.usgs.gov that has to be downloaded in small pieces.

I have imported both the vector and the raster with minimum trouble, since the first has a .prj file which GRASS recognized, and the second had a detailed import tutorial in the cygwin-GRASS install page. However, despite the fact the the two layers supposedly have the same NAD83 projection info (I did not enter anything manually), when I try to display them together on a GRASS monitor, the state layout is in the middle of the ocean compared to the elevation raster.

Is there a way to overlap the two layers or to chek for a coordinate error?

Thanks,

Luigi

Ian MacMillan wrote:

Luigi, I am not sure about your manual problems, but you can grab data of California from several location for free. Try http://data.geocomm.com/ for a vector outline of the state. You have to create an account on that site, but it is free, and they have a lot of useful data.

The best DEM of the state can be had from the USGS seamless server here: http://seamless.usgs.gov. DEMs of a variety of resolutions are available there.

If you have already imported data and can't view it, have you made sure to set your region?
g.region rast=your_rast
or
g.region vect=your_vect
then
d.rast rast=your_rast

Hope this helps,
Ian

On Nov 3, 2005, at 3:58 PM, Luigi Ponti wrote:

Markus Neteler wrote:

Luigi Ponti wrote:

Hi there,

while going through the Neteler's GRASS 6 in a nutshell tutorial I get
the message:

"Cannot open database
/home/Administrator/grassdata/spearfish60/user1/dbf/ by driver dbf"

when at section 3.1.2. I try to import vector ESRI SHAPE files form
the package tiger2000_UTM13_nad27.tar.gz, via the "import various
formats using ogr" option in the GIS manager.

Does anybody know what the problem could be?

Thanks,

Luigi,

please check if the directory
/home/Administrator/grassdata/spearfish60/user1/dbf/
exists. If not, please create it. From current GRASS 6.1-CVS onwards the
dbf/ directory in the current
mapset is created automatically.

Hope this helps

Markus

Thanks Markus,

It worked out perfectly fine.

Now I would need to create my own location with a state boundary and an elevation model for California (this would provide the basis to map the output of my ecological model --reminder). Does anybody know if there is online files available for download that would allow me to accomplish such task? I have been trying to do this by OGR importing ESRI .shp files from http://gis.ca.gov/ into a location with projection based on .prj file content, but the new layer does not show up in the monitor. Any hints?

Also when I try to access GRASS help, the following message shows up:
"help /usr/local/grass6.0.0/scripts/g.manual: line 71: iexplore: command not found".
Is it that GRASS cannot find the browser?

Finally, when trying to access the man version via "g.manual -m d.rast", the shell says:
"fopen: No such file or directory
Cannot open man page /home/Administrator/C:\cygwin/usr/local/grass6.0.0/man/man1/ d.rast.1.gz
No manual entry for C:\cygwin/usr/local/grass6.0.0/man/man1/d.rast.1".
Is there any fix to this?

Ciao,

Luigi

>
What happens if a big asteroid hits Earth? Judging from realistic simulations involving a sledge hammer and a common laboratory frog, we can assume it will be pretty bad.
- Dave Barry

-------------------------------------------------------------
This message has been scanned by Postini anti-virus software.

Ian_MacMillan1 · November 10, 2005, 12:52am

Luigi, glad that you mostly got things to work. Good to know that there are alternate sources of data. I am not sure, but I would guess that you are using a UTM projection. I bet that the vector is using one zone, and the DEM is using another. Make sure that you import them into their true projection, then reproject the data as you see fit (I would reproject the vector to save time).

Just for completeness, I think GRASS can import e00 vectors without too much trouble with v.in.e00.

-Ian

On Nov 9, 2005, at 11:10 AM, Luigi Ponti wrote:

Thanks Ian,

I got the vector .shp outline of the state from the California Spatial Information Library http://gis.ca.gov/index.epl since the one you recemmended me on http://data.geocomm.com/ was in .e00 format and GRASS would not import it directly.

The DEM I got it from NOAA/GLOBE http://www.ngdc.noaa.gov/mgg/topo/globe.html where you can customize your download area as opposed to the DEM on http://seamless.usgs.gov that has to be downloaded in small pieces.

I have imported both the vector and the raster with minimum trouble, since the first has a .prj file which GRASS recognized, and the second had a detailed import tutorial in the cygwin-GRASS install page. However, despite the fact the the two layers supposedly have the same NAD83 projection info (I did not enter anything manually), when I try to display them together on a GRASS monitor, the state layout is in the middle of the ocean compared to the elevation raster.

Is there a way to overlap the two layers or to chek for a coordinate error?

Thanks,

Luigi

Ian MacMillan wrote:

Luigi, I am not sure about your manual problems, but you can grab data of California from several location for free. Try http://data.geocomm.com/ for a vector outline of the state. You have to create an account on that site, but it is free, and they have a lot of useful data.

The best DEM of the state can be had from the USGS seamless server here: http://seamless.usgs.gov. DEMs of a variety of resolutions are available there.

If you have already imported data and can't view it, have you made sure to set your region?
g.region rast=your_rast
or
g.region vect=your_vect
then
d.rast rast=your_rast

Hope this helps,
Ian

On Nov 3, 2005, at 3:58 PM, Luigi Ponti wrote:

Markus Neteler wrote:

Luigi Ponti wrote:

Hi there,

while going through the Neteler's GRASS 6 in a nutshell tutorial I get
the message:

"Cannot open database
/home/Administrator/grassdata/spearfish60/user1/dbf/ by driver dbf"

when at section 3.1.2. I try to import vector ESRI SHAPE files form
the package tiger2000_UTM13_nad27.tar.gz, via the "import various
formats using ogr" option in the GIS manager.

Does anybody know what the problem could be?

Thanks,

Luigi,

please check if the directory
/home/Administrator/grassdata/spearfish60/user1/dbf/
exists. If not, please create it. From current GRASS 6.1-CVS onwards the
dbf/ directory in the current
mapset is created automatically.

Hope this helps

Markus

Thanks Markus,

It worked out perfectly fine.

Now I would need to create my own location with a state boundary and an elevation model for California (this would provide the basis to map the output of my ecological model --reminder). Does anybody know if there is online files available for download that would allow me to accomplish such task? I have been trying to do this by OGR importing ESRI .shp files from http://gis.ca.gov/ into a location with projection based on .prj file content, but the new layer does not show up in the monitor. Any hints?

Also when I try to access GRASS help, the following message shows up:
"help /usr/local/grass6.0.0/scripts/g.manual: line 71: iexplore: command not found".
Is it that GRASS cannot find the browser?

Finally, when trying to access the man version via "g.manual -m d.rast", the shell says:
"fopen: No such file or directory
Cannot open man page /home/Administrator/C:\cygwin/usr/local/grass6.0.0/man/man1/ d.rast.1.gz
No manual entry for C:\cygwin/usr/local/grass6.0.0/man/man1/d.rast.1".
Is there any fix to this?

Ciao,

Luigi

>
What happens if a big asteroid hits Earth? Judging from realistic simulations involving a sledge hammer and a common laboratory frog, we can assume it will be pretty bad.
- Dave Barry

-------------------------------------------------------------
This message has been scanned by Postini anti-virus software.

>
What happens if a big asteroid hits Earth? Judging from realistic simulations involving a sledge hammer and a common laboratory frog, we can assume it will be pretty bad.
- Dave Barry

-------------------------------------------------------------
This message has been scanned by Postini anti-virus software.

H_B · November 10, 2005, 5:12am

I have imported both the vector and the raster with minimum trouble,
since the first has a .prj file which GRASS recognized, and the second

had a detailed import tutorial in the cygwin-GRASS install page.
However, despite the fact the the two layers supposedly have the same
NAD83 projection info (I did not enter anything manually), when I try
to display them together on a GRASS monitor, the state layout is in
the middle of the ocean compared to the elevation raster.

Is there a way to overlap the two layers or to chek for a coordinate
error?

r.info and v.info will give you the map extents, you should then be able
to figure out which is correct.

Also try "d.where -w" to give you a lat/lon value for different parts of
your map which you can compare to the map on the wall.

Hamish

Luigi_Ponti · November 14, 2005, 10:06pm

Thanks Amish,

looking at v.info and r.info, it seems the two maps have the same projection (see below). However, when looking at the output "d.where -w" for a sample map point (the San Francisco Bay), the longitude is very different in the two maps:

                                                   WGS84 Co-ordinates
             EAST: NORTH: LON:
              -191114.6346037 -24240.66895625 143:03:58.588289W 29:24:27.1
(San Francisco Bay on the California State vector layer - WRONG)

1839139.30311877 663981.00484798 122:19:50.783112W 37:57:48.1
(San Francisco Bay on the Elevation Model raster layer - CORRECT)

Do you think this could be due to an error in the import process or should I just look for another vector map?

Thanks again,

Luigi

v.info map=state layer=1
+----------------------------------------------------------------------------+
| Layer: state Organization: |
| Mapset: luigi Source Date: |
| Location: Cal_nad83-zone3 Name of creator: |
| Database: /home/andy |
| Title: |
| Map Scale: 1:1 |
| Map format: native |
|----------------------------------------------------------------------------|
| Type of Map: Vector (level: 2) |
| |
| Number of points: 0 Number of areas: 51 |
| Number of lines: 0 Number of islands: 48 |
| Number of boundaries: 59 Number of faces: 0 |
| Number of centroids: 51 Number of kernels: 0 |
| |
| Map is 3D: 0 |
| Number of dblinks: 1 |
| |
| Projection: Lambert Conformal Conic (zone 0) |
| N: 922206.900 S: -701781.905 |
| E: 1043412.917 W: -374425.808 |
| B: 0.000 T: 0.000 |
| |
| Digitize threshold: 0.00000 |
| Comments: |
| |
+----------------------------------------------------------------------------+

r.info map=e10g@PERMANENT
+----------------------------------------------------------------------------+
| Layer: e10g Date: Tue Nov 8 17:25:24 2005 |
| Mapset: PERMANENT Login of Creator: andy |
| Location: Cal_nad83-zone3 |
| DataBase: /home/andy |
| Title: ( e10g ) |
|----------------------------------------------------------------------------|
| |
| Type of Map: raster Number of Categories: 5443 |
| Data Type: CELL |
| Rows: 7192 |
| Columns: 11628 |
| Total Cells: 83628576 |
| Projection: Lambert Conformal Conic (zone 0) |
| N: 3353466 S: -3838238 Res: 999.95884316 |
| E: 6073284 W: -5555183 Res: 1000.04016168 |
| Range of data: min = -84 max = 5443 |
| |
| Data Source: |
| |
| |
| |
| Data Description: |
| generated by r.proj |
| |
| |
+----------------------------------------------------------------------------+

Hamish wrote:

I have imported both the vector and the raster with minimum trouble, since the first has a .prj file which GRASS recognized, and the second

had a detailed import tutorial in the cygwin-GRASS install page. However, despite the fact the the two layers supposedly have the same NAD83 projection info (I did not enter anything manually), when I try
to display them together on a GRASS monitor, the state layout is in
the middle of the ocean compared to the elevation raster.

Is there a way to overlap the two layers or to chek for a coordinate
error?

r.info and v.info will give you the map extents, you should then be able
to figure out which is correct.

Also try "d.where -w" to give you a lat/lon value for different parts of your map which you can compare to the map on the wall.

Hamish

Ian_MacMillan1 · November 15, 2005, 2:31am

Luigi, if you look at r.info and v.info, you will only find the projection information of your location, not necessarily the projection/datum/etc. information of the original data (I believe). Do you have the original metadata? If not, try running gdalinfo on your original data files to determine their original projection info.

man gdalinfo
will give you the correct syntax.

G'luck
Ian

On Nov 14, 2005, at 2:06 PM, Luigi Ponti wrote:

Thanks Amish,

looking at v.info and r.info, it seems the two maps have the same projection (see below). However, when looking at the output "d.where -w" for a sample map point (the San Francisco Bay), the longitude is very different in the two maps:

                                                  WGS84 Co-ordinates
            EAST: NORTH: LON:
            -191114.6346037 -24240.66895625 143:03:58.588289W 29:24:27.1
(San Francisco Bay on the California State vector layer - WRONG)

1839139.30311877 663981.00484798 122:19:50.783112W 37:57:48.1
(San Francisco Bay on the Elevation Model raster layer - CORRECT)

Do you think this could be due to an error in the import process or should I just look for another vector map?

Thanks again,

Luigi

v.info map=state layer=1
+----------------------------------------------------------------------------+
| Layer: state Organization: |
| Mapset: luigi Source Date: |
| Location: Cal_nad83-zone3 Name of creator: |
| Database: /home/andy |
| Title: |
| Map Scale: 1:1 |
| Map format: native |
|----------------------------------------------------------------------------|
| Type of Map: Vector (level: 2) |
| |
| Number of points: 0 Number of areas: 51 |
| Number of lines: 0 Number of islands: 48 |
| Number of boundaries: 59 Number of faces: 0 |
| Number of centroids: 51 Number of kernels: 0 |
| |
| Map is 3D: 0 |
| Number of dblinks: 1 |
| |
| Projection: Lambert Conformal Conic (zone 0) |
| N: 922206.900 S: -701781.905 |
| E: 1043412.917 W: -374425.808 |
| B: 0.000 T: 0.000 |
| |
| Digitize threshold: 0.00000 |
| Comments: |
| |
+----------------------------------------------------------------------------+

r.info map=e10g@PERMANENT
+----------------------------------------------------------------------------+
| Layer: e10g Date: Tue Nov 8 17:25:24 2005 |
| Mapset: PERMANENT Login of Creator: andy |
| Location: Cal_nad83-zone3 |
| DataBase: /home/andy |
| Title: ( e10g ) |
|----------------------------------------------------------------------------|
| |
| Type of Map: raster Number of Categories: 5443 |
| Data Type: CELL |
| Rows: 7192 |
| Columns: 11628 |
| Total Cells: 83628576 |
| Projection: Lambert Conformal Conic (zone 0) |
| N: 3353466 S: -3838238 Res: 999.95884316 |
| E: 6073284 W: -5555183 Res: 1000.04016168 |
| Range of data: min = -84 max = 5443 |
| |
| Data Source: |
| |
| |
| |
| Data Description: |
| generated by r.proj |
| |
| |
+----------------------------------------------------------------------------+

Hamish wrote:

I have imported both the vector and the raster with minimum trouble, since the first has a .prj file which GRASS recognized, and the second

had a detailed import tutorial in the cygwin-GRASS install page. However, despite the fact the the two layers supposedly have the same NAD83 projection info (I did not enter anything manually), when I try
to display them together on a GRASS monitor, the state layout is in
the middle of the ocean compared to the elevation raster.

Is there a way to overlap the two layers or to chek for a coordinate
error?

r.info and v.info will give you the map extents, you should then be able
to figure out which is correct.

Also try "d.where -w" to give you a lat/lon value for different parts of your map which you can compare to the map on the wall.

Hamish

>
What happens if a big asteroid hits Earth? Judging from realistic simulations involving a sledge hammer and a common laboratory frog, we can assume it will be pretty bad.
- Dave Barry

-------------------------------------------------------------
This message has been scanned by Postini anti-virus software.

Luigi_Ponti · November 15, 2005, 11:03pm

Dear Ian,

the original metadata of my state outline shapefile is:

PROJCS["unnamed",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG",7019]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG",6269]],PRIMEM["Greenwich",0,AUTHORITY["EPSG",8901]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG",9108]],AXIS["Lat","NORTH"],AXIS["Long","EAST"],AUTHORITY["EPSG",4269]],PROJECTION["Albers"],PARAMETER["standard_parallel_1",34],PARAMETER["standard_parallel_2",40.5],PARAMETER["latitude_of_center",0],PARAMETER["Central_Meridian",-120],PARAMETER["False_Easting",0],PARAMETER["False_Northing",-4000000],UNIT["m",1]]

Based on this, following direction in the GRASS/cygwin installation page
http://geni.ath.cx/grass.html#toc15

I have created a dummy vector for reprojection of the raster DEM of western USA - e10g tile. Then I created a new location from EPSG code to match the vector .shp state outline file, and I have imported both the dummy raster DEM via "v.proj" and the vector state outline via "v.in.ogr".

Using "gdalinfo" on the original raster file (e10g), originated the following output:

gdalinfo c:/cygwin/home/andy/DEM_data/e10g
Driver: EHdr/ESRI .hdr Labelled
Size is 10800, 6000
Coordinate System is `'
Origin = (-180.000000,50.000000)
Pixel Size = (0.00833333,-0.00833333)
Corner Coordinates:
Upper Left (-180.0000000, 50.0000000)
Lower Left (-180.0000000, 0.0000000)
Upper Right ( -90.0000000, 50.0000000)
Lower Right ( -90.0000000, 0.0000000)
Center (-135.0000000, 25.0000000)
Band 1 Block=10800x1 Type=Int16, ColorInterp=Undefined
NoData Value=-500

Do you see anything wrong?

In addition, I tried to import v.in.e00 (after adding the required .exe file in the bin folder of GRASS) but I have not been successful yet.

Regards and thanks,

Luigi

On Mon, 14 Nov 2005 18:31:43 -0800
Ian MacMillan <ian_macmillan@umail.ucsb.edu> wrote:

Luigi, if you look at r.info and v.info, you will only find the projection information of your location, not necessarily the projection/datum/etc. information of the original data (I believe). Do you have the original metadata? If not, try running gdalinfo on your original data files to determine their original projection info.

man gdalinfo
will give you the correct syntax.

G'luck
Ian

On Nov 14, 2005, at 2:06 PM, Luigi Ponti wrote:

Thanks Amish,

looking at v.info and r.info, it seems the two maps have the same projection (see below). However, when looking at the output "d.where -w" for a sample map point (the San Francisco Bay), the longitude is very different in the two maps:

                                                  WGS84 Co-ordinates
            EAST: NORTH: LON:
            -191114.6346037 -24240.66895625 143:03:58.588289W 29:24:27.1
(San Francisco Bay on the California State vector layer - WRONG)

1839139.30311877 663981.00484798 122:19:50.783112W 37:57:48.1
(San Francisco Bay on the Elevation Model raster layer - CORRECT)

Do you think this could be due to an error in the import process or should I just look for another vector map?

Thanks again,

Luigi

v.info map=state layer=1
+---------------------------------------------------------------------- ------+
| Layer: state Organization:
                    |
| Mapset: luigi Source Date:
                     |
| Location: Cal_nad83-zone3 Name of creator: |
| Database: /home/andy
                                                      |
| Title:

         |
| Map Scale: 1:1
                                                           |
| Map format: native
                                                        |
|---------------------------------------------------------------------- ------|
| Type of Map: Vector (level: 2)
                                         |
|

                |
| Number of points: 0 Number of areas: 51 |
| Number of lines: 0 Number of islands: 48 |
| Number of boundaries: 59 Number of faces: 0 |
| Number of centroids: 51 Number of kernels: 0 |
|

                |
| Map is 3D: 0
                                               |
| Number of dblinks: 1
                                               |
|

                |
| Projection: Lambert Conformal Conic (zone 0)
                            |
| N: 922206.900 S: -701781.905
                                |
| E: 1043412.917 W: -374425.808
                              |
| B: 0.000 T: 0.000
                                      |
|

                |
| Digitize threshold: 0.00000
                                             |
| Comments:
                                                               |
|

                |
+---------------------------------------------------------------------- ------+

r.info map=e10g@PERMANENT
+---------------------------------------------------------------------- ------+
| Layer: e10g Date: Tue Nov 8 17:25:24 2005 |
| Mapset: PERMANENT Login of Creator: andy |
| Location: Cal_nad83-zone3
                                                 |
| DataBase: /home/andy
                                                      |
| Title: ( e10g )
                                                       |
|---------------------------------------------------------------------- ------|
|

                |
| Type of Map: raster Number of Categories: 5443 |
| Data Type: CELL
                                                      |
| Rows: 7192
                                                      |
| Columns: 11628
                                                     |
| Total Cells: 83628576
                                                  |
| Projection: Lambert Conformal Conic (zone 0)
                       |
| N: 3353466 S: -3838238 Res: 999.95884316 |
| E: 6073284 W: -5555183 Res: 1000.04016168 |
| Range of data: min = -84 max = 5443
                                 |
|

                |
| Data Source:
                                                            |
|

                |
|

                |
|

                |
| Data Description:
                                                       |
| generated by r.proj
                                                    |
|

                |
|

                |
+---------------------------------------------------------------------- ------+

Hamish wrote:

I have imported both the vector and the raster with minimum trouble, since the first has a .prj file which GRASS recognized, and the second

had a detailed import tutorial in the cygwin-GRASS install page. However, despite the fact the the two layers supposedly have the same NAD83 projection info (I did not enter anything manually), when I try
to display them together on a GRASS monitor, the state layout is in
the middle of the ocean compared to the elevation raster.

Is there a way to overlap the two layers or to chek for a coordinate
error?

r.info and v.info will give you the map extents, you should then be able
to figure out which is correct.

Also try "d.where -w" to give you a lat/lon value for different parts of your map which you can compare to the map on the wall.

Hamish

>
What happens if a big asteroid hits Earth? Judging from realistic simulations involving a sledge hammer and a common laboratory frog, we can assume it will be pretty bad.
- Dave Barry

-------------------------------------------------------------
This message has been scanned by Postini anti-virus software.

INFINITO ADSLFLAT 4 MEGA: SOLO 27,90 EURO AL MESE IVA INCLUSA IP STATICO, BANDA GARANTITA 256Kbps ANTIVIRUS E FIREWALL INCLUSI NEL PREZZO

http://adsl.infinito.it

[GRASSLIST:8825] GRASS 6 in a nutshell tutorial

============================================================================ System Information

============================================================================
System Information