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

[Heiko Schröter]

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=&TIME=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=&TIME=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
>> ---------------------------------------------------------------------
>> --
>
>