[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
Hi:
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
PROJDIR=/nas/uncch/depts/cep/emc/lran/soft/proj-4.5.0/exe/lib
PROJ_H=/nas/uncch/depts/cep/emc/lran/soft/proj-4.5.0/src/projects.h
# Library paths
LIBPATH=-L. -L$(PROJDIR)
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:
PROJ
args=+proj=lcc+datum=NAD83+lat_1=30.000000+lat_2=60.000000+lon_0=-100.000000+lat_0=40.000000+units=m
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.
Thanks,
LR
More information about the Proj
mailing list