title: CartoDb - Aggregate and Infowindow

text: I'd like to use CartoDb to map historic addresses of a city (link to demo). Individual addresses are mapped as points on the map. Detailed information for each address point should be displayed in the infowindow The geocoded addresses are stored like this:               NAME   PROFESSION ADDRESS ... (+ CartoDb fields like cartodb_id, the_geom_webmercator)     Smith  shopkeeker grand-rue 10     Martin baker      grand-rue 11     Wagner grocery    grand-rue 11      As long as there is only one person living at an address, I get the desired infowindow by simply doing a _select_ * _from_ _table_name_. When two or more persons are associated with the same address I obviously would need some aggregate function in PostgreSQL. I am experimenting with:               SELECT  address, string_agg(name || ', ' || profession, '; ' ORDER BY name, profession) FROM table_name group by address      I am running into two problems:   1. In order to map, CartoDb expects the columns 'cartodb_id' and 'the_geo_webmercator' to be included in the SELECT statement. When doing so, the result of the query is not grouped by 'address' in the desired way; instead the unique 'cartodb_id' is used for grouping. So, the following is not grouping according to 'address':              SELECT  address, string_agg(name || ', ' || profession, '; ' ORDER BY name, profession), the_geo_webmercator, cartodb_id FROM table_name group by address, the_geo_webmercator, cartodb_id        2. The output of 'string_agg' is treated as one field in the design for the infowindow.              Grand-rue 11     Martin, baker; Wagner, grocery      What I would like to have, though, is a separate treatment of 'name' and 'profession' in the infowindow, like:               Grand-rue 11     1. Martin        baker     2. Wagner        grocery      Any hint is appreciated. Thanks! Peter


title: .dbf assocaited with aerial photos

text: I have aerial photos that are .jp2 format and are compiled into a .dbf. In ArcGIS I can simple add the .dbf file and all the aerials automatically appear. Can this be done in QGIS, rather than insert each individual .jp2 file?


title: Autodesk Infrastructure Map Server - WFS 1.1.0 GetFeature Request error

text: This is a duplicate request posted on the AutoDesk MapGuide Enterprise Discussion Group. If I take a working request...               <wfs:GetFeature service="WFS" version="1.1.0" maxFeatures="10"     xmlns:wfs="http://www.opengis.net/wfs">       <wfs:Query typeName="Water:WA_LARGE_VALVE" SRSNAME="EPSG:26910">         <Filter xmlns="http://www.opengis.net/ogc">           <Intersects xmlns="http://www.opengis.net/gml">             <PropertyName>Geometry</PropertyName>             <Envelope srsName="EPSG:26910">               <lowerCorner>491784.0 5456759.0</lowerCorner>               <upperCorner>491809.0 5456805.0</upperCorner>             </Envelope>           </Intersects>         </Filter>       </wfs:Query>     </wfs:GetFeature>      And replace Envelope with Polygon it fails...               <wfs:GetFeature service="WFS" version="1.1.0" maxFeatures="1"     xmlns:wfs="http://www.opengis.net/wfs">       <wfs:Query typeName="Water:WA_LARGE_VALVE" SRSNAME="EPSG:26910">         <Filter xmlns="http://www.opengis.net/ogc">             <Intersects xmlns="http://www.opengis.net/gml">             <PropertyName>Geometry</PropertyName>             <Polygon srsName="EPSG:26910">               <exterior>                 <LinearRing>                   <posList>489093.629 5456057.158 489092.66826402                   5456047.4034839 489089.822976626 5456038.02382838                   489085.202480615 5456029.37948835 489078.984339059                   5456021.80266094 489071.407511651 5456015.58451938                   489062.763171618 5456010.96402337 489053.383516101                   5456008.11873598 489043.629 5456007.158 489033.874483899                   5456008.11873598 489024.494828382 5456010.96402337                   489015.850488349 5456015.58451938 489008.273660941                   5456021.80266094 489002.055519385 5456029.37948835                   488997.435023374 5456038.02382838 488994.58973598 5456047.4034839                   488993.629 5456057.158 488994.58973598 5456066.9125161                   488997.435023374 5456076.29217162 489002.055519385                   5456084.93651165 489008.273660941 5456092.51333906                   489015.850488349 5456098.73148062 489024.494828382                   5456103.35197663 489033.874483899 5456106.19726402 489043.629                   5456107.158 489053.383516101 5456106.19726402 489062.763171618                   5456103.35197663 489071.407511651 5456098.73148062                   489078.984339059 5456092.51333906 489085.202480615                   5456084.93651165 489089.822976626 5456076.29217162                   489092.66826402 5456066.9125161 489093.629                   5456057.158</posList>                 </LinearRing>               </exterior>             </Polygon>             </Intersects>         </Filter>       </wfs:Query>     </wfs:GetFeature>      I receive is this error:               An exception occurred in FDO component.         String does not represent a valid filter      Is there something wrong with this request or maybe both requests?


title: Selecting lines occupying identical positions in MapInfo or QGIS

text: I have two shapefiles from a client, one is their entire network (complete with unique IDs) and another a selection I need. However the selection file does not have unique IDs but the entire one does. How would I be able to select the lines that appear in both files from the entire network shapefile? Thanks, Alan


title: Openlayers -- Is there a way to delay server requests on zoom for a particular layer?

text: I have a rather complex issue where I have 10 TMS Layers (I'm using OpenLayers) that animate. When I zoom in and out quickly I can quickly rack up 1000's of hits to the same server. Is there a way I can delay the request of that layer until zooming has stopped for x amount of time? Thanks!


title: Spatial Data with Unity

text: I was wondering if anyone has used spatial data as an input source to build 3D environments and interactive applications with Unity? The geospatial 3D world is very young, and some what immature, while the gaming industry is much more robust and developed. Unity is very developed and provides a lot of great options for building 3D environments. However, the typical data sources used in the geospatial world (vectors, imagery, DEMs) are foreign to Unity. I have seen some round about ways to use this data in a Unity project, but I was wondering if anyone has experience or success with this process, and would be willing to share. Thanks...


title: check if mapquest geocoding return no values

text: I use mapquest's geocoding api. I want to check if the request is sucessfull, but returned no data. I type `vdfsbvdf54vdfd` (just a silly string) in a form field and nothing happens. I expect an alert like so `Sorry, wrong input`. Alert never happens This is my code snippet               $.ajax({         type: "POST",         url: "http://open.mapquestapi.com/geocoding/v1/address?key=012",         data: { location: ad, maxResults: 1}         })               .done(function( response ) {             var geoclng = response.results[0].locations[0].latLng.lng;             var geoclat = response.results[0].locations[0].latLng.lat;                  if (geoclng=="") {alert("Sorry, wrong input");}                   //create magic using lon/lat                               )}      If it helps, when I alert the `response` I get `object Object`. I dont know if its the API or I handle wrong the request's response Thanks in advance


title: Refering to 0° of longitude/latitude in a range

text: This is more of a question in diction... In scientific writing, is there a preferable way to write a range, bounded by 0° longitude, in a way similar to "5° W to 10° E" if there is no cardinal direction for `0`. The two ways that seem most reasonable to me are "0 to 5 degrees longitude" or "0° longitude to 10° E". Thanks!


title: Rebuilding Indexes in SDE Database throws "Could not get a lock." error

text: While trying to do some maintenance on our SDE database but i am always getting a "Could not get a lock." error on the details section of the process. Do i need to setup something for me to successfully run the index rebuilding tool? Btw, i am using ArcMap 10.1 to rebuild the indexes. Thanks.


title: Make opengeo suit online,Need some free web server

text: I have installed opengeo suit and it works fine on localhost,How can I make it online?Which kind of web hosting do I need?Do I need any server or a linux webhost is ok? Is there any documention that I can use it to learn how to do this? Do you know some free web servers to test opengeo on it? Thanks a lot for your help


title: Software needed to cross-reference geocodes to a Shapefile

text: I do not have any GIS software except Microsoft MapPoint. I have State district KML and Shapefiles of the districts. I am assuming these are for ArcGIS. I do not want to draw maps. Instead, I am looking for a program that, if I give it one of 500,000 geocodes, it will tell me from a database perspective which District it is in. I am looking for analytical information not visual.


title: Flow Mapping Lines by Magnitude

text: ![enter image description here](http://i.stack.imgur.com/dZ7V3.png)I am attempting to map the flow of remittances into Japan using the FlowMapper plug in in QGIS. When I use the Style-->Advanced--->Size Scale field--->magnitude command the lines become so thick and varied that they are impossible to distinguish from one another (in fact my map becomes one giant block of "lines"). Does anyone have any advice on how to fix this other than manually changing the size of lines classified under equal intervals? Is there a way to label the flow lines by country of origin so that the text is in and oriented with the flow line, not located on top (that is, once I'm able to get the flow lines to their desired thickness? How do I make the lines curved so they are more aesthetically pleasing? Thank you!


title: Cannot merge my map on CartoDB

text: I am trying to create a map on CartoDB using this kml file   http://cdn.journalism.cuny.edu/mapFiles/NYC_Community_Districts.kml and a CSV file. I have uploaded both to CartoDB but there the name column comes up empty so I cannot merge the two. The CSV file contains a community district column that I want to merge the map with but can't. Any ideas?


title: Batch merge large numbers of.asc grid files

text: I have >2.5k .asc grid files which I need to merge into a single overarching dataset. I would like to avoid having to convert them into esri grids first so was looking to use QGIS (v.2.0). I am aware of the command line length limitation with the QGIS merge tool so was looking to use the OSGeo4W shell to automate the process. Can someone help me with some code to cycle through the files and merge them? I was looking at this thread but the suggested solution doesn't work - it results in an "index out of range error" The files are all of the same format, in a single folder and have the name convention SDXXX.asc - where X is a single digit.


title: arcpy - ascii encoding error when updating field content in SQL Server table

text: I know there are plenty of posts and blog explanations about unicode errors, but I still can't figure out to handle it in my particular case. So here is my problem: I am writing a Python script to update records in an SQL Server table, with field content from a shapefile that has been edited in ArcPad. I use an `arcpy.da.UpdateCursor` to update existing records. I've put `# -*- coding: cp1252 -*-` on top of my script as I work with French characters. Just doing so, I get the ascii encoding error: > 'ascii' codec can't encode characters in position 0-1: ordinal not in > range(128) I've tried putting `u` before the field content:               ...     with arcpy.da.UpdateCursor(DBtable, DBFields, where_clause) as DBCur:         for DBrow in DBCur:             ...             DBrow[8] = u"{} - {}: {}".format(DBrow[8], date , AProw[8]) # AProw comes from a SearchCursor reading the shapefile.      I can then print a message with the text content, and it doesn't return an error message, but all characters are replaced (I get something like `???4???5????????>???5???` in my database field). Any accentuated or punctuation character can be present in the fields (there's a free text comment field) so I don't want to check for every possible non ascii character and replace it. I work with ArcGIS 10.2.2, the shapefile is edited in ArcPad 10.2 and the database is SQL Server 2008 R2. What am I missing? Thanks for any help. **EDIT: This only occurs if the destination database is SQL Server. No problem with a file gdb.** I have to add that the SQL Server table already contains non ascii characters.


title: How to replicate ArcGIS Intersect in PostGIS

text: I'm trying to replicate this ArcGIS process in PostGIS: http://blogs.esri.com/esri/arcgis/2012/11/13/spaghetti_and_meatballs/. This describes how to break buffered points down into polygons based on their intersections, counting the number of layers, and attributing that to the polygons in order to classify them. I'm using it to create a rough point density map with vectors, and the results were surprisingly nice for my data set in ArcGIS. However, I am struggling to come up with something workable in PostGIS where I need it for producing dynamic point density layers for a web map. In ArcGIS, I simply ran the Intersect tool on my buffered points layer to create the shapes I needed. In PostGIS, I ran this query: `CREATE TABLE buffer_table AS SELECT a.gid AS gid, ST_Buffer(a.geo,.003) AS geo FROM public.pointTable a;` `CREATE TABLE intersections AS SELECT a.gid AS gid_a, b.gid AS gid_b, ST_Intersection(a.geo,b.geo) AS geo FROM public.pointTable a, public.pointTable b WHERE ST_Intersects(a.geo, b.geo) AND a.gid < b.gid;` `DELETE FROM intersections WHERE id_a = id_b;` The output looks pretty much identical to the ArcGIS output, except that it is not breaking the polygons down to the same extent that is required for a meaningful density map. Here are screenshots of what I mean: ![ArcGIS](http://i.imgur.com/JDY0KO5.jpg) ![PostGIS](http://i.imgur.com/MBRj6dP.jpg) ArcGIS is on the left, and PostGIS is on the right. It is slightly difficult to tell, but the ArcGIS image shows the 'interior' polygon created where all 3 buffers intersect. The PostGIS output, on the other hand, does not create that interior polygon and instead it keeps its components intact. This makes it impossible to provide a classification for just that interior area with 3 layers on top of each other compared to just 1 for the outer parts. Does anyone know of any PostGIS function to break the polygon down to the extent I need? Alternatively, does anyone know of a better way to produce a point density map with vectors in PostGIS?


title: Transfering a lot of features from geodatabases

text: I'm having a problem where I need to export a large amount of content from 2 geodatabases of different nature (SDE -> File). I have three options: Spatial Select, a regular SQL select and the full feature class/table. What are the best practices to handle these large transfers? I'm having issues with the following code:                   public void ExportAll(IGisEntity entity)         {             ITable table = _Context.GetEntityTable(entity.EntityName);             IFeatureWorkspace workspace = _OutWorkspace as IFeatureWorkspace;             ITable outTable = workspace.OpenTable(entity.EntityFeatureClass);                  if (table == null)                 throw new ArgumentException("A tabela de origem não existe.");                  if (outTable == null)                 throw new ArgumentException("A tabela de destino não existe.");                  ICursor selectCursor = table.Search(null, false);             IRow selectRow = null;                  ICursor insertCursor = outTable.Insert(true);             IRowBuffer buffer = outTable.CreateRowBuffer();                  int recordCount = table.RowCount(null);             int count = 0;             int step = recordCount;             if (recordCount > 10000)                 step = Convert.ToInt32(Math.Round(recordCount * .1, 0));                  int remainder = 0;                  while ((selectRow = selectCursor.NextRow()) != null)             {                 buffer = selectRow;                 insertCursor.InsertRow(buffer);                 count++;                 Math.DivRem(count, step, out remainder);                 if (remainder == 0)                     insertCursor.Flush();             }             insertCursor.Flush();                  CustomComReleaser.ReleaseComObject(insertCursor);             CustomComReleaser.ReleaseComObject(selectCursor);         }              public void ExportSelect(IGisEntity entity, int sectorNo)         {             ITable table = _Context.GetEntityTable(entity.EntityName);             IFeatureWorkspace outWork = _OutWorkspace as IFeatureWorkspace;             ITable outTable = outWork.OpenTable(entity.EntityFeatureClass);                  if (table == null)                 throw new ArgumentException("A tabela de origem não existe.");                  if (outTable == null)                 throw new ArgumentException("A tabela de destino não existe.");                  IQueryFilter filter = new QueryFilterClass();             filter.WhereClause = String.Format(@"{0}_GR_IDENTIFICADOR LIKE '1{1}%'", entity.EntityCode.Substring(1, 3), sectorNo.ToString());                  ICursor selectCursor = table.Search(filter, false);             IRow selectRow = null;                  IRowBuffer buffer = outTable.CreateRowBuffer();             ICursor insertCursor = outTable.Insert(true);                  int recordCount = table.RowCount(filter);             int count = 0;             int step = recordCount;             if (recordCount > 10000)                 step = Convert.ToInt32(Math.Round(recordCount * .1,0));                  int remainder = 0;                  while ((selectRow = selectCursor.NextRow()) != null)             {                 buffer = selectRow;                 insertCursor.InsertRow(buffer);                 count++;                 Math.DivRem(count, step, out remainder);                 if (remainder == 0)                     insertCursor.Flush();             }             insertCursor.Flush();                  CustomComReleaser.ReleaseComObject(insertCursor);             CustomComReleaser.ReleaseComObject(selectCursor);         }      When one operating on ArcSDE systems these methods tend to hang and gives me I/O errors. I look into IGeodatabaseBridge interface and IGeoDbDataTransfer, but none of them gives me the option to export a cursor (product of a select, for example). **This is now a major issue: everytime I'm getting a RPC_ERROR_SERVERFAULT with this snippet.**                       while ((row = select.NextRow()) != null)             {                 buffer = row as IRowBuffer;                 insert.InsertRow(buffer);                 count++;                 if (count == 4999)                 {                     insert.Flush(); // this is where the error explodes at my face!                     count = 0;                 }             }      Any choices?


title: Get features in order they are selected using QGIS

text: Is there a way of establishing the order that a user selected features using QGIS and Python? By example: we have a network of line features representing a road network, and the user defines a route from point A to point B. Currently, using layer.selectedFeatures() returns the correct features, but in the order they were drawn rather than the order they were selected. I've tried a few things to do this manually:   * QgsMapToolEmitPoint(canvas) obviously doesn't work - I want the feature not the point clicked. But I guess in worst case I could hack a spatial query to do a look up of the closest feature. Not ideal.   * QgsMapToolSelect(canvas) doesn't seem to be available as part of the Python accessible GUI or CORE, although I can find mention of it existing elsewhere. Can anyone suggest a best case solution for this? All suggestions/tutorials/related questions appreciated. Thanks in advance for your time, Dan


title: With just a web browser we want our field staff to collaboratively make complex maps with secured private data--looking for options?

text: **Simple Question** With just a web browser we want our field staff to collaboratively make complex maps with secured private data. What are our options? Note: Please provide specific product examples to help me and others start the research process and ensure we dont overlook a great solution. I will give a concrete example in response to the comments ( @brad and @iant ) that my earlier version of the question above is too general. **Concrete Real-World Scenario Example** For example, consider there are 100 local people spread around a province in Afghanistan and they only have a web browser. The objective is to have them collaboratively making a map of the safety and danger zones in their province. **Requirements** They do NOT require advanced spatial analysis such a routing etc. Nevertheless, the final map-data product needs to be structured (geoJSON, shapefile, etc) so some one can later use the data for more complex spatial analysis. The key functional requirement is that these people need to upload to the collaboratively edited base map the following new features with associated labels:   * point locations of violent incidents with a labels   * polygon zones of dangerous and safe zones with labels   * line roads with labels  **More elaborate question** What are the specific products I should look into where the data will be secure, away from the public eyes --- _ideally_ something that connects with Google's Datastore, which is approved for SBU data (sensitive but unclassified data) --- _most likely_ some product where spatial data is stored and served internally ?


title: What is the difference between a Custom Layer & a plugin data source?

text: You can add custom data to ArcGIS desktop in two different ways:   * **Plug-in Data Sources**   * **Custom Layer** What is the difference between these two? I'm looking at the possibility of creating an editable layer from some custom format. I already know that it is not possible to create an editable Plug-in Data Source. Is an Editable custom Layer possible? Or is there some other way of achieving what I want?


title: How to draw ploygon in leaflet that use event?

text: Is there only one way to draw a polygon through the Leaflet.draw plugin? If you want draw anything, you must add the control ("map.addControl(drawControl);"). I want to draw a polygon on leaflet map by Event that I can trigger through any DOM event. Does anybody have a good idea how I can do this? Thank you for your time.


title: Simplify shapefile parcel layer

text: I have a very large shapefile with over a million parcels within it. The problem I am having is that the polygons representing each parcel are very jagged. What I am looking for is a simple and fast way of essentially "simplifying" these polygons and thus getting cleaner lines out of things. Naturally, in doing so I need to ensure all parcels are still linked together so there are no spaces between the parcel boundaries. What's the best way of achieving these objectives?


title: How to highlight and zoom to features after applying a definition query programmatically in Arcgis silverlight api

text: A definition query is applied to a bunch of layers after a user queries for a feature( with an ID). The def query is applied but I'm looking for a way to highlight those features that are a result of the query and zoom to them. Any ideas?


title: WGS72 to NAD83- how to find correct transformation

text: I am trying to get this raster file to convert from WGS72 to NAD83, but ArcMap10 doesn't appear to have the correct transformation listed. Where can I find this transformation, and/or how do I import ones that are not already listed as available?


title: How can I get my GPS track data into a graphics program?

text: I'm using a Garmin Edge 705 bike GPS/computer unit, which produces .gpx files. I can use the standard software it comes with to plot on Google Maps etc, but I'd like to be able to output just a vector polyline for use in Adobe Illustrator. Any ideas? I'm guessing .svg is going to be the best format, and I could use Inkscape or similar if that makes it any easier. * * * Sorry, other projects are taking up my time but I will try the suggestions out as soon as I can.


title: Incremental spatial clustering algorithm

text: I'm in search of a incremental spatial clustering algorithm. Here is my use case:   * users create entries with an initial position   * users may change the positions of existing entries I now want to implement a decoupled service that provides clustering information of these data. The service will be notified each time a new entry has been added or an existing entry has been moved. What is a good clustering algorithm therefore? Ideally, it should scale well up to high amounts of data and if there's a trade-off between cluster quality and runtime complexity, I'm okay with degrading results, and _eventually consistent_ results (stale results are okay for some time). To summarize my requirements:   * spatial clustering based on positions   * incremental modifications on changes   * add new positions    * change existing positions   * good runtime performance  Thanks in advance!


title: Need symbology that allows for symbols based on multiple fields, sometimes simultaneously (for certain features)

text: I've got a shapefile containing a large number of O&G parcel features. Each feature parcel has an associated status (things like "Held by Production," "Shut-in," "Leased, Non-producing") and an associated lease expiration date. I need to create a map that depicts this information by color-coding each parcel by its expiration date AND its status. For example, I need to color all parcels with leases expiring in 2014 red, 2015 orange, 2016 yellow, etc. I also need to identify those parcels with Status = "Paid Up" with diagonal lines overlaying the color fill designating expiration year. I hope my question is clear, and thanks in advance for any feedback you can offer! EDIT Throughout the process of working with this problem, I've come up with another tripping point: How does one assign symbology based only on the YEAR of a date field? For instance, my shapefile has a DATE field with the following format: 01/01/2014. I would like to symbolize each feature based on the YEAR only. The symbology tool allows me to assign symbology for each individual value (1/1/2014, 1/2/2014, etc.), but not based on only the YEAR portion of the DATE field. In fact, when I populate the symbology tool with all unique values of the DATE field and sort, the list sorts alphabetically rather than chronologically. This leads to weirdness like 1/1/2015 appearing before 1/2/2014.


title: To get FeatureInfo in OpenLayers

text: I am new to OpenLayers and I want to add OpenLayers.Control.WMSGetFeatureInfo to get feature   info to website.But i can't figure out right method. I also see the source code of GeoServer OpenLayers page and still I can't find the syntax of OpenLayers.Control.WMSGetFeatureInfo code? I use geoserver,apache-tomcat-6.0.35 and open layer 2.10.Should i use proxy to call OpenLayers.Control.WMSGetFeatureInfo? This is my GeoServer WMS URL:               http://localhost:8080/geoserver/sac/wms?     service=WMS&     version=1.1.0&     request=GetMap&     layers=sac:contour&     styles=&     bbox=187754.5,204980.703125,192031.40625,210022.890625&     width=434&height=512&s     rs=EPSG:5235&     format=application/openlayers      Pls give me sample code which calling OpenLayers.Control.WMSGetFeatureInfo? Thank you.


title: add cookies in WMS request

text: I need to add a cookie in each WMS request. i tried to do something like this but i don't think this is the right way.               getURL : function(bounds) {     var requestString = OpenLayers.Layer.WMS.prototype.getURL.apply(this, arguments);     window.document.cookie="referrer="+window.location.origin;     return this._encode(requestString) + "&KEY=" + this.navkey;      is it possible to do this? thanks you


title: Flow mapper plugin failed

text: I want to create a desire line with the Flow Mapper pluginwith qgis ver.2.1 . After I did it has failed. The data that I use as an attachment below. There may be error procedure I have done?. Expected clue how to do it right. Thank you for the response. ![Data](http://i.stack.imgur.com/W47rM.png)


title: GeoExt: Layer leaf nodes replicated across all tree nodes

text: This is going on with GeoExt 1.1 and OpenLayers 2.12. I need to build a layer tree dynamically, with layers and nodes that are only known at run time. The snippets I show below correctly create a layer tree with several nodes and adds the layers node leafs. The problem is that though each layer is assigned to a single tree node, it is then shown under each and every tree node created. During init the tree and its root are set up:                 treeRoot = new Ext.tree.TreeNode       ({           text        : 'Project',           id          : '0',           leaf        :false,           expandable  : true,           expanded    : true       });            layerTree = new Ext.tree.TreePanel       ({         title: 'Map Layers',         region: "west",         width: 200,         collapsible: true,         autoScroll: true,         enableDD: true,         plugins:         [{           ptype: "gx_treenodecomponent"         }],         loader:         {           applyLoader: false,           uiProviders:           {             "custom_ui": LayerNodeUI           }         },         root: treeRoot,         rootVisible: false,         lines: false       });      Later on new tree nodes can be added according to a tag, and kept in an array:               if(treeNodes[tag] == null)     {         treeNodes[tag] = new GeoExt.tree.OverlayLayerContainer         ({                 text: tag,                 leaf: false,                 expanded: true         });         treeRoot.appendChild(treeNodes[tag]);     }      Layers are added to the map and to a specific tree node:               var layer = new OpenLayers.Layer.WMS(title, serviceURL, params, options);     leftMap.addLayer(layer);     layer.events.register("visibilitychanged", this, toggleLayer);          var newNode = new GeoExt.tree.LayerNode(     {         text: title,         layer: layerName,         leaf: true,         checked: false,         children: [],         nodeType: "gx_layer"     });     treeNodes[tag].appendChild(newNode);      But then each layer appears several times, once under each tree node. Any hints on why it is happening this way are welcome. Thank you.


title: How to project a geoJSON file in OpenLayers?

text: I have a json file well formed. I want to read a json file as a base layer from OpenLayers.                   <script type="text/javascript">             var styles = new OpenLayers.StyleMap({                     "default": {                         strokeWidth: 2                     },                     "select": {                         strokeColor: "blue",                         strokeWidth: 4                     }                 });                      // add rules from the above lookup table                 styles.addUniqueValueRules("default", "RP_TYPE", {                     10: {strokeColor: "#000000", strokeWidth: 2},                     12: {strokeColor: "#222222", strokeWidth: 2},                     14: {strokeColor: "#444444", strokeWidth: 2},                     16: {strokeColor: "#666666", strokeWidth: 2},                     18: {strokeColor: "#888888", strokeWidth: 2},                     19: {strokeColor: "#666666", strokeWidth: 1}                 });                           // add rules from the above lookup table                 styles.addUniqueValueRules("default", "RP_TYPE", {                     10: {strokeColor: "#000000", strokeWidth: 2},                     12: {strokeColor: "#222222", strokeWidth: 2},                     14: {strokeColor: "#444444", strokeWidth: 2},                     16: {strokeColor: "#666666", strokeWidth: 2},                     18: {strokeColor: "#888888", strokeWidth: 2},                     19: {strokeColor: "#666666", strokeWidth: 1}                 });                      var vectors = new OpenLayers.Layer.Vector("Lines", {                     strategies: [new OpenLayers.Strategy.Fixed()],                                     protocol: new OpenLayers.Protocol.HTTP({                         url: "data/locality_region.json",                         format: new OpenLayers.Format.GeoJSON()                     }),                     styleMap: styles                 });                      map.addLayer(vectors);                 map.addControl(new OpenLayers.Control.LayerSwitcher());                 map.zoomToMaxExtent();                  }        I have another json files which are downloaded from api and those files works fine. How about my own json file that got converted from db. Do I need to work more on it. If so please guide me how to project a json on OpenLayers that created from a PostgreSQL db. How can I fix the bounds for my own map or layer; Thank you


title: Error on Geoserver Web Admin Page

text: When i click on Geoserver Web Admin Page gives Message as 'Bad Request (Invalid Hostname)' in web page. I have already data of on going project in geoserver. I have started Geoserver. My geoserver version 2.1.3. in window Xp operating System. Can any one please guide me to get out of the error and get back Geoserver web Admin Page. Please see the Page result in images below.![enter image description here](http://i.stack.imgur.com/6t8hw.png)


title: Get the vertices on a LineString either side of a Point

text: I have a shapely LineString and defined a shapely Point which lies along the LineString. How can I find the vertices of the LineString which lie either side of the point? (split the line in two). This is very similar to this blog post by the author: http://sgillies.net/blog/1040/shapely-recipes/ Except, instead of cutting by distance along the line, I want to cut given a point that lies along the line.


title: how display graduated symbol in the legend in print composer?

text: this is the problem: why the graduated symbol is not reported at all in the legend of the print composer? it is a bug or i am doing something wrong? here below the image that shows the problem ![enter image description here](http://i.stack.imgur.com/QIIGR.jpg)


title: Mapinfo and qgis file's datum problem

text: please refer to following link for the explanation of the problem I had. do you have any idea on this? http://osgeo-org.1560.x6.nabble.com/Mapinfo-and-qgis-file-s-datum-problem- td5091891.html#a5100434


title: QGIS 2.0 Plugin: Select and Identify a Single Feature

text: given is a layer with some nodes on it and every node has some attributes(name, id, etc...). For showing up the informations of a single node I can use the "Identify Feature Tool". Choose the tool, click the node and see the informations. I am trying to develop a Plugin with Python, which is doing almost the same: Start the Plugin, click on a node and get the informations about this node(output in the python console is enough for the first steps). Okay, so this is the status: I was following this tutorial: http://www.qgisworkshop.org/html/workshop/plugins_tutorial.html which works not for QGIS 2.0. The 'MouseClickEvent' is working, I get the X,Y-coordinates. But I am not able to select and to get the informations of a node. That is my selectFeature function, which is called by `result = QObject.connect(self.clickTool, SIGNAL("canvasClicked(const QgsPoint &, Qt::MouseButton)"), self.selectFeature)`:               def selectFeature(self, point, button):         #pntGeom = QgsGeometry.fromPoint(point)         #pntBuff = pntGeom.buffer( (self.canvas.mapUnitsPerPixel() * 2),0)          #rect = pntBuff.boundingBox()              layer = self.canvas.currentLayer()         selected_features = layer.selectedFeatures()         for i in selected_features:             attrs = i.attributeMap()             for (k,attr) in attrs.iteritems():                 print "%d: %s" % (k, attr.toString())      This let me see nothing. I was reading and trying out a lot of things and at the moment in feel a bit lost. I lost the overview. Do I need these geometry functions to capture the node by clicking the mouse button? What functions do I need to get the attributes of a single selected feature? Thank you :)


title: Cleaning up Line Features

text: I'm in the process of creating some routes for which I need left and right sections. I have an original streets feature class, which I used "copy parallels" on, and came up with something like the picture below. I'm wondering if there is any tool that would clean up these roads, erase some dangles and intersecting lines. I've tried using the "trim" tool, but cannot come up with a standard length to trim. Any ideas? ![Routes](http://i.stack.imgur.com/7a05c.png) Using ArcGIS 10.2.1 Advanced


title: Setting ArcGIS Server 10.2 Data Store Folder as Local vs. Network (UNC) path

text: Is there any difference? How setting ArcGIS Server Data Store Folder as local vs. network path will affect the publishing process? I saw warnings in ArcMap 10.1 when I Analyze the MXD before I publish, complaining that network path is used. (Assuming this might slow down the REST service if network access to MXD data is slow) Example: My file server has designated shared folder for ArcGIS data called "GisDataStore". It is mapped as local drive on my ArcGIS server, so I can access it both ways: \MyFileServer\GisDataStore and Z:\GisDataStore I wonder which I should use to Register the "GisDataStore" folder as ArcGIS 10.2 Server datastore. ![Local Path](http://i.stack.imgur.com/yWrCy.png) ![Network Path](http://i.stack.imgur.com/XWG8F.png) Here's the warning I am getting: SEVERITY STATUS CODE DESCRIPTION Medium Unresolved 10027 Layer's data source is referenced via a UNC path ![UNC Path Warning](http://i.stack.imgur.com/D8HaB.png)


title: How to generate JPEG tiles using gdal2tiles.py

text: I am using gdal2tiles.py to generate tiles and packing it into MBTiles using mb-util. gdal2tiles by default generates PNG tiles. When I convert PNG tiles to JPEG large black area around tiles are left (due to lack of transparency in JPEG): http://www.imagemagick.org/discourse-server/viewtopic.php?f=1&t=24048 Is there a way to directly generate JPEG tiles in gdal2tiles.py ? I am using GDAL 1.9.2 on Ubuntu 12.04.


title: Using the python shape library: pyshp - how to convert polygons in .csv file to .shp

text: I'm converting csv data to shapefiles using the pyshp python library. What I want to do is very similar to Using the python shape library: pyshp - how to convert .csv file to .shp, except that I don't know how to handle polygons. My data is in a csv file. Each row corresponds to a named rectangle, i.e. a string and the latitude or longitude coordinates of each side of the rectangle: left side (xl), right side (xr), top (yt) and bottom (yb). So my data looks like this:               name,xl,xr,yt,yb     some Name,-25.3125,22.5,47.517193,31.353634     another Name,-103.359375,-0.703125,80.87282,74.40216     ...      And my python code is pretty simple. It is only slightly modified from the points example. But when I try to import this data into google maps, it has some errors parsing it. So I think I'm doing something wrong, but I'm not sure what?               #Set up blank lists for data     name, polyPart = [],[]          #read data from csv file and store in lists     with open(in_file, 'rb') as csvfile:       r = csv.reader(csvfile, delimiter=',')       for i,row in enumerate(r):         if i > 0: #skip header                # parse data into a point array representing the bounding box           xl = float(row[1])           xr = float(row[2])           yt = float(row[3])           yb = float(row[4])           tl = [xl, yt]           tr = [xr, yt]           br = [xr, yb]             bl = [xl, yb]           parr = [tl, tr, br, bl, tl]                # array of one "part", the part is an array of points           polyPart.append([parr])           name.append(row[0])          #Set up shapefile writer and create empty fields     maxStringLength = 50     w = shp.Writer(shp.POLYGON)     w.field('name','C',maxStringLength)          #loop through the data and write the shapefile     for j, name in enumerate(name):       w.poly(parts=polyPart[j])       w.record(name)          #Save shapefile     w.save(out_file)


title: Attempting to use ST_DWithin to locate objects within a given distance of meters using SRID 4326

text: I'm a Rails developer with minimal GIS experience, working with PostGIS and the RGeo gem. I'd like to find all of the cows within 500 meters of my current location. So, I'll create a model for my cows:               class Cow < ActiveRecord::Base     set_rgeo_factory_for_column(:latlon, RGeo::Geographic.spherical_factory(:srid => 4326))        I'll then seed my database with enough cows to keep me busy:               5000.times{ Cow.new(latlon:"POINT(#{40.735767 + rand() * 0.1  } #{-77.919674 + rand() * 0.1  })" ).save }      To find the cows within 500 meters of my location, I'll use the ST_DWithin PostGIS function that all the tutorials suggest:               cows_found = Cow.find_by_sql("SELECT * FROM cows WHERE ST_DWithin(latlon, '#{search_latlon}', 500)")      The unhappy result is the following: ![enter image description here](http://i.stack.imgur.com/c3J7O.png) Moving to the East or the West, the 500 meter radius threshold appears to be observed correctly, but moving to the North or the South, I am finding cows as far as two kilometers away. Now, obviously I'm extremely naive in my approach, but I'm not completely clueless; I understand I'm populating my cows in a spheroidal projection and handing the coordinates returned from ST_DWithin to Google's Mercator. I expect distortion. However, this is larger than I anticipated and suggests to me that I've fundamentally misunderstood something. I had thought that in my use of ST_DWithin I was passing as arguments two geographies, and that distance would therefore be specified in meters. I was able to set the use_spheroid boolean argument to false without PostGIS complaining, which I do not think I should be able to do if I were in fact comparing geometries. In fact, if I use psql to look at the latlon column, it is of type geography:                latlon     | geography(Point,4326)       |      However, perusing related questions, I see several answers that state that using ST_DWithin with SRID 4326 results in the distance threshold being specified in radians. (Edit: I now understand this would only be the case if one was using 4326 geometries, not 4326 geographies.) So, given the goal of locating cows within a given number of meters, and tolerating an error of tens of meters (nobody will die if we find a few extra cows) what would be the best way to proceed? Throw out the 4326 spatial reference system, and use a different SRID? Use 4326, but calculate a conversion from meters to radians, for use with ST_DWithin? (This conversion would obviously be subject to the specific latitude and longitude from which I am searching.) Any suggestions would be appreciated. Edit: Both arguments to the ST_DWithin call are 4326 geographies. The first argument, latlon, is a 4326 geography column, and the second argument is drawn from another column of the same type. A commenter's (very perceptive) suggestion that I ensure the second argument was in fact a geography by using ST_GeogFromText led to precisely the same end result. So, I have high degree of confidence that I am using the correct version of ST_DWithin with the correct parameters. As a sanity test, I repeated exactly the same procedure (generating a dense set of latlons and performing ST_DWithin upon it) very near the equator. This is the result: ![enter image description here](http://i.stack.imgur.com/HutR0.png) As you can see, the points returned by ST_DWithin are much less distorted. More importantly, they are all within 500 meters of the central point. The results that I am able to obtain at the equator are precisely in line with my expectations. However, at roughly the 41st parallel North, I am getting the distortion observed in the first example. Now, I expect distortion from handing my 4326 coords to Google Maps (who internally transform them to their own web Mercator) and I expect some distortion from using 4326 (PostGIS In Action warns that local distance measurements will be 'bad'). I expect these distortions to become more severe as I move away from the equator. The first map example in this question is obviously subject to both of those distortions. So my question is: Is the level of distortion I am seeing in the first map example actually a correct result? (Meaning, the result one would expect given the limitations of the 4326 spheroid model and distortion from the projection to web Mercator?) If this is a correct result, I would be eager to hear suggestions regarding a PostGIS strategy to retrieve points within a given radius that was less subject to these distortions, or a strategy to ameliorate this distortion.


title: Function to determine font color of data driven page attribute in ArcMap 10.2

text: ArcGIS 10.2 I am looking for an IF/ELSE statement that I can place in the text properties of my Data Driven Page Attribute item to change the font color based on the value in that field. Here is the function I have been trying. Its a variation on a function for basic map label. I can't get this function to work w/ dynamic text. `Function FindLabel<dyn type="page" property="attribute" field="RP_Slope" domainlookup="true"/> > 2 then FindLabel = <CLR red='255'><dyn type="page" property="attribute" field="RP_Slope" domainlookup="true"/></CLR> else FindLabel = <dyn type="page" property="attribute" field="RP_Slope" domainlookup="true"/> end if End Function`


title: osm2po: import multiple .osm files for disparate regions?

text: Does osm2po support importing **multiple** OpenStreetMap files into the same Postgres/PostGIS database (i.e., to be used for routing)? If yes, can the OSM files be for geographic regions that are **not** next to each other? For example, could I import something like NewYork.osm and SanFrancisco.osm into the same routing database? Thanks!


title: Bing map v7 adding infobox to georss parser module

text: Using this plugin (source) to load pins from an xml file, pins get added correctly, but now trying to add infoboxes to each respective pin. In the function that receives the array of pins I'm trying to loop through each item and add infoboxes accordingly. Any ideas as to what I'm missing here?                     function GeoRSSImportedCallback(items, bounds) {               map.entities.push(items);                    for (var i = 0; i<items.length; i++) {                 items[i].Title = items[i].title;                 items[i].Description = items[i].description;                 Microsoft.Maps.Events.addHandler(itmes[i], 'click', displayInfoBox);               }                    map.entities.push(infoboxLayer);               map.setView({ bounds: bounds });           }              function displayInfoBox(e) {             pinInfoBox.setOptions({title: e.target.Title, description: e.target.Description, visible:true, offset: new Microsoft.Maps.Point(0,25)});             pinInfoBox.setLocation(e.target.getLocation());     }      I do initialize the following in the `GetMap()` function                       pinInfoBox = new Microsoft.Maps.Infobox(new Microsoft.Maps.Location(0, 0), { visible: false });             infoboxLayer.push(pinInfoBox);


title: scalebar doesn't print properly in IE or FF

text: I have a ruler-style scalebar in my webapp, which uses the ArcGIS JS API version 3.1. When I send the map to print using window.print(), the scalebar renders as expected in Chrome, but in both IE and Firefox, the scalebar renders as a hollow bar with no markings. From what I can tell, the scalebar=line option is not available for version 3.1 of the API, and there are other dependencies in the code that prevent me from changing the API.


title: Find index location of a lat/lon point on a raster grid in ArcPy

text: This question has been confusing me for a while, and I'd love y'all's help. I'm trying to find the index location on a raster (in terms of column and row#) of a specified lat/lon point. My first thought was to find the cell value at the lat/lon point, convert the raster to a numpy array, and find the index of the point (but that's clunky, and requires each cell value to be unique...which isn't true in the raster)               #get values, assuming the extent object is already saved     raster = "raster.shp"     x_coord = extent.XMin     y_coord = extent.YMin     point = arcpy.Point(x_coord, y_coord)     cell_val = arcpy.GetCellValue_management(raster,point,"")     raster_array = numpy.array(raster)     final_index = raster_array.index(cell_val)     print final_index      But that's ungainly, and also requires the cell values to be unique in the raster (which again isn't the case) Is there any other way to find the cell index of a pixel within a raster that corresponds to a point object? Thanks!


title: layer.visible command not working in arcpy.mapping module

text: I am using the arcpy.mapping module for the first time, and the lyr.visible command is not working as it should (or I am not using it as I should). I have a list that contains the the layer names of .img files that are loaded into a single data frame - the list is called "layernames". I am using a while loop to go though the layernames list, turn on one of the layers at a time, export a .png file with that layer visible, then turn that layer off and move on to the next layer in the list. At the beginning of each pass through the loop I specify which raster layer should be turned on by setting layernames.visible = True, and at the end of the pass through the loop I turn off that layer by setting layernames.visible = false. The first time through the loop the script works and does turn on the correct layer in the list "layernames". But for all subsequent passes through the loop, only the first layer remains visible, and none of the others are turned on or off. In other words, every .png image shows the same raster (the first one that is turned on), instead of a different raster in each .png image. I have pasted my script below. Any help would be greatly appreciated!               mxd = arcpy.mapping.MapDocument("Z:\\sharks\\Models\\Summer2013\\daily figures.mxd")     layernames = arcpy.mapping.ListLayers(mxd)     counter = 0     yr = 2005     while counter < 365:         layernames[counter].visible = True         x = str(layernames[counter])         y = x[9:12]         j = int(y)         date = datetime.datetime(yr,1, 1)+ timedelta(days=(j-1))         strDate = str(date)         if (j + 1)<= 31:             mxd.title = "January (" + strDate[0:10] + ")"         elif (j + 1) >= 32 and (j+1) <= 59:             mxd.title = "February (" + strDate[0:10] + ")"         elif (j + 1) >= 60 and (j+1) <= 90:             mxd.title = "March (" + strDate[0:10] + ")"         elif (j + 1) >= 91 and (j+1) <= 120:             mxd.title = "April (" + strDate[0:10] + ")"         elif (j + 1) >= 121 and (j+1) <= 151:             mxd.title = "May (" + strDate[0:10] + ")"         elif (j + 1) >= 152 and (j+1) <= 181:             mxd.title = "June (" + strDate[0:10] + ")"         elif (j + 1) >= 182 and (j+1) <= 212:             mxd.title = "July (" + strDate[0:10] + ")"         elif (j + 1) >= 213 and (j+1) <= 243:             mxd.title = "August (" + strDate[0:10] + ")"         elif (j + 1) >= 244 and (j+1) <= 273:             mxd.title = "September (" + strDate[0:10] + ")"         elif (j + 1) >= 274 and (j+1) <= 304:             mxd.title = "October (" + strDate[0:10] + ")"           elif (j + 1) >= 305 and (j+1) <= 334:             mxd.title = "November (" + strDate[0:10] + ")"         elif (j + 1) >= 335 and (j+1) <= 365:             mxd.title = "December (" + strDate[0:10] + ")"           mxd.save()         out_png = "Z:/sharks/Models/Summer2013/RF/ALL/SB_RF/AVpresRFP/figures/" + strDate[0:10] + ".png"         arcpy.mapping.ExportToPNG(mxd, out_png)         layernames[counter].visible = False         mxd.save()         counter +=1


title: Geoserver WMS v Geowebcache offset?

text: Having successfully used the new OpenLayers.Layer.GeoWebCache to add a local tile cache created using Geowebcache, I thought I would compare against the same data as a wms from my local GeoServer. Looking at the geowebcache.js file there is an option to specify the cache layer as a base layer, so I sent it to false and added a layer switcher to the map Both layers add to the map however there is a significant offset. ![enter image description here](http://i.stack.imgur.com/E8lIU.png) The cache layer is the one that is set to the right. In the geowebcache.js there is no mention of offset needing to be set. Any ideas on what would cause this and what I need to change? thanks ## UPDATE As requested my openlayers code               <script src="http://openlayers.org/api/OpenLayers.js"></script>     <script src="geowebcache.js" type="text/javascript"></script>     <script type="text/javascript">          var map;             function init() {                      var world = new OpenLayers.Bounds(                     0, 0, 700000, 1344000                 );                      map = new OpenLayers.Map({                     div: "map",                     projection: new OpenLayers.Projection("EPSG:27700"),                     numZoomLevels: 11,                     resolutions: [896,448,224,112,56,28,14,7,4,2,1],                     maxExtent: world,                     center: new OpenLayers.LonLat(400000, 500000),                     zoom: 11                 });                      map.addControl(new OpenLayers.Control.LayerSwitcher());                      var cache = new OpenLayers.Layer.GeoWebCache({                       url: "https://url/",                       isBaseLayer: false,                       name: "Cache Layer",                       gridSetId: "zoommap",                   });                      // Add WMS layer to our map                 map.addLayer(cache);                      var bwms = new OpenLayers.Layer.WMS( "Base Map",                   "http://localhost:8080/geoserver/wms?",                   { "transparent":"false",                     isBaseLayer: true,                     "layers":"osgb:aa_base",                     "format":"image/png" },                   { "reproject":"true" }                   );                      // Add WMS layer to our map                 map.addLayer(bwms);      Here is my GWC configuration file                 `          <?xml version="1.0" encoding="utf-8"?>         <gwcConfiguration xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://geowebcache.org/schema/1.5.0 http://geowebcache.org/schema/1.5.0/geowebcache.xsd" xmlns="http://geowebcache.org/schema/1.5.0">       <version>1.5.0</version>       <backendTimeout>120</backendTimeout>       <gridSets>         <gridSet>           <name>zoommap</name>           <srs>             <number>27700</number>           </srs>           <extent>             <coords>               <double>0.0</double>               <double>0.0</double>               <double>700000.0</double>               <double>1344000.0</double>             </coords>           </extent>           <alignTopLeft>false</alignTopLeft>           <resolutions>             <double>896.0</double>             <double>448.0</double>             <double>224.0</double>             <double>112.0</double>             <double>56.0</double>             <double>28.0</double>             <double>14.0</double>             <double>7.0</double>             <double>4.0</double>             <double>2.0</double>             <double>1.0</double>           </resolutions>           <metersPerUnit>1.0</metersPerUnit>           <pixelSize>2.8E-4</pixelSize>           <scaleNames>             <string>zoommap:0</string>             <string>zoommap:1</string>             <string>zoommap:2</string>             <string>zoommap:3</string>             <string>zoommap:4</string>             <string>zoommap:5</string>             <string>zoommap:6</string>             <string>zoommap:7</string>             <string>zoommap:8</string>             <string>zoommap:9</string>             <string>zoommap:10</string>           </scaleNames>           <tileHeight>250</tileHeight>           <tileWidth>250</tileWidth>           <yCoordinateFirst>false</yCoordinateFirst>         </gridSet>       </gridSets>       <layers/>     </gwcConfiguration>      The data is stored in a PostGIS database and the projection is definitely EPSG:27700. As you can see I have defined the projection in my OpenLayers code and included the same resolutions as my geowebcache configuration resolutions.


title: Troubleshooting why ArcGIS Model works on one PC but not another?

text: Do you have specific expertise you could share of examples where an ArcGIS Modelbuilder has worked on one PC but not on another? Are there any issues between different Windows versions? There is a related question here Model not editable on a different PC In my example I have a model that runs to completion on a Windows 8 machine. I run the models from an external hard drive. When the same model is run from the same external drive but on a different PC - Windows 7 then the model crashes. The model can be opened, edited and verified on both Pc's.It just crashes on one but not the other. I have explored possible options: ArcGIS Version and service packs - both are 10.1 SP1.   RAM - the PC on which the models work has LESS RAM than the PC on which the models crash OS type - both machines are 64 BIT.   Model build - i have re-built the model by copying and pasting into new / clean models (I did this on the Windows 8 PC though) I am not sure that the cause is the different Windows version, it is just the main difference I can see between the two PC's Any expertise, answers or references to sources of information would be much appreciated. I am using ArcGIS 10.1 SP1 for Desktop Advanced.


title: How to export all the layers from Table of Content to shp file?

text: I have got lots of different layers in the Table of Content panel in ArcMap. I want to export them all in shp format in one go. How can I do that? So far I have been doing it one by one by right clicking on the layer-->data-->export as shp. and it is taking longer time. I want to do them all together.


title: Polygonize.py: AttributeError: 'NoneType' object has no attribute 'CreateDataSource'

text: **Given a topographic raster tif** of an european country where elevations are simplyfied to **either 0, either 200**. **How to create a polygon ?** * * * **Comment** : With the following polygonize.py command:               gdal_polygonize.py level200.tif -f "ESRI_Shapefile" level200.shp      I get to following error:               Creating output level200.shp of format ESRI_Shapefile.     Traceback (most recent call last):       File "/usr/bin/gdal_polygonize.py", line 164, in <module>         dst_ds = drv.CreateDataSource( dst_filename )     AttributeError: 'NoneType' object has no attribute 'CreateDataSource'


title: Pure JavaScript Library for Geometry Operations

text: Does anyone know of a pure JavaScript library for doing geometry operations such as buffer, union, intersects, etc? I want to play with doing these operations in the browser. It looks like GeoScriptJS has some dependencies.


title: Are there any Indoor GIS platforms?

text: What sort of infrastructure is needed to allow me to write a simple mobile app for indoor navigation? Professor Mike Goodchild said: > the average American only spends 13 percent of the time outdoors. GIS-based > services help us find restaurants and hotels, but they offer almost no > support for navigating in the complex indoor spaces of shopping centers, > hospitals, mines, or airports. On my smartphone, I'd like to show a floor plan oriented to the real world, centered on my location.


title: Pure Python Library for Geometry Operations

text: Does anyone know of a pure Python (no C or Java dependencies) library for doing geometry operations such as buffer, union, intersects, etc? I'm looking for something that Google AppEngine compatible. Thanks!


title: Topology Edit Tool in ArcGIS Engine

text: Is a Topology Edit Tool available "out of the box" for ArcGIS Engine 10? Or do I have to write one myself?


title: Renaming legend items in ArcGIS for Desktop?

text: Can someone please tell if it is possible to rename legend items in the map layout in ArcGIS for Desktop? If yes, how can I do it?


title: Customize geocoding locator styles in an ArcGIS Server environment

text: Using ArcGIS Desktop, it is relatively straightforward to customize locator styles by editing the corresponding CLS files. We are struggling to do the same in a server environment. We have created a locator in ArcSDE (9.3.1) for Oracle 11g and a corresponding _Geocode Service_ in ArcGIS Server (9.3.1). How do we customize locator styles in this environment?


title: ReClassify Vector Data Between Values

text: I would like to be able to take all values between 1 and 10 and have a field define them as 1, and all values between 10 and 20 and set them to 2, how do I do that with vector data? I read how to set all building to 20, and grass to 40, but I'm pretty sure that to do that for a range of values, like 1-10, I would have to manually enter each number, and decimal in that range. Am I wrong? and if not, how would I accomplish this? Thanks


title: Is there a way for multiple users to use an SQLite/SpatiaLite db simultaneously?

text: Is there a way that multiple users can use the SQLite / SpatiaLite database simultaneously? Thank you for your feedback.


title: How can I set the output feature type name dynamically?

text: How can I set the output feature type name so that it consists of a constant part and a variable part, such as "Circle_205", where the number part could be different each time? When using a fanout, this is easy to set, but without a fanout, I do not know how to do this, unfortunately.


title: Calculating Mean Upslope slope, curvature, etc

text: I'm trying to generate rasters whose values summarize the mean value of a given parameter (slope, curvature, etc.) of the upslope terrain. I feel like this should be doable using only a flow direction raster and a raster of the parameter of interest, but I don't really know where to begin. I've seen this instruction http://forums.esri.com/Thread.asp?c=93&f=995&t=303114, but I can't seem to figure out a way to interpret it that produces a result that makes any sense. I'm familiar with ArcMAP, QGIS, SAGA, and a bit of GDAL and GRASS, so I'll gladly accept instructions for any of these platforms.


title: geowebcache tiles not showing up in folder

text: I am using geoserver 2.5.1 and am atempting to use the intergrated geowebcache.I have installed tomcat and am using the webapp process for the geoserver. I managed to add all the tiles and got the jai working. when i had tested out the first time of installing geoserver and the tiling on my local host it worked in seeding the layers and creating cached tiles. The only problem is i don't exactly remember what I did and all the things that i remember what was dont i did it. I then istalled the tomcat and geoserver on our developer server so it is now using an ip address instead of a local host, but for some reason all the things i tried with the localhost does not seem to work on the developer one. I also noticed when I either add:               <context-param>             <param-name>GEOWEBCACHE_CACHE_DIR</parm-name>             <parm-value> C:\temp</param-value>         </context-param>          then it seems to make pink tiles on my map. and when I create a folder after setting up the environmental variables at the system variables in the C:\temp, when I create a folder in my C drive that matches the environmental variables cause i thought u might have to create the folder for it to work but it just seems to make my geoserver work disapear and makes it as if has installed all over again but when I delet the folder then i get my data back. Can any one help with this problem. I have looked through all the questions that vaigley match what im trying to do but to much availe im still having problems and the geoserver and geowebcache documentation is not making too much sence. When I click on the side of the geoserver under serice Capabilities when I click the WMS section it goes to a 404 page on either one of those but when I click on my localhost it work fine and gives me an xml page.


title: vector layer, no more editing is allowed

text: I have a polygon vector layer with already 529 objects but without any attributes so far and now I get the message that shapefiles without attributes are treated as write protected. What can I do to go on working with my shapefile?


How do I turn off labels that have null values? I am using a leader with this label but I want it turned off if the values are null


title: Possible to downgrade Postgresql and PostGIS versions and still be compatible with latest OpenGeo Suite?

text: I'm trying to connect ArcGIS Desktop 10.1 to PostGreSQL/PostGIS tables. I have PostGreSQL 9.2/PostGIS 2.0, which was installed with the latest OpenGeo suite. Unfortunately ArcGIS direct connections to PostGres are only supported by ESRI for version 9.1.3. Would it be possible to downgrade my version of PostGres 9.2 to 9.1.3 while maintaining the other up-to-date components of the OpenGeo Suite (i.e. geoserver, geowebcache, etc.)?


title: GeoServer WCS geotiff style not loading in ArcMAP

text: I have published a GeoTIFF file using GeoServer's (2.4.5) WCS service. I also defined a colormap style on GeoServer. I know the style works, because, when I load the map through WMS the colormap is there and the legend image is also generated. But when I try to load it into ArcMAP or QGis using WCS, I just get a stretched gray scale. Opening the same GeoTIFF file directly from the catalog in arcmap, lets me calculate the statistics and histogram, and then I can define and edit a colormap for that layer, but this option to calculate the statistics for the WCS raster is not available. I also tried to create a classified or unique value colormap in the symbology tab at the properties dialog, but those options are not available for the WCS, only for the file opened directly. How can I make ArcMap take into account the style information I created using GeoServer or, alternatively, how can I make ArcMap let me create and use a classified colormap for this WCS raster. Thanks to anyone who tries to help.


title: Pass Queried Feature Set to geoprocessor as a parameter

text: I am new to geoprocessing services and tried to publish my first. I created a simple model with one parameter, a featureset, and I am intersecting it with another layer which is specified in the model and uploaded to the server. If I pass the required parameter,as below, it never executes. If I pass an empty parameter object, var params = {}, it runs.                   var app = {};         app.qtCounties = new QueryTask("http://myserver/arcgis/rest/services/PA_Counties/MapServer/0");         app.qCounties = new Query();         app.qCounties.returnGeometry = true;         app.qCounties.outFields = ["*"];              paCounties.on("click", getData);function getData(evt) {                 console.log("In GetData The Event graphic is " + evt.Graphic + "Point is ", evt.mapPoint.x, evt.mapPoint.y)                 var county, promises, qGeom, point, pxWidth, padding;                      point = evt.mapPoint;                 pxWidth = map.extent.getWidth() / map.width;                 padding = 2 * pxWidth;                 qGeom = new Extent({                     "xmin": point.x - padding,                     "ymin": point.y - padding,                     "xmax": point.x + padding,                     "ymax": point.y + padding,                     "spatialReference": point.spatialReference                 });                      app.qCounties.geometry = qGeom;                      county = app.qtCounties.execute(app.qCounties);                 console.log("deferreds: ", county);                 promises = new all([county]);                 promises.then(handleQueryResults);                 console.log("created d list");                           function handleQueryResults(results) {                     console.log("queries finished: ", results);                          gp = new Geoprocessor("http://myserver/arcgis/rest/services/TestIntersectModel/GPServer/TestIntersectModel");                     gp.setOutSpatialReference({ wkid: 4326 });                               var params = { "Input_Polygon": results[0].features };                     gp.execute(params, completed);                     function completed(res) {                         console.log("Completed GeoProcssing", res)                     }                               console.log("Params ", params)}      When the query runs I get a valid object in results which is set as the lone parameter. Apparently I am not passing the parameter as a featureset but not sure how else to pass it? Any help would be appreciated. I published my service using this example. Though, I did not publish it as a map service, the last step, and I am trying to access it via javascript api instead of the server web page method. The model works fine when run in desktop and providing the parameter manually.


title: How do you edit a stack of polygons simultaneously?

text: Using ArcMap 10.1. Lets say you have 100 polygons that all are in the exact same shape and in the exact same place. They appear to be one polygon. How would you edit the shape of all 100 at the same time, can you share vertices among the 100 polygons? Can you join all 100 and keep their separate attributes and Object ID's?


title: Sending coordinates via VHF Radios to Computers

text: I need to send and receive coordinates in NMEA format from a VHF radio to a computer. Any suggestions where to look in order to accomplish this task? The transmission will occur only when the user presses the PTT button. Thank you all.


title: Halo in OpenLayers (Line, Vector)

text: How to make a halo around the line of one pixel? For example, line (vector) hase blue color and its frame (one pixel) the black color. Thanks!


title: Vector Building Footprint Resources - Tokyo

text: I typically know where to look for this data in the United States, but I'm at a loss when looking for building footprint shapefiles abroad. I thought I could find what I was looking for with Open Street Map data, but when I check out http://downloads.cloudmade.com/asia/eastern_asia/japan/tokyo#downloads_breadcrumbs I don't get the files I need. Anyone know other resources, either official or osm?


title: How geographically accurate is Google Earth?

text: I'm using a KML file which I have loaded in to Google Earth. I am looking at some islands, and see items such as the following island: ![enter image description here](http://i.stack.imgur.com/nSYhr.jpg) Note the blue region is the polygon describing the area provided, and the image in the background. Sometimes I see similar shapes, but they are off a bit, such as below: ![enter image description here](http://i.stack.imgur.com/iBSpd.jpg) Bottom line is, I'm trying to get as accurate of a map of some of these island nations that I can. Is it fare to use the image that comes from Google Earth to attempt to make a better representation of such islands, or is there enough inaccuracy in the creation of Google Earth islands that it is not worth making any changes?


title: Spatial Views in SDE - What causes them to become unregistered from the geodatabase

text: Spatial views in ArcSDE allow you to join columns from a 'non-spatial' table in the database to a table with a 'geometry column'. They can also allow you to format and choose the attribute columns that you want to have included in your feature class. A spatial view is created with the SDE commandline tools to 'register' the view with the geodatabase. The view won't show up in ArcCatalog unless it is registered. I have a data maintenance process that utilizes spatial views. After running unattended for close to a year, I am now having a problem where the spatial view becomes partially unregistered from the geodatabase. This prevents the geoprocessor from seeing the feature class. It also introduces the fairly common SDE issue where you can't re-create the view because the tool tells you that it already exists, but you can't delete it because the tool tells you that it doesn't exist. Has anyone else run into this issue and come up with either a prevention or remediation strategy. Platform: ArcSDE 9.3x, Oracle10g, UNIX.


title: Python - indentation error

text: I've started learning python for future use in QGIS and using Ubuntu, installed with Python 3. I'm doing the tutorials from the great book for newbies "A byte of python". I'm having troubles running some python-programs from the chapter "Control overflow", all with the same error. For example running this file of mine, "while.py" : Input:               number = 23     running = True     while running:     guess = int(input('Enter an integer : '))     if guess == number:     print('Congratulations, you guessed it.')     running = False # this causes the while loop to stop     elif guess < number:     print('No, it is a little higher than that.')     else:     print('No, it is a little lower than that.')     else:     print('The while loop is over.')     # Do anything else you want to do here      Output(the error):                 File "if.py", line 4         print('Congratulations, you guessed it.') # New block starts here             ^     IndentationError: expected an indented block     morck@morck-N150-N210-N220:~/dokumenter$ python3 while.py       File "while.py", line 4         guess = int(input('Enter an integer :'))             ^     IndentationError: expected an indented block     morck@morck-N150-N210-N220:~/dokumenter$      Where the right output according to the tutorial should have been:               $ python3 while.py     Enter an integer : 50     No, it is a little lower than that.     Enter an integer : 22     No, it is a little higher than that.     Enter an integer : 23     Congratulations, you guessed it.     The while loop is over.     Done      Can't figure out what causes the error and need some insight understanding this. Thank you


title: How to know if shapefile is georeferenced?

text: I have a shapefile which is located in province in indonesia. it contains many sub-province if i removed sub regions inside a shapefile, and it leaves a small region, and my question is that small region's scale is as same as a previous scale into a map how to reset a new scale that i want?


title: How to save a raster and vector into a single PNG file in GRASS?

text: I am trying to save a raster and vector into a PNG file. The vector is overlayed on top of the raster map. I was only able to save the raster file into a png using r.out.png but without the vector.


title: Looking for .Net alternative to Mapxtreme

text: I'm pretty new to the GIS area. I was involved in a project that used MapXtreme server and MapXtreme SDK to query the tabs (as in SearchNearest , SearchInRadius etc...) I'm looking for an open source alternative that I can integrate in .Net web app that will use OpenLayers + Tile Cache and I'll be able to query Tab (or Shape) files after the user pin the map. Thank you in advance


title: Extract certain area of contours/satellite images

text: I have 2 satellite images (jpeg2000 format) and 2 xyz-files from the same area. I am trying to extract certain region from the area including contour lines and the satellite image in order to import these into Google Sketchup for terrain modeling. I have been able to create contour lines from xyz-files. Anyway, I would like to first merge contours and then satellite pictures and finally extract only certain part of the data (both contours and satellite pictures). Is there some kind of instruction available that would help in a problem like this?


title: spgrass6 and GIS_LOCK

text: I would like to thank's everyone first, i learn a lot here (gis.stackexchange) about GIS. I have one question about the use of GRASS within R session, i'm able to do some basic task, but i have one important issue is to use GRASS graphics through spgrass6::execGRASS. I'm aware of the capabilities of R in graphics (lattice, ggplot and others), but i would like to be able to use GRASS graphics because of some option like nviz. Here is a minimal exemple               require(spgrass6)     require(spatstat)     # to use SpatialGridDataFrame in ggplot     fortify.SpatialGridDataFrame <- function(data) as.data.frame(as.im(as.image.SpatialGridDataFrame(data)))          # initialisation and the use of spearfish60 data     initGRASS(gisBase = "/usr/lib/grass64", home = tempdir(), gisDbase = "/home/dickoa/GIS/grass/grassdata/", location = "spearfish60", mapset = "usr1", override = TRUE)          # verif     gmeta6()          execGRASS("g.list", parameters = list(type = "vect"))          parseGRASS("d.mon")          # and this command gave me errors     execGRASS("d.mon", parameters = list(start="x0"))      and the error message is :               ERREUR :Unable to get GIS_LOCK enviroment variable value      And generaly, i really like the R workflow using emacs (ESS) and the idea of keeping everything in script, it's one of the main raison for me to use grass within R, so i want also to ask for other R and GRASS users if this integration through spgrass6 (Bivand 2011) can allow me to do everything as in interactive session of GRASS. here is my R session info (i use Kubuntu 10.10 as OS)               R version 2.12.2 (2011-02-25)     Platform: x86_64-pc-linux-gnu (64-bit)          locale:      [1] LC_CTYPE=fr_FR.UTF-8       LC_NUMERIC=C                    [3] LC_TIME=fr_FR.UTF-8        LC_COLLATE=fr_FR.UTF-8          [5] LC_MONETARY=C              LC_MESSAGES=fr_FR.UTF-8         [7] LC_PAPER=fr_FR.UTF-8       LC_NAME=C                       [9] LC_ADDRESS=C               LC_TELEPHONE=C                 [11] LC_MEASUREMENT=fr_FR.UTF-8 LC_IDENTIFICATION=C                 attached base packages:     [1] grid      stats     graphics  grDevices utils     datasets      [7] methods   base               other attached packages:      [1] spgrass6_0.6-26 XML_3.2-0       rgdal_0.6-33         [4] sp_0.9-79       spatstat_1.21-5 deldir_0.0-13        [7] mgcv_1.7-4      ggplot2_0.8.9   proto_0.3-9         [10] reshape_0.8.4   plyr_1.4        MASS_7.3-11              loaded via a namespace (and not attached):     [1] lattice_0.19-17    Matrix_0.999375-48 nlme_3.1-98


title: How to get an older version of OpenStreetMap data?

text: For some reason in the map of the Dominican Republic's Capital there are a bunch of streets which doesn't have any names. I had an OSM file which I downloaded ealier this year which did have those streets name. I really don't know WHY the names has just dissapeared, but I would like to know if there is a way that I can get them back. Maybe I can download an older version and then re-upload it so the current version get the streets' names back. So can you tell me how to get an older version of the OpenStreetMap or can you give me any other alternatives to get the streets' names back? Thanks in advance.


title: Export Map Schema\Metadata

text: I've been digging around, but I can't seem to find out if it's possible. I create a series of maps for a client and I need to document all the specifications in the map. So I need to create a spreadsheet/word doc of all the properties set in each layer (i.e. symbology colors, border thickness, definitions queries, label expressions, source, coordinate system...basically all the settings in the layers properties dialogue box) Same thing for the data frame properties settings. Is there a function in ArcMap to do this? Thanks, Mike


title: Polygon to Polyline tabular data transfer

text: I am working with a Divisional database. The division is divided by districts and then the district is further sub-divided by its micro units called sub- districts. Its topological relationship was error free and in the table the polygon data were well defined in its sub-district, district and divisional columns. Now I have drawn its road network throughout the division. There are some highway with unique name that covers more than 6-7 sub-district polygon. Now I want to create a road database that will break the road at the sub- district boundary line and join the polygon data (Sub-district, district and divisional data) to the clipped road with its existing database. I think it is possible but I couldn't find any solution. Can anybody help me?


title: Client Side Cache google map

text: I need to Cache the google map on the client side through the openlayer.I havent tried the client side cache yet and have very little idea about it too.


title: Road Topology Problem solution for junction

text: I am working with a Road Network. I have digitized the lines with AutoCAD Map 3D. Then Imported it to ArcMap. After entering all kind of database I have created Routing Dataset. But unfortunately its not giving the desired route. The most probably problem may be the junction. I think there are few problems with junction between the road. Its a huge database and it cant be check manually every junction. Is there any way to create a network topology to check the junctions easily whether the roads snapped or not? Thanks


title: Multiple vertices after editing giving wrong start and end x, y co-ordinates

text: I have a small problem which I have been fighting to solve but in vain. Initially I generated some polylines across the river at various locations. I found some of them are smaller in size and not fitting the river width, Therefore, I edited those using advanced editing tool by using extend polyline tool. However,I did not have knowledge that extending would create one more vertex in the middle. When I am computing distance using my script, the script is calculating distance between two earlier vertices it did not consider the another end vertices to compute the x,y of the actual end point of the modified lines. Please suggest me how to get solution. I tried to delete middle vertices still I find wrong distance. Think that wanna to rebuild/repair the geometry, I am not sure. Thanks and appreciate your help!


title: How to set up a Coverage Layer for Atlas Plugin?

text: Has anybody utilized the Atlas plugin? I am starting to take a look and setting it up on a project. My question is how do you set up the layer that is utilized for the coverage layer. Is there a way to get a polygon the mimics the view of the map area shown in the Print Composer?


title: Passing selected polygon to INTERSECTS filter in a WFS request

text: In the code below,                  <ogc:Filter>          <ogc:Intersects>            <ogc:PropertyName>the_geom</ogc:PropertyName>            <ogc:Function name="querySingle">               <ogc:Literal>sf:restricted</ogc:Literal>               <ogc:Literal>the_geom</ogc:Literal>               <ogc:Literal>cat = 3</ogc:Literal>            </ogc:Function>          </ogc:Intersects>        </ogc:Filter>      i'd like to replace                  <ogc:Literal>sf:restricted</ogc:Literal>        <ogc:Literal>the_geom</ogc:Literal>        <ogc:Literal>cat = 3</ogc:Literal>      with a polygon selected by a user. I get the polygon from a selection event. I've tried a number of options, but none is working, yet. Any assistance will be highly appreciated. My environment: GeoServer 2.1.3, OpenLayers 2.11. **UPDATE:** I've transformed my geometry to `EPSG:4326`, the one used by my layers in GeoServer:               var geometry = event.feature.geometry.transform(new OpenLayers.Projection("900913"), new OpenLayers.Projection("EPSG:4326"));      The event.feature.geometry represents a polygon selected by a user. I've also converted the geometry to WKT format as shown:               var wktParser = new OpenLayers.Format.WKT();     var feature = new OpenLayers.Feature.Vector();     feature.geometry = geometry;          var wktGeom = wktParser.write(feature);      I've then passed the geometry in `WKT` format to the filter as shown below:               '<ogc:Filter>' +     '<ogc:Intersects>' +         '<ogc:PropertyName>the_geom</ogc:PropertyName>' +         '<ogc:Function name="querySingle">' +             '<ogc:Literal>ws:layer_2</ogc:Literal>' +             '<ogc:Literal>the_geom</ogc:Literal>' +             '<ogc:Literal>INTERSECTS(the_geom,' + wktGeom + ')</ogc:Literal>' +         '</ogc:Function>' +     '</ogc:Intersects>' +     '</ogc:Filter>';      But when I send a request, the following error is returned: `<ows:ExceptionText>Process failed during execution java.lang.ClassCastException: org.geotools.filter.AttributeExpressionImpl cannot be cast to org.opengis.filter.Filter org.geotools.filter.AttributeExpressionImpl cannot be cast to org.opengis.filter.Filter</ows:ExceptionText>`. Any idea on what I may be doing wrong.


title: Importing .svg images to use as Symbols in ArcMap 10.1

text: I have several .svg files that I want to use as symbols in ArcMap 10.1. However, I do not see any import feature or the sort to use these .svg pictures. Thanks


title: How to get an GeoExt LegendPanel with image/png instead of image/gif?

text: I am using a `GeoExt.LegendPanel` to show legends for my various layers in _OpenLayers_. The `GeoExt.WMSLegend` supports configuring which image format to use, but the LegendPanel does not utilize this causing the WMSLegend to fallback to _image/gif_. So far I have manually edited the _GeoExt.js_ minified file to make the WMSLegend use _image/png_ as default instead of _image/gif_ , but I would prefer to parameterize this. **Am I missing some way to configure my LegendPanel to use _image/png_ as my WMSLegend image format?**


title: Algorithm for estimating position of many objects and differentiating those?

text: Let's assume a crowd-sourced tree location use case where a research team asks hikers to GPS the location of a specific tree type in case they encounter it. This will lead to a large number of longitude-latitude points each identifying the measured (noisy) GPS location of a tree. The researchers would now like to compute from those measurements a map of those trees. That means they have to differentiate location point estimates for close trees and after assigning the estimates to a specific tree calculate the location of that tree. Now trees don't move and GPS measurements can be pretty inacurate. What clustering algorithm would be the canoncial choice for this application? (If you have an R package::function to suggest please tell me)


title: How to get the coordinate system from a PostGIS database?

text: My question is probably amazinlgy stupid, but I can't find the answer myself. I would like to extract from a database the latitude and longitude of points. For now, I'm getting x(the_geom) and y(the_geom), and I would like to know what is the coordinate system of the results I get. How can I know this?


title: create POV-Ray image from DEM-file in grass

text: I am working on `r.out.pov` Tool of grass 6.4 creating a POV-ray(beautiful 3D) image from an elevation raster. ![enter image description here](http://i.stack.imgur.com/wQjiB.jpg) You can see this link for more information about 3D analysis of `GRASS`.   r.out.pov manual   Can anyone help me(step by step for windows) creating such an image of my DEM file? and some help about targa file (.tga).


title: Polygon country boundaries in decimal degrees

text: Using Google Maps, I'd like to create draggable polygons (on the fly) of country boundaries (they don't need to be particularly precise). Something similar to this, but for countries. So I want the boundaries in decimal degrees, but I can only find the data in .shp or kml/kmz files via sites like: shp:   * http://www.naturalearthdata.com/   * https://hiu.state.gov/data/data.aspx   * http://www.gadm.org/ KML:   * https://www.google.com/fusiontables/DataSource?docid=1N2LBk4JHwWpOY4d9fobIn27lfnZ5MDy-NoqqRpk ... but I'm unsure about how I go about extracting decimal coordinates from those data sources. Any advice on how to best solve this problem? Thank you.


title: Can we use procedures in a Query layer (ArcGIS Desktop 10.1 & Server 10.1)?

text: I have one basic question that.. Can we use procedures in Query layer (ArcGIS Desktop 10.1)?? Currently I am using query layer for ArcGIS Desktop 10.1 and its working fine . My DBA person told me that its difficult to modify query for each time.. As per discussion he told me that for large database its difficult to use query layer he suggested to use procedures.. I searched on Google and ArcGIS help..but there is no such information.. We are using Sql Server 2008 r2 , ArcGIS Desktop/Server 10.1.. So any help/links will be highly appreciated : )


title: What is the Best Practice for ArcGIS Server File Structure?

text: My organization is trying to decide how to best organize our file structure for services we put on our ArcGIS Server. As of now we have subfolders within the root so our structure is:               http://<servername>/arcgis/rest/services/<subfolder>/<servicename>/MapServer      However, after looking around most other organizations we have found have a structure of:               http://<servername>/arcgis/rest/services/<servicename>/MapServer      Essentially just putting all the services in the root. Is there a best practice on this? I can see reasons for both, and I can't really see an issue with subfolders, but the fact that most organizations put them in the root makes me wonder.


title: Caching / Preloading data for use with WMS queries

text: Hoping someone can help me with this issue. I've got a map implemented in my web application using OpenLayers that makes WMS queries to GeoServer that in turn queries a Postgres (PostGIS) database to produce WMS tiles. It also makes use of WMS GetFeatureInfo to get information about features being moused over. Because of the amount of data that is being queried in the initial load for the tiles can take anywhere between 1-10+ minutes depending on the data that is being queried (query can be adjusted by date/time ranges and a few other properties), which is acceptable because it's a lot of data. But I've got a couple of performance related issues that I'm trying to get around.   1. Every time the zoom level is changed the database has to be re-queried. Is there a way the data can be cached so that it doesn’t have to be re-queried for each change in zoom level? The data doesn’t change depending on zoom level so the query being run is returning the exact same result each time.    2. When a user mouses over a feature, a GetFeatureInfo request is made and a popup is shown over the map detailing that feature. It once again has to run the whole query again and then it narrows it down to the feature that has been moused over.  Previously before changing to using WMS we generated a single KML file from the data we got from our database query and then this would be used for loading the map at all zoom levels and getting the feature info. We had to change from this approach to deal with bigger and more complex polygons. So simply, is there anyway to get GeoServer to use a cache result of the query data or is there any way I could perhaps generate the KML file how I was before but then pass it to GeoServer to render it as WMS tiles? Would much appreciate any input on the matter or any suggestions to get around the problems I've stated. Thanks in advance, Mark


title: Why is the dijit.popup not working correctly?

text: I take the following steps:   * Click 2 overlapping features from 2 different featLayers.   * Popup auto-selects and highlights map.infoWindow.features[first object]   * The popup content is actually from map.infoWindow.features[second object] After you navigate to the second graphic and back, the popup corrects it's self?? here's an example using esri data: http://jsfiddle.net/92Q2G/1/ Any advice? Thanks.


title: Geoserver - layer from sql server

text: I am new to Geoserver (version 2.4.3). I have configured SQL store and then tried to create a Layer over one of my views with the correct coordinate system, when I tried to preview the layer using openlayers - I got this error in the logs.               Caused by: java.io.IOException         at org.geotools.jdbc.JDBCFeatureSource.getReaderInternal(JDBCFeatureSource.java:610)         at org.geotools.data.store.ContentFeatureSource.getReader(ContentFeatureSource.java:563)         at org.geotools.data.store.ContentFeatureCollection.features(ContentFeatureCollection.java:165)         ... 108 more     Caused by: com.microsoft.sqlserver.jdbc.SQLServerException: Incorrect syntax near the keyword 'FROM'.         at com.microsoft.sqlserver.jdbc.SQLServerException.makeFromDatabaseError(SQLServerException.java:197)         at com.microsoft.sqlserver.jdbc.SQLServerStatement.getNextResult(SQLServerStatement.java:1493)         at com.microsoft.sqlserver.jdbc.SQLServerStatement.doExecuteStatement(SQLServerStatement.java:775)         at com.microsoft.sqlserver.jdbc.SQLServerStatement$StmtExecCmd.doExecute(SQLServerStatement.java:676)         at com.microsoft.sqlserver.jdbc.TDSCommand.execute(IOBuffer.java:4575)         at com.microsoft.sqlserver.jdbc.SQLServerConnection.executeCommand(SQLServerConnection.java:1400)         at com.microsoft.sqlserver.jdbc.SQLServerStatement.executeCommand(SQLServerStatement.java:179)         at com.microsoft.sqlserver.jdbc.SQLServerStatement.executeStatement(SQLServerStatement.java:154)         at com.microsoft.sqlserver.jdbc.SQLServerStatement.executeQuery(SQLServerStatement.java:611)         at org.apache.commons.dbcp.DelegatingStatement.executeQuery(DelegatingStatement.java:208)         at org.apache.commons.dbcp.DelegatingStatement.executeQuery(DelegatingStatement.java:208)         at org.geotools.jdbc.JDBCFeatureReader.<init>(JDBCFeatureReader.java:140)         at org.geotools.jdbc.JDBCFeatureSource.getReaderInternal(JDBCFeatureSource.java:580)         ... 110 more     2013-12-29 12:49:52,626 INFO [geoserver.wms] -      Request: getServiceInfo     2013-12-29 12:49:52,632 ERROR [geotools.rendering] - Error querying layer     java.lang.Exception: Error transforming bbox         at org.geotools.renderer.lite.StreamingRenderer.getStyleQuery(StreamingRenderer.java:1068)         at org.geotools.renderer.lite.StreamingRenderer.processStylers(StreamingRenderer.java:2009)         at org.geotools.renderer.lite.StreamingRenderer.paint(StreamingRenderer.java:837)         at org.geoserver.wms.map.RenderedImageMapOutputFormat.produceMap(RenderedImageMapOutputFormat.java:481)         at org.geoserver.wms.map.RenderedImageMapOutputFormat.produceMap(RenderedImageMapOutputFormat.java:251)         at org.geoserver.wms.map.RenderedImageMapOutputFormat.produceM


title: change style of newly drawn points

text: I've a vector layer with a bunch of draw controls on it. You can add polygons, points and you can modify polygon. Now I changed the style of the points so they have a nice pin-like image. But them I noticed when I select the modify polygon tool and select a polygon, the handles on the polygon are also replaced by my point-image. How do I fix that?               this.panel = new OpenLayers.Control.Panel({             allowDepress: true,             autoActivate: false         });         this.controls = {             "polygon": new OpenLayers.Control.DrawFeature(this.layer, OpenLayers.Handler.Polygon, {                 "title": "New Zone"             }),             "POI": new OpenLayers.Control.DrawFeature(this.layer, OpenLayers.Handler.Point, {                 "title": "New POI",                 "displayClass": "olControlDrawPoi"             }),             "modify": new OpenLayers.Control.ModifyFeature(this.layer, {                 "standalone": false,                 "title": "Edit Zone"             }),             "clear": new OpenLayers.Control.Button({                 displayClass: "olControlClearLayer",                 trigger: myLocations.clearLayer,                 title: "Clearlayer"             })         };         for (key in this.controls) {             this.panel.addControls(this.controls[key]);         }              myMap.map.addControl(this.panel);      I tried adding the style in the Point Handler but that has no effect. I also tried to set a feature property so I can use a symbolizer in the style, but i can't figure out how to add this property before the feature is actually drawn. (so this works, but only after a layer.redraw(). So I want to have defautl styling for the polygons and its handels when its in modify mode, but I do want a nice marker image for my newly drawn points. Any ideas?


title: How to fill the area bounded by lines in QGIS?

text: There are two types of area I want to fill color:   * The areas bounded by a line, such as a building   * The area bounded by different lines, such as the river I don't know how with QGIS. I only know how to fill colors in polygons but bot lines .May you provide some suggestions? Thanks! The image here shows that the building and rivers are composed of lines: ![The image here shows that the building and rivers are composed of lines](http://i.stack.imgur.com/3yLVT.png)