title: OpenLayers: failed to create a Ext.list that wraps a Openlayers.request.Get function

text: I am using Sencha Touch + OpenLayers to develop a demo app. One feature is to get the catalog of feature data in GeoServer. I want to use Ext.List to display the name of feature layers. I tried following codes:               App.AddFeatureLayerList = Ext.extend(Ext.List, {         createStore: function(){         Ext.regModel('WFSLayers', {             fields: ['name', 'title', 'srs', 'featureNS']         });         var data = [];         var request = OpenLayers.Request.GET({             url: root+'/geoserver/ows?service=wfs&version=1.0.0&request=GetCapabilities',             handler: function(response){                 var XMLformat = new OpenLayers.Format.XML();                 var xml = XMLformat.read(response.responseText);                 var CAPformat = new OpenLayers.Format.WFSCapabilities();                 cap = CAPformat.read(xml);                      for (var i = 0; i < cap.featureTypeList.featureTypes.length; i++) {                     var featureT = cap.featureTypeList.featureTypes[i];                     data.push({                         name: featureT.name,                         title: featureT.title,                         srs: featureT.srs,                         featureNS: featureT.featureNS                     });                 }                 return new Ext.data.Store({                     model: 'WFSLayers',                     sorters: 'title',                     data: data,                     getGroupString: function(record){                         return record.get('title')[0];                     },                 });             }         });     },     initComponent: function(){         this.store = this.createStore();         this.itemTpl = new Ext.XTemplate('{name}');         this.grouped = true;         this.listeners = {             itemtap: function(dataview, index, item, e){                       }         };         App.AddFeatureLayerList.superclass.initComponent.call(this);     }     });     Ext.reg('app_addFeatureLayerList', App.AddFeatureLayerList);      However, i got an error:               Uncaught DataView requires tpl, store and itemSelector configurations to be defined.      then if I remove the codes:               new Ext.data.Store({                 model: 'WFSLayers',                 sorters: 'title',                 data: data,                 getGroupString: function(record){                     return record.get('title')[0];                 },             });      and put them to pass to `this.store` after call "`this.createStore()`", I will get a empty list. however, I traced array `data` and it turns out it has data, but the store is empty. another funny thing is if I put '`createStore()`' as a public function, and call it in list configs, e.g, '`store: createStore()`', I will get an error message saying my IP is not allowed to access to target URL, which this is not correct since I use a proxy and I added other WFS layers from this server already... I can't find what's wrong with it. Anybody could help me give me some hints? Appreciate all help!!!


title: Map Design: Pros and Cons of using ArcGIS Halos to improve readability of text in a map

text: I have to create a map which involves 2 overlapping Layers. The underlying layer is coloured, the upper layer is hatched. On top of this I need to annotate the regions. In ArcGIS the dirtiest method is to simply attach a halo to the text. Then the text is always readable. But I find the halo covers up so much area of the map and in print looks pretty awful. I was thinking of exporting the geometries into The Gimp or Illustrator....but then my GIS Workflows become very complicated and updates are much more difficult to do. Does anyone have any good examples of annotating under similar conditions and can offer some good advice? ![Map with hatching and annotations](http://i.stack.imgur.com/z9uEa.png)


title: How to integrate SAGA Python Modules into ArcGIS Desktop?

text: I have been having a look at SAGA (http://www.saga-gis.org/en/index.html) and it seems to have some nice tools for geoprocessing. I was wondering if anyone out there has had any luck bringing SAGA Modules (python) into ArcGIS and ArcToolbox. If so what steps did you use. Thanks in Anticipation, CDB


title: ArcGIS: Maintaining the same layer styles over multiple data frames

text: Let's say you have a map with three or four data frames. In each data frame there is the same road layer. Is there a way to ensure that all the symbology/labeling/annotation is the same for that layer across all data frames. So let's say I modify the symbology of the road layer in data frame 1, is it possible to automate that change across all data frames?


title: Statistical relationships between rasters

text: I am trying to analyze how sets of rasters relate to each other using some statistical techniques. As, I don't have much experience using the spatial statistics tools in _ArcGIS_ I was exporting my rasters as _Ascii_ files, and analyzing them using _R_ (specifically the `maptools` package, and `readAsciiGrid()` ). This has been functioning ok (but as the datasets have 90,000 points it is slow to run the analysis), but I don't know if I am recreating in _R_ , existing functionality in _ArcGIS_. For example, I want to perform regressions between each of these rasters using a few different transformations (logarithmic, exponential, etc). Can this be done within _ArcGIS_? A second broader question is if there are standard statistical methods for examining this type of data? Each raster pair has matching data/no-data values and all parameters are identical, aside from the gridcell value.


title: Question on how to use the editor tool and snapping in ArcGIS 10

text: I am new to ArcGIS 10 and am having trouble in using the new editing tool and snapping interface they have now. I cannot find the Task and Target options that allow me to choose Modify feature and define which layer I need to edit (and toggle between layers). Also, I could not find the Trace tool that used to be in ArcGIS 9.x versions that could be useful for redrawing curvy boundaries. There seems to be a problem with Sanapping because now there is a small toolbar that gives me only the options (End, Vertex, Point, Edge, Tangent, etc) but without reference as to which layer I want to snap to. Yet I can not snap properly to the vertex I want. Previously the cursor would be snapped to the point I was hovering over within a reasonable proximity (defined by the tolerance). Now it seems that that sometimes does not happen even if I zoom in close enough. This of course results in problems since sometimes I got the snapping right, sometimes not. Has anybody experienced similar problems and if so how did you resolve them? I would appreciate the information very much. Thanks


title: "File not found at..." error when loading KML into Google Maps

text: I have a KML file on a public URL like http://host.com/dir/file.kml If I put that KML file on another server, and put the new URL in the search box on Google Maps, it loads the KML into google maps. If I put the above URL into google maps like so: http://maps.google.com/maps?q=http://host.com/dir/file.kml It does not load, and I get the following message:               File not found at http://host.com/dir/f....          Suggestions:              Make sure the URL is spelled correctly.         Make sure the file exists.      Therefore, I assume this is a configuration issue on the server. They are not being very helpful with debugging/configuring. Can anyone here offer advice on how to fix this?


title: How to get rid of double labels on ArcGIS 10

text: I'm displaying several point shapefiles. I want them all labeled with the field "ident" which is a 3 digit number starting at 001. Two of the three shapefiles show the labels correctly but the 3rd one shows a double label for each point. Instead of just saying "027" it says "027027". In the Layer properties window, I tried unchecking the box that says "Label features in this layer" but this removes the labels completely, then when I recheck "Label features" in the menu that comes when I right click the layer in the table of contents, the double layers re-appear. Any suggestions would be much appreciated! Molly


title: Looking for HowTos or Examples of Presentation Platforms targeting Management/Clients

text: Any good ideas about using Google Earth to present self contained geospatial data&app to a client or senior management? The broader issue being explored here is that many of our clients are not likely to maintain GIS systems themselves, or be competent enough to install a free app and then browse the data. And we would like the legacy we leave to a client to be independent of us, so we don't want them to be relying on our server or any one else's server. A single file container (yes I know you will say SpatiaLite but what app?) that can be USB'd or the like. One feature I really like in GE, is to geolocate photos and have them appear on the map, to provide a record, or some interactivity etc. I can do the same in QGIS, but not without knowing a lot about the installation. In GE, I know that if they have got the KMZ, they will find the photos. And what about senior management? Because you know they are luddites at best, won't have time to view it when you have time etc. So my thought here is to put GE up as the app, and KMZ as the container, and ask the question in this forum. What is the best example of GE (or your choice app) as a self contained presentation?


title: How to calculate the number of buildings along the left/right side of a street?

text: I'm struggling to find a way to get the line direction stored on shape file from its geometry, is there any tool can get it for me? Also, I have two layers of map: Roads layer (of type polyline), Building layer (of type polygon) I want to get the number of buildings located on the left, right sides of a street segment. Is there any way to do that?


title: Table not found ConvertFeatureClass how to project?

text: I get a Table [Streets] not found when I try to convert a shapefile to a geodatabase. I'm thinking I'll try setting the spatial reference like in this post, http://forums.esri.com/Thread.asp?c=159&f=1707&t=299127 But I'm not sure what the steps are to get a GeometryDefEdit from my converter code               //input      IWorkspaceName pInShpWorkspaceName = new WorkspaceNameClass() as IWorkspaceName;     pInShpWorkspaceName.PathName = shapeFileDirectory;     pInShpWorkspaceName.WorkspaceFactoryProgID = "esriDataSourcesFile.ShapefileWorkspaceFactory";     IFeatureClassName pFCName = new FeatureClassNameClass() as IFeatureClassName;     IDatasetName pShpDatasetName = pFCName as IDatasetName;     pShpDatasetName.Name = shapeFileName;//this is your input file      pShpDatasetName.WorkspaceName = pInShpWorkspaceName;     IName sourceName = (IName)pFCName;     IFeatureClass sourceFeatureClass = (IFeatureClass)sourceName.Open();     IFields sourceFields = sourceFeatureClass.Fields;               //output     IWorkspaceName pWorkspaceName = new WorkspaceNameClass() as IWorkspaceName;     pWorkspaceName.WorkspaceFactoryProgID = "esriDataSourcesGDB.FileGDBWorkspaceFactory";     pWorkspaceName.PathName = GDBPath; //output file      IFeatureClassName pOutputFC = new FeatureClassNameClass();     IDatasetName pDataSetName = pOutputFC as IDatasetName;     pDataSetName.WorkspaceName = pWorkspaceName;     pDataSetName.Name = "Streets";     IFeatureDataConverter pShpToFc = new FeatureDataConverterClass();          pShpToFc.ConvertFeatureClass(pFCName, null, null, pOutputFC, null, sourceFields, "", 1000, 0);      Also is there anything else I should try ? I guess it was the 5th Param into ConvertFeatureClass, but unfortunately it didn't make any difference. _**_** _*_ **** _Update_ ** _*_ ** _*_ *** So I did the following operation:               ConvertFieldName(sourceFields, "CITY", "CITY3"); //change column name from city to city3          void ConvertFieldName(IFields fields, string oldname, string newname)     {       if (string.IsNullOrEmpty(oldname)) return;       if (newname == oldname) return;       List<string> fieldnames = GetFieldNames(fields);            IFieldsEdit editfields = fields as IFieldsEdit;       int index = editfields.FindField(oldname);       if (index != -1)       {         IField field = editfields.get_Field(index);         IFieldEdit addresspd = field as IFieldEdit;         addresspd.Name_2 = newname;         addresspd.AliasName_2 = newname;       }     }      And also I try this: (Change column name, and change it back)               ConvertFieldName(sourceFields, "CITY", "CITY3");      ConvertFieldName(sourceFields, "CITY3", "CITY");       Seems insane, I know, but for some reason, doing these operations (either changing or changing back and forth) allowed the GDB file to be created. I guess the code was cleaning up a resource or something, which allowed for the GDB creation.


title: Geoserver replies HTTP Error 503

text: I´ve seen this question was asked before but I couldn´t find any solution to the problem. Since today, my Geoserver does not reply on any requests anymore, I can´t even access the Geoserver UI. I didn´t change any settings since the last time I used Geoserver and it worked fine then. This is the HTTP error message:               HTTP ERROR: 503     SERVICE_UNAVAILABLE     RequestURI=/geoserver/web     Powered by jetty://      Going further, access to my PostGreSQL/PostGIS Database is denied as well, telling me...               could not connect to server: Connection refused (0x0000274D/10061)      Is the server running on host "127.0.0.1" and accepting TCP/IP connections on port 54321?       Has anybody experienced something similar? If lines from the log are needed, let me know which ones. Didn´t want to post all the content in here. Any help is appreciated! Philippe EDIT: This is probably why it fails:                   ERROR [context.ContextLoader] - Context initialization failed     org.springframework.beans.factory.BeanCreationException: Error creating bean with name 'gwcCatalogConfiguration' defined in URL [jar:file:/C:/Program%20Files%20(x86)/OpenGeo/OpenGeo%20Suite/webapps/geoserver/WEB-INF/lib/gwc-2.1-SNAPSHOT.jar!/geowebcache-geoserver-context.xml]: Cannot resolve reference to bean 'rawCatalog' while setting constructor argument; nested exception is org.springframework.beans.factory.BeanCreationException:      Error creating bean with name 'rawCatalog' defined in URL [jar:file:/C:/Program%20Files%20(x86)/OpenGeo/OpenGeo%20Suite/webapps/geoserver/WEB-INF/lib/main-2.1-SNAPSHOT.jar!/applicationContext.xml]: Initialization of bean failed; nested exception is java.lang.RuntimeException: java.lang.IllegalArgumentException: Resource must be part of a store


title: Does mapbasic support function overloading?

text: I am creating application with mapbasic. Now that it's been some time with mapbasic I found that the language is very fundamental and does not support functionality that I took for granted in Java or Other language. > `How do I know what mapbasic supports or does not support`? I had a look at reference guide as well as user guide but could not find anything.


title: How to edit the PGDB data in QGIS

text: I have a small query regarding the Q-GIS software. I hope someone in the know would be able to answer them.   1. I am using 1.9.0-Master q gis. My datasets are inside a ESRI personal Geodatabase. I am able to load and display the vector data using the “load vector” but there doesn’t appear to be any way of loading the rasters inside the PGDB.   2. Along the same line as #1, I can only display the vector data inside the PGDB but am unable to edit it. Is it possible to edit them at all? If not, are there any other ways of storing multiple datasets in a single collection and edit them?


title: Arcmap Update field from a different table

text: I want to update a field in one table using two matching values in another table. Let's call the fields A, B and C in one table, and X, Y and Z in the other table. So, where A=X and B=Y then update C with Z


title: Is there a version of the ESRI ArcGIS Symbol Library PDF for ArcGIS 10, or some other cheat sheet for symbols/line weights?

text: ESRI has an old set of PDFs from 2004 showing all the symbols available in ArcGIS at the time and their names. I don't need it often but if I'm searching for a certain symbol and don't know its name it's sometimes easier to look at a few sheets of paper than to use the Symbol Selector. It doesn't seem like it's been updated for ArcGIS 10, though. I was wondering if anyone is aware of a new version of this, or if anyone has created one on their own they might be willing to share.


title: Redcap Shapefile Error

text: For some reason I am unable to import shape files in the form of point data to REDCAP. Do these data points need projected coordinate systems? I am aware that the sample size is 500 when double clicking the REDCAP file. I am also encountering errors when using the command line to launch the application. Java is not recognized as an internal or external command when trying to execute the code included in the manual to launch the application. There are apparently memory unless launched from Java. This enables working with larger data sets. Neither polygons or points are opening, the data have been tested in a GCS and PCS. The sample data opens without flaw. http://www.spatialdatamining.org/software/redcap Any suggestions would be much appreciated.


title: Smooth, not simplify, lines features in PostGIS

text: I'm trying to automatically generate paths to place labels on that follow the shape of rivers, 4wd tracks and other very windy lines. At first I thought the `st_simplify()` function was the answer, but it basically has exactly the opposite goal. Here's a result of `st_simplify(way,300)`: ![enter image description here](http://i.stack.imgur.com/xP0Fq.png) `st_simplify()` is trying to preserve the maximum amount of information from the underlying geometry, resulting in a very spiky line. What I'm looking for is a rough approximation of the line - something that would kind of pass down the center, without actually hitting any of the corner points. The ultimate goal is to get more labels rendered in TileMill, without having to increase the `text-max-char-angle-delta` setting, which looks pretty bad. Is there a way to achieve this? I'm looking for a script or function to call - not just an algorithm to attempt to implement. (OpenStreetMap data in PostGIS, rendered using TileMill. Happy to pre-process - performance is not a big concern.) Also, note the difference from this question - I don't want "a smoother line of roughly the same significance to a human interpreter", I want "a smoother line roughly in the centre of the old line".


title: I made a Relationship Class(1-M). Its a street X our Address table. I can´t split the street. How to fix it?

text: I have a table "Addresses" and the Feature "streets". The feature "Streets" are lines. I need to do a relation with this table with this feature. What I do? I tried built a "relationship class". The type is "1-M". It is the line of a street with the table of a adress. But happens an error. I can´t split the line of a street. Anyone knows how to fix it? That is the error: "The feature could not be split. The field is not nullable." I found the error "FDO_E_FIELD_NOT_NULLABLE" description on a ESRI site. But they dont told how to fix it. This is the link: http://edndoc.esri.com/arcobjects/8.3/ComponentHelp/esriCore/fdoError.htm I need the relationship and to split the streets! How to solve/fix this error?


title: Hindi character support in MS SQL 2012

text: I have spatial MS SQL 2012 and i need to store hindi characters in a column. I tried using nvarchar datatype but I'm getting only "`Question marks(????)`". Is it possible to store non-english characters?


title: QGIS 2.4 : How to search an Attribute Table?

text: I am learning how to use QGIS on v.2.4. When opening an attibute table, it does not have any clear "look for" search bar as seen in older versions in screenshots and videos. There is now a (SQL?) query builder that takes much longer to operate than a simple search bar. **How to efficiently search within an attribute table with QGIS 2.4+ ?**


title: svg export not keeping line widths

text: I just upgraded to 2.4 I have noticed that when I export as an svg, for opening in illustrator, the line/pen/border widths are coming out all different. For example a line width of 0.25 in QGIS is coming out larger than 1 in illustrator. Before the difference was only about 0.01 so not really noticable. My border widths are in millimeters. Has anyone else had this problem or know how to remedy it?


title: Full featured opensource webgis solution with editing capabilities based on user roles

text: we are looking for a complete opensource webgis solution with editing capabilities. Permission to read/write on layers should be based on user roles, or at least user/group acl. The system should be able to upload multimedia elements too. Is opengeo suite sufficient? other alternatives?


title: Layer with custom UI not working properly

text: I have a problem with QGIS 2.4 when introducing layers with a custom UI. More precisely, widgets such as CheckBoxes are not recognised as such (appear as Text Edit in Attributes) and the corresponding Checked value appears as NULL (although checked by default in the UI file). When I am manually changing the specific field:               instudy_id = layer.fieldNameIndex("instudy")     layer.setEditorWidgetV2(instudy_id, 'CheckBox')       Everything is normal. The corresponding part of the UI file is: ![UI file](http://i.stack.imgur.com/WGqSi.png) Sorry if my explanation is poor, I am quite new to QGIS. Hope the screenshots help. ![Layer -> Properties](http://i.stack.imgur.com/hUl0X.png) ![Attribute Table of Layer](http://i.stack.imgur.com/xgLEV.png) To conclude:   * some options of my Widgets (defined in my .ui file) are overwritten by QGIS (e.g. "Checked" or "Unchecked" state of a QCheckBox)   * others, which are introduced in the corresponding .py file (like QgsEditorWidgetWrapper.fromWidget( in study ).setEnabled(False) ) are not respected by QGIS.  What is the proper way to introduce a custom / customise an .ui file ? Is there anywhere detailed information/examples on the use of newer methods like QgsVectorLayer.setEditorWidgetV2, setEditorWidgetV2Config, QgsEditorWidgetWrapper ?


title: How to create a Shapely LineString from two Points

text: If have two points, from which I want to create a straight `LineString` object:               from shapely.geometry import Point, LineString     A = Point(0,0)     B = Point(1,1)      The Shapely manual for `LineString` states: > A sequence of `Point` instances is not a valid constructor parameter. A > `LineString` is described by points, but is not composed of Point instances. So if I have two points `A` and `B`, is there a shorter/better/easier way of creating a line `AB` than my current "best" guess...               AB = LineString(tuple(A.coords) + tuple(B.coords))      ... which looks rather complicated. Is there an easier way? **Update** With today's released Shapely 1.3.2, the above statement from the manual is no longer correct. So from now on,               AB = LineString([A, B])      _works_!


title: TypeError: L.Proj is undefined in leaflet

text: I am using leaflet 0.7.3. I am trying to load a wms tile:               var crs = new L.Proj.CRS('EPSG:2236', '+proj=tmerc +lat_0=24.33333333333333 +lon_0=-81 +k=0.999941177 +x_0=200000.0001016002 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs',     {         resolutions: [8192, 4096, 2048],     });          // pass crs into map constructor     var map = new L.map('map',     {         crs: crs     });          ...      But am getting an error: > TypeError: L.Proj is undefined Has anybody received this error, and if so, how did you resolve this issue?


title: What date was a table loaded into PostGIS

text: Is it possible to find out the date a table was loaded into a PostGIS schema using PgAdmin3? The tables are usually supplied from varying external sources.


title: QGIS API documentation in French?

text: I am looking for a document which explains what the API of QGIS is, with examples of this API but in French.


title: Population by Census Tract

text: I am looking to obtain count of population by census tract. The ideal output would be Census Tract 1919.01 Population = 7,435. Does anyone know of resources where I can get this information. I tried searching Census.Gov w/ no luck.


title: Adding shapefile elevation data to dted

text: I am a newbie to QGIS. Currently trying to augment a dted raster value file with vector shapefile(polygon) containing height of buildings. Have converted the shapefile into raster via rasterize function. WheneverI mouse over the converted shapefile(now raster), I am able to get the value of the polygonal shape (from the band value) corresponding to the old shapefile polygon elevation For the dted, the mouseover seem to give me the quantized color depth(band value). How could we get the real height value? Or display them? Lastly, how do we add the two values together and merge them into a single raster data set. I have already projected them to use the same reference coordinate system Thanks


title: How to import and project ASTER 1B data into Arcmap?

text: I recently purchased a granule of ASTER 1B data from ERSDAC, when I import any of the associated rasters into ArcGIS, the image is displayed as being only 1m across, when the actual area represented is tens of kilometers across. Aside from this, the image is also in the incorrect location. The raster came already projected in UTM zone 35N but the area the raster represents is in UTM Zone 35S. I have tried both defining and projecting the coordinate system of the the raster to WGS1984 Zone 35S to no avail. I am using ArcGIS 10.2 Does anyone have any idea of what the issue might be?


title: Saving Shapefile to File Geodatabase in QGIS 2.2?

text: I searched the StackExchange forum for an answer too my question and apologize if I missed it. I'm using QGIS 2.2 64-bit and able to read feature classes from an ESRI file geodatabase without a problem. The problem occurs when I attempt to save a layer to a file geodatabase by right-clicking on the layer and selecting "Save As." I have "ESRI FileGDB" selected as the desired format and enter a name for the exported feature class. The result is a file geodatabase created with the name of the feature class and an error message > Export to vector file failed. > > Error: creation of layer failed (OGR error:Error: Failed at creating > FeatureDataset (General function failure.)) Any idea why the feature class export to a file geodatabase fails?


title: Creating counties in QGIS

text: I have a raster map (a photocopy I took from an Atlas), and a vector map of the UK. The raster map shows pre-1890 Scottish county borders. I have geocoded the map so that the raster borders overlap the vector map, but here comes my question: how do I create a vector layer showing the separate counties? This may be especially difficult, since some of the counties are not contiguous.


title: GoogleMaps URL string to display unique marker

text: This is a topic I thought as easy as pie to solve with just a Google search and which is in fact quite complicated. When I use the following URL string to open a GoogleMaps window on a specific location, I get 2 markers (one green on the exact position of the Lat/Lng coordinates I put in the URL, and one red which is the nearest address). I would like to get just the green marker, which means opening a new GoogleMaps window with a marker matching the Lat/Lng coordinates from the URL. This is quite simple, but I can't find out on the web if it possible or not. This is possible with OpenStreetMap using the following URL.


title: Reprojected Shapefile does not align with Openlayers Plugin

text: I am based in NSW, Australia. I have a shapefile that when I underlay with the openlayer plugin (google satellite), the shapefile is skewed despite all my efforts to correct it. I changed the project properties to GDA 94 and GDA 94 MGA 56, as well as re-saving the shapefile as WGS 84 EPSG:3857 to match the openlayer plugin. Nothing seems to work. Can anyone help? Is it the original shapefile?


title: Updating an ArcPad Server Extension publication

text: If you have succesfully created a publication (axf file) that syncs between the server and multiple users, what happens if somebody wants to change the publication data, e.g. add some new fields ? Can you just republish the data and edit your original forms, or do you need to format all of the re published default-generated form controls and redevelop any event scripts again for the new publication ?


title: How can i convert .tif file into headerless (flat) binary files?

text: I am using NDVI files (value range -1 to +1) created from Landsat CDR images to produce a timeseries in TIMESAT software. My .tif file has following characteristics. File type TIFF,Data Type Double Precision, Byte Order Host (Intel), Interleave BSQ As per the requirement of TIMESAT, the image should be headerless(flat) binary file, and values can be stored as 8 bit integer(0-255), 16 bit signed integer(-32767 to +32768) or 32 bit real. Now i need to change my data type into any of these 3 format and .tif file into headerless binary files. Any idea?


title: How to create a summary statistics table with as many rows as iterated features? I only get one row

text: I have a feature layer of 25 buffers, which are the light blue circles you see. Let’s call each buffer a community. Many buffers overlap. Each buffer has a unique ID. I also have a feature layer of points (centroids). Some of the centroids fall within certain buffers, some fall outside. Every centroid has columns of information associated with it, two of which are important to me: Origins and Destinations. ![Overlapping Buffers and Centroids](http://i.stack.imgur.com/ECd6H.jpg) WHAT I WOULD LIKE TO DO: (1) I want to iterate through each feature in the buffer layer (loop through the buffers one-by-one) and ask: which centroids are completely within this buffer? (2) Then, I want to ask, of the centroids that fall within the buffer, what is the sum of the origins and the sum of the destinations? (3) Last, I want to combine the results of (1) and (2) in a table with three columns of information: The buffer ID, the sum of origins in that buffer, and the sum of destinations in that buffer. MY SETUP and MY PROBLEM SETUP   * I am using the Feature Selection Iterator in ModelBuilder to loop through the buffer file. I have chosen to group by the BufferID field.   * I use the Select Layer by Location tool to select the centroids that are ‘completely within’ each iterated feature in the buffer file. (The buffer is the selecting feature, the centroids are the input feature layer)   * The output of this geoprocess, labeled Centroids(2) in the image below, is, I believe, the collection of centroids ‘completely within’ a particular buffer. ![enter image description here](http://i.stack.imgur.com/4ayHS.jpg)   * I use the Summary Statistics tool to calculate the sum of the Origins and Destinations associated with Centroids(2). In the Summary Statistics tool interface, I select Origins and Destinations as statistics fields and I calculate their sum. The output of this tool is a table. I select the location of this output table.  ![enter image description here](http://i.stack.imgur.com/3sX0U.jpg) PROBLEM   * This is where I am stuck. If I run the model above, I get this: ![enter image description here](http://i.stack.imgur.com/ND6CM.jpg)   * This is interesting, because I see the tools turning to red as many times as there are buffers – so I do believe that model is iterating through each feature in the buffer layer, selecting centroids by location, and creating summary statistics, when run.    * Clearly, though, I have not set up the model properly to get the desired output because the resulting table has only one row. There should be as many rows as there are buffers, which in this case is 25. Moreover, I would like the OBJECTID* in the above image to be the BufferID from the buffer layer. When I grouped the Feature Selection Iterator by BufferID I thought this might carry over to the Summary Statistics.    * Although the model is iterating 25 times (as evidenced by the tools turning red 25 times) the output it generates is for one iteration only. What do I have to change in order to get (A) a table with 25 rows, and (B) an ID that corresponds to the BufferID? A very sincere thank you for reading and helping


title: ArcPy delete terrain: cannot be deleted

text: I can't get my ArcPy script to delete a terrain. If it exists, it should be deleted. However it throws an > ERROR 000317: {path} cannot be deleted My script is a lot longer than this, but this is what it boils down to. fds is a feature dataset in a file geodatabase. The feature dataset contains a terrain with 2 associated feature classes:   * a points feature class used as masspoints    * a polygon feature class used as boundary/softclip Here's my code:               import os     import arcpy          fds = r'D:\Dumfries & Borders.gdb\DEM'          arcpy.env.workspace = fds     l = arcpy.ListDatasets()     for d in l:         if arcpy.Exists(os.path.join(fds, d)):             arcpy.Delete_management(os.path.join(fds, d))      Earlier on I tried to have the script delete the whole feature dataset at once. Didn't work either. I then decided to delete the datasets in it one by one. The ListDatasets() lists only one: the terrain, which sparks the error mentioned above. It is possible to delete the terrain or the feature dataset by hand in ArcCatalog, or using the Python window in ArcCatalog. However, even after closing all other programs (to preclude any locks from obstructing) except my Python interpreter (PyCharm) it still throws the same error. So, any ideas? Will post this to the ESRI forums as well. edit follow-up: removing the spaces and & from the gdb's name does not help, as ListDatasets() won't list anything now anymore. Also scriptwise trying to delete the feature dataset now throws an > ERROR 000732: Input Data Element: Dataset {path} does not exist or is not > supported.


title: ArcGIS multiple scripts conflict

text: I am developing separate javascript widgets in different scripts and they work independently. However, when I try and put both scripts together, neither script works at all. So this brings me to my question, is there any way for me to _easily_ integrate multiple Javascript ArcGIS scripts? For example, if I only change the require statements from:               require([             "esri/map",             "esri/layers/FeatureLayer",             "esri/dijit/PopupTemplate",             "esri/request",             "esri/geometry/Point",             "esri/graphic",             "dojo/on",             "dojo/_base/array",             "dojo/domReady!"           ],      to:               require([             "esri/map",             "esri/tasks/RouteTask",             "esri/tasks/RouteParameters",             "esri/layers/FeatureLayer",             "esri/dijit/PopupTemplate",             "esri/request",             "esri/geometry/Point",             "esri/graphic",             "dojo/on",             "dojo/_base/array",             "dojo/domReady!"           ],      breaks all of my code. It gives me an error on the following line:               var popupTemplate = new PopupTemplate({               title: "{title}",               description: "{description}"             });      In this case, when I say break I mean that the functional aspects of my code no longer work. It just shows an empty map, but none of my feature layers display anything.


title: How to get ArcObjects SDK for .NET?

text: Forgive my newbness, but how do I get the .NET ArcGIS SDK? I can't access any of the ArcGIS templates when I create a new project in Visual Studio. I'm using Visual Studio 2010, SP1 and ArcGIS 10.1. I'm VERY green WRT to programming with ArcObjects and I think I'm missing something really simple! Any help would be greatly appreciated.


title: How do I insert a geometry as WKT to SQL Server Spatial with PHP?

text:
I have a web mapping application built on Bootleaf (Leaflet + Bootstrap). Using the leaflet Draw plugin the user can create a polygon on the map. I'd like to insert that polygon into a geometry field to an SQL Server 2012 using PHP. I have a field called **Boundary** with spatial type _Geography_ in a table called _tblProjects_. There are hard coded values here just to get things working. They'll come from a jQuery post once I have things on the server side ironed out.               <?php         include_once("connection.php");         $connectionInfo = array("Database"=>"$dbname", "UID"=>$user, "PWD"=>$pass);         $DBH = sqlsrv_connect($server, $connectionInfo);              if( !$DBH ) {             echo "Connection could not be established.
";             if( ($errors = sqlsrv_errors() ) != null) {                 foreach( $errors as $error ) {                     echo "SQLSTATE: ".$error[ 'SQLSTATE']."
";                     echo "code: ".$error[ 'code']."
";                     echo "message: ".$error[ 'message']."
";                 }             }         }              $sql = "INSERT INTO tblProjects (ProjectID, Boundary) VALUES (1, geography::STGeomFromText( "POLYGON (-93.30388806760311 27.994401411046173, -94.62224744260311 33.37641235124676, -79.70281384885311 31.80289258670676, -93.30388806760311 27.994401411046173)",4326));";              $stmt = sqlsrv_query($DBH, $sql);         if( $stmt === false ) {             die( print_r( sqlsrv_errors(), true));         } else {             echo "SUCCESS!";         }      It doesn't like my geography type, and I'm not sure how to include it; here's the PHP error: ![enter image description here](http://i.stack.imgur.com/0KDvF.png)


title: Getting null values using Point sampling tool in QGIS 2.0.1

text: I would like to obtain the information of a raster layer at some points. Apparently the Point sampling tool does the job. I have used this tool, and all the raster values extracted at specific points are NULL... Both vector and raster layer are in the same CRS. I use QGIS 2.0.i in Windows 8.1. Any clue ? By the way, does anybody knows how the raster value is computed by this tool = Is it interpolated ? or is it the closest raster value to the point coordinate ? )


title: OpenLayers not refresh a layer

text: i have a problem with refreshing layers.   I have WMS and WFS layers from one source. I edit WFS and after do `wms_layer.redraw(true)` but layer redraw only after zoom in or zoom out. Same with WFS. i have `wfs_layer.visibility=false` and do `wfs_layer.visibility=true` `wfs_layer.refresh({force:true})`. And layer refresh only after zoom in or zoom out.   For do this operations i gonna do a something special?


title: How to programmatically create ArcGIS Image Service and add it to map using C# / .NET?

text: How can I programmatically create ArcGIS Image Service and add that service to my MXD file programmatically? I am using C# .Net.


title: How to determine layers affected by a user-drawn polygon

text: # Scenario: I have a set of layers, each covering one country. In my OpenLayers client, I enable users to randomly draw a polygon. A user is then able to perform analysis based on the drawn polygon. Analysis include getting the population covered by the polygon, population density of the area covered by the polygon, etc. This separation of data into separate layers is for performance reasons, as there are millions of grid record data. Putting all this data in one layer is likely to degrade performance considerably. # The Problem: This works perfectly within a single country, which is represented by a single layer in GeoServer. However, I want to enable users to draw a polygon that covers sections in more than one country. This means that to perform the analysis, data must be extracted from all the affected layers. I'm looking for a pure OpenLayers+GeoServer approach of achieving this. The features extracted are then passed to a GeoServer WPS process. # My attempt so far: I've created one layer group that includes all the layers for all the countries. I wish to load as minimal data as possible to the client. Which means I can't use WMSGetFeatureInfo, or the WFS version that retrieves information based on a WMS layer. Which means that the layer group approach may not be the appropriate one. All my processing is done in WPS processes. It's only the result that is sent to the OpenLayers client. Any suggestions will be appreciated very much. My environment: GeoServer 2.1.3, PostGIS 1.5, OpenLayers 2.11.


title: Can I have Points, Lines and Polygon features in one QGIS layer?

text: I have a situation where planning is done on a "base map" consisting of rivers, contours, roads etc. There is a need to create a "overlay(s)" over the "base map" consisting of drawings, that could be point(s), line(s) or polygon(s) on one overlay. Is it possible to create point, line and polygon features on one QGis layer/overlay I have looked at the vector overlay with a memory data provider but it looked like you can only create a instance of the memory data provider for one type of feature only? I would appreciate any help, thank you.


title: How to split features?

text: I would like to split a shapefile in qgis according to how the polygons overlaps. Don't know how to though? See pic for explanation: ![enter image description here](http://i.stack.imgur.com/qRVb5.png)


title: FME- ask about interpretation of log

text: I'm running some process for about 48 hours and I'm not sure about what is going on or maybe better - how should I interpret the log. Please have a look at attached part of the log. Just after successful ? clipping (at least is something I can read from log ?) the FME executes "performing low-level intersection" and repeat that part at least about 200 times in my case. Is it correct ? Should I patiently wait to finish that part of the process or it somehow stuck in that place ? Thanks in advance for help ! * * *               Finished clipping 303150 / 525807 clippees against all clippers     Performing low-level intersection at phase #1... 100% done     Setting up lines for phase# 2 ... 100% done     Setting up lines for phase# 2 ... 100% done     Performing low-level intersection at phase #1... 100% done     Setting up lines for phase# 2 ... 100% done     Setting up lines for phase# 2 ... 100% done     Performing low-level intersection at phase #1... 100% done     ... Last line repeated 2 times ...     Setting up lines for phase# 2 ... 0.48% done     Setting up lines for phase# 2 ... 100% done     Setting up lines for phase# 2 ... 100% done     Performing low-level intersection at phase #1... 100% done     ... Last line repeated 16 times ...     Setting up lines for phase# 2 ... 0.48% done     Setting up lines for phase# 2 ... 100% done     Setting up lines for phase# 2 ... 100% done     Performing low-level intersection at phase #1... 100% done     Setting up lines for phase# 2 ... 0.48% done     Setting up lines for phase# 2 ... 100% done     Setting up lines for phase# 2 ... 100% done      Thanks for remarks. Two data sources: Geotiffs (about 13k of tiles) and water areas features (about 600k) in shapefile. The bounding boxes of geotiff's in the first part of the workflow should clip the WA features. According to info from Dale process all DEM tiles at once raster feature's don't hold any data... so I believe I'm not processing DEM data at all and in the process only them bounding boxes are used. Then it should be simply use of clipper with clippers and clippees both being vector features. This is also a place of the process that I believe (interpreting the log) the process stuck so I'm not describing the further part of the workbench. If any won't find the extra information enough I'm able to share with the workbench which will be more self-explanatory for the rest of the process (can't find any button to attach it here).


title: Opening Annotated shape files in QGIS

text: Ie've been given a number of files including some annotated shape files - Basically kilometre post markers and cable identification markers, to go with a polyline shape file. I was hoping that we could just drag 'n' drop the files over the images and lines in QGIS map window and they would show, but we get an error message (unknown file type) or the annotations just show up as blobs.![image showing annotations and route with shaded relief](http://i.stack.imgur.com/xPSea.jpg) The files are .DBF files, .PRJ files, .SBN files, .SBX file with an autoCAD LT Shape source and AutoCAD LT Compiled Shape, together with an .XML file. My question is how can I overlay the annotations and get them to show text rather than the blob? I can click on the properties in the layer browser, and see the text that is supposed to display, but I can't work out how to get it to show.


title: How to assign CRS to RasterLayer in R

text: I have a layer with `coord. ref. : NA`. How can I assign `coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0` to `coord. ref. : NA`?               class       : RasterLayer      dimensions  : 127, 146, 18542  (nrow, ncol, ncell)     resolution  : 0.01, 0.01  (x, y)     extent      : -122.39, -120.93, 36.72, 37.99  (xmin, xmax, ymin, ymax)     coord. ref. : NA      data source : C:\2012Tif\2012ASC5min_ppt_spas1306_0001_20121222_0005_UTC.asc.tif      names       : X2012ASC5min_ppt_spas1306_0001_20121222_0005_UTC.asc      values      : 0, 0.0782172  (min, max)               class       : RasterLayer      dimensions  : 180, 240, 43200  (nrow, ncol, ncell)     resolution  : 0.01, 0.01  (x, y)     extent      : -123.305, -120.905, 36.705, 38.505  (xmin, xmax, ymin, ymax)     coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +to wgs84=0,0,0      data source : C:\Users\nahmlee\Documents\1998Test.tif      names       : X1998Test      values      : 0, 0.07884484  (min, max)


title: JSONP working with Geoserver and OpenLayers

text: I'm in this scenario (fighting with the same origin police issue): I have a map initialized with OpenLayers and a WFS request trough the Protocol.Script object that is giving me this error:               Uncaught TypeError: Object function (a,b){this.metadata={};b=OpenLayers.Util.extend({},b);null!=this.alwaysInRange&&(b.alwaysInRange=this.alwaysInRange);this.addOptions(b);this.name=a;if(null==this.id&&(this.id=OpenLayers.Util.createUniqueID(this.CLASS_NAME+"_"),this.div=OpenLayers.Util.createDiv(this.id),this.div.style.width="100%",this.div.style.height="100%",this.div.dir="ltr",this.events=new OpenLayers.Events(this,     this.div),this.eventListeners instanceof Object))this.events.on(this.eventListeners)} has no method 'addOptions'       The error seems to be related with OpenLayers lib (maybe the version, actually is the 2.13), but i don't know. What i tried is:               consulados = OpenLayers.Layer.Vector("Consulados", {             strategies: [new OpenLayers.Strategy.BBOX()],             protocol: new OpenLayers.Protocol.Script({                 url: "http://"+ip+':'+port+"/geoserver/wfs",                 callback: 'handleResponse',                 params: {                     service: "WFS",                     version: "1.1.0",                     srsName: "EPSG:4326",                     request: "GetFeature",                     typeName: "topp:tazmania_roads",                     outputFormat: "text/javascript"                 },             })     });      Any help is appreciated.


title: Los Angeles Public Parking Lot Polygons?

text: So, I'm looking for some relevant data to incorporate into my map, and one of them is public parking polygons. Is there such a .shp file with these data?


title: How to write scripts that can call on both Arcpy and PyQGIS?

text: I have spent quite a bit of time looking for a clear answer this topic and have still been unable to really accomplish my goal of accessing PyQGIS in a stand alone script. Parts of answers have been provided across many different questions such as: Standalone applications using QGIS and environment variables pyQgis in stand alone script, vector loading error How to write standalone Python scripts using QGIS? The most thourough Answer that I have found was provided by @gene **here**. What I particularly am interested in, is being able to write scripts that can call on both Arcpy, and PyQGIS. I have set both PATH and PYTHONPATH system variables according to the PyQGIS cookbook documentation as:                C:\OSGeo4W\apps\qgis\bin,%PATH%;C:\Python27\ArcGIS10.2\Lib\site-packages\PyQt4;C:\Program Files (x86)\Microsoft SQL Server\100\Tools\Binn\;C:\Program Files\Microsoft SQL Server\100\Tools\Binn\;C:\Program Files\Microsoft SQL Server\100\DTS\Binn\;C:\Program Files (x86)\Microsoft SQL Server\100\Tools\Binn\VSShell\Common7\IDE\;C:\Program Files (x86)\Microsoft Visual Studio 9.0\Common7\IDE\PrivateAssemblies\;C:\Program Files (x86)\Microsoft SQL Server\100\DTS\Binn\      and               C:\OSGeo4W\apps\qgis\python; C:\OSGeo4W\apps\Python27\Lib\site-packages      Respectively. I am starting with this script:               from qgis.core import *          QgsApplication.setPrefixPath("C:\OSGeo4W\bin\qgis.bat", True)     QgsApplication.initQgis()     # or your solution     # read a shapefile      layer = QgsVectorLayer('H:\DATA\Airfields_AB.shp', 'Airfields', 'ogr')     layer.isValid()     True     # loop through layer      for elem in layer.getFeatures():         geom= elem.geometry()         attr =elem.attributes()         (processing)      I am looking to avoid batch scripts that set the appropriate settings, but rather call them from inside the Python script. The optimal scenario for me would be to use the IDLE in the ArcGIS 10.2 folder to run the preceding code. I am using a Windows 7 machine with ArcGIS 10.2 and QGIS 2.01. Currently the results of this code is:               Traceback (most recent call last):       File "F:\Users\...\PyQGIS_Test_01.py", line 1, in <module>         from qgis.core import *     ImportError: DLL load failed: The specified module could not be found.      Thank you for any thoughts or comments!


title: How to maintain many-to-many relationship between feature class and table with multiple disconnected editors?

text: I have an application coded with VB.NET (ArcGIS Engine) and its functionality is:   1. digitize a parcel and enter information for its owners, so i have a M-N relationship between them.   2. every parcel has it stones so an M-N relationship ... i have two PCs where the application is installed, and I am using them in disconnected mode. The problem is that I want to be able to export daba and import it between the 2 PCs without losing data or the relationship. Does anyone have any suggestions for how to do this?


title: Where to purchase global boundaries data on city level for commercial use?

text: Recently I was asked to deal with an app geo locating improve project. It can't be completed without any boundaries metadata on city level. But so far most of my research led me only to things for non-commercial use.


title: scaling of results: the comparing of osm-files for a few days - a statistical approach

text: i am pretty new to Geographic-information Systems so do not bear with me if i ask newby questions: i am interested in gatherin information on schools Some days before i have heard that it is interesting to have a closer look at the taginfo-site (see below): i have question: what does this mean https://taginfo.openstreetmap.org/tags/amenity=school here we have two requests on the above site: **a. yesterday: 2014-06-23 23:58 UTC**               Type count Objekts all           **all:**  592672 0.02% -          [Node] Node 325997 - 0.36%     [Way] Way 260585 - 0.11%     [Relation] Relation 6090 - 0.23%      **b. some (guess two or three) days before:** we have had the following values:               countl Objekte all: 592 443 0.02%          [Node] Node 325 974 0.36%     [Way] Way 260 384 0.11%     [Relation] Relation 6 085 0.23%      see the **details:** today all: **592672** 0.02% - some days before: countl Objekte all: **592 443** 0.02% you see: that differs: we have a differende of 229 x (entries, records or some what) **question:** **a.** is it true that we can see the diff - that means how the entries of schools grow? **b.** is there a total growth of records (!!!) of schools or does the this mean - some records have been changed!? look forward to hear from you


title: Finding coordinate pair at any given distance along calibrated route?

text: Goal: given a route as a set of ordered coordinates, measure the overall distance of the route and be able to access the coordinate pair at any given distance along the route. Application: I'm building software to analyze driving routes and I would like to be able to get data for any particular section of the route that I want. For example, if I want to know the elevation profile between 50.3 and 58.7, I need the points from the route that are included between those distances. Problem: I have tried this in Google Earth and various online route building websites like MapMyRide with absolutely horrible results. The process has been to take the file produced by the service, which is a set of points, and interact with the route based on the points in the file. When measured from start to finish, the distance is the same as the service claims. BUT, if I request a point at a distance along the route, using the same calculations for distance between points, I get a points that is anywhere from .5 to 2km off as compared to the distance from the service (basically, the km markers along the route). Other info: This is a python project and I'm using Vincenty's ellipsoid to measure distance. This question is stemming from a conversation that started at Calculated distance doesn't match Google Earth.


title: OpenLayers: How to remove custom controls?

text: I am trying to remove custom draw controls from my OpenLayers application which I create like this:                var drawControls = {          polygon: new OpenLayers.Control.DrawFeature(polygonLayerFreehand, OpenLayers.Handler.Polygon),          select: selectControl,          box: new OpenLayers.Control.DrawFeature(boxLayer,               OpenLayers.Handler.RegularPolygon, {                 handlerOptions: {                     sides: 4,                     irregular: true                 }             }         )     };      and which I add to the map like this:               for(var key in drawControls) {                 _map.addControl(drawControls[key]);                 }      I now would like to know how I can remove this custom control. I found examples for removing the standard OL controls like e.g. this               var ovMapControl = _map.getControlsBy("CLASS_NAME", "OpenLayers.Control.OverviewMap")[0];     _map.removeControl(ovMapControl);      or                var ovMapControl = _map.getControlsBy('displayClass','olControlPanZoomBar')[0];      _map.removeControl(ovMapControl);      but what is the 'displayClass' or the 'CLASS_NAME' in the case of costum controls like my 'drawControls'? Regards and thanks again Jan


title: Labels placement in ArcGIS

text: I need all labels to be placed outside of the boundary (polyline features). Most of the labels are placed outside the boundary line, while others (marked with red color) are placed inside the boundary. What is the best way place all labels on the outside of the polyline? ![Label Placement](http://i.stack.imgur.com/LzdW1.jpg)


title: Get box related info in OpenLayers

text: I just started to play a little bit around OpenLayers and I am experiencing some issues while trying to let the user draw a bounding box within a map. For this, I looked at this example and tried to modify the code accordingly without much success: whenever a new box is drawn, I'm getting the following error:               Uncaught TypeError: Object [object Object] has no method 'getBounds'      This is my code (simplified) so far (please, notice it's written in CoffeeScript):               ######################################################################     ############################## OPEN LAYERS ###########################     ######################################################################              mapOL = null         boxLayer = null         heatmapOL = null         extent = null         panel = null         drag = null         draw = null              $('body').ready () ->             fromProjection = new OpenLayers.Projection("EPSG:4326");    # Transform from WGS 1984             toProjection   = new OpenLayers.Projection("EPSG:900913");  # to Spherical Mercator Projection             extent = new OpenLayers.Bounds(2.101425, 41.309534, 2.2621, 41.45585).transform(fromProjection, toProjection);                  mapOL = new OpenLayers.Map('openlayers', {                 controls: [                     new OpenLayers.Control.Attribution(),                     new OpenLayers.Control.Navigation()                 ],                 numZoomLevels: 18,                 displayProjection: fromProjection,                 projection: toProjection             })                  mapnik = new OpenLayers.Layer.OSM.Mapnik('Mapnik', {                 displayOutsideMaxExtent: false,                 wrapDateLine: true             })                  boxLayer = new OpenLayers.Layer.Vector('Box layer');             mapOL.addLayers([mapnik, boxLayer])                  mapOL.addControl(new OpenLayers.Control.MousePosition())             drag = new OpenLayers.Control.DragPan({                 title: 'Drag map'             })             draw = new OpenLayers.Control.DrawFeature(boxLayer, OpenLayers.Handler.Box, {                 title: 'Select a bounding box'             })             panel = new OpenLayers.Control.Panel({                 defaultControl: drag             })             panel.addControls([drag,draw])             mapOL.addControl(panel)                  mapOL.setBaseLayer(mapnik)             mapOL.zoomToExtent(extent)      The map and its controls are correctly displayed (I am omitting CSS and HTML code for simplicity). Any clues on what I am missing?


title: WHERE query when creating layer in MapServer

text: Is there a way to have WHERE clause when you're querying for postgis data on mapscript ? like, `geom FROM country WHERE country = 'usa' USING unique gid`. My current workaround is by creating VIEW. Thank you. :)


title: Possible to set an attribute field to update automatically based on another field?

text: If I have 2 attribute columns and I wish to change certain values, is it possible for the second attribute column to automatically update its values based on an expression? Example:               Age |  Status     10  |  Child     14  |  Child      If I change the `Age` from _14_ to _18_ , I would like the `Status` to automatically change from _Child_ to _Adult_ , so it would look like the following:               Age |  Status     10  |  Child     18  |  Adult      I must only be able to do this within QGIS and without the use of database functionalities such as PostGIS/SpatiaLite etc. I've looked at **Value relations** but that only seems to provide drop-down box for the user to select the relating columns. I've also looked at **Relation Reference** but it won't allow me to select a "Relation" (not even sure what Relation Reference actually does). Please advise? I am using QGIS 2.2.


title: How can I create a form so the user has to be authenticated to access a project?

text: I have a login form in QGIS and when it is loaded, python gives me an error for, apparently, not recognizing the _findchild_ part. This is due to the fact that the form isn't loaded before the script is executed, so there are no "childs". How can I connect the form so that the user has to be authenticated first in order to access the project?               import psycopg2     import hashlib     from PyQt4.QtCore import *     from PyQt4.QtGui import *          nameField = None     myDialog = None     passField = None          def formOpen(dialog,layerid,featureid):     global myDialog     myDialog = dialog          global nameField     nameField = dialog.findChild(QLineEdit,"username")          global passField     passField = dialog.findChild(QLineEdit,"password")          loginButton = dialog.findChild(QPushButton,"btnLogin")     loginButton.clicked.connect(authenticate)          def authenticate():          conn_string = "host='CORFOU' dbname='addisababa' user='addisababa' password=1"     conn = psycopg2.connect(conn_string)          username = nameField.getText()     hashpassword = hashlib.md5(passField.getText())          cursor = conn.cursor()     cursor.execute('SELECT * FROM users WHERE username=%s AND password=%s', (username, hashpassword))          if cursor.rowcount == 1:     myDialog.accept()     else:     msgBox = QMessageBox()     msgBox.setText("Wrong username or password")     msgBox.exec_()


title: How to Create shapefile using Geometry field

text: I am working on ArcObject. I want to create shape file using arcObject. I am very close to it. but need some better alternative as i feel my solution may fall short in time when there are lakhs of records I have Geometry field in sql database. from that I get the string which is Linstring,Multilinestring or Point Eg. "LINESTRING (-98.71702604799998 37.746413141000062, -98.717210399999942 37.738209632000064, -98.717347586999949 37.732091870000033)" from this I manually extract the coordinates like string d = st.Substring(12); // to get rid of LINESTRING label d = d.Substring(0, d.Length - 1); // to get points as x1 y1,x2 y2,x3 y3 arr = d.Split(','); // now to get array [x1 y1, x2 y2, x3 y3] here is how my code goes private IPolyline CreatePolyline(DataSet ds) { ESRI.ArcGIS.Geometry.IPolyline pSegPoly1 = new ESRI.ArcGIS.Geometry.PolylineClass();                       pSegPoly1.SpatialReference = createSpatialReference();             ESRI.ArcGIS.Geometry.IGeometryCollection pGeoColl12 = pSegPoly1 as ESRI.ArcGIS.Geometry.IGeometryCollection;                  foreach (DataRow row in ds.Tables[0].Rows)             {                 string st = row["ShapeGeoWKT"].ToString();                 string[] arr = null;                 if (st.StartsWith("MULTI"))                 {                     string d = st.Substring(18);                     d = d.Substring(0, d.Length - 2);                     arr = d.Split(',');                 }                 else if (st.StartsWith("LINE"))                 {                     string d = st.Substring(12);                     d = d.Substring(0, d.Length - 1);                     arr = d.Split(',');                 }                 ESRI.ArcGIS.Geometry.IPoint[] pntArrayTemp = new ESRI.ArcGIS.Geometry.IPoint[arr.Length];                 int k = 0;                 ESRI.ArcGIS.Geometry.ILine[] linearray = new ESRI.ArcGIS.Geometry.LineClass[pntArrayTemp.Length];                 ESRI.ArcGIS.Geometry.ISegmentCollection[] paths = new ESRI.ArcGIS.Geometry.PathClass[pntArrayTemp.Length];                 foreach (string st1 in arr)                 {                     linearray[k] = new ESRI.ArcGIS.Geometry.LineClass();                     paths[k] = new ESRI.ArcGIS.Geometry.PathClass();                     pntArrayTemp[k] = new ESRI.ArcGIS.Geometry.PointClass();                     string[] temp1 = st1.Trim().Split(' ');                     pntArrayTemp[k].X = Convert.ToDouble(temp1[0].Replace(")", "").Replace("(", "").Trim());                     pntArrayTemp[k].Y = Convert.ToDouble(temp1[1].Replace(")", "").Replace("(", "").Trim());                     object obj1 = Type.Missing;                     paths[k].AddSegment((ISegment)linearray[k], ref obj1, ref obj1);                     k++;                      }                 for (int j = 0; j < paths.Length - 1; j++)                 {                     object obj12 = Type.Missing;                     pGeoColl12.AddGeometry((ESRI.ArcGIS.Geometry.IGeometry)paths[j], ref obj12, ref obj12);                     linearray[j].PutCoords(pntArrayTemp[j], pntArrayTemp[j + 1]);                 }                  }             pGeoColl12.GeometriesChanged();             return pGeoColl12 as IPolyline;         }      from this I get array of the co-ordinates which i use to create shape for line. It is very tedious procedure for many record[lakhs or more]. Is their any simple easy method in arcObject that I can get the coordinates from LINESTRING on fly or I can use the LINESTRING to draw lines ?


title: How to get around "name 'row' is not defined" error when row in question doesn't exist?

text: I am getting a "name 'row' is not defined" error when I run a particular piece of code (see below). This only happens when my input layer has one entry instead of two. So when the attribute query is run to find a second row there is no row to find and hence the error appears. The aim of my overall script mean this will happen from time to time, what I am looking for is a way to accommodate this error and either avoid this section of script when an error is likely to occur or continue with the script even if this error pops up. I have tried to incorporate a 'try' and 'except' statement into the loop but have been unsuccessful. Would this work? if so how? If anyone has an idea on how to get around this problem please let me know.               import arcpy     from arcpy import env          env.workspace = "C:\MyArcGIS\EDRN_Shp"          arcpy.env.overwriteOutput = True          join110 = "C:\\MyArcGIS\\EDRN_Shp\\join110.shp"     joinField = "LCC_DRN_ID"     query1 ='"FID" = 1'             arcpy.MakeFeatureLayer_management(join110, "feature_layer12")     arcpy.SelectLayerByAttribute_management ("feature_layer12", "NEW_SELECTION", query1)          rows = arcpy.SearchCursor("feature_layer12")     for row in rows:         if row.getValue("edrn_in") >= "E":             arcpy.JoinField_management("feature_layer10", joinField, "feature_layer12", joinField, "edrn_in")         if row.getValue("edrn_out") >= "E":             arcpy.JoinField_management("feature_layer10", joinField, "feature_layer12", joinField, "edrn_out")          del row, rows


title: What projection system is used in China which is aligned with WGS84

text: I have some GPS tracks in Beijing. I want to project it and convert the lat long degree values in meters for calculation in euclidian space. I was wondering what is the most suitable coordinate and projection system China uses that closely resembles WGS 84? (should be one that is geocentric, I guess)


title: loading geojson in leaflet

text: I just added the geoJSON vector layer in leaflet for my hybrid mobile app using phonegap, the problem is that when I add the geoJSON file of size larger than 2mb the application crashes. What is the problem does larger geoJSON files can not be loaded this way? The code for loading geoJSON file in my map window is below.               var roadStyle = {          color: "#00ff00",          weight: 3       };          var roads = [geojson]; // geojson data is added in roads variable          L.geoJson(roads, {      style: roadStyle     }).addTo(map);      Thanks in advance.


title: Losing mouse focus on map in safari browser

text: I am facing problem when I use SAFARI browser for GIS WEB Application. I have select tool in my application to select the features on the map. When i use Safari the point focus is not correct it drags to some other location. I have tried the same functionality in other browsers like IE, Chrome and Mozilla Firefox. It works perfectly. Any suggestion what additional needs to be done to make it work in Safari Browser? Also is Safari recommended for GIS Web Application? The snapshot1 is of Safari web browser where the mouse focus is lost and it goes out. ![snapshot1 - Safari](http://i.stack.imgur.com/zULl0.jpg) Whereas the actual result in other browser is in snapshot2. ![snapshot2 - Firefox/IE/Chrome](http://i.stack.imgur.com/TOPoP.jpg)


title: Looking for Melbourne shapefile?

text: I am looking for Melbourne shape file map (boundaries or suburbs or roads). Just for Melbourne. Any one knows where I can download it?


title: How to store route names of a stop belonging to multiple routes

text: I've been reading the ArcGIS help files on route analysis to gather more information on what attributes are stored with `Stops`. Maybe I've stared at the screen too long, but I cannot seem to understand the concept below. ![RouteName description](http://i.stack.imgur.com/7WnlM.png) When I `Load Locations` into the `Stops` layer, the default value in `RouteName` is null, as expected. I do not understand what is meant by the second sentence: > By using this property, stops within one route analysis layer can be > assigned to multiple routes. How is this done, if `RouteName` can only contain one route's name? I tried to understand the theory by using an example. Two buses `BusA` and `BusB` travel in different sections of town, on completely different routes. The only bus stop they both stop at is `Stop4`, so if someone was on BusA and needed to change routes, they could only get onto BusB at Stop4. When the network dataset was set up, the two bus lines were assigned separate connectivity groups, while the bus stops was assigned to both. When using the Route solver for this scenario, how is `Stop4` assigned a route name, since it is common to both bus lines? Will it only be assigned the route name of one of the buses, but not both?


title: Integrating ArcGIS Online Web map to my Java script file(test.html)

text: Is it possible to integrate ArcGIS Online web map to my java script API?? I check the samples but in which only web map is added..I have few layers in .html file (Like feature,dynamic etc. and I want to integrate one web map to my application.) Any help will be highly appreciated : ) thank in advance : )


title: How to visualize geolocated data on a slippy map on a website?

text: I am trying to visualize some data on a website that has been measured while driving around. My source is a CSV file with longitude, latitude and my measured data. I have been able to export it to a GPX file and to display it using OpenLayers using the instructions found on http://wiki.openstreetmap.org/wiki/DE:Openlayers_Track_example. The problem is that I would like to not only visualize the (points of the) track but also the measured data that has been recorded with them. Some variation in color or line width would be great, but I could not find an example ho to do anything like this in OpenLayers (or another JS library). I guess it would be possible to use rule based styling as explained on http://docs.openlayers.org/library/feature_styling.html#rule-based-styling but I was unable to find out how to retrieve the property from the GPX file. I found an example that partially matches my needs on geo.dianacht.de/osm/track.html but this uses different data sources for track and data. It would be more useful if I could just use one source. If this should not be possible, any other suggestion is highly appreciated.


title: What needs to be done to downgrade a ESRI addin from 10.1 to 10?

text: I have created an esri addin which targeted ArcGIS 10.1, however I need to have the addin to work on ArcGIS 10, what needs to be done within my solution to get it working? I have checked the references, I have changed the targeted version witihin the config.xml, but when I try and run the addin (from Debug mode in VS2010) I don't see the addin as I would on the 10.1 version. Thanks


title: How can I select enclosed features within the same layer?

text: My problem: I have a single shapefile from a habitat-survey and need to select all hedge-features that are enclosed in meadows.The features that are to be selected were created by filling rings. Select by attribute does not work since not all hedges are within meadows. I've tried the Spatial Query Plugin's 'within' and 'contains' option after duplicating the layer, but since the features do not overlap, the plugin finds no result. The 'touches' option gets results, but unfortunately also selects neighbouring features which I don't need. The question now is, how I can select only the enclosed objects without having to do it manually? Grateful for any help, Gandalf


title: how to calculate area of triangles formed by delaunay triangulation

text: After applying delaunay triangulation 160 triangles were formed and i need to calculate the area of each triangles. Kindly help us please.


title: QGIS Clipping and reprojecting SRTM data creates a HUGE file

text: I'm having difficulty in QGIS doing what should be a very simple thing and I haven't been able to find anyone who has had the same issue. **What I'm trying to accomplish:** I need to export a very specific clipped region of elevation data from SRTM tiles in .bin format. It needs to be in my local UTM zone projection. **Problem:** This should be simple but the process I'm using creates a runaway file size that fills my hard disk. **My process is:**   1. Download SRTM tiles (SRTM3 data version 2.1 from http://dds.cr.usgs.gov/srtm/version2_1/SRTM3/North_America/   2. Import as layers of raster data. CRS is WGS84   3. Change project CRS to UTM zone 18 (EPSG:32618) enabling "reproject on the fly". (Note - some of the data in in zone 17, but not the area I am trying to clip)   4. Raster->Misc->Merge. This works properly and creates a 23.1MB file which loads into the canvas as a layer. Remove all other layers to clean things up.   5. Raster->Extract->Clipper. Mode "Extent" Enter extent coordinates manually. The area is very small - perhaps 10 miles by 10 miles. The dialog appears to freeze but if you look at the output file it is building, it expands to gigabytes in size eventually filling the hard disk after an hour or two and generating an error. The file size is around 300 gigabytes when it errors. I suspect the problem is related to my merge/reprojection and subsequent extraction of the merged rasters. If I simply clip a small section from a single tile without reprojection or merge, it works instantly, creating a small file. However, my scenery crosses tiles and I cannot extract the area I wish without first merging the SRTM tiles together, as far as I can tell. I'm brand new to QGIS and GIS in general and am trying to use SRTM data as elevation data for a Condor soaring flight simulation scenery landscape. If anyone has a better suggestion or tool for exporting clipped SRTM data to the .bin format, I'd welcome it.


title: What are the pros/cons of using an MSSQL SDE Database rather than a File-gdb in ArcGIS Server

text: I have taken over an ArcSDE 9.3 Enterprise database setup which has around 80% of it´s vector data in an MSSQL SDE Database system and 20% in FileGDB on the application server. I am trying to understand if there is a logical reason as to why someone would set this up. Can anyone offer reasons why a MSSQL SDE Database is prefered to FileGeodatabases? Versioning is not being used in either case. The web clients which have read/write access seem to all be in the MSSQL SDE DB. There are also 6 SDE Services running although most of the connects are direct. If anyone could explaing the logics behind this (if there are any) I would be grateful.


title: Extract the spectral signature of Landsat 7 and Landsat 8

text: I have some polygons and I want to extract the spectral signature of the adjacent satellite images. The satellite images that i use are the NDVI from Landsat 7 and 8 and I want to compare the spectral signature of them. I use both ENVI 5.03 and Arcgis 10.2. Is enough to extract the values from the pixels?


title: Connecting MapInfo 12 to ArcGIS ImageServer

text: I'm trying to figure out how to connect to an ArgGIS ImageServer raster layer via MapInfo. I can successfully connect via WMS, however this is obviously rather slow due to the lack of any caching. Therefore instead of using a WMS connection, I'd prefer to use the ImageServer 'exportImage' request in the fashion of a tile cache. However, I'm not sure if this is possible with MapInfo. So for example, the URL I'd like to use to connect to is of the form:               http://<servername>/arcgis/rest/services/Basemaps/Topographic/ImageServer/exportImage?f=image&bbox=260000,5116000,632000,5662000&bboxSR=28355&size=600,600      So can I connect to this with MapInfo 12? Be great if I could just add it like another Tile Server like the standard support for Bing, etc. _NOTE:_ I did try to install ImageServer MapInfo client, however the installer reported that I did not have MapInfo installed and would not let me proceed - probably as it's only meant to support up to MapInfo 9.5.


title: Assigning ID to features with similar attributes

text: I have a polygon feature with three attribute columns. The values in each of the attribute columns are number coded. I want to generate an ID (fourth column) for the unique combination of the values of the attribute columns. Rows with similar attributes are assigned same ID ## Attrib1 Attrib2 Attrib3 ID 1 1 33125   2 2 34020   1 2 60014   1 1 40290   3 1 55041   1 3 55047   2 3 56008   1 1 56026   3 3 56032   1 1 56035   1 1 57103 There are several thousands rows that it is time consuming to identify the unique combinations and assign ID manually. Thanks in advance.


title: Geometry vs geography for line-to-polygon intersection

text: I have two ~~WGS84~~ NAD83 shapefiles, one polygonal, the other containing lines:               Geometry: Polygon     Feature Count: 1294     Extent: (-79.761059, 44.991143) - (-57.105486, 62.582865)     Layer SRS WKT:     GEOGCS["GCS_North_American_1983",         DATUM["North_American_Datum_1983",             SPHEROID["GRS_1980",6378137.0,298.257222101]],         PRIMEM["Greenwich",0.0],         UNIT["Degree",0.0174532925199433]]          Geometry: Line String     Feature Count: 427227     Extent: (-79.882212, 44.929534) - (-55.834430, 62.427856)     Layer SRS WKT:     GEOGCS["GCS_North_American_1983",         DATUM["North_American_Datum_1983",             SPHEROID["GRS_1980",6378137.0,298.257222101]],         PRIMEM["Greenwich",0.0],         UNIT["Degree",0.0174532925199433]]      which I export to PostGIS 2 using the `-G` flag, to obtain geographies:               $ shp2pgsql -G -I -s 4326 polygons | psql ...     $ shp2pgsql -G -I -s 4326 lines | psql ...      However when I perform a query to obtain their intersection:               select * from polygons p, lines l     where p.geog && l.geog      and st_intersects(p.geog, l.geog)      it runs for a very, very long time, and doesn't yield the same results as the one performed **over the equivalent tables, but with geometries instead** (i.e. `shp2pgsql` without the `-G` flag), which runs very rapidly. According to the docs, `st_intersects` should play well with both types, so why this difference?


title: OpenLayers - Introduction tutorial problem

text: I have just started looking into OpenLayers and I can't seem to even get the example provided to work, I just get a blank screen. Am I missing something really obvious? Here's the html I am trying to run. Many thanks!                           <title>OpenLayers Example</title>


title: GeoWebCache with Geoserver layer using CQL_FILTER

text: I'm trying to setup a tile layer with a CQL_FILTER parameter but without success until now. I'm using GeoWebCache 1.5.0 integrated in Geoserver 2.4.4. I've set the CQL_FILTER parameter and the geowebcache complained about the FILTER parameter that I don't use, so I've set the FILTER parameter too and now i get every time cache miss without any reason message. What can I be doing wrong? Here's the parameter filters definition:                   <regexParameterFilter>           <key>CQL_FILTER</key>           <defaultValue>INCLUDE</defaultValue>           <regex>^(.*)$</regex>         </regexParameterFilter>       <regexParameterFilter>           <key>FILTER</key>           <defaultValue>TRUE</defaultValue>           <regex>^(.*)$</regex>         </regexParameterFilter>      I've successfully managed to setup layers that don't receive any filter parameter. And by the way, setting this up in the GeoWebCache 1.5.1 standalone works just fine, why this doesn't work on the integrated version? Note that in the integrated version the cache layer has been specified in the LayerInfoImpl*.xml instead of the geowebcache.xml. Also note that the same geowebcache.xml from the standalone version doesn't work well in the integrated version, giving me this error:               2014-02-12 16:40:03,663 ERROR [geoserver.ows] -      java.lang.IllegalStateException         at com.google.common.base.Preconditions.checkState(Preconditions.java:129)         at org.geoserver.gwc.wms.CachingWebMapService.invoke(CachingWebMapService.java:92)         at org.geoserver.gwc.wms.CachingWebMapService.invoke(CachingWebMapService.java:54)         at org.springframework.aop.framework.ReflectiveMethodInvocation.proceed(ReflectiveMethodInvocation.java:172)         at org.geoserver.ows.util.RequestObjectLogger.invoke(RequestObjectLogger.java:54)         at org.springframework.aop.framework.ReflectiveMethodInvocation.proceed(ReflectiveMethodInvocation.java:172)         at org.springframework.aop.framework.JdkDynamicAopProxy.invoke(JdkDynamicAopProxy.java:202)         at $Proxy58.getMap(Unknown Source)         at sun.reflect.GeneratedMethodAccessor159.invoke(Unknown Source)         at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)         at java.lang.reflect.Method.invoke(Unknown Source)         at org.geoserver.ows.Dispatcher.execute(Dispatcher.java:774)         at org.geoserver.ows.Dispatcher.handleRequestInternal(Dispatcher.java:272)         at org.springframework.web.servlet.mvc.AbstractController.handleRequest(AbstractController.java:153)         at org.springframework.web.servlet.mvc.SimpleControllerHandlerAdapter.handle(SimpleControllerHandlerAdapter.java:48)         at org.springframework.web.servlet.DispatcherServlet.doDispatch(DispatcherServlet.java:923)         at org.springframework.web.servlet.DispatcherServlet.doService(DispatcherServlet.java:852)         at org.springframework.web.servlet.FrameworkServlet.processRequest(FrameworkServlet.java:882)         at org.springframework.web.servlet.FrameworkServlet.doPost(FrameworkServlet.java:789)         at javax.servlet.http.HttpServlet.service(HttpServlet.java:643)         at javax.servlet.http.HttpServlet.service(HttpServlet.java:723)         at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:290)         at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:206)         at org.geoserver.filters.ThreadLocalsCleanupFilter.doFilter(ThreadLocalsCleanupFilter.java:27)         at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:235)         at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:206)         at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:74)         at org.geoserver.filters.SpringDelegatingFilter.doFilter(SpringDelegatingFilter.java:45)         at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:235)         at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:206)         at org.geoserver.platform.AdvancedDispatchFilter.doFilter(AdvancedDispatchFilter.java:49)         at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:235)         at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:206)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:311)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.invoke(FilterSecurityInterceptor.java:116)         at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.doFilter(FilterSecurityInterceptor.java:83)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.access.ExceptionTranslationFilter.doFilter(ExceptionTranslationFilter.java:113)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerAnonymousAuthenticationFilter.doFilter(GeoServerAnonymousAuthenticationFilter.java:53)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.authentication.www.BasicAuthenticationFilter.doFilter(BasicAuthenticationFilter.java:150)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.geoserver.security.filter.GeoServerBasicAuthenticationFilter.doFilter(GeoServerBasicAuthenticationFilter.java:82)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.context.SecurityContextPersistenceFilter.doFilter(SecurityContextPersistenceFilter.java:87)         at org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1.doFilter(GeoServerSecurityContextPersistenceFilter.java:52)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.springframework.security.web.FilterChainProxy.doFilter(FilterChainProxy.java:173)         at org.geoserver.security.GeoServerSecurityFilterChainProxy.doFilter(GeoServerSecurityFilterChainProxy.java:134)         at org.springframework.web.filter.DelegatingFilterProxy.invokeDelegate(DelegatingFilterProxy.java:346)         at org.springframework.web.filter.DelegatingFilterProxy.doFilter(DelegatingFilterProxy.java:259)         at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:235)         at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:206)         at org.geoserver.filters.LoggingFilter.doFilter(LoggingFilter.java:75)         at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:235)         at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:206)         at org.geoserver.filters.GZIPFilter.doFilter(GZIPFilter.java:48)         at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:235)         at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:206)         at org.geoserver.filters.SessionDebugFilter.doFilter(SessionDebugFilter.java:47)         at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:235)         at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:206)         at org.geoserver.filters.FlushSafeFilter.doFilter(FlushSafeFilter.java:43)         at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:235)         at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:206)         at org.vfny.geoserver.filters.SetCharacterEncodingFilter.doFilter(SetCharacterEncodingFilter.java:109)         at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:235)         at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:206)         at org.apache.catalina.core.StandardWrapperValve.invoke(StandardWrapperValve.java:233)         at org.apache.catalina.core.StandardContextValve.invoke(StandardContextValve.java:191)         at org.apache.catalina.core.StandardHostValve.invoke(StandardHostValve.java:127)         at org.apache.catalina.valves.ErrorReportValve.invoke(ErrorReportValve.java:103)         at org.apache.catalina.core.StandardEngineValve.invoke(StandardEngineValve.java:109)         at org.apache.catalina.connector.CoyoteAdapter.service(CoyoteAdapter.java:293)         at org.apache.coyote.http11.Http11AprProcessor.process(Http11AprProcessor.java:879)         at org.apache.coyote.http11.Http11AprProtocol$Http11ConnectionHandler.process(Http11AprProtocol.java:617)         at org.apache.tomcat.util.net.AprEndpoint$Worker.run(AprEndpoint.java:1760)         at java.lang.Thread.run(Unknown Source)      Thanks!


title: How to determine and send alert SMS for next bus stop in a GPS tracking system?

text: In my gps tracking system (PHP/MySql) I have got 5 Stops for a vehicle. Whenever the bus reaches stop 1, the passengers at stop 2 should get an SMS alert. And in the same way for the sops followed by. I have given stops as points. How can I determine whether the bus reached a particular stop? Should always monitor all 6 stop whether the bus reached once the trip is started? Can anybody give some idea? (Map used: Google Map/Openlayers) See the image. ![Image showing 6 Bus Stops and start & stop points](http://i.stack.imgur.com/mSjmT.png)


title: Force an ArcGIS API for JavaScript WMSLayer to use EPSG: 900913 (instead of 102100)

text:
I have a use case where an external WMS server can serve maps in the Web Mercator projection using EPSG code 900913, but won't recognize ESRI's 102100 code. My question is: **can I force all WMS requests to use EPSG code 900913?** To illustrate, I have a simple map page:               <!DOCTYPE html>                           <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">             <title>Test</title>                  <link rel="stylesheet" href="http://js.arcgis.com/3.8/js/esri/css/esri.css">

                    Where test.js creates a simple map and adds a WMSLayer to it:               var wmsAddress = "http://some.server.com/wms";     var map, srs, extent;          require(["esri/map",              "esri/layers/WMSLayer",              "esri/layers/WMSLayerInfo",              "esri/geometry/Extent",              "esri/geometry/Point",              "esri/SpatialReference",              "dojo/domReady!"],             function (Map, WMSLayer, WMSLayerInfo, Extent, Point, SpatialReference) {              srs = new SpatialReference({wkid:900913});         extent = new Extent(                 1362550.0288221193, 7616458.979152972,                 2281987.5791366613, 9525359.127428867,                 srs);         map = new Map("mapDiv", {             center: new Point(1800000, 7900000, srs),             extent: extent         });              var wmsService = new WMSLayer(wmsAddress, {             resourceInfo: {                 extent: extent,                 layerInfos: [new WMSLayerInfo({name:"layerName",title:"TestMap"})]             },             visibleLayers: ["layerName"]         });         map.addLayers([wmsService]);     });      When loading the application, the request that is sent to the WMS server looks like this (viewed in a browser console): > > http//some.server.com/wms?SERVICE=WMS&REQUEST=GetMap&FORMAT=image/png&TRANSPARENT=TRUE&STYLES=&VERSION=1.3.0&LAYERS=layerName&WIDTH=1904&HEIGHT=400& > **CRS=EPSG:102100** > &BBOX=-2743182.35289663,6945549.925862053,6343182.35289663,8854450.074137948 Which results in this response from the server:               <ServiceExceptionReport xmlns="http://www.opengis.net/ogc" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" version="1.3.0" xsi:schemaLocation="http://www.opengis.net/ogc http://some.server.com:80/geoserver/schemas/wms/1.3.0/exceptions_1_3_0.xsd">         <ServiceException code="InvalidCRS">             Error occurred decoding the espg code urn:x-ogc:def:crs:EPSG:102100 No code "EPSG:102100" from authority "European Petroleum Survey Group" found for object of type "IdentifiedObject".         </ServiceException>     </ServiceExceptionReport>


title: How can I create a hexbin layer using the SQL and CartoCSS parameters in cartodb.js?

text: CartoDB has a great hexbin visualization wizard baked into the client dashboard: ![enter image description here](http://i.stack.imgur.com/wkeGB.jpg) However, I like using cartodb.js to customize events and interactions; in particular, I'd like to make styles dependent on map zoom level. Recognizing that the hexbin visualization is a combination of SQL and CartoCSS parameter configuration, what are the ideal functions to pass to CartoDB to return hexes with cartodb.js?


title: How to change layer transparency without changing border lines?

text: I have a raster image of elevations. I have added a vector shapefile that contains local political divisions. I want to keep the elevation colors, but I want to keep the black lines representing the political divisions as well. I have reduced the transparency of the political divisions shapefile, but I can't figure out how to have this NOT affect the political division lines. Advice for reducing the fill transparency while keeping the black border colors at full?


title: How to identify line intersection in QGIS when I have more than 2 lines?

text: I am trying to identify where overlaps exist between multiple shapefiles. I would use the intersect tool, but it only allows me to intersect two layers. I have over one hundred individual line shapefiles and I want line segments that overlap. Any ideas? I'm using QGIS 2.0.1 and GRASS for the creation of the shapefiles.


title: How to change the style of the north arrow from “ESRI North 4” to “ESRI North 1”?

text: I have a Python assignment where I have to write a script to change the style of the north arrow from “ESRI North 4” to “ESRI North 1”. I wrote the following code:               import arcpy     mxd = arcpy.mapping.MapDocument(r"F:/Lab9/Georgia.mxd")          for elm in arcpy.mapping.ListLayoutElements(mxd, "MAPSURROUND_ELEMENT"):              if elm.name == "North Arrow":             elm.style = "ESRI North 1"          mxd.save()          del mxd      It gave me no error, but It didn't change the existing style to the one I required.


What are the valid values for "band" in a GeoTiff file? I've looked through the Geotiff Spec at http://www.remotesensing.org/geotiff/spec/geotiffhome.html and I am having trouble finding information about multi-band imagery. I'm looking for documentation that answers these questions:   * Do the band values always start at 1?    * Is band 0 a valid band?   * What is the max band allowed? I am looking at the Band information that is reported from running gdalinfo on a geotiff file.


title: What's the correct csv format for MMQGIS Join by attribute?

text: When trying to join attributes from CSV-file with MMQGIS, I currently get the error message: * * * File inaccessible or in unrecognized CSV format /Netzplanung/Truck_park/trkbuspk15_2.csv * * * Tried with semicolons and commas as seperators, didn't work. Using Qgis master 1.9.0, MMQGIS 2012.12.07. Anybody got an idea? Attached sample csv file (semicolons as separators) sample file


title: Importing shapefile in grass

text: I have very simple shapefile, which consists of couple of polygons, and I wanted to snap boundaries automatically. Unfortunately qgis doesn't have such function, but grass has according this answer: QGIS :Check Topology I followed provided walk-through, i.e. I made new grass mapset, then imported shapefile layer from qgis, but when I wanted to apply `v.clean.break` or `v.clean.snap` command, I couldn't select input vector map as there was not option to select one (FYI I also did refresh in Browser tab): ![enter image description here](http://i.stack.imgur.com/DCAXi.png) Browser shows that there is one, but all values are zero, so I wonder if I made a mistake somewhere: ![enter image description here](http://i.stack.imgur.com/xvB2x.png) I then deleted created mapset and closed qgis, then opened grass to do the same task but now by just using grass. I imported shapefile to vector mapset but it's just same issue, it's like mapset is empty or something like that. I use up to date OSGeo4W environment. What could be issue here? * * * **Update** :               (Tue Dec 11 17:09:11 2012)                                                           v.in.ogr dsn=C:\Temp\geo\1.shp layer=1 output=x      Projection of input dataset and current location appear to match     Layer: 1     Counting polygons for 11 features...     Importing map 11 features...     -----------------------------------------------------     Building topology for vector map <x_tmp>...     Registering primitives...     0 primitives registered     0 vertices registered     Number of nodes: 0     Number of primitives: 0     Number of points: 0     Number of lines: 0     Number of boundaries: 0     Number of centroids: 0     Number of areas: -     Number of isles: -     Building topology for vector map <x>...     Registering primitives...     0 primitives registered     0 vertices registered     Building areas...     0 areas built     0 isles built     Attaching islands...     Attaching centroids...     Number of nodes: 0     Number of primitives: 0     Number of points: 0     Number of lines: 0     Number of boundaries: 0     Number of centroids: 0     Number of areas: 0     Number of isles: 0     (Tue Dec 11 17:09:13 2012) Command finished (1 sec)