[GRASSLIST:5459] geogrid in ps.map

Ian MacMillan ian_macmillan at umail.ucsb.edu
Thu Jan 20 14:44:32 EST 2005


Hi all, I am using grass5.7 on osx (May 2004 or so) and I am trying to 
make a world mercator map using ps.map.  I would like to add a geogrid, 
however the output only gives latitudinal lines, no longitudes.  No 
errors were returned.  The command I gave was

ps.map input=- output=test.eps -e
geogrid 30 d
color black
numbers 1
end
end

g.region -p yields

projection: 99 (Mercator)
zone:       0
datum:      ** unknown (default: WGS84) **
ellipsoid:  wgs84
north:      11714654.84
south:      -11714654.84
west:       -20037508.34
east:       20037508.34
nsres:      3697.80771465
ewres:      3710.64969259
rows:       6336
cols:       10800


Alternatively, I have also tried to import (using v.proj) a worldwide 
vector grid with 30 degree spacing, however every attempt to import 
yields in errors such as (using grass 5.3):

Input Projection Parameters: +proj=latlong +a=6378137 
+es=0.006694379990141316 +towgs84=0.000,0.000,0.000

Output Projection Parameters: +proj=merc +lat_ts=0.0000000000 
+lon_0=180.0000000000 +k_0=1.0000000000 +a=6378137 
+es=0.006694379990141316

Creating dig file...
pj_transform() failed
cause: latitude or longitude exceeded limits
Error in pj_do_proj

OR (using 5.7)

v.proj input=grid30 location=latlong_wgs84

Input Projection Parameters: +proj=latlong +a=6378137 +rf=298.257223563 
+no_defs +towgs84=0.000,0.000,0.000
Input Unit Factor: 1

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

Creating vector file...
Building topology ...
Registering lines:

At this point v.proj stops running.  Trying to display the vector with 
yields:

GRASS 5.7.0:~ > d.vect grid30
WARNING: coor files of vector 'grid30 at ianmacmillan' is larger than it
          should be (14312 bytes excess).
WARNING: Cannot display areas, topology not available

I think there is a problem with worldwide mercator projections, but I 
can't figure what it is.  Displaying other worldwide data (I can't 
remember how I got it into this location in the first place) shows that 
the region is appropriate for a world map.  It is hard to tell for sure 
though since units are measured in meters.

Any ideas?

Thanks,
Ian




More information about the grass-user mailing list