[GRASSLIST:4782] trouble with v.proj
bfdi533
bfdi533 at cbacc-security3.bus.utexas.edu
Thu Nov 11 15:01:22 EST 2004
I have a TIGER line file that was imported into a LL location with an
output from g.projinfo here:
-----------------------------------------------------------
PROJ_INFO file:
name: Latitude-Longitude
datum: wgs84
dx: 0.000000
dy: 0.000000
dz: 0.000000
proj: ll
ellps: wgs84
-----------------------------------------------------------
PROJ_UNITS file:
unit: degree
units: degrees
meters: 1.0
-----------------------------------------------------------
I then input another roads file into a location in TX State Plane format
with output from g.projinfo here:
-----------------------------------------------------------
PROJ_INFO file:
name: State Plane
proj: lcc
a: 0.63782064e+07
es: 0.6768657997291094e-02
x_0: 0.6096012192024384e+06
y_0: 0
lon_0: 100d20'w
lat_0: 29d40'n
lat_1: 31d53'n
lat_2: 30d7'n
-----------------------------------------------------------
PROJ_UNITS file:
unit: meter
units: meters
meters: 1.0
-----------------------------------------------------------
I then tried v.proj from the LL location INTO the State Plane location
and did a d.vect on both sets of roads. They do not line up.
Here is the region from the STP imported vector file:
GRASS:~ > g.region vect=streetcl
GRASS:~ > d.erase
GRASS:~ > g.region -l
long: 48.34520 lat: 86.33708 (north/west corner)
long: 51.00098 lat: 85.50330 (north/east corner)
long: 41.70189 lat: 85.31135 (south/east corner)
long: 38.11301 lat: 86.14387 (south/west corner)
rows: 512
cols: 514
Center Longitude: 44:47:24.969634E [44.79027]
Center latitude: 85:49:26.046164N [85.82390]
And here is the region that I get from running the same commands on the
vector file that was imported with v.proj:
GRASS:~ > g.region vect=tx_travis
GRASS:~ > d.erase
GRASS:~ > g.region -l
long: -98.17640 lat: 30.64542 (north/west corner)
long: -97.35039 lat: 30.62892 (north/east corner)
long: -97.36951 lat: 30.00828 (south/east corner)
long: -98.19023 lat: 30.02468 (south/west corner)
rows: 150
cols: 147
Center Longitude: 97:46:17.870455W [-97.77163]
Center latitude: 30:19:36.571973N [30.32683]
It seems quite likely that I have something wrong with the STP location
but I have recreated it more than once using the right info and then
imported data (arcinfo shape files) that say they are in the STP coords
that I created. But, the above just does not line up right and the
coords look to be correct for the projected vector file, actually.
Any ideas on what I am doing wrong?
Ed
More information about the grass-user
mailing list