Am Dienstag 13 Oktober 2009 12:30:32 schrieben Sie: Hello,
the mapfile i've tested looks similar to your setup. I'am just trying to strip out everything which is not needed. But still it gets stuck in the EPSG:4326 projection. From strace i can see that libgdal is touched but the search fails for '/usr/lib/gdal16plugins' (if that is important). Regards Heiko Call: /usr/lib/cgi-bin/mapserv -nh 'QUERY_STRING=MAP=/home/schroete/sciadaten/raster_test/raster.map&LAYERS=Bild&SRS=epsg:3575&mode=map' > out.png Mapfile: MAP PROJECTION # Does not make any diff here "init=epsg:3575" END WEB METADATA "wms_title" "Raster Test" "wms_onlineresource" "http://kahlo1/cgi-bin/mapserv?map=/home/schroete/sciadaten/raster_test/raster.map&" "wms_srs" "epsg:3575 epsg:4326" END END LAYER NAME "Bild" DATA "/home/schroete/sciadaten/raster_test/wmap_work.png" TYPE RASTER STATUS ON # PROJECTION # Comment or not won't change a thing # "init=epsg:3575" # END END END > Hi > > The projection of your Bild layer is > PROJECTION > "proj=ortho" > END > > I am not familiar with ortho projection, but are you sure that this is > enough? In your gdalwarp example you are using also "lat_0=90". > There is no difference in reprojecting raster or vector data. This is a > working example of one orthophoto layer: > METADATA > "wms_title" "orthophotos 92-96" > "wms_srs" "EPSG:2393 EPSG:3067 EPSG:4326 EPSG:900913" > END > PROJECTION > "init=epsg:2393" > END > > -Jukka- > > > Heiko Schröter wrote: > > > > > > Am Dienstag 13 Oktober 2009 11:50:13 schrieben Sie: > > > > Hello, > > > > sorry, no that doesn't work. I have tried with epsg:3575 > > (wich is somewhat polar and recognized by GDALWARP) but that > > fails as well. > > Does anyone has a working example for reprojecting raster > > images via gdallib ? > > > > > > Regards > > Heiko > > > > > > > > > Hi, > > > > > > Perhaps it would work by making "+proj=ortho +lat_0=90" as > > the default projection of your map by adding MAP level > > PROJECTION block as: > > > > > > PROJECTION > > > "proj=ortho" > > > "lat_0=90" > > > "no_defs" > > > END > > > > > > If you want to get the image through WMS I suppose you > > should give it your own EPSG number and add corresponding > > line to epsg file, just like it has been done with EPSG:900913. > > > > > > -Jukka Rahkonen- > > > > > > > Heiko Schröter wrote: > > > > > > > > Hi, > > > > > > > > is there a small howto recipe of how to reproject Raster > > Images via > > > > MapServer (CGI call) with the Gdal Lib ? > > > > I need to go from EPSG:4326 to "+proj=ortho +lat_0=90" (Polar > > > > View) and i cannot figure out how to pass the Projection. > > > > The output image stays in EPSG:4326. > > > > This is not a proj4 issue. WMS calls for vector maps to be > > > > reprojected are ok (in a different scenario). > > > > > > > > Any pointer to some infos or sample code are highly appreciated. > > > > > > > > Thanks and Regards > > > > Heiko > > > > -------------------------------------------------------------- > > > > -------------------------------- > > > > > > > > Reprojecting with gdalwarp works fine: > > > > gdalwarp -of GTiff -s_srs "epsg:4326" -t_srs "+proj=ortho > > > > +lat_0=90" -wo SOURCE_EXTRA=50 wmap_work.png wmap_out.png > > > > > > > > /usr/lib/cgi-bin/mapserv -v > > > > MapServer version 5.4.0-rc1 OUTPUT=GIF OUTPUT=PNG OUTPUT=JPEG > > > > OUTPUT=WBMP OUTPUT=SVG SUPPORTS=PROJ SUPPORTS=FREETYPE > > > > SUPPORTS=ICONV SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT > > > > SUPPORTS=THREADS INPUT=EPPL7 INPUT=POSTGIS INPUT=OGR INPUT=GDAL > > > > INPUT=SHAPEFILE > > > > > > > > The Test mapfile: > > > > MAP > > > > NAME WELT > > > > CONFIG "MS_ERRORFILE" "/tmp/ms_raster_error.txt" > > > > DEBUG 5 > > > > STATUS ON > > > > SIZE 800 400 > > > > IMAGETYPE PNG24 > > > > IMAGECOLOR 240 240 240 > > > > SHAPEPATH "data" > > > > EXTENT -180 -90 180 90 > > > > UNITS DD > > > > WEB > > > > IMAGEPATH "images" > > > > IMAGEURL "images" > > > > MINSCALE 50000 > > > > MAXSCALE 400000000 > > > > METADATA > > > > # "wms_title" "Example WMS Server" > > > > # "wms_onlineresource" > > > > "http://kahlo1/cgi-bin/mapserv?map=heiko5.map&" > > > > "wms_srs" "ortho epsg:4326 epsg:900913" > > > > # or 'proj=ortho' doesn't matter > > > > END > > > > END > > > > > > > > LAYER > > > > NAME "Bild" > > > > DATA "/home/schroete/sciadaten/raster_test/wmap_work.png" > > > > TYPE RASTER > > > > STATUS ON > > > > PROJECTION > > > > "proj=ortho" > > > > END > > > > END > > > > END > > > > > > > > The MapServer Test call: > > > > /usr/lib/cgi-bin/mapserv -nh > > > > 'QUERY_STRING=MAP=/home/schroete/sciadaten/raster_test/raster. > > > > map&LAYERS=Bild&TRANSPARENT=true&SERVICE=WMS&VERSION=1.1.1&FOR > > > MAT=image/png&SRS=ortho&mode=map' > out.png > > > > > > > > > > > > Debug output: > > > > [Tue Oct 13 10:29:35 2009].656759 msDrawRasterLayerLow(Bild): > > > > entering. > > > > [Tue Oct 13 10:29:35 2009].845281 msDrawMap(): Layer 0 (Bild), > > > > 0.189s [Tue Oct 13 10:29:35 2009].845376 msDrawMap(): > > > > Drawing Label Cache, 0.000s [Tue Oct 13 10:29:35 2009].845397 > > > > msDrawMap() total time: 0.193s [Tue Oct 13 10:29:35 > > > > 2009].958674 msSaveImage() total time: 0.113s [Tue Oct 13 > > > > 10:29:35 2009].958837 mapserv request processing time > > (loadmap not > > > > incl.): 0.307s [Tue Oct 13 10:29:35 2009].958856 > > > > msFreeMap(): freeing map at 0x1729dc0. > > > > > > _______________________________________________ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users