Re: [mapserver-users] zoom extent
Template example (but other demo on the web are still very useful for basic commands): I have a div which must contain something created by Mapserver, so I use a particular CGI variable and Mapserver constructs that element (the legend, the reference map, ...) Use square brackets [ ] for defining the cgi-variables. div [cgi-variable] /div or if you want to view the actual scale input name=scale value=[scale] / Mapserver itself puts the real actual scale value in the value attribute of the input. Hope this helps Bye, Carlo john tse ha scritto: thanks for the reply, however, i've been playing around with it the last few days and still no luck for a customized zoom button. It would be most useful if anyone can actually show me a sample script of how to put something similar together so that i can learn from it. I don't get how the template fits with the CGI. On Sun, Apr 13, 2008 at 11:20 PM, Carlo Pelliconi [EMAIL PROTECTED] mailto:[EMAIL PROTECTED] wrote: I think you should try strarting with mapserver basic, not javascript. Mapserver CGI documentation (for functions) here: http://mapserver.gis.umn.edu/docs/reference/cgi/controls Then, for the template (i.e. the html which embeds the map): http://mapserver.gis.umn.edu/docs/reference/templatereference/templatereference/view I've never used javascript for functions such as zoom in, because mapserver itself has them already defined. Good work. Carlo john tse ha scritto: Hi,Guys I'm currently doing a volunteer project and i need to get a zoom to button to work. This is the code i have right now html body head script type = text/javascript // zoom to Region function GoRegion() { //alert(document.mapserv.mapzoomarea.value); document.mapserv.target=; document.mapserv.mode.value=browse; document.mapserv.imgext.value = (document.mapserv.mapzoomarea.value); document.mapserv.submit(); /script /head select name=mapzoomarea onChange=GoRegion(); optionZoom to Region/option option/option option value=5457203 487104 5451873 482304Reserve 2/option option value=5435308 491431 5434498 490403Reserve 4/option option value=5425919 541,907 5502970 444168Musqueam Statement Of Intent/option option value=5435189 475830 5457751 504449Richmond Boundary/option option value=5448859 481843 5463269 500122Vancouver Boundary /option option value=5432692 482228 5457268 515426North Arm Fraser River /option /select /body /html it doesn't do anything, and i have no idea what's wrong with it, any help will be much appreciated. thanx ___ mapserver-users mailing list mapserver-users@lists.osgeo.org mailto:mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users -- ing. Carlo Pelliconi [EMAIL PROTECTED] mailto:[EMAIL PROTECTED] SIS.TER S.p.A. www.sis-ter.it http://www.sis-ter.it via Mentana, 10 40026 Imola (BO) Italy ATTENZIONE! Questo documento è indirizzato esclusivamente al destinatario. Tutte le informazioni ivi contenute, compresi eventuali allegati, sono soggette a riservatezza secondo i termini del D.Lgs. 196/2003 in materia di privacy e ne è proibita l'utilizzazione da parte di altri soggetti. Si ricorda che la diffusione l'utilizzo e/o la conservazione dei dati ricevuti per errore costituisce violazione di carattere penale ai sensi dell'art. 616 Codice penale (Violazione, sottrazione e soppressione di corrispondenza). Se avesse ricevuto per errore questo messaggio, La preghiamo cortesemente di contattare il mittente al più presto e di cancellare il messaggio subito dopo. Grazie. ATTENTION! This message is being sent from Sis.Ter S.p.A. and may contain information which is confidential or privileged. If you are not the intended recipient, please advise the sender immediately by reply e-mail and delete this message and any attachments without retaining a copy. Any unauthorized use of the content of this message is a breach of your duty to respect the confidentiality of the correspondence between other persons and can expose the responsible party to civil and/or criminal penalties (D.Lgs. 196/2003 Tutela della privacy e art. 616 Codice penale Violazione, sottrazione e soppressione di corrispondenza), and may constitute a more serious offence. Thanks.
[mapserver-users] Connexion between Mapserver MySql
Hi! I would like use Database data to print point on my map. I searched on the web and found the method to that. So, I have the following code in a .ovf file: - MYSQL:gpvfr,user=gpvfr2,password=,host=localhost,port=3306,tables=terrain SELECT a.oaci, lat, 'long', nom FROM terrain t,airfields a WHERE a.oaci = t.oaci EPSG:27582 wkbPoint - And when I run ogrinfo on this file, I obtain: --- ERROR 1: MySQL connect failed for: gpvfr,user=gpvfr2,password=31p87qr2,host=localhost,port=3306,tables=terrain;airfields Access denied for user 'gpvfr2'@'localhost' (using password: YES) FAILURE: Unable to open datasource `var/www/html/ms/map_files/mysql.ovf' with the following drivers. - ESRI Shapefile - MapInfo File - UK .NTF - SDTS - TIGER - S57 - DGN - VRT - REC - Memory - BNA - CSV - GML - GPX - KML - GeoJSON - GMT - SQLite - MySQL - AVCBin -- I'm sure of the connexion informations and MySql is installed on my system and GDAL library has been built with ogr option. Where is the problem? Thanks for your help and sorry for my English, I'm frenchy :) Sincerely yours Clement _ Découvrez les profils Messenger de vos amis ! http://home.services.spaces.live.com/___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] Connexion between Mapserver MySql (2)
hi! Sorry, my Mail box eat a part of my .ovf code. I use the code that I give you like that: http://www.forumsig.org/showthread.php?p=144045#post144045 Thanks Clement _ Avec Hotmail, vos e-mails vous suivent partout ! Mettez Hotmail sur votre mobile ! http://www.messengersurvotremobile.com/?d=hotmail___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] Connexion between Mapserver MySql
Clement, look on the mapserver homepage, I've sent a mysql/ovf howto document there last year. It describes step by step what to do. http://mapserver.gis.umn.edu/docs/howto/debian-odbc-mysql There is another one: http://mapserver.gis.umn.edu/docs/howto/ogrmysql Hope that helps! Regards, Till Hi! I would like use Database data to print point on my map. I searched on the web and found the method to that. So, I have the following code in a .ovf file: - MYSQL:gpvfr,user=gpvfr2,password=,host=localhost,port=3306,tables=terrain SELECT a.oaci, lat, 'long', nom FROM terrain t,airfields a WHERE a.oaci = t.oaci EPSG:27582 wkbPoint - And when I run ogrinfo on this file, I obtain: --- ERROR 1: MySQL connect failed for: gpvfr,user=gpvfr2,password=31p87qr2,host=localhost,port=3306,tables=terrain;airfields Access denied for user 'gpvfr2'@'localhost' (using password: YES) FAILURE: Unable to open datasource `var/www/html/ms/map_files/mysql.ovf' with the following drivers. - ESRI Shapefile - MapInfo File - UK .NTF - SDTS - TIGER - S57 - DGN - VRT - REC - Memory - BNA - CSV - GML - GPX - KML - GeoJSON - GMT - SQLite - MySQL - AVCBin -- I'm sure of the connexion informations and MySql is installed on my system and GDAL library has been built with ogr option. Where is the problem? Thanks for your help and sorry for my English, I'm frenchy :) Sincerely yours Clement _ Découvrez les profils Messenger de vos amis ! http://home.services.spaces.live.com/___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users -- -- terrestris GmbH Co. KG Irmintrudisstrasse 17 53111 Bonn Germany Till Adams Geschäftsführung Tel: +49 (0)228 / 962 899-52 Fax: +49 (0)228 / 962 899-57 [EMAIL PROTECTED] http://www.terrestris.de Amtsgericht Bonn, HRA 6835 - Komplementärin: terrestris Verwaltungs GmbH vertreten durch: Hinrich Paulsen, Till Adams ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] How to tile a large TIF image?
Hi there, I would like to use instead of single large tif image smaller tiles. Although I am well aware of mapserver's and gdal's possibilities to create the shapes for it, I first need to split the large tif into 20 or 50 or 100 tiles. Can anyone give me a hint with what kind of software this is possible? There is an ArcGIS script, but it doesn't work on my machine ( http://arcscripts.esri.com/details.asp?dbid=13978 ). Thanks for any hints, Stef Stefan Schwarzer Lean Back and Relax - Enjoy some Nature Photography http://photoblog.la-famille-schwarzer.de Appetite for Global Data? UNEP GEO Data Portal: http://geodata.grid.unep.ch ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] How to tile a large TIF image?
Thanks for the info. Did it like this: Library/Frameworks/GDAL.framework/Versions/1.5/Programs/gdal_translate -outsize 10'' 10'' -co TILED=YES gebco/bathymetry.tif gebco/ bathymetry_tiled.tif But get the message: Segmentation fault Anything that I did wrong? Or should do differently? Thanks for a hint, Stef Hi, Gdal_translate program has options for you: c:\FWToolsgdal_translate Usage: gdal_translate [--help-general] [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/ CInt16/CInt32/CFloat32/CFloat64}] [-strict] [-of format] [-b band] [-outsize xsize[%] ysize[%]] [-scale [src_min src_max [dst_min dst_max]]] [-srcwin xoff yoff xsize ysize] [-projwin ulx uly lrx lry] [-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value] [-gcp pixel line easting northing [elevation]]* [-mo META-TAG=VALUE]* [-quiet] [-sds] [-co NAME=VALUE]* src_dataset dst_dataset By playing with -srcwin or possibly with -outsize and -projwin you should be able to split your image as you wish. Read more from http://gdal.org/gdal_translate.html -Jukka Rahkonen- Lähettäjä: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] ] Puolesta Stefan Schwarzer Lähetetty: 18. huhtikuuta 2008 12:05 Vastaanottaja: mapserver-users@lists.osgeo.org Aihe: [mapserver-users] How to tile a large TIF image? Hi there, I would like to use instead of single large tif image smaller tiles. Although I am well aware of mapserver's and gdal's possibilities to create the shapes for it, I first need to split the large tif into 20 or 50 or 100 tiles. Can anyone give me a hint with what kind of software this is possible? There is an ArcGIS script, but it doesn't work on my machine ( http://arcscripts.esri.com/details.asp?dbid=13978 ). Thanks for any hints, Stef Stefan Schwarzer Lean Back and Relax - Enjoy some Nature Photography http://photoblog.la-famille-schwarzer.de Appetite for Global Data? UNEP GEO Data Portal: http://geodata.grid.unep.ch ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] How to tile a large TIF image?
Stefan, I don't think you can specify 10'', you need something in pixels. But the command is also probably not what you really need to do. More likely, you should be doing the following: gdal_translate -co TILED=YES gebco/bathymetry.tif gebco/ bathmetry_tiled.tif This will create a single tif file that has an internal block size of 256x256 - you can think of this as having broken your tif up into 256x256 tiles but keeps them all within the same file. This is equivalent to breaking the tif up into individual files, creating a shapefile that has rectangles for each individual files' extent, and using that as a tile index in mapserver. Next, you want to do this: gdaladdo gebco/bathymetry.tif 2 4 6 8 16 this will pre-compute smaller versions of the tif image called overviews at 1/2, 1/4 1/6 1/8 and 1/16 of the original size of the tif - this makes it much more efficient for gdal to return exactly what mapserver is asking for at any given scale. Depending on the range of scales that you need to display your raster image at, you can add more overview levels (or take some away). With this two commands, you can make your rasters much more efficient for mapserver and probably avoid the need to split them up into many files. It does make the file somewhat larger. If you are working with very large files ( 4GB ) then you may run into some problems with tiff and may need to investigate another format or a compressed format like ECW or MrSID (both requiring commercial licenses I believe). There are also some built in compression schemes for tiffs in gdal which you can apply when running gdal_translate (see http://www.gdal.org/frmt_gtiff.html for creation options), for instance: gdal_translate -co TILED=YES -co COMPRESS=JPEG -co JPEG_QUALITY=80 gebco/bathymetry.tif gebco/bathmetry_tiled.tif would compress the tif using JPEG compression (lossy) set to 80% (low compression, small loss). Cheers Paul On 18-Apr-08, at 5:47 AM, Stefan Schwarzer wrote: Thanks for the info. Did it like this: Library/Frameworks/GDAL.framework/Versions/1.5/Programs/ gdal_translate -outsize 10'' 10'' -co TILED=YES gebco/bathymetry.tif gebco/bathymetry_tiled.tif But get the message: Segmentation fault Anything that I did wrong? Or should do differently? Thanks for a hint, Stef Hi, Gdal_translate program has options for you: c:\FWToolsgdal_translate Usage: gdal_translate [--help-general] [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/ CInt16/CInt32/CFloat32/CFloat64}] [-strict] [-of format] [-b band] [-outsize xsize[%] ysize[%]] [-scale [src_min src_max [dst_min dst_max]]] [-srcwin xoff yoff xsize ysize] [-projwin ulx uly lrx lry] [-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value] [-gcp pixel line easting northing [elevation]]* [-mo META-TAG=VALUE]* [-quiet] [-sds] [-co NAME=VALUE]* src_dataset dst_dataset By playing with -srcwin or possibly with -outsize and -projwin you should be able to split your image as you wish. Read more from http://gdal.org/gdal_translate.html -Jukka Rahkonen- Lähettäjä: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] ] Puolesta Stefan Schwarzer Lähetetty: 18. huhtikuuta 2008 12:05 Vastaanottaja: mapserver-users@lists.osgeo.org Aihe: [mapserver-users] How to tile a large TIF image? Hi there, I would like to use instead of single large tif image smaller tiles. Although I am well aware of mapserver's and gdal's possibilities to create the shapes for it, I first need to split the large tif into 20 or 50 or 100 tiles. Can anyone give me a hint with what kind of software this is possible? There is an ArcGIS script, but it doesn't work on my machine ( http://arcscripts.esri.com/details.asp?dbid=13978 ). Thanks for any hints, Stef Stefan Schwarzer Lean Back and Relax - Enjoy some Nature Photography http://photoblog.la-famille-schwarzer.de Appetite for Global Data? UNEP GEO Data Portal: http://geodata.grid.unep.ch ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users __ Paul Spencer Chief Technology Officer DM Solutions Group Inc http://www.dmsolutions.ca/ ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] Displaying points labels
Hi! I'd like displaying points and labels. My points are on MySql database and I use them with OGR. The connexion is OK and I think my mapfile too. But I don't see the points and the labels. My codes can be see on: http://www.forumsig.org/showthread.php?p=144085#post144085 I hope find someone to help me. Have a nice day. Clement PS: Excuse me for my English, I'm a frenchy! _ Créez votre disque dur virtuel Windows Live SkyDrive, 5Go de stockage gratuit ! http://www.windowslive.fr/skydrive/___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
RE: [mapserver-users] MapServer v. 4.10.0 featureid problem
Jacob is bang on (I'm no expert on ArcGIS :) If I invoke a tool like ogrinfo on the shapefile, I see no attribute called FID. I do see ID and DISTRICT attributes. Try setting gml_featureid ID in your mapfile, which should pick up the field and process accordingly. ..Tom -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Delfos, Jacob Sent: 17 April, 2008 11:55 PM To: mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Hi Guys, I'm no expert on WFS, but I do know ArcGIS. The FID field is not really an attribute. It is more like an internal ESRI unique row ID (which gets updated after edits to be sequential). If you open up the attributes in another application (if you use shapefiles, open up the DBF in OpenOffice Calc or Excel), then you will find there is no FID field. Within ArcGIS, however, the FID field is presented as if it were a normal attribute. Hope this helps, Jacob -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Brian Mixon Sent: 18 April 2008 01:34 To: Kralidis,Tom [Burlington]; mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Tom, As soon as the Brenton on this list gets in to work we'll try to put something together for you. I'm working remotely and need his GIS knowledge to get something for you. Brian -Original Message- From: Kralidis,Tom [Burlington] [mailto:[EMAIL PROTECTED] Sent: Thursday, April 17, 2008 10:51 AM To: Brian Mixon; mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Hi Brian, What's your underlying data? If there's anyway you can package a small testcase, I can do a local test. ..Tom -Original Message- From: Brian Mixon [mailto:[EMAIL PROTECTED] Sent: 17 April, 2008 12:24 PM To: Kralidis,Tom [Burlington]; mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Hi Tom, Yes, it does. It's an Arc field that is present in all our data. Thanks, Brian -Original Message- From: Kralidis,Tom [Burlington] [mailto:[EMAIL PROTECTED] Sent: Thursday, April 17, 2008 10:11 AM To: Brian Mixon; mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Does the column/field fid exist in your data? ..Tom From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Brian Mixon Sent: 17 April, 2008 11:51 AM To: mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: [mapserver-users] MapServer v. 4.10.0 featureid problem Hello, I'm a MapServer newbie and am having trouble with a WFS v. 1.0.0 request. I'm trying to have all the features for a particular layer returned (which happens successfully) with the fid (feature id) included as well (does not happen). I've included the gml_featureid 'fid' tag in the layer's metadata section of the mapfile but only receive the following warning the the xml/gml response:!-- WARNING: FeatureId item 'fid' not found in typename 'senate_districts'. -- I've searched and read the various tickets that claim this problem is fixed so I'm wondering what I'm doing wrong? Any information to point me in the right direction would be greatly appreciated. Regards, Brian Mixon ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users This email and any attachments are intended solely for the named recipient and are confidential. The copying or distribution of them or any information they contain by anyone other than the named recipient is prohibited. If you have received this document in error, please notify the sender and delete all copies from your computer system and destroy all hard copies. It is the recipient's responsibility to check this email and any attachments to this email for viruses before use. ___ 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
Re: [MAPSERVER-USERS] Multiple GetFeatureInfo on Raster Data
Matthew Doyle wrote: Hi MapServer users, We are trying to do a GetFeatureInfo request which will query multiple layers. These are both WMS Raster data layers, namely 'rain_jan' and 'rain_jun'. We are constructing the URL as follows: http://www./latest?SERVICE=wmsVERSION=1.1.1REQUEST=getfeatureinfoLAYERS=rain_jan,rain_junQUERY_LAYERS=rain_jun,rain_janx=200y=150INFO_FORMAT=text/plainsrs=EPSG:4283BBOX=120,-50,180,0FORMAT=image/pngwidth=400height=300stylesFEATURE_COUNT=2 Now, this will only ever query the first layer, it will never pick up multiple layers, no matter how many we put in the LAYERS or QUERY_LAYERS fields, or no matter how high our FEATURE_COUNT var is. All we ever get is the first layer, and other 'features' which are close in proximity (lat/lon) to the point selected within the same layer. Matt, I presume this is a flaw in my implementation of raster query. I would appreciate your filing a ticket with a small, easy to reproduce test case for me. ie. a map file with two small raster files as layers, and the mapserv url (like above) to invoke the request. Assign the ticket directly to me (warmerdam). Best regards, -- ---+-- I set the clouds in motion - turn up | Frank Warmerdam, [EMAIL PROTECTED] light and sound - activate the windows | http://pobox.com/~warmerdam and watch the world go round - Rush| President OSGeo, http://osgeo.org ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] Raster problems
Hello! I'm new in Mapserver, and after testing and testing, searching and searching... I have a problem! I'm doing a map file with two raster layers. I defines the MAXSACLe and MINSCALE for this layers but at the end i have a white screen :S What I'm doing wrong? It's possible to create a raster based map with diferent raster levels for teh same area?(I think yes but i couldn't do it :S) Thanks for teh answers ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] Displaying points label
Hi, you are using a layer of type annotation. This type of layer usually does not display the actual point symbol, it only renders the label at that position. Thus you could try to use type point, to really draw the points and the labels... I just realized that you have two TYPE declarations in your layer: point and annotation. That could also be a problem ;-) Cheers, Michael 2008/4/18 Clemiou [EMAIL PROTECTED]: Hi! I'd like displaying points and labels. My points are on MySql database and I use them with OGR. The connexion is OK and I think my mapfile too. But I don't see the points and the labels. My codes can be see on: http://www.forumsig.org/showthread.php?p=144085#post144085 I hope find someone to help me. Have a nice day. Clement PS: Excuse me for my English, I'm a frenchy! _ Créez votre disque dur virtuel Windows Live SkyDrive, 5Go de stockage gratuit ! http://www.windowslive.fr/skydrive/___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users -- --- Michael Schulz [EMAIL PROTECTED] in medias res Gesellschaft für Informationstechnologie mbH In den Weihermatten 66 79108 Freiburg Tel +49 (0)761 556959-5 Fax +49 (0)761 556959-6 http://www.webgis.de / http://www.zopecms.de --- ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] How to tile a large TIF image?
Thanks a lot for these details. That should be added to the existing documents! But unfortunately I still get a Segmentation fault message... Stef Stefan, I don't think you can specify 10'', you need something in pixels. But the command is also probably not what you really need to do. More likely, you should be doing the following: gdal_translate -co TILED=YES gebco/bathymetry.tif gebco/ bathmetry_tiled.tif This will create a single tif file that has an internal block size of 256x256 - you can think of this as having broken your tif up into 256x256 tiles but keeps them all within the same file. This is equivalent to breaking the tif up into individual files, creating a shapefile that has rectangles for each individual files' extent, and using that as a tile index in mapserver. Next, you want to do this: gdaladdo gebco/bathymetry.tif 2 4 6 8 16 this will pre-compute smaller versions of the tif image called overviews at 1/2, 1/4 1/6 1/8 and 1/16 of the original size of the tif - this makes it much more efficient for gdal to return exactly what mapserver is asking for at any given scale. Depending on the range of scales that you need to display your raster image at, you can add more overview levels (or take some away). With this two commands, you can make your rasters much more efficient for mapserver and probably avoid the need to split them up into many files. It does make the file somewhat larger. If you are working with very large files ( 4GB ) then you may run into some problems with tiff and may need to investigate another format or a compressed format like ECW or MrSID (both requiring commercial licenses I believe). There are also some built in compression schemes for tiffs in gdal which you can apply when running gdal_translate (see http://www.gdal.org/frmt_gtiff.html for creation options), for instance: gdal_translate -co TILED=YES -co COMPRESS=JPEG -co JPEG_QUALITY=80 gebco/bathymetry.tif gebco/bathmetry_tiled.tif would compress the tif using JPEG compression (lossy) set to 80% (low compression, small loss). Cheers Paul On 18-Apr-08, at 5:47 AM, Stefan Schwarzer wrote: Thanks for the info. Did it like this: Library/Frameworks/GDAL.framework/Versions/1.5/Programs/ gdal_translate -outsize 10'' 10'' -co TILED=YES gebco/ bathymetry.tif gebco/bathymetry_tiled.tif But get the message: Segmentation fault Anything that I did wrong? Or should do differently? Thanks for a hint, Stef Hi, Gdal_translate program has options for you: c:\FWToolsgdal_translate Usage: gdal_translate [--help-general] [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/ CInt16/CInt32/CFloat32/CFloat64}] [-strict] [-of format] [-b band] [-outsize xsize[%] ysize[%]] [-scale [src_min src_max [dst_min dst_max]]] [-srcwin xoff yoff xsize ysize] [-projwin ulx uly lrx lry] [-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value] [-gcp pixel line easting northing [elevation]]* [-mo META-TAG=VALUE]* [-quiet] [-sds] [-co NAME=VALUE]* src_dataset dst_dataset By playing with -srcwin or possibly with -outsize and -projwin you should be able to split your image as you wish. Read more from http://gdal.org/gdal_translate.html -Jukka Rahkonen- Lähettäjä: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] ] Puolesta Stefan Schwarzer Lähetetty: 18. huhtikuuta 2008 12:05 Vastaanottaja: mapserver-users@lists.osgeo.org Aihe: [mapserver-users] How to tile a large TIF image? Hi there, I would like to use instead of single large tif image smaller tiles. Although I am well aware of mapserver's and gdal's possibilities to create the shapes for it, I first need to split the large tif into 20 or 50 or 100 tiles. Can anyone give me a hint with what kind of software this is possible? There is an ArcGIS script, but it doesn't work on my machine ( http://arcscripts.esri.com/details.asp?dbid=13978 ). Thanks for any hints, Stef Stefan Schwarzer Lean Back and Relax - Enjoy some Nature Photography http://photoblog.la-famille-schwarzer.de Appetite for Global Data? UNEP GEO Data Portal: http://geodata.grid.unep.ch ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users __ Paul Spencer Chief Technology Officer DM Solutions Group Inc http://www.dmsolutions.ca/ ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] Raster problems
David Martinez Morata wrote: Hello! I'm new in Mapserver, and after testing and testing, searching and searching... I have a problem! I'm doing a map file with two raster layers. I defines the MAXSACLe and MINSCALE for this layers but at the end i have a white screen :S What I'm doing wrong? It's possible to create a raster based map with diferent raster levels for teh same area?(I think yes but i couldn't do it :S) Thanks for teh answers David, Yes, this is possible but a white (presumably empty) map can be caused by many things. First, ensure that each layer is working without minscale/maxscale settings. Then make sure you set minscale/maxscale such that one or the other layer is always going to be enabled. It would be helpful for you to post the layer definition for the two layers in question. It is also helpful to have some understanding of what does work, and what issues you have ruled out. Debugging a map file is a logical process, just like debugging a program. Best regards, -- ---+-- I set the clouds in motion - turn up | Frank Warmerdam, [EMAIL PROTECTED] light and sound - activate the windows | http://pobox.com/~warmerdam and watch the world go round - Rush| President OSGeo, http://osgeo.org ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] How to tile a large TIF image?
Sounds like a problem with your gdal install - I'm running a mac also and have no problems. I recommend going back to the kyngchaos site and downloading the latest GDAL and installing it. Cheers Paul On 18-Apr-08, at 9:02 AM, Stefan Schwarzer wrote: Thanks a lot for these details. That should be added to the existing documents! But unfortunately I still get a Segmentation fault message... Stef Stefan, I don't think you can specify 10'', you need something in pixels. But the command is also probably not what you really need to do. More likely, you should be doing the following: gdal_translate -co TILED=YES gebco/bathymetry.tif gebco/ bathmetry_tiled.tif This will create a single tif file that has an internal block size of 256x256 - you can think of this as having broken your tif up into 256x256 tiles but keeps them all within the same file. This is equivalent to breaking the tif up into individual files, creating a shapefile that has rectangles for each individual files' extent, and using that as a tile index in mapserver. Next, you want to do this: gdaladdo gebco/bathymetry.tif 2 4 6 8 16 this will pre-compute smaller versions of the tif image called overviews at 1/2, 1/4 1/6 1/8 and 1/16 of the original size of the tif - this makes it much more efficient for gdal to return exactly what mapserver is asking for at any given scale. Depending on the range of scales that you need to display your raster image at, you can add more overview levels (or take some away). With this two commands, you can make your rasters much more efficient for mapserver and probably avoid the need to split them up into many files. It does make the file somewhat larger. If you are working with very large files ( 4GB ) then you may run into some problems with tiff and may need to investigate another format or a compressed format like ECW or MrSID (both requiring commercial licenses I believe). There are also some built in compression schemes for tiffs in gdal which you can apply when running gdal_translate (see http://www.gdal.org/frmt_gtiff.html for creation options), for instance: gdal_translate -co TILED=YES -co COMPRESS=JPEG -co JPEG_QUALITY=80 gebco/bathymetry.tif gebco/bathmetry_tiled.tif would compress the tif using JPEG compression (lossy) set to 80% (low compression, small loss). Cheers Paul On 18-Apr-08, at 5:47 AM, Stefan Schwarzer wrote: Thanks for the info. Did it like this: Library/Frameworks/GDAL.framework/Versions/1.5/Programs/ gdal_translate -outsize 10'' 10'' -co TILED=YES gebco/ bathymetry.tif gebco/bathymetry_tiled.tif But get the message: Segmentation fault Anything that I did wrong? Or should do differently? Thanks for a hint, Stef Hi, Gdal_translate program has options for you: c:\FWToolsgdal_translate Usage: gdal_translate [--help-general] [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/ CInt16/CInt32/CFloat32/CFloat64}] [-strict] [-of format] [-b band] [-outsize xsize[%] ysize[%]] [-scale [src_min src_max [dst_min dst_max]]] [-srcwin xoff yoff xsize ysize] [-projwin ulx uly lrx lry] [-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value] [-gcp pixel line easting northing [elevation]]* [-mo META-TAG=VALUE]* [-quiet] [-sds] [-co NAME=VALUE]* src_dataset dst_dataset By playing with -srcwin or possibly with -outsize and -projwin you should be able to split your image as you wish. Read more from http://gdal.org/gdal_translate.html -Jukka Rahkonen- Lähettäjä: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] ] Puolesta Stefan Schwarzer Lähetetty: 18. huhtikuuta 2008 12:05 Vastaanottaja: mapserver-users@lists.osgeo.org Aihe: [mapserver-users] How to tile a large TIF image? Hi there, I would like to use instead of single large tif image smaller tiles. Although I am well aware of mapserver's and gdal's possibilities to create the shapes for it, I first need to split the large tif into 20 or 50 or 100 tiles. Can anyone give me a hint with what kind of software this is possible? There is an ArcGIS script, but it doesn't work on my machine ( http://arcscripts.esri.com/details.asp?dbid=13978 ). Thanks for any hints, Stef Stefan Schwarzer Lean Back and Relax - Enjoy some Nature Photography http://photoblog.la-famille-schwarzer.de Appetite for Global Data? UNEP GEO Data Portal: http://geodata.grid.unep.ch ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users __ Paul Spencer Chief Technology Officer DM Solutions Group Inc http://www.dmsolutions.ca/ __ Paul Spencer Chief Technology
[mapserver-users] MapServer logs
Debugging a map file is a logical process, just like debugging a program. Thus spake Warmerdam in nuther thread. Ok then. I have the following situation -- Client OpenLayers calling MapServer via WMS. MapServer getting data out of a db (homegrown SQLiteGIS) via an OGR connection. CONNECTIONTYPE OGR CONNECTION /Users/punkish/Data/sqlitegis/wi.db DATA SELECT wkt_geometry FROM sg_zip I want MapServer to write every call that it makes to the database via OGR. So, I added the following line to my map file WEB .. LOG /Users/punkish/Sites/sqlitegis/tmp/ms.log END tmp directory is writeable. The map is rendered just fine, thank you, however, no logfile. What is my culpa? -- Puneet Kishor ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] Displaying points lab
Hi Clement, please leave the thread also on the list... Ok, that is not the real debug output, merely the log. Did you try it from the command line with shp2img: shp2img -m aero.map -l AERO -all_debug 5 this should print some debug information to stdout. Are you using mapserver 5? Then you can set a special error log for each mapfile, with the directive CONFIG MS_ERRORFILE /path/to/mserror.log at the map level. Cheers, Michael 2008/4/18 Clemiou [EMAIL PROTECTED]: Hi The one thing that I can see in the mapserver log is: Fri Apr 18 15:18:51 2008,18284,127.0.0.1,MS,3,1.241143 46.540217 1.445917 46.744991,-1.00 -1.00,FRANCE ,normal execution ^^ Very funny isn't it? Clement Date: Fri, 18 Apr 2008 15:03:59 +0200 From: [EMAIL PROTECTED] To: [EMAIL PROTECTED]; mapserver-users@lists.osgeo.org Subject: Re: [mapserver-users] Displaying points labe Hi, could you have a look at the debug output? Turn debugging on with DEBUG on on the map and on the layer level, then you should get debug output, either to your configured mapserver log or the webserver (error) log. Another way woudl be to use shp2img. Cheers, Michael 2008/4/18 Clemiou : Hi, I have modifief my layer for just have the points, as you can see: # # Le layer des aérodromes . # # LAYER NAME AERO TYPE POINT STATUS ON PROJECTION init=epsg:27582 END TEMPLATE ttt_query.html CONNECTIONTYPE OGR CONNECTION mysql.ovf #TOLERANCE 5 #TYPE ANNOTATION #LABELITEM nom DATA my_layer CLASS NAME my_class #LABEL # ANGLE auto #TYPE TRUETYPE #ANTIALIAS TRUE #MINFEATURESIZE 100 #COLOR 255 255 255 #OUTLINECOLOR 255 0 0 #SIZE 20 #FONT arial #FORCE TRUE #END # Fin LABEL AERO STYLE SYMBOL 'circle' SIZE 15 COLOR 0 255 0 #OUTLINECOLOR 255 0 0 END # Fin du STYLE AERO END # Fin de CLASS AERO END # Fin du LAYER AERO But it doesn't work. Thakns for your answer Clement Date: Fri, 18 Apr 2008 14:51:44 +0200 From: [EMAIL PROTECTED] To: [EMAIL PROTECTED] Subject: Re: [mapserver-users] Displaying points label CC: mapserver-users@lists.osgeo.org Hi, you are using a layer of type annotation. This type of layer usually does not display the actual point symbol, it only renders the label at that position. Thus you could try to use type point, to really draw the points and the labels... I just realized that you have two TYPE declarations in your layer: point and annotation. That could also be a problem ;-) Cheers, Michael 2008/4/18 Clemiou : Hi! I'd like displaying points and labels. My points are on MySql database and I use them with OGR. The connexion is OK and I think my mapfile too. But I don't see the points and the labels. My codes can be see on: http://www.forumsig.org/showthread.php?p=144085#post144085 I hope find someone to help me. Have a nice day. Clement PS: Excuse me for my English, I'm a frenchy! _ Créez votre disque dur virtuel Windows Live SkyDrive, 5Go de stockage gratuit ! http://www.windowslive.fr/skydrive/___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users -- --- Michael Schulz [EMAIL PROTECTED] in medias res Gesellschaft für Informationstechnologie mbH In den Weihermatten 66 79108 Freiburg Tel +49 (0)761 556959-5 Fax +49 (0)761 556959-6 http://www.webgis.de / http://www.zopecms.de --- _ Créez votre disque dur virtuel Windows Live SkyDrive, 5Go de stockage gratuit ! http://www.windowslive.fr/skydrive/ -- --- Michael Schulz [EMAIL PROTECTED] in
Re: [mapserver-users] MapServer logs
On 4/18/08, Daniel Morissette [EMAIL PROTECTED] wrote: P Kishor wrote: I want MapServer to write every call that it makes to the database via OGR. So, I added the following line to my map file WEB .. LOG /Users/punkish/Sites/sqlitegis/tmp/ms.log END IIRC, LOG only logs some info on the mapserv request status at the end of its execution. I don't use it much and don't know much about it. To get debugging output, with MapServer 5.0+, set: CONFIG MS_ERRORFILE /Users/punkish/Sites/sqlitegis/tmp/ms.log ... and then set DEBUG level (ON, or number between 1 and 5) at the top-level in the mapfile and in each layer for which you want debugging output. Very nice. tail tmp/ms.log shows me [Fri Apr 18 08:46:54 2008].505608 msOGRLayerClose(/Users/punkish/Data/sqlitegis/wi.db). [Fri Apr 18 08:46:54 2008].505626 msOGRFileClose(/Users/punkish/Data/sqlitegis/wi.db,-1). [Fri Apr 18 08:46:54 2008].505647 msConnPoolRelease(sg_zip,/Users/punkish/Data/sqlitegis/wi.db,0x1126de0) [Fri Apr 18 08:46:54 2008].505654 msConnPoolClose(/Users/punkish/Data/sqlitegis/wi.db,0x1126de0) [Fri Apr 18 08:46:54 2008].506741 msDrawMap(): Layer 0 (sg_zip), 0.257s [Fri Apr 18 08:46:54 2008].216458 msOGRFileOpen(/Users/punkish/Data/sqlitegis/wi.db)... [Fri Apr 18 08:46:54 2008].228475 OGROPen(/Users/punkish/Data/sqlitegis/wi.db) [Fri Apr 18 08:46:54 2008].231712 msConnPoolRegister(sg_zip,/Users/punkish/Data/sqlitegis/wi.db,0x1126de0) Unfortunately, it doesn't log the actual SQL calls being made to the db. Perhaps that is an OGR logging issue. So, the next step is to find a way to log that. More details are available in RFC-28: http://mapserver.gis.umn.edu/development/rfc/ms-rfc-28 Daniel -- Daniel Morissette http://www.mapgears.com/ ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] How to tile a large TIF image?
Dejan, I'm not really sure that there would be a noticeable performance difference, and I definitely don't feel sufficiently educated to comment on why one approach would be superior to another, except that I find it much easier to manage single files and to set up a map file to point to them. The process of creating multiple files and setting up a tileindex can be automated, certainly, but I think it is a little more complicated especially for the beginner. The multiple files approach requires mapserver to open the tileindex file, find the polygons that intersect the requested extent, extract the references to the actual disk files, open the disk files using gdal to get the content and composite those into the memory copy of the map. The single file approach does something very similar except that GDAL handles it and there is only one open file handle. I'm not sure that there is a significant overhead to the first approach just on the basis of files needing to be opened. I'm hoping others more qualified will jump in here (Frank, Ed ...) and provide more educated opinions. Cheers Paul On 18-Apr-08, at 9:40 AM, [EMAIL PROTECTED] wrote: Paul, I am just interested in one thing: If you say breaking a single tiff file with TILED=YES is equivalent to breaking into individual files and using tileindex - does it mean it is ALWAYS better to create a single file with internal tiling than create many files and using tileindex? I suppose the answer is NO, but I would very like to know why? For example, if you need to display a region that covers several blocks, then in the first case you open just one file, in the second you open several files. So the first method has less overhead right? What are the opposite situations then? thanks very much regards, dejan [EMAIL PROTECTED] wrote on 18.04.2008 13:11:47: Stefan, I don't think you can specify 10'', you need something in pixels. But the command is also probably not what you really need to do. More likely, you should be doing the following: gdal_translate -co TILED=YES gebco/bathymetry.tif gebco/ bathmetry_tiled.tif This will create a single tif file that has an internal block size of 256x256 - you can think of this as having broken your tif up into 256x256 tiles but keeps them all within the same file. This is equivalent to breaking the tif up into individual files, creating a shapefile that has rectangles for each individual files' extent, and using that as a tile index in mapserver. Next, you want to do this: gdaladdo gebco/bathymetry.tif 2 4 6 8 16 this will pre-compute smaller versions of the tif image called overviews at 1/2, 1/4 1/6 1/8 and 1/16 of the original size of the tif - this makes it much more efficient for gdal to return exactly what mapserver is asking for at any given scale. Depending on the range of scales that you need to display your raster image at, you can add more overview levels (or take some away). With this two commands, you can make your rasters much more efficient for mapserver and probably avoid the need to split them up into many files. It does make the file somewhat larger. If you are working with very large files ( 4GB ) then you may run into some problems with tiff and may need to investigate another format or a compressed format like ECW or MrSID (both requiring commercial licenses I believe). There are also some built in compression schemes for tiffs in gdal which you can apply when running gdal_translate (see http://www. gdal.org/frmt_gtiff.html for creation options), for instance: gdal_translate -co TILED=YES -co COMPRESS=JPEG -co JPEG_QUALITY=80 gebco/bathymetry.tif gebco/bathmetry_tiled.tif would compress the tif using JPEG compression (lossy) set to 80% (low compression, small loss). Cheers Paul On 18-Apr-08, at 5:47 AM, Stefan Schwarzer wrote: Thanks for the info. Did it like this: Library/Frameworks/GDAL.framework/Versions/1.5/Programs/ gdal_translate -outsize 10'' 10'' -co TILED=YES gebco/ bathymetry.tif gebco/bathymetry_tiled.tif But get the message: Segmentation fault Anything that I did wrong? Or should do differently? Thanks for a hint, Stef Hi, Gdal_translate program has options for you: c:\FWToolsgdal_translate Usage: gdal_translate [--help-general] [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/ CInt16/CInt32/CFloat32/CFloat64}] [-strict] [-of format] [-b band] [-outsize xsize[%] ysize[%]] [-scale [src_min src_max [dst_min dst_max]]] [-srcwin xoff yoff xsize ysize] [-projwin ulx uly lrx lry] [-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value] [-gcp pixel line easting northing [elevation]]* [-mo META-TAG=VALUE]* [-quiet] [-sds] [-co NAME=VALUE]* src_dataset dst_dataset By playing with -srcwin or possibly with
Re: [mapserver-users] How to tile a large TIF image?
Paul, I am just interested in one thing: If you say breaking a single tiff file with TILED=YES is equivalent to breaking into individual files and using tileindex - does it mean it is ALWAYS better to create a single file with internal tiling than create many files and using tileindex? I suppose the answer is NO, but I would very like to know why? For example, if you need to display a region that covers several blocks, then in the first case you open just one file, in the second you open several files. So the first method has less overhead right? What are the opposite situations then? thanks very much regards, dejan [EMAIL PROTECTED] wrote on 18.04.2008 13:11:47: Stefan, I don't think you can specify 10'', you need something in pixels. But the command is also probably not what you really need to do. More likely, you should be doing the following: gdal_translate -co TILED=YES gebco/bathymetry.tif gebco/ bathmetry_tiled.tif This will create a single tif file that has an internal block size of 256x256 - you can think of this as having broken your tif up into 256x256 tiles but keeps them all within the same file. This is equivalent to breaking the tif up into individual files, creating a shapefile that has rectangles for each individual files' extent, and using that as a tile index in mapserver. Next, you want to do this: gdaladdo gebco/bathymetry.tif 2 4 6 8 16 this will pre-compute smaller versions of the tif image called overviews at 1/2, 1/4 1/6 1/8 and 1/16 of the original size of the tif - this makes it much more efficient for gdal to return exactly what mapserver is asking for at any given scale. Depending on the range of scales that you need to display your raster image at, you can add more overview levels (or take some away). With this two commands, you can make your rasters much more efficient for mapserver and probably avoid the need to split them up into many files. It does make the file somewhat larger. If you are working with very large files ( 4GB ) then you may run into some problems with tiff and may need to investigate another format or a compressed format like ECW or MrSID (both requiring commercial licenses I believe). There are also some built in compression schemes for tiffs in gdal which you can apply when running gdal_translate (see http://www. gdal.org/frmt_gtiff.html for creation options), for instance: gdal_translate -co TILED=YES -co COMPRESS=JPEG -co JPEG_QUALITY=80 gebco/bathymetry.tif gebco/bathmetry_tiled.tif would compress the tif using JPEG compression (lossy) set to 80% (low compression, small loss). Cheers Paul On 18-Apr-08, at 5:47 AM, Stefan Schwarzer wrote: Thanks for the info. Did it like this: Library/Frameworks/GDAL.framework/Versions/1.5/Programs/ gdal_translate -outsize 10'' 10'' -co TILED=YES gebco/bathymetry.tif gebco/bathymetry_tiled.tif But get the message: Segmentation fault Anything that I did wrong? Or should do differently? Thanks for a hint, Stef Hi, Gdal_translate program has options for you: c:\FWToolsgdal_translate Usage: gdal_translate [--help-general] [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/ CInt16/CInt32/CFloat32/CFloat64}] [-strict] [-of format] [-b band] [-outsize xsize[%] ysize[%]] [-scale [src_min src_max [dst_min dst_max]]] [-srcwin xoff yoff xsize ysize] [-projwin ulx uly lrx lry] [-a_srs srs_def] [-a_ullr ulx uly lrx lry] [-a_nodata value] [-gcp pixel line easting northing [elevation]]* [-mo META-TAG=VALUE]* [-quiet] [-sds] [-co NAME=VALUE]* src_dataset dst_dataset By playing with -srcwin or possibly with -outsize and -projwin you should be able to split your image as you wish. Read more from http://gdal.org/gdal_translate.html -Jukka Rahkonen- Lähettäjä: [EMAIL PROTECTED] [mailto: [EMAIL PROTECTED] ] Puolesta Stefan Schwarzer Lähetetty: 18. huhtikuuta 2008 12:05 Vastaanottaja: mapserver-users@lists.osgeo.org Aihe: [mapserver-users] How to tile a large TIF image? Hi there, I would like to use instead of single large tif image smaller tiles. Although I am well aware of mapserver's and gdal's possibilities to create the shapes for it, I first need to split the large tif into 20 or 50 or 100 tiles. Can anyone give me a hint with what kind of software this is possible? There is an ArcGIS script, but it doesn't work on my machine ( http://arcscripts.esri.com/details.asp?dbid=13978 ). Thanks for any hints, Stef Stefan Schwarzer Lean Back and Relax - Enjoy some Nature Photography http://photoblog.la-famille-schwarzer.de Appetite for Global Data? UNEP GEO Data Portal: http://geodata.grid.unep.ch
Re: [mapserver-users] How to tile a large TIF image?
Dejan - I NEVER say always... But no, I do not think one could state that one approach (a single TIFF with overviews vs. multiple TIFFs at different resolutions) is ALWAYS better than the other. There are a number of factors involved, some of which depend on your filesystem and disk hardware characteristics. There isn't usually much overhead in opening a file by itself. What's expensive is a random seek from one place on a disk to another - that's VERY slow compared to almost anything else you can do on your machine except typing. If you have a single very large TIFF file with overviews, GDAL will need to read the TIFF directory, locate the overview, and seek to it to read it. That's really not very much different from reading the disk directory, locating the file, and seeking to it to read it in the other case. A very large file is more likely to be fragmented into different pieces across the disk, which will require more seeks and slow things down. And a large, fragmented directory will slow down opening one file due to lots of disk seeks. I'm sure there are scenarios in which one approach is faster than the other, but I would not worry about this factor right now. The most important things to do are: 1. Ensure there is a pre-built overview that matches the resolution MapServer will request. If possible, constrain the UI so only a predefined set of view scales are available, and build overviews at all of those scales. 2. Divide large files into tiles (either internally or externally) that are roughly comparable to the size of the output map image, or somewhat larger. It is inefficient to read a lot of small files to create one map, and it is also inefficient to parse through a single large (untiled) file to extract a small area to create one map. 3. Don't put a large number of files into a single directory. On many filesystems this will cause the directory to be stored in separate pieces across the disk, and this can turn just open the file into lots of disk seeks. If you have more than a few thousand files in one directory or so (that's a very rough guess) think about reorganizing the data. - Ed On 4/18/08 10:01 AM, Paul Spencer [EMAIL PROTECTED] wrote: Dejan, I'm not really sure that there would be a noticeable performance difference, and I definitely don't feel sufficiently educated to comment on why one approach would be superior to another, except that I find it much easier to manage single files and to set up a map file to point to them. The process of creating multiple files and setting up a tileindex can be automated, certainly, but I think it is a little more complicated especially for the beginner. The multiple files approach requires mapserver to open the tileindex file, find the polygons that intersect the requested extent, extract the references to the actual disk files, open the disk files using gdal to get the content and composite those into the memory copy of the map. The single file approach does something very similar except that GDAL handles it and there is only one open file handle. I'm not sure that there is a significant overhead to the first approach just on the basis of files needing to be opened. I'm hoping others more qualified will jump in here (Frank, Ed ...) and provide more educated opinions. Cheers Paul On 18-Apr-08, at 9:40 AM, [EMAIL PROTECTED] wrote: Paul, I am just interested in one thing: If you say breaking a single tiff file with TILED=YES is equivalent to breaking into individual files and using tileindex - does it mean it is ALWAYS better to create a single file with internal tiling than create many files and using tileindex? I suppose the answer is NO, but I would very like to know why? For example, if you need to display a region that covers several blocks, then in the first case you open just one file, in the second you open several files. So the first method has less overhead right? What are the opposite situations then? thanks very much regards, dejan [EMAIL PROTECTED] wrote on 18.04.2008 13:11:47: Stefan, I don't think you can specify 10'', you need something in pixels. But the command is also probably not what you really need to do. More likely, you should be doing the following: gdal_translate -co TILED=YES gebco/bathymetry.tif gebco/ bathmetry_tiled.tif This will create a single tif file that has an internal block size of 256x256 - you can think of this as having broken your tif up into 256x256 tiles but keeps them all within the same file. This is equivalent to breaking the tif up into individual files, creating a shapefile that has rectangles for each individual files' extent, and using that as a tile index in mapserver. Next, you want to do this: gdaladdo gebco/bathymetry.tif 2 4 6 8 16 this will pre-compute smaller versions of the tif image called overviews at 1/2, 1/4 1/6 1/8 and 1/16 of the original size of the tif - this
Re: [mapserver-users] Displaying points la
Hi Clement, please answer not directly to me, but to the mapserver list. sorry i forgot you have to add the -o output option to shp2img: shp2img -m aero.map -l AERO -all_debug 5 -o aero.gif Then it should save the image to the specified file. Still no points? Michael 2008/4/18 Clemiou [EMAIL PROTECTED]: Hi! I resolved the last error but now, I don't have error! Just this log and always no points labels. [Fri Apr 18 16:58:43 2008].559339 msOGRFileOpen(mysql.ovf)... [Fri Apr 18 16:58:43 2008].569987 OGROPen(/usr/local/mapserver-5.0.2/../../../var/www/html/ms/map_files/mysql.ovf) [Fri Apr 18 16:58:51 2008].160106 msConnPoolRegister(AERO,mysql.ovf,0x8802e20) [Fri Apr 18 16:58:51 2008].332909 msOGRLayerClose(mysql.ovf). [Fri Apr 18 16:58:51 2008].332975 msOGRFileClose(/usr/local/mapserver-5.0.2/../../../var/www/html/ms/map_files/mysql.ovf,0). [Fri Apr 18 16:58:51 2008].332989 msConnPoolRelease(AERO,mysql.ovf,0x8802e20) [Fri Apr 18 16:58:51 2008].333003 msConnPoolClose(mysql.ovf,0x8802e20) [Fri Apr 18 16:58:51 2008].333276 msDrawMap(): Layer 1 (AERO), 7.953s [Fri Apr 18 16:58:51 2008].30 msDrawMap(): Drawing Label Cache, 0.000s [Fri Apr 18 16:58:51 2008].48 msDrawMap() total time: 7.976s GIF87a� � ��,� � ��ŋ���ڋ���H�扦�ʶ� ��L ËçôºýŽÏë÷ü¾ÿ (8HXhx��� )9IYiy *:JZjz�� K[k{ ,N^n~ /?O_⎺ Ÿ¯¿Ïßïÿ0��� ]m}�� ¼øñäË›? ½úõìÛ» ?�ۿ扤��`� H◆ ˆ◆‚~ .È◆ƒ▒„ NH▒… ^ˆ▒† ┼È▒‡ ␉ˆŽH␉‰žˆ␉Š*(R)È␉‹.␌Œ2ÎH␌ 6Þˆ␌Ž:îÈ␌ ␍ BI␍'F ‰␍'J.É␍N␊RNI␊•V^‰␊–Z┼É␊—^· °˜␉ŽI°™°ž‰°š┘(R)É°›┼¾±œ⎼ÎI± ┴Þ‰±ž≥îɱŸ·þ  ‚J¡† Š¢Š.Ê£Ž ␋¤'NJ␋¥–^Š␋¦š┼Ê␋§ž· ┘¨¢ŽJ┘(c)¦žŠ┘ªª(R)Ê┘«(R)¾ ┐¬²ÎJ┐¶ÞŠ┐(R)ºîÊ┐¯¾þ ┌°ÂK┌±Æ ‹þ┌²Ê.Ë┌³Î └´ÒNK└µÖ^‹└¶Ú┼Ë└·Þ· ┼¸âŽK┼¹æž‹┼ºê(R)Ë┼»î¾ ⎺¼òÎK⎺½öÞ‹⎺¾úîË⎺¿þþ ⎻À Lp� �p [Fri Apr 18 16:58:51 2008].616500 msFreeMap(): freeing map at 0x87297f8. === Thanks for your help Clement Date: Fri, 18 Apr 2008 15:45:17 +0200 From: [EMAIL PROTECTED] To: [EMAIL PROTECTED]; mapserver-users@lists.osgeo.org Subject: Re: [mapserver-users] Displaying points lab Hi Clement, please leave the thread also on the list... Ok, that is not the real debug output, merely the log. Did you try it from the command line with shp2img: shp2img -m aero.map -l AERO -all_debug 5 this should print some debug information to stdout. Are you using mapserver 5? Then you can set a special error log for each mapfile, with the directive CONFIG MS_ERRORFILE /path/to/mserror.log at the map level. Cheers, Michael 2008/4/18 Clemiou : Hi The one thing that I can see in the mapserver log is: Fri Apr 18 15:18:51 2008,18284,127.0.0.1,MS,3,1.241143 46.540217 1.445917 46.744991,-1.00 -1.00,FRANCE ,normal execution ^^ Very funny isn't it? Clement Date: Fri, 18 Apr 2008 15:03:59 +0200 From: [EMAIL PROTECTED] To: [EMAIL PROTECTED]; mapserver-users@lists.osgeo.org Subject: Re: [mapserver-users] Displaying points labe Hi, could you have a look at the debug output? Turn debugging on with DEBUG on on the map and on the layer level, then you should get debug output, either to your configured mapserver log or the webserver (error) log. Another way woudl be to use shp2img. Cheers, Michael 2008/4/18 Clemiou : Hi, I have modifief my layer for just have the points, as you can see: # # Le layer des aérodromes . # # LAYER NAME AERO TYPE POINT STATUS ON PROJECTION init=epsg:27582 END TEMPLATE ttt_query.html CONNECTIONTYPE OGR CONNECTION mysql.ovf #TOLERANCE 5 #TYPE ANNOTATION #LABELITEM nom DATA my_layer CLASS NAME my_class #LABEL # ANGLE auto #TYPE TRUETYPE #ANTIALIAS TRUE #MINFEATURESIZE 100 #COLOR 255 255 255 #OUTLINECOLOR 255 0 0 #SIZE 20 #FONT arial #FORCE TRUE #END # Fin LABEL AERO STYLE SYMBOL 'circle' SIZE 15
RE: [mapserver-users] Displaying points labels
Hi! Still no points! But I found perhaps something... My datas are in WGS84 format and my map in Lambert II. So, in my mapfile, my Map object has WGS84 projection, the layer of my Map, Lambert II format and the layer of my points the WGS84 format. But still no points. I don't know why. ^^ Thanks for your help. Clement _ Découvrez les profils Messenger de vos amis ! http://home.services.spaces.live.com/___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
Re: [mapserver-users] Displaying points label
Hi, in your ovf as well as in the layer definition you claim your data to be in EPSG:27582, if it is actually in epsg:4326 then it should say so in both places. The layer projection defines in which system your data is in, the map projection defines in what system your image is rendered (and your data projected to on the fly if necessary). Cheers, Michael 2008/4/18 Clemiou [EMAIL PROTECTED]: Hi! Still no points! But I found perhaps something... My datas are in WGS84 format and my map in Lambert II. So, in my mapfile, my Map object has WGS84 projection, the layer of my Map, Lambert II format and the layer of my points the WGS84 format. But still no points. I don't know why. ^^ Thanks for your help. Clement _ Découvrez les profils Messenger de vos amis ! http://home.services.spaces.live.com/___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users -- --- Michael Schulz [EMAIL PROTECTED] in medias res Gesellschaft für Informationstechnologie mbH In den Weihermatten 66 79108 Freiburg Tel +49 (0)761 556959-5 Fax +49 (0)761 556959-6 http://www.webgis.de / http://www.zopecms.de --- ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] variable substitution
Hello, i want to use variable suistitution in a map but i dont understand hoe i could replace teh varaible i n the template file, im using postgis for load the shape but i wnt to change the expression via variable sustitution. Help me please and thanks _ Connect to the next generation of MSN Messenger http://imagine-msn.com/messenger/launch80/default.aspx?locale=en-ussource=wlmailtagline___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
RE: [mapserver-users] MapServer v. 4.10.0 featureid problem
Tom, Perhaps I should tell you why I'm after the (or any) unique AND sequential id field for a service's features. I'm trying to page the feature data in a UI and started designing my solution after having read the WFS v 1.0.0 spec. I assumed (yes, I know what happens when one does that (;)) that I'd just use feature ids. I should have asked what mechanism MapServer employs to allow me to do something like this (the partial retrieval of features from a large dataset on some service somewhere). I'm going to bring back a certain number of features with each call to a WFS and as the user scrolls through this data I will make subsequent calls to the same service and ask for the next batch of features. Does this make sense and can you point me in the right direction for getting MapServer to do this? Thanks for your patience, Brian -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Kralidis,Tom [Burlington] Sent: Friday, April 18, 2008 5:36 AM To: Delfos, Jacob; mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Jacob is bang on (I'm no expert on ArcGIS :) If I invoke a tool like ogrinfo on the shapefile, I see no attribute called FID. I do see ID and DISTRICT attributes. Try setting gml_featureid ID in your mapfile, which should pick up the field and process accordingly. ..Tom -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Delfos, Jacob Sent: 17 April, 2008 11:55 PM To: mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Hi Guys, I'm no expert on WFS, but I do know ArcGIS. The FID field is not really an attribute. It is more like an internal ESRI unique row ID (which gets updated after edits to be sequential). If you open up the attributes in another application (if you use shapefiles, open up the DBF in OpenOffice Calc or Excel), then you will find there is no FID field. Within ArcGIS, however, the FID field is presented as if it were a normal attribute. Hope this helps, Jacob -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Brian Mixon Sent: 18 April 2008 01:34 To: Kralidis,Tom [Burlington]; mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Tom, As soon as the Brenton on this list gets in to work we'll try to put something together for you. I'm working remotely and need his GIS knowledge to get something for you. Brian -Original Message- From: Kralidis,Tom [Burlington] [mailto:[EMAIL PROTECTED] Sent: Thursday, April 17, 2008 10:51 AM To: Brian Mixon; mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Hi Brian, What's your underlying data? If there's anyway you can package a small testcase, I can do a local test. ..Tom -Original Message- From: Brian Mixon [mailto:[EMAIL PROTECTED] Sent: 17 April, 2008 12:24 PM To: Kralidis,Tom [Burlington]; mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Hi Tom, Yes, it does. It's an Arc field that is present in all our data. Thanks, Brian -Original Message- From: Kralidis,Tom [Burlington] [mailto:[EMAIL PROTECTED] Sent: Thursday, April 17, 2008 10:11 AM To: Brian Mixon; mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: RE: [mapserver-users] MapServer v. 4.10.0 featureid problem Does the column/field fid exist in your data? ..Tom From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Brian Mixon Sent: 17 April, 2008 11:51 AM To: mapserver-users@lists.osgeo.org Cc: Brenton Wonders; Russell Huff Subject: [mapserver-users] MapServer v. 4.10.0 featureid problem Hello, I'm a MapServer newbie and am having trouble with a WFS v. 1.0.0 request. I'm trying to have all the features for a particular layer returned (which happens successfully) with the fid (feature id) included as well (does not happen). I've included the gml_featureid 'fid' tag in the layer's metadata section of the mapfile but only receive the following warning the the xml/gml response:!-- WARNING: FeatureId item 'fid' not found in typename 'senate_districts'. -- I've searched and read the various tickets that claim this problem is fixed so I'm wondering what I'm doing wrong? Any information to point me in the right direction would
Re: [mapserver-users] Investigating OpenLayers - MapServer (WMS) - OGR - db
On Fri, Apr 18, 2008 at 12:31:24PM -0500, P Kishor wrote: SELECT wkt_geometry FROM table WHERE extent of features is in the BBOX That is what I want to investigate. What exactly is OGR/MapServer executing? What is the form of WHERE extent of features in the BBOX? How can I log this level of detail? You'll likely need to modify OGR to log this level of detail. In general, this type of logging is done by the 'server' in DB environments: Since SQLIte doesn't *have* a server, I don't think there's any way to do this external to OGR, and I don't believe that it's an existing feature of the OGR sqlite driver. In any case, MapServer is the wrong list for this: gdal-dev is probably the right place to go. Regards, -- Christopher Schmidt MetaCarta ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] Experiences with CartoWeb
Hi there, we've used mapserver a lot over the last years, in some minor as well as some major applications. Now, we're thinking to replace a mapping module in a major portal application with cartoweb, in order to avoid a duplication of the development of code, of having a more modular approach, to facilitate update and maintenance etc. What are your experiences in using Cartoweb? What are, upon your opinion, the major advantages or drawbacks? Anything we should be aware of before undergoing such a change in our application? Thanks for any advices, Stef ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] redrawing the map
right now, i have a nearly functional code. This code is for zooming to a predetermined area. the problem is that the map doesn't get redrawn when i use this script. I found this script in the archive, and would like to know is there a line i have to add in order to get it to redraw or refresh the map. More importantly if this script is actually correct, any help is much appreciated. John html head script language=javascript function AreaSelected() { alert(document.myform.ViewArea.value); document.myform.target=; document.myform.mode.value=browse; document.myform.imgext.value = (document.myform.ViewArea.value); document.myform.submit(); } /script /head body form method =POST name=myform select name=ViewArea onChange=AreaSelected(); optionZoom to Region/option option/option option value=5457203 487104 5451873 482304Indian Reserve 2/option option value=5435308 491431 5434498 490403Indian Reserve 4/option option value=5425919 541907 5502970 444168Musqueam Statement Of Intent/option option value=5435189 475830 5457751 504449Richmond Boundary/option option value=5448859 481843 5463269 500122Vancouver Boundary /option option value=5432692 482228 5457268 515426North Arm Fraser River /option /select input name=mode input type=hidden name=imgxy value=299.5 299.5 input type=hidden name=imgext value=[mapext] /form /body /html ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
RE: [mapserver-users] Displaying points label
Hi, Now I have the good projections I think but I don't see the points.I have: WGS84 for my map.EPSG:27582 for my map Layer WGS84 for my label layer. What do you think about that? Is it Ok with my datas (WGS84 for database datas and Lambert II for my Map). Thanks for you answer Clement Date: Fri, 18 Apr 2008 18:25:01 +0200 From: [EMAIL PROTECTED] To: mapserver-users@lists.osgeo.org Subject: Re: [mapserver-users] Displaying points label Hi, in your ovf as well as in the layer definition you claim your data to be in EPSG:27582, if it is actually in epsg:4326 then it should say so in both places. The layer projection defines in which system your data is in, the map projection defines in what system your image is rendered (and your data projected to on the fly if necessary). Cheers, Michael 2008/4/18 Clemiou [EMAIL PROTECTED]: Hi! Still no points! But I found perhaps something... My datas are in WGS84 format and my map in Lambert II. So, in my mapfile, my Map object has WGS84 projection, the layer of my Map, Lambert II format and the layer of my points the WGS84 format. But still no points. I don't know why. ^^ Thanks for your help. Clement _ Découvrez les profils Messenger de vos amis ! http://home.services.spaces.live.com/___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users -- --- Michael Schulz [EMAIL PROTECTED] in medias res Gesellschaft für Informationstechnologie mbH In den Weihermatten 66 79108 Freiburg Tel +49 (0)761 556959-5 Fax +49 (0)761 556959-6 http://www.webgis.de / http://www.zopecms.de --- ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users _ Créez votre disque dur virtuel Windows Live SkyDrive, 5Go de stockage gratuit ! http://www.windowslive.fr/skydrive/___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] How best to serve very large (300GB) raster data?
I've been using ESRI products for quite a while. To serve a very large raster dataset I have (about 300 GB), or really any raster or vector, I usually just load the data into ArcSDE installed on Oracle and server through ArcIMS. However, I want to move away from ESRI/Oracle and into a more open environment, with all components being open source if possible. How best to serve a large raster dataset of about 300 GB? My data is in 2010 tiles, each a 3-band TIF of about 135 MB. From my initial research, I like what MapServer can do and all of the potential clients (OpenLayers, ka-Map, and many more.) I also really like PostGIS for vectors...and have already started a project with PostGIS points, being maintained by Quantum GIS and mapped through Google Maps Javascript API and some AJAX/PHP code I wrote. My biggest worry is dealing with the very large rasters. Ideally, a tile server would be nice but I don't want to create the tiles through an ESRI product (like Arc2Earth or HTMLImageMapper.) I really like the possibility of TileCache (http://www.tilecache.org/) working with MapServer. The immediate question is what would be the raster data source to MapServer? Should I create tiles inside each TIF using gdal_translate and then a shapefile tileindex? Maybe go with compressed JPGs instead of TIFs since they usually have much smaller file sizes? I'm not sure where to start. Thanks for any guidance you can offer. - John John Callahan Geospatial Application Developer Delaware Geological Survey University of Delaware 227 Academy St, Newark DE 19716-7501 Tel: (302) 831-2833 Email: [EMAIL PROTECTED] http://www.dgs.udel.edu ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users