[mapserver-users] Rewarping/Reprojection with GDAL and -wo EXTRA_SOURCE parameter

[Eichner, Andreas - SID-NLKM]

Can you provide a simple test case, probably a MAP file with inline data? I've created some simple test data but I am unable to reproduce the problem. 

> -----Ursprüngliche Nachricht-----
> Von: Heiko Schröter [mailto:schroete at iup.physik.uni-bremen.de]
> Gesendet: Dienstag, 17. Dezember 2013 21:01
> An: Eichner, Andreas - SID-NLKM; mapserver-users at lists.osgeo.org;
> Heiko Schröter
> Betreff: Re: [mapserver-users] Rewarping/Reprojection with GDAL and -
> wo EXTRA_SOURCE parameter
> 
> Seems not to work. Failure can be reproduced by changing BBOX params.
> I've setup a test site with three layers so you may have a look.
> Blue Marble = RASTER
> Grid = LINE
> Coloured Rect = POSTGIS
> 
> 
> Warping ok from epsg4326->ortho:
> 
> http://www.iup.uni-bremen.de/warehouse/cgi-bin/laura?BBOX=-6500000,-
> 6500000,6500000,6500000&FORMAT=IMAGE/PNG&FROMDAT=2006-08-
> 06%2000:00&HOEHE=-
> 1&LAYERS=KARTE1,WILLI,GRID&MAP=/var/www/localhost/mapserver/laura.map
> &PRODUKT=stro3_21&REQUEST=GETMAP&SERVICE=WMS&SRS=EPSG:0815&STYLES=&TI
> ME=2006-10-06%2000:00/2006-10-08%2023:59&TODAT=2006-08-
> 07%2023:59&TRANSPARENT=FALSE&VERSION=1.1.1&WIDTH=600&HEIGHT=600
> 
> Warping distorted:
> 
> http://www.iup.uni-bremen.de/warehouse/cgi-bin/laura?BBOX=-4500000,-
> 4500000,4500000,4500000&FORMAT=IMAGE/PNG&FROMDAT=2006-08-
> 06%2000:00&HOEHE=-
> 1&LAYERS=KARTE1,WILLI,GRID&MAP=/var/www/localhost/mapserver/laura.map
> &PRODUKT=stro3_21&REQUEST=GETMAP&SERVICE=WMS&SRS=EPSG:0815&STYLES=&TI
> ME=2006-10-06%2000:00/2006-10-08%2023:59&TODAT=2006-08-
> 07%2023:59&TRANSPARENT=FALSE&VERSION=1.1.1&WIDTH=600&HEIGHT=600
> 
> 
> 
> Am 17.12.2013 17:57, schrieb Eichner, Andreas - SID-NLKM:
> > Can you try if PROCESSING "POLYLINE_NO_CLIP=YES" helps on TYPE LINE
> layers?
> >
> >> -----Ursprüngliche Nachricht-----
> >> Von: Heiko Schröter [mailto:schroete at iup.physik.uni-bremen.de]
> >> Gesendet: Dienstag, 17. Dezember 2013 16:46
> >> An: Eichner, Andreas - SID-NLKM; mapserver-users at lists.osgeo.org
> >> Betreff: Re: [mapserver-users] Rewarping/Reprojection with GDAL
> and -
> >> wo EXTRA_SOURCE parameter
> >>
> >> Thanks for the hint, but sorry no, does not change rewarping
> >> behaviour
> >> of POLYGON or LINE Layers.
> >> Distortion and missing parts do remain.
> >>
> >> The odd thing:
> >> The projection of all layers is ok if the BBOX equals or is
> greater
> >> than
> >> max extension of the map.
> >> i.e. BBOX=-6500000,-6500000,6500000,6500000 [proj=ortho units=m]
> >> By "zooming" in i.e. reducing the BBOX the RASTER layer rewarps
> ok,
> >> POLYGON and LINE layers do not.
> >>
> >>
> >> Am 17.12.2013 16:11, schrieb Eichner, Andreas - SID-NLKM:
> >>> Looks like that causes the same effect as setting
> >>>     PROCESSING "LOAD_WHOLE_IMAGE=YES"
> >>>
> >>>> -----Ursprüngliche Nachricht-----
> >>>> Von: mapserver-users-bounces at lists.osgeo.org [mailto:mapserver-
> >> users-
> >>>> bounces at lists.osgeo.org] Im Auftrag von Heiko Schröter
> >>>> Gesendet: Dienstag, 17. Dezember 2013 14:09
> >>>> An: mapserver-users at lists.osgeo.org
> >>>> Betreff: [mapserver-users] Rewarping/Reprojection with GDAL and
> -
> >> wo
> >>>> EXTRA_SOURCE parameter
> >>>>
> >>>> To make use of it in mapserver a hint has been given in an old
> >>>> mailing
> >>>> thread (sorry, lost the source) to adjust the following code
> >> snippet
> >>>> in
> >>>> resample.c:
> >>>>
> >>>> Version: mapserver-6.4.0
> >>>> resample.c
> >>>> <snip>
> >>>> /* -------------------------------------------------------------
> --
> >> ---
> >>>> -- */
> >>>>      /*      Project desired extents out by 2 pixels, and then
> >> strip
> >>>> to      */
> >>>>      /*      available
> >>>> data.                                                 */
> >>>>      /*
> >>>> ----------------------------------------------------------------
> --
> >> --
> >>>> */
> >>>>      memcpy( &sOrigSrcExtent, &sSrcExtent, sizeof(sSrcExtent) );
> >>>>      sSrcExtent.minx = floor(sSrcExtent.minx-1.0);
> >>>>      sSrcExtent.maxx = ceil (sSrcExtent.maxx+1.0);
> >>>>      sSrcExtent.miny = floor(sSrcExtent.miny-1.0);
> >>>>      sSrcExtent.maxy = ceil (sSrcExtent.maxy+1.0);
> >>>> <snap>
> >>>>
> >>>> Changing the -/+1.0 to -/+1250.0, for example, warps the RASTER
> >> image
> >>>> correctly.
> >>>>
> >>>
> >>
> >>
> >> --
> >> ------------------------------------------------------------------
> ---
> >> --
> >> 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-
> >> 62070
> >> University of Bremen (FB1)
> >> P.O. Box 330440               email:  schroete at iup.physik.uni-
> >> bremen.de
> >> Otto-Hahn-Allee 1
> >> 28359 Bremen
> >> Germany
> >> ------------------------------------------------------------------
> ---
> >> --
> >
> >