[GRASSLIST:10027] projection problems

Kirk R. Wythers kwythers at umn.edu
Fri Jan 27 12:52:04 EST 2006


It seems that I could use some very basic explanation concerning  
locations and projections. I am trying to read several shape files  
into a single location. However, I am getting a projection error on  
some of them.

g.region -p for the location is:

GRASS 6.1.cvs (minnesota_utm):~ > g.region -p
projection: 1 (UTM)
zone:       15
datum:      nad83
ellipsoid:  grs80
north:      5472415.17999963
south:      4816310.42
west:       189788.954
east:       761663.27300056
nsres:      0.99999963
ewres:      1.00000056
rows:       656105
cols:       571874

and v.info for the state outline shape file (which read in fine with  
v.in.ogr) is:

GRASS 6.1.cvs (minnesota_utm):~ > v.info mn_state
+----------------------------------------------------------------------- 
-----+
| Layer:    mn_state                        
Organization:                     |
| Mapset:   pls                            Source  
Date:                      |
| Location: minnesota_utm                  Name of  
creator:                  |
| Database: /Volumes/disk1/home1/kwythers/ 
grassdata                          |
|  
Title:                                                                   
    |
| Map Scale:   
1:1                                                            |
| Map format:  
native                                                         |
|----------------------------------------------------------------------- 
-----|
|   Type of Map:  Vector (level:  
2)                                          |
|                                                                        
      |
|   Number of points:       0               Number of areas:       
1          |
|   Number of lines:        0               Number of islands:     
1          |
|   Number of boundaries:   1               Number of faces:       
0          |
|   Number of centroids:    1               Number of kernels:     
0          |
|                                                                        
      |
|   Map is 3D:               
0                                                |
|   Number of dblinks:       
1                                                |
|                                                                        
      |
|   Projection: UTM (zone  
0)                                                 |
|            N: 5472414.182    S:  
4816310.425                                |
|            E: 761662.273    W:  
189788.954                                  |
|            B: 0.000         T:  
0.000                                       |
|                                                                        
      |
|   Digitize threshold:  
0.00000                                              |
|    
Comments:                                                                
  |
|                                                                        
      |
+----------------------------------------------------------------------- 
-----+

Although I do notice that the projection reads UTM (zone 0) which  
seems odd to me... it should be zone 15.

Now for the really fun part...

If I try and read another shape file into the location, I get the  
projection error about not matching the current location:

GRASS 6.1.cvs (minnesota_utm):~ > v.in.ogr dsn=~/projects/mn_dnr/ 
shape_files/pls_corners_btrees/county/cook/ output=cook_test
ERROR: Projection of dataset does not appear to match current location.

        LOCATION PROJ_INFO is:
        name: UTM
        datum: nad83
        towgs84: 0.000,0.000,0.000
        proj: utm
        ellps: grs80
        a: 6378137.0000000000
        es: 0.0066943800
        f: 298.2572221010
        zone: 15

        cellhd.proj = 0 (unreferenced/unknown)

        You can use the -o flag to v.in.ogr to override this check.
        Consider to generate a new location with 'location' parameter  
from
        input data set.

If I use the -o flag (or location=some_new_location) the file comes  
in as projection (0) x,y

I asked one of our ARC/GIS folks downstairs to see if he could import  
the shape files and lay them on top of each other and he did not have  
any trouble. I am wondering if I have set something very basic up wrong.

Any suggestions for diagnosing this problem and getting these shape  
files to match up would be much appreciated.

Thanks




More information about the grass-user mailing list