title: get linestring and polygon into leaflet from geoserver

text:
I can load points from my geoserver to leaflet via jsonp. Now I would like to know how can I load linestring or polygon with similar method into leaflet. Here is the example for points:               var owsrootUrl = 'http://217.8.255.188:8080/geoserver/opengeo/ows';          var defaultParameters = {         service : 'WFS',         version : '2.0',         request : 'GetFeature',         typeName : 'opengeo:csg_bt',         outputFormat : 'text/javascript',         format_options : 'callback:getJson',         SrsName : 'EPSG:4326'     };                     var geoJsonLayer = new L.GeoJSON();                function getJson(data) {               geoJsonLayer.addData(data);           }               var parameters = L.Util.extend(defaultParameters);     var URL = owsrootUrl + L.Util.getParamString(parameters);          var markers = L.markerClusterGroup();     var WFSLayer = null;     var ajax = $.ajax({         url : URL,         dataType : 'jsonp',         jsonpCallback : 'getJson',         success : function (response) {            WFSLayer = L.geoJson(response, {                 style: function (feature) {                     return {                         stroke: false,                         fillColor: 'FFFFFF',                         fillOpacity: 0                     };     },                 onEachFeature: function (feature, layer) {                     popupOptions = {maxWidth: 600};                     layer.bindLabel('
'+feature.properties.musno+'
');                     //sidebar.setContent('
'+feature.properties.musno+'

'+'
'+feature.properties.exchange_name+'

'+feature.properties.pcp, popupOptions);                     layer.bindPopup('
'+feature.properties.musno+'

'+'
'+feature.properties.exchange_name+'

'+feature.properties.pcp                         ,popupOptions);                 }             }).addTo(markers.addTo(map));         }     });


title: TileMill for Windows

text: After hours and days of trying to build TileMill from source on Windows I'm more than frustrated. Is there anyone how is willing to share his up-to-date TileMill for Windows with me? I tried to build a version from a few weeks ago. But basically I'm only interested in a version which supports setting a transparent color for raster data (which current 0.10.1 for Windows does not). It's "image- filters:colorize-alpha", that I'm interested in. Thanks in advance, Daniel


title: How can i turn off the minus to zoom control in OpenLayers while using OpenLayers.Control.KeyboardDefaults?

text: I have a map using OpenLayers.Control.KeyboardDefaults but I also have a form next to the map that acepts coordinates. As is stands pressing '-' for a negative coordinate does not type but zooms the map out. Is there a parameter I can add between the brackets of `new OpenLayers.Control.KeyboardDefaults()` that will allow me to turn off the listener for '-' while retaining all the other controls? Thanks


title: How do I specify an extent to my map without disabling zooming in?

text: I'm a bit new to GIS in general so apologies for the newbish question but in ArcGIS 10 I'm trying to get it so my map is only for my local city government (where all my data is) and a small surrounding area. I'd like to have it set up so that by dragging the scroll bars I'm confined to that area instead of the current scrolling around the entire earth. It renders the scroll bars useless besides a quick arrow nudge here or there. I've tried different extent settings on my data frame from this page from ESRI to no avail. I want to still be able to zoom in to inspect features closer but all viewing extent settings seem to disable zooming in and out. Is my base map causing a problem since it's worldwide imagery? All my other data is confined to my city and some small surrounding area. Am I out of luck in doing this? Thanks in advance!


title: Cannot import Grass in Console of Qgis (1.9.0 Master)

text: I can't import Grass in Console of Qgis (1.9.0 Master). I am guessing it is because of some Grass environment variables which are not set up correctly. However I don't understand why it is working when I start my script directly from the script editor (new feature in QGIS 1.9 it seems). How come environment variables could be set up differently for the console and for the editor? Thanks in advance for any explanation and help with this. I would like my script start from the console (with command execfile) > import grass.script as grass               Traceback (most recent call last):       File "<input>", line 1, in <module>       File "C:/Program Files (x86)/QGIS Weekly/apps/qgis-dev/./python\qgis\utils.py", line 453, in _import         mod = _builtin_import(name, globals, locals, fromlist, level)     ImportError: No module named script


title: How can i display the color of my labels into the content table?

text: The thing is that i need to display labels based on colors but without the symbol. In this case i'm working with points in Arcgis 10 and the maps needs to display only labels with color-based categories and in the content table the layers display only the symbol. My question is if its possible to, somehow, the color of my labels to be displayed in the content table? I hope you can help me on this.


title: How to create "two-dimensional" color ramps in QGIS?

text: Is there a way in QGIS to create a "two-dimensional" color ramp for the graduated symbology; i.e. a color ramp that determines a color based on two different attributes? My situation is that I want to control the hue of the color of a polygon using one attribute and the saturation of the color using another attribute.


title: postgis: problems creating flow direction polyline; not all needed connections are drawn

text: I am using post-gis to create a polyline that describes the flow direction of water in a landscape/catchment, to use in a hydrological model. To do that I start with a regular grid (as polygon) and connect the centre of each cell to the centre each of the 8 cells around it, using ST_Intersect. (afterwards I delete the irrelevant lines, based on the elevation of the grid). When I do this, most connections are made as I want, except that only half of the diagonal connections are made: ![enter image description here](http://i.stack.imgur.com/Lo3G3.png) In the figure the green squares is my grid, and the white lines are the output polyline here is a code sample I use:               INSERT INTO test_olaf.olafmain(id,id_in,id_out,ngf_in,ngf_out,the_geom,geom_in,geom_out SELECT DISTINCT ON (ST_Intersection(parcel_a.the_geom,parcel_b.the_geom) CAST(CAST(parcel_a.id as TEXT)||CAST(parcel_b.id as TEXT) as INTEGER),      nextval('olaf_s2s_id'),     parcel_a.id,      parcel_b.id,     parcel_a.elevation,     parcel_b.elevation,     ST_MakeLine(ST_centroid(parcel_a.the_geom),ST_centroid(parcel_b.the_geom)),     parcel_a.the_geom,     parcel_b.the_geom           FROM test_olaf.grid as parcel_a, test_olaf.grid as parcel_b -- test_olaf.grid is my grid (based on a DEM, converted into a polygon)     WHERE ST_Intersects(parcel_a.the_geom,parcel_b.the_geom)          )     ;      Does anyone know how I can connect my cell-centre to all the cells around it, or why the ST_intersect function does not work for all the cells? thanks! Annette


title: ArcMap AddIn Build Error when referencing external ActiveX Interface in VS2010

text: If the "Register for COM interop" is unchecked the projects builds without errors. ![VS2010 Project Properties Build Tab](http://i.stack.imgur.com/uQPjq.png) But I have to reference external ActiveX Interface in my project so I have to have that check-box checked. I am getting the build error: Cannot register assembly "...\MyCustomArcMapAddIn.dll". Could not load file or assembly 'ESRI.ArcGIS.Desktop.AddIns, Version=10.1.0.0, Culture=neutral, PublicKeyToken=8fc3cc631e44ad86' or one of its dependencies. The system cannot find the file specified. Is there workaround here? Thanks


title: Transit webapp development (journey planner)

text: I am going to develop a journey planner (like Google transit), i was just wondering that if there any opensource API available that can help me, or i have to do every thing on my own. And where can i find some help to get started. kindly help me out . Thanks.


title: How to get Esri *.gdb file into QGIS on Mac?

text: I have tried the above options, but not being too computer savy, some I didnt get. I tried the homebrew option, but terminal says an error occurs when trying to install. I need this 1 file for a work project. Anything will help. Thank you!


title: Search in Bounding Boxes and Locations

text: Description of the problem (what I need): I have a set of bounding boxes (maybe the actual data is a location[longitude, latitude] in which case the bounding box will have identical corners). I have a search box (it is a bounding box too). Now I want to find which bounding boxes in my set intersect with my search box. Question: Which databases/tools (I prefer open source ones and among them I prefer free ones) can do this? Note: Currently I am using an in-memory RTree that I have developed. But my data is increasing and maintaining it; is becoming hard.


title: Is it better to define object properties or layer properties?

text: I have a question, but it is difficult to phrase. I am picking up GIS methods as I go, and I have run into this problem frequently. Please excuse me if I am in the wrong place. When organizing information, is it better to group Features together and display/arrange accordingly? Or is it better to attach specific information to each Feature individually? **For a real-world example** , I have lots of Features (>10000). Each Feature is given an icon, where the icon represents aspects of the feature. The map consists of layers that organize and display the Features. Take, for example, the question of where to store the icon information? Should I let the layers dictate the icons that are displayed? Or should I let the Features render themselves in whatever layers they are in? This is coming from a programmer, so processing speed/memory are the main issues. Is there a general design principle for this sort of behavior? If my question indicates a deeper misunderstanding, please fill me in. I've been running in circles.


title: Dynamically Change Annotation Font Size According To Parcel Area

text: I'm looking for a solution to change annotation font size in an automated/dynamic way. The annotation data represents civic addresses and are displayed over a parcel fabric. What we would like to do is have the font size change in proportion to the area of the parcel. The current dataset has several classes and this kind of lets us do what we want, but still, it doesn't cut it. In the past we have done this manually, but our quickly growing municipality makes this task rather arduous and cumbersome. I can't foresee a simple solution to this problem, but I'm willing to try anything. I have access to FME 2012 and ArcGIS 9.3.1. Scripting solutions are fair game. Any hints/suggestions are welcome. Cheers.


title: How do I select * rows of multiple columns in one query using Select by Attributes within ArcMap?

text: I am trying to convert a psql query string into arcpy formatting for an UpdateCursor where expression parameter:               '''SELECT \"NAME_0\",\"NAME_1\",\"NAME_2\",\"NAME_3\",\"NAME_4\",\"NAME_5\",\"WIKI_URL\", \"LEVEL_DEPT\" FROM      AdminBoundaries WHERE \"WIKI_URL\" IN(SELECT \"WIKI_URL\" FROM AdminBoundaries GROUP BY \"WIKI_URL\" HAVING      (COUNT (\"WIKI_URL\") > 1)) ORDER BY \"WIKI_URL\";     '''


title: How to convert a data frame grid to a shapefile?

text: I have a map layout set up in ArcMap 9.3.1 with an index grid on the data frame. What I would like to do is convert this index grid to a shapefile that includes the grid reference information. I've tried the Convert Graphics to Features button on the Drawing toolbar, but it is greyed out. Is anyone aware of a method or script to perform this operation?


title: How to automatically convert line features into polygons using QGIS/PostGIS?

text: I have a master table of line features and a number of tables that I use to join (in the SQL sense) to produce various layers. Some of theses layers consist of boundaries of administrative zones (think counties ;) Using SQL I can either show the outline of all the zones or subsets of the zones. Now I want to have a layer that consist of a set of polygons, one for each zone -- so I can use it with the atlas generator. I would prefer a solution that can be automated, preferably in an sql query. Some points:   * Currently the line segments end points _don't_ match exactly (how would one get them to?)   * I know about ST_BuildArea but I think it expects the lines to form perfect closed polygons. When I try and use it I get back n empty polygons where n is the number of lines in the input. What I am looking for is something like ST_BuildArea which will allow some tolerance in the line end points. Another way of coping with this would be to have a method of standardising endpoints of lines given some tolerance Does such a thing exist? Before I go and start re-inventing the wheel. I figure it should be a fairly easy plugin to write.


title: Plugin QGIS Error "QTextEdit.setPlainText(QString): not enough arguments"

text: I try to make save as button but i got this error message when try to save file in my directory               File "C:\Users\Mr.Pakde/.qgis2/python/plugins\condition\conditiondialog.py", line 80, in saveAs         self._save( _filename )     File "C:\Users\Mr.Pakde/.qgis2/python/plugins\condition\conditiondialog.py", line 68, in _save         f.write( 'Rekomendasi :' + str(self.ui.textRec1.setPlainText() ) )     TypeError: QTextEdit.setPlainText(QString): not enough arguments      this is my code                       Smpan = self.ui.btnSave             QtCore.QObject.connect(Smpan, QtCore.SIGNAL('clicked()'),self.saveAs)                     def Compare1(self):             input1 = self.ui.lineInput1.text()             input2 = self.ui.lineInput2.text()             compare = ''             if input1 == input2:                 compare = 'Data dapat digunakam'             else:                 compare = 'Data tidak cocok'             self.ui.textRec1.setPlainText(compare)              def _save(self, simpan):                     f = open( simpan, 'w' )             f.write( 'Input 1 :' + str(self.ui.lineInput1.text()) )              f.write( 'Input 2 :' + str(self.ui.lineInput2.text()) )             f.write( 'Rekomendasi :' + str(self.ui.textRec1.setPlainText()) )                    f.close()         def savefile(self):             if self.simpan:                   self._save( "%s" % self.simpan )             else:                   self.saveAs()              def saveAs(self):                 tulis = QtGui.QFileDialog(self).getSaveFileName()                 if tulis !="":                     _filename = "%s" % tulis                     self._save( _filename )


title: Can GeoProcessor use credentials?

text: I've set up a web service that uses ESRI's REST API spec. This allows me to make requests using the (Silverlight)/WPF SDK. If I turn on windows authentication though, I get a 401 error when using GeoProcessor.ExecuteAsynch. I can successfully call the same url with a WebClient with credentials set:               WebClient c = new WebClient();     c.Credentials = CredentialCache.DefaultCredentials;      If I use WebClient, then I have to deserialize ESRI's json, which is non- trivial, so I'd really like to use Geoprocessor.ExecuteAsynch instead. Is there some way to send credentials with Geoprocessor?


title: How to split a line at designated points and vertices(confluences)?

text: I have a line that already has its own segments and attribute table： ![enter image description here](http://i.stack.imgur.com/UuvNH.png) ![enter image description here](http://i.stack.imgur.com/lOj80.png) I've used **line dissolve** tool in **Geoalgorithms** to make the previous attribute table look like this：so I make the segmented line become an integrated one. ![enter image description here](http://i.stack.imgur.com/Okv8B.png) And now I want to **re-clip** this line with：   1. its confluences   2. designated points the designated points are locations on the line and distribute like this： ![enter image description here](http://i.stack.imgur.com/ElN7A.png) How can I use these two features to clip this line and get re-arranged segments, new attribute table？ I am working on Win 8.1 64-bits, Qgis 2.0. Thanks for giving me advices.


title: Is there a way to get ST_ClosestPoint for Geography type?

text: Using PostGIS 2.0.0 I'd like to find the point on a LINESTRING that is closest to a given point. The LINESTRING represents a great circle line (ie. geography type). ST_ClosestPoint appears to do exactly what I want, however, I find that the returned point does not lie on the great circle line:               WITH points AS (         SELECT ST_GeomFromEWKT('SRID=4326;POINT(41 12)') AS point,             ST_GeomFromEWKT('SRID=4326;LINESTRING(40 5,41 15)') AS border     )     SELECT ST_AsText(ST_ClosestPoint(border, point)), -- POINT(40.7029702970297 12.029702970297)         ST_DWithin(ST_ClosestPoint(border, point)::geography, border::geography, 700) -- false     FROM points;      It's actually over 700 metres from the line. I believe this happens because ST_ClosestPoint operates on a plane. How would I do the same thing on the spheroid (or at least on the sphere)? **Edit:** I tried MerseyViking's function (ported to plpgsql) and it mostly works very well, but in some cases returns a point that is not on the line. A simple example:               WITH points AS (         SELECT ST_GeomFromEWKT('SRID=4326;POINT(10 3)') AS point,         ST_GeomFromEWKT('SRID=4326;LINESTRING(10 5,10 15)') AS border     ),     calc AS     (         SELECT ST_ClosestPoint(border, point)::geography AS plane_closest,             closest_point_to_line(border, point) AS sphere_closest         FROM points     )     SELECT ST_AsText(plane_closest) AS plane_closest, -- POINT(10 5)         ST_DWithin(plane_closest, border::geography, 1), -- true         ST_AsText(sphere_closest) AS sphere_closest, -- POINT(10 3)         ST_DWithin(sphere_closest, border::geography, 1), -- false     FROM points, calc;


title: What kind of high level semantic free GIS layers out there?

text: I'm looking for free GIS data that contains some high level information such as type of vegetation or type of soil or some other visual characteristics of the neighborhoods. I need more information than Google/Bing/Yahoo Maps. Do you know any GIS server that has this kind of information in it? Thanks,


title: ArcGIS 9.3 Geoprocessor - AttributeError: Object: Environment <workspace> not found

text: So, this could be the most puzzling python geoprocessor error I have ever dealt with... This problem is occurring in ArcGIS Server 9.3, Python 2.5.2 I am exposing a typical python script as a geoprocessing service. In this script, I frequently use the gp.ListFeatureClasses() method, which searches the current workspace as shown in snippet below.               import arcgisscripting     gp = arcgisscripting.create(9.3)     gp.workspace = "C:/someworkspace/mygdb.gdb"      featureClasses = gp.ListFeatureClasses()      In the initial portion of the script, everything works as expected. But in later portions of the script, the service crashes when I call gp.workspace reporting the following error: > `AttributeError: Object: Environment <workspace> not found` My initial guess was that maybe the workspace was being set to a non-existent workspace or gdb, but I have been able to verify that the workspace is always infact valid. Digging around a bit has shown that others have experienced the same problem. However, none have been able to set the workspace AT ALL. Has anyone ever experienced such issues?


title: Delete points outside a polygon with PostGis

text: I have to use PostGis because I have a very large database (more than 7 000 000 lines), but I am not used to work with PostGis; and since days, I would like to do that: Select and then delete all points which are geographically outside from a polygon (the points are in a difference table than the polygon table). I read many forums etc, but I don't find what I want. Did somebody have an advice or an example command which I can use? I tried to do it through QGis, but the treatment is too large and I have to wait a long time to finally obtain an error message! Cheers


title: Trouble trying to zoom to graphic (ArcGIS 9.3.1, VBA)

text: I am writing an application where one of the things it does is create a point then creates buffer (which is a graphic) at the specific point, what I want to do is when the buffer has been created, is to zoom to the extent of the buffer. I have tried the following code but it doesn't seem to work               Dim penv2 As IEnvelope     Set penv2 = New Envelope     Set penv2 = pActiveView.extent     Dim pArea As IArea     Set pArea = pElement.Geometry.Envelope     penv2.CenterAt pArea.Centroid                 pActiveView.Refresh      I've also tried               Dim penv2 As IEnvelope     Set penv2 = pElement.Geometry.Envelope     pActiveView.Extent = penv2      pActiveView.Refresh       But nothing happens. Any ideas on where I am going wrong? Thanks


title: How to fix "Error, query must return columns 'id', 'source', 'target' and 'cost'"?

text: I have one spatial table with **gid** , **FNODE_** , **TNODE_** and **LENGTH** column in my Postgresql/PostGIS database. I am using pgRouting. Problem is, when i try with shortest path query it shows following error. But it should not be. I double check my geographical database there is a connection between 505793 and 496514 nodes. But, whenever i query it shows following error.               SELECT * FROM shortest_path('                     SELECT "gid" as id,                              "FNODE_"::integer,                              "TNODE_"::integer,                              "LENGTH"::double precision as cost                             FROM roads_dist',                     505793, 496514, false, false);                    ERROR:  Error, query must return columns 'id', 'source', 'target' and 'cost'               ********** Error **********          ERROR: Error, query must return columns 'id', 'source', 'target' and 'cost'     SQL state: XX000      Thanks, in advance.


title: ArcGIS Symbology Transparency

text: Using ArcGIS 10.1, I am interested in symbolizing a series of polylines to generate a pseudo-heatmap like one of the answers to this question. The essence of the solution idea is to symbolize each polyline as red center, green border, and set the _individual line transparency_ to 20%. Using ArcGIS, it is trivial to set the entire layer transparency, but then the cumulative, low alpha, overlaps never sum to an opaque level. Is it possible to set opacity / transparency at the symbol level and not the layer level? Similar to How to display transparent overlapping polygons in ArcGIS?, but I am interested in polyline and using color, instead of hatching. Edit * * * I have tested the ability to do this in two ways:   1. If I buffer the polyline and add a transparency field, I can directly set transparency for each layer. Unfortunately, this still appears to 'flatten' the symbology, i.e. the cumulative alpha of overlapping polylines is lost.   2. By exporting 4 of the polylines, known to overlap, it was possible to manually set the transparency ( _on the layer_ ) to 80%. When all layers are viewed the cumulative alpha is displayed as anticipated, i.e. two overlapping lines are 'darker' (40% opaque). As an aside, I tested this using a google fusion table and the result is as expected, i.e. alpha is per line and cumulative for overlapping symbols. This can also be done in OpenLayers, but that removes the step from my ArcGIS workflow (and more importantly, my users ArcGIS workflow). ![enter image description here](http://i.stack.imgur.com/gqKw2.png) ## Edit 2 I was able to accomplish this manually in ArcGIS using the following workflow. This is non-ideal as it requires a duplication of the layer _n_ times, where _n_ is the number of features.   1. Copy the layer and paste to the same dataframe.   2. Define a definition query on a unique field. In my case I used a grouping: `"id" >= 20 AND "id" < 30.   3. Repeat 2 until all features are in their own layer...   4. Add the features to a New Group Layer   5. In Group Layer Properties select Symbol Levels   6. Check 'Draw this layer using the symbol levels specified below'   7. Manually check Merge for each layer I suspect that the above workflow can be automated as a python add-in. I will need to explore that next. Result with black background: ![enter image description here](http://i.stack.imgur.com/pG53Q.png)


title: Failed to convert multilinestring to linestring

text: I tried to convert some multilinestring data into linestring using:               SELECT ST_AsText(     ST_LineMerge(         ST_SnapToGrid(             ST_GeomFromText(the_geom),             0.001)         )     );      and               ST_AsText(ST_LineMerge(     ST_GeomFromText(the_geom)     )     )       function, but when executing a database some data showing the same multistring values below is one of them:               "MULTILINESTRING((76.4505052014918 29.3605577579601,76.4506218636754 29.3607478108812,76.4511480538787 29.3613277796464,76.4518103152086 29.3620711781354,76.4527672628892 29.3631772485864,76.4531572408356 29.3635752267738,76.4539256694916 29.364440449285,76.4547439604311 29.365369956029,76.4550240785647 29.3656618043844,76.4552426759001 29.3658855962542,76.4555650340103 29.3662556751861,76.4556087620657 29.3663556730666,76.455767298346 29.3668617755123,76.4557673174369 29.3669897943356,76.4558918607226 29.3675826328776,76.4560706995576 29.368422784851,76.4562862936152 29.3694406538537,76.4566033421227 29.3709822940838,76.4568395123292 29.372037249194,76.4570259582372 29.3729957081483,76.4571197426559 29.3734350731336,76.4572698354037 29.3739352007255,76.4573190500858 29.3741392890071,76.4573682490733 29.3742469314174,76.457344669698 29.3743867963434),(76.4573092877099 29.3746523684764,76.4572972911226 29.3747446298991,76.4571720791717 29.3751492332069,76.4570081738852 29.3753587826844,76.4567056895133 29.3756411728579,76.4564101418548 29.3757981626948,76.4561806358876 29.3758648580016,76.4559546435516 29.3759683021827,76.4558235061862 29.3760683197683,76.4552352631393 29.3766506122008,76.4546155360981 29.377308054817,76.4537679226558 29.378150072491,76.4537570213123 29.3782127134124))"          "MULTILINESTRING((80.8962792489216 26.8237892138329,80.8962792489216 26.8237892138329,80.8960534536316 26.8238581478127,80.8959042776709 26.8238877309904),(80.8959404040242 26.823976632038,80.8960326566971 26.8239310587069,80.8962703359495 26.8238515195471,80.896404030529 26.8238037960245,80.896404030529 26.8238037960245))"      Above there are two examples. How to convert those values?


title: Workflow Scripting in QGIS

text: I have created a hexagonal grid shapefile using MMQGIS and want to use this for hexagonal binning. I have another point shapefile which contains patient records for three different types of hospital activity (let's call them A,B and C) for two different years. What would be the best way of quickly generating point counts (via Points in Polygon) for each subset of the activity shapefile so that I ended up with 9 additional columns in my hexagonal grid shapefile?   1. Year 1, activity A   2. Year 2, activity A   3. Difference between years 1 and 2 for activity A   4. Year 1, activity B   5. Year 2, activity B   6. Difference between years 1 and 2 for activity B   7. Year 1, activity C   8. Year 2, activity C   9. Difference between years 1 and 2 for activity C Is it possible to script this using pyqgis? I've started reading up on pyqgis but I'm a long way off being able to tell if this is possible. If it is possible, are there any tutorials out there that would cover this sort of processing or at least point me in the right direction to start figuring it out? I really want to learn scripting in QGIS and this feels like something which would lend itself to getting started.


title: IOError trying to change PICTURE_ELEMENT sourceImage

text: I recently ran into this issue after I published a Geoprocessing Task to a production server. It appears to be the same problem discussed in this question, but there is no answer posted there yet. Hopefully, sharing the details of my circumstances will allow someone to provide an answer. My geoprocessing python code is an enhanced Web Map Print service that loads a map via the `arcpy.mapping.ConvertWebMapToMapDocument` function and either sets the `sourceImage` of a `PICTURE_ELEMENT` or moves it off of the page. The published GP Service works correctly on a development server, but returns the following error when running on the production server.               Unable to complete operation.     Error executing tool.:       Traceback (most recent call last):         File "<ags_folder>\PythonWebPrint.GPServer\extracted\v101\my_toolboxes\service.py", line 643, in <module> main()        File "<ags_folder>\PythonWebPrint.GPServer\extracted\v101\my_toolboxes\service.py", line 634, in main (message, output_file, image_files) = execute(web_map_json, output_format, layout_template, layout_template2, georef_info, quality)        File "<ags_folder>\PythonWebPrint.GPServer\extracted\v101\my_toolboxes\service.py", line 482, in execute staticLegend.sourceImage = legendPath        File "c:\program files\arcgis\server\arcpy\arcpy\arcobjects\_base.py", line 87, in _set return setattr(self._arc_object, attr_name, cval(val))       IOError: Cannot set path based on current settings.       Failed to execute (PythonWebPrint). Failed to execute (Python Web Print).      I have verified that the correct file path is being used, and that the file exists and that the folder and file permissions are identical to the MXD files being loaded by the `arcpy.mapping.ConvertWebMapToMapDocument` function. Any ideas as to the source of or how to fix this issue, or other diagnostics that I can perform?


title: Arcpad checkout to an area

text: I would like to check out data from an SDE database but for performance sake, i prefer to delete all features outside of my survey area. i've tried to copy features in area from sde to a file GDB then checkout , but the checkin could not be perfect, because i can only synchronize new added features to sde database. so i cannot synchronize modified or deleted features. i've tried to checkout directly from sde Database then delete all features outside area using Bounding box query , but it take much time once i start arcpad after checkout. i'm using Arcpy to checkin and checkout. Any help please.


title: How to extract an adjacency matrix from OSM?

text: I would like to use graph algorithms from the Boost graph library on OSM data. What are my choices to use Openstreetmaps data from within a C++ program ? I need at least a flat file with a list of nodes and edges to build the graph using Boost graph. I am not looking for code that computes anything from the data. I'm thinking of something like an XML parser in C++ that I could integrate to my code.


title: Adding Rules to Topology with Arcpy: getting ERROR 000714

text: I want to make a very basic python script, however I keep running into error 000714 (general/unknown script error). The topology is created, the features are added to it, however the rule is not. I've tried different rules, same problem; I will want to add more rules and validation if this error is overcome. I have tried doing the same process in modelbuilder but the I get error 999999, the syntax is the same when I export the code from it. Anyone run into this problem or find solution? I am using ArcGIS 10.1. Code (The script tool's only parameter is a Dataset):               import arcpy, sys     from arcpy import env          env.workspace = sys.argv[1]          topolo = "GBA_Topology"          transLn = "TRANSPORT_LINE"     gbaTransLn = "GBA_TRANSPORT_LINES"          if arcpy.Exists(topolo):         arcpy.Delete_management(topolo)          arcpy.CreateTopology_management(env.workspace, topolo, "0.001")          arcpy.AddFeatureClassToTopology_management(topolo, transLn, "5", "1")     arcpy.AddFeatureClassToTopology_management(topolo, gbaTransLn, "5", "1")          arcpy.AddRuleToTopology_management(topolo, "Must Not Self-Intersect (Line)", transLn, "", "", "")


title: Using geojson hosted on GitHub from Leaflet: it doesn't work

text: I'm trying to use a geojson hosted on GitHub (rif. https://raw.githubusercontent.com/cesaregerbino/DatiGeo/master/ScuoleTorino.geojson), from Leaflet. here you are my code .... My first Leaflet map               <meta name="viewport" content="width=device-width, initial-scale=1.0">          <link rel="stylesheet" href="http://localhost/Leaflet/072/leaflet.css" />                                        The page doesn't show the markers of my geojson file: it seems that no return come from the request. Here you are the situation ![enter image description here](http://i.stack.imgur.com/k52wn.png) If I put the same geojson file under my apache it works fine. The new situation ... ![enter image description here](http://i.stack.imgur.com/xhGgE.png) I've tried to use the same geojson form QGIS desktopanche it works too. The image .... ![enter image description here](http://i.stack.imgur.com/M2fFk.jpg) Any suggestions? Thank you very much in advance ... Bye Cesare


title: mapnik verbose?

text: I've got a PostGIS NHD layer thats not rendering in mapnik (using toposm code scripts) and I don't know why. I've confirmed that the table data is there, positioned correctly and has the proper projection via QGIS, and the styles seem to have the table names correctly, but I'm still getting a features layer without my nhd water features. Is there a way I can tell mapnik.render to spit out the queries its running against my PostGIS db? Is there another way to debug this? TIA!


title: Merge Error: Not Enough Memory

text: I'm attempting to make a flood plain map for a school project, and am using FEMA Flood Hazard Symbology for the legend. I've dissolved by a Flood Zone Field, but am left with two attributes both named "X" they are slightly different, but for what I need it does not matter. I want to consolidate these attributes into one. So far I've tried to merge them using the editor toolbar. The process starts to run, but ArcMap locks up and eventually says the merge cannot be completed: Not enough memory. I'm running Arc on a computer with 8Gb or RAM. Is there another way to consolidate these into one, or a way to not make it lock up? ![enter image description here](http://i.stack.imgur.com/ZGbqM.jpg)


title: Buffer from Field using Different Unit

text: Is it possible to create a buffer polygon using an attribute field which is in different units than the feature class coordinate system? Example: I have a field that stores distances in miles (double data type) and the feature class is in state plane (units are feet). Is there a way to do this? I'm using Esri ArcGIS 10.2. Thanks in advance.


title: JSAPI 2.8 and Dojotoolkit 1.7.2? May?

text: I realize no one can give a definitive answer but is it reasonable to think the next release of ESRI JSAPI will be built on top of dojotookit 1.7? And will be released in the next 2-3 months? Any documention for preparing 1.6 code for the move to 1.7 that takes advantage of some of the loader features introduced in 1.6?


title: Help with Geographic Weighted Regression: Condo Prices

text: I was working on a GWR Problem in ArcGIS to figure out the potential real estate values of Condos in a city. I thought I was following everything correctly but I keep receiving the same error over and over again. I have attached a link to the dataset below.   The dependent variable is Finalprice. The independent (explanatory) variables are all the variables that are KM_* (distances to desirable features) and the DIST_TRECC (km to city center). It keeps throwing the error that explains that there is a major model problem with multicollinearity. I also did the OLS regression analysis and a link to those results is below as well. Everything I read said that the numbers from the OLS analysis show that the GWR should be working. I am at a loss and not sure what to do next. LINKS   DATASET: https://db.tt/9EUXfJy4   OLS Results: https://db.tt/LwS8pSIm


title: Update ArcMap links to ArcSDE Geodatabase features

text: I've seen questions about attemting to repair a .mxd file with broken links. What I would like to do is get ahead of this issue and "fix" it before it becomes a problem. I have my data in an ArcSDE geodatabase. The SDE is being navigated to a new machine in the near future so the database connection properties will change. Is there anything I can do ahead of time to avoid broken links and slow loading maps? Here is a similar but different question: How do you deal with broken MXD\LYR paths ? Thanks for any advice.


title: Can the concentration of values in a raster map be known?

text: I have a generated raster map containing different flood levels. Is there a grass script to determine the percentage of values in a raster map? I'm thinking of using r.mapcalc to solve this, but I'm not yet sure how to do it.


title: identifying single overlaping features from a GeoServer layer

text: I'm trying to use GeoServer to display features from a PostGis table the url will look similar to this. Example URL Localhost:8080/geoserver/cite/wms?service=WMS&version=1.1.0&request=GetMap&layers=cite:Layer1, cite:Layer2& styles=&bbox……..width=512,heaight=430&srs=EPSG:4326&application/openlayers               &FILTER= (<Filter xmlns="http://www.opengis.net/ogc">                    <FeatureId fid=”Layer1.1234”></FeaturId>                    <FeatureId fid=”Layer1.5678”></FeaturId></Filter>)                 (<Filter xmlns="http://www.opengis.net/ogc">             <FeatureId fid=”Layer2.345”></FeaturId></Filter>)      When I click on a polygon I need to highlight a single polygon if two (or more) polygons overlap i need to be able to step through each polygon in turn showing just it's attributes. At present I click on a feature and the attributes for the entire layer are shown in my example I use PostGis I've also tried this using a shapefile and the same thing happens. Any help would gratefully received ThanX Holly


title: Neighbours list extracted out of polygon regions

text: I've got a SQL database which contains some coded polygon structures. Those can be extracted as follows               poly <- data.frame(sqldf("SELECT ST_astext(geometry) FROM table"))      The data.frame 'poly' contains strings that now can be converted to real 'SpatialPolygons' objects as follows (for the first string)               realWKT(poly[1,1])      I can do the previous for each string, and save it in a vector               list <- c()     for (i in 1:100){        list <- c(list, readWKT(poly[i,1])     }      The last thing I want to do, is to create a neighbourhood list, based on all the SpatialPolygons by making use of the following function               poly2nb(list)      But sadly, this command results in the following error               Error: extends(class(pl), "SpatialPolygons") is not TRUE      I know that the problem has something to do with the classtype of the list, but I really don't see a way out.. Any help will be appreciated!


title: Why are there both square and rectangle world maps?

text: When i search for images of earth divides to longitude ans latitudes seems that the are 2 major types: ![one](http://i.stack.imgur.com/HenqM.png) is a perfect square and the latitudes in it are grow a different rate from the equator to up (and down) [mean, seems that the length from 0 to 66 is shorter than the way from 60 to the pole (90)] and the ![second](http://i.stack.imgur.com/4UlN4.gif) is rectangle (oblong) and it's latitudes are grow in the same rate all the way up (and down).


title: What tools are available for uploading gis data to a database?

text: I would like a nice graphical tool that allows me to take standard gis data such as shp files and kml files and upload them to a database such as PostGIS, MySql or Oracle. Are there any such tools? What do other people use to move gis data into a database?


title: How would I create my own geocoding web service

text: I would like to write a geocoding web service for genealogy in Java. Because it is for genealogy, I don't need to go down to the street level. I only need to go down to the county. Looking at genealogy data, misspelled names are very common. People also use a lot of non standard abbreviations. For instance for Baltimore County might look like:   1. Baltimore   2. Baltimore County   3. Baltimore Co   4. Baltimore Cty


title: ArcGIS Javascript: upgrade ComboBox to FIlteringSelect using ItemFileReadStore

text: I would like to upgrade a web app I am working on from v2.1 to v3.1 of the Javascript API (I have tried 3.2 but I am getting some very odd rending problems, so I'm settling with 3.1 for now...) One issue I have encountered is that a ComboBox populated with a store does not seem to fully function with 3.1. I was told I can fix this by replacing the ComboBox with a FIlteringSelect widget, but so far I haven't succeeded. Below is my code. I was hoping I wouldn't need to do much more than replace the ComboBox object with the FilteringSelect object, but that apparently isn't sufficient. Can somebody help me see what I am missing? Thanks, Jason JS (old code:)                     function populateSchoolMenu(results)      {         var schoolName;         var values = [];         var testVals = {};              var features = results.features;         dojo.forEach(features, function (feature) {             if ( (schoolLayerFlag == 0) || (schoolLayerFlag == 2) )             {                 schoolName = feature.attributes.Name;                  }             else if (schoolLayerFlag == 1)             {                 schoolName = feature.attributes.Site_name;                  }                       if (!testVals[schoolName])              {                 testVals[schoolName] = true;                 values.push({ name: schoolName });             }         });              var dataItems = {             identifier: 'name',             label: 'name',             items: values         };       var menuStore = new dojo.data.ItemFileReadStore({ data: dataItems });            dijit.byId("schoolMenu").store = menuStore;        }      HTML (old:)                           <select id = "schoolMenu"                 dojotype="dijit.form.ComboBox"                 value = "Select School"                 autoComplete="true" class="comboBoxClass"                 forceValidOption="false"                           onchange="executeSelectedQuery(this)">                 </select>


title: How to specify WGS84 CRS when importing a Shapefile into Idrisi?

text: Could you help me, please, to open shapefile from QGIS to Idrisi. I'm having problem answering to the Reference system, the project is using WGS84 but that is not in the Pick list. Idrisi is asking for a .ref file but I have .prj; .qpj; .shp and .shx files.


title: Ideas for Kids GIS Workshop

text: I'm working on getting an 8 hour workshop made for Middle and High Schoolers for GIS. I know it can be hard to grab their attention so I was wondering if you guys had any ideas about what we can do to keep them interested. We should be able to have access to ArcMap as well as ArcGIS online. If you guys have any ideas let me know.


title: Layer vs. Color Transparency with ArcGIS Javascript API

text: When creating services with ArcGIS Server, I have paid attention to the warning that Layer Transparency will make response time slower and that I should change to color transparency instead so I have been changing it when I "share" the service. Maybe I'm not doing this correctly because I'm not getting the result that I want or maybe I need to ignore the warning. Currently all my layers are in the same service and I can set opacity to some value 0.5 (and/or can set it individually for different layers e.g. just the layer on top) but the layers on top still cover layers on the bottom. I'd like to see the parcel boundaries underneath the other polygon layers.               mylayer[2].setOpacity(0.1);      Even if I set layer 2 to something almost invisible (like 0.1 as shown in code above; and as rendered below) it covers the parcels underneath (see parcel 0353 in the screenshot below). How should I be doing this instead? Do I need to ignore the warning? ![enter image description here](http://i.stack.imgur.com/InJlw.jpg)


title: Why does a z factor of 1 get turned into 0.999998 in my output

text: Just built a model that uses 3D Analyst - Surface Volume. I am trying to use a z factor of 1. Whenever I run the model my output table shows that a Z factor of 0.999998 was used instead of 1. I know this is a very small difference but even so, it is skewing my results somewhat. Really what I would like to know, is there a good explanation of why it gets changed to 0.999998 or is there a way to make it stay at 1? * * * Thanks again for your input. When I've tried changing the type of raster, the tool grays out the z factor all together and locks it in at 0.999998. So... I'm guessing there is really no way around this using integer or floating point data. I even tried converting to integer and using the raster calculator to multiply the data by 1000 to preserve the decimals in integer form, but was still stuck with the 0.999998 z factor. No dice.


title: USA Historical Satellite Imagery for an Online Application

text: I am looking for different options to add historical satellite imagery to our online-based application. We currently incorporate Bing Maps into our application, but cannot get historical imagery. I am essentially looking for something like Google Earth with its historical imagery slide bar. I am currently looking into using services from Google Earth and Digital Globe. I have also looked into NAIP imagery. I tried downloading the .sid files for one state and one year, but am having trouble. I downloaded the data from Geospatial Data Gateway in county files. I was trying to clip the county .sid files because of the noise around each .sid file. The noise that I am referring to is not the black border, but pixels that received none zero values in the border. In the end I need to have a seamless image. I am not sure where I would go from this spot of having clean county imagery. I need an opensource way to create a mosaic dataset. If I stick with this option, I think I would be looking at Terabytes of imagery. What are my options? What else is out there? What could I be doing differently?


title: How to add an attribution on a GeoJSON Layer from Leaflet?

text: I need to use a GeoJSON Layer on my Leaflet map. Here is a sample of my code:               function onEachFeature(feature, layer) {                           if (feature.properties && feature.properties.popupContent) {                              layer.bindPopup(feature.properties.popupContent);                           }                          }          myGeoJsonLayer = L.geoJson(data, {                                          pointToLayer: function (feature, latlng) {                                                        return L.circleMarker(latlng, geojsonMarkerOptions);                                          },                                          onEachFeature: onEachFeature                   });     myGeoJsonLayer.addTo(map);                              TOC.addOverlay(myGeoJsonLayer, "My GeoJSON Layer");      All is working. Now I'd like to add an attribution on my layer but how? Any suggestions are appreciated.


title: How to display multiple row counts with some varying information for one city?

text: My table has two rows for the same city, latitude and longitude but there are different counts, office names and office number designations. For example: Raleigh, NC; Long -78.6389, Lat 35.7806; Count 745; Office body shop; Office number 1 Raleigh, NC; Long -78.6389, Lat 35.7806; Count 3445; Office glass shop; Office number 2 I went to map view to display the information and turned on display count, office and office number. When I click the city it only shows one row's information and the other office is missing. What adjustments need to be made to display both offices that share the same location? I will also need to do this for other locations that share the same location but are designated with different office names and office numbers.


title: Is there a way to visualize TIN surfaces for web-mapping purposes?

text: I have created a TIN data set of electricity consumption. It is visually appealing, and often times TINS help the map reader understand geographical relationships better. Are there methods to display data in a 3D manner within Leaflet JS?


title: What can cause Geoserver to deliver a WMS tile so slowly?

text: I have a PostGIS database with a table containing about 300 000 simple polygons (table has spatial index). I have a Geoserver installed on the same server delivering this table with WMS. I have a response time of about 6 sec to receive a single WMS tile from that dataset. (I'm using firebug to get the getMap query url and parameters including the BBOX for a single tile) When querying that postgis database directly for extracting the same data using the same BBOX with the following query I receive the result in less than 1sec...               SELECT         a.*,         ST_Transform(a.geom, 3857)     FROM         role.my_table a     WHERE         a.geom && ST_Transform(ST_MakeEnvelope(-8079088.141629988,5721158.693088872,-8077865.1491774265,5722381.685541436, 3857), 4326);      So my question is what can be causing that slow response from Geoserver while PostGIS can extract and transform the data very fastly? Geoserver need 5 sec just to render it? It may worth mentionning that the same Geoserver deliver tiles very fastly from others tables with fewer rows... Thanks a lot


title: How to store Python Add-In menu customization/configuration in *.esriaddin file?

text: There is a gap in my understanding of Python Add-Ins. When I create a Python Add-In it is easy to see where buttons, tools and extensions are implemented as Classes in the *_addin.py script of the Install folder, and how these get written into the shareable *.esriaddin file. However, I find the ArcGIS 10.2 Online Help for **Creating an add-in menu** to be more than minimalist. By some trial and error I have figured out, within ArcMap, how to drag add-in menus, buttons and tools from the Customize dialog onto an add-in toolbar to achieve a menu like below. ![enter image description here](http://i.stack.imgur.com/RfW57.png) What I am unable to figure out is how this "layout" of the add-in toolbar, with its add-in menu and add-in button and add-in tool gets written into the *.esriaddin file. Is the way to get add-in menu customization/configuration written into the *.esriaddin file by editing the config.xml, and if so, is the syntax used to do that documented anywhere?


title: ArcGIS Server: Execute Export Web Map Task with custom image service rendering rule

text: I am trying to use the out-of-the-box _Print_ ==> _Export Web Map Task_ for ArcGIS Server (10.1+) a web map containing an image service layer with a custom rendering rule. The web map json spec mentions the ImageService Layer has : ArcGIS ImageService Layer                 {         "id": "<webmapOperationalLayer_imageService>",         "url": "<url1>,         "token": "<tokenString1>,         "title": "<title1>",         "opacity": <opacity1>,         "visibility": <true | false>,         "minScale": <minScale1>,         "maxScale": <maxScale1>,         "noData": <no data value>,         "format": "<jpgpng | png | png8 | png24 | jpg | bmp | gif | tiff>",         "interpolation": "<RSP_BilinearInterpolation | RSP_CubicConvolution |                                  RSP_Majority | RSP_NearestNeighbor>",         "compressionQuality": <compression quality>,         "bandIds": [           <bandId1>,           <bandId2>         ],         "mosaicRule": "<mosaic rule 1>",         "**renderingRule**": **-rendering rule 1-**        },      When the renderingRule is included, it receive a 404 error. When I do not set renderingRule, the task works fine. I use the same rendering rule as when I call exportImage on the image service itself. Has anybody seen an example of Export Web Map Task **handling a custom renderingRule** for an Image Service Layer? Thanks in advance for your help.


title: Dissolve multiple rows straight into one layer and add OIDs into table

text: I am currently running this script to direct multiple dissolves of a polygon layer               FC = "BLPUsCavPlace"     TABLE = "BLPUsCavPlaceContigBLPUs"     FIELDS = ["blpuid", "contigid"]     with arcpy.da.SearchCursor(TABLE, FIELDS) as cursor:         for row in cursor:             OID1 = str(int(row[0]))             OID2 = str(int(row[1]))             SQL = "OBJECTID = %s OR OBJECTID = %s"%(OID1, OID2)             arcpy.MakeFeatureLayer_management(FC, "TEMP_LYR", SQL)             OUTPUT_NAME = "Dissolve_%s_%s"%(OID1, OID2)             arcpy.Dissolve_management("TEMP_LYR", OUTPUT_NAME)             arcpy.Delete_management("TEMP_LYR")      However, this knocks out individual dissolved layers. I can always append them afterwards, but I would like to keep it clean and write them straight into another polygon layer called BLPUsCavPlaceDissolved. I would also like each row to contain the two OIDs that were dissolved, so there would be two extra columns called OID1 and OID2 that record them in the layer. Is this possible? Or will dissolve only ever create individual layers? Thanks


title: Projection 27700 in the OpenGeo Suite SDK GXP

text: I have some data stored as UK grid references, they display just slightly off alignment on the map after I followed the OpenGEo EXP tutorial using these settings (below) however when I change the EPSG to 27700 for UK grid reference, my map view loads with no layers, or tools visible. I know it displays ok using 4326 however I thought It is best to follow good practice and amend the projection to the correct EPSG.                 // map and layers     map: {         id: "mymap", // id needed to reference map in portalConfig above         title: "Map",         projection: "EPSG:4326",         center: [358, 54],         zoom: 6,                  layers: [{             source: "ol",             type: "OpenLayers.Layer.WMS",             args: ["World map", "http://vmap0.tiles.osgeo.org/wms/vmap0", {layers: 'basic'}],             group: "background"     },


title: ArcView: Trim polyline, that is shorter than specific dangle length (alternative to tool: Trim Line)?

text: I am looking for a way to trim dangles from a polyline shape, that are shorter than a specific length like the Trim Line Tool for ArcEditor and ArcInfo, but with ArcView-License. I already discover, that the ET geowizards-Toolbox has a similar tool called "Clean Dangling Nodes". Unfortunately this Tool also create new features where other features are closer than the specified dangling tolerance. You can see this in the following picture: Result (orange) of the tool Clean Dangling Nodes compared to the input data (blue) Where the tool works like I want to is marked with green circles. But there are several situations where the tool creates new features (what I don't want) that are marked with red circles. Any help is appreciated!


title: Why ModelBuilder cannot see Summary Statistics field?

text: I'm having some problems adapting a few models inherited by me from a colleague. Inside this model, we attempt to calculate Summary Statistics, but the model can't find the field that it needs to use. ModelBuilder refuses to show it to me, despite the fact that the Feature Class in this calculation possesses the field in question. Has anyone seen weird problems like this one? The original model was for 9.3 and I'm rebuilding it in 9.3.1. Any pitfalls I'm not aware?


title: Several questions about using spatial data?

text: 1. If you have spatial data for SQL server, can you use the spatial data (both geometry and geography) in another application for instance ESRI ArcGIS?   2. Is spatial data a standard that every GIS application can use it? For instance same spatial data can be used for SAGA-GIS, ESRI ArcGIS, Oracle Spatial, PostGIS without doing some convert or some edition in the data.   3. What system does spatial data use in SQL server? I heard that it use GPS? True? If yes, is there different version?


title: Modifying the start/end points of a LineString and maintaining connectivity in OpenLayers

text: I'm developing an application where the user can view and edit a street grid of a city. The streets are presented as features on a vector layer, and the user has to be able to modify the elements (both the attributes and the geometry - which is LineString - of the feature). The editing of a geometry is limited to dragging the end points of the LineString, at least for now. The trick is that the new end points must intersect with the current geometry of the feature, so that the user cannot drag a road to start from outside of the street grid. Also, modifying the geometry should affect the neighbouring geometries, if any. That is, if another road begins where one ends, both should have their start/end point changed. I don't suppose OpenLayers supports this out of the box so how should I begin to implement this?


title: How Does Directed and has_reverse_cost Arguments when set to [true, false] and [true, true] Affect Routing?

text: What will be the difference on the routing behaviour when I try to play with the directed and has_reverse_cost arguments on all the three routing algorithms? Example: the following SQL queries give me a different result each time. > SELECT rt.gid, ST_AsGeoJSON(rt.the_geom) AS geojson, at_2po_4pgr.gid FROM > at_2po_4pgr, (SELECT gid, the_geom FROM dijkstra_sp_delta_directed( > 'at_2po_4pgr', 66126, 9013, 0.1, **false** , **false** ) ) as rt WHERE > at_2po_4pgr.gid=rt.gid; > > SELECT rt.gid, ST_AsGeoJSON(rt.the_geom) AS geojson, at_2po_4pgr.gid FROM > at_2po_4pgr, (SELECT gid, the_geom FROM dijkstra_sp_delta_directed( > 'at_2po_4pgr', 66126, 9013, 0.1, **true** , **false** ) ) as rt WHERE > at_2po_4pgr.gid=rt.gid; > > SELECT rt.gid, ST_AsGeoJSON(rt.the_geom) AS geojson, at_2po_4pgr.gid FROM > at_2po_4pgr, (SELECT gid, the_geom FROM dijkstra_sp_delta_directed( > 'at_2po_4pgr', 66126, 9013, 0.1, **true** , **true** ) ) as rt WHERE > at_2po_4pgr.gid=rt.gid; I can understand if both are set to false but what if one is true and the other is false? Any explanation please? Thanks.


title: Finding the point closest to a point on a near feature

text: I have a dataset _osmMajor_ containing road data (from OSM; 670,780 features) and another dataset _points_ containing some points of interest (500,000 points). For each feature in my point layer, I want to find the coordinate of the closest point located on the road network. ![enter image description here](http://i.stack.imgur.com/33JLO.png) I have two options to do this: either ArcMap's near tool (Analysis tools/Proximity/Near), or Spatialite. I tried using ArcMap, but it seems a very long task to perform (I had to stop the process after 3 hours). For this reason, I was thinking about using Spatialite to do this, but it seems that the process is less straightforward than what I thought. My question is thus: **how should I write my Spatialite query to get what I want?** For Spatialite, I think that I should follow these steps:   1. For each point feature, calculate the distance between the point and the closest road    2. Around each point, create a buffer using the distance to the closest feature   3. Get the intersection between the buffer and the road data (it should be a point)   4. Get the coordinates (x,y) of the resulting point Trying to implement these 4 steps on a sample dataset with only 10 points, and I encountered several issues. For point (1.), I can find the overall smallest distance between the points and the roads, but I cannot have it for each point:               SELECT test.ID, min(Distance(test.Geometry, osmMajor.Geometry)) AS distance     FROM test, osmMajor     WHERE osmMajor.ROWID IN (     SELECT ROWID      FROM SpatialIndex     WHERE f_table_name = 'osmMajor'     AND search_frame = test.Geometry)     GROUP BY test.ID      I get a distance for only three points instead of the 10 that are in the _test_ table. I suspect it has something to do with the way the spatial index works, because the query only returns the results for points 2, 3 and 4, which are the closest to the road network. ![enter image description here](http://i.stack.imgur.com/Bp8xF.png) Can anyone give me advice on the best way to get the coordinate of the point on a feature which is closest to another point? Thank you!


title: Have there been any geospatial applications run in OpenStack by non-governmental entities?

text: This whitepaper about NASA's Nebula project, describes how cloud computing was used by NASA's Short-term Prediction Research and Transition Center to identify damage scars left by tornadoes. Since then, Nebula was open-sourced to create OpenStack. Have there been any geospatial applications run in OpenStack by non- governmental entities? ![enter image description here](http://i.stack.imgur.com/8QPlO.png) ![enter image description here](http://i.stack.imgur.com/vxNGG.png)


title: Can a mask layer in advanced drawing options be easily reversed?

text: I wonder if there is an solution to my crude workaround. I use often use the Advanced Drawing Options to mask items such as line underneath a contour label, etc. I am continually running into situations where I need to **reverse** this mask. Instead of masking a line features that intersects several large lakes for example I would like to mask the portion of that line feature that does not cross intersect those lakes: a reverse mask. The crude workaround is to create a large polygon feature and clip all the "lakes" out then use this feature as the Masking layer. It works but it is not easily maintained especially when the "masked" feature changes. _**Example:** I have a map of a province with shaded relief and 20 layers of various spatial information. I have a polygon of an area of interest (AOI) and I need to show all spatial information within the AOI but only basic information outside this AOI (i.e. shaded relief, hydrology). If I create a layer that contains the province polygon with the AOI clipped out, i can use it as a Mask in the Advanced drawing options. This creates redundant data and extra work. Is there a better workaround?_


title: Problem with mcc-lidar

text: I am running mcc-lidar on a win7 64bits with 16gb of ram. Nevertheless when I try to process a las file with 1, 800 000 points using -s 0.5 and -t 0.3 the program tells me: Error: insuficient memory; if applicable use a larger post spacing. I am able to process the data using -s 1.5 -t 0.3, so I figure the problem is not with the number of points. The program is installed at c:\program files (x86), which I understand is for 32 bits programs. Is there some specific instructions to installit as a 64 bit aplication in order to use my 16 gb of RAM? Thanks in advance Gerardo


title: render kml 10mb KML file cordinates data on OSM map throws stop running script error in OpenLayer.js file

text: I am trying to render the kml file having 18000 line string data on OSM map using openlayer HTTP strategy.   To render OSM map i am using GeoExt.MapPanel and ExtJS 3.4 viewport. The problem is I am facing the _stop running script error_ on page and in IE8 hang the PC and browser.


title: Update SDE Feature Class from Shapefile

text: I work for a department that get's a regular **City Council District Shapefile** showing the boundaries of each Councilors District. In our SDE we have a Council_District table with an exact duplicate of the previous **City Council District Shapefile**. In SDE we have added other fields, links to websites for each councilor, etc. We utilize Flex Maps to make these maps public. Everytime a new updated shapefile comes out, we copy over the Shapefile and write-over the SDE version with the new data. Then we have to re-create the added fields and their data. This requires that we shut down every web map using this file which can disrupt both in-house and external users and is not efficient at all. We are currently running version 9 servers, we will be updating to 10.1 in January. We will NOT be updating to 10.2 (yet!) due to data issues from our data providers. # Question(s) Is there any way to just update the SDE from the Shapefile to avoid thisor do we always have to shut down webmap services?


title: CartoDB Runtime Visualization on Local Instance

text: I'm trying to create a runtime visualization on a local instance of CartoDB. Here's the relevant code:               cartodb.createLayer(map, {                 user_name: 'testuser2',                 type:'cartodb',                 tiler_host:'', //??                 tiler_port: '', //??                 tiler_protocol: '', //??                 sublayers: [{                     sql: "SELECT * FROM metro_stop AS pop",                     cartocss: "#pop {marker-fill: #0fb16a; marker-width: 10;}"                 }]             }, {no_cdn: true}).addTo(map);      The problem is, I don't know the values of the three missing parameters: tiler_host, tiler_port, and tiler_protocol. Does anyone know how I would find those out? Edit to add: Here's the error I'm receiving when I set tiler_host to 'localhost.lan', tiler_port to '8181', and tiler_protocol to 'http': http://testuser2.cartodb.com:8181/api/v1/map?config=%7B%22version%22%3A%221.0.0%22%2C%22stat_tag%22%3A%22API%22%2C%22layers%22%3A%5B%7B%22type%22%3A%22cartodb%22%2C%22options%22%3A%7B%22sql%22%3A%22SELECT%20*%20FROM%20metro_stop%20AS%20pop%22%2C%22cartocss%22%3A%22%23pop%20%7Bmarker- fill%3A%20%230fb16a%3B%20marker- width%3A%2010%3B%7D%22%2C%22cartocss_version%22%3A%222.1.0%22%7D%7D%5D%7D&callback=jQuery17209372722392436117_1401785836614&_=1401785836710 Failed to load resource: Could not connect to the server. For whatever reason, this is pointing to CartoDB.com to try to get the info.


title: Problems running grass.mapcalc

text: I am using the GRASS Python Scripting Library on Ubuntu 14.04 with the packages available from the Universe repository. I have the following block of code:                       try:                 grass.mapcalc("$weightLine = $seg_rast * $cost",                             weightLine = "weightLine",                             seg_rast = "seg_rast",                             cost = "cost")             except:                 self.logger.debug("Exception running mapcalc: " +                                    str(sys.exc_info()))                  kwargs = {}             kwargs['stdout'] = grass.PIPE             kwargs['stderr'] = grass.PIPE             ps = grass.start_command("r.sum", rast="weightLine", **kwargs)             mesg = ps.communicate()             self.logger.info("r.sum: %s", str(mesg))      Which is logging the following:               [2014-06-03 09:57:35,532] DEBUG: Exception running mapcalc: (<type 'exceptions.SystemExit'>, SystemExit(1,), <traceback object at 0x7ffcbab1de18>)     [2014-06-03 09:57:35,538] INFO: r.sum: ('', 'ERROR: Raster map <weightLine> not found\n')      The employment of `sys.exit()` is not very helpful, even though it can be worked around. Without the `except:` block it exits my script and any other processes that may have invoked it. Note also that `except Exception, e:` does not catch `exceptions.SystemExit`. But worst of all, the `mapcalc` instruction is not succeeding. It is either failing - in which case an error message is lacking - or it is simply not being executed. How can I know the cause of this? **Update I** : I have tried to run `grass.mapcalc` in different ways:               grass.mapcalc("weightLine = seg_rast * cost")      And:               weightLine = "weightLine"     seg_rast = "seg_rast"     cost = "cost"                 grass.mapcalc("%s = %s * %s" % (weightLine, seg_rast, cost))      But the result is always the same: `exceptions.SystemExit` and the command is not executed.


title: GPservice cannot connect to sde

text: Recently we updated from ArcMap 10.1 to 10.2.2. ArcGIS Server is still on 10.1. We have numerous GP and Mapservices that were published on ArcMap 10.1 and are still running without problems. I also published a Map service from ArcMap 10.2, with data in an Enterprise GDB that is registered, without any problems. Now I have a python script that uses data from the same DB and works perfectly in ArcMap 10.2. I can publish it but when I run it, it returns false on the first               arcpy.Exists(sdeconnection\\FeatureClass).      According to this link I should however be able to publish form 10.2.2 to 10.1 without too many problems. Any ideas on what could be wrong? I guess it should be the data connection. But my Map service uses the same SDE, so my data should still be accessible. EDIT: In the original file, my workspace was defined as follows:               arcpy.env.workspace = "\\\\v057GIS\\Data\\Connection to v057gis.sde"      After publishing, this line changed suddenly to:               arcpy.env.workspace = arcpy.env.packageWorkspace      I can open the "C:\arcgisserver\directories\arcgissystem\arcgisinput\GPservice.GPServer\extracted\v101" file in python IDE, change the line and run it correctly. I'm however not able to save it. I read on this topic how to do it, but am a bit afraid of doing it since we only have production server. Is there a workaround for that? Maybe directly in packageWorkspace?


title: Most suitable open source GIS platform

text: I'm looking for an open source GIS server or desktop package that supports as many file formats as possible, is fast, runs on Linux or Mac and can be controlled via an API. Planning to do very basic editing such as intersects so I guess most GIS platforms will be able to do this. Does anyone have any recommendations as to a suitable package?


title: QGIS layers not adding properly

text: using QGIS 1.8.0 lisboa on windows 7 pc adding vector layers and they are not adding properly they don't show up in the layers list but they are there when I happen to merge the the layers altogether, I can see them listed, they are not hidden. it seems to happen when i try to add layers and have control layer render order on, because when i turn it off the layers adds fine, but i cannot delete the previous layer that didn't add correctly, i can't see them to select and remove them ![enter image description here](http://i.imgur.com/nuSANEH.jpg)


title: How to create a search box and SQL statement for searching an attribute field?

text: How can I add a HTML search bar that the user can use to search an attribute field for a polygon layer, and then tie that input into the SQL statement? So far I am able to hard code the query of sending a SQL statement and returning a answer using cartodb. My problem is figuring out how to take the input from the search box and automatically inserting the input into the sql statement. I am using Cartodb api.                var sql = new cartodb.SQL({ user: 'jmann111' });        sql.getBounds("select * from stations where station_nm like ='%" + user_input +"'").done(function(bounds) {       map.fitBounds(bounds)       });


title: How to uniquely identify instances of same Compatible Unit (CU) type in work location?

text: I'm attempting to synchronize a design with edits made in our customer's external Work Management System. To do this, I need to be able to distinguish between two compatible units (CUs) of the same type (e.g. 40ft wooden poles) in the same work location so that I can update the extended data of the correct pole. There seems to be a unique identifier for each instance of a CU in the PROCESS.MM_WMS_COMPATIBLE_UNIT table called ID. Also, there seems to be changes in the schema between the version we're using currently (9.3.1) and version 10.1. That said, I am hoping the API provides access to these, or some equivalent, unique identifiers. Initially, I thought IMMCompatibleUnit2.DesignerId and IMMGisUnit.GisUnitId might be what I was looking for. From the Developer help: DesignerID Property - "Gets the DesignerID associated with the CU." I couldn't find anything in the help about GisUnitId. Also, I couldn't find where a description of the Design ID was given, so I'm not even sure what it is. When I attached the debugger and examined the DesignerIds and GisUnitIds, they were always 0, which just added to my confusion. Does anyone know how to uniquely identify CUs in a DesignTree? Thanks for your help!


title: Downloading an unzipped File Geodatabase from USGS

text: I am attempting to download a file geodatabase from the pubs.usgs.gov site. The file db is DDS69ff.gdb The browser sees the file as a directory so I am having trouble downloading the complete file as a geodatabase accessible in Arcgis 10. I've tried to retrieve the file in wget and I've had no luck.


title: Issues ST_MakePoint in PostGIS with Geoserver

text: I am trying to setup a GIS system to use PostGIS and Geoserver. I have loaded all of the data into a table. The data that is received has long/lat coordinates in Degrees, Minutes, Seconds, Direction, each of which is a separate column. What I did is created a column (location) of type geography and I'm trying to give proper info for each record. The command that I am using is this (of course this is only in the Western Hemisphere):               UPDATE *tablename* SET location = ST_SetSRID(ST_MakePoint((lat_degrees+(lat_minutes+lat_seconds/60)/60), (-1*(long_degrees+(long_minutes+long_seconds/60)/60))), 4326);      My DMS -> DD conversions should be correct, and I've verified that by creating separate columns to test each individually. But the problem is, my geoserver doesn't detect the proper bounds. When I manually set the bounds (I set them to -180,-180,180,180) and check my data in Google Earth, it only zooms to 0,0 and I can find some horribly warped data, in Antarctica. One thing to note... not all of my records include location data, could that be an issue? Or is ST_MakePoint simply not creating the right location for my data? Thanks for the help!


title: Formating geojson file with leaflet ajax-plugin

text:
First of all I apologize for my basic english. I have to say, too, that I'm really a newby in coding. I'm trying to create an interactive map using leaflet (you can see it here). I was able to show my geojson file using ajax- plugin. Now I want to change the style of my geojson file, but I really don't know how I can do that. Can someone explain me? This is my code:               <!DOCTYPE html>               <meta charset=utf-8 />     <title>Geojson data</title>     <meta name='viewport' content='initial-scale=1,maximum-scale=1,user-scalable=no' />                                  <link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.7.2/leaflet.css" />

                     Thanks


title: How do I parse the file format for attributes in Raster Product Format (RPF) maps and imagery?

text: I'm trying to implement an open-source NITF reader (using only java, no external libraries such as GDAL or NITRO), and its mostly going OK. However I'm having trouble understanding the RPF extensions, in particular the attributes in the RPFDES TRE. I know that they are documented in MIL-STD-2411-2, and I've read that and the main references (MIL-STD-2411 and MIL-STD-2411-1). I still don't understand how to parse it out though. As a concrete example, http://download.osgeo.org/gdal/data/nitf/bugs/bug3337.ntf (which is NITF 2.0) has the RPFDES TRE looking like:               01efe0 20 02 00 20 02 00 20 02 00 20 02 44 45 52 65 67  > .. .. .. .DEReg<     01eff0 69 73 74 65 72 65 64 20 45 78 74 65 6e 73 69 6f  >istered Extensio<     01f000 6e 73 20 20 20 20 30 31 55 20 20 20 20 20 20 20  >ns    01U       <     01f010 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20  >                <     *     01f0a0 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 55  >               U<     01f0b0 44 49 44 20 20 30 30 31 30 30 30 30 52 50 46 44  >DID  0010000RPFD<     01f0c0 45 53 30 30 36 32 31 00 21 00 02 00 00 00 00 00  >ES00621.!.......<     01f0d0 08 00 01 01 00 00 00 01 08 00 02 01 00 00 00 01  >................<     01f0e0 10 00 03 01 02 00 00 01 18 00 04 01 02 00 00 01  >................<     01f0f0 20 00 04 02 02 00 00 01 2a 00 04 03 02 00 00 01  > .......*.......<     01f100 32 00 04 04 02 00 00 01 36 00 05 01 02 00 00 01  >2.......6.......<     01f110 3d 00 05 02 02 00 00 01 3f 00 05 03 02 00 00 01  >=.......?.......<     01f120 43 00 05 04 02 00 00 01 47 00 06 01 01 00 00 01  >C.......G.......<     01f130 4b 00 07 01 01 00 00 01 4f 00 08 02 01 00 00 01  >K.......O.......<     01f140 53 00 09 02 01 00 00 01 55 00 0a 02 01 00 00 01  >S.......U.......<     01f150 57 00 0b 02 01 00 00 01 59 00 0c 01 02 00 00 01  >W.......Y.......<     01f160 5b 00 0e 01 02 00 00 01 5e 00 10 01 02 00 00 01  >[.......^.......<     01f170 60 00 10 02 02 00 00 01 64 00 11 01 02 00 00 01  >`.......d.......<     01f180 66 00 11 02 02 00 00 01 6a 00 12 01 02 00 00 01  >f.......j.......<     01f190 6c 00 12 02 02 00 00 01 70 00 13 01 02 00 00 01  >l.......p.......<     01f1a0 72 00 13 02 02 00 00 01 76 00 14 01 02 00 00 01  >r.......v.......<     01f1b0 78 00 14 02 02 00 00 01 80 00 14 03 02 00 00 01  >x...............<     01f1c0 82 00 14 04 02 00 00 01 8a 00 19 01 02 00 00 01  >................<     01f1d0 92 00 19 02 02 00 00 01 94 31 39 38 34 30 33 31  >.........1984031<     01f1e0 30 32 30 30 31 31 30 31 32 31 39 38 34 30 33 31  >0200110121984031<     01f1f0 30 47 4e 43 20 20 20 20 20 20 20 30 31 20 20 20  >0GNC       01   <     01f200 20 20 20 55 4e 44 20 20 20 20 20 20 20 37 54 43  >   UND       7TC<     01f210 00 00 00 00 40 40 00 00 3f 7f e5 c9 4d 53 4c 20  >....@@..?...MSL <     01f220 57 47 45 00 00 00 00 00 00 00 00 00 49 4e 54 75  >WGE.........INTu<     01f230 64 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00  >d...............<     01f240 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00  >................<     *     01f260 00 00 00 00 00 00 00 00 00 00 08 00 52 00 10 00  >............R...<     01f270 05 40 52 f6 db ba 46 e8 05 c0 60 e0 00 00 00 00  >.@R...F...`.....<     01f280 00 40 52 f6 db ba 46 e8 05 c0 60 e0 00 00 00 00  >.@R...F...`.....<     01f290 00 40 52 f6 db ba 46 e8 05 40 46 80 00 00 00 00  >.@R...F..@F.....<     *     01f2b0 00 40 52 f6 db ba 46 e8 05 c0 60 e0 00 00 00 00  >.@R...F...`.....<     01f2c0 00 00 05 40 56 80 00 00 00 00 00 c0 66 80 00 00  >...@V.......f...<     01f2d0 00 00 00 40 56 80 00 00 00 00 00 40 66 80 00 00  >...@V......@f...<     01f2e0 00 00 00 40 50 ff f6 ab 93 ae fe 40 66 80 00 00  >...@P......@f...<     01f2f0 00 00 00 40 50 ff f6 ab 93 ae fe c0 66 80 00 00  >...@P.......f...<     01f300 00 00 00 40 56 80 00 00 00 00 00 c0 66 80 00 00  >...@V.......f...<     01f310 00 00 00 00 00 00 08 00 01 00 19 30 30 30 30 30  >...........00000<     01f320 30 32 33 2e 47 4e 39 30 30 30 30 30 30 31 33 2e  >023.GN900000013.<     01f330 47 4e 39 33                                      >GN93<     01f334      I'm OK with the parts up and including the RPFDES tag name and length. After 01f0c7, I'm lost. From using GDAL, I can see that there are some dates starting at 0x01f1d9, and some other things that are probably the CRS and elevation references, but I'm mystified as to how I'm meant to extract them. Can someone work through the data structure parsing to get to the first date, or at least provide a step- through of how I'm meant to use MIL-STD-2411-2 to work it out? _For those wondering why it has to be only java - its a deployment thing for OSGI. For those wondering where the source code is, seehttps://github.com/bradh/nitf-file-reader._


title: How to Overlay the WMS Layer on Google map

text: I want to overlay WMS layer on Google layer of OpenLayers with my map server. but My WMS not Overlayed on it. my code is here:                function init()         {             var map;                           map = new OpenLayers.Map('map', {                                    projection: new OpenLayers.Projection("EPSG:4326"),                     displayProjection: new OpenLayers.Projection("EPSG:900913")                 });                 map.addControl(new OpenLayers.Control.LayerSwitcher());                 var gmap = new OpenLayers.Layer.Google(                          "Google Streets",                     { numZoomLevels: 20},                      { sphericalMercator: true }                           );                 var ghyb = new OpenLayers.Layer.Google(                     "Google Hybrid",                     { type: google.maps.MapTypeId.HYBRID, numZoomLevels: 20 }                 );                      var gsat = new OpenLayers.Layer.Google(                 "Google Satellite",                 { type: google.maps.MapTypeId.SATELLITE, numZoomLevels: 22 }                  );             map.addLayers([gmap,ghyb, gsat]);                      var wmsm = new OpenLayers.Layer.WMS("wms", "http://188.136.237.210:8090/cgi-bin/mapserv.exe?map=C:/ms4w/Apache/htdocs/maps/nikgis.map&",           {                    layers: 'GIS_Division',               format: 'image/png',               transparent: 'true',               visibility: true,                },                     {                         "reproject": "true"                          }                    );                wmsm.isBaseLayer = false;                 map.addLayer(wmsm);                 map.setCenter(new OpenLayers.LonLat(0, 0), 2);                      //this line is for become SRS:EPSG=4326                gmap.projection = map.projection;                             }      and the scripts that imported is:                <script src="http://maps.googleapis.com/maps/api/js?key=AIzaSyDY0kkJiTPVd2U7aTOAwhc9ySH6oHxOIYM&sensor=false"></script>      after running the project this URL is sent: mapserverURL&LAYERS=GIS_Division&FORMAT=image%2Fpng&TRANSPARENT=true&REPROJECT=true&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&SRS=EPSG%3A4326&BBOX=-90,-45,5009287.085,5009332.085&WIDTH=256&HEIGHT=256 (this URL should be sent with EPSG:4326.) the status code is OK. without google layer , the WMS Layer is shown nice and the URL that sent is: mapserverURL&LAYERS=GIS_Division&FORMAT=image%2Fpng&TRANSPARENT=true&REPROJECT=true&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&SRS=EPSG%3A4326&BBOX=39.375,33.75,45,39.375&WIDTH=256&HEIGHT=256 I think the problem is for BBox. is there any way to show My WMS Layer.


title: Postgresql Postgis Query is taking long time

text: I have a table that included only 472 rows and geometry of table is Multipolygon. I am sending query to databse via **PgAdmin III**               select * from mytable // this query takes 99413 ms      this is vey long time. what do you think about long queryies on postgresql? Server resurce, disc, cpu etc.


title: Is it possible to import arcpy module to PyPy interpreter?

text: I'm trying to use PyPy for my code because of speed. But my code needs arcpy module to run and I cant find out how to install or import this module to arcpy. I'm using PyPy 2.3, Python 2.7 and Win 7. Any suggestions?


title: Adding attributes to a hexagonal grid generated with MMQGIS

text: I generated a hexagonal grid with MMQGIS, and I need to add some attributes to each hexagon. Is there a convenient way to do this or I have to work with the dbf file? In that case, what is the meanign of the attributes stored in the generated dbf file?               left,N,24,16     top,N,24,16     right,N,24,16     bottom,N,24,16      Thanks for help!


title: OpenLayers and Geoserver. Using multiple Servers for tiles

text: OpenLayers allows the wms layer to specify up to 4 servers to permit the browser to connect to a server multiple times Has anyone had any experience in mapping those for servers to four Physical geoserver Instances. The Geoservers (v2.2.5) are sharing a common data directory, different cache directories for the Integrated GeoWebCache Would it be reasonable to expect it to work ok? I believe the functionality is to point four different dns entries at the SAME Server, not Four DNS Entries to Four DIFERENT Servers sharing a common database. The reason is we are trying to speed up tile delivery and the performance on the scenario of pointing 4 dns entries to one server is not what we would expect. We have done all the tuning on the Java side that we can find information on and want to see what impact this makes, if it is at all even possible. Platform is a public facing windows IIS7.7 and geoserver 2.2.5 on it own machine on a private network We are also using a proxy that distributes the tile request to the Geoserver Instance on a private network. If we proceed, we could also use 4 proxy servers as the end points for the public Geoserver Request. Comments and suggestions welcome.


title: Importing dxf as a vector layer or using dxf2shp function (Differences?)

text: I was trying to import a DXF file in QGIS. First tried using dxf2shp function but it didn´t work. I browse the dxf file and i put a name to the output and nothing was created in directory i chose. I couldn´t convert a DFX file in order to import it with this method, i dont know what i did wrong. I tried in other way with "add vector layer" selecting (Source: file,code:UTF-8; and Browse my dxf) and it works. what is the difference between the 2 ways, is there a disadvantage with the vector layer option? Do you know why first method didn´t work for me? Regards


title: Writing shapefile metadata with geotools

text: Using GeoTools, I'm able to write features to a shapefile by essentially passing a `SimpleFeatureCollection` to a `ShapefileDataStore`. The feature collection provides both a schema and features including their attribute values. In addition to having the row-per-feature attribute table in the shapefile, I'd like to include file-level metadata - either structured fields or at least free-text title/description fields. Is there support for doing this in GeoTools? Or do I need to serialise my own `*.shp.xml` file?


title: Informations about API in GDAL for geocoding

text: I write to ask you if did you ever used the new API in GDAL for geocoding. (Link for description: http://www.gdal.org/ogr/ogr__geocoding_8h.html#a2b0382ceb7e3935d7fe1da2adada564e ). Do you know what is the method with gdal can associate a location with specific coordinates? I ask you because I don't understand what does it means : "It runs a geocoding request.If the result is not found in cache, a GET request will be sent to resolve the query". I found a simple code to require a location coordinates:               /* Create a session with default options */          OGRGeocodingSessionH hSession = OGRGeocodeCreateSession(NULL);          /* Now query the service */          OGRLayerH hLayer = OGRGeocode(hSession, "Paris", NULL, NULL);          /* Use the result layer (the first feature is generally the most relevant) */          /* Cleanup */          OGRGeocodeFreeResult(hLayer);          OGRGeocodeDestroySession(hSession);      Do you confirm this is the code to require a location coordinates? And if I want to find coordinates of a street adress?


title: QGIS Extend Toolbar (GUI) Button Functionality via Python Plugin?

text: Is it possible to piggyback onto the QGIS split method via a python plugin? Are toolbar buttons accessible via the python API? This looks like the class that I need to access. Access ideas? Proposed Workflow to Give an Implementation Example:   1. Initialize editing   2. Split the line as per How to Split Features in QGIS?   3. Launch a custom GUI or otherwise call a set of functions to start to build a list that I push to the DB table on saving edits.


title: Cached Featured in WFS - can I change default

text: QGIS currently sets the default to cached download for adding WFS - is it possible anywhere to set the default to non cached mapping (I know you can untick the box, but its for a larger user group)?


title: Problem with Create Fishnet on a Landsat footprint

text: I have a Landsat footprint polygon shapefile that I simply want to create a 6 x 6 fishnet on. I thought this would be simple, but every time I run the tool I get a fishnet output that extends longer than the original footprint and is shorter in height. I think the problem may stem from the fact that the footprint is slanted and the y-axis may be a different angle than the x-axis. I used Feature Vertices to Points and then Add XY Coordinates to get the exact coordinates of the footprint polygon's corners. Then within the Fishnet tool, I used the footprint's extent as the 'Template Extent', the bottom left corner as the 'Origin Coordinate', the top left corner as the 'Y-Axis Coordinate', and the top right corner as the 'Opposite corner of Fishnet.' Coordinate system is WGS_1984_UTM_Zone_18S. Any ideas for how I can fix this and generate a properly placed fishnet? (I can't post images yet as this is my first post on here). Thank you. EDIT: Can post images now. Here's the strange looking fishnet output: ![Fishnet](http://i.stack.imgur.com/7NDh1.jpg) And here are the parameters I passed into the fishnet tool, along with the proper corner coordinates in the top left of the image (Point_X and Point_Y): ![Fishnet tool screenshot](http://i.stack.imgur.com/89Y43.jpg)


title: Running OpenLayers in .NET 4.5 breaks

text: So, I'm writing an web mapping application using OpenLayers in WebForms in .NET 4.5. It was working wonderfully until my browser upgraded to IE9. Now the JavaScript fails. After researching I found that the problem seems to be that IE9 changed how it handles certain kinds of JavaScript. So the solution, until the community can rewrite the JavaScript libraries, seems to be to have the browser run in compatibility mode. Unfortunately, I can't seem to get compatibility mode to work correctly. Currently my application has one page and this is the header:               <%@ Page Language="C#" AutoEventWireup="true" CodeBehind="default.aspx.cs" Inherits="AgPlaceHolder._default" %>          <!DOCTYPE html>          <html xmlns="http://www.w3.org/1999/xhtml">     <head runat="server">     <meta http-equiv="X-UA-Compatible" content="IE=EmulateIE8" />     <title>Agrilogic Place Holder</title>     <script src="OpenLayers.js"></script>     <script src="AgInit.js"></script>     <script src="AgLabels.js"></script>          <style>         html, body         {             width: 100%;             height: 100%;             margin: 0;             padding: 0;         }     </style>     </head>      It didn't work. So I found this alternative to add to the web.config file:               <system.webServer>       <httpProtocol>         <customHeaders>           <add name="X-UA-Compatible" value="IE=EmulateIE8"/>         </customHeaders>       </httpProtocol>     </system.webServer>      This also didn't work. I've tried some variations I've seen used around, but none of them seem to be working. Does anyone see what I might have done wrong? Or has another solution to this problem? Thanks. On Edit: I'm using OpenLayers 2.12. The error is:               Unhandled exception at line 670, column 407 in http://localhost:57639/OpenLayers.js          0x800a138f - Microsoft JScript runtime error: Unable to get value of the property 'resources': object is null or undefined      The code for setting up the map is:               map = new OpenLayers.Map({         div: "map",         layers: layerArray,         controls: [             new OpenLayers.Control.LayerSwitcher({ roundedCorner: true }),             new OpenLayers.Control.PanZoomBar(),             new OpenLayers.Control.Navigation(),             new OpenLayers.Control.ScaleLine({ geodesic: true }),             new OpenLayers.Control.KeyboardDefaults(),             new OpenLayers.Control.MousePosition({                 displayProjection: new OpenLayers.Projection("EPSG:4326"),                 numDigits: 5,                 emptyString: '',                 prefix: "x: ",                 separator: " y: "             }),             //new OpenLayers.Control.EditingToolbar(editableLayer),             //new OpenLayers.Control.SelectFeature(editableLayer, { toggle: true, clickout: true }),             //new OpenLayers.Control.DragFeature(editableLayer),             overviewMapControl         ]     });      The java script works fine, except for another issue I have posted a question on over in GISStackExchange, but that's not preventing the maps from showing. The problem there is I can't drag the label-points the program creates. But now with the upgrade of the browser I can't even display the maps anymore. The page just fails. The newly built page has the same errors. When I remove Bing:               0x800a138f - Microsoft JScript runtime error: Unable to get value of the property 'wrapDateLine': object is null or undefined      When I include Bing:               0x800a138f - Microsoft JScript runtime error: Unable to get value of the property 'resources': object is null or undefined      I am at the mercy of your expertise.


title: Do FeatureLayers reproject their features to adapt them to the map projection?

text: I'm working with the ArcGIS Flex API. For example, ArcGISDynamicMapLayers do reproject their data in order to make it in the same projection as the map. Do FeatureLayers have the same behaviour?


title: Retrieve polygon smallest MBR

text: I have a postgis table containing polygons. I want to transform my polygons to 4 points polygons : retrieve the smallest MBR. Here comes an illustration of what I want : ![alt text](http://i.stack.imgur.com/H6qTG.png) My original polygon is in pink color. I guess, it will be quite easy to get the MBR (with ST_Envelope and BBox function) But is there an easy way to get what I want? Any ideas? Many thanks


title: How to interpolate values sampled at regular time intervals?

text: I want to numerically integrate a set of discrete values that keep arriving at constant time intervals. Which interpolation method do I choose? how to reduce the ringing artifacts if I use cubic spline interpolation?


title: How to PGRoute from OpenLayers via PostGIS/PostgreSQL

text: I need to route a shortest distance for a map which I have loaded from a PostgreSQL by means of JSON file. I can project the map in my browser and I can identify the linestrings. For a particular linestring I wanted to check whether distance is calculated or not. Soon after, I would like to check among the different linestrings(nearest node). All I want is after marking start point and end point that route would calculate a distance. Here I could get vertices of a point when clicked. Data I am having is a JSON file and PostgreSQL Server/DB, JScript, OpenLayers.. And moreover, I would like to display attributes over the linestring from the JSON file properties. Is this possible. Please help me out. I got stuck for the next step. Thank you


title: Visualization libraries for rendering 3D GIS world on iOS/Android

text: Are there any good 3D libraries/frameworks for rendering 3D worlds in geospatial coordinates on mobile platforms? Use case: I have some known 3D geospatial features (eg. buildings) with their corresponding lat,long,elevation values as well as an underlying DEM. I would like to render such a "scene" from a given camera viewpoint on either an iOS or an Android device. I could use opengl-es but it would take some effort to convert the geospatial coordinates to cartesian coordinates. Is there a "geospatial-aware" soft of 3D library for mobile platforms that I could use? Something like the AGI insight3D libraries but for mobile. Thanks!


title: How to display local file's geometries as rubber bands in QGIS

text: I'd like to get geometries from a local file (e.g. GML) displayed as rubber bands (QgsRubberBands) in QGIS. However, I wouldn't like to load the layer into the map because in the real scenario there could be a large number of files to get geometries from. I've been attempting it with the following code (you can copy it directly to the QGIS Python Console), but the problem is that only the last feature is displayed:               from qgis.core import QgsFeature, QgsVectorLayer     from qgis.gui import QgsRubberBand     def loadGMLAsRubberBands(fileName):         featList = []            rb = QgsRubberBand(qgis.utils.iface.mapCanvas(), True)         vlayer = QgsVectorLayer(fileName, "layer", "ogr")         allAttrs = vlayer.dataProvider().attributeIndexes()         vlayer.select(allAttrs)                 f = QgsFeature()         while vlayer.dataProvider().nextFeature(f):             featList.append(f)                       rb.addGeometry(featList[-1].geometry(), None)             loadGMLAsRubberBands("/tmp/gml/sample_three_features.gml")      A sample GML file can be found here. I've already attempted several things with no success, namely to store layers (vlayer), features, geometries and/or rubber bands in a list (as suggested in some QGIS mailing-list threads). I'm using QGIS 1.8.0-Lisboa (rev. a1255fc), Python v.2.7.2, GDAL/OGR v.1.7.3 on Ubuntu/Linux. Any help would be appreciated.


title: Eclipse can't find Driver.class from postgresql-9.2-1002.jdbc4.jar file nor AbstractJdbc4Connection class

text: I am using the jar file in my configuration path after following the example from postgresql/jdbc tutorial. But after I added the jar file and configured it for my reference library when I run it I am getting an error that says it can't find the source attachment for that jar file despite my retrieving that jar file from the jdbc postgres website. I have tried a couple of different configuration methods, but it the jar file and its classes have mostly already been added with the exception of a couple classes like Driver.class. Should I use a different version or is it a configuration issue? I don't believe its something in my code yet because no code error other than not finding the driver class would cause that. I am connecting to pgadminIII to update some rows.


title: How to geocode Guam addresses with their latitude/longuitude?

text: I am attempting to get lat/longs for a table of 25,000 addresses for Guam. Guam's addressing style doesn't follow the typical address number range. I've exhausted all of my geocoding resources. ESRI, Google, Bing, Yahoo, HERE, and various educational facilities. Google, thus far, has been most promising, but I am only get street centerline values returned. I was curious whether anyone out there had attempted this before, and if anyone had any resources/ideas about how I can go about this. Thank you! Using ArcGIS 10.2 Advanced Sample Addresses:               103 Ciriaco Loop, 96910-2070     143 Taluba Rd, 96910-2078     176 Chorito Blvd, 96910-1600     310 Ladera Ln, 96913-5868     141 Jolene Leon Guerrero St, 96913-1324