title: ExtJS+GeoExt LayerTree (Map and Tree in separate views) Cannot get layers store to work

text: I have an Ext application i'm setting up and one of the page panels is a GeoExt map. The layerout of the map area is set up like this:               MapTabs      - Arctic Map Tab         - Map View         - Layer Tree View      - Antarctic Map Tab         - Map View         - Layer Tree View      So essentially, each tab on the MapTabs panel contains a map and a collapsible layer tree. Following is the code for the Arctic Tab. The key lines to the issue I am having are in ArcticMap.js and ArcticTree.js Basically ArcticMap.js is an ExtView (extended GeoExt.panel.Map) with id: 'ArcticMap'. In ArcticTree I need to get the layers object of the GeoExt map panel ('ArcticMap') and do so using the Ext compoonent query (see ArcticTree.js). When all is executed the result from the Ext component query looks identical to the result from the GeoExt2 LayerTree example but when I look at the page the layertree folders do not contain any layers. MapTabs.js               Ext.define('OPS.view.MapTabs.MapTabs' ,{              extend: 'Ext.tab.Panel',         alias: 'widget.maptabs',              requires: [             'OPS.view.MapTabs.ArcticTab.ArcticTab',             'OPS.view.MapTabs.AntarcticTab.AntarcticTab'         ],              initComponent: function() {                  this.items = [                 {                     xtype: 'arctictab',                     collapsible: false,                     header: false,                 },                 {                     xtype: 'antarctictab',                     collapsible: false,                     header: false,                 }             ]                  this.callParent(arguments);         }     });      ArcticTab.js               Ext.define('OPS.view.MapTabs.ArcticTab.ArcticTab' ,{              extend: 'Ext.Panel',              layout: 'border',         defaults: {             collapsible: false,             split: true,             bodyStyle: 'padding:0px'         },              alias: 'widget.arctictab',         title: 'Arctic',              requires: [             'OPS.view.MapTabs.ArcticTab.ArcticMap',             'OPS.view.MapTabs.ArcticTab.ArcticTree',         ],              initComponent: function() {                  // ARCTIC MAP TAB ITEMS             this.items = [                 {                     xtype: 'arcticmap',                     region: 'center',                     collapsible: false,                     header: false,                 },                 {                     xtype: 'arctictree',                     region: 'east',                     collapsible: true,                     width: 200,                 }                    ];                  this.callParent(arguments);         }     });      ArcticMap.js               // OPENLAYERS WMS URL     var arcticWms = OPS.Global.baseUrl.concat('/geoserver/arctic/wms');          // DEFINE AN EXT (INNER GEOEXT) OBJECT     Ext.define('OPS.view.MapTabs.ArcticTab.ArcticMap' ,{              id: 'ArcticMap',         extend: 'GeoExt.panel.Map',         alias: 'widget.arcticmap',         zoom: 3,         center: [118200, -1800000],         map: {             allOverlays: false,             projection: "EPSG:3413",             units: 'm',             maxExtent: new OpenLayers.Bounds(-8125549,-6101879,8186727,3197247),             controls: [                 new OpenLayers.Control.Navigation({dragPanOptions: {enableKinetic: true}}),                 new OpenLayers.Control.Zoom(),                 new OpenLayers.Control.MousePosition({prefix: '' +'EPSG:3413: '}),                 new OpenLayers.Control.ScaleLine(),                 new OpenLayers.Control.LayerSwitcher({'ascending':false}),                 //new OpenLayers.Control.Click({'active':true})             ],         },         layers: [                  // BASE LAYERS             new OpenLayers.Layer.WMS(                 "Natural Earth I",                 arcticWms,                 {layers: 'arctic:arctic_naturalearth'},                 {isBaseLayer:true}             ),                  new OpenLayers.Layer.WMS(                 "Landsat 7 Natural 90m",                 arcticWms,                 {layers: 'arctic:greenland_natural_90m'},                 {isBaseLayer:true}             ),                  new OpenLayers.Layer.WMS(                 "Landsat 7 B/W 45m",                 arcticWms,                 {layers: 'arctic:greenland_bw_45m'},                 {isBaseLayer:true}             ),                  // OVERLAYS             coastline = new OpenLayers.Layer.WMS(                 "Greenland Coastline",                 arcticWms,                 {layers: 'arctic:greenland_coastline',transparent:true},                 {isBaseLayer:false,visibility:false}             ),                  rdsFlightlines = new OpenLayers.Layer.WMS(                 "Radar Depth Sounder Flightlines",                 arcticWms,                 {layers: 'arctic:arctic_rds_line_paths',transparent:true},                 {isBaseLayer:false,visibility:true}             ),                  accumFlightlines = new OpenLayers.Layer.WMS(                 "Accumulation Radar Flightlines",                 arcticWms,                 {layers: 'arctic:arctic_accum_line_paths',transparent:true},                 {isBaseLayer:false,visibility:false}             ),                  snowFlightlines = new OpenLayers.Layer.WMS(                 "Snow Radar Flightlines",                 arcticWms,                 {layers: 'arctic:arctic_snow_line_paths',transparent:true},                 {isBaseLayer:false,visibility:false}             ),                  kubandFlightlines = new OpenLayers.Layer.WMS(                 "KUBand Radar Flightlines",                 arcticWms,                 {layers: 'arctic:arctic_kuband_line_paths',transparent:true},                 {isBaseLayer:false,visibility:false}             ),                  // VECTOR LAYERS             polygonBoundary = new OpenLayers.Layer.Vector(                 '',                 {displayInLayerSwitcher:false}             ),                  selectedLine = new OpenLayers.Layer.Vector(                 '',                 {displayInLayerSwitcher:false}             )              ]          });      ArcticTree.js               Ext.define('OPS.view.MapTabs.ArcticTab.ArcticTree' ,{              extend: 'GeoExt.tree.Panel',         alias: 'widget.arctictree',         title: 'Map Layer Selection',              initComponent: function() {                  this.rootVisible = false,             this.lines = false,                  this.store = Ext.create('Ext.data.TreeStore', {                 model: 'GeoExt.data.LayerTreeModel',                 root: {                          expanded: true,                     children: [                         {                             plugins: [                                 {                                     ptype: 'gx_layercontainer',                                     store: Ext.ComponentQuery.query('#ArcticMap').layers                                 }                             ],                             expanded: true                         },                          {                             plugins: ['gx_baselayercontainer'],                             expanded: true,                             text: "Reference Layers"                         },                          {                             plugins: ['gx_overlaylayercontainer'],                             expanded: true,                             text: "Data Layers"                         }                     ]                 }             });                  this.callParent(arguments);         }     });      Here is a preview of the portal (the layertree w/o layers is on the right): ![enter image description here](http://i.stack.imgur.com/f3AEa.jpg)


title: Issues with importing a netcdf file to arcgis as a point feature layer

text: I have a netcdf file of my computational grid. It is a structured curvilinear grid. When I make a feature layer from a it with the Multidimension Tool toolbox, I cannot get the point locations to plot properly. If I don't tell arcgis what netcdf dimensions to use (i.e. row dimensions option empty) all I get is a single point at the 0,0 origin. If I try to use the xi and yi dimensions, which are the indices of of the netcdf variable that I want to plot (i.e. with row dimensions as xi and yi) I get them to plot with the coordinates as the xi and yi indices instead of latitudes and longitudes. I have the set of lon - lat variables in my netcdf files, but I cannot select them to be used as the dimensions or the dimension values (the tool box will not see them at all, although they are in the netcdf file). How can I fix this issue? Using ArcGIs 10.0.


title: Edited features not appearing in ArcGIS

text: In ArcGIS 10.2 I created multi-segment lines using the Editor. Once a line was drawn, I selected "Finish Sketch" and then "Saved Edits". My new line appears in the attribute table but not on the map. However, if I double click said line in the attribute table, the line appears as selected in the map but disappears upon removing my selection. Why don't my edits appear permenantly on the map?


title: Not a valid data source of kml created from image overlay by Google Earth

text: I did image overlay in Google Earth and saved it as kml file. But it was showing me that it is not a valid data source when add it as vector(kml) in Qgis 1.8.0. Any solutions? Thanks!


title: Minimizing bandwidth while refreshing an ArcGIS dynamic map layer

text: I've been tasked with creating a web app where multiple users can draw on a map and have that reflected to other users in real time. Persistence of the drawing if users leave and come back is important. The server will be ArcGIS Server and the client will be done with the JavaScript API. I'm pretty sure the best way to do this is with a dynamic map service layer set with a short refresh time, but I'm worried about bandwidth usage. Are there any strategies for minimizing the volume of data sent such as having the server only send what's changed after the first query?


title: attribute table in qgis wont save custom values

text: I am trying to edit a point shapefile in qgis and I have editing on for that layer, I have added a new column for the custom data, I have set it to text and have 8 characters allowed in the column but every time I save and toggle editing off all of the values go back to null.


title: Is there a Select Feature by Attribute tool for Feature Classes when performing selection ModelBuilder?

text: I want to perform a Select by Attribute selection on a Feature Class in ModelBuilder, but I know the Select by Attribute tool does not work in this situation. I need to perform this function as I want to avoid creating an additional dataset. This task is the second to last one, the last one being a Delete Rows.


title: IObjectInspector and ClassExtensions

text: I have a little bit of dilema here. I have a hierachical collection of classExtensions that I very much would like to keep them as they are now: 1) Timestamper (master concrete class Extension) 1.1) BaseCadastralExtension (derives from Timestamper) 2) Derived classes (from BaseCadastralExtension) But I want to implement custom IObjectInspectors for each of the derived classes. I can't keep them deriving from BaseCadastralExtensions, because the IObjectInspector needs to derive from UserControl (am I right?). So, I cannot do this with multiple inheritance :P Is it possible to implement an IObjectInspector without trashing the class structure? Or do I need to rewrite? I wonder that composition in this matter (baseCadastralExtension to have a "funcionalities" member to take care of all the events, etc) and implement for each a custom IObjectInspector would do better? I'm not sure if I'm being clear, but I would like suggestions from you expert ESRI devs. Thanks, George


title: Corrupted rasters after import into ArcGIS FGDB

text: I have been seeing a number of corrupted rasters recently after importing them into file geodatabases (FGDB). The image below left shows a .tif prior to the import, while the image to the right shows the image after importing into the FGDB. Could these errors be due to NoData handling? What could be causing these errors? **Source TIFF info:**   * Columns and Rows: 5985, 7571   * Bands: 1   * Cellsize: 1,1 m   * uncompressed Size: 21.61 MB   * Source Type: continuous   * Pixel Type: Unsigned integer   * Pixel Depth: 4 Bit   * Number of classes: 4   * Storage Location: Google Drive ![enter image description here](http://i.stack.imgur.com/dlNhk.jpg)


title: FIELD CALCULATOR VB scripting

text: In my attribute tables field. I am trying to get rid of the "&" in between names i.e. "John Adams & wife" and want it to say "John Adams and Wife" in all the attributes with "&". What is the VB for that? I can't seem to find it in the help page of ArcGIS. Thanks


title: Plotting KML files onto a Google Map

text: I would like to show a map on my website. I have 2 different kinds of KML files that need to be plotted. One defines region boundaries. The other points within those boundaries. The region boundaries are in standard google map format. The points are in "EPSG:4326" format. Is google maps the right tool for the job or should I use something else? I heard openlayers is an alternative.


title: ARCGis LayerToKML Skips Feature Layers and Doesn't Process Layers in Group Layers

text: I am using ARCGis desktop and have written a quick conversion program to convert my layers to kml but I ran into some issues. It seems the only layers it can process correctly are group layers. When it processes the group layers I am missing all the feature layers included in the group layer in the resulting doc.kml. It's also failing on feature layers with an error that says it's not a layer with the following message: ERROR 000840: The value is not a Layer. When I know it's a layer, I can pull it up in ARCMap and, indeed, ARCMap agrees with me that it's a layer. I can use the MapToKML function without issue but I need to control my layers individually so I was hoping I can generate KML for each individual layer on a map. Thanks in advance for any help you can provide. Here is my source code:                 static void Main(string[] args)       {        IMapDocument mapDocument = new MapDocumentClass();        IPage page;             string mapFilename = @"c:\maps\BaseMap_NAD83_SD_StatePlane_Plumes_Scott.mxd";             Console.WriteLine("Initializing GEO...");        Geoprocessor GP = new Geoprocessor();             if (mapDocument.get_IsMapDocument(mapFilename) && mapDocument.get_IsPresent(mapFilename))        {         Console.WriteLine("Loading map...");         mapDocument.Open(mapFilename, "");              page = mapDocument.PageLayout.Page;              IMap map;         for (int iMapIndex = 0; iMapIndex < mapDocument.MapCount; iMapIndex++)         {          map = mapDocument.get_Map(iMapIndex);               Console.WriteLine("Converting layers for map {0}...", map.Name);               ILayer layer;          for (int iLayerIndex = 0; iLayerIndex < map.LayerCount; iLayerIndex++)          {           layer = map.get_Layer(iLayerIndex);                     ConvertLayerToKML(GP, layer as ILayer2);                Marshal.ReleaseComObject(layer);          }               Marshal.ReleaseComObject(map);         }        }             mapDocument.Close();        Marshal.ReleaseComObject(mapDocument);       }            private static string[] ConvertLayerToKML(Geoprocessor GP, ILayer2 layer)       {        string [] messages = new string[] {};             string layerName = layer.Name;             LayerToKML layerToKML = new LayerToKML();        layerToKML.layer = layer;        layerToKML.layer_output_scale = 10000;        layerToKML.out_kmz_file = string.Format("{0}.kmz", layerName);             Console.WriteLine("Converting layer {0}", layerName);             object result = GP.Execute(layerToKML, null);             messages = new string[GP.MessageCount];        for (int i = 0; i < GP.MessageCount; i++)        {         messages[i] = GP.GetMessage(i);        }             return messages;       }


title: datum conflict when creating buffer in arcgis 10

text: A buffer I created successfully executed, but I received the following message: "WARNING 000632: Datum conflict between map and output." How can I fix this?


title: ArcGIS 10.2 Restrict Map Service to specific IP Addresses

text: I've been asked if it is possible to restrict access to certain map services so that only a few whitelisted IP addresses are allowed to connect. What is the easiest and most straightforward way of accomplishing this? I was thinking about putting an .htaccess file in the service directory and allowing and blocking IPs with that, but I'm not sure what the actual service directory is located. (So as example I want to block http://devmymapservices.mydomain.local/dev/rest/my_map_service/MapServer?f=jsapi so that is only accessible by my IP address) how would I do that?


title: What is the inverse of QgsGeometry.asGeometryCollection() in qgis python?

text: What is the inverse operation of QgsGeometry.asGeometryCollection() as documented here: http://qgis.org/api/classQgsGeometry.html#a7cb1cbf8cfa441456067a66ad3a25a94 I'd like to do something like:               from qgis.core import (QgsFeature, QgsGeometry,                        QgsVectorLayer, QgsMapLayerRegistry,                        QgsField)     from PyQt4.QtCore import QVariant     from qgis.utils import iface     import math               def myFunction(arg):         return arg          layer=iface.mapCanvas().currentLayer()     feature = layer.selectedFeatures()[0]     geom = feature.geometry()     if geom.isMultipart():         newgeom=QgsGeometry.INVERSE([myFunction(part) for part in geom.asGeometryCollection ()])     else:         newgeom=myFunction(geom)     changeGeometry(feature.id(),newgeom)      Things that do not work:               # doesn't work:     newgeom=QgsGeometry.combine([part for part in geom.asGeometryCollection ()])     # doesn't work:     newgeom=([part for part in geom.asGeometryCollection ()])     # doesn't work:     newgeom = [part for part in geom.asGeometryCollection()]     # doesn't work: argument 1 has unexpected type 'list'     newgeom = QgsGeometry.fromPolyline([part for part in geom.asGeometryCollection()])     # QgsGeometry has no attribute 'MultiPolyline':      newgeom = QgsGeometry.MultiPolyline([part for part in geom.asGeometryCollection()])     # MuliPolyline is not defined:      newgeom = MultiPolyline([part for part in geom.asGeometryCollection()])     newgeom = QgsGeometry.fromMultiPolyline([part for part in geom.asGeometryCollection()])     newgeom = QgsGeometry.fromMultiPolyline([part.asPolyline for part in geom.asGeometryCollection()])          newgeom = QgsGeometry.fromMultiPolyline([part.asPolyline for part in geom])     newgeom = QgsGeometry.fromMultiPolyline([part.asPolyline for part in geom.asMultiPolyline()])      Works with Polylines / LineStrings:               newgeom = QgsGeometry.fromMultiPolyline([part for part in geom.asMultiPolyline()])          newgeom = QgsGeometry.fromMultiPolyline([userFunction(part,*args).asPolyline() for part in geom.asGeometryCollection()])      http://docs.qgis.org/testing/en/docs/pyqgis_developer_cookbook/geometry.html has some documentation on accessing and creating the various types. http://qgis.org/api/classQGis.html lists the enumerations of WKB* types and plain geometry types. It appears that each multipart type must be treated separately, and the multi-* types can't be disassembled and reassembled generally. In other words, you need to spell out and process the individual geometries separately like:               parts = ([userFunction(part,*args) for part in geom.asGeometryCollection()])          if geom.wkbType() == Qgis.WKBMultiLineString:          newgeom= QgsGeometry.fromMultiPolyline([x.asPolyline() for x in parts])     elif geom.wkbType() == QGis.WKBMultiPolygon:         newgeom= QgsGeometry.fromMultiPolygon([x.asPolygon() for x in parts])     elif geom.wkbType() == QGis.WKBLineString25D:         newgeom= QgsGeometry.fromMultiPolyline25D([x.asPolyline25D() for x in parts])     elif geom.wkbType() == QGis.WKBMultiPolygon25D:         newgeom= QgsGeometry.fromMultiPolygon25D([x.asPolygon25D() for x in parts])     ...          layer.changeGeometry(feature.id(),newgeom)      Is there a function or combination of functions that would simplify the above? Something like:               newgeom = CAST_AS_GEOMETRY_TYPE([part for part in geom.asGeometryCollection()],geom.WKBType())


title: Assigning field value based on location without join

text: I'm interested in updating a field in geodatabase feature classes based on their location in relation to another feature class in the geodatabase, but without having to perform a join (and therefore creating another feature class). I've got a polygon feature class ("boundaries") and then multiple point and line features that fall within those boundaries ("fieldData," "streams," "origins"). All of the features, once loaded to the database, have a common field labeled "siteCode" that I would like to populate programmatically rather than manually doing a Select by Location and calculating the fields (or creating extra files like performing a join). Better yet, I would like this to occur as the point and line files are loaded to the geodatabase, but that's not necessary. None of the files have this information stored in their attributes prior to being loaded to the database, but the siteCode is the first two letters in the name of each source shapefile. As it stands now, I select a feature in the "boundaries" feature class, Select by Location on each of the other files one at a time, and populate the fields using Calculate Field. I am open to any suggestions for workflow improvement here! My experience is limited with Python but I am interested in learning. If people respond with Python suggestions I would just ask that they include as much information as possible! I feel like this must have been answered a million times but I'm seeing a lot about using a join and not many ways to avoid it. Also: There is this posting that appears to be accomplishing something similar but I'm too much of a novice to understand exactly what it's doing and if it will do what I need. If someone can confirm that some variation on that code will work for me I will dive in and figure it out!                   #Make a feature layer for the neighourhoods     arcpy.MakeFeatureLayer_management(neighbourhood, "neighbourhoodLayer")     arcpy.MakeFeatureLayer_management(parcel, "parcelLayerTest1")     #Setup cursor to update new NBHD field     with arcpy.da.UpdateCursor(parcel, (parcelIDField,nbhdField)) as parRows:         for par in parRows:             # Create a query string for the current patrol zone              #parcelIDstring = par[0]             #queryString = '"' + str(parcelIDField) + '" =' + "'" + str(parcelIDstring) + "'"                  #Select parcels based on location contained within neighbourhood             arcpy.SelectLayerByLocation_management("parcelLayerTest1", "HAVE_THEIR_CENTER_IN", neighbourhood)                  nbhdTable = arcpy.SearchCursor(neighbourhood)             for nbhd in nbhdTable:                 nbhdValue = nbhd.getValue("Neighbourh")                 if nbhdValue:                     par[1] = nbhdValue                     parRows.updateRow(par)      Thanks in advance for any assistance.


title: What does Error 000885: Input Signature file does not have file extension indicate?

text: I can't run a Maximum Likelihood Classification because an Error 000885 says my input signature file does not have a file extension but it actually does as you can see in the picture. What am I doing wrong? ![enter image description here](http://i.stack.imgur.com/HYFah.jpg)


title: Web Editing Versioned Enterprise GDB

text: I would like to implement a web editing environment where edits are posted to a version of a database (delta tables) not too default. The idea behind this is so that any edits made through a web app would be looked over and posted to default by a GIS admin. I have tried to publish a version of the data, but I get a "Feature service requires a registered database" error (code: 00090).


title: Removing last point in polygon geometry using ArcPy?

text: This polygon was created from a grid of 500 x 500 meters, then Union, Generalize and Dissolve was used to arrive at this geometry, the laspoint not help me, I just want to have points on each vertex. How I can eliminate this point (lastpoint) that I can not enumerate my vertex. ![laspoint from polygon's geometry](http://i.stack.imgur.com/JvPDQ.gif) First thank you all for your attention to my problem. (Sorry for my bad English). I'm working with ArcGIS 10.1 and according to my workflow, it is performed as follows:   * From a grid of 500 x 500 meters, I select an area and generated a polygon with `arcpy.Union_analysis` as shown in the graphs.   * Using `arcpy.Dissolve_management` and `arcpy.Generalize_edit`, try to delete the intermediate points.   * From this polygon Figure **A** , generated a Feature (point type). With `arcpy.FeatureVerticesToPoints_management` to label then using the ID field and deploy the polygon to see a result like Figure **B**. Sometimes it works without problems, sometimes it fails. ![enter image description here](http://i.stack.imgur.com/nISTQ.gif) Thanks again to everyone for their comments and suggestions. Well, the first graph posted has a wrong result, and the second graph has a correct result, was only posted to show that sometimes it works and sometimes it does not and I do not know why. In fact, unnecessary points that show the outcome of this third graph, I can remove the very fast and easy way to "Start Editing" ArcGIS, but the fact is that I am creating automated processes with python, then put them as services in ArcGIS Server. I have tried using many methods like `Dissolve, Generalize, SimplifyPolygon_cartography, geometry Modify, Repair Geometry`, I tried to recreate the polygon convert to line `UnsplitLine_management`. The most important thing for me is to label the vertices of the polygon with the corresponding geometry (would like to use some expression), as not doing so, I convert the polygon to get points and the geometry of these points as a label. ![enter image description here](http://i.stack.imgur.com/VUcnO.gif)


title: Get the all coordinates (points ) inside a polygon boundery

text: I work on the following map in arcgis10.2 and i get vertices in bounderies ![enter image description here](http://i.stack.imgur.com/i7cGQ.jpg) I need to get all coordinates (points )inside this polygon(map).How can i get it ?


title: UpdateCursor 'is not defined' Python

text: I have several files (shapefiles). Table 1 has only one row with several fields whereas Table 2 has several rows. What I want is to take a value from Table 1 and populate all rows in Table 2. I used a SearchCursor to find the value I want and an UpdateCursor to populate the rows.               import arcpy          arcpy.env.workspace = "C:\\Users\\Documents\\Scoring\\Ecodis_clip_split"      #set workspace to folder with values          rows = arcpy.SearchCursor ("Table1.shp", "", "", "area_ha", "")      #search for the value in attributetable          arcpy.env.workspace = "C:\\Users\\Documents\\Scoring\\Split_wetlandt"      #set workspace to folder with files which need to be updated          rows = UpdateCursor ("Table2.shp", "", "", "area_eco", "")      #Update all rows of file with value from SearchCursor above      Running this, I keep getting this Errormessage:               Traceback (most recent call last):     File "<module1>", line 17, in <module>     NameError: name 'UpdateCursor' is not defined      If something is not defined, this means that Python sees it as a variable? Even if I could get rid of the errormessage: would it populate all the rows in Table 2 with the value or do I need to do something else? Note that I just started to learn Python and never scripted anything what so ever... Edit 6/17/14 11:32 am @ Aaron, would this replace the values in Table 2 with my searched values from table 1?               import arcpy     arcpy.env.workspace = "C:\\Users\\Documents\\Scoring\\Ecodis_clip_split"      #set workspace to folder with values          rows = arcpy.SearchCursor (table 1", "", "", "area_ha", "")      #search for the value in attributetable          arcpy.env.workspace = "C:\\Users\\Documents\\Scoring\\Split_wetlandt"      #set workspace to folder with files which need to be updated          with arcpy.UpdateCursor (Table 2 ["area_eco"])      #Update all rows of file with value from SearchCursor above     for row in cursor:         area_eco.replace(SearchCursor)      EDIT 6/27614 12:04pm I updated my script to this. Somehow it deleted my whole attributetable in Table 2 but I got it restored. I run it again and now nothing happens. It does not update the rows, the value in the rows stays 0.               import arcpy     arcpy.env.workspace = "C:\\Users\\Documents\\Scoring\\Ecodis_clip_split"      #set workspace to folder with values          rows = arcpy.SearchCursor ("Table1", "", "", "area_ha", "")      #search for the value in attributetable          arcpy.env.workspace = "C:\\Users\\Documents\\Scoring\\Wetland"      #set workspace to folder with files which need to be updated          rows = arcpy.UpdateCursor ("Table2", "", "", "area_ha", "")     del rows


title: Openlayers Select Feature: avoiding the on/off aspect of selecting a feature: want to be able to click, then click again

text: Map here: http://tinyurl.com/potjxlu I can click an orange icon, popup appears, hit the X button to exit out of the popup.That is OK. But if I want to bring that popup again for the same point, I have to click it twice, not once. Also, sometimes the X button in the popup becomes disabled. It's not a consistent error so not sure what's going on there. Looking at the OL docs for SelectFeature, I thought `highlightOnly` may be the solution but it is not. The current code for my select feature control:               var selector = new OpenLayers.Control.SelectFeature([wfs, vectors_origin, vectors_destination], {                 clickout: false,                 toggle: true,                 hover: false,                 autoActivate:true         });      Any input on how to achieve my desired functionality is appreciated. Edit: If I just use the default vectors (the little orange circles), what's basically happening is: -I'll click the orange circle, the orange circle becomes highlighted with a dark purple: the popup now appears. -The popup has an X button. I click this to exit out and that aspect works. -Now at this point: the orange circle is still highlighted the dark purple: in other words it is still selected. -So what happens is I have to click WITHIN the orange circle (currently highlighted dark purple) to UNSELECT it. Now that dark purple has turned back to the initial orange.  -Now I can click it again and the popup should appear correctly: which it does. -So how do I skip this middle step? Note that my issue is only for points that you want to click again: if I click the Austin, TX feature then the Los Angeles, CA feature, there are no issues with the popups. But I basically want to cut out this middle deselection step as it seems unncessary for what I want to achieve.


title: Error adding a delimited text layer

text: Im working with QGIS 1.8.0 and I have been trying to upload a delimited text layer but keep getting the following error : The layer file:///C:/Users/fernando/Desktop/ATLAS/TAMAULIPAS/Estaciones%20Meteorol%C3%B3gicas%20Tamaulipas%20SMN.csv?delimiter=%5B,%5D&delimiterType=regexp&xField=X&yField=Y is not a valid layer and can not be added to the map I attach you a picture of the settings Im using when uploading the CSV file. Do you have any ideas why could it be wrong? Thanks in advance. ![enter image description here](http://i.stack.imgur.com/a9fkq.png)


title: Is it possible to automatically update a Buffer FeatureClass when creating new features?

text: Is there a way I can update the results of my Buffer analysis? I have a roads FeatureClass, and I've built a buffer for it. I will still make edits to the road FeatureClass. Let's say I add a new segment of road, how can I automatically buffer the new segment without creating a new buffer FeatureClass? I use ArcGIS 10.1 Thank you for your help. PS: Apparently when using Geomedia if your information changes in your data source, so does the analysis in the workspace. I want to know if ArcMap can do the same.


title: How to get JSON response from submitJob operation?

text: I'm trying to get the JSON response of the submitJob operation launched by a restful call: http:///submitJob. I'm using Java7. The response I'm getting is an HTML response, using the code below:               URL url = new URL("http://<gp-task-url>/submitJob");          HttpURLConnection connection = (HttpURLConnection) url             .openConnection();     connection.setRequestProperty("Content-Type", "application/json");      connection.setRequestMethod("GET");          InputStream is = connection.getInputStream();      Do you have any idea on how to retrieve a JSON response?


title: How to add a MSSQL spatial layer using Python?

text: I'm having difficulty adding a MSSQL spatial data layer using Python. The database is password protected, and I can add the layer without any problems using the built in feature, but I can't make it work with a Python script. I have been trying to follow the advice on this website without any success. Using the following script (except with the actual variables obviously) gives me a 'Invalid layer' error every time and I have tried with various different syntax:               uri = "MSSQL:server=xx.xx.xx.xx,database='MYDBNAME',user='MYUSERNAME',PWD='MYPASSWORD',tables=dbo.MYTABLENAME"     qgis.utils.iface.addVectorLayer(uri,'newlayer','ogr')      Any help would be greatly appreciated. I'd also like to be able to run SQL queries on this as well so if anyone has some advice that would be great. The advice on the qgis website explains how to load just about every other type of layer except MMSQL (http://www.qgis.org/en/docs/pyqgis_developer_cookbook/loadlayer.html)


title: OpenLayers doesn't show <Point>'s label

text:
I want to display a KML file via OpenLayers using the example given at this page. Displaying a KML File on a Vector Layer in OpenLayers For your information I've used the same code samples shown there to get the results discussed here. An excerpt of the KML file that I want to be displayed is:                 <Placemark>       <name>Tropical Depression EMONG</name>       <description>
Actual Position

June 17, 2013 02:00 AM

Source: PAGASA
       </description>              <Point>       <coordinates>127.5,13.5</coordinates>       </Point>       </Placemark>      where "Tropical Depression EMONG" is the label. The fact is that when the KML file is opened in Google Earth all labels are shown as they should, whereas when I use OpenLayers code only the icons are shown without the associated label. Is there a limitation in OpenLayers that won't displayd POINT elements, or do I need to provide special parameters to make the labels to show up? Thanks for your help. Robert Shoferbridge


title: how to project a flat image to spherical mercator EPGS:900013 to use in a Image Layer in openlayers

text: I have a flat 400x400 pixels image which represents an area of 400x400 kilometers. I need to draw this image as an overlay on my OpenLayers mapping system. I use OpenStreetMap/Google Maps as base layers. I want to draw the overlay image using an Image Layer and I want to center my overlay image at a specific lat/lon point. The projection of the base layer is EPSG:900913 which is, in fact, the projection used by OpenStreetMap and Google Maps. What I have understood (please correct me if I'm wrong), is that I can't use a flat, 'unprojected' image as an Image Layer, but I have to do a sort of 'reprojecting' to adjust the image to be compliant with the projection of the map, which is EPSG:900913. Searching here and the net, it seems I should use a library/tool called GDAL to do this, but whereas I am completely new to this I ask you:   * Can anyone address me on how to use it?   * Could you possibly tell me if what I have just said is correct and the way I described here is valid?   * How could I use gdal to translate my flat unprojected 400x400 pixels image to an image projected in EPSG:900913, and suitable to be put as an Image Layer as an overlay? Thanks.


title: QGIS contour plugin exception struck help

text: I am a fisheries biologist and I am just getting started with QGIS. I took a few GIS classes in college so my GIS knowledge is not extremely high. On to my question. I am trying to make a depth map of one of the lakes that I manage. I recorded my data using a lowrance depth finder. I have imported the data into QGIS and now have all the points that were recorded with depth and Lat. and Long. When I try to use the contour plugin I get an exception struck #### error. The error displays a number but when I delete that line from the attribute table the number changes. I keep deleting the number of the error that is returned but it seems like an endless loop. Any help with what I am doing wrong would be Greatly Appreciated. Data: ![Example of data](http://i.stack.imgur.com/9mdcD.png) Error: ![This is the error I am getting.](http://i.stack.imgur.com/fyGLP.png) I can go through and manually delete points and after deleting a certain amount I no longer get the exception struck error. But by the time I do this I have lost many points and the map will no longer be accurate.


title: Cell Size of Raster Explained

text: I'm not quite understanding the cell size values that a given raster has. I been doing some kriging interpolation on some very dense point data and if I leave the cell size parameter blank, it's cell size is like an extremely small number, like e-006. What does this extremely small cell size number mean? I always that the cell size of a raster is the size of the cell that each value of the raster represents spatially. But with such a tiny number, doesn't seem to make sense. When running the Kriging tool and if I enter a whole number, like 1 for the cell size, it just errors out. If I zoom into the raster and measure the size of the cell pixels, they look like there around 3 ft by 3 ft in size. Am I confusing cell size with pixel size here? If so, what determines the size of the pixels?


title: Clip in QGIS 2.4.0 resulting in empty output files

text: I'm seeing the same problem to Clip or Intersection tool - QGIS giving empty files trying to use the Clip tool since upgrading to QGIS 2.4.0, but have not been able to fix the problem based on the responses there. I tried creating a simple example to ensure was not my data but am seeing the same problem. I have a series of points & a polygon. I want to retain those points within the polygon only. ![Example of points and polygon](http://i.stack.imgur.com/vznTz.jpg) I use Vector -> Geoprocessing Tools -> Clip ![Commands used](http://i.stack.imgur.com/w8LKm.jpg) This results in an empty attribute table. Options I've tried:   * I have tried switching the Input and Clip layer in case I had these round the wrong way    * Disable on the fly reprojections (confirmed the layers do overlap, they are in the same projection)   * Checked through Properties to confirm both layers have CRS EPSG:4326 - WGS 84   * I left this simple example running for approximately 40 minutes with no change in outcome in case I was just not leaving enough time (in previous versions this was successful in a short time span).   * I have tried using other Geoprocessing tools to see if it is just clip - I do not receive any output with any of them.


title: What does 8GK100K mean?

text: I have a lot of mosaic datasets:               Zone8GK100k, Zone9GK100k, ... , Zone7GK25k, ... , Zone7Gk50K, ... , Zone7GK200K.      As I know, 7GK, 8GK... are utm zone numbers. 1) What does 25K, 50k, 100k, 200k mean? 2) What is the relation between these numbers and the ArcGis scales (1:1 000, 1:10 000, 1:24 000, 1:100 000, 1:250 000 ...)? _**Update:_** Thanks. I have layer Zone9GK200k. So, I set: "In beyond 1:200 000". What minimum scale(out beyond) for this layer do I set?


title: How to remove white parts in ArcMap 10.1?

text: This is the part of my map. ![enter image description here](http://i.stack.imgur.com/Mp25N.jpg) On this picture you can see white strips. Some of this white parts are appeared on the intersection of 2 Mosaic Dataset (for eaxample, there is the real map in 1 Mosaic dataset, there is the empty scan in 2 Mosaic dataset). **How to set option "no-data" of Mosaic Layer?** I want to display map without white parts. _**UPDATE_** On this picture you can see different MosaicDatasets. On the other forum people tell me that it is impossible to delete white parts in this case. ( **because it is not a single mosaic dataset** ) Is it true?


title: What is the projection of WGS84 in GIS software?

text: WGS84 is a reference ellipsoid that determines coordinates of every point on Earth using latitude, longitude and height about the its surface. This ellipsoid is often used as a base to make maps. Though, to make a map on a flat paper or a flat digital screen you need a projection too. An example of a popular projection would be UTM. Mainstream GIS software give plethora of projection alternatives to display your data layers. Among them, they also let you display your data using WGS84 only. My question is: how is that possible? Of course they are using a certain projection but what projection is that? Why they do not tell the user about that?


title: GRASS in QGIS not working (windows XP)

text: I am unable to use GRASS algorithms in QGIS 2.0.1 (on windows xp). I've already configured the GRASS path. The algorithms seems to work correctly (I cannot see errors in the log), however in the end the following windows appears: > Oooops! The following output layers could not be open Output layer: > C:\DOCUME~1\MCORDI~1\IMPOST~1\Temp\processing\ebc399ac72364b68aad151d232a1379a\output.shp > The above files could not be opened, which probably indicates that they were > not correctly produced by the executed algorithm. [...] It seems that GRASS > is not correctly installed and configured in your system. Please install it > before running GRASS algorithms. Here below the log for the **v.distance** command > Algorithm v.distance - Finds the nearest element in vector map 'to' for > elements in vector map 'from'. started g.proj -c proj4="+proj=utm +zone=32 > +ellps=intl +towgs84=-87,-98,-121,0,0,0,0 +units=m +no_defs" v.in.ogr > min_area=0.0001 snap=-1 dsn="E:/Tavolo Tecnico AUSL-Pr/GIS" > layer=asma_DOM_select output=tmp13811305945511 --overwrite -o v.in.ogr > min_area=0.0001 snap=-1 dsn="E:/GIS DataBase" > layer=eurostreets_EmiliaRomagna_FRC04 output=tmp13811305945512 \--overwrite > -o g.region n=4999266.578 s=4848836.96069 e=801201.169359 w=516043.246024 > res=1 v.distance from=tmp13811305945511 to=tmp13811305945512 upload=cat > column=mindist output=output6f83c8f8869a4a41bdbc595d37c12891 --overwrite > v.out.ogr -e input=output6f83c8f8869a4a41bdbc595d37c12891 > dsn="C:\DOCUME~1\MCORDI~1\IMPOST~1\Temp\processing\ebc399ac72364b68aad151d232a1379a" > format=ESRI_Shapefile olayer=output type=auto > > E:\PROGRA~1\QGISDU~1\bin>set HOME=C:\Documents and Settings\mcordioli > > E:\PROGRA~1\QGISDU~1\bin>set GISRC=C:\Documents and > Settings\mcordioli.qgis2\processing\processing.gisrc > > E:\PROGRA~1\QGISDU~1\bin>set GRASS_SH=E:\Programmi Lavoro\Quantum GIS > Lisboa\apps\msys\bin\sh.exe > > E:\PROGRA~1\QGISDU~1\bin>set PATH=E:\Programmi Lavoro\Quantum GIS > Lisboa\apps\msys\bin;E:\Programmi Lavoro\Quantum GIS > Lisboa\apps\msys\lib;E:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\bin;E:\PROGRA~1\QGISDU~1\apps\msys\bin;E:\PROGRA~1\QGISDU~1\bin;C:\WINDOWS\system32;C:\WINDOWS;C:\WINDOWS\System32\Wbem;E:\PROGRA~1\QGISDU~1\apps\msys\bin;E:\PROGRA~1\QGISDU~1\apps\Python27\Scripts;E:\PROGRA~1\QGISDU~1\apps\qgis\bin;E:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\lib > > E:\PROGRA~1\QGISDU~1\bin>set WINGISBASE=E:\Programmi Lavoro\GRASS GIS 6.4.3 > > E:\PROGRA~1\QGISDU~1\bin>set GISBASE=E:\Programmi Lavoro\GRASS GIS 6.4.3 > > E:\PROGRA~1\QGISDU~1\bin>set GRASS_PROJSHARE=E:\Programmi Lavoro\GRASS GIS > 6.4.3\share\proj > > E:\PROGRA~1\QGISDU~1\bin>set GRASS_MESSAGE_FORMAT=gui > > E:\PROGRA~1\QGISDU~1\bin>if "" == "" set PATH=E:\Programmi Lavoro\GRASS GIS > 6.4.3\bin;E:\Programmi Lavoro\GRASS GIS 6.4.3\lib;E:\Programmi > Lavoro\Quantum GIS Lisboa\apps\msys\bin;E:\Programmi Lavoro\Quantum GIS > Lisboa\apps\msys\lib;E:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\bin;E:\PROGRA~1\QGISDU~1\apps\msys\bin;E:\PROGRA~1\QGISDU~1\bin;C:\WINDOWS\system32;C:\WINDOWS;C:\WINDOWS\System32\Wbem;E:\PROGRA~1\QGISDU~1\apps\msys\bin;E:\PROGRA~1\QGISDU~1\apps\Python27\Scripts;E:\PROGRA~1\QGISDU~1\apps\qgis\bin;E:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\lib > > E:\PROGRA~1\QGISDU~1\bin>if not "" == "" set PATH=E:\Programmi Lavoro\GRASS > GIS 6.4.3\bin;E:\Programmi Lavoro\GRASS GIS 6.4.3\lib;;E:\Programmi > Lavoro\GRASS GIS 6.4.3\bin;E:\Programmi Lavoro\GRASS GIS > 6.4.3\lib;E:\Programmi Lavoro\Quantum GIS Lisboa\apps\msys\bin;E:\Programmi > Lavoro\Quantum GIS > Lisboa\apps\msys\lib;E:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\bin;E:\PROGRA~1\QGISDU~1\apps\msys\bin;E:\PROGRA~1\QGISDU~1\bin;C:\WINDOWS\system32;C:\WINDOWS;C:\WINDOWS\System32\Wbem;E:\PROGRA~1\QGISDU~1\apps\msys\bin;E:\PROGRA~1\QGISDU~1\apps\Python27\Scripts;E:\PROGRA~1\QGISDU~1\apps\qgis\bin;E:\PROGRA~1\QGISDU~1\apps\grass\grass-6.4.3\lib > > E:\PROGRA~1\QGISDU~1\bin>set GRASS_VERSION=6.4.0 > > E:\PROGRA~1\QGISDU~1\bin>if not "" == "" goto langset > > E:\PROGRA~1\QGISDU~1\bin>FOR /F "usebackq delims==" %i IN (`"E:\Programmi > Lavoro\GRASS GIS 6.4.3\etc\winlocale"`) DO @set LANG=%i > > E:\PROGRA~1\QGISDU~1\bin>set > PATHEXT=.COM;.EXE;.BAT;.CMD;.VBS;.VBE;.JS;.JSE;.WSF;.WSH;.PY > > E:\PROGRA~1\QGISDU~1\bin>set PYTHONPATH=;E:\Programmi Lavoro\GRASS GIS > 6.4.3\etc\python;E:\Programmi Lavoro\GRASS GIS 6.4.3\etc\wxpython\n > > E:\PROGRA~1\QGISDU~1\bin>g.gisenv.exe set="MAPSET=PERMANENT" > > E:\PROGRA~1\QGISDU~1\bin>g.gisenv.exe set="LOCATION=temp_location" > > E:\PROGRA~1\QGISDU~1\bin>g.gisenv.exe set="LOCATION_NAME=temp_location" > > E:\PROGRA~1\QGISDU~1\bin>g.gisenv.exe > set="GISDBASE=C:\DOCUME~1\MCORDI~1\IMPOST~1\Temp\processing\grassdata" > > E:\PROGRA~1\QGISDU~1\bin>g.gisenv.exe set="GRASS_GUI=text" > > E:\PROGRA~1\QGISDU~1\bin>g.proj -c proj4="+proj=utm +zone=32 +ellps=intl > +towgs84=-87,-98,-121,0,0,0,0 +units=m +no_defs" > > E:\PROGRA~1\QGISDU~1\bin>v.in.ogr min_area=0.0001 snap=-1 dsn="E:/Tavolo > Tecnico AUSL-Pr/GIS" layer=asma_DOM_select output=tmp13811305945511 > --overwrite -o > > E:\PROGRA~1\QGISDU~1\bin>v.in.ogr min_area=0.0001 snap=-1 dsn="E:/GIS > DataBase" layer=eurostreets_EmiliaRomagna_FRC04 output=tmp13811305945512 > --overwrite -o > > E:\PROGRA~1\QGISDU~1\bin>g.region n=4999266.578 s=4848836.96069 > e=801201.169359 w=516043.246024 res=1 > > E:\PROGRA~1\QGISDU~1\bin>v.distance from=tmp13811305945511 > to=tmp13811305945512 upload=cat column=mindist > output=output6f83c8f8869a4a41bdbc595d37c12891 --overwrite > > E:\PROGRA~1\QGISDU~1\bin>v.out.ogr -e > input=output6f83c8f8869a4a41bdbc595d37c12891 > dsn="C:\DOCUME~1\MCORDI~1\IMPOST~1\Temp\processing\ebc399ac72364b68aad151d232a1379a" > format=ESRI_Shapefile olayer=output type=auto > > E:\PROGRA~1\QGISDU~1\bin>exit Converting outputs Loading resulting layers > Algorithm v.distance - Finds the nearest element in vector map 'to' for > elements in vector map 'from'. finished


title: QGIS 2.0 Atlas Generation: Can I force the selected feature to be the last one rendered

text: I have a layer with a number of features that represent zones which share some boundaries with others. I am using Atlas generator to produce a set of maps, one for each zone extent. The zones are rendered in different colours and I want the current feature rendered last so it is on top. That is in the red zone map I want the outline of the red zone (rendered in red ;) drawn last so that _all_ the red zone is red and not blue where the blue and red share a boundary. Alternatively can I suppress the other features altogether?


title: How to measure areas using an official CRS which is not available in QGIS?

text: I am very new to QGIS and doing a project on Glacial Lake Outburst & Flooding (GLOF). I am trying to find the area of a lake using images from Google Earth 7.0.3.8542. I am using the WGS84 CRS from my Google Earth image for the area calculation. I am getting some messy area values like 0.000103145e-5. I can't even guess in what unit it's being displayed. I have an official Local Cordinate System for our country in .prg format but I can't find it in the list of available CRSs in QGIS 2.0 or how to load it as a custom CRS. **How can I load the local CRS in my computer into QGIS?**. The contents of the .prg file are: > PROJCS["DRUKREF 03 local grid > GASA",GEOGCS["GCS_DRUKREF03",DATUM["D_DRUKREF03",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",250000.0],PARAMETER["False_Northing",-2500000.0],PARAMETER["Central_Meridian",90.0333333333333],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]


title: How do I import Bitmap symbol images to an annotation feature class in Arc GIS 10.1?

text: I am currently working on some map layers and I am trying to annotate a buildings layer with BMP symbols ( exits, fire hydrant, manhole etc..). So far I have created an annotation feature class but have been unable to import these BMP symbols, which I have saved in a folder to the annotation symbol collection.


title: selecting polygons that share XX percents from another polygon (all are in the same layer)

text: I have a layer with many overlapping polygons. I want to select / delet all polygons that another single polygon contain at least 90% (or any other value) of them. I'm using ArcGIS 10 (can get access to other programs)


title: problem in spline with barriers in ArcGIS 9.3?

text: I am trying to use the script tool "Spline with Barriers" from arcgis 9.3 3D Analyst toolboxes. I am not getting expected output. The output is just only one value over the surface. We are using ArcGIS 9.3 in XP Remote Desktop and it is a German Version of Windows, so the system represents the decimal place as a comma. Is it possible to edit this tool to store the field values as 3.40 instead of 3,30?


title: What are license options to share and distribute ArcGIS ModelBuilder toolkit?

text: Do you have expertise or experience to share on what software license can be used to share a model created in ArcGIS ModelBuilder to ensure that it is free to use, or used for non commercial purposes. I am currently looking at the GPL v3 license, but since the model is based on using ArcGIS software I am not sure if this is appropriate. UPDATE I should have added for clarity. We wish that although the Models are free to use that there is the ability to restrict commercial use (which may then be under a different license) We wish to be able to ensure that any edits or improvements to the Models are also made available and free to use.


title: Draw line between points from multiple datasets

text: I have 2 point datasets and want a tool or script to draw a line between all points in one dataset to all the points in the 2nd dataset. Is it possible? I have tried ESRI "Point Distance" and "Generate near Table" but these just give me a table of points. I need to also calculate the distance between the points and this is between MGA Zone 55S and 54S ![enter image description here](http://i.stack.imgur.com/M7QJa.jpg) I checked How to write a script that will automatically draw lines from one point to multiple destination points? (ArcMap 9.3) and tried the arcgis 10 tool but it gives me write errors in ArcGIS 10.2 An ESRI or QGIS solution is fine. \-- The MMQGIS looks perfect but I have a python error. I am using Ver 2.0.1 on a Win 7 Enterprise.               An error has occured while executing Python code:          Traceback (most recent call last):       File "C:\Users\georgec/.qgis2/python/plugins\mmqgis\mmqgis_menu.py", line 299, in hub_lines         dialog = mmqgis_hub_lines_dialog(self.iface)       File "C:\Users\georgec/.qgis2/python/plugins\mmqgis\mmqgis_dialogs.py", line 1397, in __init__         self.set_hub_attributes(self.hublayer.currentText())       File "C:\Users\georgec/.qgis2/python/plugins\mmqgis\mmqgis_dialogs.py", line 1411, in set_hub_attributes         for index, field in enumerate(layer.dataProvider().fields()):     AttributeError: 'NoneType' object has no attribute 'dataProvider'          Python version:     2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)]               QGIS version:     2.0.1-Dufour Dufour, d94c044          Python path: ['C:/OSGEO4~1/apps/qgis/./python/plugins\\processing', 'C:/OSGEO4~1/apps/qgis/./python', 'C:\\Users\\georgec/.qgis2/python', 'C:\\Users\\georgec/.qgis2/python/plugins', 'C:/OSGEO4~1/apps/qgis/./python/plugins', 'C:\\INTERG~1\\ERDASI~1\\usr\\lib\\Win32Release\\python', 'C:\\OSGEO4~1\\bin\\python27.zip', 'C:\\OSGEO4~1\\apps\\Python27\\DLLs', 'C:\\OSGEO4~1\\apps\\Python27\\lib', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\plat-win', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\lib-tk', 'C:\\OSGEO4~1\\bin', 'C:\\OSGEO4~1\\apps\\Python27', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\PIL', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\win32', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\win32\\lib', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\Pythonwin', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\Shapely-1.2.18-py2.7-win-amd64.egg', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\six-1.3.0-py2.7.egg', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\xlrd-0.9.2-py2.7.egg', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\xlwt-0.7.5-py2.7.egg', 'C:\\OSGEO4~1\\apps\\qgis\\python\\plugins\\fTools\\tools', 'C:\\Users\\georgec\\.qgis2\\python\\plugins\\mmqgis/forms']      I have moved the projection issue to Calculate length of line across zones with at least m accuracy


title: I have a WKT How can i put it to MAPINFOW.prj file

text: I have this WKT file:               PROJCS["PRS92 / Philippines zone 3"     ,GEOGCS["PRS92"     ,DATUM["D_Phillipine_Reference_System_1992"     ,SPHEROID["Clarke_1866",6378206.4,294.9786982138982]]     ,PRIMEM["Greenwich",0]     ,UNIT["Degree",0.017453292519943295]]     ,PROJECTION["Transverse_Mercator"]     ,PARAMETER["latitude_of_origin",0]     ,PARAMETER["central_meridian",121]     ,PARAMETER["scale_factor",0.99995]     ,PARAMETER["false_easting",500000]     ,PARAMETER["false_northing",0]     ,UNIT["Meter",1]]      How can I use it in mapinfo .prj file so that i can use the CRS? I'm very new in this kind of matter. I hope someone can help me. Cheers,


title: QGIS select features from attribute table using the "LIKE" operator

text: I created an expression with the help of the GUI in the QGIS attribute table of a vector (shapefile) layer of Native American ("Indian") reservations in the USA state of New Mexico. "NAME" LIKE 'Picuris' No row became selected. The full "NAME" in the desired row is 'Picuris Indian Reservation' . I don't mind typing the whole name, but in principle, is not the LIKE function supposed to find partial matches? That is, if I want 'State Hwy 68' (which is in the FULLNAME field) (of a different layer), will the LIKE operator not find it if I use "FULLNAME" LIKE 'Hwy 68' ?


title: Is rid=2 significant in PostGIS raster structure?

text: It seems that in most, if not all, examples with raster appear to include the `WHERE rid=2` statement in the query. http://postgis.net/docs/RT_ST_Value.html Is there something about the rast structure that requires this, or is this just an artifact of previous examples used to reduce the sample dataset?


title: Is it possible to cluster markers based on city name in Google Maps API

text: I have a requirement to display all the customers in map. Since the customer list is too large it looks messy when I place markers. I have tried the marker cluster functionality but it groups based on distance and that's not what I prefer. Is there a way to cluster markers based on city names or custom logic?? Thanks in advance


title: nearest facility with pgrouting

text: I have a postgres table of students with lat, lon and geom, , and a table of schools with the same, and a table of roads that has been converted to nodes. I'm trying to find the driving distance to the nearest school for each student. But I'm new to postgres and pgrouting, but I'm not seeing a great way to structure the multiple joins to work with the pgr_drivingDistance function. Ideally I would end up with a new table that would have the student id, the id of the closest school via driving distance, and the distance. Any thoughts on how to approach this would be great. Update addressing Jendrusk's questions:   1. what do you mean by converting roads to nodes? Maybe you meant network graph? select pgr_createTopology('roads', 0.0001, 'geom', 'id')   1. Are you able to find at least one way and the trouble is with many ways? The challenge is finding the shortest way for multiple students. I could take the lat/lon for one student, find the nearest node and then calculate distances that way, but I need to do it for hundreds.   1. Nearest school euclidean or by the roads? Driving distance.   1. What was the work-flow of this conversion form roads to nodes? See above.   1. What are the names of tables and column names in this tables? Three tables: Kids, roads, schools. Kids and schools have a GID, lat, lon, geom, and a bunch of other non-spatial info. Roads has GID, geom, source, target and length.   1. How big are your tables with students and schools? 100k students, several dozen schools.   1. Is it OK to select nearest school euclidean No, interested in driving distance. > > > Also function st_DrivingDistance is not finding route, but it's creating > set of edges reachable in specified time (e.g. how far you can go in 15 > minutes). >> >> Also 2 you need distance from student to nearest school but which is nearest if you don't know the distance? Do you have some kind of relation between schools and students? If not you have to count distance from one student to every school to know which is nearest. In most cases nearest euclidean (in straight line) will be also nearest by the roads, but if you have some special cases in your graph such as river with one bridge you cant be sure.<<< Again, I do not need to find the route. I need to identify the closest school based on driving distance. My impression is that using pgr_drivingDistance it is possible to create a table containing driving distance length from a specified node to all other nodes in a network. You can scope this operation to nodes within a particular distance of the index node, and obviously, in my case, only consider nodes that are the closest node to a school. I'm basing this info on the following post: http://anitagraser.com/2011/05/13/catchment-areas-with-pgrouting- driving_distance/ I believe pgr_drivingDistance can limit the query based on a number of different costs, including driving time and distance, depending on how the cost is specified both in the query and in the data.


title: How to narrow the result from Nominatim with only cities and countries, no additional areas with the same name?

text: I'm trying to get WKT or GeoJson polygons of cities and countries from Nominatim (for example, Paris and France). I send a `GET` request but it returns me several results named "Paris" or "France". How may I narrow results only to a capital and a country. Here's my request url:               http://nominatim.openstreetmap.org/search/fr/paris?format=json&polygon_geojson=1      And I get an array of 10 objects where only one is the capital:               class: "boundary", type: "administrative".      The rest are:               class: "place", type: "isolated_dwelling"     class: "place", type: "hamlet",     class: "place", type: "neighbourhood"     and so on...      How to select only the desired city and only the desired country ignoring rest of results with the same name?


title: Hyperlink broken path - python check

text: I've recently come across an issue where some of our hyperlinks aren't working. I'm familiar with python and have made similar scripts before, but I would first like to ask for your assistance to maybe something already pre- existing or another tool to use. I would like to cycle through all the hyperlink paths for a certain shapefile and search the C:\Scans\ for the path. If the path doesn't exist than I want the path to be written to an excel file. Can anyone help me with this? It would be greatly appreciated, thanks! Final code so far               import arcpy     import os          fc = r"G:\GIS\Lauren\Shapes\watermains.shp"     fields = ["Truck_Path"]          with arcpy.da.SearchCursor (fc,fields) as cursor:         for row in cursor:             Roll = str(row[0])             if os.path.exists(Roll):                 pass             else:                 f = open(r"C:\Scans\BrokenLinks.txt","a")                 f.write(Roll + os.linesep)                 f.close()


title: jQuery tabs into OpenLayers switcher layers

text: How can I put jQuery UI tabs in openlayers layerswitcher? I can not imagine how this can be done ...


title: How to zoom to the extent of Qgsvectolayer qgis?

text: Anyone know how to zoom to the extent of a Vectorlayer in qgis using python plugin?


title: Calculate x,y,z distance between 3D point and multipatch closed polygon in ArcScene

text: I have a multipatch polygon and a pointZ FeatureClass containing a single point. I would like to calculate the distance in X,Y,Z from the point to the closest point on the surface of the multipatch polygon. From reading the documentation, I cannot figure out if the ArcScene "3D Near" function is the right tool or not. Has anyone else tried something similar? What tool did you end up using?


title: Difference between Heat Map and Hot Spot Analysis?

text: What is the difference between "Heat Map Analysis" and "Hot Spot Analysis"? Are we perhaps using hot spot analysis to get the required relations between features, and then using the results as an input for a heat map?


title: GPS tracking in ArcGIS Engine application

text: Is there an ArcGIS api for ArcObjects /.Net that I can use to build a function in an ArcGIS Engine 10 application that allows me to track a GPS signal and display the movement plus lets me export and import GPS-tracks?


title: How to make GRASS and SAGA tools capable to use in QGIS？

text: As my title, I installed Qgis, GRASS and GDAL by osGeo4w package, and as I finish the installation, I want to use GRASS tool, I face this problem： ![enter image description here](http://i.stack.imgur.com/m1Ttz.png) But I already installed GRASS and SAGA, how to make those tools in QGIS to be used？


Is it possible to render POIs/landmarks from Postgres/postgis DB using google Map like icons? Any readily available icon sets?


title: How to use a polygon grid to calculate statistics on line and point features

text: What seems a simple task but I can't find the right tool to do it. I have a polygon grid cells and say a point feature class. I want to be able to create a table that contains the grid cell ID and for example the sum of points within that cell or the average of the value for those points. For a line feature I would like to simply have a 0/1 value if the line falls into a grid cell. I am doing this in ArcMap 10. This is like Zonal statistics but on a vector layer instead of a raster.


title: How to disable all mouse event in openlayers' map div?

text: version: OpenLayers-2.13.1 ## brower : IE8 (I Hit it!!! But I need to work on it) I create a div with id 'map' to display openlayer map. Sometimes, I need disable map div catch any mouse events, let me enable "drag" items into map (not really drag and drop event). But I don't know how to disable the map receives mouse events, so, I dynamic create a new div into html body to cover above the map, make the new div receive all mouse events and without pass these events to the map, by this way, Google Chrome works well, but in IE8 is not!!!! So..... sad!! In IE8, the created new div only receive mousemove event at itself's border, when mouse inner blank area of the new div, the mousemove events be received by the map div, upper level's new div receive nothing. In Chrome, when mouse move onto the new div, no matter on the border or inner blank area, all mouse event pass to upper level's new div first, the map div would not receive! Could anybody told me how to make openlayers map div avoid receive mouse events and how to restore receive mouse events ability? If use new div cover upper the map, how to let upper div receive all mouse events in IE? Thank You very much! \---TingSky---


title: Can tilecache create UTFgrids?

text: Am i able to create utf grids in tilecache or do i have to move to tilestache? I ask because i am using, and would prefer to continue to use, utmZ10 and not wgs84 or spherical mercator.


title: pure OpenLayer is good or GeoExt

text: pure **OpenLayer** is good or **GeoExt** WHich have more options w.r.t GIS client-side


title: Can ArcGIS make a join with selected attributes being added into the outcome?

text: My question is when making a spatial join, is it possible to only assign selected attributes of source layer to the target layer instead of getting all attributes from source layer? I'm asking is because each time when I do a join, I need to clean up the table to remove useless attributes. Thanks!


title: osm2pgrouting: Import multiple OSM files?

text: I want to import multiple OSM data files for a few different countries into postgres using osm2pgrouting. Can I just run osm2pgrouting on each file with the same database settings (e.g., import countryA.osm, then countryB.osm, etc.) or would this cause problems?


title: Cause of ImportError: No module named mapscript in QGIS mapserver export?

text: I'm trying to expor project to mapfile in new QGIS Chugiak but a get this error below, someone have the same error? How can I fix this?thnks.               An error has occured while executing Python code:               Traceback (most recent call last):       File "C:/Users/Artur/.qgis2/python/plugins\rt_mapserver_exporter\plugin.py", line 65, in run         from .mapfileexportdlg import MapfileExportDlg       File "C:/PROGRA~1/QGISCH~1/apps/qgis/./python\qgis\utils.py", line 454, in _import         mod = _builtin_import(name, globals, locals, fromlist, level)       File "C:/Users/Artur/.qgis2/python/plugins\rt_mapserver_exporter\mapfileexportdlg.py", line 30, in          import mapscript       File "C:/PROGRA~1/QGISCH~1/apps/qgis/./python\qgis\utils.py", line 454, in _import         mod = _builtin_import(name, globals, locals, fromlist, level)     ImportError: No module named mapscript               Python version:     2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)]               QGIS version:     2.4.0-Chugiak Chugiak, 8fdd08a          Python path: ['C:/PROGRA~1/QGISCH~1/apps/qgis/./python/plugins\\processing', 'C:/PROGRA~1/QGISCH~1/apps/qgis/./python', u'C:/Users/Artur/.qgis2/python', u'C:/Users/Artur/.qgis2/python/plugins', 'C:/PROGRA~1/QGISCH~1/apps/qgis/./python/plugins', 'C:\\INTERG~1\\ERDASI~1\\usr\\lib\\Win32Release\\python', 'C:\\PROGRA~1\\QGISCH~1\\bin\\python27.zip', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\DLLs', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\plat-win', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\lib-tk', 'C:\\PROGRA~1\\QGISCH~1\\bin', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\PIL', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\jinja2-2.7.2-py2.7.egg', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\markupsafe-0.23-py2.7-win-amd64.egg', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\pytz-2012j-py2.7.egg', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\win32', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\win32\\lib', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\Pythonwin', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\Shapely-1.2.18-py2.7-win-amd64.egg', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\six-1.3.0-py2.7.egg', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\xlrd-0.9.2-py2.7.egg', 'C:\\PROGRA~1\\QGISCH~1\\apps\\Python27\\lib\\site-packages\\xlwt-0.7.5-py2.7.egg', 'C:\\PROGRA~1\\QGISCH~1\\apps\\qgis\\python\\plugins\\fTools\\tools', 'C:/Users/Artur/Documents/GIS DataBase']


title: Lambert93 to WGS84, IGN algorithms, longitude problem

text: I'm looking for a Java-based solution to convert coordinates from Lambert93 to WGS84. I found this, "lambert-java", which is an implementation of the algorithms (PDF) published by IGN (French government agency in charge of geographic information). However, when testing it with data from North of France, the resulting longitude is off by something like 50km. Here's an example:   * Source data (in Lambert93): X = 668832.5384, Y = 6950138.7285   * Geofree conversion to WGS84: lat = 49.64961, lng = 2.56865 (this is correct)   * Output of the above Java program: lat = 49.65262, lng = 1.90503 (this is incorrect) As far as I can tell, the program is "faithful" in that it faithfully implements the algorithms from IGN; so the fault either lies in the algorithms themselves, or in some parameter or some constant that needs to be adjusted. How can one sort this out? (There are other Java-based solutions capable of doing all sorts of GIS manipulations, but since I only need to do this conversion (on lots of data), I'd rather use a very simple, very light package, and "lambert-java" fits this bill nicely; however of course it also needs to be correct...) Edit 2014-03: this was solved by the original author five months ago and has been working wonderfully ever since.


title: Why do valid polygons repeat the same start and end point?

text: In the GIS world, for most modern formats and Standards (e.g., Shapefiles, WKB/WKT, GML, KML, etc.), _valid_ polygons need to have closed linear rings, which is a linestring of coordinates where the first point is a repeat of the last point. For example, a triangle requires _four_ points (not three).               POLYGON ((10 20, 30 60, 50 20, 10 20))      ![Triangle](http://i.stack.imgur.com/BuiqS.png) Who started this convention, and why? Is it some legacy of pre-Shapefile storage? (like how MS Windows still uses 2-byte CR+LF newlines?) Other non-GIS standards (e.g. SVG) _don't_ require this repetition to encode polygons.


title: Profile plot trouble in Grass GIS

text: I'm new to Grass and GIS environments in general, and I recently faced a problem when trying to get the profile plot from a DTM. When I draw the transect, instead of the plot, I get the following message:               Traceback (most recent call last):       File "C:\Program Files (x86)\GRASS     6.4.2\etc\wxpython\gui_modules\profile.py", line 468, in     OnCreateProfile          self.DrawPlot()       File "C:\Program Files (x86)\GRASS     6.4.2\etc\wxpython\gui_modules\profile.py", line 490, in     DrawPlot          marker=self.properties['marker']['type'])       File "C:\Program Files (x86)\GRASS 6.4.2\Python27\lib     \site-packages\wx-2.8-msw-unicode\wx\lib\plot.py", line 342,     in __init__          PolyPoints.__init__(self, points, attr)       File "C:\Program Files (x86)\GRASS 6.4.2\Python27\lib     \site-packages\wx-2.8-msw-unicode\wx\lib\plot.py", line 145,     in __init__          self._points =     _Numeric.array(points).astype(_Numeric.Float64)     ValueError     :     could not convert string to float: *      Does anyone out there have an idea about my problem? Thanks in advance for your support!


title: File IMG from OSM

text: Is it possible to convert or open an IMG file from OSM? I don't find a solution for the visualization and conversion. Thanks


title: ArcMap is exporting to extremely large pdf files

text: My organization has one user who recently has noticed very large pdf files being exported from her mxds. When I recreate and do the export, I am getting a 4mb file (normal). When she exports the file we are seeing a 2.4 GB file. I changed some settings and was able to get it down to 1GB, but that wasn't satisfactory (obviously). Things we have recently tried: -Deleted the normal.mxt (reduced pdf size to 100mb) -Turned off the setting that saves the spatial information in the pdf -DPI is at 300 It isn't an overly large mxd. Normal vector data with some labeling. The pdf size is 42" x 60". Any ideas on what to try next? She is on an ArcEditor 9.3.1 license.


title: QGIS : copy/paste with "Avoid Intersection" option

text: I need to copy and paste polygons from one layer to another but with adjustments in the pasted geometries: I need them to be cut in the same way they would have been cut by a single edit with the "Avoid Intersection" option in the project snapping options. In Arcmap I can nearly do this with the "construct feature" tool (topology toolbar) user the target layer existing objects. But in my QGIS project I need the use a third layer to cut the pasted polygons. Tanks


title: How applications like Global Mapper, ArcMap and ArcScene support showing high resolution geotiff images?

text: Considering my question in Game Developement: How to create textures of big images in opengl (bigger than the MAX_TEXTURE_SIZE) it seems that, there's no way to show an image bigger than the MAX_TEXTURE_SIZE in OpenGL (I'm not sure yet, looking forward to more answers in Game Developement) In my case MAX_TEXTURE_SIZE = 8192 but my image is a 16997x15931 which can be shown in openGL based on my current knowledge about this API. Just tobe curious, and maybe they're using a toolkit or API which is not free **What way programs like Global Mapper, ArcMap and ArcScene manage to show high resolution geotiff images for example bigger than 15000x15000 with high performance?** **What toolkit, API, 3rd party library do they use for`display tasks` and `2D or 3D rendering` in their programs?** I highly doubt that image-pyramids is the concept related to my question?   if it's true beyond answering my question about APIs please introduce me a link in GDAL documentation that has an excellent tutorial about pyramids?


title: How to create sidebar help information for custom arcgis script tools?

text: I have been creating custom gis tools with python via the arcpy module. I can set up the tool successfully with parameters and all, but I cannot figure out how to create help text that shows up on the sidebar of the tool like the standard arcgis tools have.


title: Layer package duplication

text: I create a layer package from desktop 10 and the double click to decompress it on my machine (to see what the customer is going to see).   I get the initial folder (same as my lpk file) and then another one duplicated inside. c:\documents and settings\users\myusername\My Documents\ArcGIS\Packages\mylpkfile2-17-11\mylpkfile2-17-11   There are data files in both folder structures (with shape files in both).   There is a lyr file in the root folder that points to the first dataset.   I see a lyr for v10 in the v10 folder.   However this lyr uses the same data as the other lyr.   ![Layer structure](http://i.stack.imgur.com/v7Wes.png) It looks to me that there is no reason for this data structure.   I don't know if it is just my machine or ... Does anyone else see this also? The above occurs if I double click on the lpk. If I run the extract package tool from within arcmap it extracts only one folder structure to the location I instruct it to.


title: postprocessing GPS coordinates with RINEX data to improve precision?

text: Im doing some mapping for OSM in my region and I am curious if the precision of tracks from a GPS Receiver can be improved if they are post processed with RINEX data. The tracking is done with a Garmin GPSMAP 60CSX or alternatively my Android Phone and OSMTracker for Android. The RINEX data is available from http://igs.bkg.bund.de/filestatus/dailyrinex cheerio steve


title: Using python lookup table to compare records with existing geodatabases

text: I am trying to use python to look through an access table containing facilities that need to have geodatabases with geodatabases that already exist in order to generate a list of facility geodatabases that still need to be created. The script is working but the loop is stopping at the first record of the table and not continuing to page through the rest of the records. Here is what I have:               lookuptbl = r"P:\\USACE Mobile 0610-01 2010 63D RSC NR Contract\\03-Source Documents\\GIS\\SITE_ID_LOOKUP.mdb\\new_facility"     q = '[include_in_gis] = ' + "'" + "X" + "'"               rows = gp.SearchCursor(lookuptbl, q)     row = rows.next()          while row:         facil_id = row.GetValue("facil_id")         lookuptbl2 = r"P:\\USACE Mobile 0610-01 2010 63D RSC NR Contract\\03-Source     Documents\\GIS\\SITE_ID_LOOKUP.mdb\\gdbs"         q2 = '[facil_id_in_gdbs] = ' + "'" + facil_id + "'"         #IF q is satisfied, the gdb that is currently being processed was found in lookup table;         #If q is NOT satisfied, the script will skip the gdb alltogether; nothing will be done              rows = gp.SearchCursor(lookuptbl2, q2)         row = rows.next()              while row:             print facil_id + "--already exists"             row = rows.next()                   row = rows.next()              I'm new to scripts so any help would be appreciated. Thanks!


title: How do I check if Mapserver is installed or not?

text: I have installed MS4W (full) and choose port 8080 because on port 80 i have already installed IIS7. Now installation of MS4W is successful , but how i check ? I put http://localhost:8080/ in my browser, but i could not see the page. any idea how to fix it ?


title: How to use python to turn a .dbf into a shapefile

text: I have been scouring the internet trying to find a pythonic (sp?!) way to process this data.. Everyday we will recieve a load of data in .dbf format (hopefully) - we then need to save this data as a shapefile. Does anyone have any links or any suggestions as to my process?


title: In OpenLayers zoom to extent of all overlays

text: In Openlayers I can zoom to the bounding box of all layers like this:               map.zoomToExtent();      I can also zoom to the bounding box of a single layer:               map.zoomToExtent(layer_name.getDataExtent());      Is there an easy way of zooming to the extent of all layers except the base layer? In other words showing the extent of all overlays.


title: How to Convert a Number to Percentage in Field Calculator

text: Using ArcGIS 10.2 and Personal geodatabase,I have a field with values of Numbers like **0.461538** and Type of " **Double** " can you please let me know how I can Export them to new Field in Percentage Format to be shown like 46% or at least round the number to looks like 46 then I can add "%" on label to them (Honestly! I do not think this is a good idea for further Analysis!). Thanks


title: Viewshed calculation using SQL

text: Can i calculate viewshed in MS-SQL server? Viewshed analysis is Raster operation. I have elevation value vector observation as point vectors. Are there any methods to calculate Viewshed (using algorithm or something like that) ?


title: Adding and Editing Excel table to the data frame

text: 1. I have few excel tables of almost 5-6 columns and i need to paste them at different locations in the same map.   2. The values in the excel table gets revised every quarter and i need to edit the table. Please don't confuse it with the attribute table. I have attached the way it look in the arcmap and the excel table. Help me in figuring out a solution for the same. Suggest me way by which i dont have to copy/paste the table each time. Is there a way to edit the table![enter image description here](http://i.stack.imgur.com/LYj34.png) ![enter image description here](http://i.stack.imgur.com/PsN1o.png)


title: Does Y mean latitude and X mean longitude in every GIS software?

text: Now it's been 3 days I am using Mapinfo and it has Y as latitude and X as longitude. Is that the same case for all mapping software? As for any country their respective value is multiple of 1 or -1. So for Nepal can I say it is on positive side +1 for both latitude and longitude? And for USA to be +1 Y and -1 X.


title: How to change QGIS project color space

text: Is it possible to change the color model of an entire QGIS project from RGB to CYMK? I know that it's possible to change the color models of individual symbols/features (through the style manager), but I can't seem to find an option for the whole project. While my workflow usually sees me polishing my maps in Inkscape/Illustrator (where it's easy to switch between RGB/CMYK), it would be nice to have the option in QGIS as well.


title: How can you insert a custom symbol?

text: How can you use a custom symbol for either dot density or proportional symbol maps in QGIS? Is there a way to import a symbol of your own such as a drawing or image that you can create or download from somewhere else?


title: Zooming on point cluster in openlayers

text: I have clustering enabled on an OpenLayers map and have used code from the first answer in: Openlayers zoom in on cluster to zoom in on a cluster.               function onFeatureClick(event) {         if (event.feature.cluster.length > 1) {             // *** Change this so that the map zooms on cluster location ***             map.setCenter(event.feature.geometry.getBounds().getCenterLonLat());             map.zoomIn();         } else {             // popup code         }     }      **Q:** Is it possible to zoom on the cluster location so that it stays in the same place? Right now the cluster gets centered on the map. In other words the cluster should not move when zoomed in.


title: Why ArcGIS turn Front Rd to Frontage Rd

text: I did a google search and a search in the ArcGIS forum. Does anyone know why when ArcGIS Desktop standardize an address such as "285 LAKE FRONT RD" and turn it into "Frontage Road"? It's a fairly small problem but it really shouldn't happen. Right? Edit: I'm using ArcGIS Desktop but I'm not sure why that would make a difference. I now realize that I should say that when I was standardizing the address field, I was using the US Single House locator. The problem IS that "Lake Front Rd" is converted to "Frontage Road" after using the Standardize Addresses tool. Frontage is not recognized by USPS as a street suffix.


title: Most reliable way to import large dataset with osm2psql

text: I know my question is similar to slow import via osm2pgsql to postgresql database and Optimizing osm2pgsql imports for OSM data but as we currently have quite some problems importing a large set of OSM data I open a new one... What is the best way to import a large dataset (OSM Export of Europe) into a postgres DB Our computer has 32 GB of RAM... so it could use all of that :) We tried a couple of params, but had no success... last try we used the               osm2pgsql -c -S /usr/share/osm2pgsql/default.style --slim -d osm-europe -U postgres -C 25000 europe-latest.osm.pbf      But we ran out of memory even though our sever has 32 GB of RAM available.               pending_ways failed: out of memory for query result     (7)     Error occurred, cleaning up      Any help on how to improve out import command is highly appreciated! Even if it takes longer... But we need to have the data imported into our postgres DB. Would you recommend using an EC2 for the task or should our setup work with different parameters? Thanks


title: Network Analyst in ArcGIS: Multiple vehicles, multiple trips, finite vehicles in a vehicle routing problem

text: I am working on a problem where I have 30 vehicles who can all make at most 6 stops in a day. They have to service more than 500 shops. I can get Network Analyst to do a Vehicle Routing Problem easily, but I can't figure out how to make the vehicles do multiple trips. The Network Analyst only lets me do one trip for each vehicle so that they only serve 180 shops. How can I get ArcGIS to allow me to do multiple trips? I am interested in distance of routes. I can measure this if I create enough vehicles to visit all the shops, but it is extremely tedious to create 75+ vehicles by hand ("Add Item") and will be impossible in my target case of more than 2000 shops. Please help! Thank you.


title: Unserialize socklist error when using rasterEngine in r

text: Basically, I am running the "rasterEngine" function in R. I want to use parallelization to change values in a raster image (I have around 400 images that I need to work on). However, the rasterEngine from the spatial.tools package throws an error and does not execute on my machine. Below I will outline the code, system info, and error message.               Sys.info()     sysname release version nodename      "Windows" "7 x64" "build 7601, Service Pack 1" "CDV0267"      machine  login    user   effective_user      "x86-64"  "x"  "x"   "x"       This the code I am employing. I took the function from the rasterEngine tutorial on this website (http://publish.illinois.edu/jgrn/software-and- datasets/rasterengine-tutorial/). I also used a function that I wrote myself but I gave me the same error:                require(spatial.tools)      require(doParallel)           # Begin a parallel cluster and register it with foreach:      cpus = 2 # The number of nodes/cores to use in the cluster      cl <- makeCluster(spec = cpus, type = 'PSOCK', methods = FALSE)      # Register the cluster with foreach:      registerDoParallel(cl)           name_r <- "C://test//tcv.tif"      start_raster <- raster(x=name_r, values = T)           meters_to_feet <- function(elevation_in_meters,...)       {      feet <- elevation_in_meters * 3.28084      return(feet)      }           tahoe_lidar_highesthit_feet <- rasterEngine(elevation_in_meters=start_raster,     fun=meters_to_feet)      Here is the error message, and the traceback:               Error in unserialize(socklist[[n]]) : error reading from connection     In addition: Warning messages:     1: closing unused connection 4 (<-CDV0267.CDV.ORG:11840)      2: closing unused connection 3 (<-CDV0267.CDV.ORG:11840)           > traceback()     12: unserialize(socklist[[n]])     11: recvOneData.SOCKcluster(cl)     10: recvOneData(cl)     9: recvOneResult(cl)     8: dynamicClusterApply(cl, fun, length(x), argfun)     7: clusterApplyLB(cl, argsList, evalWrapper)     6: e$fun(obj, substitute(ex), parent.frame(), e$data)     5: foreach(chunkID = chunkID, .packages = c("rgdal", "raster", "spatial.tools",        "mmap", .packages), .verbose = verbose) %dopar%      focal_hpc_pixelChunkFunction(chunkID,        tr, x, chunk_format, fun, args, layer_names, outbands, filename)     4: focal_hpc_pixel_processing(tr, chunkArgs)     3: focal_hpc(x, fun = focal_hpc_multiRaster_function, args = c(list(fun = fun,        debugmode = debugmode), args), window_dims = window_dims,        window_center = window_center, filename = filename, overwrite = overwrite,        outformat = outformat, additional_header = additional_header,        processing_unit = processing_unit, chunk_format = chunk_format,        minblocks = minblocks, blocksize = blocksize, outbands = outbands,        outfiles = outfiles, setMinMax = setMinMax, debugmode = debugmode,        .packages = .packages, verbose = verbose)     2: rasterEngine(elevation_in_meters = start_raster, fun = meters_to_feet)     1: system.time(tahoe_lidar_highesthit_feet <- rasterEngine(elevation_in_meters =    start_raster,        fun = meters_to_feet))      I also turned off as much of the firewall as I can -- the computer is still inside of the company's network (so that might be the problem). Does anyone know what to do? Why does this process fail? Thank you.


title: How to grab user's mouse click (coordinates) in Javascript

text: I am just getting to grips with using ArcGIS Server 9.3.1. I am developing in .NET with a very basic Web app in C# - it has the standard ESRI web ADF ASP controls in the web page (Default.aspx) including a map and toolbar with some out of the box nav buttons. I have also made some custom tools on the toolbar. _I am interested in working with the Web ADF Javascript libraries not the ArcGIS JavaScript API - a subtle but important differnce I believe._ How would I go about doing the following (bearing in mind I want to be firing the JavaScript via the ClientAction property - custom):   1. I need some javascript to work out where on the Map control the user has clicked and get a hold of the xy coordinates.   2. Send these coordinates to a web service (I already have a web service set up to accept coordinates and send XML back).   3. Parse the XML from the web service and put in a table. Cheers


title: Script to either rotate existing text element or insert graphic element

text: I have about 300 maps that contain a text element that needs to be rotated. I am pretty new to python but I have heard how great it is. This seemed like the perfect task for python but I have had no luck finding any examples on how to do this. The text is "Figure 1-1" and it just needs to be rotated 90 degrees and moved to a set location in the corner. If rotating is not an option I thought I could just manually rotate the text and save it as a graphic and then write a script that would insert the graphic in the desired location. This also doesnt seem possible. Any help would be much appreciated.


title: Openlayers Vector Feature Popup

text: im trying to create a popup in Openlayers that would open when someone clicks a marker. And the marker is setup like this: `new OpenLayers.Feature.Vector()` and not with the ordinary marker layer. The marker is then placed inside another layer `vesselLayer.addFeatures([markerLayer])` which is then added to the map: `map.addLayer([vesselLayer])` Now, How can i get a popup to show when someone selects one of the markers? Let it be noted that each vesselLayer has only 1 marker. But there can be a 100 dynamically added vesselLayers to the map. (Made it that way for easy access trough the LayerSwitcher). I tried reading the docs: http://openlayers.org/dev/examples/select-feature- openpopup.html But i could not figgure it out. The error that comes out of the console is saying that it cannot read the property "events" of null.. I double checked the object instantiation and it does get instantiated and the property is a part of the name. Heres the object logged out right after instantiation. ![Object Overview](http://i.stack.imgur.com/Z2aGr.png)


title: Get Spatial Table Information Oracle Spatial

text: How can I get the tablename, Srid and GeometryType of all spatial tables in Oracle(I use 11g)? In postgis, I can do something like               select f_table_name,srid,type from public.geometry_columns order by f_table_schema, f_table_name;      Does anything similar exist in Oracle?


title: How can you parameterize ColorMapEntry in a style?

text: Is is possible to parameterize ColorMapEntry in a style. I want to be able to specify different values for "quantity" in ColorMapEntry when I request a raster layer. Thanks!


title: Raster operation with grass.mapcalc

text: I need to iteratively add each raster map of viewshed analysis obtained from a series of single observer locations. For this i used grass.mapcalc with Python (GRASS 7 for Windows). However, my problem is that in `grass.mapcalc` you cannot allow functions like the following that we do apply in normal programming.               grass.mapcalc("viewshedN = viewshed + ViewshedN", overwrite = True, quiet = True)      Therefore, i have a trouble in adding each raster iteratively. Is it possible to add an incremental number to the output each time a new raster created? I just tried this in a loop with an integer variable but failed (parse error) in the following code. Is there any other work around/ method or alternative to add each raster result iteratively in GRASS with Python?               formula = 'viewshed' + str(i)               grass.mapcalc(formula, overwrite = True, quiet = True)


title: deploying geoserver on linux or windows: performance implications?

text: I am deploying my first geoserver in production. I develop on a Macbook Pro retina, and my client uses Windows, so I deploy on Windows servers. My first impression was that geoserver seemed a lot slower on the Windows server than my development machine, but I was able to improve this greatly using all the performance tips found all over the net. In short: install Native JAI, tune jvm, enable gwc. But still. In my search for performance tips, I stumbled accross the following results: http://www.slideshare.net/gatewaygeomatics.com/wms-performance-shootout-2011 In this slideshow is shown that MapServer on linux is considerably faster than Mapserver on windows. I tried to look for a similar results or benchmark for Geoserver, but did not find it. Aside of that I also wonder if the container has any impact on the performance of GeoServer (e.g. jetty/tomcat/jboss/...).


title: How to colour polygons with a completely transparent fill using MapServer?

text: Is it possible to fill polygons so that they are completely transparent in MapServer? I want to have a black background that is slightly opaque, and wherever there is a polygon I want that to be completely clear with no fill. I set my background in the following way:                  IMAGECOLOR     0 0 0        TRANSPARENT OFF      I've tried creating a class that is white to fill the polygon, and experimented with opacity using `TRANSPARENCY` with values from 0 - 100 but it does not result in a clear polygon:               CLASS           NAME "light"             EXPRESSION ([max_cost] == 2000)                STYLE                    OUTLINECOLOR 255 255 255             COLOR 255 255 255             SIZE       1             WIDTH      2.0                   END     END       This is results in the following type of image: ![enter image description here](http://i.stack.imgur.com/vWcYD.png) Instead of the polygon having a white fill; I want it to have a fill with no color. I found one suggestion online to define the color like this:               COLOR -1 -1 -1      but it also did not work. In this part of the documentation, the following is described: > “PALETTE=/path/to/palette.txt” is used to define the absolute path where > palette colors can be found. This file must contain 256 entries of r,g,b > triplets for RGB imagemodes, or r,g,b,a quadruplets for RGBA imagemodes. The > expected format is one triplet (or quadruplet) per line, each value > separated by commas, and each triplet/quadruplet on a single line. If you > want to use transparency with a palette, it is important to have these two > colors in the palette file: 0,0,0,0 and 255,255,255,255. After reading this, I also tried using `COLOR 0 0 0 0` but this throws an error: > Unknown identifier. Parsing error near (0):(line 71) I don't understand how I refer to a define color in `PALETTE=/path/to/palette.txt` but this seems like one possible way to achieve this.


title: Shortest path to draw all polyline routes in geometric network without repeating

text: My problem is related to the "shortest path" problem, however it's slightly different in that I don't want to draw one path but to draw all paths without repeating any routes (as much as is possible). I'm trying to work out what this geometric problem is called (if it has been covered at all) in order to find an algorithm to work with. Say there are 6 points in a network that are attached via a web of 10 edges (non-weighted). I must make a route starting form any one single point and move down every edge in any order but without repeating any of these edges. When I find myself in a position where there are no possible new edges forward I can either take one of the previously travelled edges back or I can go "off- road" directly to a new point, however these are considered "lost time". What I want is to find an algorithm (or even the name of this geometric problem) that will minimise "lost time" As I attempt to move down all routes. An easier way to explain might be a printer that draws out all of these paths: while the printers pen is down (i.e. drawing) we are not wasting energy, however as the printers pen will not repeat any drawn lines it may have to move between points with the pen up, which we consider as wasted energy. I want to minimise energy waste. Can anyone help point me in the right direction?


title: What is the unit the shapely length attribute?

text: I'm doing a very simple calculation of the length of a polyline using shapely:               from shapely.geometry import LineString     ...      xy_list = [map(float,e) for e in xy_intm]     line = LineString(xy_list)     s = '%s,%s,%s' % (fr,to,line.length)      My coordinates are in WGS84. I can't seem to find any information about shapely's length attribute. What is the unit of the length attribute? Is there an easy way to convert to to km or meters?


title: How to find matching pixel values in seperate rasters

text: I am looking to scan two separate rasters, capture the pixel information, then draw a polygon of the matching pixels. any thoughts? I am at a loss. Edit: I have two aerials, both are in different geographic areas (same projection and within about 1 km of eachother). I am looking to determine, if between the images, there are matching pixel grids.


title: Match bicycle routes to street lines

text: I have a KML layer of bicycle routes. I have street line layers. Using either ArcGIS or QGIS, how can I select the streets that correspond to the bicycle routes?