[Proj] Problem in calling pj_transform with retutn x and y = nan

Limei Ran lran at unc.edu
Mon Nov 19 09:52:44 EST 2007


Could someone help me on solving the problem I have in calling 
pj_transform?  I am calling pj_transform to do the projection from our 
GIS software in C.  I obtained PROJ.4 and datum shift grids files.  I 
installed the software on our linux server and tested the software using 
cs2cs.  It works fine with the right output. 

I compiled our software with the built PROJ.4 library using the lines in 
our Makefile:
#set PROJ.4 libary

# Library paths
LIBS = $(LIBPATH) -lspatial -lm -lproj -static
allocator.exe: $(ALLOC)
        $(CC) -o $@ $(ALLOC) $(LIBS)

Our program is compiled without errors with PROJ.4 software. 
In my calling scripts, I set:
setenv PROJ_LIB /nas/uncch/depts/cep/emc/lran/soft/proj-4.5.0/exe/share/proj

The program gives the right projected x and y when from and to 
projections are different, but the datum systems are the same (both 
NAD27, NAD83, ot sphere).  But, if I give different datum for from and 
to projections, the print out of projected x and y = nan.  Here are the 
output lines:
Input projection:
PROJ args=+proj=lcc+datum=NAD27+lat_1=33+lat_2=45+lat_0=40+lon_0=-97

Output projection:

In Projection: x=1100.000 y=-11200.000    Out: x=nan y=nan
In Projection: x=-25000.000 y=25000.000  Out: x=nan y=nan

Basically, the data transformation parts did not give me the right 
result even though the program ran without any error.  There are many 
libary files under ./proj/lib directory.  I am wondering whether LIBS in 
my Makefile is correct.  Or, I need to set PROJ_LIB somewhere else.  I 
appreciate any help you might have.



More information about the Proj mailing list