Hi,
I note the following in MAP_INSTALL, so was not sure if this should be a bug or feature request:
81 NOTE: Map projection support in GDL is a work in progress. There are
82 BUGS!! These will be fixed as time permits. Additional projections
83 will also be added in the same manner. Bug reports and projection
84 requests are appreciated.
We are trying to access the MAP_PROJ_INIT UTM GCTP projection from GDL on Linux, but are getting an error regarding the loading of MAP_PROJ_GCTP_FORINIT, which appears to only be internally accessible to IDL.
To the best of my knowledge, all instructions in MAP_INSTALL have been followed, however, the list of required *.pro files has been copied over from an IDL 8.1 installation given that the link provided to downloadable files in the MAP_INSTALL instructions only contains a subset of the files.
We have the following example function:
pro maptest
zoneID=-54
maxLongi=140.96964
minLongi=139.718946172
maxLat=-36.9494155233
minLat=-38.3109611892
ctrLong=(minLongi+maxLongi)/2
ctrLat=(minLat+maxLat)/2
UTM_Struct=MAP_PROJ_INIT('UTM',/GCTP, ELLIPSOID=12,$
CENTER_LONGITUDE=ctrLong,CENTER_LATITUDE=ctrLat,ZONE=zoneID)
end
When called in GDL we get:
GDL> maptest
% Compiled module: MAPTEST.
% Compiled module: MAP_PROJ_INIT.
% Compiled module: MAP_PROJ_INIT_COMMON.
% Compiled module: STRSPLIT.
% Compiled module: UNIQ.
% MAP_PROJ_INIT: Procedure not found: MAP_PROJ_GCTP_FORINIT
% Error occurred at: MAP_PROJ_INIT 1378 /apps/gdl/0.9.3/share/gnudatalanguage/lib/gshhs/map_proj_init.pro
% MAPTEST 9 maptest.pro
% $MAIN$
% Execution halted at: MAPTEST 9 maptest.pro
Doing a recursive grep on the GDL install directory indicates that only map_proj_init.pro contains the string "MAP_PROJ_GCTP_FORINIT".
Doing the same on the 8.1 IDL install directory shows that the string is present in the libidl.8.1.so file, but I cannot see it listed in the symbols list when running nm libidl.8.1.so, so I'm assuming IDL links this in at compile time or something similar.
Is there anything we might be doing wrong? Or is this path of functionality not currently supported in GDL. Are there any thoughts on if, how, or when this projection might be made available.
Regards,
Luke.
Hi,
The support for map projections is still limited to a subset of what is in IDL6, and that may even be broken in places! GDL is far from competing with idl8.1 especially since it introduced all the GCTP projections. All of them are in the libproj library but the connexion between the libproj and gdl syntax must be done manually with some educated knowledge of the projection's particulars. Not terribly difficult but one has to carefully read both manuals (idl and libproj) and test all cases.
for the record, only the following projections have a chance today to work partially:
Stereographic Projection: (1)
Orthographic Projection: (2)
Lambert Conformal Conic: (3)
Lambert Equal Area Conic: (4)
Gnomonic: (5)
Azimuthal Equidistant: (6)
Satellite: (7) (add 28/07/2008)
Cylindrical Equidistant: (8)
Mercator: (9)
Aitoff: (12)
Robinson: (17) (add 28/07/2008)
Help welcome, of course, do not hesitate to propose patches!
Thanks heaps. My colleagues implemented a work around by SPAWNing out to another GCTP application and marshaling the result data into IDL/GDL representation. But they will be pleased they won't have to do that again.
Well test it and only Rey if there are problems.
Regards,
Luke Domanski
On 14/06/2014, at 7:17 AM, "giloo" gilles-duvert@users.sf.net wrote:
Related
Bugs:
#588OK, everything works now. Will be fixed in 0.9.6 and in cvs soon.
fixed in cvs. see MAP_INSTALL for new instructions for mapping