[GRASS-dev] Need to update g.setproj

Sun Oct 1 01:46:22 EDT 2006

Either this or Hamish's suggestion would work from my vantage--whichever is
easier to do with the g.setproj code.

There would need to be some tables with user-readable names from which to
parse correct values. But, as you show below, that should be quite doable.

Michael
__________________________________________
Michael Barton, Professor of Anthropology
School of Human Evolution & Social Change
Center for Social Dynamics & Complexity
Arizona State University

phone: 480-965-6213
fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

> From: Glynn Clements <glynn at gclements.plus.com>
> Date: Sun, 1 Oct 2006 03:56:59 +0100
> To: Michael Barton <michael.barton at asu.edu>
> Cc: Paul Kelly <paul-grass at stjohnspoint.co.uk>, GRASS-DEV
> <grass-dev at grass.itc.it>
> Subject: Re: [GRASS-dev] Need to update g.setproj
> 
> 
> Michael Barton wrote:
> 
>>> Is there any real reason why GRASS needs to comprehend the projection
>>> parameters, i.e. to have lib/gis/geo_init.c or all of that boilerplate
>>> in g.setproj?
>>> 
>>> IOW, can projection parameters not be treated as opaque key/value
>>> pairs which are simply passed straight through to PROJ, with any UI
>>> consisting of generic code using an external database[1] of projection
>>> parameters?
>>> 
>>> [1] By which, I mean a text file in an easy-to-parse format.
>>> 
>>> It seems to me that all of that boilerplate could be replaced by a
>>> list of projections (+proj= name, description), a database of
>>> parameters (+<opt>= name, description, bool/int/float), and a matrix
>>> of "projection X accepts/requires parameter Y with default value Z".
>> 
>> This seems like it would be a big help in designing a GUI wrapper for these
>> functions.
> 
> Converting lib/gis/geo_init.c and g.setproj/main.c to text files was
> simple enough; those are attached.
> 
> AFAICT, it should be straightforward to replace the boilerplate in
> g.setproj with code which uses the tables.
> 
> In the process, I noticed this in g.setproj/main.c:
> 
> case LAT3:
> sprintf(tmp_buff, "%.10f", LLSTUFF[i]);
> G_set_key_value("lat_2", tmp_buff, out_proj_keys);
>                 ^^^^^
> 
> -- 
> Glynn Clements <glynn at gclements.plus.com>
> 
> ALPHA:alpha:Azimuth angle at Cartesian Origin
> AZIM:azi:Azimuth Angle of Tilt in Decimal degrees
> HEIGH:h:Height of Viewing Point in Meters
> KFACT:k_0:Scale Factor at the Central Meridian
> LAT0:lat_0:Central Parallel
> LAT1:lat_1:First Standard Parallel
> LAT2:lat_2:Second Standard Parallel
> LAT3:lat_3:Third Standard Parallel
> LATB:lat_b:Angular Distance from Tangency Point
> LATTS:lat_ts:Latitude of True Scale
> LON0:lon_0:Central Meridian
> LON1:lon_1:First Standard Meridian
> LON2:lon_2:Second Standard Meridian
> LON3:lon_3:Third Standard Meridian
> LONC:lon_c:Longitude of Cartesian Origin
> LOTSA:lotsa:LOTSA
> MFACT:m:m factor
> MSFACT:M:M factor
> NFACT:n:n factor
> NOCUT:no_cut:Both Hemispheres
> NOROT:no_rot:Suppress Rotation
> NOSKEW:ns:Suppress Skew
> NOUOFF:no_uoff:Suppress Offset from Pre-Rotated Axis
> OLATP:o_lat_p:Latitude of New Pole
> OLONP:o_lon_p:Longitude of New Pole
> QFACT:q:q factor
> ROTCONV:rot_conv:Origin Convergence Angle
> SNUM:lsat:Satellite Number
> SOUTH:south:South Hemisphere
> SPATH:path:Satellite Path Number
> THETA:theta:Theta Angle
> TILT:tilt:Tilt Angle in Decimal Degrees
> WFACT:W:W factor
> X0:x_0:False Easting
> Y0:y_0:False Northing
> ZONE:zone:Projection Zone
> AEA:lat_0=ask,23.0;lat_1=ask,29.5;lat_2=ask,45.5;lon_0=ask,-96.0;x_0=ask,0.0;y
> _0=ask,0.0
> AEQD:lat_0=ask,0.0;lon_0=ask,20.0
> AIRY:lat_0=ask,0.0;lat_b=ask,90.0;lon_0=ask,20.0;no_cut=ask,nodfl
> AITOFF:lat_0=ask,0.0;lon_0=ask,20.0
> ALSK:lat_0=noask,64.0;lon_0=noask,-152.0
> APIAN:lat_0=ask,0.0;lon_0=ask,20.0
> AUGUST:lat_0=ask,0.0;lon_0=ask,20.0
> BACON:lat_0=ask,0.0;lon_0=ask,20.0
> BIPC:lat_0=ask,0.0;lon_0=ask,-90.0;ns=ask,nodfl
> BOGGS:lat_0=ask,0.0;lon_0=ask,20.0
> BONNE:lat_1=ask,40.0;lon_0=ask,20.0
> CASS:lat_0=ask,0.0;lon_0=ask,20.0
> CC:lat_0=ask,0.0;lon_0=ask,20.0
> CEA:lat_0=ask,0.0;lat_ts=ask,0.0;lon_0=ask,20.0
> CHAMB:lat_0=ask,55.0;lat_1=ask,40.0;lat_2=ask,20.0;lat_3=ask,35.0;lon_0=ask,20
> .0;lon_1=ask,5.0;lon_2=ask,55.0;lon_3=ask,65.0
> COLLG:lat_0=ask,0.0;lon_0=ask,20.0
> CRAST:lat_0=ask,0.0;lon_0=ask,20.0
> DENOY:lat_0=ask,0.0;lon_0=ask,20.0
> ECK1:lat_0=ask,0.0;lon_0=ask,20.0
> ECK2:lat_0=ask,0.0;lon_0=ask,20.0
> ECK3:lat_0=ask,0.0;lon_0=ask,20.0
> ECK4:lat_0=ask,0.0;lon_0=ask,20.0
> ECK5:lat_0=ask,0.0;lon_0=ask,20.0
> ECK6:lat_0=ask,0.0;lon_0=ask,20.0
> EQC:lat_0=ask,0.0;lat_ts=ask,0.0;lon_0=ask,0.0
> EQDC:lat_0=ask,40.0;lat_1=ask,20.0;lat_2=ask,60.0;lon_0=ask,20.0
> EULER:lat_0=ask,55.0;lat_1=ask,45.0;lat_2=ask,65.0;lon_0=ask,20.0
> FAHEY:lat_0=ask,0.0;lon_0=ask,20.0
> FOUC:lat_0=ask,0.0;lon_0=ask,20.0
> FOUC_S:lat_0=ask,0.0;lon_0=ask,20.0
> GALL:lat_0=ask,0.0;lon_0=ask,20.0
> GINS8:lat_0=ask,0.0;lon_0=ask,20.0
> GNOM:lat_0=ask,0.0;lon_0=ask,20.0
> GN_SINU:lat_0=ask,0.0;lon_0=ask,20.0;m=ask,1.0;n=ask,1.0
> GOODE:lat_0=ask,0.0;lon_0=ask,20.0
> GS48:lat_0=noask,45.0;lon_0=noask,-120.0
> GS50:lat_0=noask,45.0;lon_0=noask,-120.0
> HAMMER:M=ask,1.0;W=ask,0.5;lat_0=ask,0.0;lon_0=ask,20.0
> HATANO:lat_0=ask,0.0;lon_0=ask,20.0
> 
IMW_P:lat_0=ask,0.0;lat_1=ask,20.0;lat_2=ask,60.0;lon_0=ask,20.0;lon_1=ask,20.>
0
> KAV5:lat_0=ask,0.0;lon_0=ask,20.0
> KAV7:lat_0=ask,0.0;lon_0=ask,20.0
> LABRD:azi=noask,18.9;k_0=noask,0.9995;lat_0=noask,18.9;lon_0=noask,46.43720833
> 3;x_0=noask,400000.0;y_0=noask,800000.0
> LAEA:lat_0=ask,55.0;lon_0=ask,20.0;x_0=ask,0.0;y_0=ask,0.0
> LAGRNG:W=ask,2.0;lat_0=ask,0.0;lat_1=ask,0.0;lon_0=ask,20.0
> LARR:lat_0=ask,0.0;lon_0=ask,20.0
> LASK:lat_0=ask,0.0;lon_0=ask,20.0
> LCC:lat_0=ask,23.0;lat_1=ask,33.0;lat_2=ask,45.0;lon_0=ask,-96.0;x_0=ask,0.0;y
> _0=ask,0.0
> LEAC:lat_0=ask,55.0;lat_1=ask,0.0;lon_0=ask,20.0;south=ask,nodfl
> LEE_OS:lat_0=noask,-10.0;lon_0=noask,-165.0
> LOXIM:lat_0=ask,0.0;lat_1=ask,45.0;lon_0=ask,20.0
> LSAT:lat_0=ask,0.0;lon_0=ask,20.0;lsat=ask,1;path=ask,1
> MBTFPP:lat_0=ask,0.0;lon_0=ask,20.0
> MBTFPQ:lat_0=ask,0.0;lon_0=ask,20.0
> MBTFPS:lat_0=ask,0.0;lon_0=ask,20.0
> MBT_FPS:lat_0=ask,0.0;lon_0=ask,20.0
> MBT_S:lat_0=ask,0.0;lon_0=ask,20.0
> MERC:k_0=ask,1.0;lat_ts=ask,0.;lon_0=ask,-96.0
> MILL:lat_0=ask,0.0;lon_0=ask,20.0
> MIL_OS:lat_0=noask,18.0;lon_0=noask,20.0
> MOLL:lat_0=ask,0.0;lon_0=ask,20.0
> MPOLY:lat_0=ask,0.0;lat_1=ask,-20.0;lat_2=ask,20.0;lon_0=ask,20.0;lotsa=ask,no
> dfl
> MURD1:lat_0=ask,0.0;lat_1=ask,-20.0;lat_2=ask,20.0;lon_0=ask,20.0
> MURD2:lat_0=ask,0.0;lat_1=ask,-20.0;lat_2=ask,20.0;lon_0=ask,20.0
> MURD3:lat_0=ask,0.0;lat_1=ask,-20.0;lat_2=ask,20.0;lon_0=ask,20.0
> NELL:lat_0=ask,0.0;lon_0=ask,20.0
> NELL_H:lat_0=ask,0.0;lon_0=ask,20.0
> NICOL:lat_0=ask,0.0;lon_0=ask,20.0
> NSPER:h=ask,50000000.0;lat_0=ask,55.0;lon_0=ask,20.0
> NZMG:lat_0=noask,-41.0;lon_0=noask,173.0;x_0=noask,2510000.0;y_0=noask,6023150
> .0
> OB_TRAN:lat_0=ask,0.0;lon_0=ask,0.0;o_lat_p=ask,90.0;o_lon_p=ask,0.0
> OCEA:lat_0=ask,0.0;lat_1=ask,-45.0;lat_2=ask,45.0;lon_0=ask,20.0;lon_1=ask,-20
> .0;lon_2=ask,60.0
> OEA:lat_0=ask,0.0;lon_0=ask,20.0;m=ask,1.0;n=ask,1.0;theta=ask,0.0
> OMERC:k_0=ask,1.0;lat_0=ask,0.0;lat_1=ask,-45.0;lat_2=ask,45.0;lon_0=ask,20.0;
> lon_1=ask,-40.0;lon_2=ask,40.0;no_rot=ask,nodfl;no_uoff=ask,nodfl;rot_conv=ask
> ,nodfl
> ORTEL:lat_0=ask,0.0;lon_0=ask,20.0
> ORTHO:lat_0=ask,0.0;lon_0=ask,20.0
> PCONIC:lat_0=ask,0.0;lat_1=ask,33.0;lat_2=ask,45.0;lon_0=ask,20.0
> POLY:lat_0=ask,0.0;lon_0=ask,-90.0
> PUTP1:lat_0=ask,0.0;lon_0=ask,20.0
> PUTP2:lat_0=ask,0.0;lon_0=ask,20.0
> PUTP3:lat_0=ask,0.0;lon_0=ask,20.0
> PUTP3P:lat_0=ask,0.0;lon_0=ask,20.0
> PUTP4P:lat_0=ask,0.0;lon_0=ask,20.0
> PUTP5:lat_0=ask,0.0;lon_0=ask,20.0
> PUTP5P:lat_0=ask,0.0;lon_0=ask,20.0
> PUTP6:lat_0=ask,0.0;lon_0=ask,20.0
> PUTP6P:lat_0=ask,0.0;lon_0=ask,20.0
> QUA_AUT:lat_0=ask,0.0;lon_0=ask,20.0
> ROBIN:lat_0=ask,0.0;lon_0=ask,20.0
> RPOLY:lat_0=ask,0.0;lat_ts=ask,0.0;lon_0=ask,20.0
> SINU:lat_0=ask,0.0;lon_0=ask,20.0
> SOMERC:k_0=noask,1.0;lat_0=noask,46.952405556;lon_0=noask,7.4395833333;x_0=noa
> sk,600000.0;y_0=noask,200000.0
> STERE:k_0=ask,1.0;lat_0=ask,0.0;lat_ts=ask,0.0;lon_0=ask,20.0
> TCC:lat_0=ask,0.0;lon_0=ask,20.0
> TCEA:k_0=ask,1.0;lat_0=ask,0.0;lon_0=ask,20.0
> TISSOT:lat_0=ask,0.0;lat_1=ask,-30.0;lat_2=ask,45.0;lon_0=ask,20.0
> TMERC:k_0=ask,1.0;lat_0=ask,23.0;lon_0=ask,-96.0;x_0=ask,0.0;y_0=ask,0.0
> TPEQD:lat_0=ask,0.0;lat_1=ask,-45.0;lat_2=ask,45.0;lon_0=ask,20.0;lon_1=ask,-2
> 0.0;lon_2=ask,60.0
> TPERS:azi=ask,0.0;h=ask,10000.0;lat_0=ask,0.0;lon_0=ask,20.0;tilt=ask,0.0
> UPS:south=ask,nodfl
> URM5:alpha=ask,0.0;lat_0=ask,0.0;lon_0=ask,20.0;n=ask,1.0;q=ask,1.0
> URMFPS:lat_0=ask,0.0;lon_0=ask,20.0;n=ask,1.0
> UTM:south=ask,nodfl;zone=ask,nodfl
> VANDG:lat_0=ask,0.0;lon_0=ask,20.0
> VANDG2:lat_0=ask,0.0;lon_0=ask,20.0
> VANDG3:lat_0=ask,0.0;lon_0=ask,20.0
> VANDG4:lat_0=ask,0.0;lon_0=ask,20.0
> WAG1:lat_0=ask,0.0;lon_0=ask,20.0
> WAG2:lat_0=ask,0.0;lon_0=ask,20.0
> WAG3:lat_0=ask,0.0;lon_0=ask,20.0
> WAG4:lat_0=ask,0.0;lon_0=ask,20.0
> WAG5:lat_0=ask,0.0;lon_0=ask,20.0
> WAG6:lat_0=ask,0.0;lon_0=ask,20.0
> WAG7:lat_0=ask,0.0;lon_0=ask,20.0
> WEREN:lat_0=ask,0.0;lon_0=ask,20.0
> WINK1:lat_0=ask,0.0;lat_ts=ask,0.0;lon_0=ask,20.0
> WINK2:lat_0=ask,0.0;lat_1=ask,0.0;lon_0=ask,20.0
> WINTRI:lat_0=ask,0.0;lat_1=ask,0.0;lon_0=ask,20.0