[GRASS-user] Difficlty with projections - newbie

pmarc paulomarcondes at gmail.com
Tue Jun 20 22:40:12 EDT 2006


Hi,
before anyone question my newbie status, I wish to reaffirm it!
I've been reading this list since a long time but never got my hands
dirty working with GRASS.

well, back to the point. or polygon as is the case.

So, The following vector is supposed to be in polyconic projection, as
used in Brasil. The following info came with the vectors:

- Ellipsoid: UGGI 67
- Semi-major axis= 6378.160 km
- Semi-minor axis = 6356.775 km
- Excentricity = 0.0818201799960599
- Flattening = 0.00335289186923715
- Planimetric Datum: South American 1969 - SAD69
- Altimetric Datum: Imbituba

- Origin Latitude:  0º ( Equator);
- Origin Longitude : - 54º (54W).

>From what OGR tells me, I believe there is no projection information
on the shapefile:

paulo at onix:~/gis_data/ibge/mapa_indice$ ogrinfo -so
MALHA_MUNICIPAL_2000.shp layer MALHA_MUNICIPAL_2000
INFO: Open of `MALHA_MUNICIPAL_2000.shp'
using driver `ESRI Shapefile' successful.

Layer name: MALHA_MUNICIPAL_2000
Geometry: Polygon
Feature Count: 5726
Extent: (-2205855.200000, -3736283.700000) - (2615330.400000, 586413.000000)
Layer SRS WKT:
(unknown)
CODIGO_UF: String (25.0)
NOME_MUNIC: String (64.0)
GEOCODIGO_: String (7.0)
AREA_2000: Real (24.5)
LATITUDE_S: Real (24.5)
LONGITUDE_: Real (24.5)
ID: Integer (10.0)
mslink: Integer (10.0)
mapid: Integer (10.0)

What I want to do: add the projection information and then re-project
this to lat-lon in order to  overlay this on the GLCF shapefiles.

I need some help to glue this all together.

thanks you all.
-- 
Paulo Marcondes = PU1/PU2PIX
-22.915 -42.229 = GG87vc (http://www.amsat.org/cgi-bin/gridconv)
Debian GNU/Linux = http://rj.debianbrasil.org = http://www.debian.org




More information about the grass-user mailing list