title: XY event source layer problem (converts into feature class after running python script ) in ArcMap 10.1

text: I am using ESRI 10.1 suit (ArcGIS Desktop & ArcGIS Server, No ArcSDE) at my end and after running python script n number of times event layer converts into feature class. Following steps I am following : -I am using query layer in ArcMap & fetching records from MS sql server, table type is plane table.  -Created XY event source layer from query layer and displaying points on the map. -layer is classified using unique value symbology. -My query layer is updating five times in day and using python script (example2) to automate the unique value symbology (using windows scheduler) five times in day. -Overwriting mxd on ArcGIS Server 10.1 using ArcPy in specific time interval. -After 10 times my XY event source is automatically converts to feature class then I am not able to use ArcPy script to update the legend at my end So is there any bug related to XY event source layer or am I doing something wrong?? Researching on Google and ArcGIS web help not getting any reason/limitation for that. Any help will be great !! thanks in advance !! **EDIT :** Legend update python script Map service overwrite script ![Screen-shot](http://i.stack.imgur.com/d4YVy.png)


title: Is there a way to pass arguments to qbrowser?

text: is there a way to pass arguments to qbrowser? I was thinking that maybe one could use the QGIS browser as a simple shapefile viewer, but it would be nice if one could specify the folder/file it should open when launched I am on Windows 7 64-bit, and I have QGIS 2.2 installed through OSGeo4W 64-bit.


title: Changing default parser in ArcGIS Field Calculator to Python?

text: Is there a way to change the default parser in the Field Calculator to Python? So annoying to always change now that I'm trying to move away from using VBScript... I'm using ArcGIS 10.1 for Desktop.


title: Density slice/Colour slice in QGIS?

text: Can you do Raster Color Slices (like in ENVI) in QGIS? I.e. Where the histogram can be divided up into discrete units to represent different classes based on the values of the pixels.


title: Using OpenStreetMap for simulation of commuter routes

text: I'm seeking some advice on how to build a simulation based on OpenStreetMap Data for my masterthesis. My knowledge of spatial data is very limited, since I come from another field of science, namely information systems. I have experience with object oriented programming, as well as mySQL as relational database. The research goal is to evaluate if applications which inform the user of the current gasoline price actually have the potential to save money. So the idea is to simulate the routes people travel between their home and work by car and when the car needs gasoline to find the nearest gasoline station and refill. As Database I will use data from the www.zensus2011.de which can tell me from which administrative district people leave and to which administrative district they travel. Within OSM I have the following Tags that can identify the administrative regions/districts which I need   * de:amtlicher_gemeindeschluessel => integer   * de:regionalschluessel => integer   * admin_level => 6   * boundary => administrative So my idea is to use this data to generate random routes with a certain length based on other theses which describe how far people travel between home and work. One of the main problems which I could identify from reading the OSM wiki and other online posts is that the OSM data per se is not routable and needs to be transformed for that purpose. I for that I found two programs osm2pgrouting and OSM2PO. With the latter one I managed to import the map of Germany into a PostgreSQL Database. I also imported the OSM data using osm2pgsql with the "default.style" settings into another PostgreSQL Database. Now I'm a little overwhelmed by all the data and the way how to use this data. But let's continue with the idea. What I need to understand now is how the map data with all the information and the routing data can be connected to generate the routes. What I thought of was: I use the administrative boundary to get a list of streets within this boundary, then randomly choose one as a starting point and do the same for the destination, calculate the route and save the route to a database for further use in the simulation. Still I need to find out a way to determine the start and destination points better, but for the beginning and understanding it would be fine with the above. Some final words on the implementation of the simulation. Currently my plan is to use python and code the whole simulation, which I hope is feasible. I found two related works which I'm currently investigating namely MATSim and www.roaf.de. Also I'm still browsing other questions here which might help me understand the topic. To summarize my questions:   * Is the way I choose feasible?    * How would SQL statements look like to connect the two Data sets?    * Is it possible to save generated routes to a database? Thank you for reading and I hope you guys can give me some advice on where to look at.


title: Displaying click results in a side panel

text: I am trying to do something similar to the interactive choropleth map (http://leafletjs.com/examples/choropleth.html) on leaflet's website. In my case I want to use a click instead of a hover and want to show the results in a side panel (i.e. separate div). Can someone point me to an example or suggest how I can do this? UPDATE: Based on Bronco's reply: That part I understand. It is this snippet that I need to change but I don't know how:               var info = L.control();          info.onAdd = function (map) {     this._div = L.DomUtil.create('div', 'info'); // create a div with a class "info"     this.update();     return this._div;     };          // method that we will use to update the control based on feature properties passed     info.update = function (props) {     this._div.innerHTML = 'US Population Density' + (props ? '' + props.name + '     ' + props.density + ' people / mi2'     : 'Hover over a state');     };     info.addTo(map);


title: Extracting shape point coordinates from polylines in MapInfo

text: My question follows from Coordinate extraction in a different projection in MapInfo . I now have WGS84 begin and end coordinates for the polylines(which are in british grid coordinates) How do I now extract the begin and end coordinates of each shape point associated with each polyline in WGS84?


title: Error 000865 when using arcpy.ASCIIToRaster_conversion

text: I don't know what I'm doing wrong. I always get the error 000865 when using arcpy.ASCIIToRaster_conversion. I'm new and do not have much Python experience right now – so sorry for that.               >>> import arcpy     import os     from arcpy import env     env.workspace = r"D:\Julia_T\projekt\erste_schritte\dritter_tag_python\testdaten"          #folder to list     ascFileList = arcpy.ListFiles("*.asc")          for ascFile in ascFileList:     # geoprocessing steps     ascFileName = os.path.splitext(ascFile)[0]     # define the output file     rastFile = ascFileName + "_output.img"     # run the tool     arcpy.ASCIIToRaster_conversion(ascFile,rastFile, "INTEGER")      > Runtime error Traceback (most recent call last): File "", line 15, in File > "c:\programme\arcgis\desktop10.2\arcpy\arcpy\conversion.py", line 2193, in > ASCIIToRaster raise e ExecuteError: ERROR 000865: Eingabe-ASCII-Raster- > Datei: 2001002_cloud_free.asc ist nicht vorhanden. Do I have to use the CatalogPath? I'm working with ArcMap 10.2.1. Thanks


title: Incorrect display of national characters from Spatialite in QGIS, Windows

text: I used this tutorial by Underdark (thanks much :-) ) to convert my area of interest from OSM file (whole Czech Republic downloaded from Geofabrik.de server). Used OSGEO4W shell (saying "GDAL 1.11.0, released 2014/04/16") and following command: `C:\GIS\CZ_OSM_data_2014-07-22T20-21-02>ogr2ogr -f "SQLite" -dsco SPATIALITE=YES -spat 12.50004 50.46600 12.93001 50.15486 nejdek_OSM_2014_07.db czech-republic-l atest.osm.pbf` I loaded the exported DB in QGIS and found out that the places labels are not displayed correctly in QGIS 2.4.0-Chugiak (GDAL/OGR 1.11.0, SpatiaLite 4.1.1) running on Windows 7 64bit while the same project is ok when opened in 2.4.0-Chugiak (GDAL/OGR 1.10.0, SpatiaLite 4.1.1) running on 32bit Kubuntu linux on my laptop. Tried to find and encoding settings or something similar but no success... Any advice or idea what am I doing wrong?


title: Using update cursor with list of shapefiles

text: I'm still pretty new to using cursors. What I'd like to do is create a list of shapefiles. Within those shapefiles I would like to populate a field "NAME" with the shapefile name less the .shp extension. I am using 10.1 and my code is below:               bufferCalc = arcpy.ListFeatureClasses("*.shp")     for calc in bufferCalc:     NOTE: This space contains calculations for other fields     with arcpy.da.UpdateCursor(calc, ("PROJNAME")) as cursor:         for shp in arcpy.ListFeatureClasses():             cursor.updateRow(os.path.splitext(shp))      The other calculations I need run smoothly however, the error I get when the cursor is ran states "StopIteration: iteration not started". This has got to be something simple I'm doing wrong but can't seem to figure it out. Thanks for any clues you can give me.


title: Publicly available ArcIMS server?

text: Just..don't...ask. It's a dark tale. :-P Does anyone know of a publicly / externally available ESRI ArcIMS server / service / services with the full US base map? Doesn't have to be free, but that's always good. I've seen some city-specific ones, but I need a USA one with street-level detail. Thank you.


title: Attribute information not showing in Silverlight viewer popup

text: I have feature service which I am accessing through a Silverlight viewer. The popups contain all the available fields (about 40), most of which have coded value domains. Looking at the popup, these fields appear blank. If I click on the `Edit Values` button, the fields' existing attribute information appears (which I can see in ArcMap). This might be an issue with version 3.1 of the viewer; similar issues are mentioned on the ESRI forums. I have not had this problem before, and I could see the values in the viewer three days ago. The only thing that's different is I had to republish the service after making cosmetic changes to the map. The data, field aliases etc all stayed the same.


title: Adding attributes to line feature according to adjacent polygon

text: I'm trying to add attributes to a field in a polylines attribute table according to directly adjacent polygons.... ![Copy attributes adjacent polys](http://i.stack.imgur.com/0Ke7W.jpg) In the pic, the purple line has an empty field labelled 'FG' and I'd like to fill it with whatever is in the 'FVG' field in the adjacent polygon that shares the boundary ie. 'Rf'.. 'Mit'...'Con'. The purple polyline is split into segments at each polygon boundary. This should be quite simple but can't for the life of me figure it out. Is it spatial join I'm after? Haven't had much experience with that


title: Building java based desktop map application using embedded database such as spatialite

text: I am relatively very new to building GIS based applications and have a very limited understanding in this area. So I need some help in building a java based desktop application which would generate map tiles on the fly for a given tile coordinates (zoom level, x and y coordinate of the tile). I could build a reasonable solution using the below procedure using postgresql/postgis and geotools APIs. But now my client wants me to use the embedded database instead of postgresql/postgis since that eliminates the database administration. Here is what I have done so far:   1. Downloaded osm data and used osm2pgsql tool to populated the required osm tables in posgresql/postgis data.   2. I have used mapnik's SLD files to style the map objects (openstreetmap.org serves the tiles generated by mapnik). I have used mapnik to generate the SQL script that populates the map objects tables in posgresql/postgis database using osm tables in step 1 and the SLD files that can be used by geotools API to syle the objects. By map objects, I mean objects like roads, buildings, addresses etc which are map features that get loaded into each layer of the map.   3. geotools has APIs (StreamingRenderer) to generate the raster image of the bounding box that corresponds to a given tile coordinate. The generated raster image gets served as the tile. I have started exploring spatialite (spatial extension to sqlite) since I found that geotools has a plugin for spatialite which will load the map objects into the map. I am able to load the OSM data into spatialite tables using the spatialite-tool "spatialite_osm_map". Looking at the tables generated in spatialite, it seems like all map objects got populated. But I am not sure how to render these map objects of spatialite. I am not able to use the SLD files that were meant for postgis map objects. I would want to leverage the styling done for postgis, but that does not seem to be possible. I have googled quite a bit and did not find anyone using SLD files to style the spatialite map objects. Typically one uses desktop application such as QGIS to show the map objects of spatialite. The spatialite 4.1.0 has started supporting SLD styling, but that version is not compatible with geotools plugin for spatialite which relies on version 2.4.0. So with all these constraints, I guess I need to adopt entirely different approach of using embedded database. May be some other embedded databases. I looked at this discussion in this site, but among all of them I found spatialite to be more popular and being actively used and hence started exploring it. Has anyone used any other embedded database to solve this kind of problem? Can someone please guide me through this? That would be really helpful and highly appreciated. Any alternate approaches, thoughts are most welcome. I guess one way could be storing the already generated tiles in sqlite (conforming to MBtiles spec)? But that may not be acceptable by my client since he wants tiles to be generated on the fly. Thanks and Regards, Janaki @BradHards - Ok, I will try to be more specific. I would like use a embedded database (for e.g. spatialite) in my java desktop application to store the map data (osm data) instead of postgis/postgresql and be able to generate the tile on the fly for a given tile coordinate using the stored map data. The generated tile should look similar to what one sees in openstreetmap.org. Any suggestions to achieve this workflow? It seems QGIS does not out of the box support for java based on this discussion. I did look at spatiliate link that you mentioned - The SLD styling is supported in spatialite 4.1.0 and the geotools-spatialite plugin which I am presently using is compatible with spatialite 2.4.0. However I can go ahead and use spatialite 4.1.0, but I did not find any java framework that can possibly help me in achieving the workflow of reading the map data and generate tiles similar to openstreetmap. @blah238 - Thanks for the pointer. I will take dig at uDig and see if Neo4j spatial can be used in my case. Regarding on-the-fly generation of tiles - My client does not want to store all the tiles for various zoom levels (required zoom level is around 20) in the database or file system since that would be a huge data to be stored and is impracticable for him.


title: How to access Bing Map Metadata in ArcGIS Desktop?

text: Is there any way to access the metadata of Bing Maps in ArcGIS desktop? Any good plug-in or extension?


title: How to add multiple markers on map and popup different information using OpenLayers?

text: Can anybody give me some examples of adding multiple markers on the map?


title: QGIS - delete polys based on vector line layer - possible?

text: I have a layer of polygons and a layer of vector lines. I want to delete any poly that is touched/crossed/intersects a vector line. Is this possible? I have tried the Geoprocessing Tools. Looks like I want 'Difference' but the best I can do is delete segments of lines that touch polys. I want the opposite.


title: How to sum up all the single points in a state

text: I am new to using ArcMap. I have a list of locations with lat/long. I'm wondering if it is possible for ArcGIS to compute the sum of the locations by state and then shade in the state depending on the total amount of locations.


title: vector file is used by two mapsets at a time in grass gis

text: i want to study about vector processing.if a vector file is used by two mapsets at a time, which one will be saved first and why.if a team want to work on the same set of shape files, what is the procedure. will it be saved in dbf format?


title: open current selection's tab in tabbed infoWindow container: ArcGIS Javascript API

text: I am updating a webApp using ArcGIS Javscript API 2.1 I have an infoWindow with tabs which shows all of the features in my map. I built it based on this example: http://resources.esri.com/help/9.3/arcgisserver/apis/javascript/arcgis/help/jssamples_start.htm#jssamples/identify_drilldown.html I would like to have the infoWindow open with the tab that matches the selected feature on top (or, selected). Currently, in my map and in the above example, the tab that opens is the tab last clicked during the preceding instance of the infoWindow. I see there is a 'selected' property for tab containers. I can set it for the initial layout of the infoWindow onload in HTML, but setting it in the Javascript code does not seem to work.. Does anyone have an idea how I may do this? Thanks,


title: st_intersection; What Im doing wrong?

text: Im trying to run this query :               create TABLE data.agri_grc as select     id, code_00,     ST_Intersection(data.agri.geom,data.grc_bor.geom) as geom from     data.agri;      But im getting this error:               ERROR:  missing FROM-clause entry for table "grc_bor"     LINE 3: ST_Intersection(data.agri.geom,data.grc_bor.geom) as geom fr...      Whats I am missing? In postgis2 db


title: Error : keyword can't be an expression

text: I am getting "keyword can't be an expression" error while running the following code:               >>> import arcpy     >>> arcpy.MakeFeatureLayer_management("Point","P1","Layer" = 'WM-VALVE-ID')      I ran "Make Feature Layer" tool from the toolbox which worked. Below is the python snippet:               arcpy.MakeFeatureLayer_management("Point","Point_Layer","/\Layer/\ = 'WM-VALVE-ID'","#","OBJECTID OBJECTID VISIBLE NONE;Shape Shape VISIBLE NONE;Entity Entity VISIBLE NONE;Handle Handle VISIBLE NONE;Layer Layer VISIBLE NONE;LyrFrzn LyrFrzn VISIBLE NONE;LyrLock LyrLock VISIBLE NONE;LyrOn LyrOn VISIBLE NONE;LyrVPFrzn LyrVPFrzn VISIBLE NONE;LyrHandle LyrHandle VISIBLE NONE;Color Color VISIBLE NONE;EntColor EntColor VISIBLE NONE;LyrColor LyrColor VISIBLE NONE;BlkColor BlkColor VISIBLE NONE;Linetype Linetype VISIBLE NONE;EntLinetype EntLinetype VISIBLE NONE;LyrLnType LyrLnType VISIBLE NONE;BlkLinetype BlkLinetype VISIBLE NONE;Elevation Elevation VISIBLE NONE;Thickness Thickness VISIBLE NONE;LineWt LineWt VISIBLE NONE;EntLineWt EntLineWt VISIBLE NONE;LyrLineWt LyrLineWt VISIBLE NONE;BlkLineWt BlkLineWt VISIBLE NONE;RefName RefName VISIBLE NONE;LTScale LTScale VISIBLE NONE;Angle Angle VISIBLE NONE;ExtX ExtX VISIBLE NONE;ExtY ExtY VISIBLE NONE;ExtZ ExtZ VISIBLE NONE;DocName DocName VISIBLE NONE;DocPath DocPath VISIBLE NONE;DocType DocType VISIBLE NONE;DocVer DocVer VISIBLE NONE;ScaleX ScaleX VISIBLE NONE;ScaleY ScaleY VISIBLE NONE;ScaleZ ScaleZ VISIBLE NONE;ID ID VISIBLE NONE;ID_ ID_ VISIBLE NONE;B_P B_P VISIBLE NONE")      So, when I copied the Python snippet into my code and tried running it, I got another error "ERROR 000230: Failed selecting with /\Layer/\ = 'WM-VALVE-ID' Failed to execute (MakeFeatureLayer)." Any help is very much appreciated. Thanks Sam


title: Is there a way to get accurate daily U.S. populations?

text: I suppose you could interpolate from the Census data, but I am interested in a data feed that has accurate DAILY total population of the U.S. The finer the granularity, the better (so if there was daily populations by state or county that would be great). I know this is a somewhat weird/obscure question, but thanks.


title: In ArcMap, what is the difference between Copy Parallel and Buffer?

text: To build an offset of a linear feature in ArcMap, what is the difference between Copy Parallel and Buffer? In terms of the resulting features and their linear differences from the source linear feature? I am building an easement of a centerline feature. I want to build a feature that is 25feet either side of my linear feature. Where the centerline is going without direction changes there seems to be no difference. Where the centerline changes direction the differences start to appear. Thanks.


title: How to find qgis2 folder to delete in ubuntu

text: How can i find the qgis2 in qgis and where it is installed, i had one error on when qgis open this error will occurred "[ERROR] Can not make qgis.db private copy" some links told delete the folder qgis2, i can't find this folder how to solve this problem help me in advance.


title: SQL to get points within a certain distance

text: I am new to using spatial data types in sql server. I have managed to create a column of Approx 50,0000 geometry points (from x,y data) and would like to find out which (count number of) points are within a radius of 50m from another point. I have read a ton of blog posts and advice pages and think that I need to do this by creating a polygon layer of circles of radius 50m from each point, create a spatial index on the polygon layer and then do a count of all points in each polygon. The problem is that I can't work out how to this using my point data to create the new polygon layer. The columns of data I have are 'ref','x', 'y', 'geompointloc' I would also also be interested to know if this is the most efficient way to complete this task Thanks


title: MapServer/Mapcache Getting "This XML File does not appear to have any style associated with it" when trying to open map

text: I have installed mapserver and mapcache using the guide here: http://trac.osgeo.org/mapserver/wiki/RenderingOsmDataUbuntu I have successfully generated a map, and verified that it works using the shp2img command. When I navigate to:               http://localhost:9282/mapcache/?MAP=/gisDirectory/htdocs/mapserver/osm-outlined,google,usshields.map&SERVICE=WMS&REQUEST=GetCapabilities      I simply see get the error:               This XML file does not appear to have any style information associated with it. The document tree is shown below.      The document tree is as follows:               <WMT_MS_Capabilities version="1.1.1">         <Service>             <Name>OGC:WMS</Name>             <Title>no title set, add some in metadata</Title>             <OnlineResource xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://localhost:9282/mapcache/?"/>         </Service>         <Capability>             <Request>                 <GetCapabilities>                     <Format>application/vnd.ogc.wms_xml</Format>                     <DCPType>                     <HTTP>                         <Get>                             <OnlineResource xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://localhost:9282/mapcache/?"/>                     </Get>                    </HTTP>                    </DCPType>                </GetCapabilities>                 <GetMap>                     <Format>image/png</Format>                     <Format>image/jpeg</Format>                     <DCPType>                          <HTTP>                             <Get>                                 <OnlineResource xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://localhost:9282/mapcache/?"/>                            </Get>                         </HTTP>                     </DCPType>                 </GetMap>                 <GetFeatureInfo>                     <Format>text/plain</Format>                     <Format>application/vnd.ogc.gml</Format>                         <DCPType>                             <HTTP>                                 <Get>                                     <OnlineResource xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://localhost:9282/mapcache/?"/>                                 </Get>                             </HTTP>                         </DCPType>                 </GetFeatureInfo>             </Request>             <Exceptions>                 <Format>text/plain</Format>             </Exceptions>             <VendorSpecificCapabilities>                 <TileSet>                     <SRS>EPSG:900913</SRS>                     <BoundingBox SRS="EPSG:900913" minx="-20037508.342789" miny="-20037508.342789" maxx="20037508.342789" maxy="20037508.342789"/>                     <Resolutions>                        156543.03392804099712520838 78271.51696402048401068896 39135.75848201022745342925 19567.87924100512100267224 9783.93962050256050133612 4891.96981025128025066806 2445.98490512564012533403 1222.99245256282006266702 611.49622628141003133351 305.74811314070478829308 152.87405657035250783338 76.43702828517623970583 38.21851414258812695834 19.10925707129405992646 9.55462853564703173959 4.77731426782351586979 2.38865713391175793490 1.19432856695587897633 0.59716428347793950593                     </Resolutions>                     <Width>256</Width>                     <Height>256</Height>                     <Format>image/png</Format>                     <Layers>osm</Layers>                     <Styles/>                 </TileSet>             </VendorSpecificCapabilities>         <Layer>             <Title>no title set, add some in metadata</Title>             <SRS>EPSG:900913</SRS>             <SRS>EPSG:4326</SRS>             <SRS>EPSG:3857</SRS>             <Layer cascaded="1" queryable="0">                 <Name>osm</Name>                 <Title>OSM MapServer served map</Title>                <Abstract>                     see http://trac.osgeo.org/mapserver/wiki/RenderingOsmDataUbuntu                </Abstract>                 <LatLonBoundingBox minx="-180.000000" miny="-90.000000" maxx="180.000000" maxy="90.000000"/>                 <BoundingBox SRS="EPSG:900913" minx="-20037508.342789" miny="-20037508.342789" maxx="20037508.342789" maxy="20037508.342789"/>                 <SRS>EPSG:900913</SRS>                 <SRS>EPSG:3857</SRS>                 <BoundingBox SRS="EPSG:4326" minx="-180.000000" miny="-90.000000" maxx="180.000000" maxy="90.000000"/>                 <SRS>EPSG:4326</SRS>             </Layer>         </Layer>     </Capability>


title: uDig window to change between historical dates?

text: I'm new to uDig, but I'm experienced Java developer. I have points with temperature measurements that change overtime in SQL database. I would like to give user ability to change between historical data. Is there window in uDig which allows this? If not: Can you point me where to start to write a plugin allowing this?


title: Syntax for using ArcGIS for Server layer with WMS service?

text: I am having some issues with getting an image/png format from an ArcGIS Server. The console returns:               Resource interpreted as Image but transferred with MIME type application/vnd.ogc.se_xml:     "http://meshekgis.com/ArcGIS/services/mustang/Mustang_Map/mapserver/WMSServe…3A900913&BBOX=15028131.255,-5009377.085,20037508.34,0&WIDTH=256&HEIGHT=256".       Here is the javascript I am using:               var arcserver_wms =   "http://myserver.com/ArcGIS/services/mustang/Mustang_Map/mapserver/WMSServer";          var arcserver_wms = new OpenLayers.Layer.WMS("ArcServer WMS", arcserver_wms,                          {                            layers:"0",                            FORMAT: "image/png-8"                         },                         {                             isBaseLayer: false,                             opacity: 0.5                         });


title: How to import ESRI Ascii format into GRASS?

text: I need to take a set of datafiles, generated by our sensor system that reads GPR (Ground penetrating Radar) deviations over a surface as we drive along. I have discrete readings from a set of sensors and can accurately place them relative to the GPS unit spitting out location values once a second, along with track and speed in meters/second. I want to create an Esri ASCII data file for each second, with each row representing the sensor values along the rail, with one row representing all 16 sensors (spaced equally apart), read every 1/25th of a second. I can figure out the offsets to accurately place the sensors, in meters, relative to the GPS unit and as each sensor is read. It is relatively easy to then figure out the location for each reading along the track the GPS unit is measuring in Meters/second, thus arriving at the separation, in meters, between rows of generated sensor data. What I am having a hard time doing is converting this information into ESRI ASCII format, specifically the CELL size. My sensor rack is about 4 meters long and the GPS unit is centrally located with a small offset forward of the sensor array. How can I determine the cell size for my many rows of sensor data? I will locate this data with the GPS output, and have no problem calculating the boundaries for the one second time period. What has me stumped is the cell size. Example data: Sensors 1 - 16 232.2 233. 244.4 256.8 260.1 ..... 211.5 next row 0.1 meter later next row 0.1 meter later 25 rows in all Sensors separated 0.25 meter from each other Any ideas folks? What does that pesky cell size have to be and more importantly, how did you figure out that number?


title: How to export tracks of a specific polygon in overpass-turbo?

text: I would like to export some features (in this case `route=hiking`) that are inside a specific boundary. So not the boundary box that I can choose through the pan tool, but all the routes that for example belong to a national park. Can I do this with the wizard? I've founded this example in the OSM wiki. I'm trying to perform something like that.


title: Using GIS / Map inside c# WPF

text: How could i create e mapping applications that displays a custom map from this URL: (Capabilties with a mapview) http://openwms.statkart.no/skwms1/wms.sjo_hovedkart2?request=GetCapabilities&Service=WMS inside c# with WPF controls. It would be nice if there was some full feature map libraries so i could add pinpoints to certain Long,lats e.e (Like Openlayers for Javascript) Google / bing maps API cannot be used as you are not allowed to replace their maps.


title: Help needed with the rendering of QGIS Heatmaps

text: Please could someone give me some advice on creating heatmaps in QGIS? I have tried to create these through both QGIS and GRASS, and fail miserably each time. In QGIS, I either receive an error that not all the fields have been filled in, which I most certainly have done, or it fails to render anything remotely resembling a heatmap. In GRASS, it takes hours to render and the final result is just a bunch of colourful, concentric circles. What am I doing wrong? * * *   1. I'm using QGIS 1.8.0-Lisboa, using the heatmap plugin via Raster ---> Heatmap ---> Heatmap    2. The data I am using is fire point data from NASA's Aqua and Terra satellites. As far as workflow is concerned, I've tried the 3 methods shown in the tutorial:  a) Open vector point shapefile, then select 'Heatmap' as per above, provide an output raster name, set the radius as 100000 (100km), decay ratio as 0, click on advanced, x=100, y=100. This usually produces a blank or a grey rectangle (with a lot of holes in it). b) Using the GRASS v.kernel module: Open the mapset (created previously). Go to GRASS Tools and select v.kernel, then select the input vector and enter 1000 as the Standard deviation, but when I run it, an error appears which says the input vector and standard deviation have not been entered. c) From GRASS tools, I select GRASS shell, then enter v.kernel in the commandline interface. This takes hours and is what delivered the concentric circles. Not sure how to share the data? Still new on the forum and wasn’t allowed to include a screenshot.


title: Openlayers loop through json array and output result with as javascript

text: Connection.php file which stored database connection code with json:               <?php     $username = "root";     $password = "";     $database = "roadmap";     $url = "localhost";          $con = mysql_connect($url, $username , $password);     if (!$con)     {     die('Could not connect:' . mysql_error());     }          mysql_select_db($database, $con) ;     $data = array();     // query your DataBase here looking for a match to $input     $query = mysql_query("SELECT * FROM students");     while ($row = mysql_fetch_assoc($query)) {     $json = array();     $json['Id'] = $row['Id'];     $json['Lon'] = $row['Lon'];     $json['Lat'] = $row['Lat'];     $data[] = $json;     }     header("Content-type: application/json");     echo json_encode($data);     mysql_close($con);         ?>      In separate file index.php:               $.getJSON("connection.php", function(data) {                 addMarker( data[0].Lon,data[0].Lat ,data[0].desc);                 addMarker(-0.13264,51.50918 , );                 addMarker( -0.12498,51.50807 , );                 center = bounds.getCenterLonLat();                 map.setCenter(center, map.getZoomForExtent(bounds) - 1);                 zoom = map.getZoom();                  });       this currently pulls the json data through and I have to use it for each individual assigning of a part from the database such as data[0].lon will get the longitude entry for the first record with in the database. As I am going to be working with quite a few entry's. I want to loop through the json array and then out put all records as so:               addMarker( RECORD 1 );          addMarker( RECORD 2 );          addMarker( RECORD 3 );      Can anyone help me with this.


title: Remove first n characters in label using ArcGIS 10 label expression builder with VBscript

text: I have a planning cadastre dataset. It has the type of plan (DP, RP or SP) and then the plan number, thusly:               DP123456     RP1248590000     SP3855     etc...      I want to make a label in the expression builder in ArcGIS 10 which skips the first two characters and only has the plan number. In full, this is what I have done:               left([planlabel],2)&vbnewline&left([planlabel],len([planlabel]-2))&vbnewline& [lotnumber]       You can see here that I am making a stacked label which will put the type of plan on the first row (it is always the first two characters of [planlabel] and then the second row has the plan number, followed by [lotnumber] on the last line. In Excel, string indexing of this sort is apparently achieved by CELL,len(CELL-2) and I was hoping the same would hold for VBscript in ArcGIS but apparently this is not the case. Any clues?


title: Can GDAL Perform Vertical Transformations?

text: Can GDAL perform vertical transformations? If the answer is yes, please indicate which amongst these **listed pages** can explain how exactly to accomplish vertical transformations? If there is some other page discussing how exactly to accomplish vertical transformations, please indicate URL. This question arises after reading Jason's question Vertical Datum Conversion Utilities from a couple of years ago and the partial answer given.


title: Geoservices REST JSON to SQL Server Geometry using Python

text: I'd like to hit a Geoservices REST endpoint (below) to return json and push this into a SQL Server Geometry table using Python. Ok so I can achieve this via ogr using below via command line: ogr2ogr -overwrite -nln "LCCGIS.CAD_Parcels_JSON" -f "MSSQLSpatial" "MSSQL:s erver=SQL2008P1\MSSQLSERVER32BIT;database=GISTEST;Uid=lccgis;pwd=password" "json" -lco "GEOM_NAME=Shape" -lco "SRID=28355" -lco "DIM=2" Can someone help with how I do this using **python bindings for GDAL** or point me to a forum etc?... table_driver = ogr.GetDriverByName( 'MSSQLSpatial' ) table_ds = shp_driver.CreateDataSource( "MSSQL:server=SQL2008P1\MSSQLSERVER32BIT;database=GISTEST;Uid=lccgis;pwd=password") out_layer_name = "LCC.CAD_Parcels_JSON" lco = ["GEOM_NAME=Shape","DIM=2"] table_layer = shp_ds.CreateLayer( out_layer_name, geom_type = in_defn.GetGeomType(), srs = in_layer.GetSpatialRef(), options = lco ) geom_type and srs seem ok so I suspect options is the culprit. My table_layer ends up being a NoneType


title: How do you get all the results from a Selection made by Attribute or Location using ArcPy?

text: Using arcpy I'm trying to get the results of a selection, or what is currently selected as a list. Basically I'm attempting to build an addin which on click will either perform an operation where I invoke SelectLayerByAttribute_management (or ByLocation) to select some features, or the user performs a selection ByAttribute (or ByLocation), and while they have their selection hightlighted executes my snippet by pressing a button on the toolbar. In each case I want to get a list of results to perform a specific task on only certain selected features. I can't find out how to get a list of the currently selected features, or find a listener for selection having been made so I can enable the button, and ensure I can perform my operation on an existing selection. Thanks


title: Arcpy - What's the best way to check whether feature class has selection?

text: In a sentence, **how do you test in Arcpy whether a feature class has any selected records or not**? More info below to explain: I am doing a batch process in which a sampling layer (10km grids) is used to select polygon records (of which there are about 100,000 in total). My script returns the list of polygon IDs that intersect each 10km grid square. I then do some further processing/summarizing using that list of polygon IDs. In instances where there are polygons within the sample square, ArcGIS performs a selection on the polygons, and it's relatively easy to iterate through those records using a searchcursor. Note I'm using v10.0. But in instances where there are no selected records (ie there are no polygons intersecting my 10km sample square), the searchcursor would return ALL records in the database. So what's the most robust way to test programatically whether a feature class has a selection or not? I know you can use arcpy.GetCount_management to return a number that tells you how many records are selected, and if that number = the number of records in the entire layer then you can infer that either none, or all are selected. But that's not the same as determining if there is a selection or not. What's the best way to do that in Arcpy? Thank you


title: Is there an event handler for a Selection being made using ArcPy?

text: I want to be able to capture an arcpy.SelectByAttribute() or arcpy.SelectByLocation() event being triggered. There doesn't seem to be an extension that is able to handle this. Is there a way to do this?


title: Create a simple WPS client

text: It looks like the only way to get flat files/data from a WPS service is to use a client. I'm looking to create a super-simple WPS client in Python 2.6 (or other language) without any dependencies on non-default modules. Right now I'm gathering resources - I haven't done a lick of coding yet. Can I use a POST HTTP request with a proper header to request a WPS file? Pointers in the right direction would be appreciated. # Additional Info:   * There's an OWS LIB example client, but that requires an OWSlibrary.   * There's a Javascript/Open Layers client from PyWPS  - this looks like its the most independent of external dependencies sans the openlayers.js , but then ... well, openlayers is its own client.    * I have the Firefox HTTP Resource Test for POST and GET requests   * I have FireBug


title: Custom python datasource to QGIS

text: I'm working on a homebrew spatial feature data format, and using Python to write it's API. It occoured to me that use some tool like QGIS to display my data would be very nice to have, for debug/test purpose. Though based on QGIS's mail archive this and this, QGIS's Python binding to data provider is access-only, not something can be extended. And I have not dig deep into OGR DataSource Implementation, but it seems extending that is not very hack-friendly either, and limited to c++. Then I found Web Feature Service (WFS) is something can be done in Python, but the spec is complicated. The most decent Python WFS Server I found is featureserver, for now my best try would be hacking featureserver. Is there any other shortcut that can achieve my goal (custom python datasource to qgis layer) out there that I missed? Please don't tell me to dump my data to postgis/GeoJson/Shapefile, I'm aware of that and that's not my point... Thanks.


title: Openlayers Css problem?

text: I have succeeded to show `Openlayers` map in my application, when i show css popup(not Openlayers Popup) map has to blank as the rest of the application. But `Openlayers` map is not getting blank. What could be the problem? Instead of `Openlayers` if i use `Google Maps` its getting blank when i show popup. I am not using any css related to `Openlayers`, just i included this http://openlayers.org/api/OpenLayers.js


title: How to reset MacOS QGIS to defaults?

text: I can't find a way to clear away all imported QGIS styles, even completely removing all QGIS files, after new install - old styles are still there. Where to look? Removing one by one is the only option from within the application, but there are hundreds that I have imported. I found the "first run state" question, but that answer works for Windows, but how about MacOS?


title: Using os.walk to update broken links and save mxd

text: I've got code to walk through an entire directory and sub-directories repair broken links in each mxd and save the mxd. It works for the mxd's in the top directory but fails when trying to enter any sub. My code is as follows:               import arcpy, os     from arcpy import env     arcpy.env.overwriteOutput = True     ws = "C:\Users\mittcla\Desktop\New folder"     oldpath = r"C:\NewProjects"     oldbase = r"C:\Base Data"     newpath = r"C:\Users\mittcla\Documents\NewProjects"     newbase = r"C:\Users\mittcla\Documents\Base Data"     for root, dirs, files in os.walk(ws):         for f in files:             if f.endswith('.mxd'):                 fullpath = os.path.join(ws,f)                 mxd = arcpy.mapping.MapDocument(fullpath)                 print "Replacing path for " + f +"..."                 mxd.findAndReplaceWorkspacePaths(oldpath, newpath)                 mxd.findAndReplaceWorkspacePaths(oldbase, newbase)                 mxd.save()                 del mxd      I get an error saying > Runtime error Traceback (most recent call last): File "", line 14, in File > "c:\program files > (x86)\arcgis\desktop10.2\arcpy\arcpy\arcobjects\mixins.py", line 609, in > **init** assert (os.path.isfile(mxd) or (mxd.lower() == "current")), > gp.getIDMessage(89004, "Invalid MXD filename") AssertionError: Invalid MXD > filename. What is causing it to break when entering the new directory?


title: Trouble downloading Geodatabase wetlands from FWS.gov

text: Have tried using every browser (chrome, firefox, safari, ie, opera) to download wetland data from: FWS.gov wetland data download with no luck. It either chokes around 400 MB or finishes up with a corrupted zip file which cannot be opened. It is only an issue with the larger .gdb files (over 400 MB), like Texas which is 1.7 GB. Any ideas on how to download these files without error?


title: Shapefiles of Central America and Mexico coastlines?

text: Does anybody know/has a shapefile of Central America + Mexico (WSG84) (with no administrative boundaries) at a low resolution that could share with me? I'm looking for a very simple shapefile without the small islands and smoothed out coastlines. Any help is greatly appreciated!


title: "Can't extract geo keys from object, malformed geometry?" Possible self intersection?

text: My '2dsphere' index creation fails at a document. The command i write for creation of index is:               db.collection.ensureIndex({'geometry' : "2dsphere"})      It returns the following error:               {         "err" : "Can't extract geo keys from object, malformed geometry?:...",         "code" : 16755,         "n" : 0,         "singleShard" : "localhost:27026",         "ok" : 1     }      The document at which the index creation fails is:               {         "_id" : ObjectId("53978003e6919576808d12ad"),         "type" : "Feature",         "properties" : {             "NumPoints" : 185         },         "geometry" : {             "type" : "Polygon",             "coordinates" : [                  [                      [                          -123.751691277,                          39.26661433100003                     ],                      [                          -123.75229746,                          39.26815238500006                     ],                      [                          -123.751289221,                          39.26838356800005                     ],                      [                          -123.751592313,                          39.26915259600003                     ],                      [                          -123.752600561,                          39.26892141100007                     ],                      [                          -123.75229746,                          39.26815238500006                     ],                      [                          -123.753305692,                          39.26792119500004                     ],                      [                          -123.753002587,                          39.26715217100008                     ],                      [                          -123.756027215,                          39.26645855800007                     ],                      [                          -123.755724088,                          39.26568954100003                     ],                      [                          -123.757740455,                          39.26522709700004                     ],                      [                          -123.7580436009999,                          39.26599610900007                     ],                      [                          -123.760059962,                          39.26553362800007                     ],                      [                          -123.760363134,                          39.26630263500005                     ],                      [                          -123.761371315,                          39.26607137900004                     ],                      [                          -123.761674503,                          39.26684038300004                     ],                      [                          -123.762682687,                          39.26660911700003                     ],                      [                          -123.762985891,                          39.26737811700008                     ],                      [                          -123.7639940779999,                          39.26714684000007                     ],                      [                          -123.764297299,                          39.26791583800008                     ],                      [                          -123.765305489,                          39.26768455000007                     ],                      [                          -123.765608726,                          39.26845354500006                     ],                      [                          -123.76661692,                          39.26822224700004                     ],                      [                          -123.766920173,                          39.26899123800007                     ],                      [                          -123.76792837,                          39.26875992900006                     ],                      [                          -123.768231639,                          39.26952891700006                     ],                      [                          -123.7712562199999,                          39.26883493500003                     ],                      [                          -123.771559525,                          39.26960391500006                     ],                      [                          -123.772567715,                          39.26937256900004                     ],                      [                          -123.772871036,                          39.27014154600005                     ],                      [                          -123.77387923,                          39.26991018900003                     ],                      [                          -123.774182567,                          39.27067916300007                     ],                      [                          -123.776198955,                          39.27021641900006                     ],                      [                          -123.776502318,                          39.27098538800004                     ],                      [                          -123.7785187,                          39.27052260700003                     ],                      [                          -123.778822089,                          39.27129157000007                     ],                      [                          -123.781846642,                          39.27059733000004                     ],                      [                          -123.782150067,                          39.27136628400007                     ],                      [                          -123.784166423,                          39.27090341200005                     ],                      [                          -123.784469873,                          39.27167236100007                     ],                      [                          -123.785478051,                          39.27144091000008                     ],                      [                          -123.785781518,                          39.27220985600008                     ],                      [                          -123.7867896989999,                          39.27197839500008                     ],                      [                          -123.787093182,                          39.27274733800005                     ],                      [                          -123.788101367,                          39.27251586500006                     ],                      [                          -123.788404866,                          39.27328480500006                     ],                      [                          -123.7904212349999,                          39.27282183100004                     ],                      [                          -123.79072476,                          39.27359076600004                     ],                      [                          -123.791732945,                          39.27335926400008                     ],                      [                          -123.792036486,                          39.27412819500006                     ],                      [                          -123.793044674,                          39.27389668300003                     ],                      [                          -123.793348231,                          39.27466561100005                     ],                      [                          -123.794356423,                          39.27443408800008                     ],                      [                          -123.794659996,                          39.27520301300007                     ],                      [                          -123.795668191,                          39.27497147900004                     ],                      [                          -123.79597178,                          39.27574040100006                     ],                      [                          -123.796979978,                          39.27550885700003                     ],                      [                          -123.797283584,                          39.27627777600003                     ],                      [                          -123.796275376,                          39.27650932200004                     ],                      [                          -123.796578979,                          39.27727824300007                     ],                      [                          -123.795570755,                          39.27750978400007                     ],                      [                          -123.795874354,                          39.27827870600004                     ],                      [                          -123.7948661139999,                          39.27851024200004                     ],                      [                          -123.79516971,                          39.27927916600004                     ],                      [                          -123.794161454,                          39.27951069600005                     ],                      [                          -123.794465046,                          39.28027962200008                     ],                      [                          -123.793456774,                          39.28051114600004                     ],                      [                          -123.793760363,                          39.28128007500004                     ],                      [                          -123.79174378,                          39.28174310300005                     ],                      [                          -123.792047356,                          39.28251203600007                     ],                      [                          -123.790030727,                          39.28297503700003                     ],                      [                          -123.79033429,                          39.28374397400006                     ],                      [                          -123.777225449,                          39.28675272000004                     ],                      [                          -123.7769220119999,                          39.28598375100006                     ],                      [                          -123.759778205,                          39.28991614500006                     ],                      [                          -123.759171682,                          39.28837811900007                     ],                      [                          -123.758163184,                          39.28860935800003                     ],                      [                          -123.757859942,                          39.28784034100005                     ],                      [                          -123.755842946,                          39.28830279000005                     ],                      [                          -123.756146169,                          39.28907181100004                     ],                      [                          -123.755137652,                          39.28930302600003                     ],                      [                          -123.754834439,                          39.28853400200006                     ],                      [                          -123.7538259249999,                          39.28876520600005                     ],                      [                          -123.753522728,                          39.28799617900006                     ],                      [                          -123.750497178,                          39.28868973400006                     ],                      [                          -123.749890862,                          39.28715166400008                     ],                      [                          -123.748882351,                          39.28738282800003                     ],                      [                          -123.747972955,                          39.28507571000006                     ],                      [                          -123.7469644669999,                          39.28530685800007                     ],                      [                          -123.746055158,                          39.28299972700006                     ],                      [                          -123.745046693,                          39.28323086000006                     ],                      [                          -123.744743613,                          39.28246181300005                     ],                      [                          -123.7427266829999,                          39.28292404800004                     ],                      [                          -123.742423629,                          39.28215499600003                     ],                      [                          -123.741415164,                          39.28238609900006                     ],                      [                          -123.741112126,                          39.28161704300004                     ],                      [                          -123.739095197,                          39.28207922000007                     ],                      [                          -123.738792184,                          39.28131015800005                     ],                      [                          -123.73778372,                          39.28154123200005                     ],                      [                          -123.737177733,                          39.28000310200008                     ],                      [                          -123.735160824,                          39.28046521500005                     ],                      [                          -123.73485786,                          39.27969614400007                     ],                      [                          -123.733849406,                          39.27992718600007                     ],                      [                          -123.733546458,                          39.27915811200006                     ],                      [                          -123.732538007,                          39.27938914300006                     ],                      [                          -123.732235075,                          39.27862006600003                     ],                      [                          -123.730218173,                          39.27908209900005                     ],                      [                          -123.729915267,                          39.27831301700007                     ],                      [                          -123.727898359,                          39.27877501200004                     ],                      [                          -123.7275954789999,                          39.27800592400007                     ],                      [                          -123.726587024,                          39.27823690700006                     ],                      [                          -123.726284161,                          39.27746781600007                     ],                      [                          -123.718216374,                          39.27931536500006                     ],                      [                          -123.717913595,                          39.27854625300006                     ],                      [                          -123.716905103,                          39.27877715800003                     ],                      [                          -123.716602339,                          39.27800804300006                     ],                      [                          -123.71559385,                          39.27823893700003                     ],                      [                          -123.715291103,                          39.27746981900003                     ],                      [                          -123.714282617,                          39.27770070200006                     ],                      [                          -123.713979886,                          39.27693158100004                     ],                      [                          -123.711962915,                          39.27739331900006                     ],                      [                          -123.711357511,                          39.27585506500003                     ],                      [                          -123.710349035,                          39.27608591600006                     ],                      [                          -123.7100463529999,                          39.27531678600008                     ],                      [                          -123.70903788,                          39.27554762700004                     ],                      [                          -123.7087352139999,                          39.27477849400003                     ],                      [                          -123.706718269,                          39.27524014500006                     ],                      [                          -123.706415628,                          39.27447100700005                     ],                      [                          -123.703390192,                          39.27516341600006                     ],                      [                          -123.703087587,                          39.27439426900003                     ],                      [                          -123.7020791049999,                          39.27462505300008                     ],                      [                          -123.701776516,                          39.27385590300003                     ],                      [                          -123.700768038,                          39.27408667700007                     ],                      [                          -123.700465465,                          39.27331752400005                     ],                      [                          -123.69945699,                          39.27354828700004                     ],                      [                          -123.699154434,                          39.27277913100005                     ],                      [                          -123.695120507,                          39.27370209100008                     ],                      [                          -123.694817997,                          39.27293292500008                     ],                      [                          -123.693809509,                          39.27316364200004                     ],                      [                          -123.693507014,                          39.27239447200003                     ],                      [                          -123.6965324289999,                          39.27170230300004                     ],                      [                          -123.696229911,                          39.27093314000007                     ],                      [                          -123.698246803,                          39.27047165900007                     ],                      [                          -123.697944273,                          39.26970250000005                     ],                      [                          -123.698952699,                          39.26947175000004                     ],                      [                          -123.69925524,                          39.27024090600003                     ],                      [                          -123.70228051,                          39.26954859800003                     ],                      [                          -123.7019779409999,                          39.26877944900008                     ],                      [                          -123.702986342,                          39.26854866600007                     ],                      [                          -123.703288921,                          39.26931781200005                     ],                      [                          -123.707322499,                          39.26839458800004                     ],                      [                          -123.707625124,                          39.26916372400007                     ],                      [                          -123.710650269,                          39.26847121200007                     ],                      [                          -123.709742327,                          39.26616382700007                     ],                      [                          -123.710750667,                          39.26593298100005                     ],                      [                          -123.710145384,                          39.26439472600004                     ],                      [                          -123.719219978,                          39.26231679100005                     ],                      [                          -123.719522704,                          39.26308589700005                     ],                      [                          -123.72053097,                          39.26285497200007                     ],                      [                          -123.719925506,                          39.26131676600005                     ],                      [                          -123.720933746,                          39.26108583800004                     ],                      [                          -123.720631013,                          39.26031673600005                     ],                      [                          -123.72668026,                          39.25893101200006                     ],                      [                          -123.7269830499999,                          39.25970009900004                     ],                      [                          -123.730007616,                          39.25900711900005                     ],                      [                          -123.730310442,                          39.25977619800005                     ],                      [                          -123.7313186279999,                          39.25954518600003                     ],                      [                          -123.73162147,                          39.26031426100008                     ],                      [                          -123.7326296589999,                          39.26008323900004                     ],                      [                          -123.733235382,                          39.26162138300003                     ],                      [                          -123.734243584,                          39.26139034700003                     ],                      [                          -123.734546465,                          39.26215941600003                     ],                      [                          -123.740595599,                          39.26077301600003                     ],                      [                          -123.740898545,                          39.26154206900003                     ],                      [                          -123.741906721,                          39.26131097100006                     ],                      [                          -123.7422096829999,                          39.26208002100003                     ],                      [                          -123.743217863,                          39.26184891300005                     ],                      [                          -123.743520841,                          39.26261796000006                     ],                      [                          -123.744529024,                          39.26238684100008                     ],                      [                          -123.744832018,                          39.26315588500006                     ],                      [                          -123.749872884,                          39.26200015500007                     ],                      [                          -123.751085121,                          39.26507627500007                     ],                      [                          -123.752093313,                          39.26484509500006                     ],                      [                          -123.7526994899999,                          39.26638314600007                     ],                      [                          -123.751691277,                          39.26661433100003                     ]                 ]             ]         }     }      I projected this polygon on a map and I see that two points of this polygon overlap. They do not essentially intersect though but they do overlap. Is there any way to make this polygon work in MongoDB? I can surely modify the polygon to remove the overlapping point, but does anyone has a different idea?


title: find and Calculate centroid within polygon either in ArcMap or in MapInfo

text: How can i calculate single/multiple centroid(point) within polygon either in Arcmap or in MapInfo. Both object(point, polygon) are in defferent layer. Pls suggest..... Thanks


title: How to display Heatmap in Internet Explorer 7 using ArcGIS JavaScript API?

text: Trying to find a solution to heat maps on IE7/IE8 (i.e. not using Canvas/HTML5/heatmapapi). Any ideas?


title: QGIS linking census shapefiles with geocoded addresses

text: I'm trying to merge a census shapefile of New Orleans with a health outcomes dataset by address with QGIS. I have geocoded the addresses and added them as a layer in QGIS but don't know on which variable within the street level city shapefile to merge them. The dataset has long/lat but the shapefile does not... any help would be appreciated! Thanks!


title: How to determine number of layers in a given .mpk file and then enumerate it using ArcGis Runtime SDK for WPF 10.2

text: I need to determine number of layers in a given .mpk file and then enumerate it using ArcGis Runtime SDK for WPF 10.2. Can anyone provide me some link or code/snippet..thankyou


title: How do I only use a subset of features from Arcgis feature services in leaflet?

text: Is there a possibility to only use a selection of features based on an attribute from an arcgis feature service in leaflet? As you can set a query to a service, could this be done easily in leaflet js, using vector layer plug-in or esris plugin? Now I made all objects i don't want to show transparent which is heavy, dirty and makes the transparent ones clickable... Cheers Henrik


title: <#> doesn't give a useful distance

text: I'm new to postgresql/postgis and trying to measure a distance. So `st_distance()` may be the most accurate, it's the slowest too (as I've read). I read that one can use geometric (instead of geographic) measuring with less accuracy, but the numbers I get don't make sense. This is what I've got:               (all points are geography(Point,4326))          select       st_astext(point),       st_astext(point2),       round(st_distance(point, point2)),       round(point::geometry <#> point2::geometry) -- cast to geometry and use distance operator     from testgeo2 limit 50          "POINT(53.887 54.043)","POINT(70.764 80.922)","3053406","32" -- ok     "POINT(93.262 59.546)","POINT(84.331 84.334)","2774883","26" -- ok     "POINT(99.624 67.098)","POINT(95.718 85.252)","2027908","19" -- ok     "POINT(54.406 84.531)","POINT(80.27  83.385)", "326308","26" -- hm?!      So I don't get how this can be right. Any suggestions? One more example               -- new york, tokyo     insert into testgeo2(point, point2) values (        st_geographyfromtext('point(-74.007873 40.717602)'),         st_geographyfromtext('point(139.754333 35.675853)'))          -- get distance     select       st_astext(point),       st_astext(point2),       -- km, most accuracy       round(st_distance(point, point2)) / 1000,       -- km, with less accuracy       round(st_distance(point, point2, false)) / 1000,       -- km, Lambert       round(st_distance(         st_transform(point::geometry, 3587),         st_transform(point2::geometry, 3587))) / 1000,       -- km, from degree       round(point::geometry <#> point2::geometry) * 111     from testgeo2          -- results     POINT(-74.007873 40.717602) POINT(139.754333 35.675853),         10871, 10847, 11224, 23754      is geometry just not working on a larger scale?


title: How to change the default gray color of the QgsComposition background?

text: I'm building a map, with elements like QgsComposerMap of course, but also QgsComposerLabel, QgsComposerLegend ,etc... Each of these elements has a boolean "setBackgroundEnabled" property, to make them white backgrounded or not. But when aggregating these white backgrounded elements in the QgsComposition, we can see gray bands between them, caused by the deepest gray background. And it's printed when saved as image ! As for canvas, ( i mean canvas.setCanvasColor(Qt.white)), is there a way to set the color of the overall QgsComposition background (the deepest one) ? Is there something like :               comp = QgsComposition(mapRenderer)     comp.setCompositionColor(Qt.white)      ? Thank you


title: How to export data from sql database in the decimal format in QGIS?

text: Currently I am working with QGIS and fetching data from SQL server and export it in the format of shapefile in QGIS. While importing the data from sql server attribute table showing data with geography having decimal value. ![enter image description here](http://i.stack.imgur.com/zYyWI.jpg) But after exporting in the format of shapefile data in the attribute table not showing the decimal value. ![enter image description here](http://i.stack.imgur.com/ntCIF.jpg) according to article as follows decimal degree data erase itself between edit session in qgis in the description they mentioned to add new column and make decimal places precision as 5 or more than that, but attribute table showing for data imported from sql doesn't have option to make decimal places and add the precision as 5 or more than that. If we check this option in the attribute table of shapefile we can found that. But it not open properly. It got crashed while editing the column. Another question in front of me is how to copy all huge data to new column? ![enter image description here](http://i.stack.imgur.com/W7GYf.jpg) ![Attribute editing for sql data table](http://i.stack.imgur.com/jWmIv.jpg) I wanted to know:   1. How data got automatically converted in to integer type after exporting it from database   2. Does any new way to add decimal precision for adding new column in sql database imported attribute table and copy all data in the form of decimal places?   3. Does any tools, plug in to convert data type from float to float or any mapping type option for data type available?


title: How to run a new Service Area analysis

text: I will run a new service area analysis with Default Breaks at 30.60 and 90 minutes. I will try to do this in a run with ModelBuilder. Does anyone please help me out here. Also wonder how I select facilities?


title: Check if property exists on features in layer

text: I was wondering if there is a good way to check to see if a property exists on a particular feature layer in GeoServer? My use case is that I have a bunch of features I want to display on a map. Some of the feature layers will have a date property associated with them. If they have a date property I want to be able to display a control to filter the features shown by date. I noticed that WFS has an operation called `GetPropertyValue` which will return an exception if the property doesn't exist. However, if it does it (as expected) returns the value for that property for all features which can be very slow. My general algorithm is as follows:               var wfsLayer = new OpenLayers.Layer.Vector('MyVectorLayer', {         strategies: [new Openlayers.Strategy.BBOX()],         protocol: new OpenLayers.Protocol.WFS({             url: 'http://mygeoserver/geoserver/wfs',             featureType: 'MyVectorLayer',         })     });          // if layer has date property (this is the part I need help with)         var dateFilter = new OpenLayers.Filter.Comparison({             type: OpenLayers.Filter.Comparison.BETWEEN,             property: 'date',             lowerBoundary: startDateVal,             upperBoundary: endDateVal         });         wfsLayer.filter = dateFilter;     // end if


title: how to call the method of DescribeFeatureType in owslib

text: I want to execute `DescribeFeatureType` operation of wfs services. However, I can’t find the method to fulfill this operation in the owslib. Can anyone kind enough to tell me how to do this. The following is my code snippet:               from owslib.wms import  WebMapService     from owslib.wfs import  WebFeatureService as wfsServericelayer          wfslayer=wfsServericelayer('http://***:8080/geoserver/test/ows')     for i in wfslayer.operations:         print i.name          operationT=wfslayer.getOperationByName('DescribeFeatureType')  # how to execute this method???


title: SDE connection error

text: In my company we have a number of databases in an MS SQL instance. This morning only 1 database i generating following error when I try to do a direct-connect in ArcCatalog: Failed to connect to database. Underlying DBMS error[sde.DEFAULT][STATE_ID = 1] I have increased the number of allowed connections to 10 000 and that's is more than enough. Still the problem exists. Does anyone have a clue? :)


title: Merging multipolygon data into a polygon

text: I have a MultiPolygon representation of some shapes in POSTGIS, I am trying to convert them to Polygon data. I tried, `st_geometryn(r.geom, 1)` but that didnt work because the way polygons are embedded I think. Just returns a partial. Tried `ST_ConcaveHull(ST_Collect(r.geom), 0.99)`, no luck returns nothing Tried:               SELECT f.gid, ST_AsGeoJSON(ST_Union(f.geom)) as singlegeom FROM (SELECT gid,             (ST_Dump(geom)).geom As geom FROM public.nyc_boroughs ) As f GROUP BY gid"       but getting back to Multi Polygon, I feel like a little bit of direction could help me figure out. Data is here http://www.sendspace.com/file/2savyw


title: How to clip from a point layer using postgis?

text: I've loaded a LIDAR xyz file into a PostGIS2.0 database and created a 3D point geometry column ("tblCLOUD") from the data. I've also loaded into the database a shapefile that contains six (non- overlapping) polygons ("tblZONES"). What I'm trying to do now is extract the points covered by each polygon to an individual table. So the end result would be six tables containing N points each. I've hit my limits with PSQL! Thanks! EDIT: My two tables - tblcloud & tblpoly are stuctured as follows:               CREATE TABLE tblcloud     (       "pointID" integer NOT NULL DEFAULT nextval('"tblcloud_pointID_seq"'::regclass),       x double precision,       y double precision,       z double precision,       geom geometry(PointZ,2193),       CONSTRAINT pk_point_id PRIMARY KEY ("pointID" )     )     WITH (       OIDS=FALSE     );     ALTER TABLE tblcloud       OWNER TO postgres;          -- Index: tblcloud_geom_gist          -- DROP INDEX tblcloud_geom_gist;          CREATE INDEX tblcloud_geom_gist       ON tblcloud       USING gist       (geom );               CREATE TABLE tblpoly     (       gid serial NOT NULL,       minx numeric,       miny numeric,       maxx numeric,       maxy numeric,       cntx numeric,       cnty numeric,       area numeric,       perim numeric,       height numeric,       width numeric,       geom geometry(MultiPolygon,2193),       CONSTRAINT tblpoly_pkey PRIMARY KEY (gid )     )     WITH (       OIDS=FALSE     );     ALTER TABLE tblpoly       OWNER TO tom;          -- Index: tblpoly_geom_gist          -- DROP INDEX tblpoly_geom_gist;          CREATE INDEX tblpoly_geom_gist       ON tblpoly       USING gist       (geom );


title: How to check the resolution of converted raster file

text: I have a raster file whose resolution I need to change to 0.5. I do this using this: `g.region rast=my_raster res=0.5` After running this command, how do I verify whether the resolution has actually changed to 0.5?


title: Popup positioning issue in OpenLayers 2.10

text: We were using OpenLayers 2.7 before and everything was working fine. Then we decided to upgrade our system to use OpenLayers 2.10. We have popup positioning issue with this new version. The popup is getting displayed away from the actual marker. It is working fine if I roll back to use OpenLayers 2.7. I can not insert a snap here but you can have a look at the examples below. Ref. http://dev.openlayers.org/releases/OpenLayers-2.10/examples/popupMatrix.html I guess there is not an integration issue in my system/code as it's working fine with OpenLayers 2.7. Ref. http://dev.openlayers.org/releases/OpenLayers-2.7/examples/popupMatrix.html I checked it with all browsers, the behavior is same for all.


title: CalculateDefaultGridIndex

text: Has anyone found that using CalculateDefaultGridIndex in arcGIS9.3 does not return the ';' delimited string suggested in the documentation?


title: Using projections to obtain perspective views

text: I understand that a perspective view is essentially a type of projection. Is it possible to create projection definitions for a typical GIS system that will mimic a perspective view from a particular location? If so, how feasible would it be to use common 2D GIS tools to create perspective renderings of data, similar to what can be achieved in 3D tools such as Arc Scene and Google Earth?


title: How to find Clusters of N Points Within X Distance?

text: We have tried two methods of cluster analysis:   * Cluster and Outlier   * Hot Spot ... but neither meets our specification which has these 3 parameters:   * a Set of Point features   * N points within a Cluster   * X distance band to look in For example "find me the points where each cluster has 10 points within 50 meters" or something to that effect. ArcGIS 10 has the two tools I’ve linked above but neither seems to address the N points within a cluster part. Can our requirements be met with other tools in ArcGIS? Is what I am calling cluster analysis known as a particular clustering type? We did have a version of this in another application but I really wanted to use the tools within ArcGIS.


title: Creating new project on QT Creator using ArcGIS Engine controls and C++?

text: Can anyone explain me on how to create a new project on QT Creator and I need to call ArcGIS Engine controls into the Qt creator? I have looked into Getting started with Qt and ArcGIS Engine help. But I need to develop map control on QT Creator not using ArcGIS Runtime SDK for QT. I have installed ArcGIS Engine for C++.


title: Export directions from Google to KML?

text: I used Google to draw directions from one location to another location, but can't find a way to export it as KML so I can import it into my Android smartphone. Is there a way, or must directions be drawn in Google (Maps | MapsEngine)? Thank you.


title: Exporting selected Features to CAD drawing using ArcObjects and C#?

text: I'm new to ArcObjects, and I have to export selected features to CAD Drawing using just ArcObjects and C#. The code below works well when I convert to a shapefile `ShapeFileWorkspaceFactory`, but with the `CadWorkspaceFactory` it does not work. Any help will be greatly appreciated.               public void featureToCad(){                  IFeatureLayer fl = search_layer2(axMapControl1.Map, "parcelle") as IFeatureLayer;                  IFeatureClass pFC;                  IFeatureClassName pInFeatureClassName;             IDataset pDataset;             IDatasetName pInDsName;             IFeatureSelection pFSel;             ISelectionSet pSelSet;             IFeatureClassName pFeatureClassName;             IDatasetName pOutDatasetName;             IWorkspaceName pWorkspaceName;                  // *******************************             // Get the featureclass from the layer             // *******************************             pFC = fl.FeatureClass;             // ***********************************************             // Get the Featureclassname from the featureclass             // ***********************************************             pDataset = pFC as IDataset;             pInFeatureClassName = (IFeatureClassName)pDataset.FullName;             pInDsName = (IDatasetName)pInFeatureClassName;             // ***************************             // Get the selected features             // ***************************             pFSel = fl as IFeatureSelection;             pSelSet = pFSel.SelectionSet;             IWorkspace pScratchWorkspace;             IScratchWorkspaceFactory pScratchWorkspaceFactory;             pScratchWorkspaceFactory = new ScratchWorkspaceFactory();             pScratchWorkspace = pScratchWorkspaceFactory.DefaultScratchWorkspace;                  // ***************************************************             // Define the output featureclass             // give it the name of the input feature class + exp             // ***************************************************             pFeatureClassName = (IFeatureClassName)new FeatureClassName();             pOutDatasetName = (IDatasetName)pFeatureClassName;             pOutDatasetName.Name = (pFC.AliasName + "exp");             pWorkspaceName = (IWorkspaceName)new WorkspaceName();             pWorkspaceName.PathName = "D:\\";             pWorkspaceName.WorkspaceFactoryProgID = "esriDataSourcesFile.CadWorkspaceFactory";//it works with "ShapeFileWorkspaceFactory"             pOutDatasetName.WorkspaceName = pWorkspaceName;             // *************************************************************             // Give the output shapefile the same props as the input dataset             // *************************************************************             pFeatureClassName.FeatureType = pFC.FeatureType;             pFeatureClassName.ShapeType = pFC.ShapeType;             pFeatureClassName.ShapeFieldName = pFC.ShapeFieldName;             // *************************************************************             // Export selected features             // *************************************************************             pExportOp = new ExportOperation();             pExportOp.ExportFeatureClass(pInDsName, null, pSelSet, null,pFeatureClassName,0);          }


title: ArcGIS map service layer not loading properly

text: I'm at a total loss. My initial problem was an issue with the layer not loading at all when the page was set to the full extent, despite having changed the Min/Max Scale settings generously when publishing the map service. I corrected this by rectifying the Scale settings in JavaScript:               landBaseLayer.on("load", function() { landBaseLayer.minScale = 0; landBaseLayer.maxScale = 0; });       I now have a problem with the layer not loading when entering the page. It loads when I zoom in, and the loaded subset will remain visible when I zoom back out, but everything around it will remain blank. What am I doing wrong? Sorry, I'm totally new to this. I'm using a cached map service, with a Max Record Count of 10,000 (the layer contains just over 7,000 points). The following is the entire code.               <!DOCTYPE html>                      <meta http-equiv="Content-Type" content="text/html; charset=utf-8">                  <meta name="viewport" content="initial-scale=1, maximum-scale=1,user-scalable=no">         <title>Identify with Popup</title>          <link rel="stylesheet" href="http://js.arcgis.com/3.8/js/esri/css/esri.css">


title: How to style graduated using size-scale field

text: Is it possible to style based on size-scale field using graduated symbol but instead of using color ramp, I would like to display its size. I would also like to display this in the legend in the map composer. It would be nice if I could also apply this if I'll use svg marker (example 2nd image below). From this: ![enter image description here](http://i.stack.imgur.com/pjl5K.jpg) to something like this: ![enter image description here](http://i.stack.imgur.com/ntpMm.jpg)


title: eVis-Plugin centers mapwindow on coordinate origin

text: I have a QGIS-project with layers loaded from a postgis-database and also some shapefiles. I use the eVis-Plugin to display pictures.   The displaying is working fine but after clicking on the geometry to open the eVis-Browser the mapwindow section always jumps to the origin of the coordinatesystem (0/0). Also the asterisk from eVis showing which geometry is selected is located at (0/0)-but I don't have any geometries located there. Because of that I always have to zoom back to the map section with the geometries. How can I avoid this change of the mapwindow section? I have this problem no matter if the path of the photo comes from a shapefile- layer or from the database. The QGIS-project is defined with the coordinatesystem EPGS:21781 (CH1903, Switzerland). Because of that the geometries have coordinates like 600'000/200'000.


title: What dem data are Flex Viewers using to create elevation profiles?

text: I asked a similar question a year or two ago here: Can Elevation Profile data in Arc's Flex Viewer 3.0 be exported? which appears to not have an answer. Another question regarding the Flex Viewer (3.3 now) is where does it grab the dem data for these static elevation profile images it creates? Is this dem data a WMS that can be utilized in ArcGIS desktop ? ![screen shot of a Flex Viewer elevation profile](http://i.stack.imgur.com/dUx6i.png)


title: Creating simple maps to view on Mobile devices?

text: My organization wants to create a few maps that will be interactive and viewable via phones (android, iphones,etc) and ipads. We already have a solution using Geocortex (Silverlight) for browser which cannot display on most devices. I've been looking into HTML5, Javascript and ArcGIS Online. What I need is something that displays nicely for mobile use and that is easy to set up. Looking at Javascript is slightly intimidating to me and ArcGIS Online maps don't display well in the phone's browser. I do not come from a **Comp Sci or development background** so if there's already an existing interface that makes it easily to display for mobile I'd prefer that rather than messing with code. The map would contain point features with hopefully a pop-up containing a picture and a short description of the property. I'm pretty computer savvy but I'd like this to be as painless as possible. Are there any other interfaces that exist that make this easier? **edit: The mobile map will be only for viewing purposing, no editing or capture** I'm open to looking at other platforms outside of ArcGIS as well. As an added bonus (but not mandatory), I've been looking into Story Maps and plan on using one in the future, specifically the Short List. However, I don't see this working out for a mobile device as its very picture heavy and would crowd the display. If anyone has any other alternative ideas I'm open to them! Thanks **edit** I've been looking at Leaflet and Mapbox as a possible solution for mobile use - there is plenty of documentation - but still a lot of code! **edit2** For now we've gained access to HTML5 viewer for geocortex and will be using it for our mobile sites. If there are any other alternatives I haven't already listed above - that are about on-par with ArcGIS Online for ease of use for non-programmers, I would be interested to hear about them.


title: Auto increment a new column using Model Builder Field Calculator in QGIS?

text: I am making a model which requires creating a grid and clipping various layers to it. Using the **Create grid** function does not provide an ID column, instead it creates 2 columns: "longitude" and "latitude". Because of this, I would like the Model Builder to insert a new ID column using either **Field Calculator** or **Advanced Python Field Calculator** but I have no idea what the expression is. I am using QGIS 2.2.0 Valmiera. Thank you very much. NOTE: I have noticed this topic: How can you auto-increment in a Python Script (not Field Calculator)? If using the **Advanced Python Field Calculator** , would the python coding be inserted to the _Global expression_ or to the _Formula_?


title: Spatial Join Attribution for Report Generation

text: I generated a grid index for a very large area of land in order to create a mapbook to keep track of some environmental sites for monitoring. I was able to generate the pages of maps without issue, however, I now wish to tabulate the sites (small polygons) contained within each individual map (grid index cell) in a separate file to append to the map book. I spatially joined the small polygons to the grid index cells so that I could index the sites by page number. When I generate a report (extracting only the PageNumber and Site ID in order to create a checklist for each map page), I cannot get the report to group this data into separate pages. It just generates a running spreadsheet in ascending order. What I am aiming for is to create a document where every page contains only the sites for each respective PageNumber. I thought it would be straight forward but I can't figure it out and don't want to export to an excel file to manually separate 75 pages worth of sites.


title: How to merge dbf files in Python or arcpy?

text: I have over one hundred dbf files to merge. Each of files have two columns: `ID` and another field which differs from one file to another. I want to merge those fields based on `ID` such that each raw corresponds to an unique `ID` containing all merged fields. For example:               file 1:     ID FirstName     1  Alan     2  Bell          file 2:     ID LastName     1  R     2  A          merged file:     ID FirstName LastName     1  Alan      R     2  Bell      A      **Solution:** I have a series of files named as 'rawfile200101', 'rawfiles200102',...,'rawfile201104'. Each of files have two columns: `ID` and another field same as the file name. I want to merge those fields based on `ID` such that each raw corresponds to an unique `ID` containing all merged fields. Different from Stata, the "merge" process is called "Join Field" in ArcGIS. Since I have lots of these files, I will do this in batch. Take the following as an example               import arcpy     from arcpy import env     from arcpy.sa import *     arcpy.CheckOutExtension("Spatial")     env.overwriteOutput = True          yrmn = ['200102', '200103', '200104']          for list in yrmn:         inFeatures = "I:/rawfile200101.dbf"         joinField = "VALUE"         joinTable = "I:/rawfile{}.dbf".format(list)         fieldList = ["rawfile{}".format(list)]          arcpy.JoinField_management (inFeatures, joinField, joinTable, joinField, fieldList)


title: need help to automatize report generation

text: i have shapefile that contains a lot of parcels "parcel has like fields (n° parcelle,proprietaire, cin) see picture : ![enter image description here](http://i.stack.imgur.com/9MXVY.png) and i want for all parcels to generate this plan : take a parcel and generate this. (red text is clarification) ![enter image description here](http://i.stack.imgur.com/311pW.jpg) software or language programming not important, i'm open to any idea or answer that can solve this problem(product presentation format for all parcels with the format in picture 2 ).


title: Improve performance of Ext Grid with GeoExt WMSCapabilitiesStore

text: i have a question concerning the performance of the Ext Grid when loading big datasets from the GeoExt WMSCapabilitiesStore. My application loads about 2500 layer records in an Ext Grid Panel. Once the records are loaded the application works very well, thanks to the Buffered Grid View plugin. But the problem is, that the initial load takes a very long time (about 10 seconds) which is not very responsive. I understand that it takes quite a while to generate layer records from a XML file. But since it is not possible to receive WMS Capabilities as JSON (which is recommented for huge datasets in Ext), maybe there is an other way. I assume that most of the time is used to provide an array from the XML request. It seems not very clever, that this process should be repeated every time someone opens the page, especially because the data is static. Maybe there is some way to create an array with all the layer records once and load it into the grid, instead of transforming the data every time? I'd be very grateful if someone has a hint or an advise. Thank you!


title: QGIS 2.0 Problem with loading a Layer Style file that contains a transparent pixel list

text: I work with QGIS 2.0.1. For my raster file I have created a Layer Style with a transparent pixel list. When I import a rasterfile in a new project and load the layer-style file, the map changed to the correct representation. But the transparent pixel list is empty. **When i want to apply the style settings, the presentation changed back to the original shape style.** * * * I have looked in the style file and the transparent pixel list are stored. Code:               <rasterTransparency>         <singleValuePixelList>           <pixelListEntry min="0" max="0" percentTransparent="100"/>           <pixelListEntry min="1" max="1" percentTransparent="40"/>         </singleValuePixelList>     </rasterTransparency>      * * * Second workflow:   * Open the Properties from the raster file    * Select the the Transparency menu   * Click the + Button (add values button)    * Enter the raster values and adjust the transparency    * Export the list to .txt file After this I load the transparency pixel list textfile into the new Project additonally to the layerstyle . This works. Why is the transparent pixel list stored in the layer style file but you must store the pixel list as a extra text file? What can you do, that you dont have to save 2 files (Layerstyle and transparent pixel list .txt)? Thanks for your help.


title: Is it a good idea to use GRASS Algorithms for Geoprocessing in a Server Environment?

text: I´m designing a web-service to run in a server environment that will do some geoprocessing actions to datasets uploaded by users and send the results back. I´m programming in Python, and so far I´ve been using different GRASS algorithms to do the geoprocessing. Which means, I´m using a lot of code that looks like this:               grass.run_command("v.buffer", input="a", output="b", type="line", distance=c)      Does anyone have experience in using GRASS-algoriths in a server environment? Is it a good idea to use GRASS or should I better be using other libraries (e.g. because they are faster or more stable?). It has to be OpenSource, since we don´t have the budget for ArcGIS Server, etc.