title: config parser error raised by TileCache and Mapserver with OpenLayers

text: Just a forewarning, I will try to be as extensive with my question as possible. Due to a number of reasons why this might not be working, I have posted all sorts of portions of my debugging below. Web map located here Tile cache is being setup as a test run for the "hii_5" layer in OpenLayers. MapServer is the backend. I am trying to get TileCache to cache images of my MapServer WMS. However, no images are produced. The TileCache is "trying" to load the tiles (as you can see on the live page), however, no images are returned. When I follow the path of one of the image tiles, I see this error:               An error occurred: File contains parsing errors: C:\ms4w\Apache\cgi-bin\tilecache.cfg         [line 82]: '  type=WMS\n'         [line 83]: '  url=http://141.117.104.183/cgi-bin/mapserv.exe\n'         [line 84]: '  mapfile=C:/ms4w/Apache/htdocs/moot2_dev/hii_railroads.map\n'         [line 85]: '  extension=png\n'         [line 86]: '  layers=HII_railroads\n'       File "C:\ms4w\Apache\cgi-bin\TileCache\Service.py", line 85, in _load         config.read(files)       File "C:\Python27\lib\ConfigParser.py", line 297, in read         self._read(fp, filename)       File "C:\Python27\lib\ConfigParser.py", line 538, in _read         raise e        * I've run a GetCapabilities request and there are no warnings   * but GetMap request returns a blank page   * and the template return blank Here are the relevant mapfile, tilecache.cfg, and openlayers source: Not sure how much help I will get here, but the questions out there for anyone to take a stab at it! * * * **hii_railroads.map**               MAP                                                     #-> BEGIN MAP         OUTPUTFORMAT                                        #-> BEGIN OUTPUT FORMATTING DEFINITION             NAME "png"                                      # Name of format             DRIVER AGG/PNG                                  # Image driver             MIMETYPE "image/png"                            # Type of image             IMAGEMODE RGBA                                  # Image mode             EXTENSION "png"                                 # Image filetype             TRANSPARENT ON                                  # Transparency         END                                                 #-> END OUTPUT FORMATTING DEFINITION              WEB                                                 #-> BEGIN WEB DEFINITION                  METADATA                 "wms_title"           "HII_railroads"                 "wms_onlineresource"  "http://rinnerweb.arts.ryerson.ca/cgi-bin/mapserv.exe?map=/ms4w/Apache/htdocs/moot2_dev/hii_railroads.map&"                 "wms_srs"             "EPSG:3857"                 "ows_enable_request"  "*"             END                  TEMPLATE "template.html"             IMAGEPATH "/ms4w/tmp/ms_tmp/"             IMAGEURL "/ms_tmp/"         END                                                 #-> END WEB DEFINITION          LAYER                                                   #-> BEGIN LAYER DEFINITION: Railroads <-#              METADATA             "wms_title" "HII_railroads"         END              NAME    HII_railroads                               # Layer identifier         DATA    "data/wgs84/railroads_360.tif"              # Name of the data (in this case, string identifying a raster)         STATUS  ON           * * * **tilecache.cfg**               [HII_railroads]       type=WMS       url=http://141.117.104.183/cgi-bin/mapserv.exe       mapfile=C:/ms4w/Apache/htdocs/moot2_dev/hii_railroads.map       extension=png       layers=HII_railroads      * * * **index.html**               hii_5 = new OpenLayers.Layer.WMS("Railroads","../../cgi-bin/tilecache.cgi?",{         layers:"HII_railroads",         format:"image/png"},{         isBaseLayer:false,         transparent:true,         format:"image/png",         alpha:true,         displayInLayerSwitcher:false         });


title: How to handle "illegal data type (Float64)" error when trying to convert .tif to .bil?

text: Each time I try to convert (in QGIS 2.2) an elevation model from .tif to .bil format a get this error: ERROR 1: Attempt to create ESRI .hdr labelled dataset with an illegal data type (Float64). The same error occurs when I use the clipper with .bil as output file format. any hint or tipps what the error message means? thank you and best regards, Flo


title: How to make a popup containing links?

text: I have a list of links taken from a GeoJSON file. What I want to do, is program a popup with all these links inside, with OpenLayers. So, I wrote a function which returns a string containing HTML code for my popup, and that string is used by the "htmlContent" property of the Popup class constructor. I was asking myself if there was a better way to accomplish this task, such a callback in the popup, which calculates popup content, or templates ? Thanks


title: C# Library for converting Shapefile Projection

text: Is there any C# library which can take a shapefile as an input with projection information defined (e.g, EPSG 3011) and return a new shapefile with user specified projection (e.g, EPSG 4326) ?


title: how to load my local host(Mapserver) on my website

text: My localhost is http://localhost:8080/ and i am using mapserver 3.0.3 on this local host for making the application as Hello world, single map, more than one layer map, and so on. Now my folder is ms4w in C drive like:- c/ms4w and i want to upload this to my own server on internet, I thing it would be a simple application, and not running application as like my localhost, but in my server not have mapserver, it should be if yes how to run mapserver on my website and i want to run application as that I create in my local host. In my local host it's working fine, because its run on mapserver on my localhost. Hence everybody known this problem please please help me how can i solve the problem, its urgent. Thanks a lot in advanced. Ganesh


title: How to compile pgRouting 1.05 binaries for Windows?

text: Anybody with a compiled 1.05 binary for (Windows / PostgreSQL 8.4 / PostGIS 1.5.3) or a detailed description on how to compile?


title: How to render a table with mixed geometry types in QGIS?

text: I created a table with geometry type "GEOMETRY" AND inserted data of point, line and polygon type. My question is how to show this spatial data Qgis not supporting for geometry type "GEOMETRY".


title: How to set the font size of ITextElement if only know Height of the text?

text: When I read text elements from a .dgn file. I only can get the "height", "length" attribute of the text (it doesn't have a "size" attribute) (example: total height = 0.9; total length = 8.6). I need to add these text elements into a PageLayout of Arcgis. I used the **IFontDisp** to identify font attributes for text. It requires to identify the "size" attribute. How can I know that ?               ITextElement pTextElement = new TextElementClass();     IFontDisp pFont = new StdFontClass() as IFontDisp;     pColor = new RgbColorClass();     pColor.Red = 0;     pFont.Name = fontName;     pFont.Size =  **???**     pFont.Weight = weight;


title: How can i spatially add text to feature in argis 10?

text: I have alot of text that write with drawing tools and same feature . i want to find a method to add texts to features spatially like spatial join for feature to feature . is there any method to add these texts? ![enter image description here](http://i.stack.imgur.com/IU2FA.png)


title: Is there a way to draw a straight line on a map which also represents a straight line on the earth surface?

text: I'm trying to understand more on projections. ![enter image description here](http://i.stack.imgur.com/qkOgu.png) In this example, why must one calculate the curvature of the earth (great circle) if the map is already projected? Is there a way to draw a visibly straight line on that map which also represents a straight line on the surface of the earth? Why or Why Not?


title: OpenLayers.Control.DrawFeature Accuracy problem

text: When i draw a polygon using the OpenLayers.Control.DrawFeature then i export that polygon to KML using OpenLayers.Format.KML. Though, when I open it in google earth there is a big difference between the areas the polygon crosses even though the vertices of that polygon are exactly the same. **WHY?** How would i go about getting an accurate output to KML? ![enter image description here](http://i.stack.imgur.com/B1Fu7.png)


title: Large format Plotter Purchase

text: We are looking to purchase a second large format plotter for plotting land use planning maps and aerial imagery for field use. Specifically an HP DesignJet T790ps. Anyone have an opinion on this plotter? Also is there any issues that may have developed since your purchase. We currently use an HP DesignJet 3500CP for all the large plots from ArcGIS and some other graphics functions. We like the quality of this plotter but would like to have one that plots a little faster. Also what are your opinions on extended warranties? Are they really worth it?


title: ArcGIS javascript API - How To Use Tiled Map Services with Different Spatial References?

text: I have the following code fragment in an esri page to try and change basemaps, the following code works, **unless the next basemap layer uses a different spatialreference.**               if (baseMap)       _map.removeLayer(baseMap);     baseMap = new esri.layers.ArcGISTiledMapServiceLayer(mapSource);     _map.addLayer(baseMap);      It seems like my only option is to destroy the map instance and start from scratch every time the user desires to change the basemap. Am I going about this the wrong way... anyone have any insight? I have seen samples of the basemapgallerywidget, but haven't used it. So the basic question is: To change basemaps,not knowing the future basemap spatial reference, is my only option to destroy the map instance and re-create it? Thanks in advance!


title: How can I add users to the agsadmin and agsusers groups (ArcServer)?

text: I have a problem with adding users to the agsadmin and agsusers groups in ArcServer. How can I do the step? I've tried but I can't find it: http://help.arcgis.com/en/arcgisserver/10.0/help/arcgis_server_dotnet_help/index.html#/Configuring_a_one_machine_deployment/0093000000m6000000/


title: Cutting (Exporting) TIF images from WMS in ArcGIS 10 to use in AutoCad

text: The "ArcGIS for AutoCad" tool with WMS keeps giving FDO errors and is unreliable so we have procedure for putting cropped raster (TIF images) from ArcGIS into AutoCAD that has been used for years. We do not use ArcGIS's export world file function since ArcGIS doesn't allow for precise specific rectangle exports in Data View. We do all exporting from ArcGIS in Layout View. Basically will build a world file through several complex manual tasks in ArcGIS and AutoCAD. (i.e. shape ArcGIS's data frame to what is needed, write down coordinates from 2 corners of data frame in ArcGIS, draw a rectangle in AutoCAD based on these coordinates, insert the tif image, scale, etc.) Just wondering if anyone knows of a better way to accomplish this?


title: Is there a way to quickly reduce precision in a number of DOUBLE fields?

text: I was wondering if there is a way to reduce the precision for display of a DOUBLE field of a number of fields, in a shapefile. I have 25 fields and would like to do it as automatically as possible but I'm not certain where to start. I don't particularly want to create new fields and transfer over the values manually. I am displaying the text using dynamic text on data driven pages . I just want to display it with two rather than 5-6 decimal places. This is what the attribute table looks like, along with the relevant portion of layout: ![enter image description here](http://i.stack.imgur.com/gCjQ3.jpg) Here is the field properties box: ![enter image description here](http://i.stack.imgur.com/EhEH3.jpg)


title: Replicating ArcMap's Point Sketch Geometry Editing (Feature Templates) in ArcObjects

text: I'm looking for a way to place, select, move, delete features **inside a custom C# ArcObjects application** (using AxMapControl). I found the exact behavior in ArcMap, but I want to save the user some steps and avoid all of the dialogs, etc. It appears to be creating a point sketch geometry. In ArcMap, I can do the following steps:   1. Open a new map doc   2. Choose start editing from the Editor toolbar   3. Then I get a dialog asking me to choose the layer to edit   4. After selecting a layer, a panel appears that is titled "Create Features". Under it, there is another panel named "Construction Tools"   5. When I click on a feature in the upper panel, I get a "Point" option in the lower panel   6. Now when I move the mouse over the map, the cursor has the same icon as the feature I selected.    7. Clicking will place a feature, then I can drag or delete. So ... how can I replicate this behavior in a custom app using ArcObjects. I've already been able to replicate the behavior of other tools by doing something like this:               ICommandPoolEdit _commandPool = new CommandPoolClass();     ICommand solveRouteToolCmd = new ControlsNetworkAnalystSolveCommand();     _commandPool.AddCommand(solveRouteToolCmd, null);          ICommand command = commandPool.get_Command((int)solveRouteToolCmd );          //pass the MapControl to the oncreate method     command.OnCreate(_axMapControl.Object);     //run the tool     command.OnClick();      However, I can't seem to find a tool class object that corresponds to this behavior. In this example ControlsNetworkAnalystSolveCommand class is the hook to the same behavior that happens when you click the solve route button in ArcMap. I'm looking for that same hook for the scenario described above. Note that I found the tools by looking at the customize menu dialog. They are in the commands tab, in the editor category. The specific command is called Point and the description says, "Creates a point sketch geometry".


title: If clause within if clause VBA Arcgis10

text: I'm using ArcGIS10. I do not know what's wrong with this **VBA code** in the Field Calculator: *All the fields are String type               If [StartY]  <  [EndY]  Then       If [DOT_OSM] ="1" Then         Output = "F"       Else         Output = "T"       End if     Else if      If [DOT_OSM] ="1" Then         Output = "T"       Else         Output = "F"       End if     End If      ![Field calculator](http://i.stack.imgur.com/o3TsH.png) The Error from ArcGIS: "Syntax error" Thanks!


title: Can Mapnik 2.2 and Qgis 2.0.1 be installed together?

text: I'm using LM15 (== Ubuntu 13.04 Raring) with the Launchpad repos for Qgis Unstable 2.0.1 and Mapnik 2.2.0. Quantum wants libgdal1h; Mapnik wants libgdal1. As soon as I install one, the other gets uninstalled. Is it possible to install them both together? I'm willing to try compiling from scratch if that's the only way to get things working.


title: Big data rastering using ArcGIS Platform?

text: I have the daunting task doing raster analysis on a record set that runs into quite a few millions, and worse still is that the data is in **Hive**. So far I have managed to use record set of 2 million and create a raster but, in doing so, I consumed about 15Gb of my server's RAM and took around 5 minutes to complete the whole process. I am expecting the volume of data to increase a few fold. I have gone through the **big data analysis pages of ArcGIS** but that did not help me. Can somebody suggest a better option of doing big data analysis (perhaps by not fetching the data to the code, rather the reverse)?


title: How to call Geowebcache tile openlayers

text: Can anyone provide the solution, how to call GeoWebCache tiles in openlayers as i have 20GB of tiles created in gwc folder.


title: How can I import KML/KMZ files to Google Maps Engine Lite?

text: MapsEngine seems to import only CSV, XLS files for creating a map. How can I import my KML/KMZ files with Maps Engine Lite


title: FME Desktop 2013 Rastermosaiker: Why do I get strange edge patterns?

text: I have just mosaiked around 200 ecw tiles to a single ecw file. On inspecting the result, I see that I fme has produced a strange pattern around the edge. The Pixel values vary from R-G-B: 1 -1 -1 to 2-1-1 to 0-1-1 to 0-0-1and similar combinations very reaching over 2. ![Edge Pattern](http://i.stack.imgur.com/tkzfm.png) Can anyone offer an explaination and a solution so that this doesn´t happen? Thanks, Rob


title: Create subsetted multiband image from multiple files using gdal or python?

text: Given a time series of an area (single band), how you subset them and combine the subsets into a single multiband file (eg. vrt) using gdal or gdals' bindings for python?


title: What is the easiest way to let a user upload overlayers to openlayers (Geoserver/Mapserver)?

text: I need to let a user upload vector layers he created with QGIS or MapInfo professional to my website (using OpenLayers as client and either Geoserver or Mapserver as WMS). What will be the easiest way for a non programmer user to draw vector layer in program (one of the two I mentioned above) and be able to upload it to my server and add it to the map?


title: Openlayers KML Layer Styling

text: I am trying to change styling of Overlay KML layer in Openlayers, but I am not able to do so. here is my code:               <!DOCTYPE html>                    <title>OpenLayers KML as Base Map Example</title>


title: Load and display a large number of points using openlayers

text: I am using OpenLayers.Layer.Vector to load and display KML files. Some of the KML files have a large number of points in them (lets say around 20,000). I am surprised they load at all (in chrome), but it unfortunately does take too long. Would there be any way to speed it up (for example, parse the KML in parts).


title: Plot points using ArcGIS javascript api from co-ordinates

text: I would like to map co-ordinates that I got from the exif in photo's using my own command line tool that I made. As opposed to uploading your private photo's for everyone to see on flickr. How to I get the coordinates plotted on a map like this: http://help.arcgis.com/en/webapi/javascript/arcgis/jssamples/#sample/map_graphics Using a file generated from the stdout of my command line tool, here is an example: http://dpaste.com/880696/plain/ I can't find any info on how esri.request works in the Javascript API for arcGIS. I also tried looking at this loading text sample in arcgis javascript api's sample section, the example doesn't seem to work. Basically all I'm asking is what is the simplest method to plot co-ordinates on the map as they are doing in the example with flickr, but from a text file. Perhaps I should first turn them into the format that is being used in the example with flickr? I tried using this format: on georss dot org but using sed to replace the values with my own values but it didn't work.


title: How to set up a maximal zoom for polygons in Leaflet?

text: How to set up a maximal zoom for some layer, for example state polygons? I need that my US polygons which represent states have maxZoom 6, and when zoomed higher than 6, that US counties show up instead of states. Is that possible in Leaflet? All states are defined like Leaflet polygons. Thanks in advance.


title: Google Map layer loads tiles even if this layer is disabled in OpenLayers

text: Sorry for cross posting I've asked this question in OL mailing list, but not get response yet. It is very important for me to resolve this issue. I've made a little example that reflect my problem. It's not clear for me why GM layer loads tiles even if this layer is disabled (it can be observed against firebug in Firefox) and how I can to prevent this?


title: Read specific column and row value from a text file using Python

text: I have a text file as following:               1.25        177.             1.50        220.             2.00        272.             5.00        427.         10.00       556.             25.00       763.             50.00       916.              and i want to read second last value from second column. I can move onto line of interest but somehow can't figure out a way to read required value. Any help would be appreciative.               a = open("test.txt","r")     b = a.readlines()     a.close()          count = 0     for line in b:         count += 1         if count > 6:             print line,


title: How to retrieve the mapPoint behind a PictureMarkerSymbol?

text: I am using ArcGis Javascript API 3.2. After instantiating a Point, PictureMarkerSymbol, attribite and an InfoTemplate, I use this code to plot a point on the map.               var graphic= new esri.Graphic(candidatePoint,candidateSymbol,attr,infoTemplate);     map.graphics.add(graphic);      I want to retrieve the point behind the graphic when it is being clicked. How will I do it? Thanks, Prasanna


title: Problem with server jsp geoserver

text: I want run jsp code in geoserver 2.2.3! I put page jsp in data_dir/www but it not run, html page run fine in data_dir/www! Please help me! Thanks!!


title: How to convert a finished project in QGIS to PDF

text: I have finished my project and need to convert it into a PDF file. How can I do it?


title: ArcGIS for INSPIRE, build Protected sites Inspire layer

text: Can anyone add some infomation about building protected sites inspire layer? I cannot understand the whole specification schema and I didn't know how to use it in a real database.


title: GeoServer renders shape incorrectly/distorted

text: I am having trouble displaying country borders properly in GeoServer, as you can see from the link below (map of Italy). http://trideg.server.de/tridec/wms?service=WMS&version=1.1.0&request=GetMap&layers=tridec:static_hasc_regions&styles=&bbox=489196.9,4109254.4,2426416.9,5987770.7&width=775&height=330&srs=EPSG:900913&format=application/openlayers&viewparams=codes%3A%27IT%27 The problem occurs with country boundaries which have very sharp edges (I assume), e.g. the shape of Albania shows correctly: http://trideg.server.de/tridec/wms?service=WMS&version=1.1.0&request=GetMap&layers=tridec:static_hasc_regions&styles=&bbox=489196.9,4109254.4,2426416.9,5987770.7&width=775&height=330&srs=EPSG:900913&format=application/openlayers&viewparams=codes%3A%27AL%27 I obtained the ESRI shapefile from naturalearthdata.com and used ogr2ogr to transfer them in a MySQL database, which GeoServer uses as a datasource and created a layer in it. Anyone has a lead? Either I missed something in the layer configuration (the projections are alright though, as the shapes are positioned correctly in the map), or ogr2ogr messed up the data during the conversion.


title: QGIS saving analysis output as Postgis table

text: i am using **QGIS and Postgis** together. I can load data in Postgis and than i import it to QGIS. But i am trying to **save an analysis output (example: buffer analysis, voronoi polygon or SEXTANTE analysis tools) directly in Postgis**. Is there anyway to make this? For instance in ArcGIS you can save an analysis output in a **.GDB .MDB geodatabase** but in QGIS firstly i save an output as .SHP than i load it to Postgis or in OpenJump after analysis you must save layer in Postgis table. So there is one more step !


title: Dubai - Sector/Community GIS SHP Data?

text: Where can I download free or open access GIS/SHP file datasets for Dubai, including their respective spatial boundaries? If these might not be available, where can I find a map outlining the Sectors or Communities in Dubai? I double checked the sources I found from Wikipedia ( http://en.wikipedia.org/wiki/List_of_communities_in_Dubai ), but wasn't able to find anything as none of the links worked. I also tried OpenStreetMap, but only was able to get the road network. This doesn't show/depict community/sector outlines. I want to match up the population statistics from this document (listed by sector level and community level) with a GIS file that shows each of these respective sectors/communities - http://dsc.gov.ae/Publication/Population%20Bulletin%20Emirate%20of%20Dubai%202012.pdf


title: animation contour with Time Manager in QGIS?

text: I have a temperature log for a number of town for a day: Extract from log: Logitude Latitude DateTime temperature 26.418131 -0.6014717 2014-02-03 00:00:04 14.3 26.418131 -0.6014717 2014-02-03 00:00:09 14.6 26.418131 -0.6014717 2014-02-03 00:00:14 14.6 26.418131 -0.6014717 2014-02-03 00:00:19 14.5 26.418131 -0.6014717 2014-02-03 00:00:24 15.0 26.418131 -0.6014717 2014-02-03 00:00:29 14.3 the time between two measurement is 5 seconds. For every point wants to build contour intervals (isopleth map) :10, 15, 20 degrees my goal to show how it is changing with time is to say make a movie, I've read that I should use plugin Time Manger ,have nickel install it works but i can't manage to couple them with ,if anyone can help me or clear me the way i am very grateful![enter image description here](http://i.stack.imgur.com/6Tqcd.png) here is my initial layer (initial data longitude-latitude)


title: What in-depth spatial database systems tutorials exist?

text: Is there a good tutorial that explains in depth the internals of a GIS and a spatial database system such as PostGIS with examples and without any background assumption? I'm specifically looking for answers to questions such as:   1. What is a geometry in GIS?    2. How is it represented?    3. Given a lat/long coordinate, what are the operations I need to perform on it to bring it to a state where I can call something like ST_contain PostGIS to check if the polygon contains the lat/long coordinate? I have some experience with general relational database systems but no background with GIS or the spatial database paradigm.


title: ArcGIS 10 export metadata to XML programmatically

text: I have some code that exports metadata from an iMetadata instance to a standalone XML file. This code uses the ExportXML type and the iMetadataExport interface (see below).               Private Sub ExportMetadatatoXML(ByVal Metadata As IMetadata, ByVal strXMLdocName As String, ByVal strOutputXmlFolder As String)         Try                  Dim ExportXML As ExportXML             Dim xml_Export As IMetadataExport             ExportXML = New ExportXML             xml_Export = ExportXML             If System.IO.Directory.Exists(strOutputXmlFolder) Then                 xml_Export.Export(Metadata, strOutputXmlFolder & "\" & strXMLdocName & ".xml")                 Logging.WriteToLog("Metadata Exported to " & strOutputXmlFolder & "\" & strXMLdocName & ".xml", Serverity.SeverityLevel.Information)             Else                 Logging.WriteToLog("Output XML Directory " & strOutputXmlFolder & " cannot be found", Serverity.SeverityLevel.Fatal)                 Result = Result - 1             End If                  Logging.WriteToLog("Metadata was not exported to " & strOutputXmlFolder & "\" & strXMLdocName & ".xml", Serverity.SeverityLevel.Information)              Catch ex As Exception             Dim strRoutine As String = "Error thrown in: " & System.Reflection.MethodBase.GetCurrentMethod.Name             Logging.WriteToLog(strRoutine & " - " & ex.Message, Serverity.SeverityLevel.Fatal)             Result = Result - 1         End Try          End Sub      This code worked a treat for 9.x versions of ArcGIS, but with the new metadata changes for 10+, it no longer works. The ExportXML type seems to no longer exist and I can't seem to figure out a replacement. Does anyone know how I can achieve this with 10+? I've noticed there's a metadata toolbox but would rather avoid using geoprocessing. Thanks for your time. Edit: I have since tried using geoprocessor to export to XML, but i'm getting a hresult error on gp.execute. Can anyone suggest what the problem might be?                       Dim GP As ESRI.ArcGIS.Geoprocessor.Geoprocessor = New ESRI.ArcGIS.Geoprocessor.Geoprocessor             Dim parameters As IVariantArray = New VarArray             parameters.Add("C:\Temp\Temp.gdb\BLOCKS")             parameters.Add("C:\Program Files\ArcGIS\Desktop10.1\Metadata\Translator\ARCGIS2FGDC.xml")             parameters.Add("C:\Temp\Metadata\blocks.xml")             GP.Execute("ExportMetadata_conversion", parameters, Nothing)


title: Map automation, create and populate inset maps for densely-labelled areas

text: I would like to know if it is possible to automate the example in the picture below, unlike the picture below i'm using vector data. I would like for a fixed scale generate map, but when labels are too dense, we could have some little maps (inset maps) on the side of the Main map. i have found some tips using arcpy. but the **problem is how to detect densely labelled areas** ![enter image description here](http://i.stack.imgur.com/Q0Hti.png)


title: Reprojecting from MGRm to WGS84 using open-source tools

text: I would like to work with some of the Michigan state data (for example, this shapefile of cities ), but it is not in a common projection. However, when I try to reproject the file using ogr2ogr and any one of the possible EPSG codes, I get a map that is anywhere from rotated & skewed to flat:               ogr2ogr -s_srs EPSG:102121 -t_srs EPSG:4326 city_mi_84.shp city_mi.shp       Any suggestions on how I can make this work? I don't have access to Arc at the moment.


title: Algorithm to find polygons enclosing points

text: I'm trying to find an algorithm that can determine the smallest possible polygons to cover a number of points. I know how to get the convex hull around all the points, but say that the points are located on different islands, is it possible to determine that there is a gap between different groups and get separate polygons for each group?


title: gdalwarp results all black

text: I'm trying to convert a bmp image in BNG to WGS84 using gdalwarp (via QGIS). It seems to be running without errors but the output is all black. What am I doing wrong?


title: Why is my controlling shell script not working?

text: **I updated my post after I resolved the previous errors, but additional errors were produced.** I want to create a controlling shell script that contains water.mapcalc within it. The problem is, I do not know how to save a code in a file called water.mapcalc. I'm trying out the hydrologic simulation. It is stated in page 12: "For those wishing to repeat the experiment with GRASS enter the model code in a file called water.mapcalc then create a controlling shell script." I understood the instruction as I should save the following in a raster map called water.mapcalc:               r.mapcalc water.mapcalc = "water = water + eval(x=elev + water,         if(x>(y=elev[-1,0]+water[-1,0]),             -.15*if(elev>y, water, x-y),             .15*if(elev[-1,0]>x, water[-1,0], y-x))+              if(x>(y=elev[1,0]+water[1,0]),             -.15*if(elev>y, water, x-y),             .15*if(elev[1,0]>x, water[1,0], y-x))+              if(x>(y=elev[0,-1]+water[0,-1]),             -.15*if(elev>y, water, x-y),             .15*if(elev[0,-1]>x, water[0,-1], y-x))+              if(x>(y=elev[0,1]+water[0,1]),             -.15*if(elev>y, water, x-y),             .15*if(elev[0,1]>x, water[0,1], y-x))+              if(x>(y=elev[-1,1]+water[-1,1]),             -.10*if(elev>y, water, x-y),             .10*if(elev[-1,1]>x, water[-1,1], y-x))+              if(x>(y=elev[1,1]+water[1,1]),             -.10*if(elev>y, water, x-y),             .10*if(elev[1,1]>x, water[1,1], y-x))+              if(x>(y=elev[1,-1]+water[1,-1]),             -.10*if(elev>y, water, x-y),             .10*if(elev[1,-1]>x, water[1,-1], y-x))+              if(x>(y=elev[-1,-1]+water[-1,-1]),             -.10*if(elev>y, water, x-y),             .10*if(elev[-1,-1]>x, water[-1,-1], y-x)))"      **Should I save the model code in a water.mapcalc raster map as what I did above? or should I save a water.mapcalc file in the /bin folder?** I tried both but it still produced errors. I already did the said map (water.mapcalc). and I created the following shell controlling script that needs the file (water.mapcalc). This is the shell script (myscript5.sh) and I saved it in the \bin folder of GRASS:               #!/bin/sh     export GRASS_PNGFILE=/tmp/grass6output.png          r.mapcalc water = 120     d.rast water          r.mapcalc < water.mapcalc     d.rast water          r.mapcalc < water.mapcalc     d.rast water          g.copy rast=water,water.1          r.mapcalc < water.mapcalc     d.rast water          r.mapcalc < water.mapcalc     d.rast water          g.copy rast=water,water.2      Then, I run the `sh myscript5.sh` in the CMD of GRASS and it produced this message:               PNG: GRASS_TRUECOLOR status: TRUE     PNG: collecting to file: C:/Users/Rhen/AppData/Local/Temp/grass6output.png,     GRASS_WIDTH=768, GRASS_HEIGHT=508     /c/Program Files/GRASS-64/bin/myscript5.sh: line 7:     water.mapcalc: No such file or directory     PNG: GRASS_TRUECOLOR status: TRUE     PNG: collecting to file: C:/Users/Rhen/AppData/Local/Temp/grass6output.png,     GRASS_WIDTH=768, GRASS_HEIGHT=508     /c/Program Files/GRASS-64/bin/myscript5.sh: line 10:     water.mapcalc: No such file or directory     PNG: GRASS_TRUECOLOR status: TRUE     PNG: collecting to file: C:/Users/Rhen/AppData/Local/Temp/grass6output.png,     GRASS_WIDTH=768, GRASS_HEIGHT=508     Copy raster <water@PERMANENT> to current mapset as <water.1>     /c/Program Files/GRASS-64/bin/myscript5.sh: line 15:     water.mapcalc: No such file or directory     PNG: GRASS_TRUECOLOR status: TRUE     PNG: collecting to file: C:/Users/Rhen/AppData/Local/Temp/grass6output.png,     GRASS_WIDTH=768, GRASS_HEIGHT=508     /c/Program Files/GRASS-64/bin/myscript5.sh: line 18:     water.mapcalc: No such file or directory     PNG: GRASS_TRUECOLOR status: TRUE     PNG: collecting to file: C:/Users/Rhen/AppData/Local/Temp/grass6output.png,     GRASS_WIDTH=768, GRASS_HEIGHT=508     Copy raster <water@PERMANENT> to current mapset as <water.2>      the error is: **water.mapcalc: No such file or directory**. I think there is a problem on how I saved the water.mapcalc file. I do not know what format (raster/.txt/.sh file) and where will I save it (should I also save it in the /bin folder together with my script). I tried all of those and it still did not work. I cannot figure out where the problem is. Thanks.


title: How to reduce the resolution of satellite imagery data for use in MATLAB

text: When I am using the 4-band multispectral data for NDVI calculation, after successful completion of runtime, it shows some warning message like "Loading 17%..." How should I rectify it? Should I do any filtering for reducing the resolution effect before processing? Thanks in advance.


title: Dissolving polygons by attribute and in a certain distance (ArcGIS)

text: in a FeatureClass I have a lot of polygons which contain attributes of different landuse-types. Now I am trying to dissolve those polygons of the same landuse-type, but only if the distance between them is less than 100m. I'm trying this with ArcGis/ArcInfo/Version 10. Thanks a lot for suggestions!


title: How to add offset to layer in OpenLayers?

text: As the title suggests,I have GeoJSON vectors on OpenLayers map where user should add X/Y offset (i.e 5 meters north,5 meters east).I plan to add some control for it. Since vectors are drawn map with **OpenLayers.Rule** object,as I searched for some offset option,then all I found was graphicXOffset & graphicYOffset which eludes normal vectors but only used for external graphics. Is there any easy way to accomplish this feature on client ?


I'm having issues getting my `Esri.FileGDBAPI.DLL` to work. I have the other two files, `FileGDBAPI.dll, FileGDBAPID.dll` in my `bin` folder. I can only add the Esri dll to my references, the other two I get an error saying that they can't be added because they either aren't accessible, not COM components or valid assemblies. This is the error I get when running it:               A procedure imported by 'Esri.FileGDBAPI.DLL' could not be loaded.      By the way, I also have the `FileGDB_DotNet.dll` file in my `Windows/system` directory and I'm using .Net4.0. Thanks!


title: Create raster distance around a single point feature class

text: Okay, so I'm relatively new to ArcMap and I'm trying to use the Fuzzy membership for suitability test. I have a single point feature class that I want to use for this analysis. I know that in order to be able to apply the Fuzzy membership you need to have a raster input. My question is how do I create a raster distance around my point on order to apply the fuzzy membership after? I tried all the Distance (Spatial Analyst) tools, but none seems to work for this situation. Help is very much appreciated :)


title: Writing length of line along the line

text: I sometimes have to write out the length along a line for every couple of hundred meters or kilometers. This applies to pipes, when the construction team wants to know where along the pipe they are or talk about.   So far I have done this in a very unefficient way:   * dissolving line segments (if there are multiple along the way, and there usually are)   * split them one piece at the time at specific distances using Editor --> Split   * Convert endpoints to point features   * Generate list of numbers in Excel to paste in a field in the point FC   * label the points without symbol (or with a cross symbol as below) This works, but there has to be more suitable ways of doing this. I'm not specifically looking for a programmatic/Python solution (even though it is welcome), anything that is a tad less complicated that the outlined workflow will do.. :) ![enter image description here](http://i.stack.imgur.com/zGFSX.png)


title: Will there be an iOS version of QGIS?

text: I saw that QGIS now works on Android, any chance there will be an iOS version?


title: How to Zoom and Pan using scroll wheel in ArcObjects?

text: is there any way to zoom using scroll wheel ? i have tried `AutoMouseWheel` of map control, it works well for the first used custom tool, but once i change the tool, it stops working , i would like to know the reason of this weird behavior. i don't know if it is possible to catch mouse wheel events in custom tool ? that's seems to be an alternative solution. Any suggestion or help is welcome :)


title: Postgis ST_AsGeoJSON() to openlayers security issues?

text: Im about to implement a python script serving my openlayers code with GeoJson data from PostGIS. However, here How to connect OpenLayers to PostGIS data? I see " Unfortunately you cannot connect a web page directly to a database because of security concerns " Could you please tell me about those concerns? Isnt using a proxy /server script good enough? PS. I would ideally like to post in the above tread but lack the 50 reputation neccesary.. sorry about that.


title: Error while installing PostGIS on ubuntu 12.04

text: This is my first practice on Ubuntu machine. Here i am trying to install PostGIS 2.1 on postgresql server 9.3.2. When i am following the procedure mentioned here Help I am getting the following error messages on the terminal window.               You might want to run 'apt-get -f install' to correct these:          The following packages have unmet dependencies:      pgagent : Depends: libkrb53 (>= 1.4.2) but it is not going to be installed                Depends: libssl0.9.8 (>= 0.9.8a-1) but it is not going to be installed      postgis : Depends: libgdal1 (>= 1.9.0) but it is not going to be installed                Depends: libgeos-c1 (>= 3.3.3) but it is not going to be installed                Depends: liblwgeom-2.1.1 (>= 2.0.0) but it is not going to be installed     E: Unmet dependencies. Try 'apt-get -f install' with no packages (or specify a solution).      Expecting some hints from anybody to resolve the issue.


title: How to clip a raster with vector boundaries?

text: I am using QGIS. I would like to clip a raster precipitation layer using an admin boundary layer that is vector data. However the geoprocessing tools seem to be usable only for vector data. How can I clip this precipitation layer? Thanks for any help!


title: I would like to convert trimble ssf format to shp format

text: I would like to convert a small amount of data from trimble ssf format to shp format for inclusion in an ArcGis database. I have been trying to do this for a few months and not success. I have a small amount of data so dont want to purchase a license for pathfinder. Any help or if anyone can do this appreciated. Thanks


title: How to show labels in arcglobe for extruded features

text: I am trying to visualize lithology data in ArcGlobe, i.e, the representation of different types of soil data under the earth. This I have done properly by creating Z aware featureclass and extrusion properties. I have attached a picture, a representation in Globe. In this diagram each color in the 3D bars is one feature representing the type of soil. Here I want to show the labels from the attribute table having the depth information of each feature. But when I am applying extrusion automatically, the Show Labels menu gets disabled. Is there any way to add labels programmatically (vertically next to each feature) like 3D text elements etc? Give me some ideas please. Thank you.... ![screen shot of 3D representation](http://i.stack.imgur.com/sYdba.jpg)


title: IEditSketch and IEditSketch2

text: What is the difference between IEditSketch and IEditSketch2? Also why there are same interfaces with the same name only with a difference "2" at the end of the interface's name ?


title: Defining a scale number from the Helmert transformation

text: I know this is more bit of a cartographic issue, but maybe somebody would be able to help. I have an old map georeferenced to the openstreetmap with 50 identical points. I want to get the scale of the old map. I am using Helmert 4-parameter tranformation (similarity) with the equations:   _X=ax-by+X0_   _Y=ay+bx+Y0_   where _x_ and _y_ are the old map coordinates in meters (same as coordinates used from openstreetmap), _X0_ and _Y0_ are translations and then there is one scaling and rotation. Scaling parameter _m_ could be obtained as:   _m=sqrt(a^2+b^2)_   Rotation is:   _alpha=arctan(b/a)_ The resolution of the old map is 400 dpi (0.0635 mm per pixel). I have all the parameters counted now with the least squares method.   Does anybody know how to get the scale of the old map from these parameters? I want the scale in a common form 1 : M (1 : 600 000).   _EDIT_ : Further on I will use iteratively reweighted least squares and I will be omitting those identical points with the biggest spatial error (I am working with very old maps which have mistakes). I will do this to make the definition of scale more accurate. That's why I need to use the transformation parameters for the map scale definition and not only to compare the distances.   Thanks a lot.


title: pivot table output

text: In an attempt to summarise landclasses, the columns of my pivot table output are a bit mixed up with the field info regarding the three input fields that I've selected. As shown in snapshot 3 the description field have the info of the grid_slo field, the grid_slo field have the info of the grid_alt field, and grid-alt showing zero. The original descriptions that should show under the descriptio field has been lost. Numerous attempts to fix this problem has been unsuccessful. Can anyone help please? ![attribute table](http://i.stack.imgur.com/3ta5U.png) ![pivot table input](http://i.stack.imgur.com/gueFu.png) ![snapshot showing the output, with the description field having the info of the grid_slo field, the grid_slo field having the info of the grid_alt field, and grid-alt showing zero. The original descriptions that should show under the descriptio field has been lost](http://i.stack.imgur.com/3nwD5.png)


title: Missing date column from interactivity callback

text: I can't seem to get "date" columns to return for popups/infoWindows. I can retrieve other data types without an issue. I have the interactivity set correctly from what I can tell, but the featureClick callback only returns one column. My code relies on v2 of the CartoDB API. Here's my call to cartoDB.createLayer: cartodb.createLayer( map, { type : "cartodb", options : layerData } ) where layerData is               interactivity: "name,admin_date,size"     query: "SELECT size, name, admin_date, the_geom_webmercator FROM book1"     table_name: "book1"     tile_style: "#book1{marker-width:5;marker-opacity: 0.9;marker-line-color: white;marker-fill:#e6e6e6;marker-line-width: 1;marker-line-opacity: 0.9;marker-placement: point;marker-type: ellipse;marker-allow-overlap: true;} #book1{[size>=0]{marker-fill:#b8cde9;}[size>=18229430]{marker-fill:#6693cc;}[size>=22061458]{marker-fill:#2f5f9d;}}"     user_name: "c2es"      The featureClick callback only return the first and last columns. Any ideas what's wrong?


title: How to close spaces within a polygon in QGIS 2.2?

text: I have a shapefile displaying the borders of several city units. One of its polygons, however, had free space consisting of four vertices in it. I created another polygon inside the free space and then used `Dissolve` to close the free space. The result was the following layer which still contains some errors (see the close-ups below). What is the best way to resolve this in QGIS 2.2? Using some kind of `v.clean`? Thank you for your help. ![enter image description here](http://i.imgur.com/lqqUSoe.jpg) ![enter image description here](http://i.imgur.com/GtecT6e.jpg) ![enter image description here](http://i.imgur.com/tKzyup1.jpg) ![enter image description here](https://i.imgur.com/dkLf9h0.png) ![enter image description here](https://i.imgur.com/FxdGrAO.png) ![enter image description here](https://i.imgur.com/tT12UPE.png) Edit: The entire data seems to be corrupted. Using `Dissolve` on the entire shape results in what the fourth screenshot shows. Edit 2: Added fifth screenshot. Edit 3: Added sixth screenshot.


title: QGIS export ESRI Shapefile to SpatiaLite format issue

text: I am trying to export a shapefile of 3221 mixed POLYGON and MULTIPOLYGON geometry elements and every time, I get this same error: ![enter image description here](http://i.stack.imgur.com/Bv6K2.png) Does anyone have any idea what's going on? And no, this is not an exact duplicate of http://stackoverflow.com/questions/15203403 as I am working solely within QGIS. Any help is greatly appreciated. Thank you! - Brandon


title: Is there a way to get centroid of geographic object in MSSQL?

text: I have a table of OSM geographic roads in MSSQL 2012. I want to create a new column with the centroid of each road (So I can set center on OpenLayers map for each road). How can I do it ?


title: How to get the best accuracy of coordinates for web GIS- application?

text: I would like to get advice from you about one thing. I would like to develop my own web GIS application. It would be about some sport places in my hometown; it would be good practice for my education - I'm still a student. First, I need to get coordinates of the features, then I think to put it in AutoCad, then export to .shp file, then work with it in QGIS and use Postgres and PostGis for databases and than use Mapserver and so on. But I need help with the first step, getting the coordinates. I don't have any high accurate instrument like some good GPS, so is there any chance to get coordinates of the places by Google Earth? How accurate is that and is it even possible to do it like that? I'm looking forward to someone's answer and please don't get mad if I said something stupid because as I said, I'm in the beginning of getting to know all this GIS web development and things like that. Thanks! * * * Thank you all for your answers, it all helped me in some way. The only thing that I'm not sure is accuracy of osm and bing maps. I was editing some stuff in JOSM recently and I've realized that bing maps don't really fits in the right place if you compare it to gpx data. (there is at least 2m difference) I'm not still quite sure is it accurate enough for me. I would like to achieve 1m maximum level of accuracy.. :/


title: v.dissolve to multipolygon

text: Using v.dissolve in GRASS, I'd like to get a multipolyon as a result. However, I end up with distinct polygons with the same 'dissolve value.' For example; lets assume I have two polygon features which are not touching each other spatially. Both have a 'DISS" value of 1. If I run v.dissolve on this DISS value, I get the same thing in the output as the input (i.e. two seperate polygons with a DISS value of 1). What I would like to be able to get is an output that contains one single multipolygon. How do I go about this?


title: Will QGIS fulfill my requirements? moveable, rotatable, draggable, complex, interactive icons

text: I'm trying to compile a list of mapping software that will fit my requirements for a desktop application I want to build. I'm trying to build something similar to MissionPlanner: https://code.google.com/p/ardupilot- mega/wiki/Mission Here is a quick moc-up I did for my application: http://i.imgur.com/52DaJtQ.png Requirements:   * I want to be flexible with what I can draw on the map: lines, rectangles, polygons    * I want to be flexible with icons: they should be able to move, rotate   * I want flexibility with icon interaction: for example, being able to easily drag a waypoint would be great   * I want to be able to use my own data for maps: for example an aeronautical chart made from a folder of .tiff images   * I want to sell my application commercially   * I don't have a preference on programming language yet Is QGIS right for me?


title: How to create a catchment area using LSOA and point file data

text: I have to create a catchment area (based on LSOA) for the 19 underground stations within Camden borough and it has to be in a polygon vector file. I know I have to somehow group the SOAs but don’t know which geoprocessing tool to use.


title: Arcpy CalculateField using formatting syntax

text: I have 2 feature classes (FC) that I am merging together and then performing some field calculations on. The 1st FC has 100k+ records in it, and the second FC is a blank FC that I'm using as a template. I need to merge the two FC's together and then field calc the values from FC1 into FC2.   The problem with FC2 is that the fields are required to often be of a different type. For example, in FC1 field "PID" is a double and in FC2 "PIN" is a string. I need to calculate the values of "PID" into "PIN" with a prefix of "139-". This pattern is repeated several times with different fields. My simple code example:               data = r"C:\Temp\Test.gdb\PAR_PROPINFO_SCO_VW2"         template = r"C:\Temp\Test.gdb\Template"         out = r"C:\Temp\Test.gdb\Scott_CO"         arcpy.Merge_management([data,template],out)         qry= '"139-"' + "+'PID'"         arcpy.CalculateField_management(out,"PIN",qry)      In this particular instance, qry threw a syntax error. If I just use:               arcpy.CalculateField_management(out,"PIN","PID")      I get a type mismatch error. Since I'm on Arc 10.0, I was trying to avoid a cursor with the amount of records involved. This process will run daily and the # of records will only increase. I was thinking I could just get away with a brute force merge and calc method. However, I run into type mismatches. I thought field mappings, but that is beyond my current level.


title: Unusual results for geoprocessing speed tests

text: I have been observing unusual performance with a Python geoprocessing script. The (attached) script performs the following actions:   1. Use a search cursor to look up UTM zone corresponding to polygon features   2. Create spatial reference object based on search cursor results   3. Convert .csv to feature layer and then to a point feature class I've noticed markedly different processing times based on how the script is run:   * **32-bit** processing using **IDLE** = 203 seconds   * **32-bit** processing foreground **script tool** = 91 seconds   * **64-bit** processing background **script tool** = 206 seconds _Why would this script perform so differently given the above conditions?_ I certainly would not expect the 32-bit script tool run in the foreground to be 2X as fast as the other methods. * * *               import arcpy, os, time          ###IDLE Parameters     ##fc = r'C:\path\to\polygon\fc\with\utm\zones\and\features'     ##outws = r'C:\out\location'     ##arcpy.env.workspace = r'C:\workspace'          ####################     ## Script tool parameters     fc = arcpy.GetParameterAsText(0)    # Feature class     outws = arcpy.GetParameterAsText(1) # Folder     arcpy.env.workspace = arcpy.GetParameterAsText(2)   # Workspace     ####################          # Tables are .csv     tables = arcpy.ListTables()          start = time.clock()          # Look up which UTM zone .csv features are in     for t in tables:         quad = t[7:17]         print quad         whereClause = """ "QUADID" LIKE '%s' """ % quad         with arcpy.da.SearchCursor(fc, ("QUADID","ZONE"), whereClause) as cursor:             for row in cursor:                 if row[0] == quad:                     utmZone = row[1]                     if utmZone == 10:                         sr = arcpy.SpatialReference(26910)  # NAD_1983_UTM_Zone_10N                     elif utmZone == 11:                         sr = arcpy.SpatialReference(26911)  # NAD_1983_UTM_Zone_11N                     elif utmZone == 12:                         sr = arcpy.SpatialReference(26912)  # NAD_1983_UTM_Zone_12N                     elif utmZone == 13:                         sr = arcpy.SpatialReference(26913)   # NAD_1983_UTM_Zone_13N                     else:                         print "The UTM Zone is outside 10-13"                 else:                     pass              # Convert .csv to feature class         try:             outLayer = "in_memory"             # Now with the sr defined, create the XY Event Layer             arcpy.MakeXYEventLayer_management(t, "x", "y", outLayer, sr, "z")             arcpy.FeatureClassToFeatureClass_conversion(outLayer, outws, t[7:17])             arcpy.Delete_management("in_memory")             end = time.clock()             print "In_memory method finished in %s seconds" % (end - start)              except:             # Print any error messages             print arcpy.GetMessages(2)          print "Processing complete"


title: how to open a layer's table in a map document

text: Trying to create a map document programmatically. Added a featureclass as a layer. When I open the map document in ArcMap, I want the table for that layer to be open. Found `TableView` (doc) and `AddTableView` (doc), but these seem to be only for standalone tables. I also looked at the properties for Layer (doc) and for DataFrame (doc) and didn't see methods that looked right. What I'm shooting for is to double-click on the mxd and have something like image below appear. ![enter image description here](http://i.stack.imgur.com/jaEzR.png) Any ideas? Thx


title: How to create a multi-colored raster image based on long,latt,color information?

text: I want to create a multi-colored raster layer to overlay on my map in OpenLayers I have the following information in my DataBase:               D_Cell_ID          Long            Lat         ColorCode     055113             121.4698        40.8417     FF6600     057113             121.4698        40.8325     FF6600     .     .     .     063125             121.0344        40.4658     FFCC00     063127             121.0344        40.4567     FFCC00      I would like to have the output something similar like this:Example ![raster overlay](http://i.stack.imgur.com/I3Kye.jpg) Hence I need to write some code to create a raster image based on long,latt,color information. Do anyone know how to do this? or any library tool which can import into Java environment to do this? Method and direction is needed. Thank you!


title: How to specify fields to export from QGIS?

text: I'm using QGIS 2.0.1 and am connected to a MS SQL Server spatially enabled database (not ESRI). I am able to connect to the layer and load it into my TOC just fine. I have hidden a number of fields that are not important for visibility sake. Now I want to export that data to another data type (we'll use shapefile for this example). When I right-click the layer, choose "Save As" and define the output file, I get the following error: > Export to vector file failed. Error: unsupported type for field GlobalID _**Note 1** : If I check "Skip attribute creation" in the dialog box, the shapefile is created with no errors._ _**Note 2** : ArcGIS does not have this problem, it just converts the GUID field to "TEXT" type._ Looking at the table, one of the fields that I have hidden is a GUID of SQL type <uniqueidentifier>, which is in all likelyhood the culprit. In an attempt to resolve this, I made the field "hidden". I would have assumed that because the field was hidden, QGIS would not try to export this field, but that turned out to not be the case. **Can you define the fields you want to export from QGIS?** Is it possible for QGIS to export just the "Visible" fields? Or is there another alternative way of doing this from QGIS (maybe define the fields via an ogr2ogr query)? This seems like a shortcoming of QGIS if I cannot define which fields to export. **UPDATE:** I did a test using a sample shapefile to see if I could limit the output fields. The following ogr2ogr command line worked and eliminated all fields except for "name":               ogr2ogr -f "ESRI Shapefile" test.shp city.shp -sql "select name from city"      but using QGIS and adding:               -sql "select name from city"      to the "layer" section of the "Save As" dialog did not result in a shapefile with just the single field. So I don't really see how I can manipulate the file output to show just a select few fields. Any help is appreciated.


title: draw multipart polygon with one surrounding border

text: As you can see on the image below, the waterbody is divided by the border line because it consists of several polygons. ![enter image description here](http://i.stack.imgur.com/73vLm.png) I could merge the polygon, but I dont want to change the shapefile. Is there a option in QGIS to draw te waterbody just with a border sourrounding its "absolute extend"?


title: How to filter big datasets in Geoserver WMS?

text: Is there a good way to filter datasets with approx 50,000 features in Geoserver? I'd have to select between 100 and 200 features based on ids. These ids are different from the FeatureIds. Couldn't find a "IN" operator similar to SQL in OGC Filter Encoding Standard. Combining multiple PropertyIsEqualTo operators performs really poorly.


title: Cannot run the Arc samples

text: When I try to compile I get Error 1 Cannot register assembly "C:\Users\Patrick\Desktop\DynamicBiking\CSharp\bin\Debug\DynamicBiking.dll" - access denied. Please make sure you're running the application as administrator. Access to the registry key 'HKEY_CLASSES_ROOT\DynamicBiking.BikingTrackModeCmd' is denied. DynamicBiking2010 Trying to compile this one: http://help.arcgis.com/en/sdk/10.0/arcobjects_net/SampleZips/101efef5-9071-49ac-9ce9-c5a5fe5cc896CSharp.zip I have administrator priveleges


title: using shapely: translating between Polygons and MultiPolygons

text: [EDIT: the solution to this was simply to use OGR to read shapefiles. See geographika's example.] In an ESRI shapefile, there is no distinction between Polygons and MultiPolygons. Furthermore, there is no explicit distinction between interior holes and exterior rings (besides the "handedness" of a given polygon). So after reading a shapefile, I have a list of coordinate sequences describing rings, but without some more intensive processing, I cannot distinguish which of these rings are exterior rings, interior holes, or additional polygons. It appears that for shapely's Polygon and MultiPolygon constructors, there must be a clear distinction between exterior and interior rings, so how should I move from an unclear list of rings to an ordered set of separated polygons, with clearly designated interior and exterior rings? To summarize: if I have a list of polygon rings, but I don't know which rings are holes in the interior or are separate polygons, how should I best sort them into separate polygons with designated interior holes? I'm looking for a simple algorithmic solution that I can implement in python, can be used to process hundreds of polygons in ~a minute or less, and I'm doing this in order to perform a large number of intersections.


title: Counting Events Per Polygon Feature

text: I would like to calculate the total number of events within each feature of a polygon using a dbf file and a polygon shapefile. Is there an easy way to automate this so I can then produce charts per feature area? Thanks, Haskett


title: Out of range error when converting from geometry to geography in PostGIS

text: I am attempting to convert a geometry column in the postgres database to a geography column. I created a new column like this:               alter table scene add column footprint_new geography(POLYGON, 4326);      and I'm trying to populate the new column like so:               update scene set footprint_new = footprint::geography;      and I get the following error:               ERROR:  Coordinate values are out of range [-180 -90, 180 90] for GEOGRAPHY type      Now, it isn't too surprising that there are some out of range values in this database. I want to skip those out of range geometries during the update and identify them later so they can be fixed. How can I do that?


title: How to use more than one CPU?

text: I'm creating a real time traffic application. I read the AVL car data from a web service and insert the records in a PostgreSQL table Then a trigger calculates the closest link and if there is a previous position for the same car I calculate the route to get the distance travelled in that time frame, then I get the average speed for that route. The thing is I get around 3000 points each 5 min, so I need to process ~600 records every min or 10 each second. I'm very close to those numbers because I filter the road network before calling the routing function and the routes are very small. I understand PostgreSQL doesn't support multithread in one connection but support multiple connections so could use the other CPU cores. I can make a different table for each State (25), but how can I partition the problem so each trigger creates a new connection and use more CPU cores instead of only one?


title: Using primary keys with ArcSDE (Oracle)?

text: Is it possible, more advisable to assign a primary key to the OBJECTID column in a featureclass in SDE. I am using Oracle and some of my flat tables (junction tables) need to SEE a primary key in the featureclass table in order to model a M-M relationship correctly. By default ,no primary keys are made by Oracle when you insert and SDE layer (which makes sense as you are doing this all through ESRI tools anyway) although I know SDE creates an index in SDE, but I did not want to cause any kind of corruption if I start meddling with things on the Oracle side - is there any danger of this? Before submitting this question I also did a little more reading around the subject, and I seem to be finding that it is bad practice to use OBJECTID as a form of primmary key in any case - so would I be better off creating my own primary key field and then let Oracle maintain this with triggers and sequnces? If I do this, effectively behind SDE's back, will I upset SDE and all it's metadata etc? **UPDATE** I have found this nugget of info from Neil Clemmons - what do you think? > Using the OID may not be a good choice. For instance, let's say you have a > road centerline feature class related to a maintenance table using the OID > field. If you decide to merge two centerlines into the same feature, ArcMap > does this by deleting the two original features and creating a new one. > Depending on how you have the relate setup, the records in the related table > could be deleted. At the very least, the new feature would not be related to > any records in the table because it has a new OID. If you based the relate > on another unique field, the new feature would retain the attributes of one > of the original features and the relationship with the table would be > maintained. and more Neil Clemmons advice... > You shouldn't use the OID as a foreign key. This field is created and > handled by the ArcGIS software and is present so that a unique identifier is > always available for the ArcObjects functions that need one. You have no > control over the value of this field. If you need such a field for use as a > foreign key, you should create one and populate it yourself. This way, the > value will stay with a feature throughout the lifetime of that feature.


title: ArcObjects ITextSymbol.GetTextSize in Python revisited: usage with StartDrawing and FinishDrawing methods

text: I am using ArcGIS 10.0/ArcObjects and Python 2.6. I am trying to get the usage of the ITextSymbol.GetTextSize method correct. The syntax, as established in an earlier question, is:               dblx, dbly = pTextSymbol.GetTextSize(pDisplay.hDC, pDisplay.DisplayTransformation, txt)      and the outparameters are returned just fine. My problem is in figuring out how to use this method so that it does not screw up ArcMap's screen cache. ESRI documentation indicates that it should be surrounded by calls to IDisplay.StartDrawing and IDisplay.FinishDrawing so that the input hDC will be correct. (Indeed, calling GetTextSize without the start and finish drawing results in an invalid hDC.) The syntax that I am using for StartDrawing is as follows:               pDisplay.StartDrawing(0, esriDisplay.esriNoScreenCache)      I have tried both 0 (as discussed in documentation) and pDipsplay.hDC for the first parameter of StartDrawing, but the ArcMap screen still gets stuck in one place. (Panning just causes a redraw from cache right back to the same place.) Either input for the first parameter results in valid outparams but a stuck screen. The problem could also be with the esirNoScreenCache constant, I suppose, but I don't see how. I am not trying to draw anything to the screen at this point. I'm just operating on a TextSymbol, not a TextElement or the like. pDisplay and txt are both valid, in-process objects, so the problem, I think, is with the hDC or the screen cache. I'm at a loss to figure out which. Any insights are greatly welcomed.


title: Spatial Classification

text: I am not sure what to call my problem. From this image, I have a set of classified polygons (colored). I would like a way to classify the white polygons in the image. I don’t have any attributes to base the classification on other than proximity to a classified polygon. I have tried a naïve approach by using the near tool in Arc. This works fine enough but it becomes sensitive to outliers. For example, a set of unclassified polygons may be surrounded by polygons of category "A" but if there is a single “island” polygon of category "B" nearby, they’ll be classified the same as the island polygon. I am looking for a way to rank the likelihood that an unclassified polygon belongs to a particular category based on proximity to arbitrary k-polygons. For example, if there is a polygon of category "B" immediately adjacent to an unclassified polygon but the next five closest polygons are of category "A", I'd like the system to rank the unclassified polygon as category "A". I am wondering if there’s anything in R that would solve my problem. Really, any off the shelf method would work (python, qgis, postgis, arc). ![enter image description here](http://i.stack.imgur.com/ZBg9e.jpg)


title: Unable to Install External QGIS Plugins

text: I'm using: QGIS version 1.8.0-Lisboa on Windows XP, installed using OSGeo4W I am behind a proxy but I have configured it in Options > Network. When I go to Plugins > Fetch Python Plugins.. All of the the Repositories are "Status: Connected" however on the Plugins tab I can only see the 6 installed plugins (Repository: only locally available) - See images below: ![Repositories ](http://i.stack.imgur.com/l1qgS.png) ![Plugins](http://i.stack.imgur.com/Y67kx.png) Any ideas why I cannot see any plugins from the "connected" repositories? Many thanks. Update - As suggested by vinayan I can install plugins by manually downloading them and putting them in the relevant folder and using the QGIS Plugin Manager. However still no joy using the Fetch Python Plugins.


title: issue with WEB ADF based Web application

text: I have a WEB ADF based application that I am having some issues with. When I start app in debug mode, it shows me a blank page but I can right click and choose the option view in browser and it works fine and loads properly. I am trying to figure out what I need to do to load this page propely in a debug mode. Appreciate any ideas. Thanks Jay


title: Does anyone know a good way to handle z ordering when using GraphicTracker in ArcEngine?

text: From the GraphicTracker Documentation: > The first graphic added will appear underneath subsequent graphics if they > have overlapping display areas. This is not true when using ArcEngine 10.0 with polygon fill symbols on picture marker symbols - sometimes it works (1st time for sure) and sometimes it does not. Has anyone figured out a proper way to manage z ordering of symbols when using ArcEngine and using GraphicTracker?


title: Polygon shapefile with null values

text: Does anyone know where I can find a sample shapefile that has null entries in it (preferably a polygon shapefile)? I'm not sure how to create one with null entries. Thanks, Jeff


title: PostGIS - Intersects feature to split into linestring

text: I would like to know if it's possible to intersect feature datasets (multipolygon, polygon, multilinestring and linestring) to get a result as a set of linestrings being split at each intersection point. ![enter image description here](http://i.stack.imgur.com/kR4n5.png) Output should give:   * 4 linestrings for the square    * 3 linestrings for the line   * 2 linestrings for the polygon


title: Lost qgis_core.dll

text: I´ve been working with QGIS 2.4 since two months ago. After doing some changes in my computer, I think I've lost something important for running QGIS 2.4 The message is: > qgis.exe - System error. > > The program can´t start because qgis-core.dll is missing from your computer. > Try reinstalling the program to fix the problem I reinstalled the 2.4 32 and 64, the 2.2 version and the problem is still there. I can't open QGIS. It shows the first window and crashes. The QGIS browser works fine.... This is the path of my computer: > C:\Program > Files\Gimp-2.8\bin;C:\Python34\;C:\Python34\Scripts;C:\Python34\Lib\site- > packages;C:\Program Files (x86)\Intel\iCLS Client\;C:\Program > Files\Intel\iCLS > Client\;%SystemRoot%\system32;%SystemRoot%;%SystemRoot%\System32\Wbem;%SYSTEMROOT%\System32\WindowsPowerShell\v1.0\;C:\Program > Files\Hewlett-Packard\SimplePass\;C:\Program Files (x86)\Windows > Live\Shared;C:\Program Files (x86)\ATI Technologies\ATI.ACE\Core- > Static;C:\Program Files\Intel\Intel(R) Management Engine > Components\DAL;C:\Program Files\Intel\Intel(R) Management Engine > Components\IPT;C:\Program Files (x86)\Intel\Intel(R) Management Engine > Components\DAL;C:\Program Files (x86)\Intel\Intel(R) Management Engine > Components\IPT and the python path is: > C:\Python34\Lib;C:\Python34\DLLs;C:\Python34\Lib\lib- > tk;C:\Python34\Scripts;C:\Python34\Lib\site-packages;C:\Python34\projects I don´t know what to do to fix the problem and I can´t find anything on the web.


title: How to measure correct distances in a project with OpenLayers background?

text: I have 3 lat/long points as follows:               A = {35.257246,139.721256}     B = {35.256541,139.721491}     C = {35.256326,139.7203271}      * * * 1) Measurement result by Movable Type Scripts: By using the distance calculation method at Movable Type Scripts, the distance results are:               distance(A,B)= 80.88745771063131(m)     distance(B,C)= 101.58773310547808(m)      I believe these results are correct. * * * 2) Measurement result by QGIS measure line tool: The QGIS desktop 2.0.1 and the Projection CRS is set to `[+proj=longlat +datum=WGS84 +no_defs]`. I use the Plugins/OpenLayers plugin/Add Google Satelite and import 3 lat/long points into QGIS. The results are:               distance(A,B) = 100.374(m)     distance(B,C) = 135.114(m)      I need to generate a map from a large set of lat/long points as above and Google Satellite background on QGIS. Could you show me how to solve this problem?


title: Error during visualization of many layers on Chrome (ExtJS, GeoExt, OpenLayers)

text: I've noticed that on the latest version of Chrome (36.0.01985) my application crashes when a lot of layers are displayed (and/or the browser is maximized and the map is zoomed out). In particular, the panels are not shown at all and the map is only partially rendered. No matter the kind of layer (raster, vector), the problem is always present when there is a lot of data to render. The problem doesn't occur on other browsers. I've also tested it also on an older version of Chrome (26) and interestingly, there was no error. I'm actually using ExtJS 3.4, GeoExt 1.1 and OpenLayers 2.13.1 (I know they are oldies, but till now they worked very well). I'm also wondering if an upgrade of these tools could solve the problem or only make me feel younger. Tips?


title: Python Syntax Error when running Python Snippet of successful geoprocess

text: I am having an issue with a simple python snippet. Here's the background: I wanted to include an additional step in a nightly update script that replaces characters in a field that simply do not play nice with a web application that allows users to query the data. These characters included angle brackets ("<" and ">") which I was trying to replace with "Less Than" and "Greater Than". So when I was testing the Field Calculation, I finally got it to work the way I wanted to, and as I have done with so many other snippets, I simply copied the Python Snippet of the geoprocess, and pasted into my script. Below is the snippet:               arcpy.CalculateField_management("path/to/my/layer","field","RemoveBad( !field!)","PYTHON_9.3","def RemoveBad(x):/n   x = x.replace(\"<=\", \" Less Than or Equal To \" )/n   x = x.replace(\">=\", \" Greater Than or Equal To \")/n   x = x.replace(\"<\", \" Less Than \")/n   x = x.replace(\">\", \" Greater Than \")/n   x = x.replace(\"=\", \" Equal To \")/n   x = x.strip()/n   return x")      Now in my script, I do have a log that reports geoprocessing errors, so when I run the script, I am getting the following:               Executing: CalculateField path\to\my\layer field RemoveBad( !field!) PYTHON_9.3 "def RemoveBad(x):/n   x = x.replace('<=', ' Less Than or Equal To ' )/n   x = x.replace('>=', ' Greater Than or Equal To ')/n   x = x.replace('<', ' Less Than ')/n   x = x.replace('>', ' Greater Than ')/n   x = x.replace('=', ' Equal To ')/n   x = x.strip()/n   return x"     Start Time: Mon Aug 25 09:35:32 2014     Failed to execute. Parameters are not valid.     ERROR 000989: Python syntax error: Parsing error SyntaxError: invalid syntax (line 1)     Failed to execute (CalculateField).     Failed at Mon Aug 25 09:35:32 2014 (Elapsed Time: 0.00 seconds)      So I started combing through the script but I cannot find the issue, which leads me to believe it has something to do with the angle brackets. To reiterate, this python snippet is from a successful Field Calculator geoprocess. Any assistance in pointing out where the syntax error is occurring would be appreciated. **UPDATE** : Per Alex's feedback I created the following variable:               RemoveBad = """def RemoveBad(x):        x = x.replace(\"<=\", \" Less Than or Equal To \" )        x = x.replace(\">=\", \" Greater Than or Equal To \")        x = x.replace(\"<\", \" Less Than \")        x = x.replace(\">\", \" Greater Than \")        x = x.replace(\"=\", \" Equal To \")        x = x.strip()        return x"""      Then I updated the snippet:                arcpy.CalculateField_management("path/to/my/layer","field","RemoveBad( !field!)","PYTHON_9.3",RemoveBad)      The error is now gone and script is working as expected.


title: Using auto increment code with counter

text: I was wondering if it's possible to use the standard auto increment code to work on groups of records:               FeatureClass = "newFC"      recs = collections.defaultdict(int) autoIncrementEXP = "autoIncrement()"      expCalc = """def autoIncrement(!Group!): recs[!Group!] += 1 return recs[!Group!]"""     arcpy.CalculateField_management(FeatureClass, "Number",autoIncrementEXP,"PYTHON",expCalc)      So instead of having the auto increment portion of the script stop when it reaches the end of the records, have it create sequential numbers using field value groups if that makes sense. You would start with something like the `Group` field below and end with the `Number` field containing the increment values:               Group   Number     A       1     A       2     A       3     B       1     B       2     C       1     C       2     C       3     C       4     C       5      Thanks for any insight you can provide.


title: Can ArcGIS Desktop 10 Add-in Use IronPython?

text: Coming from the same background as the original poster in this thread, I wonder if anyone has had any success creating an Arc10 Add-In using IronPython. I have been able to access ArcObjects via IronPython using the info in the linked thread, assorted blog postings, and a few posts on the ESRI forums. These neophyte attempts have all been run from the command line and I would to be able to gain access through either a toolbar or a button. One thought that I had was to create the add-in using the .Net 3.5 template, add a python file, and then call a function in the python file from the button's C# file. It appears though that add-ins are not yet supported in .Net4? Here is the code from the .cs file which controls button functionality.               using System;     using System.Collections.Generic;     using System.Text;     using System.IO;          //I need to then add references to IronPython     using IronPython;     using IronPython.Hosting;     using Microsoft.Scripting;     using Microsoft.Scripting.Hosting;          namespace PyButtonTest     {         public class ZoomToLayer : ESRI.ArcGIS.Desktop.AddIns.Button         {             public ZoomToLayer()             {             }                  protected override void OnClick()             {                 var ipy = Python.CreateRuntime();                 dynamic test = ipy.UseFile("ZoomToLayer.py");                 test.ZoomToActiveLayerInTOC(); //Call this function in the py file                      ArcMap.Application.CurrentTool = null;             }             protected override void OnUpdate()             {                 Enabled = ArcMap.Application != null;             }         }          }      As far as I can tell the Add-in format takes care of getting the application hook. I should then be able to, in the python file, after referencing and importing all of the necessary libraries get the map document with (I think?)               mxDoc = ArcMapUI.IMxDocument.ActiveView      From here I am stumped. Google has not turned up anything helpful. Has anyone else tried this? Got something to work? Do I need to wait until add-ins are supported in .Net4 to go this route?


title: QGIS Contour plugin: Creating complex isolines leads to isolated areas

text: I use the contour plugin (v1.1) within QGIS 2.4 to create contour lines. I have streets with different travel costs. I observe problems for streets with low travel costs if they are surrounded by streets with high travel costs. The low cost street is "piercing" into the low speed area (cf. picture). The points are created as the travel cost markers from Dijkstra algorithm using pgrouting 2.0.5. The Dijkstra algo. calculates correct values for every street. The problem I observe is now, that due to low speed areas left and right to the highway, the isolines are connecting both areas over the highway (yellow). Based on values, the correct behavior would be the one shown in green. **Can someone recommend me a way of preventing such results?** Before I used GDAL TIN interpolation and contour mapping, which leads to better results if it comes to complex networks. On the other hand it adds transformation errors (rasterizing) and extra steps to the calculations. ![Section of my map showing two 60 sec. isolines without connection. Green: The likely correct isoline. Yellow: The wrong connections for the isolines.](http://i.stack.imgur.com/iyVhw.png) ![Original result](http://i.stack.imgur.com/x0yki.png)


title: Counts the number of points in a polygon in R

text: I have two classes sharing the same CRS (Latitute and Longitude):   1. `bolognaQuartieriMap`: a `SpatialPolygonDataFrame` containing data of a city boroughts.   2. `crashPoints`: a `SpatialPointsDataFrame` containing data of accidents. They are well plotted using:               plot(bolognaQuartieriMap)     title("Crash per quartiere")     plot(crashPoints, col="red",add=TRUE)      What I need is to get the number of points (`crashPoints`) in each polygon that constitute `bolognaQuartieriMap`. I was suggested to use `over()` but I did not succeed.


title: OpenLayers: Switch layer off when switching another layer on

text: As the title says, I want to turn layers off when turning another layer on. I wrote a function that works. However, I have to click the layer twice in the LayerSwitcher to turn the layer on that I want to see. The first click turns the other layers off but only a second click turns the actual layer tho show on. Can somebody tell me where I am wrong in my code?                   function getVLayers(showLayer){             for (var i=0; i < map.layers.length; i += 1){                 if (map.layers[i].visibility){                     if (map.layers[i].name == showLayer && map.layers[i].name == "OpenStreetMap"){                          } else {                         map.layers[i].setVisibility(false);                     }                 }             }         };         map.events.register("changelayer", this, function(e){             getVLayers(e.layer.name);         });


title: QGIS GRASS Missing Dependancy

text: I want to perform a view shed analysis in QGIS using GRASS 7 r.viewshed. When I double click on r.viewshed in GRASS GIS 7 commands in the processing toolbox of QGIS I get the following message: > _Missing dependency.This algorithm cannot be run :-( This algorithm requires > GRASS GIS 7 to be run. Unfortunately, it seems that GRASS GIS 7 is not > installed in your system, or it is not correctly configured to be used from > QGIS Click here to know more about how to install and configure GRASS GIS 7 > to be used with QGIS_ GRASS commands seem to work such as r.los but the GRASS 7 commands don't. I am using a MAC OS X Mavericks with the latest version of QGIS 2.4.0 Chugiak.


title: raster2pgsql does nothing

text: I try to load rasters through raster2pgsql. I try to call it through pgAdmin command line, but when I call any `raster2pgsql` command (even `raster2pgsql -G`), it does nothing, I just get another line to write more commands. No output, no time doing something. Non-gdal psql commands work as expected. `\raster2pgsql` leads to "invalid command" error. I am new to the command line and psql, so I'm stuck here. I have PostgreSQL 9.3.4, PostGIS 2.1.3 and pgAdmin III 1.18.1 on WIndows XP 32 bit. I made sure that `raster2pgsql.exe` is in its place in `bin` folder, but I don't know of other ways to validate whether it is installed properly.


title: How to include data from other tables in the info-window in CartoDB

text: I have 3 layers(table_1,table_2,table_3) of data. In the info-window, which was created for the second layer. I would like to display one column of data from table_1 . When I used `{{table_1.name}}`to display the data, I am getting the error , `"Error in line 40: Cannot read property 'table_1' of undefined` Is there any way to display data from one table in another layer of data?


title: Formatting Cells in ArcMap

text: I am trying to join tables together and need the PLSS to line up correctly. In my table with geometries I have "CO060020N0650W0SN330" in a field. I want to make a new column and make it say 2N65W33 so I can join a table. Is there a formula I can use in the Field Calculator to get this result?


title: When styling with data defined properties or Labeling, Is there a way to replace an ID representing field with the actual values depicted by the ID?

text: I am currently using data defined properties and labeling to present pipes based on diameter. I now have a dictionary table that has id's for every combination of diameter and material and I'm using it with the value relation widget. I must find a way to access the value represented by the key in the diameter || material field. An ugly way that solves this would be to use:               CASE        WHEN id_material_diam in (id1, id2) THEN result_formula_containing_diameter        WHEN id_material_diam diameter in (id3...) THEN ...     END      Is there an elegant way that I can access the actual values for which the ID is the key? Thank you PS. QGIS 2.2 Valmiera


title: GDAL equivalent of ArcGIS Mosaic

text: The word 'mosaic' is often used for quite different things within the worlds of ArcGIS and GDAL. With GDAL and it's downstream applications like Qgis it is common to see mosaic used in the context of _"assemble all these tiled images so they appear as one big image"_ , with gdalbuildvrt and `gdalmerge` put forth as a convenient means of getting there. ArcGIS has the same starting point, _"assemble all these tiled images so they appear as one big image"_ , and then goes on to add _"while doing so [blend, min, max, mean, sum, ...] the overlapping pixels"_ and a few other things. To put in uncharitably, in this regard ArcGIS is smart and GDAL is stupid. However it's pretty evident GDAL is _not_ stupid, so my question is: in the GDAL and/or Qgis tool chain what is the equivalent workflow to ArcGIS mosaic? _(note, I'm not talking about the Mosaic Dataset, which is a different beast altogether)_


title: R - Dealing with missing data in SpatialPolygonDataFrames for Moran.test

text: I'm trying to use Moran.test on a SpatialPolygonDataFrame consisting of 7194 elements in R. I know that there is around 150 polygons with NA values. First I generate a spatial weights matrix: `W_cont_msoa <- poly2nb(all.data, queen=T) W_cont_msoa_mat <- nb2listw(W_cont_msoa, style="W", zero.policy=TRUE)` Then when I carry out Moran.test the code below returns:               > moran.test(residuals, listw=W_cont_msoa_mat, zero.policy=T)     Error in moran.test(residuals, listw = W_cont_msoa_mat, zero.policy = T) :      objects of different length      I've also tried to this but it hasn't worked:               > moran.test(residuals, listw=W_cont_msoa_mat, zero.policy=T, na.rm=TRUE)     Error in moran.test(residuals, listw = W_cont_msoa_mat, zero.policy = T,  :      unused argument(s) (na.rm = TRUE)      I've tried to go back to the original data and remove missing data using complete cases but I received the following:               > all.data[complete.cases(all.data)]     Error in all.data[complete.cases(all.data)] :      error in evaluating the argument 'i' in selecting a method for function '[': Error in       complete.cases(all.data) : invalid 'type' (S4) of argument      The structure of the data is as follows:               `Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots     ..@ data       :'data.frame':   7194 obs. of  23 variables:     .. ..$ MSOA01CD   : Factor w/ 7194 levels "E02000001","E02000002",..: 1 2 3 4 5 6 7 8    9 10 ...        .. ..$ MSOA01NM   : Factor w/ 7194 levels "Adur 001","Adur 002",..: 1468 125 126 127     128 129 130 131 132 133 ...          .. ..$ averagesp  : num [1:7194] 11.6 20.9 20.3 21.4 19.1 ...        .. ..$ nga        : int [1:7194] 41 112 177 97 120 NA 142 161 122 113 ...        .. ..$ mediansp   : num [1:7194] 11.5 13.1 12.6 13.5 10 ...          .. ..$ mbps2      : int [1:7194] 125 74 104 48 83 NA 64 67 76 25 ...         .. ..$ msoa.total : int [1:7194] 2372 114 191 97 126 NA 157 208 130 119 ...          .. ..$ p.nga      : num [1:7194] 1.73 98.25 92.67 100 95.24 ...          .. ..$ p.mbps2    : num [1:7194] 5.27 64.91 54.45 49.48 65.87 ...        .. ..$ residents  : Factor w/ 4057 levels "10,003","10,004",..: 2222 1748 27 1261       3216 NA 3380 403 3093 3256 ...        .. ..$ area       : Factor w/ 6908 levels "1,002.23","1,002.58",..: 3707 2869 2840      3304 839 NA 1899 1327 6811 748 ...          .. ..$ density    : num [1:7194] 25.5 31.3 46.9 24.8 72.1 NA 50.6 81.5 87.4 76.8 ...          .. ..$ micro.bc   : num [1:7194] 36 32 83 33 39 NA 74 60 44 32 ...           .. ..$ small.bc   : num [1:7194] 49 14 64 3 14 NA 50 14 14 73 ...        .. ..$ medium.bc  : num [1:7194] 850 0 15 0 0 NA 10 0 0 0 ...        .. ..$ large.bc   : num [1:7194] 245 0 0 0 0 NA 0 0 0 0 ...           .. ..$ primary    : num [1:7194] 60 0 0 0 0 NA 0 0 0 0 ...          .. ..$ secondary  : num [1:7194] 135 0 20 0 0 NA 0 0 5 5 ...         .. ..$ tertiary   : num [1:7194] 16600 95 335 100 145 NA 300 200 160 80 ...          .. ..$ kibs       : num [1:7194] 11805 20 80 20 30 ...           .. ..$ edu        : num [1:7194] 2205 2681 811 2658 116 ...          .. ..$ p.non.white: num [1:7194] 21.4 35 45.4 19 33.7 NA 32.8 30.4 38.8 30.8 ...        .. ..$ claimants  : num [1:7194] 117 245 316 254 295 ...         ..@ polygons   :List of 7194         **Can anyone suggest how I overcome this problem please?** Thanks, Ed


title: ArcGIS Publish query layer as Feature Service fails

text: I have a query layer that's simply a DB view call in PostgreSQL, the view is as such: ` CREATE OR REPLACE VIEW daily_query AS SELECT b.id, b.first_name, b.last_name, c.assignment_date::date AS assignment_date, d.case_id, d.case_sequence, d.status, a.address_complete, a.the_geom_point FROM universe a, staff b, enumerator_assignments c, enumerator_case_assignments d WHERE b.id = c.enumerator_id AND c.objectid = d.enumerator_assignments_id AND a.id = d.case_id;` The query layer is simply a call to the view `select * from daily_query` The layer renders fine in ArcMap - but then I try to publish the layer as a Service to ArcGIS Server, I want this to be a FeatureService so I enable "Feature Access" Capability and I only allow "Query" operation. and I don't allow "Geometry Updates". When I run "Analyze" it returns with "00134: Layer's data source is not supported" The Service Analyzes and Publishes fine if I don't enable "Feature Access" on it. I'm running 10.2.2 Any ideas?


title: How do I change the targeted .NET 4.0 Framework to version 3.5 in an existing VS 2010 EXPRESS project?

text: I am upgrading from ArcGIS desktop 10.0 to 10.1 and I have a bunch of addins developed in VS 2008 Express using ArcObjects 10.0. I would like to upgrade everything **to 10.1 using VS 2010 Express**. I understand that for ArcGIS Desktop **.NET Framework 3.5 Service Pack 1 is still required.** I converted my first existing addin to 10.1 in VS 2010 Express but it seems to be targeting .Net 4.0. I only found instructions on changing the target version in full version of VS but I can't find this feature in the Express edition. Can someone shed some light on this? EDIT: I am adding a screenshot from Project -> Project Properties -> Application to show where I would expect to change the Target Framework but I don't see where it could be changed. ![enter image description here](http://i.stack.imgur.com/C9nwd.jpg) **EDIT 2:** I found an option to change Target Framework (All Configurations) under **Project -> Project Properties -> Compile -> Advanced Compile Options -> Target Framework (All Configurations)** Can someone confirm that his is where I change the target framework in the 2010 Express edition?


title: MapInfo License Server - Can It Serve Multiple Versions of MapInfo

text: I have a MapInfo License Server Utility that serves 30+ MI v10.5 concurrent licenses across my organisation. v11.5 has now landed on my desk (with new serial numbers and access codes), so can I activate my v11.5 licenses within the same instance of License Server utility, whilst retaining the v10.5 licenses, in order to gradually migrate users onto v11.5 over the coming months?


title: How to calculate the shortest path between two points characterized by latitude and longitude with pgrouting?

text: I have to calculate the **length** of the shortest path between A (lat1, long1) and B (lat2, long2). But all the solutions I've found do not get a correct result because they only calculate the distance between the nearest vertex of A to the nearest vertex of B.


title: How can I select features by map clicks

text: I'm developing a flex web application using the ArcGIS flex API 3.1. I have a map with a polyline feature layer. I'd like to be able to select polylines from the map by clicking on the features on the map. Could you guide me to a sample? Regards.


title: Algorithm for finding irrregular polygon centroid (label point)

text: I need to find a centroid (or label point) for irregularly shaped polygons in Google Maps. I'm showing InfoWindows for parcels and need a place to anchor the the InfoWindow that is guaranteed to be on the surface. See images below. ![alt text](http://i.stack.imgur.com/sGcEX.png) ![alt text](http://i.stack.imgur.com/pQDX0.png) In reality I don't need anything Google Maps specific, just looking for an idea of how to automatically find this point. My first idea was to find the "false" centroid by taking the average lat and lngs and the randomly placing points out from there until I find one that intersects the polygon. I already have the point-in-polygon code. This just seems awfully "hacky" to me. I should note that I don't have access to any of the server side code outputting the geometry so I can't do anything like ST_PointOnSurface(the_geom).


title: How to display only one pie chart for polygon features in QGIS?

text: I'm working with Qgis on a project and I have a small problem. I have to do some demographic research and therefore I have a .shp file from a given country with 10 districts for example. Now, I know how to insert pie charts but the problem is that I only want to insert one pie chart for one district because this area is of importance to me for example. So I added a column in the attribute table with data for the district that I want to display and left all the other values 'NULL'. However, if I insert pie charts, they appear for all the districts, even though I've not inserted any values in the attribute table for the other 9 districts that are of no interest to me. Does anybody know how I can only display one pie chart for 1 feature of the polygon instead of having 1 good pie chart and 9 empty on the map?


title: QGIS Cross Tab of Towns and Routes

text: Using QGIS, I have a layer with 100 towns. I have another layer with 200 cycle routes, each with a different name. I want to create (automatically natch) a table that will help the people in each town, for example, get a list of all of the routes that pass within 1Km of their town centre.


title: Understanding web mapping tools

text: Before I start, this is going to be a general question but with specific cases. I've been trying to understand GIS tools and applications to achieve the following:   1. Build my own basemap   2. Publish the basemap online   3. Add my own data layers to the online basemap   4. Perform analysis (on-the-fly)    5. Add interactivity to the map So, In my first attempt, I built a spatially enabled database on PostGIS, queried this using PHP, overlayed it on google maps using their API and added interactivity using their api itself. I carried out analysis like area, denisty etc. on the fly using postGIS functions. For my next attempt, I downloaded TileMill and rendered some really beautiful maps and added zoom level functions to it. I also exported my MBtile files. When I came across trying to publish this online, I hit quite a few names, but got lost and confused:   1. TileStache   2. TileCache   3. TileStream   4. Mapserver   5. Mapnik   6. GeoServer   7. Mapbox So, I understand these are mapservers. But what exactly does a mapserver do? Why is it used? Why can't i server my tiles from a normal web Apache server? And for MBtiles, how many servers are out there which lets you publish it directly? I can also export my basemap as a png/vsvg format. Can these formats be used to make a basemap through a mapserver? When researching interactivity I came across:   1. Leaflet   2. ModestMaps   3. OpenLayers   4. Mapstraction I understand these are Javascript Mapping API's. But on what base map can I use these? Can it be used on Google maps? Can I use them on my own base maps too? I basically want to understand the differences between these tools and which would be the best software stack to build a fully customized, interactive map with my own basemap. Sorry if my question sounds to general. I will definitely be here to clarify my questions


title: Change Vector projection in Openlayers

text: I have added vector layer to map                   vactorTest = new OpenLayers.Layer.Vector("Editable Features", {     strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],     projection: new OpenLayers.Projection("EPSG:900913"),     protocol: new OpenLayers.Protocol.WFS({         srsName: "EPSG:900913",     url: "http://localhost:8080/geoserver/wfs",     featureType: "ROADS",     featureNS: "http://localhost:8080/MIA",     featurePrefix: "MIA",     geometryName: "GEOM",     version: "1.1.0"     })     });          map.addLayer(vactorTest);      And than added polygon draw control               var drawPolygon = new OpenLayers.Control.DrawFeature(     vactorTest, OpenLayers.Handler.Polygon,     {     title: "Draw Polygon",     displayClass: "olControlDrawFeaturePolygon",     multi: true     }     );      with drawn polygon I want to filter layer               testLayer.mergeNewParams({         "CQL_FILTER": mergeFilter     });      I use within filter: WITHIN(GEOM, POLYGON((4693721.49445 672633.89194, 4560244.10812 322255.75282, 4703732.29843 422363.79257, 4693721.49445 672633.89194))) target layer is in Georgia, and if I draw polygon near Italy than filter works.... How can I change vector projection like my layer to get right result?


title: Is there a way to figure out what files an MXD points to?

text: Let's say I load up two layers, and then save an .mxd file. Now I want to get the path names of my layers from the .mxd file. Is there a way to do this in ArcEngine? I'm using ArcEngine 10 with VS2010 & C#