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

2013-12-19 Thread Eichner, Andreas - SID-NLKM
Have you set the PROCESSING directive at the grid layer? With the patch 
applied, adding PROCESSING FORCE_FULL_EXTENT=YES at the GRID layer fixed 
it for me, ie:

 LAYER
 TYPE LINE
 STATUS DEFAULT
 EXTENT -180 -90 180 90
 NAME GRID
 OPACITY 50

 # the following line enables the hack:
PROCESSING FORCE_FULL_EXTENT=YES

I used the BBOXes from your examples below to generate the attached images.

 The POSTGIS (i.e. POLYGON) and RASTER layers are fine in any zoom or
 pan
 level. The LINE layer is ok if the BBOX has not been paned to far of
 center (for the ortho projection) and if extremly zoomed in the
 poles
 are cut.
 
 It looks like that POLYGON and LINE Layers are treated a little bit
 differently.
 
 If you like here is the test site with the patch included.
 
 Zoomed in, plus poles cut out in LINE layer:
 http://www.iup.uni-bremen.de/warehouse/cgi-bin/laura?BBOX=-250,-
 250,250,250FORMAT=IMAGE/PNGFROMDAT=2006-08-
 06%2000:00HOEHE=-
 1LAYERS=KARTE1,WILLI,GRIDMAP=/var/www/localhost/mapserver/laura.map
 PRODUKT=stro3_21REQUEST=GETMAPSERVICE=WMSSRS=EPSG:0816STYLES=TI
 ME=2006-10-06%2000:00/2006-10-08%2023:59TODAT=2006-08-
 07%2023:59TRANSPARENT=FALSEVERSION=1.1.1WIDTH=600HEIGHT=600
 
 Half a LINE layer missing when BBOX to far off center:
 http://www.iup.uni-bremen.de/warehouse/cgi-bin/laura?BBOX=-650,-
 650,250,250FORMAT=IMAGE/PNGFROMDAT=2006-08-
 06%2000:00HOEHE=-
 1LAYERS=KARTE1,WILLI,GRIDMAP=/var/www/localhost/mapserver/laura.map
 PRODUKT=stro3_21REQUEST=GETMAPSERVICE=WMSSRS=EPSG:0816STYLES=TI
 ME=2006-10-06%2000:00/2006-10-08%2023:59TODAT=2006-08-
 07%2023:59TRANSPARENT=FALSEVERSION=1.1.1WIDTH=600HEIGHT=600
 
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


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

2013-12-19 Thread Heiko Schröter

Oh Andreas, blind me.
Copied it into the POSTGIS layer but not the LINE Layer.
Sorry for that.
It works now as described, panning and zooming is ok.

Btw, where to send the bugreport ?

Thanks and sorry for the additional noise.

Regards
Heiko


Am 19.12.2013 10:01, schrieb Eichner, Andreas - SID-NLKM:

Narf. Send before attaching the images. So here it goes again, sorry.

Have you set the PROCESSING directive at the grid layer? With the
patch applied, adding PROCESSING FORCE_FULL_EXTENT=YES at the GRID layer
fixed it for me, ie:
  
   LAYER

   TYPE LINE
   STATUS DEFAULT
   EXTENT -180 -90 180 90
   NAME GRID
   OPACITY 50
  
   # the following line enables the hack:

   PROCESSING FORCE_FULL_EXTENT=YES
  
I used the BBOXes from your examples below to generate the attached

images.


The POSTGIS (i.e. POLYGON) and RASTER layers are fine in any zoom or
pan
level. The LINE layer is ok if the BBOX has not been paned to far of
center (for the ortho projection) and if extremly zoomed in the
poles
are cut.

It looks like that POLYGON and LINE Layers are treated a little bit
differently.

If you like here is the test site with the patch included.

Zoomed in, plus poles cut out in LINE layer:
http://www.iup.uni-bremen.de/warehouse/cgi-bin/laura?BBOX=-250,-
250,250,250FORMAT=IMAGE/PNGFROMDAT=2006-08-
06%2000:00HOEHE=-
1LAYERS=KARTE1,WILLI,GRIDMAP=/var/www/localhost/mapserver/laura.map
PRODUKT=stro3_21REQUEST=GETMAPSERVICE=WMSSRS=EPSG:0816STYLES=TI
ME=2006-10-06%2000:00/2006-10-08%2023:59TODAT=2006-08-
07%2023:59TRANSPARENT=FALSEVERSION=1.1.1WIDTH=600HEIGHT=600

Half a LINE layer missing when BBOX to far off center:
http://www.iup.uni-bremen.de/warehouse/cgi-bin/laura?BBOX=-650,-
650,250,250FORMAT=IMAGE/PNGFROMDAT=2006-08-
06%2000:00HOEHE=-
1LAYERS=KARTE1,WILLI,GRIDMAP=/var/www/localhost/mapserver/laura.map
PRODUKT=stro3_21REQUEST=GETMAPSERVICE=WMSSRS=EPSG:0816STYLES=TI
ME=2006-10-06%2000:00/2006-10-08%2023:59TODAT=2006-08-
07%2023:59TRANSPARENT=FALSEVERSION=1.1.1WIDTH=600HEIGHT=600



--
---
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:  schro...@iup.physik.uni-bremen.de
Otto-Hahn-Allee 1
28359 Bremen
Germany
---

___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


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

2013-12-19 Thread Eichner, Andreas - SID-NLKM

 Oh Andreas, blind me.
 Copied it into the POSTGIS layer but not the LINE Layer.
 Sorry for that.
 It works now as described, panning and zooming is ok.
 
 Btw, where to send the bugreport ?


No problem. I just found https://github.com/mapserver/mapserver/issues/4483 
which is tagged with MapScript. But reading through the bug description it 
sound like it is strongly related to your case. Since your problem seems to 
start in msDrawVectorLayer() this turns into a more general problem.
I'm currently wondering why msProjectRect() transforms your request BBOX 
(-650,-650,250,250) into 
(-179.464412, 0.002066, 178.470560, 66.995391) and returns MS_SUCCESS.

It seems to be a good idea to create a new issue at 
https://github.com/mapserver/mapserver/issues, tag it with Bug, Proj 
Support and MapServer C Library, reference the issue mentioned above and 
attach the map file, proj string and request that triggers the issue. You 
should also mention the found cause: msProjectRect() used by 
msDrawVectorLayer() doesn't create a proper query BBOX to fetch the features in 
such cases.
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users


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

2013-12-19 Thread Eichner, Andreas - SID-NLKM
Looks like it has been found already 2 years ago, so you might hook onto this

https://github.com/mapserver/mapserver/issues/4035


 -Ursprüngliche Nachricht-
 Von: mapserver-users-boun...@lists.osgeo.org [mailto:mapserver-users-
 boun...@lists.osgeo.org] Im Auftrag von Eichner, Andreas - SID-NLKM
 Gesendet: Donnerstag, 19. Dezember 2013 14:16
 An: 'Heiko Schröter'; mapserver-users@lists.osgeo.org
 Betreff: Re: [mapserver-users] Rewarping/Reprojection with GDAL and -
 wo EXTRA_SOURCE parameter
 
 
  Oh Andreas, blind me.
  Copied it into the POSTGIS layer but not the LINE Layer.
  Sorry for that.
  It works now as described, panning and zooming is ok.
 
  Btw, where to send the bugreport ?
 
 
 No problem. I just found
 https://github.com/mapserver/mapserver/issues/4483 which is tagged
 with MapScript. But reading through the bug description it sound like
 it is strongly related to your case. Since your problem seems to
 start in msDrawVectorLayer() this turns into a more general problem.
 I'm currently wondering why msProjectRect() transforms your request
 BBOX (-650,-650,250,250) into
 (-179.464412, 0.002066, 178.470560, 66.995391) and returns
 MS_SUCCESS.
 
 It seems to be a good idea to create a new issue at
 https://github.com/mapserver/mapserver/issues, tag it with Bug,
 Proj Support and MapServer C Library, reference the issue
 mentioned above and attach the map file, proj string and request that
 triggers the issue. You should also mention the found cause:
 msProjectRect() used by msDrawVectorLayer() doesn't create a proper
 query BBOX to fetch the features in such cases.
 ___
 mapserver-users mailing list
 mapserver-users@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/mapserver-users
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/mapserver-users