[GRASS5] PROJ_INFO and pj_get_kv() and r.in.gdal

Markus Neteler neteler at itc.it
Sat Oct 6 02:54:46 EDT 2001

Hi all,

while working with/on r.in.gdal I have realized that three
entries are not written to PROJ_INFO when creating a new
location from a geocoded dataset. These entries are
"name", "a" and "es". Example:

name: Transverse Mercator
datum: rome40
dx: -225.000000
dy: -65.000000
dz: 9.000000
proj: tmerc
ellps: international
a: 6378388.0000000000
es: 0.0067226700
f: 297.0000000000
lat_0: 0.0000000000
lon_0: 9.0
k_0: 0.9996000000
x_0: 1500000.0000000000

If "name", "a" and "es" the r.proj, v.proj and s.proj fail because
the function pj_get_kv() fails. This function is in

In my opinion there is no need to change r.in.gdal ( see wkt_to_grass()
function there which defines the proj_info and proj_units struct,
which is later used by G__make_location() ) because
"name", "a" and "es"
are redundant. The needed values are "proj" and "ellps".

Question: Can we simplify pj_get_kv() or do we have to update
wkt_to_grass() in src/raster/r.in.gdal/main.c

If you look into src/raster/r.in.gdal/main.c you see my additions
to optionally re-project GCPs in a SAR image from lat/long to 
a user defined target projection. But since pj_get_kv() fails
here due to the missing values, I had to hard-code lat/long as
input "projection" (which might not be true for all SAR data).

Comments would be welcome. I have spend some time to update wkt_to_grass()
but didn't manage (a bit too complex for me...)

Thanks for hints,


More information about the grass-dev mailing list