[GRASS-user] Mercator projection

Steve Smith ssteve52 at yahoo.com
Mon Jan 14 23:18:27 EST 2008


I am attempting to reproject data from a lamber conformal conic (lcc) projection to a mercator projection. The reproject appears to go well. The output is listed below. However, when I attempt to display the data I get the message: The bounding box of the map is outside the current region. So I use g.region to change the bounds of the region to: 50, 20, -40, -100, but I still get the same error message. As a test, I created another Location, this time as Global LatLon, and after reprojecting the data, from my lcc projection, I was able to display the data, after changing the bounds of the region to the values above.  This makes me think that there is something wrong about my mercator projection. When I run v.info on my data the projection info is:

Projection: x,y (zone 0)                       
           N: 4370846.189    S: 3988531.316                                
           E: -8414703.985    W: -9288507.105                              
           B: 0.000         T: 0.000 

It does not say I have a mercator projection?  When I created the location, I selected merc as my projection, declined to select a datum(wgs84 was used by default), and then when it asked for an ellipsoid I selected sphere, with the following settings:

lat_ts=0, lon0=0, k=1, a=6378137.0, north=85.05112878, south=-85.05112978, west=-180, east=180, nsres=0.00000904, ewres=0.00000898

Steve
Below is the output from my reproject:
v.proj input=urbnarea location=NorthCarolina mapset=PERMANENT dbase=/var/www/html/mapdata/grassdata output=NCurbnarea

Input Projection Parameters: +proj=lcc +lat_0=33.7500000000 +lat_1=34.3333330000 +lat_2=36.1666600000 +lon_0=-79.0000000000 +x_0=609601.2200000000 +y_0=0.0000000000 +a=6378137 +rf=298.257222101 +no_defs +towgs84=0.000,0.000,0.000
Input Unit Factor: 1

Output Projection Parameters: +proj=merc +lat_ts=0.0000000000 +lon_0=0.0000000000 +k_0=1.0000000000 +a=6378137 +b=6378137 +no_defs
Output Unit Factor: 1

Creating vector file...
Building topology ...
Registering lines:    1000
1968 primitives registered      
Building areas: 


666 areas built      
183 isles built
Attaching islands: 

Attaching centroids: 

Topology was built.
Number of nodes     :   1485
Number of primitives:   1968
Number of points    :   0
Number of lines     :   0
Number of boundaries:   1311
Number of centroids :   657
Number of areas     :   666
Number of isles     :   183
Number of areas without centroid :   9

       
---------------------------------
Be a better friend, newshound, and know-it-all with Yahoo! Mobile.  Try it now.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20080114/2ba1f325/attachment.html


More information about the grass-user mailing list