[Mapserver-users] Globe View with Orthographic Projection

[Harald Wehr]

Thanks for the changes you made to the source code. I've made a cvs 
checkout for proj.4, gdal and mapserver but I fail to compile mapserver:

gcc -O2  -Wall -DIGNORE_MISSING_DATA  -DUSE_EPPL -DUSE_PROJ 
-DUSE_PROJ_API_H -DUSE_WMS_SVR -DUSE_WMS_LYR -DUSE_WFS_SVR -DUSE_WFS_LYR 
-DUSE_MING_FLASH   -DUSE_JPEG -DUSE_GD_PNG -DUSE_GD_JPEG -DUSE_GD_WBMP 
-DUSE_GD_FT   -DUSE_OGR -DUSE_GDAL -DUSE_POSTGIS 
-I/home/harald/programme/gd/include 
-I/home/harald/programs/proj_test/include 
-I/home/harald/programs/gdal_test/include 
-I/home/harald/programs/gdal_test/include 
-I/home/harald/programme/postgres/include 
-I/home/harald/programme/curl/include 
-I/home/harald/programme/ming/include  shp2img.o  -L. -lmap 
-L/home/harald/programme/gd/lib -lgd -ljpeg -lfreetype -lpng -lz 
-ljpeg -lfreetype -lpng -lz   -L/home/harald/programs/proj_test/lib 
-lproj -ljpeg  -L/home/harald/programs/gdal_test/lib -lgdal 
-L/home/harald/programs/gdal_test/lib -lgdal 
-L/home/harald/programme/ming/lib -lming 
-L/home/harald/programme/postgres/lib -lpq 
-L/home/harald/programme/curl/lib -lcurl -lz -lssl -lcrypto -ldl -lz 
-lm -lstdc++   -o shp2img
./libmap.a(mapresample.o): In function `msResampleGDALToMap':
mapresample.o(.text+0x1091): undefined reference to `msGetGDALGeoTransform'
collect2: ld returned 1 exit status
make: *** [shp2img] Fehler 1

It does not find msGetGDALGeoTransform. Where is this function defined? 
Can you help me here?

Thanks for your help!

Harald


> 
> Sorry for not responding to this message or the form on the osrs-proj list.
> There are definately challenges with using MapServer for global level 
> mapping
> in projections like orthographic, but with some caveats it can work.  I 
> would
> like to dig into this but haven't really had time.
> 
> I did alot of investigation, and have come to the conclusion that the
> PROJ.4 pj_transform() function is not good at returning only a subset of
> points successfully transformed.  Also, the mapresample.c code, the code 
> used
> to reproject *raster* data input to MapServer does not have any reasonable
> handling for output map areas that cannot completely be transformed into 
> the
> input space.
> 
> I have made some corrections to mapresample.c in MapServer and to the
> pj_transform() function in PROJ.4 and was able to produce the attached png
> with the following .map file:
> 
> MAP
> NAME  "w2qpgi5qsg4k"
> DEBUG on
> SIZE  400 400
> EXTENT  -6500000 -6500000 6500000 6500000
> 
> LAYER
> NAME  "world"
> TYPE  Raster
> STATUS  default
> DATA "/home/warmerda/public_html/world.tif"
> PROCESSING "DITHER=YES"
> 
> PROJECTION
> "proj=latlong"
> END
> END
> 
> 
> PROJECTION
> "proj=ortho"
> "lat_0=0.0"
> "lon_0=0.0"
> "x_0=0.0"
> "y_0=0.0"
> END
> 
> END
> 
> Success!
> 
> Now there are still "issues" with global scale mapping in MapServer, but it
> seems we have worked out some significant issues for raster input with
> reprojection at least.
> 
> The updates are all in CVS, but if you are not comfortable with building
> from CVS it may be a while before you can get things working at your end.
> Keep in mind that updates are required to the core MapServer as well as
> the PROJ.4 libraries.
> 
> Best regards,
> 
> 
> ------------------------------------------------------------------------
>