Has anybody done performance benchmarks between C++, C#, and Java ArcGIS Server Object Extensions? I'm wondering writing lower-level code translates to better performance.


title: How to create a custom enterprise QGIS installer with plugins and symbols?

text: I am looking at installing QGIS at a company level. I would like my installer to install specific plugins, and add our company map symology. Are there any ways to do this? I have seen post about the custom GUI installer, but I didn't feel it was the same as my request. Thanks...


title: Viewshed / Line of Sight analysis- as a batch process?

text: I wish to carry out a line of site analysis similar to the subject in here (preferably within QGIS): Visual Impact Assessment in QGIS or GRASS **Scenario:** I have two points, A and B. Point A is a signal provider and Point B is a signal receiver. Point A is at **X** altitude and Point B is at **Y** altitude. There may be hills between Points A and B (represented by either vector contour lines or a raster DTM) that block Point B from receiving signal. **Data:**   * A vector shapefile containing hundreds of points (the viewing position / Point A) ![Viewing Position data](http://i.stack.imgur.com/M3xqH.jpg)   * A vector shapefile containing several thousands of points (each point in this shapefile represents a position that needs to be viewed / Point B)   * A 50 m resolution DTM in raster format (currently looking for higher resolution free data for the UK) ![DTM, Points A\(black dots\) and Points B \(blue dots\)](http://i.stack.imgur.com/gG0mO.jpg) **Question:** How can I establish that Point As are able to view all of Point Bs (via say, a mask?). And how do I run this as a batch process? At the moment I can only run the line of sight (LOS) analysis for one coordinate at a time: ![r.los](http://i.stack.imgur.com/lKqsZ.jpg)


title: raster-derived KMZ cannot display in google map api v3

text: I got some raster(tiff)-generated KMZ files, but they cannot be displayed at all in googe map api v3. However, all of my shapefile-generated KMZ files can be displayed properly. Anyone faced the similar situation before?


title: how to do change detection in ERDAS

text: I'd like to ask a really simple question. I want to obtain the size of unchanged, decreased, and increased area etc (just like shown in image below) from 2 classification images in 2000 and 2007. How to do this in ERDAS? ![result of land cover change detection](http://i.stack.imgur.com/7p6TK.png). Thank you so much


title: Which commercial off-the-shelf GIS server products work on Linux?

text: I am looking for suitable COTS GIS server product for Linux platform. The product should allow 500 Clients concurrently and should have Java / C++ development APIs. So far I found ESRI ArcGIS server which supports for Linux. I need few more products which have good performance like ESRI. Kindly share the list of available COTS (Commercial off the shelf) GIS server products.


title: Basemap Gallery Error - ESRI JS API

text: I'm getting an error when trying to switch basemaps in my application(right now only works in IE9). **TypeError: Unable to get value of the property '_basemapGalleryLayerType': object is null or undefined** I'm using standard code provided in the ESRI samples. The full code is a mashup of various samples, so my error may be in the order I've placed things. However, I've tried reordering the code and haven't had any luck. When I remove `mapReady` from the `dojo.connect(map, "onLoad", mapReady);` listener, the basemap gallery works, however it causes my initial basemap to fail and the mapReady function to fail as well (obviously). There is no link that I can determine between the initial basemaps and my mapReady function so I'm not sure what the deal is there... Also, if I clear the browser cache the entire app will work, 1 time only, then on refresh, it fails as described above. Very inconsistent behavior. Any help would be greatly appreciated.


title: OpenLayers OSM tile boundary

text: I have a local tile source of OSM tiles. I'd like to show the tile zoom, x, y on the tile so that I can optimize and create small offline groups of tile for specific areas(polygons). I can get the tile boundary to show by setting a negative gutter, but I need the tile url info (zoom, x, y) to help with the manual optimization. I saw an online resource a while back that showed this exact information on a map, but I've lost it :( I know this is probably a simple problem to solve, but I haven't found any solutions after a while of playing around with it and many searches. Here is my layer               var newLayer = new OpenLayers.Layer.OSM("GC-h", "tiles/GC/composite_h/${z}/${x}/${y}.png", {numZoomLevels: 19, gutter: -1, isBaseLayer: false});     map.addLayer(newLayer);      TIA!


title: Is the PostGIS documentation correct for ST_Project?

text: In example of the page ST_Project the select is               SELECT ST_AsText(ST_Project('POINT(0 0)'::geography, 100000, 45));                       st_astext         ------------------------------------------          POINT(0.63523102912532 0.63947233472882)         (1 row)      But, in my enviroment this code give me an error like "Azimuth must be between -2PI and 2PI". After a few tryings I solve with :               SELECT ST_AsText(ST_Project('POINT(0 0)'::geography, 100000, radians(45) ))      The sample in page tells the data type expected is in degrees, and the function expects radians. May I using wrong PostGIS version? Thanks. BTW : My PostGIS version : 2.0


title: File naming for Enterprise GIS database?

text: I'm designing the guidelines for an Enterprise GIS database. Part of this database will include customized geoprocessing tools using a variety of software. I'm interested in knowing whether anyone can foresee a problem with including '_' as a separator for filenames if GRASS GIS is being used as the main geo- processing software? Most of the tools will be coded in python, SAGA GIS, and R stats will also likely be used as background processess. Thanks for any advice!


title: How to geocorrect in qgis?

text: I need to fit a utm wgs84 tiff raster to a topomap in the same projection. I use GCP and it invariably fails. The first raster uncorrected rater superimpose well, the "geocorrected" raster is projected far away. I use the same projection for the project. What's happening? * * * Ok, sorry and thanks for your answer. I have a utm wgs84 image in tiff and toposheet in utm wgs84 in tiff. The image need geocorrection. I open the toposheet in the main qgis screen and I import the image in the georeferencer. I create GCP with the "Enter map coordinates" tool, "From map canvas". This fails: srcX will be almost the same as dstX, the absolute value of srcY will also be very close to the dstY, but srcY is negative. The new geocorrected map is projected in a totally different area.


title: How to get elem_info array from WKT or JTS Geometry in Java?

text: I have Well Known Text (WKT) string with geometry. Now I want to create Oracle Mapviewer SdoGeometry, which needs sdo_elem_info array as one of arguments which is an array of doubles. How can I extract elem_info array from WKT String? I've read that geometry into Geometry object from Vivid Solutions JTS library, but I can't find a way to get elem_info from that either. Anyone can help?


title: Is it possible to lock overview map?

text: Issue: I've got openlayers application with some overview map which has just one scale so the size of it doesn't change even though I zoomed my map. Now i would like to make overview map locked - blocked, and in the same time to make possible to move the extent rectangle on the overview map, so it's not always in the center. Summing up I want to switch properties between the extent rectangle on the overview map and overview map. I will be grateful for your help


title: Why do Print composer templates not save my maps in their entirety?

text: I have been using QGIS and its various updates in a limited way to make maps for a few years now, but I have always had the same problem with the Print Composer. I make my layered maps with grouped annotations and save them as templates, which I'd always assumed would allow me to recall them at another time, when the map would pop up just as I'd left it and if I wanted to make minor changes or adjustments, I could then do so. This has never happened. What does happen is that when I open Print Composer files, the annotations appear, but not the original map images. I have looked at various QGIS manuals but cannot find the answer to this. Would someone be able to explain to me if there is a way to save the complete content of maps for retrieval and/or editing at a future date? To me, this would appear to be the main function of a GIS, so I must be misunderstanding how it works.


title: ArcGIS Network Analyst tool for Android?

text: I am trying to include the Network Analyst tool into an Android application. I have published the tool on ArcGIS Server and I integrate its URL address to use it. On the same time, I import a MapServer which display my network with all the layers needed to make NA work. I tried to create a GeoProcessor and to input the tool in it and then execute it, but it doesn't work and I haven't found any clues on the web... Does anyone has an idea about how to deal with this ?


title: How can I convert KML files to shapefiles? (ArcView license)

text: I only have an ArcView 9.3 license at work. What's the best way to convert a large KML file to SHP format? I tried the DNRGarmin tool to import, but the file was too large for that particular program. I've also tried XTools and received an error right as the process was starting. Does anyone have a reliable method?


title: Flood records/database in Europe

text: Anyone knows where I can find a complete* database with the historic flood occurrence in Europe, at least since the 70's ? I am looking for some information on the date, number of floods* _, amount precipitation_ *, as well as the location*** (the region and/or the affected city will be nice). *complete in time(from the 70's or 80's until now) * _both informations are not extremely important *_ *a name will be enough I am not looking for a database with the geographic coordinates of the flood extent, just a reference for the name of the region/city


title: Replacing header in Flexviewer

text: Will the index.html swf need to be re-compiled in order to embed the flexviewer app in between a html header and footer wrapper? Or it this perhaps handled in the HeaderControllerWidget.mxml?


title: ArcGIS Server Rest API Export with Class Breaks

text: I'm trying to export a map of a dynamic layer with class breaks. I can export the dynamic layer like this:               $mapURL = 'http://sampleserver6.arcgisonline.com/arcgis/rest/services/Census/MapServer/export?bbox=-183.78001472868405,16.300709121618663,-61.406854669684265,74.03030803096895&bboxSR=&layers=&layerDefs=&size=&imageSR=&format=png&transparent=false&dpi=&time=&layerTimeOptions=&dynamicLayers=[{%22id%22:101,%22source%22:{%22type%22:%22mapLayer%22,%22mapLayerId%22:3},%22drawingInfo%22:{%22renderer%22:%20{%22type%22:%22simple%22,%20%22symbol%22:%20{%22type%22:%22esriSFS%22,%22style%22:%22esriSFSSolid%22,%22color%22:[255,0,0,255],%22outline%22:{%22type%22:%22esriSLS%22,%22style%22:%22esriSLSSolid%22,%22color%22:[0,255,0,255],%22width%22:1}}}}}]&gdbVersion=&f=image';      // $mapURL = 'http://sampleserver6.arcgisonline.com/arcgis/rest/services/Census/MapServer/export?bbox=-183.78001472868405,16.300709121618663,-61.406854669684265,74.03030803096895&bboxSR=&layers=&layerDefs=&size=&imageSR=&format=png&transparent=false&dpi=&time=&layerTimeOptions=&dynamicLayers=[{"id":101,"source":{"type":"mapLayer","mapLayerId":3},"drawingInfo":{"renderer": {"type":"simple", "symbol": {"type":"esriSFS","style":"esriSFSSolid","color":[255,0,0,255],"outline":{"type":"esriSLS","style":"esriSLSSolid","color":[0,255,0,255],"width":1}}}}}]&gdbVersion=&f=html' This works great. However I am having a difficult time applying a classBreaksDef with "classificationField": "POP2000", "classificationMethod": "esriClassifyNaturalBreaks", "breakCount": 5 . Does anyone have an example I can peek at to look at the syntax? Thanks, Tyler


title: Why is the dynamic text on data driven pages not displaying a value from the grid attribute table?

text: Relating to my question here, I've been trying to return a value from the attribute table using the following code: > `<dyn type="page" property="lu1p"/>` However, even if I try to use a different attribute heading (i.e. lu2p), it only displays "[empty]". Is there something which I'm missing? The grid index page is actually a set of irregular study areas (i.e. not perfectly square), but it hasn't effected the use of Data Driven Pages. I have looked at the help files for this but have not found any solution. Edit: As requested, here is a screenshot: ![enter image description here](http://i.stack.imgur.com/SudPM.jpg) (http://i.imgur.com/8FXwS.jpg) available here also Edit: I'm still having this issue so if anyone could provide an alternative answer, I'd really appreciate it. I tried again in a new mxd and it worked once but then it stopped worked when I re-opened it.


title: Proportional Symbol Formatting

text: I need to create a proportional symbol map for a GIS class. I am wondering, how can I manually change the symbology to make the dots grow in size depending on their magnitude represented, and also how can I manually change the color scheme?


title: Adding Descriptions to custom ArcToolbox 9.3.1

text: I looked http://webhelp.esri.com/arcgiSDEsktop/9.3/index.cfm?TopicName=Setting_script_tool_parameters and couldn't find how to add Descriptions for each parameter in the script tool. I see how to add a description for the overall script tool but not for each parameter.


title: How to update data from external table in real time?

text: I have been trying for a long time to connect a external table data[csv] through Excel (for exmaple) with QGis with the objective to each time i change any data in table(for example data location - coordinates [X ; Y]),this changes occur in QGis real time.Is there any procedure to get it?. Thank you very much


title: How to create a dynamic, interactive map with QGIS?

text: How can we create an interactive map with QGIS? I would keep it in the bottom of aerial photographs. When the mouse moves over elements of a layer, I wish they will be highlighted and a window shows the information contained in the attribute table of my shp layers. The export permit is done in a good resolution to zoom in and zoom out dynamically. Thank you in advance for your help ! * * * Comment peut-on créer une carte interactive avec Qgis ? je souhaiterais conserver en fond des photos aériennes. Lorsque la souris passe sur des éléments d'une couche, je souhaiterais qu'ils apparaissent en surbrillance et qu'une fenêtre indique les informations contenues dans la table attributaire de mes couches shp. Que l'export permette se fasse dans une bonne résolution pour pouvoir zoomer et dé-zoomer de manière dynamique. Merci d'avance pour votre aide!


title: Problem with extent of mapserver wms/layer

text: I use mapserver with MS SERVER 2008. The problem with the correctly computing layer extend of Map server. If map file like this                MAP           NAME "test map"           IMAGETYPE      GIF           EXTENT         19.63 41.10 180.00 81.6           SIZE           1900 800           IMAGECOLOR     255 255 255           CONFIG "PROJ_LIB" "C:\inetpub\wwwroot\MapServ\proj\nad"               WEB         IMAGEPATH "img_tmp/"         IMAGEURL "img_tmp/"         METADATA             "wms_enable_request" "*"             "wms_title" "test"             "wms_onlineresource"  "******"             "wms_srs"  "epsg:4326"             "wms_abstract" "test"             "wms_encoding" "windows-1251"             "wms_extent" "19.63 41.10 180.00 81.6"         END     END          PROJECTION         "init=epsg:4326"     END          LAYER         TEMPLATE "template/template.html"         NAME "okato"         METADATA             "wms_title" "OKATO"             "wms_srs"  "epsg:4326"             "wms_include_items" "all"             "wms_abstract" "test"         END         TYPE POLYGON         STATUS ON       ####         CONNECTIONTYPE PLUGIN         PLUGIN "C:/inetpub/wwwroot/MapServ/specialplugins/msplugin_mssql2008.dll"         CONNECTION "Server=... ; uid=... ; pwd=... ; Database=... ; Integrated Security=true"         DATA "geom from dbo.okato USING UNIQUE ID USING SRID=4326"       ####         PROJECTION             "init=epsg:4326"         END      Map server generate XML with the following BoundingBox of layer:               <CRS>epsg:4326</CRS>     <EX_GeographicBoundingBox>     <westBoundLongitude>-3.40282e+038</westBoundLongitude>     <eastBoundLongitude>3.40282e+038</eastBoundLongitude>     <southBoundLatitude>-3.40282e+038</southBoundLatitude>     <northBoundLatitude>3.40282e+038</northBoundLatitude>     </EX_GeographicBoundingBox>     <BoundingBox CRS="epsg:4326" minx="-3.40282e+038" miny="-3.40282e+038" maxx="3.40282e+038" maxy="3.40282e+038"/>      If I add extent layer like this:               LAYER     ...         EXTENT         19.63 41.10 180.00 81.6     ...     END      or               LAYER     ...         METADATA             ......             "wms_extent" "19.63 41.10 180.00 81.6"         END     ...     END      It works correctly               <CRS>epsg:4326</CRS>     <EX_GeographicBoundingBox>     <westBoundLongitude>19.63</westBoundLongitude>     <eastBoundLongitude>180</eastBoundLongitude>     <southBoundLatitude>41.1</southBoundLatitude>     <northBoundLatitude>81.6</northBoundLatitude>     </EX_GeographicBoundingBox>     <BoundingBox CRS="epsg:4326" minx="41.1" miny="19.63" maxx="81.6" maxy="180"/>      **How to make that mapserver dynamically to calculate the extent (and BBOX) of layer, and not add it manually?**


title: ESRI World Light Gray Base Map JavaScript - Prevent zoom out no more than the full extent

text: It is my first attempt to use JavaScript and the ESRI products. I am using the World Light Gray Base map for my first web GIS application. I have a question about the map full extent method. I would like to set the full extent size so users can't zoom out no more than my full extent setting. So the map always fully covers the map window. At the moment, if I zoom out the World Light Gray Base Map, the full size of the map is getting smaller so it covers a small part of my map window. Could anyone helps me how to set up the maximum zoom out size so my map window is always covered by the world light gray map? Best,


title: Why does QGIS hang when I try to create a vector grid?

text: I am working in open layers plugin.i need to project the satellite data with respect to vector data shape file i am having .for that Andre Joost suggested to create vector grid with vector data reference system .But when I am creating vector grid the system hangs this the problem?


title: Adding vector layers on top of gdal rasterized image

text: I used gdal2tile.py on a big image (actually its screenshots of image map) and now I have a nice openlayers-powered map I can pan and zoom around. But now I want to add a vector layer on top so I can highlight certain features on the map. I tried following the OpenLayers recipes but have failed. Vector Formats at best shows nice image on top zoom, shows only 3 tiles when zoomed in once and breaks down completely later on. GML Layer just shows empty map. I imagine the problems I'm having is because those examples are expecting real Earth coordinates. I only want to add a bunch of square boxes all over the map, how can I do this? This is the JS code that was generated by gdal:               var map;     var mapBounds = new OpenLayers.Bounds(0.0, -4722.0, 4150.0, 0.0);     var mapMinZoom = 0;     var mapMaxZoom = 5;     var emptyTileURL = "http://www.maptiler.org/img/none.png";     OpenLayers.IMAGE_RELOAD_ATTEMPTS = 3;          function init()         {         var options = {             div: "map",             controls: [],             maxExtent: new OpenLayers.Bounds(0.0, -4722.0, 4150.0, 0.0),             maxResolution: 32.000000,             numZoomLevels: 6         };         map = new OpenLayers.Map(options);              var layer = new OpenLayers.Layer.TMS("TMS Layer", "",                 {                     serviceVersion: '.',                     layername: '.',                     alpha: true,                     type: 'png',                     getURL: getURL                 });              map.addLayer(layer);         map.zoomToExtent(mapBounds);              map.addControls([new OpenLayers.Control.PanZoomBar(),             new OpenLayers.Control.Navigation(),             new OpenLayers.Control.MousePosition(),             new OpenLayers.Control.ArgParser(),             new OpenLayers.Control.Attribution()]);         }          function getURL(bounds)         {         bounds = this.adjustBounds(bounds);         var res = this.getServerResolution();         var x = Math.round((bounds.left - this.tileOrigin.lon) / (res * this.tileSize.w));         var y = Math.round((bounds.bottom - this.tileOrigin.lat) / (res * this.tileSize.h));         var z = this.getServerZoom();         var path = this.serviceVersion + "/" + this.layername + "/" + z + "/" + x + "/" + y + "." + this.type;         var url = this.url;         if (OpenLayers.Util.isArray(url))             {             url = this.selectUrl(path, url);             }         if (mapBounds.intersectsBounds(bounds) && (z >= mapMinZoom) && (z <= mapMaxZoom))             {             return url + path;             } else             {             return emptyTileURL;             }         }          function getWindowHeight()         {         if (self.innerHeight) return self.innerHeight;         if (document.documentElement && document.documentElement.clientHeight)             return document.documentElement.clientHeight;         if (document.body) return document.body.clientHeight;         return 0;         }          function getWindowWidth()         {         if (self.innerWidth) return self.innerWidth;         if (document.documentElement && document.documentElement.clientWidth)             return document.documentElement.clientWidth;         if (document.body) return document.body.clientWidth;         return 0;         }          function resize()         {         var map = document.getElementById("map");         var header = document.getElementById("header");         var subheader = document.getElementById("subheader");         map.style.height = (getWindowHeight() - 80) + "px";         map.style.width = (getWindowWidth() - 20) + "px";         header.style.width = (getWindowWidth() - 20) + "px";         subheader.style.width = (getWindowWidth() - 20) + "px";         if (map.updateSize)             {             map.updateSize();             }         ;         }          onresize = function ()         {         resize();         };


title: How can I create a legend with the new data defined symbols functionality?

text: I have a large geology data set that I need to display using an industry standard set of symbols and colours. The symbols and colours are attributes of the data. I have been using the data defined symbols capability of the dev build of QGIS (code revision b805894) to display the data according to its attributes. This works perfectly in terms of displaying the data correctly in the map window. The problem arises when I try to create a legend in the composer. The data- defined colours displayed on the map don't pull through to the legend. Here are some screens shots. The first shows the map window. You can see the problem manifesting in the Layers box on the left. The symbols are all coloured black. The second screen shot shows the composer window with a legend. ![enter image description here](http://i.stack.imgur.com/9MNw9.png) ![enter image description here](http://i.stack.imgur.com/doOhc.png) Thanks in advance


title: QGis creating shape file of several points out of two Excel sheets

text: I am about to create a Qgis file on which it should be visible which companies are working together. As source I have got two Excel Tables.   1. Table = Distribution Centers   2. Table = Point of sales I like to connect The point of sales (always several) with the associated distribution center. In the tables of the distribution centers you have the ID, and some (x/y) coordinates. In the Point of sale table you can find a ID, also (x/y) coordinates + a ID and which allows you to connect the Point of sales and the distribution centers (the second ID is the same like any of the distribution centers). I want to separate the group of the distribution center and the related point of sales as 1 shape file without searching manually in my excel tables and creating 1 shape manually for every group... I hope you will be able to understand my problem although my english is horrible. * * * My Lines are like this...               TDLinx-DC|TDLinx_-POS|Geometrie      00421362|00421362|LINESTRING (2.30567 48.7321)      00421362|00611434|LINESTRING (0.89378 46.4288)     00421362|00409507|LINESTRING (6.45527 46.359952)       Now it is like this. Might it be impossible to have a mixed point AND polylines layer?               TDLinx-DC|TDLinx-POS|Geo     00421102|00617429|LINESTRING (6.06043 48.70638, 7.34101 47.7495)     00421102|00256888|LINESTRING (6.06043 48.70638, 6.1812365 48.6876425)     00421102|00264459|LINESTRING (6.06043 48.70638, 5.0379298 47.3211782)     00421102|00264472|LINESTRING (6.06043 48.70638, 6.4946471 48.591343)     00421102|00264483|LINESTRING (6.06043 48.70638, 6.1618809 48.6763621)     00421102|00264501|LINESTRING (6.06043 48.70638, 7.7577037 48.5677807)     00421102|00264502|LINESTRING (6.06043 48.70638, 5.5904536 47.4484859)     00421102|00264503|LINESTRING (6.06043 48.70638, 6.1550176 47.6202079)     00421102|00264758|LINESTRING (6.06043 48.70638, 6.4474578 48.1758391)     00421102|00269360|LINESTRING (6.06043 48.70638, 6.0215766 47.2392361)     00421102|00269382|LINESTRING (6.06043 48.70638, 5.384374 49.1602425)     00421102|00269399|LINESTRING (6.06043 48.70638, 7.3564671 48.0791666)     00421102|00269478|LINESTRING (6.06043 48.70638, 6.8588071 47.6378244)     00421102|00429683|LINESTRING (6.06043 48.70638, 7.7465095 48.5839322)     00421102|00421102|MULTIPOINT (6.06043 48.70638, 7.34101 47.7495, 6.1812365 48.6876425)


title: use a selection of features in ArcMap in Python script

text: I am trying to write a python script and am stuck on the first obstacle. I want the user to select features in ArcMap and then use these selected features in a python script. Is that at all possible and does someone have a code example? The code that I have now uses all the features instead of just the selected ones. Please see code snippet below.               arcpy.MakeFeatureLayer_management(pipes,"pipeslyr")     arcpy.SetParameter(0,"pipeslyr")     pipesnew = os.path.join(ws,"pipesnew")          arcpy.CopyFeatures_management("pipeslyr", pipesnew)


title: MapServer Mapfile: One layer, disable the label but not the point for first two zoom levels

text: This map is not using OpenLayers, just MapServer with a mapfile. There is a layer that displays both the point and the label. Looking through the documentation, I'm not entirely sure what the best way is. Edit: Ideally (not even sure if this is possible), I would like to disable the label for the first two zoom levels but display the associated point. Is this even possible if they are tied together as one layer? Edit2: Looks like LABELMAXSCALEDENOM is what I'm looking for, but how to configure correctly? I know it's                LABELMAXSCALEDENOM integer      but how do you exactly know the scale of your map? Currently:               MAP       SIZE 128 768         EXTENT -180 -90 180 90      I tried LABELMAXSCALEDENOM equal to some different integers. Didn't notice any difference: the labels were just gone completely, no matter what the zoom. Any input is appreciated.


title: How can I solve 'Handle bad layers' error when using OpenStreetMap in Qgis?

text: I can easy download Portion of maps from OpenStreetMap, save them etc. Once I reopen the same map I get a message that says "Handle bad layers" for all the OSM files. What I am doing wrong?


title: ESRI imagery service that is WMS Compliant

text: I want to use the ESRI free online imagery service http://services.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer But I want a service that is WMS compliant. Is there a version of this service that is WMS compliant?


title: How to overlay geographic data from multiple sources?

text: Being a newbie to GIS, I am trying to create a map using shapefiles and addresses found online, as an exercise. The steps are below:   1. I downloaded a shapefile of Connecticut from the following link (Connecticut State House Districts; CT House Districts; 2010)   2. I picked up a list of hospitals in Connecticut from this link   3. I use the R script from this link to obtain a long/lat for each hospital   4. I followed this discussion here to create a shapefile for hospitals; hospitals are presented as points   5. I used QGIS to import the 2 shapefiles, the Connecticut boundary file, and the hospital location file But I find that the hospitals are not located within Connecticut: the dots are displaced to other empty space, not on the Connecticut map. I wonder if I missed anything when I was creating the shapefile for hospitals. Can anyone help? Thanks. Below is my R-script for creating the hospital shapefile. I guess the `CRS` is the problem, but I'm not sure.               hosp.shp <- SpatialPointsDataFrame(coords=hosp.list[,c("coord_x","coord_y")],data=data.frame(hosp_name=hosp.list[,"hosp_name"]),proj4string=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"))          writeOGR(hosp.shp,dsn="hosp_shp",layer="hosp_shp",driver="ESRI Shapefile")      ### update 01 I am using wgs84 version of Connecticut State House Districts, CT House Districts, 2010. I don't have specific reason to use this shp file, I just see this has quite a few of polygon, which may be easier for me to continue my GIS exercise. Some elaboration on step 2. From the link, I have the names of quite a few of hospital, I used R to access google map to extract the long/lat for each hospital, e.g. "Advanced Center for Rehabilitation Medicine" has the long/lat of 41.1175966,-73.4078968. I used google map to help me do the geocoding, and now I have a table with the name of the hospital, its long and lat. Then this table goes to step 3, using R again to produce a shp file. I import the self-made and the Connecticut shp file into QGIS, both layer can be show, and I can see the hospital seems to cluster to some location, but all of the hospitals are not located within the boundary of Connecticut, they are FAR away. For the hospital shp file, I used `CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")`, which is copied direct from the Connecticut file, hoping that this will fix the "misalignment" problem, but no dice.


title: How to unproject coordinates on a Google map with Openlayers?

text: Below is my code and a pic. When I finish making my circle on a Google map with Openlayers, I get the following coordinates but I want coordinates like E.G: bottom: -27.762124575762 left: 25.828500424238 right: 27.402749575762 top: -26.187875424238 But instead I get: ![enter image description here](http://i.stack.imgur.com/yR0rs.png) Here is my code:                           function init() {                      // Add the map                 var options = {                 projection: new OpenLayers.Projection("EPSG:900913"),                 displayProjection: new OpenLayers.Projection("EPSG:4326"),                 units: "m",                 maxExtent: new OpenLayers.Bounds(-20037508.34, -20037508.34,                                                  20037508.34, 20037508.34)             };             map = new OpenLayers.Map('map', options);                  // Add the OpenStreetMaps Layer             var osm = new OpenLayers.Layer.OSM.Mapnik("OpenStreetMap (Mapnik)", {                 displayOutsideMaxExtent: true,                 wrapDateLine: true,                 buffer: 0             });                  // Add the google maps Layer             var gmap = new OpenLayers.Layer.Google("Google", {sphericalMercator:true});             map.addLayer(gmap);                  // http://bingmapsportal.com/ and use that instead.             var apiKey = "AqTGBsziZHIJYYxgivLBf0hVdrAk9mWO5cQcb8Yux8sW5M8c8opEC2lZqKR1ZZXf";                  // BING Maps             var veroad = new OpenLayers.Layer.Bing({             key: apiKey,             type: "Road",             wrapDateLine: true             });                  map.addLayer(veroad);                  // This is the call to create the OpenLayers Vector Framework Circle Layer and Edit Toolbar             var wmsLayer = new OpenLayers.Layer.WMS( "OpenLayers WMS",                  "http://vmap0.tiles.osgeo.org/wms/vmap0?", {layers: 'basic'});                  // Now we must define our vector style             OpenLayers.Feature.Vector.style['default']['strokeWidth'] = '2';                       // Styles for Default and Select states                 var myStyles = new OpenLayers.StyleMap({                 "default": new OpenLayers.Style({                     fillColor: "#000",                     fillOpacity: 0.3,                     strokeColor: "#000",                     strokeOpacity: 0.3,                     strokeWidth: 1,                     graphicZIndex: 1                 }),                 "select": new OpenLayers.Style({                     fillColor: "#ccc",                     fillOpacity: 0.6,                     strokeColor: "#ccc",                     strokeWidth: 1,                     label: "Label for Mouseover",                     labelAlign: "cc",                     fontColor: "#333333",                     fontOpacity: 0.9,                     fontFamily: "Arial",                     fontSize: 14,                     graphicZIndex: 2                 })             });                  // Add a vector Layer which holds the markers             vectors = new OpenLayers.Layer.Vector("Vector Layer");             map.addLayer(vectors);                  vlayer = new OpenLayers.Layer.Vector( "Editable" , {                 styleMap: myStyles,                 rendererOptions: {zIndexing: true}             } );                     map.addLayer(vlayer);              // Options for the sides wanted in the default circle e.g 3 sides = triangle             polyOptions = {sides: 40};                  // At this point we must make a control to draw a circle on the map after a button to activate is clicked             polygonControl = new OpenLayers.Control.DrawFeature(vlayer,                                             OpenLayers.Handler.RegularPolygon,                                             {handlerOptions: polyOptions} );                  // Now we add the an alert for an option to delete once the user has drawn a circle or vector on the map             vlayer.events.on({             featuresadded: function (event) {             var confirmPolygon = function () { return confirm("Do you want to keep this polygon?") };                  if (!confirmPolygon()) {                 vlayer.removeFeatures(event.features);               }              }             });                  // Now we call an alert to get the bounds or coordinates from a circle or vector we have drawn              vlayer.events.on({             featuresadded: onFeaturesAdded             });                  function onFeaturesAdded(event){             var bounds = event.features[0].geometry.getBounds();             var answer = "bottom: " + bounds.bottom  + "\n";             answer += "left: " + bounds.left  + "\n";             answer += "right: " + bounds.right  + "\n";             answer += "top: " + bounds.top  + "\n";             alert(answer);             }                   }


How can I create polygon of gaps/holes if I have several shapefiles. Do I really have to merge all layers, it it big amount of data (NEW YORK roadbed, buildings, water,parks etc) I have ET geowizards to clean polygons or create gaps but I can only choose one layer. Merging all layers takes so long time. maybe it is another way? erasing, clip with an background layer?. Union with no gaps allowed worked but i need to be able to symbolize the gap polygon and i need to be able to select individual gap polygons. any ideas?


title: How to handle multipoint symbology catalogue (as in mil2525c, tactical graphics appendix b)

text: Is there a general format (standards based) to store a vector based drawing description for multipoint symbology (means having several anchor points reflecting symbol point <-> geographic point association)? OGC defined the Styled Layer Descriptor which uses symbolizers (point/line/polygon) but one problem is that they are better suited to be created for a single point representation (to place e.g. a raster based symbol at specific location, maybe with offset) or the have a stronger dependency to the geometry (following a line with an offset, placing somthing at the line, ...) like a line/polygon symbolizer. If you look at the mil2525c standard there are several symbols where the graphical representation has nothing in common with the stored geometry except the anchor point association to "stretch" the symbol. So I expect those problems using OGC conformant aproach using sld:   1. I don't see any default "multipoint" symbolizer, only a point symbolizer, a line symbolizer and a polygon symbolizer. Line and polygon symbolizers both don't fit because they only styles the line/polygon geometry ( => abstaction to mil2525c multipoint geometries and styling them using sld would be very complex... not to say impossible.)   2. Point symbolizer can call a mark factory (e.g. see GeoTools interface/implementation docs at http://docs.geotools.org/latest/javadocs/org/geotools/renderer/style/MarkFactory.html) to get a shape instead of a raster based symbol but there is no general approach to handle anchor points in a default manner.   3. If using the 2nd. approach using a point symbolizer: The problem of storing vector based data for symbology generation is not solved. ("Hand-made" generation by building the geometry seems to be an option, but is difficult to analyse or edit. What type is shape in the .net world or what would be the "universal" equivalent (System.Windows.Shapes.Path)? Are there other standards for multipoint symbols I didn't noticed? How would you solve this problem?


title: OpenLayers with Google, Bing & MapQuest - Does it violate the terms of service?

text: to my understanding, you are not allowed to use the Google Maps API to display Google and Bing maps (Layer Toggle) (From the terms of service). However, I'm unclear if I can have an OpenLayers map that allows the users to switch between Google, Bing and MapQuest. I want to let the user be able to choose the best Satellite/Aerial data available for their location. I'm really unclear if this is a Terms Of Service violation. Any ideas? Thanks, Nick,


title: How to rotate labels using ArcPy?

text: I'm writing Python code for ArcGIS using the arcpy library and now I'm experiencing a little problem with labels. In fact I'm not able to rotate them in layout via arcpy (I should rotate all labels of 90 degrees). I found a way to display correct field, round the number, but not to position the label correctly.


title: openlayers set wfs layer

text: # error Could not determine geoserver request from http request org.geoserver.platform.AdvancedDispatchFilter$AdvancedDispatchHttpRequest@74f148 # wfs layer initalize                                         //wfs-t editable overlay                     var wfs_layer = new OpenLayers.Layer.Vector("Editable Features", {                     strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],                     //projection: new OpenLayers.Projection("EPSG:2322"),                     protocol: new OpenLayers.Protocol.WFS({                     version: "1.1.0",                     // loading data through localhost url path                     url: "http://192.168.1.236:8080/geoserver/Kaski/wfs",                     featureNS :  "http://192.168.1.236:8080/Kaski",                     maxExtent: mapextent,                     // layer name                     featureType: "PIN",                     // geometry column name                     geometryName: "the_geom",                     schema: "http://192.168.1.236:8080/geoserver/Kaski/wfs?service=WFS&version=1.0.0&request=DescribeFeatureType&typeName=Kaski%3APIN"                     })                     });      it gets data and show but now set new data. # my proxy.ashx.cs               namespace KaskiGIS     {         public class Proxy : System.Web.IHttpHandler         {             public void ProcessRequest(System.Web.HttpContext context)             {                 var request = (System.Net.HttpWebRequest)System.Net.WebRequest.Create(context.Request["url"]);                 request.Method = context.Request.HttpMethod;                      if (request.Method == "POST")                 {                     var content = new System.Text.UTF8Encoding().GetBytes(new System.IO.StreamReader(context.Request.InputStream).ReadToEnd());                          request.ContentLength = content.Length;                     using (var stream = request.GetRequestStream())                     {                         stream.Write(content, 0, content.Length);                         stream.Close();                     }                 }                      var response = request.GetResponse().GetResponseStream();                 if (response == null) return;                      context.Response.Write(new System.IO.StreamReader(response).ReadToEnd());             }                  public bool IsReusable { get { return false; } }         }     }


title: How to make autoIncrement number start with (001....009, 011......099)?

text: i have a challenge with my time in my work, i have feature class called "Parcel" storing many data, i generated an autoIncrement number by using field calculator, in field "parcel code" `Expression: autoIncrement() Expression Type: PYTHON_9.3 Code Block: rec=0 def autoIncrement(): global rec pStart = 1 #adjust start value, if req'd pInterval = 1 #adjust interval value, if req'd if (rec == 0): rec = pStart else: rec = rec + pInterval return rec ` that i asked here before How to add automate values to attribute field? i want to start number to 001......009, and 010.........099,   thank you


title: Merge shapefiles to one tool not working?

text: When using the data management tool "Merge shapefiles to one" in QGIS why is it copying one of my shapefiles instead of merging the two together as it should?


title: How to display labels - created using LabelLayer only at certain scales

text: I am using ESRI Javascript 3.7 framework to develop a map application, in it I am trying to display labels created using LabelLayer at certain scales only, say from 1:1 to 1:500 scale only, how can I achieve this? I played a little with ScaleDependentRenderer but it results in 2 firebug exceptions. Not sure if the ScaleDependentRenderer is the right answer or not. So, if anyone can help, I would really appreciate.               var occupancyLabelLayer = new esri.layers.LabelLayer({ id: "labels" });     // add the label layer to the map     map.addLayer(occupancyLabelLayer);          var labelField = "FLOOR";     // create a text symbol to define the style of labels          var spaceLabelColor = new dojo.Color("red");     var spaceOccupancyLabel = new esri.symbol.TextSymbol().setColor(spaceLabelColor);     //spaceOccupancyLabel.font.setSize("14pt");     spaceOccupancyLabel.font.setFamily("arial");     var spaceOccupancyLabelRenderer = new esri.renderer.SimpleRenderer(spaceOccupancyLabel);          var spaceLabelScaleDependentRenderer = new ScaleDependentRenderer({         rendererInfos: [{             renderer: spaceOccupancyLabelRenderer,             minScale: 500,             maxScale: 1         }]     });          occupancyLabelLayer.addFeatureLayer(featureLayer, spaceLabelScaleDependentRenderer, "${" + labelField + "}");      Firebug Exceptions:               1)exception in animation handler for: onEnd     http://js.arcgis.com/3.7/ (line 45)     2)TypeError: b is undefined     http://js.arcgis.com/3.7/js/esri/layers/LabelLayer.js    Line 21


title: Expert (Recommender) System for Agriculture?

text: I would like to know is there any recommender system for agriculture ?? for example if i input weather and other attributes it will recommend me crop (wheat is better to be considered in these conditions) ??


title: Python OGR: convert a point geometry into polygon

text: I'm writing a routine for creating the convex hull envelop around shapefiles of islands. I could read the data, extract polygons vertices into a geometry and then compute the convex hull. Now I need to save it to a shapefile of polygons. Is there a way to transform the geometry (named `convexHull` in my code) returned by the convexHull method into a polygon geometry; or am I forced to parse all vertices of `convexHull` and create a ring of a polygon?               def doConvexHull(infile, outfile):         inH = ogr.Open(infile, 0)         if inH is None:             usage("Could not open file {0}. Exit.".format(infile))         layer = inH.GetLayer()              # get all polygons         thisGeometry = ogr.Geometry(ogr.wkbPoint)         for index in xrange(layer.GetFeatureCount()):             feature = layer.GetFeature(index)             geometry = feature.GetGeometryRef()             ring = geometry.GetGeometryRef(0)             points = ring.GetPointCount()             for p in xrange(points):                 lon, lat, z = ring.GetPoint(p)                 thisGeometry.AddPoint(lon, lat)              convexHull = thisGeometry.ConvexHull()              drv = ogr.GetDriverByName( "ESRI Shapefile" )         ds = drv.CreateDataSource( outfile )         if ds is None:             usage("Could not create file {0}".format( outfile) )              lyrname = "convexHull_${0}".format( layer.GetName() )         lyr = ds.CreateLayer( lyrname, layer.GetSpatialRef(), ogr.wkbPolygon )         thisFeature = ogr.Feature( layer.GetLayerDefn() )         thisFeature.SetGeometry( convexHull )         lyr.CreateFeature( thisFeature ): error from geometry


title: How to draw polygons on an existing map?

text: If you want a map with some circles, squares, etc. representing cities or something else you can use the following code. Note that you can change the map and the coordinates with yours, and you'll get your own map in a java application!               import java.awt.Color;     import java.net.URL;     import org.geotools.data.DataUtilities;     import org.geotools.data.shapefile.ShapefileDataStore;     import org.geotools.data.simple.SimpleFeatureCollection;     import org.geotools.data.simple.SimpleFeatureSource;     import org.geotools.factory.CommonFactoryFinder;     import org.geotools.feature.FeatureCollections;     import org.geotools.feature.simple.SimpleFeatureBuilder;     import org.geotools.geometry.jts.JTSFactoryFinder;     import org.geotools.map.FeatureLayer;     import org.geotools.map.MapContent;     import org.geotools.renderer.lite.StreamingRenderer;     import org.geotools.styling.FeatureTypeStyle;     import org.geotools.styling.Graphic;     import org.geotools.styling.Mark;     import org.geotools.styling.PointSymbolizer;     import org.geotools.styling.PolygonSymbolizer;     import org.geotools.styling.Rule;     import org.geotools.styling.Style;     import org.geotools.styling.StyleBuilder;     import org.geotools.styling.StyleFactory;     import org.geotools.swing.JMapFrame;     import org.opengis.feature.simple.SimpleFeature;     import org.opengis.feature.simple.SimpleFeatureType;     import org.opengis.feature.type.GeometryDescriptor;     import org.opengis.filter.FilterFactory;     import com.vividsolutions.jts.geom.Coordinate;     import com.vividsolutions.jts.geom.GeometryFactory;          public class MyGIS{         public static void main(String[] args) {             try {                 usingFeatureCaching();             } catch (Exception e) {                 // TODO Auto-generated catch block                 e.printStackTrace();             }         }              public static void usingFeatureCaching() throws Exception {                  //Shape-file's path             URL url = new URL("file://C:/PE/Italia/reg2011_g.shp");             ShapefileDataStore shapefile = new ShapefileDataStore(url);                  // Creates a map and adds the shapefile             MapContent map = new MapContent();                  //Set's windows title             map.setTitle("Italy");                  //Creates the map style             StyleBuilder styleBuilder = new StyleBuilder();             PolygonSymbolizer restrictedSymb = styleBuilder.createPolygonSymbolizer(Color.LIGHT_GRAY, Color.BLACK, 0);                  //Sets opacity             restrictedSymb.getFill().setOpacity(styleBuilder.literalExpression(0.5));             org.geotools.styling.Style myStyle = styleBuilder.createStyle(restrictedSymb);                  //Adds another layer to the map                  FeatureLayer layer = new FeatureLayer(shapefile.getFeatureSource(), myStyle);             map.addLayer(layer);                  //-------------------------- BUILDS THE CIRCLE ON THE MAP ------------------------------------------//                  StyleFactory sf = CommonFactoryFinder.getStyleFactory(null);             FilterFactory ff = CommonFactoryFinder.getFilterFactory2();                  SimpleFeatureSource fs = shapefile.getFeatureSource();                  SimpleFeatureType pointtype = DataUtilities.createType("Location", "the_geom:Point," + "name:String");                  SimpleFeatureBuilder sfb = new SimpleFeatureBuilder(pointtype);                  double longitude = 537319.7867d;             double latitude = 5062350.2792d;                  GeometryFactory geometryFactory= JTSFactoryFinder.getGeometryFactory(null);             com.vividsolutions.jts.geom.Point point = geometryFactory.createPoint(new Coordinate(longitude, latitude));             sfb.add(point);                  SimpleFeatureCollection col = FeatureCollections.newCollection();             SimpleFeature feature1 = sfb.buildFeature(null);             col.add(feature1);                       org.geotools.styling.Stroke stroke2 = sf.createStroke(                     ff.literal(new Color(0xC8, 0x46, 0x63)),                     //circle thickness                     ff.literal(1)             );                  org.geotools.styling.Fill fill2 = sf.createFill(                     ff.literal(new Color(0xff, 0xC8, 0x61)));                  map.setTitle("Italia");                       //Type of symbol             Mark mark = sf.getCircleMark();                  mark.setFill(fill2);             mark.setStroke(stroke2);                  Graphic graphic = sf.createDefaultGraphic();             graphic.graphicalSymbols().clear();             graphic.graphicalSymbols().add(mark);                  //circle dimension on the map             graphic.setSize(ff.literal(5));                  GeometryDescriptor geomDesc = fs.getSchema().getGeometryDescriptor();             String geometryAttributeName = geomDesc.getLocalName();             PointSymbolizer sym2 = sf.createPointSymbolizer(graphic, geometryAttributeName);                  Rule rule2 = sf.createRule();             rule2.symbolizers().add(sym2);             Rule rules2[] = {rule2};             FeatureTypeStyle fts2 = sf.createFeatureTypeStyle(rules2);             Style style2 = sf.createStyle();             style2.featureTypeStyles().add(fts2);                  map.addLayer(new FeatureLayer(col, style2));                       StreamingRenderer renderer = new StreamingRenderer();             renderer.setMapContent(map);                  //Shows the map              JMapFrame.showMap(map);         }     }      This is the result of the code: ![The code's result](http://i.stack.imgur.com/ULYkF.png) PS: Here you can find the complete project, enjoy!


title: Applying depth to Shapefile

text: Apologies if this question has been answered before. I'm digitizing a historic coastline and would like to add a single depth figure to it in order to merge it with a current contour model. Is there and easy to add depth to the polyline shapefile I am creating?


title: Get X, Y coordinates in the center of a raster cell

text: I am trying to get the X, Y coordinates of all the cells of a raster file, using Python 2.6 and ArcGIS 10.0. I would like to get the X, Y coordinates in the center of the cell. Is there any way to do so?


title: How to Calculate North?

text: I am working for field survey & client requirement is differentiate various north,like magnetic north,true north & grid north.So therefore my question is what is the methodology of calculating north & how can we differentiated various type of north as above mentioned.


title: Common class range for different attributes in a same polygon shapefile

text: I have a polygon shapefile which has population attributes pop90,pop91 and so on for 10 years. They have different minimum and maximum value for each of the years. Now I need to do classification for each attribute(pop90,pop91) but with the same number of classes and the range for each class should be same. Same number of classes is fine but I am not sure on how can I use the same class range for each of the attribute. Any idea how can this be done. Thank You.


title: Get values from select by location to append to table

text: I'm trying to calculate the road frontage of polygons by converting polygons to lines and splitting in their vertices, then using select by location to find lines that are within 7.5 meters of the road. I then want to retrieve the road label and append it to the lines table. This is what I have, which find the lines correctly and writes the road information to the lines table, but it doesn't get the correct road info. It doesn't seem to find the road that was used in the select by location tool. Can anyone help?               arcpy.PolygonToLine_management("BLPUs", "BLPUs_PolyLine")     arcpy.SplitLine_management("BLPUs_PolyLine", "BLPUsPolyLineSplit")     arcpy.MakeFeatureLayer_management("BLPUsPolyLineSplit", "BLPU_PolyLineSplit")     arcpy.SelectLayerByAttribute_management("BLPU_PolyLineSplit", "NEW_SELECTION", '"LEFT_FID" <> -1 ')     arcpy.DeleteRows_management("BLPU_PolyLineSplit")     in_layer = "BLPU_PolyLineSplit"     out_layer = "BLPUs_PolyLineSplit"     arcpy.CopyFeatures_management(in_layer, out_layer)     arcpy.AddField_management("BLPUs_PolyLineSplit", "RoadLabel", "TEXT")     arcpy.AddField_management("BLPUs_PolyLineSplit", "DescTerm", "TEXT")     arcpy.AddField_management("BLPUs_PolyLineSplit", "Nature", "TEXT")     arcpy.MakeFeatureLayer_management("BLPUs_PolyLineSplit", "BLPU_PolyLineSplitt")     with arcpy.da.UpdateCursor("BLPU_PolyLineSplitt", ["RoadLabel", "DescTerm", "Nature"]) as cursor:         field1 = "RoadLabel"         field2 = "DescTerm"         field3 = "Nature"         for row in cursor:             arcpy.SelectLayerByLocation_management("BLPU_PolyLineSplitt", "WITHIN_A_DISTANCE", "C:\Users\Ant\Documents\ArcGIS\OSMaster.gdb\ITN_RoadLinks", "7.5 Meters")             TABLE = "ITN_RoadLinks_Copy"             searchRows = arcpy.da.SearchCursor(TABLE,["RoadLabel", "DescTerm", "Nature"])             for searchRow in searchRows:                 searchRow1 = str(searchRow[0])                 searchRow2 = str(searchRow[1])                 searchRow3 = str(searchRow[2])             row[0] = searchRow1             row[1] = searchRow2             row[2] = searchRow3             cursor.updateRow(row)     arcpy.CopyFeatures_management("BLPU_PolyLineSplitt", "BLPUs_FrontageLines")


title: Can I import a .gml file into Spatialite in QGIS?

text: I have some Ordnance Survey Address Base data files which have a .gml extension. I have imported them to GGIS as vector layers along with some mapping data but would like to explore them with PostgreSQL. I see that QGIS comes with SpatiaLite installed. Can I import the gml files into SpatiaLite or for that matter PostgreSQL? I suspect the format and content of my question says more about what I don't know so any comments or even answers would be appreciated.


title: Why do text boxes and arrows shift position when exported from print composer to jpeg QGIS 2.2

text: I have created a map using a 1:20000 WMS layer. On top of this I created a shapefile layer tracing a mapped stream. In print composer I added a couple of text boxes with arrows pointing to certain sections of the stream. I then exported this as a jpeg to insert in to a report in MS Word. Unfortunately, the exported version shows the arrows and text boxes displaced from the position in print composer. I have checked my composition tab in composer and this matches the page dimensions used in Word. Does anyone have an idea why this offset is happening?


title: Geoserver Postgres Configuration

text: My imported postgres database size for plane.osm is more than 300GB under slim mode. It added four tables for Tiles Creation. _planet_osm_line, planet_osm_road, planet_osm_polygon, planet_osm_point_ Each Table contains attributes which are: > osm_id   >  access   >  addr:housename   >  addr:housenumber   >  addr:interpolation admin_level   >  aerialway   >  aeroway   >  amenity   >  area   >  barrier   >  bicycle   >  brand   >  bridge   >  boundary   >  building   >  construction   >  covered   >  culvert   >  cutting   >  denomination   >  disused   >  embankment   >  foot   >  generator:source   >  harbour   >  highway   >  historic   >  horse   >  intermittent   >  junction   >  landuse   >  layer   >  leisure   >  lock   >  man_made   >  military   >  motorcar   >  name   >  natural   >  office   >  oneway   >  operator   >  place   >  population   >  power   >  power_source   >  public_transport   >  railway   >  ref   >  religion   >  route   >  service   >  shop   >  sport   >  surface   >  toll   >  tourism   >  tower:type   >  tracktype   >  tunnel   >  water   >  waterway   >  wetland   >  width   >  wood   >  z_order   >  way_area   >  way Now I want to configure it to Geoserver for a specific attribute. Example: in polygon, I need to import `building` alone. I am following this GeoServer and OpenStreetMap and Publishing a PostGIS Table but it's not specificly clear about the above. I need to know if it's possible in Geoserver? And how ?


title: Adding/Removing Leaflet GeoJSON layers

text: I'm trying to show different GeoJSON layers at different zoom layers using the Leaflet API. I can load and display all three layers at once (though I don't actually want them all to show at once). I can load and display them at different zoom levels. I have the code set up so that at Zoom levels 1-6, the map will show one GeoJSON layer. At levels 7-10, it will show another, and at levels 11+ it will show a third. Displaying them works. What I'm trying to get working now is turning the others off if one is displayed. Going from 1-6 to 7-10 works (meaning that it turns off the 1-6 layer correctly), but not from 7-10 to 11+ (meaning that the 7-10 layer sticks around) and I can't figure out why (it uses the same code). Here's the ajax for the GeoJSON layers:               function getJson(defaultStyle, map, simp, geojsonLayer){     var url = 'file' + simp + '.json';     map.removeLayer(geojsonLayer);     geojsonLayer.clearLayers();     $.getJSON(url, function(data){         geojsonLayer = L.geoJson(data, {             style: defaultStyle,             onEachFeature: onEachFeature         });         geojsonLayer.addTo(map);     });     }      And here's the main function that calls the ajax depending on the zoom. simpCounter is set to 0 initially.               map.on('zoomend', function(e) {     if (map.getZoom() >= 7 && map.getZoom() <= 10) {         if (simpCounter == 0 || simpCounter == 2) {         getJson(defaultStyle, map, 60, geojsonLayer);         simpCounter = 1;         }     } else if (map.getZoom() >= 11) {         if (simpCounter == 0 || simpCounter == 1) {         getJson(defaultStyle, map, 35, geojsonLayer);         simpCounter = 2;         }     }     });      So again, the first transition turns off the old layer correctly, the but second transition does not. Thanks for the help.


title: How to make heat maps?

text: It might beginner question regarding heat maps..I want to prepare the heat maps so what parameters are required for that?? I gone through following posts about heat map..but I did not get exact idea.. I have large point data (I want to represent the density on the map) I gone through following posts : How to build effective heat-maps? http://forums.arcgis.com/threads/43069-How-to-create-a-heat-map-from- shapefile-data still unclear..any help will be great..thanks in advance : ) software's ( ArcGIS Desktop 10.1 and ArcGIS Server 10.1 + ArcGIS Java Script API 3.0)


title: How to merge SeaWiFS and MODIS data or use other merged products?

text: I hope to produce a time-series for sea surface temperature (SST) and chlorophyll-a (Chl-a) for use in MaxEnt without the use of arcgis, preferring Python tools to others. I am new to satellite data processing, and there appears to be quite a lot to consider. What must be taken into account for merging SeaWiFS and MODIS datasets (e.g. discrepancies between the data, projections, etc.)? Are there other products available that provide such merged data that are generally recommended? I have found the following:   * NOAA - CoRTAD (SST)   * ESA - MEDSPIRATION (SST)   * ESA - GlobColour (SST & Chl-a) Ultimately I need all of the data gridded to the same grid size/area, so I would to write a script where I could acquire the data, re-project if necessary, re-sample to the grid if necessary, and then save to desired format (ESRI's ASCII .asc) for running in the model.


title: How to extract POINT from a [LINE|MULTILINE]STRING

text: I'd like to create a vectorLayer in `QGIS` to display all the points that form a `LINESTRING` or a `MULTILINESTRING` stored in a `PostGIS` DB. I think that I need to transform all points of `LINESTRING` (or `MULTILINESTRING`) to `POINT`. **First Question** : Is there any `PostGIS` function to do that? **Second Question** : If I want to exclude duplicate points from the result `POINT` table, how can I do that?


title: Keeping attribute data when exporting KML to feature class

text: I am working on updating kml files into an mxd. I am able to use the conversion tool to change each KML to a layer file, but I am unable to keep the original kml attribute data. When I export point files, it exports as one single point instead of multiple points which does not make sense to me. Is there a way convert KML/KMZs to a feature class and keep all of the original data?


title: GRID TO TXT in saga

text: i have a problem in SAGA. I want to export several raster files to a .txt file. It means integration raster files to a .txt file. Can I do this in SAGA. Please explain? Thanks


title: Create hachures from DEM in ArcGIS 10.1

text: I have created contour lines from DEM and now I need to create hachures/bergstrich for them. I believe that spatial analyst should solve this, but I'm not used to it and I didint found any similar tutorials. Any solutions will help. thanks in advance.


title: How do I include proj4js in openlayers?

text: in my Openlayers app, since the base map projection is not epsg900913 or 4326, I need to include proj4js in order to use geolocation feature. but i got such a error notice:                GET http://svn.osgeo.org/metacrs/proj4js/trunk/lib/projCode/null.js 404 (Not Found)      i am wondering how come this null.js is needed here. the funny thing is i can't find this file from proj4js lib files.. anybody can help me out? thanks!! Code:               <!DOCTYPE html>                      <link rel="stylesheet" href="theme/default/style.css" type="text/css">      <link rel="stylesheet" href="style.css" type="text/css">


title: How to use Google Maps geocoding with Openlayers?

text: I'm making a Openlayers map with some shapefile and I have Google Maps base layer. I want to create a text area and a button. When I input a address and click the button, the map should show me the google maps address in the map. Can you help me please ?


title: Drawing pie charts with vectors on OpenLayers?

text: Is anyone got any examples of drawing pie charts with vectors as markers on OpenLayers ? Cheers,


title: leaflet layerGroup is rendered behind baselayer

text: I have several baselayers and overlays. Everything is fine except the lauerGroup is rendered behind the baselayer. First I define the baselayers and add them to a 'baselayer' variable.               // base layers     naturalEarth = L.tileLayer('http://{s}.tiles.mapbox.com/v3/mapbox.natural-earth-2/{z}/{x}/{y}.png'),          grump = L.tileLayer('http://localhost:8000/v2/gdal_urban_rural_pop/{z}/{x}/{y}.png');          var baselayer = {         'Natural Earth' : naturalEarth,         'GRUMP Population' : grump     };      Then I do the same for some overlays:               // overlay layers     osm = L.tileLayer('http://{s}.tiles.mapbox.com/v3/mapbox.world-bank-borders-en/{z}/{x}/{y}.png'),          df_admin0 = L.tileLayer('http://localhost:8000/v2/dwelling_fractions/{z}/{x}/{y}.png', {             wax: 'http://localhost:8000/v2/dwelling_fractions.json'         }),          df_port = L.tileLayer('http://localhost:8000/v2/PRT_dwelling_fractions/{z}/{x}/{y}.png', {             wax: 'http://localhost:8000/v2/PRT_dwelling_fractions.json'         });          //add layer group     hazard_curve_group = L.layerGroup();          hazard_curve1 = L.tileLayer('http://localhost:8000/v2/hazard_curve/{z}/{x}/{y}.png', {             wax: 'http://localhost:8000/v2/hazard_curve.json'         }).addTo(hazard_curve_group);          hazard_curve2 = L.tileLayer('http://localhost:8000/v2/hazard_curve2/{z}/{x}/{y}.png', {             wax: 'http://localhost:8000/v2/hazard_curve.json2'         }).addTo(hazard_curve_group);          var overlays = {         "OpenStreetMap" : osm,         "Dwelling Fractions Admin 0" : df_admin0,         "Dwelling Fractions Portugal" : df_port,         "Hazard Curves test" : hazard_curve_group,     };          // create the map     var map = L.map('map', {         ...     });          //add it all to the map     L.control.layers(baselayer, overlays, {collapsed: true, position: 'topleft'}).addTo(map);


title: How to import 3D geometries from DXF, VRML or X3D into PostGIS 2

text: we plan to develop a webbased 3D GIS using PostGIS 2 and GeoServer with a W3DS. The architecture works fine, but i have problems to import 3D geometries into PostGIS. The geometries come from AutoCAD and consist of 3DFACEs and 3DSOLIDs. I tried to import it with ogr2ogr, but I think it does not support these geometries. The idea is to convert all this 3D objects into 3D-POLYGONs or Polyhedral Surfaces. Does anyone have an idea? Thanks Frank


title: Managing a large number of layers in OpenLayers/GeoExt and Geoserver

text: I am building a web application to facilitate the handling of approximately 2500+ historic maps. I think it would be too difficult for our users to search for specific maps in an ordinary layer tree. Is there some way to implement a layersearch in a GeoExt/Geoserver stack? I'd be grateful for any idea how to deal with that even if it's not an "out- of-the-box"-solution.


title: Drop Marker with 3857 projection in leafletjs

text: Is it possible in leaflet to drop marker with 3857 projection on map? I fond a solution where first translate point to 4326 using proj4js and than drop it on map.               var source = new Proj4js.Proj('EPSG:3857');         var dest = new Proj4js.Proj('EPSG:4326');         var p = new Proj4js.Point(-12307.601047734, 6711475.1283642);         Proj4js.transform(source, dest, p);      new L.Marker([p.y, p.x], {bounceOnAdd: true}).addTo(map); is there any other way in core leaflet or leaflet plugin to drop marker without using proj4js library. Thanks.


title: Convert FeatureLayerClass to shape file error : Item not found in collection

text: so following this link http://edndoc.esri.com/arcobjects/9.2/NET/91e3e8a0-d355-43d2-a201-0a8ceb719334.htm I am able to convert my featureLayer into a shape file and export it. However, some how for 1 specific layer that i tested with i always get this error: **Item NOT found in collection** during this piece of code:               IFeatureDataConverter featureDataConverter = new FeatureDataConverterClass();     IEnumInvalidObject enumInvalidObject =       featureDataConverter.ConvertFeatureClass(((IDataset)featureLayer.FeatureClass).FullName as IFeatureClassName,       null, null, targetFeatureClassName, targetGeometryDef, targetFields, "", 1000, 0);          // Check for errors.     IInvalidObjectInfo invalidObjectInfo = null;     enumInvalidObject.Reset();     while ((invalidObjectInfo = enumInvalidObject.Next()) != null)     {         // Handle the errors in a way appropriate to the application.         Debug.WriteLine("Errors occurred for the following feature: {0}, {1}",           invalidObjectInfo.InvalidObjectID, invalidObjectInfo.ErrorDescription);         result = false;     }      It does not crash, but the enumInvalidObject always gives me the error Item not found in this collection. Does anyone have an idea on what may cause this, I know saying that ONE of my layer don't work is very vague, but i really dont know what is special about that layer. It does export the shape file, but no features inside. Thanks and Regards, Kev84


title: How to fix "Access to database denied" when accessing PostgreSQL?

text: I got a tiny problem regarding connection to my server. What should we do when this message appears in front of you? ![enter image description here](http://i.stack.imgur.com/bCeoG.jpg) I've done several changes in the ***pg_hba.conf*** entry but yet to succeed. Any careless mistake from my side? Thanks in advance. I love this sentence.


title: Qualifying field names in field calculator

text: I am trying to set up a quick model that:   1. Takes a point feature layer of development projets (a table in an access database with xy coordinates),   2. Spatially joins polygon neighborhood designations to the feature layer,   3. Join the output of 2 to the source feature layer (sitting as a table in access)   4. Updates the neighborhood field based on the neighborhood field from the attributes from the spatial join The point is that I really don't want to create a new file every time I do this, but maintain and add to the existing table of development projects. My question is how to properly qualify field names in the field calculator so I can update values in the source table (1) based on the values of the spatial join featureclass with the neighborhood designation (2). I join the output of (2) to (1) so the attributes should all be in the same name space, I imagine. I currently use this as my update expression (VB parser): `[PipeRealCurrent.PlanningDistrictsCombo]=[PipeRealCurrent_x_PDist.PlanningDistrictsCombo]`               A field name was not found or there were unbalanced quotation marks.     ERROR 999999: Error executing function.     Item not found in this collection.     Failed to execute (Calculate Field).      ![enter image description here](http://i.stack.imgur.com/NGJHW.png)


title: Length of a line within a surface

text: I need to find out the lengt of each line within a surface. I have approx. 20'000 diferent lines and around 10'000 different surfaces. what kind of SQL- sentences do I need to use? I've tried things like:               select name.line, ST_length(line.thegeom, surface.thegeom)     from x.line, y.surface     where st_intersects (line.thegeom, surface.thegeom)      but my progragm (QGIS) just dies all the time or says something like "geos intersects () threw an error" I'd love to have a sql query that gives me following GID(=number of the line), GIS_NR(=number of the surface), total length of the line within the surface. I'm quite new to GIS and postgresql I'd apreciate your help.


title: Convert DEM to STL (3D Printer Format)

text: I have a DEM in ArcGIS 10.1 that I would like to get into STL format for a 3D printer. Does anyone know some software that will allow me to make this conversion? I have to use the STL format because that is the format required by my university's 3D printer.


title: What is the approximate error of the Pythagorean Theorem vs. Haversine Formula in measuring distances on a sphere at various scales?

text: Many people when first trying to calculate distances between two longitude / latitude pairs ask if Pythagorean theorem works as an appropriate distance function. Most often people answer "no, the Pythagorean theorem only works on a 2D Euclidean plane." Rarely, however, do people mention the effect of scale and location on the sphere on how inaccurate the Pythagorean theorem is. The basic idea being at very small scales, the surface of a sphere looks very much like a plane. At very large scales, it distances along the surface are more curved and therefore the difference between the incorrect Pythagorean Theorem and the correct Haversine Formula is greater. Does anyone know a formula or rule of thumb that tells you the difference between the two distance measures based on the scale of the distance you are trying to measure? I think having this explicitly would help in [1] explaining why the Pythagorean Theorem isn't perfect and [2] in letting people who are looking for more "rough" distances know when Pythagoras actually will serve their purposes.


title: I can't export in QGIS a Shape file to PostGis

text: When I try export a pipe network shape to postgis i get this error:               Problema inserindo feições do arquivo:     C:/PP/gvSIG_Projects/Shapes/Abastecimento/QGIS_VATROCOADUTOR.shp          Ocorreu um erro no banco de dados ao executar este SQL:     SELECT AddGeometryColumn('public','QGIS_VATROCOADUTOR','the_geom',-1,'NULL',2)     O erro foi:     ERROR:  Invalid type name "NULL(2)" - valid ones are:         POINT, MULTIPOINT,         LINESTRING, MULTILINESTRING,         POLYGON, MULTIPOLYGON,         CIRCULARSTRING, COMPOUNDCURVE, MULTICURVE,         CURVEPOLYGON, MULTISURFACE,         GEOMETRY, GEOMETRYCOLLECTION,         POINTM, MULTIPOINTM,         LINESTRINGM, MULTILINESTRINGM,         POLYGONM, MULTIPOLYGONM,         CIRCULARSTRINGM, COMPOUNDCURVEM, MULTICURVEM         CURVEPOLYGONM, MULTISURFACEM, TRIANGLE, TRIANGLEM,         POLYHEDRALSURFACE, POLYHEDRALSURFACEM, TIN, TINM         or GEOMETRYCOLLECTIONM     CONTEXT:  SQL statement "SELECT AddGeometryColumn('',$1,$2,$3,$4,$5,$6,$7)"     PL/pgSQL function addgeometrycolumn(character varying,character varying,character varying,integer,character varying,integer,boolean) line 5 at SQL statement      I can't interpret this error can you help me, thanks


title: Store Graphics in Local when browser refresh

text: Using arcgis **javascript api** i am **adding** different **graphics** to the **map**. Later i will save them in server. If i forgot to save or due to some reasons **page** is **refreshed** before saving i am **losing all the graphics**. My requirement is replace the graphics to that particular session. Is there any possible way to save the graphics...?? Thanks in advance.


title: How to make digitization more detailed afterwards with QGIS?

text: I digitized already two areas, bordering each other on one side. Than I realized, more detail was needed! What do you suggest as a fast way to do so?   * Adding each time new dots/verteces takes so much time.   * But, adding a new polygon and merge both (as I used to do for a not bordering polygon) does not seem that easy in this case.   * Using the split feature causes an invalid area which I cannot merge with the other area afterwards. Thanks!


title: layer from sqlserver database table

text: I'm new in ArcGIS. Is it possible to create a layer from a SqlServer table instead of a shapefile? in other words, is it possible to store layer data in a table in Sql Sever instead of a shapfile? because I want to make a crude application with these data.. if yes.. please show me how or direct me to any tutorial, thanks for any help


title: Overlay lat/long graticules on projected map QGIS

text: I am using an equal areas projection of North America and would like to overlay a lat/long grid. I can do this in the print composer, but only for the coordinates. Can QGIS do a distorted grid for lat/long coordinates over a projected map? If so, how? Thank you very much


title: Mapping points from total station data in Angle, Horizontal Distance, Vertical Distance

text: I am working on creating a contour map from total station topo points in ArcGIS 10. The points were collected using H angle, HD, VD instead of XYZ points. I need to create points from these measurements with elevations. I have tried using the Bearing Distance to Line Tool and then creating points from the vertices of lines (see http://gis.stackexchange.com/a/31956/18682) but the Bearing Distance To Line Tool does not recognize the angles in my input file (I've tried .csv, excel, and attribute table from feature class created from the excel file). My angles are in Decimal Degrees, which is the default for the tool. What might I be doing wrong? Does anyone have any suggestions for a better method? [EDIT] Also, I am trying to use COGO to create these points. Can anyone help by explaining how to use COGO in ArcGIS Desktop 10 with a table of Distance and Bearings? Thank you very much!


title: Control view of Multi polygon in Google map

text: We have multipolygon WKT. We have no issue to draw all of them in Google map. Below are the codes. The only issue is that we want only the first polygon to be visible but the rest to transparent. How to control on that? Any idea please?               function drawPoly(data, pl) {         var pointsData = data.split(",");         var cords = [];              // iterate over each points data and create a latlong         // & add it to the cords array         var len = pointsData.length;         for (var i = 0; i < len; i++) {             var xy = pointsData[i].split(" ");                  var pt = new google.maps.LatLng(xy[1], xy[0]);             cords.push(pt);         }              // now add it to the map         // Construct the polygon         alert("I=" + pl);         var polygonFirst = "";         if (pl == 0) {             var polygon1 = new google.maps.Polygon({                 paths : cords,                 strokeColor : '#1E90FF',                 strokeOpacity : 0.8,                 strokeWeight : 2,                 fillColor : '#1E90FF',                 fillOpacity : 0.35,                 zIndex : 1             });             // polygon1.setMap(map);         } else {             var polygon2 = new google.maps.Polygon({                 paths : cords,                 strokeColor : '#1122AA',                 strokeOpacity : 1,                 strokeWeight : 2,                 zIndex : 1             });             // polygon2.setZIndex(1);             polygon2.setMap(map);         }          }


title: OpenLayers Circle not visible

text: I'm trying to draw a circle with OpenLayers with a Google Maps Layer. I'm able to transform a simple marker and have it draw at the correct location, but I can't figure out the circle. I'm trying to draw a 50 meter radius circle:               function InitializeMap() {          // create map     map = new OpenLayers.Map({         div: "map_canvas",         theme: null,         units:'m',          projection: "EPSG:900913",         displayProjection: new OpenLayers.Projection("EPSG:4326"),         controls: [             new OpenLayers.Control.Attribution(),             new OpenLayers.Control.TouchNavigation({                 dragPanOptions: {                     enableKinetic: true                 }             })         ]     });          // load layers     var gmap = new OpenLayers.Layer.Google(         "Road", // the default         {numZoomLevels: 20}     );     map.addLayer(gmap);     if (map.layers.length > 0) {         var gphy = new OpenLayers.Layer.Google(             "Terrain",             {type: google.maps.MapTypeId.TERRAIN}         );         map.addLayer(gphy);         var ghyb = new OpenLayers.Layer.Google(             "Hybrid",             {type: google.maps.MapTypeId.HYBRID, numZoomLevels: 20}         );         map.addLayer(ghyb);         var gsat = new OpenLayers.Layer.Google(             "Satellite",             {type: google.maps.MapTypeId.SATELLITE, numZoomLevels: 22}         );         map.addLayer(gsat);         map.addControl(new OpenLayers.Control.LayerSwitcher);     }     var blankmap = new OpenLayers.Layer("Blank", {         isBaseLayer: true,         maxExtent: new OpenLayers.Bounds(             -128 * 156543.03390625,             -128 * 156543.03390625,             128 * 156543.03390625,             128 * 156543.03390625         ),         sphericalMercator: true,         maxResolution: 156543.03390625,         units: "m",         projection: "EPSG:900913"});     map.addLayer(blankmap);          // Set default layer if not offline     if (map.layers.length > 1) {       switch (default_map_type) {         case 1: map.setBaseLayer(blankmap); break;           case 2: map.setBaseLayer(gmap); break;         case 3: map.setBaseLayer(gsat); break;         case 4: map.setBaseLayer(ghyb); break;         case 5: map.setBaseLayer(gphy); break;        }     }          // Add Layers to map         var currentLocationLayer = new OpenLayers.Layer.Markers("Current Location");         var currentAccuracyLayer = new OpenLayers.Layer.Vector("Current Accuracy");         map.addLayer(currentLocationLayer);     map.addLayer(currentAccuracyLayer);              // Add marker         var current_point = new OpenLayers.LonLat(-81,37.transform(map.displayProjection,map.getProjectionObject());         map.setCenter(current_point);         map.zoomTo(17);         current_loc_marker = new OpenLayers.Marker(current_point,image_GPS);         currentLocationLayer.addMarker(current_loc_marker);              // Add circle         var hA_circle_geometry = new OpenLayers.Geometry.Polygon.createRegularPolygon(current_point, 50, 50, 0).transform(map.displayProjection,map.getProjectionObject());         hA_circle = new OpenLayers.Feature.Vector(hA_circle_geometry);     currentAccuracyLayer.addFeatures([hA_circle]);     }      I tried the following for the circle and it still does nto show. Does it need to be transformed due to the Google Layer?:               var hA_circle = new OpenLayers.Control.DrawFeature(currentAccuracyLayer,                                             OpenLayers.Handler.RegularPolygon,                                             {sides: 40, radius: 50, angle: 0});     map.addControl(hA_circle);


title: Why does gdalwarp return incorrect corner coordinates for my geotiff?

text: I have a png source image that's in Lambert Azimuthal Equal Area and I have the long/lat extent. So I ran gdal_translate to create tie points and convert to gtif. Then I ran gdalwarp to reproject to Web Mercator. gdalinfo shows that the bounding coordinates are ok for the gdal_translate output, but not for the gdalwarp output. Que pasa? Here's my gdal_translate operation on the source file:               gdal_translate -of GTiff -a_srs "+proj=laea +lat_0=40. +lon_0=-105 +x_0=0 +y_0=0 +a=6370997 +rf=298.257" ^     -gcp 0 0 -119.344674 47.148697 ^     -gcp 525 0 -90.655326 47.148697 ^     -gcp 525 431 -93.609897 31.438244 ^     -gcp 0 431 -116.390103 31.438244 ^     -a_ullr -119.344674 47.148697 -93.609897 31.438244 ^     src.png tmp.gtif      (I included the bounding box long/lats for good measure but it seems to make no difference in the outcome of gdalwarp.) Here's my gdalwarp operation on the gdal_translate output (tmp.gtif):               gdalwarp -of GTiff -s_srs "+proj=laea +lat_0=40. +lon_0=-105 +x_0=0 +y_0=0 +a=6370997 +rf=298.257" -t_srs EPSG:3857 tmp.gtif final.gtif      Here's the metadata from 'gdalinfo tmp.gtif':               Driver: GTiff/GeoTIFF     Files: tmp.gtif     Size is 525, 431     Coordinate System is `'     Origin = (-119.344674000000000,47.148696999999999)     Pixel Size = (0.049018622857143,-0.036451167053364)     Image Structure Metadata:       INTERLEAVE=PIXEL     **Corner Coordinates:     Upper Left  (-119.3446740,  47.1486970)      Lower Left  (-119.3446740,  31.4382440)      Upper Right ( -93.6098970,  47.1486970)      Lower Right ( -93.6098970,  31.4382440)      Center      (-106.4772855,  39.2934705)**      Band 1 Block=525x3 Type=Byte, ColorInterp=Red       Mask Flags: PER_DATASET ALPHA      Band 2 Block=525x3 Type=Byte, ColorInterp=Green       Mask Flags: PER_DATASET ALPHA      Band 3 Block=525x3 Type=Byte, ColorInterp=Blue       Mask Flags: PER_DATASET ALPHA      Band 4 Block=525x3 Type=Byte, ColorInterp=Alpha      Corner coordinates look good. But here's the metadata from 'gdalinfo final.gtif':               Driver: GTiff/GeoTIFF     Files: final.gtif     Size is 579, 355     Coordinate System is:     PROJCS["WGS 84 / Pseudo-Mercator",         GEOGCS["WGS 84",             DATUM["WGS_1984",                 SPHEROID["WGS 84",6378137,298.257223563,                     AUTHORITY["EPSG","7030"]],                 AUTHORITY["EPSG","6326"]],             PRIMEM["Greenwich",0],             UNIT["degree",0.0174532925199433],             AUTHORITY["EPSG","4326"]],         PROJECTION["Mercator_1SP"],         PARAMETER["central_meridian",0],         PARAMETER["scale_factor",1],         PARAMETER["false_easting",0],         PARAMETER["false_northing",0],         UNIT["metre",1,             AUTHORITY["EPSG","9001"]],         EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs"],         AUTHORITY["EPSG","3857"]]     Origin = (-11688702.286414048000000,4866004.054033085700000)     Pixel Size = (0.057992644981843,-0.057992644981843)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=PIXEL     **Corner Coordinates:     Upper Left  (-11688702.286, 4866004.054) (105d 0' 5.04"W, 40d 0' 1.53"N)     Lower Left  (-11688702.286, 4865983.467) (105d 0' 5.04"W, 40d 0' 1.02"N)     Upper Right (-11688668.709, 4866004.054) (105d 0' 3.95"W, 40d 0' 1.53"N)     Lower Right (-11688668.709, 4865983.467) (105d 0' 3.95"W, 40d 0' 1.02"N)     Center      (-11688685.498, 4865993.760) (105d 0' 4.49"W, 40d 0' 1.28"N)**     Band 1 Block=579x3 Type=Byte, ColorInterp=Red       Mask Flags: PER_DATASET ALPHA      Band 2 Block=579x3 Type=Byte, ColorInterp=Green       Mask Flags: PER_DATASET ALPHA      Band 3 Block=579x3 Type=Byte, ColorInterp=Blue       Mask Flags: PER_DATASET ALPHA      Band 4 Block=579x3 Type=Byte, ColorInterp=Alpha      Also notice that the pixel size is suspect in the final image - both dx and dy are the same.


title: rgb2pct.py increases size of file 20 fold!

text: I wanted to reduce the number of bands of a tif file from 3 to 1... I used this command               rgb2pct.py -n 256 fnp_braunschweig.tif fnp_bs.tif      and it increased the size from 31mb to 778mb I actually want to decrease the size of the file. can anyone tell me why this is happening? cheers, Rob


title: ArcGIS Javascript API - IE9 Quirks Mode

text: I'm having difficulty understanding why ArcGIS Javascript API isn't rendering SVG when running in IE9 quirks mode. When I visit the ESRI sample site at http://help.arcgis.com/en/webapi/javascript/arcgis/demos/graphics/graphics_add.html and switch to IE 9 + Quirks mode I notice dojox/gfx/vml.xd.js loads but vector rendering does not work. Point, Multipoint, Extent all work fine but anything gfx related (polyline, freehand, line) does not work. Is it possible to use dojo/ArcGIS Javascript API in IE9 + Quirks mode? The reason I ask is because the pages I'm embedding the map in must run in this mode.


title: Splitting lines once at specified distance from each end?

text: This is related to an old problem of mine, regarding arrowheads on lines, which doesn't look very nice unless the line is very thin. It's even an ArcGIS idea with quite a lot of votes (although no apparent attention from ESRI). My current workaround consists of splitting the lines at a specified distance from the end, distance depending on the intended scale. Then change the symbology on the end part to contain arrows. Hopefully they will connect nicely to the rest of the line, which then can be a lot thicker. However, this time I need arrows on more lines than I care to edit manually, and I need arrows in both ends. So I figured out I need to split my lines at about 110 m from each end. How can I accomplish that? There are no vertices, points or anything at that location which I can use. All my lines are of different lengths. Every function I have tried (native ArcGIS functions, ET Geowizards) can only split the lines multiple times in equal segments, and only count from one side. I only need one segment at each end, and it has to be 110 m. I guess I could use the split function twice on each line, select the segments I need and be happy with that. It seems highly ineffecient though. How should I go about fixing this, or are there any other workarounds to get my arrows look nice? Python solutions could work as well. I have a Basic license, Python 2.7, Spatial analyst and 3D analyst to work with. I have a number of lines that looks like this: ![enter image description here](http://i.stack.imgur.com/aW6oc.png)


title: Which horizontal unit is used on the table of the Profile Tool?

text: Does anyone knows what horizontal unit is used on the table of the "Profile Tool". Vertical is in meter but I don't know about the horizontal axis. Knowing the length of my profile, I tried to divide that length by the number of extracted elevation but I noticed that the gap between each elevation data seems to change. At first I have one elevation data each 0.00453 for a hundreds of data. Then its one data per 0.0060 then 0.0105 etc etc. Would that mean that elevation data are not extracted at equally spaced points ? Thank you very much in advance ! Note: I am using as raster a .tif DEM and a line built in shapefile.


title: Projected coordinate system for South America

text: Does anybody have any suggestions for a projected coordinate system for the Northern part of South America? I have a data set that reaches from Venezuela to Southern Peru and I'd like to do some spatial autocorrelation on it. However, the current WGS84 is obviously a geographic coordinate system and not suitable for this type of analysis. I could just stick it in UTM 19 or 20 (north or south I guess given it spans the Equator) but I just wanted to know if anyone had any better suggestions.


title: How to save project with data source path as relative in QGIS?

text: Is their an option in QGIS to save the project with the data source being relative path. I Dont see the option here.QGIS Wroclaw1.7.3![enter image description here](http://i.stack.imgur.com/silhE.jpg)


title: Selection not displaying all results in QGIS

text: I have a point layer. Of which I have just selected 15 points from the attribute table manually and say show selected attributes. Of these 15 one point does not get highlighted the rest do fine. When I click on that feature with that identify option it does show me the attributes. What might the possible reasons ?


title: ArcGIS - Export raster keep color scheme

text: I have a *.flt raster loaded in ArcMap. I would like to export this raster as a georeferenced tiff file keeping the color scheme I have provided in the symbology. I right click in `TOC, Export Data` and export that raster data as a tiff file. The problem is the color scheme does not carry over. I have tried variations of "use rendered" and "force RGB" with no luck.


title: How do i create menu in base toolbar?

text: i'd like to create menus in a base toolbar like arcHdro tools or spatial analyst . how can i create menus with arcengine 9.3 in my toolbar ( the toolbar shows in arcmap )?


title: Is there any model builder in QGIS or other open source GIS?

text: Like the ModelBuilder of ArcGIS for Desktop to run simultaneous geoprocessing tools, is there any model builder available in QGIS or any other open source GIS software?


title: Backup several SDE database connections

text: This question stems from my previous question on how to list databases based on SDE. Now I want to copy them and have it written as a function so I can list out each database connection below it so it's easy to manage and update. Here's the code that lists the databases. I added the CopyFeature_Management, which may be different that what I want.               import arcpy     import os          arcpy.env.workspace = r"C:\Users\me\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Survey.sde"     def listFcsInGDB():         ''' set your arcpy.env.workspace to a gdb before calling '''         for fds in arcpy.ListDatasets('','feature') + ['']:             for fc in arcpy.ListFeatureClasses('','',fds):                 yield os.path.join(arcpy.env.workspace, fds, fc)          #Set workspace environment to geodatabase     #arcpy.env.workspace = r"C:\Users\me\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Survey.sde"          #Get list of feature classes in geodatabase     #FCs = arcpy.ListFeatureClasses()     FCs = listFcsInGDB()           #Loop through feature classes in list     for FC in FCs:               #Print the feature class name         print FC         print os.path.basename(FC)[0]         print os.path.splitext(FC)         output = "C:/Users/me/Desktop/"         arcpy.CopyFeatures_management(FC, output + os.path.basename(FC).split('.')[-1] + ".shp")              #os.path.basename(MXDPath).split('.')[0] + ".pdf"   .split('.')[0] + ".shp"      Any advice how to get the workspace and output as different variables listed below the function? I want them outputted to different folders based off of their database name. It works right now, but only for this one database connection. * * * **UPDATE 2** So now I have this working where it creates a GDB based on the SDE connection name. It will copy the contents and put them in the GDB. The only problem I get is it always leaves one FC out of the file geodatabase. if there are 19 features in survey it only copies 18 to the file GDB               import arcpy     import os          #Set workspace environment to the GDB locations.     arcpy.env.workspace = "C:/Users/me/AppData/Roaming/ESRI/Desktop10.1/ArcCatalog/"          #get every SDE connection we have. We can narrow it down by replacing the wildcard "*" with something else.     workspaces = arcpy.ListWorkspaces("*", "SDE")          #workaround to get feature classes within and outside feature datasets     def listFcsInGDB():         ''' set your arcpy.env.workspace to a gdb before calling '''         for fds in arcpy.ListDatasets('','feature') + ['']:             for fc in arcpy.ListFeatureClasses('','',fds):                 yield os.path.join(arcpy.env.workspace, fds, fc)               for workspace in workspaces:         arcpy.env.overwriteOutput = True         #list the current workspace we are looping through         arcpy.AddMessage(workspace)         #set the current environment to the SDE or GDB instance we want to loop through         arcpy.env.workspace = workspace         FCs = listFcsInGDB()          #Loop through feature classes in SDE or GDB         for GDB in FCs:             output = arcpy.CreateFileGDB_management("C:/Users/me/Desktop/", os.path.basename(GDB)[:14] + time.strftime('_%y_%m_%d_%H%M'))                  for FC in FCs:                 arcpy.env.overwriteOutput = True                 ##output = arcpy.CreateFileGDB_management("C:/Users/me/Desktop/", os.path.basename(FC)[:14])                 #output = "C:/Users/me/Desktop/" + os.path.basename(FC)[:14]+"/"                 #os.makedirs(output)                      arcpy.FeatureClassToGeodatabase_conversion(FC, output)                 print str(FC) + " Moved to " + str(output)      Any ideas why its moving them all but one from each? One of the GDBs has only one fc in it. It creates the GDB, but the database is actually empty. * * * **UPDATE 3**               import arcpy     import os     timestr = time.strftime("%Y_%m_%d_%H:%M")          #Set workspace environment to the GDB locations.     arcpy.env.workspace = "C:/Users/me/AppData/Roaming/ESRI/Desktop10.1/ArcCatalog/"          #get every SDE connection we have. We can narrow it down by replacing the wildcard "*" with something else.     workspaces = arcpy.ListWorkspaces("*", "SDE")          #workaround to get feature classes within and outside feature datasets     def listFcsInGDB():         ''' set your arcpy.env.workspace to a gdb before calling '''         for fds in arcpy.ListDatasets('','feature') + ['']:             for fc in arcpy.ListFeatureClasses('','',fds):                 yield os.path.join(arcpy.env.workspace, fds, fc)               for workspace in workspaces:         #arcpy.env.overwriteOutput = True         #list the current workspace we are looping through         arcpy.AddMessage(workspace)         #set the current environment to the SDE or GDB instance we want to loop through         arcpy.env.workspace = workspace         FCs = listFcsInGDB()          #Loop through feature classes in SDE or GDB              output = arcpy.CreateFileGDB_management("C:/Users/me/Desktop/", os.path.basename(workspace)[:14] + time.strftime('_%y_%m_%d_%H%M'))              for FC in FCs:             arcpy.env.overwriteOutput = True             ##output = arcpy.CreateFileGDB_management("C:/Users/me/Desktop/", os.path.basename(FC)[:14])             #output = "C:/Users/me/Desktop/" + os.path.basename(FC)[:14]+"/"             #os.makedirs(output)                  arcpy.FeatureClassToGeodatabase_conversion(FC, output)             print str(FC) + " Backed up to: " + str(output)


title: How can i use python QGis network analysis library to calculate short path based on speed and distance?

text: I want to calculate and show short path from starting point to an end point in a qgis python plugin. in my research i stopped by an article about QGis network analysis library; http://www.qgis.org/en/docs/pyqgis_developer_cookbook/network_analysis.html, and i found this: "To calculate edge properties programming pattern strategy is used. For now only QgsDistanceArcProperter strategy is available, that takes into account the length of the route. You can implement your own strategy that will use all necessary parameters. For example, RoadGraph plugin uses strategy that compute travel time using edge length and speed value from attributes." I'm wondering if i can define my own strategy based on (speed&distance) meaning time cost using QgsDistanceArcProperter. if somebody can explain to me if it's possible and how, i'll appreciate that.