[mapserver-users] MapServer CGI GDAL Raster reprojection

[Rahkonen Jukka]

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.
> 
> 
> 
> 
> --
> --------------------------------------------------------------
> ---------
> !!! NEUE TELEPHONNUMMER !!!
> !!! NEW PHONE NUMBER !!!
> 
> Dipl.-Ing. Heiko Schröter
> Institute of Environmental Physics (IUP)    phone: 
> ++49-(0)421-218-62092
> Institute of Remote Sensing (IFE)           fax:   
> ++49-(0)421-218-4555
> University of Bremen (FB1)
> P.O. Box 330440               email:  
> schroete at iup.physik.uni-bremen.de
> Otto-Hahn-Allee 1           
> 28359 Bremen                
> Germany
> --------------------------------------------------------------
> ---------
> _______________________________________________
> mapserver-users mailing list
> mapserver-users at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/mapserver-users
>