AW: [mapserver-users] visualisation of image with overviews
Hi Toni, I think gdal_translate is not able to save overviews once you produced them with gdaladdo before. Therefore you do have to use gdaladdo after you are done with all the other stuff. But as far as I know you can't do both: Use Internal Tiling and overviews, because gdaladdo does not support the creation option for internal tiling. I don't know what gdal_translate does with the internal overviews once they are created, maybe the one with the lowest resolution is saved. ;-) I did a lot of downsampling aerial images using gdalwarp with -rc. for creating 'external overviews'. This way I did always get very good results and I can use internal tiling too. For the mapfile I calculate proper scales depending on input image resolution and the resolution of the presentation (e.g. max scale = 20cm/dot input image x 28 dpc (72dpi) /or 38 dpc (96dpi). For sure in any case of resampling arial images you do need 24bit images. This might cause more effort, but in my opinion the results are good and fast. Best regards Andreas Mit freundlichen Grüßen Im Auftrag Andreas Fischer Kreis Unna - Der Landrat Zentrale Datenverarbeitung Friedrich-Ebert-Straße 17 59425 Unna Fon 02 3 03 27-44 16 Fax 0 23 03 27-28 96 andreas.fisc...@kreis-unna.de www.kreis-unna.de -Ursprüngliche Nachricht- Von: mapserver-users-boun...@lists.osgeo.org [mailto:mapserver-users-boun...@lists.osgeo.org] Im Auftrag von Toni Pignataro Gesendet: Freitag, 11. Februar 2011 15:51 An: mapserver-users@lists.osgeo.org Betreff: [mapserver-users] visualisation of image with overviews Hi list, I have a nonsatisfying visualisation of an geotiff image with overviews in mapserver. The different overviews look really poor zooming in the map. What i tried until now is, resampling the original image with gdal_addo (gdaladdo -r nearest , gdaladdo -r average_mp, gdaladdo -r average_magphase, gdaladdo -r mode) and i also made an gdal_translate (gdal_translate -of GTiff -co TILLED=YES), but the results did not look even better. Any one an suggestion how to handle the orginalimage? Thanks, Toni ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users Diese E-Mail wurde beim Ausgang auf Viren geprueft. Wegen der potentiellen Gefahr auf den Uebertragungswegen wird zu einer Vireneingangskontrolle geraten. Eine Haftung für Virenfreiheit wird ausgeschlossen. ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
AW: [mapserver-users] Spatial Filter with WFS
Hi Jukka, I tried another getFeature request using your syntax. I only replaced the 'EPSG:4326' to 'EPSG:25832' since the coordinates are in this reference system. Unfortunately I got the same result. It doesn't matter if I ask for VERSION=1.1.0 or VERSION=1.0.0 - I tried out both. Request http://192.168.2.117/cgi-bin/mapserv?map=/mapserver/mapfiles/u218768/ku_testwfs.mapREQUEST=getFeaturetypename=folie001SERVICE=WFSVERSION=1.0.0; Filter= ogc:Filter ogc:Intersects xmlns:gml='http://www.opengis.net/gml' ogc:PropertyName xmlns:topp=http://www.openplans.org/topp;topp:the_geom/ogc:PropertyName gml:Polygon srsName=EPSG:25832 gml:outerBoundaryIs gml:LinearRing gml:coordinates cs=, decimal=. ts= 409350.0,5709890.0 409350.0,5709910.0 409370.0,5709910.0 409370.0,5709890.0 409350.0,5709890.0/gml:coordinates /gml:LinearRing /gml:outerBoundaryIs /gml:Polygon /ogc:Intersects /ogc:Filter Response ?xml version=1.0 encoding=ISO-8859-1 ? - wfs:FeatureCollection xmlns:ms=http://mapserver.gis.umn.edu/mapserver; xmlns:gml=http://www.opengis.net/gml; xmlns:wfs=http://www.opengis.net/wfs; xmlns:ogc=http://www.opengis.net/ogc; xmlns:xsi=http://www.w3.org/2001/XMLSchema-instance; xsi:schemaLocation=http://mapserver.gis.umn.edu/mapserver http://192.168.2.117/cgi-bin/mapserv?map=/mapserver/mapfiles/u218768/ku_testwfs.mapSERVICE=WFSVERSION=1.1.0REQUEST=DescribeFeatureTypeTYPENAME=folie001OUTPUTFORMAT=text/xml; subtype=gml/3.1.1 http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd; numberOfFeatures=4 - gml:boundedBy - gml:Envelope srsName=EPSG:25832 gml:lowerCorner-1.00 -1.00/gml:lowerCorner gml:upperCorner-1.00 -1.00/gml:upperCorner /gml:Envelope /gml:boundedBy /wfs:FeatureCollection The only difference I could detect is the coordinates within the Envelope. With PostGIS datasource it is: gml:lowerCorner-1.00 -1.00/gml:lowerCorner gml:upperCorner-1.00 -1.00/gml:upperCorner With MapInfo it is: gml:lowerCorner-1.00 -1.00/gml:lowerCorner gml:upperCorner409388.231116 5709944.987211/gml:upperCorner Even if the upperCorner is set, it is not the upper Corner of the featureCollection (4 features), that should be selected by the filter condition. And in both cases I did not change the EXTENT of the MAP und LAYER sections (EXTENT 389965 5693582 420177 5732815, wfs_extent 389965 5693582 420177 5732815). So please let me know, if you have another idea to solve this problem Tanks again! Best regards Andreas Mit freundlichen Grüßen Im Auftrag Andreas Fischer Kreis Unna - Der Landrat Zentrale Datenverarbeitung Friedrich-Ebert-Straße 17 59425 Unna Fon 02 3 03 27-44 16 Fax 0 23 03 27-28 96 andreas.fisc...@kreis-unna.de www.kreis-unna.de Von: Rahkonen Jukka [mailto:jukka.rahko...@mmmtike.fi] Gesendet: Dienstag, 23. November 2010 15:01 An: Yewondwossen Assefa; Fischer, Andreas Cc: mapserver-users@lists.osgeo.org Betreff: Re: [mapserver-users] Spatial Filter with WFS Hi, Andreas, you may be without a wsf:Query element. Here is a WFS 1.0.0 request that works with Geoserver. ?xml version=1.0 encoding=ISO-8859-1? wfs:GetFeature xmlns:ogc=http://www.opengis.net/ogc; xmlns:gml=http://www.opengis.net/gml; xmlns:wfs=http://www.opengis.net/wfs; service=WFS version=1.0.0 maxFeatures=1000 outputFormat=GML2 wfs:Query xmlns:topp=http://www.openplans.org/topp; typeName=topp:states ogc:Filter ogc:Intersects xmlns:gml='http://www.opengis.net/gml' ogc:PropertyName xmlns:topp=http://www.openplans.org/topp; topp:the_geom/ogc:PropertyName gml:Polygon srsName=EPSG:4326 gml:outerBoundaryIs gml:LinearRing gml:coordinates cs=, decimal=. ts= 409350.0,5709890.0 409350.0,5709910.0 409370.0,5709910.0 409370.0,5709890.0 409350.0,5709890.0/gml:coordinates /gml:LinearRing /gml:outerBoundaryIs /gml:Polygon /ogc:Intersects /ogc:Filter /wfs:Query /wfs:GetFeature Y.A, a question for you: I have thought that also the filters should be presented as gml3 with WFS 1.1.0 and thus posList should be used instead of coordinates. Am I wrong with this belief? By the way, this is how Geoserver 2.1 beta 1 writes out on multipolygon as gml3 topp:the_geom gml:MultiSurface srsDimension=2 srsName=urn:x-ogc:def:crs:EPSG:4326 gml:surfaceMember gml:Polygon gml:exterior gml:LinearRing gml:posList 37.510991 -88.071564 37.4762736 -88.087883 37.442852 -88.311707 37.409308 -88.359177 -Jukka Rahkonen- Lähettäjä: mapserver-users-boun...@lists.osgeo.org [mailto:mapserver-users-boun...@lists.osgeo.org] Puolesta Yewondwossen Assefa Lähetetty: 23. marraskuuta 2010 15:43 Vastaanottaja: Rahkonen Jukka; Fischer, Andreas Kopio: mapserver-users@lists.osgeo.org Aihe: Re: [mapserver-users] Spatial Filter with WFS Hi all, As noted MapServer does not support gml:Surface. The original request with gml:Polygon seems to be correct
[mapserver-users] Filter Encoding with Spatial Filter within WFS-Request
Dear MapServer-User, I configured MapServer version 5.6.5 as WMS- and WFS-Server to offer information about land and properties. I use PostgreSQL/PostGIS (8.4.4/1.5.1) as datastore in this case. The spatial reference system is EPSG:25832 und the type of the geometry is MULTIPOLYGON. OS is Suse Linux Enterprise Server 11, Servicepack1. WMS and WFS work fine until I use a spatial filter within the getFeature-Request of the WFS. In this case the response doesn't deliver any features. Nevertheless the number of features (e.g. numberOf Features='4') is correct. I tried this out with a couple of different boundingboxes and the known number of features that intersect these boundigboxes. In all cases I do not get any features back, but always the right number. So it seems that at least the intersection of the boundingbox and my data works correctly. Please take a look at the following request and response that illustrate that problem: Request: http://myserver/cgi-bin/mapserv?map=/mapserver/mapfiles/u218768/ku_testwfs.mapREQUEST=getFeaturetypename=folie001SERVICE=WFSVERSION=1.1.0; Filter=ogc:Filter Intersects ogc:PropertyNamemsGeometry/ogc:PropertyName gml:Polygon srsName='EPSG:25832' gml:outerBoundaryIs gml:LinearRing gml:coordinates409350,5709890 409350,5709910 409370,5709910 409370,5709890 409350,5709890/gml:coordinates /gml:LinearRing /gml:outerBoundaryIs /gml:Polygon /Intersects /ogc:Filter Response ?xml version=1.0 encoding=ISO-8859-1 ? wfs:FeatureCollection xmlns:ms=http://mapserver.gis.umn.edu/mapserver; xmlns:gml=http://www.opengis.net/gml; xmlns:wfs=http://www.opengis.net/wfs xmlns:ogc=http://www.opengis.net/ogc; xmlns:xsi=http://www.w3.org/2001/XMLSchema-instance; xsi:schemaLocation=http://mapserver.gis.umn.edu/mapserver http://myserver/cgi-bin/mapserv?map=/mapserver/mapfiles/u218768/ku_testwfs.mapSERVICE=WFSVERSION=1.1.0REQUEST=DescribeFeatureTypeTYPENAME=folie001OUTPUTFORMAT=text/xml; subtype=gml/3.1.1 http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd; numberOfFeatures=4 gml:boundedBy gml:Envelope srsName=EPSG:25832 gml:lowerCorner-1.00 -1.00/gml:lowerCorner gml:upperCorner-1.00 -1.00/gml:upperCorner /gml:Envelope /gml:boundedBy /wfs:FeatureCollection Indeed of this behavior the WFS returns features when I do not use a spatial filter. Therefore in my opinion the WFS should be configured right in general. To me the spatial filter seems to be the problem. Unfortunately I do not have any ideas what is going wrong. For more information you can take a look at my MapFile and a part of the capabilities document returned from the WFS below: MAP NAME ku_testwfs STATUS ON SIZE 400 520 EXTENT 389965 5693582 420177 5732815 UNITS METERS PROJECTION init=epsg:25832 END SYMBOLSET /mapserver/symbols/symbols_alk.sym FONTSET /mapserver/fonts/fonts.fnt IMAGECOLOR 255 255 255 WEB IMAGEPATH /mapbender/mb262/http/tmp/ IMAGEURL http://myserver/mapbender/tmp/; METADATA ows_onlineresource http://myserver/cgi-bin/mapserv?map=/mapserver/mapfiles/u218768/ku_testwfs.map; ows_title Liegenschaftskarte, Kreis Unna ows_accessconstraints Nutzung nur für den Dienstgebrauch. ows_abstract Flurstücksfolie der Liegenschaftskarte. ows_contactperson Andreas Fischer ows_contactorganization Kreis Unna ows_contactposition Zentrale Datenverarbeitung, Geodatenmanagement wms_srs EPSG:25832 EPSG:31462 EPSG:31463 EPSG:31466 EPSG:31467 EPSG:4326 wms_feature_info_mime_type text/html wfs_srs EPSG:25832 END #METADATA END #WEB LAYER NAME folie001 STATUS ON TYPE POLYGON DUMP TRUE CONNECTIONTYPE postgis CONNECTION user=dbuser password=dbpassword dbname=dbname host=dbhost port=dbport DATA the_geom from ku_view_folie001 USING UNIQUE objnr USING SRID=25832 METADATA ows_title Flurstücke wms_srs EPSG:25832 EPSG:31462 EPSG:31463 EPSG:31466 EPSG:31467 EPSG:4326 wms_feature_info_mime_type text/html wfs_srs EPSG:25832