title: ArcFM v/s Network Engineer

text: Which software would you prefer for a telecom GIS project and why? I have experience with both, but I personally prefer Telvent ArcFM over Telcordia Network Engineer because of it's functionality and tighter integration with ESRI products, especially ArcGIS Server. Any other opinions? Edit:   I would like to have your views based on (but not limited to) these factors:   1. Performance   2. Scalability   3. Database Support   4. Customization options/ ease of customization   5. Interoperability with other (GIS) systems   6. Data exchange to/from other non-GIS systems (like Sharepoint, middlewares etc.)   7. Availability of a web-based offering  I am purposefully ignoring cost at this moment as I would like a completely technical aspect for the comparison.


title: Why CQL DWITHIN returns different result than PostGIS DWITHIN?

text: I have 2 tables in PostGIS that I publish to geoserver. both have same SRS (2039). When I do the following query in PostGIS , I get 69 results :                SELECT a.gid,b.gid      FROM Table_A a INNER JOIN Table_B b      ON ST_DWithin(a.geom, b.geom,1);      But when I do the same in CQL , I get error that result contains more than 1000 features. This is my CQL :               http://localhost:8080/geoserver/wms?     LAYERS=a,b&STYLES=&FORMAT=image/png&SERVICE=WMS&VERSION=1.1.1&     REQUEST=GetMap&EXCEPTIONS=application/vnd.ogc.se_inimage&     SRS=EPSG:2039&     BBOX=589042.42768447,4914010.3926913,609134.78143081,4928059.3431623&     WIDTH=512&HEIGHT=358&     CQL_FILTER=     DWITHIN(geom,collectGeometries(queryCollection('b','geom','INCLUDE'),1,meters);     EXCLUDE      The Exclude at the end is because I don't want to see layer b as result. Any explanations on this ? _**_* UPDATE _*_** Layer A and Layer B are published to Geoserver , The dimensions in the request I get by default when I use PREVIEW LAYER (Layer A are lines and B are polygons)


title: Visual edit geographic data with geotools

text: With JMapFrame I have options for to visualize my geographic datas, but I don´t know how to add tools for to edit this data, equals conventional GIS. It´s possible to do this? att, Augusto


title: QGIS startup error (out of range pointer)

text: I downloaded QGIS (for Windows) and installed the software. However, when I start QGIS, I get an error message:               "The Instruction at 0x04c20d55 referenced memory at 0xffffffff. The memory could not be read."      Any ideas how to solve this problem??


title: Drive time analysis using Network Analyst with ArcInfo 9.1

text: What is the most common SETTING for the Layer Properties when creating time- travel polygons? I am in the process of creating 4, 5, 6 minute polygons using Network Analyst. Every time I generate these polygons they have different shapes. I would like to be consistent from one year to another based on the number of streets that are being added to our Center Line Street Network. Using ArcInfo – 9.1 SP #1


title: represent multiple geometry as single feature

text: There is a school. it is on two disjoint areas, one is primary wing and other is secondary. I need to make single feature from the two areas such that, there is only one column in gis table, i.e. only one record. Also is it possible for them to have separate values but still show up as one record.


title: Reset Object ID After Truncating Table In SQL Server

text: I am working on a program to export data from a view containing a geometry type in Sql Server 2008, to an ArcSDE table on another server which also stores data using the geometry type. The source database is not running ArcSDE. This will be a nightly job. Following the documentation at the link below, I am able to generate an Object ID with every new insert: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/in_SQL_Server/006z00000017000000/ My problem (if it is indeed a problem) is before I copy records between the two servers, I first delete all items in target table using a SQL statement:               DELETE FROM target_table_name      However, when I insert the first record after the delete, the ObjectID starts off at the next highest value from before when I deleted the table. For example, lets assume that the highest ObjectID in the table before I deleted it was 1000. After I delete the contents of the table, and start to load new data, the new ObjectID the dbo.i_get_ids stored procedure returns is 1001. Each load has around 100,000 records and since this runs nightly, the ObjectID can become a pretty big value pretty quickly. My question is, is there a way I can reset the value? If I use the following ArcObjects code to truncate the table the ObjectID is reset to zero:               ((ITable)targetTable).DeleteSearchedRows(null);      However, if there is a way to do this without ArcObjects, it would be preferred.


title: Persistent settings in MapInfo programs

text: I have a couple of MapInfo programs that use external EXEs for some functionality (e.g. ogr2ogr for layer conversion and OV2 Tools for satnav data output). In my programs I would like to have a settings dialog where you set the location of these EXEs and then _somehow_ store these so that they are remembered even after Mapinfo has been closed down. What's the best way to do this?


title: Looking for open source component to create PDF from WMS in a web application

text: I want to create a PDF from WMS in a web map application. Please can anyone recommend a robust component? I'd like free open-source under an relatively business-friendly licence (no GPL, BSD/MIT etc. would be ideal, LGPL would be OK).   * Our WMS layers are produced by ncWMS and UNM MapServer. We will also want to overlay some layers from external WMS providers.    * _EDIT_ : We are developing our own JavaScript client based on OpenLayers.    * _EDIT_ : We have some backend code in Java. We would also consider PHP or other platforms to create the PDF. MapFish would be excellent but it is GPL.


title: Is it possible to export a GeoJSON with an array field from a joined layer?

text: I'm using QGIS to 1) load a shapefile, 2) join the shapefile data with a CSV- file and finally 3) exporting the layer as GeoJSON. Hence, I end up with something like this, where the `properties` object contains the original CSV data.               {       "type": "FeatureCollection",       "features": [         {           "type": "Feature",           "geometry": {              "type": "Point",              "coordinates": [102.0, 0.5]            },           "properties": {              "prop0": "value0",              "prop1": "value1",              "prop2": "value2"           }         },...      What I really want to do is to store some of the properties in an array, as such:               "properties": {         "prop0": "value0",         "arr": [             "value1",             "value2"     ]      Is it possible to export GeoJSON in such a format using QGIS? If not, does anyone have any other ideas as to how to edit GeoJSON data like that?


title: binding nested arrays as geojson popups in leaflet

text:
I am struggling with geojson popups in leaflet.js I have nested arrays as properties for my features - an example would be counts of fish caught at a location (not my actual project data!)               {"type":"FeatureCollection",     "features":[{"geometry":{"type":"Point",     "coordinates":[-122.679416,     45.516632]},     "type":"Feature",     "properties":{"fisher_id":"401",     "river":"Willamette",     "fishcaught":[{"type":"Salmon",     "count":5},     {"type":"Sturgen",     "count":1},     {"type":"Bass",     "count":2}]}},     {"geometry":{"type":"Point",     "coordinates":[-121.318584,     44.050165]},     "type":"Feature",     "properties":{"fisher_id":"401",     "river":"Deschutes",     "fishcaught":[{"type":"Trout",     "count":10},     {"type":"Whitefish",     "count":1}]}},     {"geometry":{"type":"Point",     "coordinates":[-124.102962,     43.983057]},     "type":"Feature",     "properties":{"fisher_id":"401",     "river":"Siuslaw",     "fishcaught":[{"type":"Steelhead",     "count":2},     {"type":"Salmon",     "count":6}]}},     {"geometry":{"type":"Point",     "coordinates":[-121.776289,     42.220382]},     "type":"Feature",     "properties":{"fisher_id":"401",     "river":"Klamath",     "fishcaught":[{"type":"Trout",     "count":3},     {"type":"Steelhead",     "count":2},     {"type":"Salmon",     "count":3}]}}]}      I am trying to put a popup on each feature (onEachFeature) that states river and a count of each type caught. I'm binding using a function like this               function onEachFeature(feature, layer) {             if (feature.properties.fishcaught) {                 var popupcontent;                 popupcontent = "River: "+feature.properties.river                 for (i=0;i<feature.properties.fishcaught.length;i=i+1){                     var fishstring = "
"+feature.properties.fishcaught[i].type+": "                     +feature.properties.fishcaught[i].count                     popupcontent = popupcontent+fishstring;                 }                 layer.bindPopup(popupcontent);             }         }          geojson = L.geoJson(inData, {         onEachFeature: onEachFeature     }).addTo(map);      The problem is that the onEachFeature is caught in an endless loop. If I remove the for loop from the function, everything works fine. Any thoughts on this?


title: Time-aware layers WCS Service

text: I am currently developing an application that requires to display raster images over a period of time (I'm looking how Hurrican Sandy progressed). While investigating, people have suggested that I could use a WCS service and view it on a viewer. We are currently using Google Maps API. My question is, how can I publish these rasters in a WCS service (or relevant format) that can allow a viewer to loop through all images. Is there a standard practice everyone follows to accomplish this? I have knowledge using ArcGIS 9.3 Desktop/Server, ArcGIS Flex Viewer, and QGIS. I appreciate any information regarding this topic or tutorials (haven't found any good ones yet). The images are geotiffs. Thanks in advance!


title: What is the future of ArcGIS web development ? Is it JavaScript and HTML5?

text: I want to know about the future of ArcGIS web development; is it HTML5? I used to work with ArcObjects and .NET Web ADF so what is the best technology I must follow, especially when I know that Silverlight version updates have been stopped? So I think that the best technology that must follow .NET in the future is JavaScript with HTML5. Is it true?


title: How to generalize a point feature class

text: I have a point featureclass with a point for every 25cm - in Total 400,000 Points. I would like to reduce this resolution to one point every Meter. How could I instruct ArcGIS Advanced(ArcInfo) to do this? I can´t see a thinning tool or generalization tool for Point feature classes and don´t want to create a raster from the data yet to resample.


title: MapFish PDF: Geoserver WMS layers not generated as tree structure

text: I am using the mapfish printing module in my web map application.In that everything is working fine but layers not generated as tree structure in output PDF file and "ODC_Graphics_Data" wms layer not printed. I have used below code to generate pdf. System: Windows XP / Tomcat7 / Geoserver 2.4.0/printing-2.1-SNAPSHOT.jar Kindly help me to fix this. Thanks in advance.                 http://localhost:8080/geoserver/pdf/print.pdf?spec={       "layout":"A3 portrait",       "srs":"EPSG:404000",       "units":"m",       "dpi":300,       "outputFilename": "ggg",       "layers":[         {             "baseURL":"http://xxx.xxx.xx.xx:8080/geoserver/cite/wms",             "opacity":1,             "singleTile":false,             "type":"WMS",             "layers":["Oxford_Map"],             "format":"image/png",             "styles":[],             "overview":true         },         {             "baseURL":"http://xxx.xxx.xx.xx:8080/geoserver/cite/wms",             "opacity":1,             "singleTile":true,             "customParams":{"CQL_FILTER": "User_id = 234"},             "type":"WMS",             "layers":["ODC_Graphics_Data"],             "format":"image/png",             "styles":[],             "overview":false         },         {                "baseURL":"http://xxx.xxx.xx.xx:8080/geoserver/cite/wms",             "opacity":1,             "singleTile":true,             "type":"WMS",             "layers":["Plymouth_Buildings"],             "format":"image/png",             "styles":[],             "overview":false         },         {             "baseURL":"http://xxx.xxx.xx.xx:8080/geoserver/cite/wms",             "opacity":1,             "singleTile":true,             "type":"WMS",             "layers":["Plymouth_Hard_Surface_Block_Paving"],             "format":"image/png",             "styles":[],             "overview":false         },         {             "baseURL":"http://localhost:8080/geoserver/cite/wms",             "opacity":1,             "singleTile":true,             "type":"WMS",             "layers":["Plymouth_QUERIES"],             "format":"image/png",             "styles":[],             "overview":false         },         {             "baseURL":"http://xxx.xxx.xx.xx:8080/geoserver/cite/wms",             "opacity":1,             "singleTile":true,             "type":"WMS",             "layers":["Plymouth_Street_Lights"],             "format":"image/png","styles":[],             "overview":false         }     ],      "pages":[         {             "center":[248744.93342196,55350.945],             "mapTitle":"",             "company":"testuser2_portrait",             "comment":"",             "scale":439,             "rotation":0         }     ] }      Please find screenshot for output PDF: ![enter image description here](http://i.stack.imgur.com/5eMOD.png) ![enter image description here](http://i.stack.imgur.com/f88A3.png)


title: How to store the CHANGES to a map database

text: I'm looking into providing a map amendment service for a particular map series. The service will help map owners look into the changes made to the operational features of the map, after the map was published. There are three basic types of changes: 1) An object has been added. 2) An object has been deleted. 3) An object has changed, either an attribute or its geometry. For object changes, the user should get both the old and the new value of the attribute or geometry. Now, I'm looking into a clever way to implement this. The basic idea is that the database administrator edits the map data in an interface, just as he would edit a normal map, and the rest is taken care of behind the scenes. As I see it, there are two basic approaches: 1) Store the changes in a separate table, and keep the main database current. 2) Timestamp the records in the main database. In case of changes, I cannot change a record, but timestamp the old one as expired and add a new urrent one. My main concern is not performance (not huge amounts of data), but ease of implementation. Is there a way I can make the DBMS (say, PostGIS) handle this automatically (for example by using triggers), Or can I build my own interface on top af an exisiting GIS (say, QGIS or ArcGIS) to handle the storage of records in the various tables as the database administrator makes changes to the map. Preferrably, I'd like a Web-interface. Any Ideas anyone?


title: using Map Service or Feature Service for Feature Layer

text: Both, Dynamic Map Services and Feature Services allow for querying data (through FeatureLayer) and bringing FeatureSets to the browser. Which of the above would yield better performance at querying attributes and geometries of large datasets?


title: Extract Street, City, And State/Region from Planet.osm with Osmosis

text: How can I extract street, city, and state/region data from Osmosis, and nothing else? I know using Xapi I could query under ways and use name=* to get streets, but I need to do this with Osmosis, in addition to city and state. My end goal is to query with xapi to get streets cities and regions in a given bounding box


title: Rendering custom OpenStreetMaps style (land=white, water=black): do I need a dedicated computer just to do this?

text: EDIT: My goal is NOT to use OpenStreetMap or cloudmade as a web map. I described exactly what I want ("a full, rasterized set of tiles for the entire world at zoom=10 (around 68 Gigapixels, water=black, land=white, no labels"). I only meant to use the cloudmade.com link as an example to illustrate what I'm shooting for. (As it is, the cloudmade style is only 90% of the way there.) I need the tiles offline. I need the image (68 Gpix) I described as a mask in a larger raster map task for a custom map I'm building using NASA and USGS data. * * * **As an illustration** I made a custom style that sort of shows what I want:   * http://maps.cloudmade.com/?styleId=22688# (No labels, no roads, no features of any kind except all water features (rivers, lakes, oceans, you name it) as black with land being solid white, and the nice Mapnik antialiasing of values in between.) * * * **What I really want** I would like a full, rasterized set of tiles for the entire world at zoom=10 (around 68 Gigapixels). I considered downloading and installing Planet.osm (11 GB compressed) and running Mapnik (for days/weeks?) to get this data. However, Planet.osm looks like it will take 100GB or maybe even 1TB once the database and index are built. I don't have that kind of disk space on my laptop and since I don't need the full dataset, is there a smaller subset of the data I could download (it says here that the polygon data is only 700 MB)? Actually from what I can tell, the data has doubled in size in the last year so I would need a beefy machine to deal with this. Is there an easier, direct way to get these tiles (level 10, black and white, just land on water) batched, or do I basically need a dedicated computer to do this?


title: DrawPolyLine not drawing completely on the right layer

text: I zoom into a map, then do a screenDisp.DrawPolyline(polyline); Now the line is in the correct place, but the map that is the background of the polyline object is actually the initial map and not the zoomed in map. Here's a screenshot: http://i11.photobucket.com/albums/a173/iteration_x/LineBehavior.jpg What you are seeing is a really zoomed in map of minnesota with most of the layers turned off, and that interior picture is part of the West Coast from the initial map load.


title: How do I disable ArcGIS Server web adf calls when zooming with a mouse wheel?

text: When zooming using the mouse wheel on a web adf based map, every time the map wheel clicks, requests are sent to ArcGIS Server for maps etc. (as evident in Fiddler) Is there any way to harvest an event that would disable the requests until you got to zoom level you want?


title: Issues with geoprocessing with .NET

text: There are some nice features in ArcToolbox we can use, but for some reason, this is NOT working properly. It doesn't even throw me an error. My software is running inside ArcMap, so no need to AoInitialize again, corret?                   public void Execute()         {             InitializeProduct();             try             {                 Geoprocessor gp = new Geoprocessor();                 gp.OverwriteOutput = true;                      FeatureToPoint featureToPoint = new FeatureToPoint();                      string outputPathName = CurrentWorkspace.PathName + "\\teste_centroide";                      featureToPoint.in_features = InputFeatureClass;                 featureToPoint.out_feature_class = outputPathName;                 featureToPoint.point_location = "INSIDE";                      IGeoProcessorResult result = (IGeoProcessorResult)gp.Execute(featureToPoint, null);                      if (result == null)                 {                     for (int i = 0; i <= gp.MessageCount - 1; i++)                     {                         Console.WriteLine(gp.GetMessage(i));                     }                 }                      IGPUtilities gpUtils = new GPUtilitiesClass();                 this.OutputFeatureClass = gpUtils.OpenFeatureClassFromString(outputPathName);             }             catch (Exception ex)             {                 MessageBox.Show(ex.Message + "\r\n");             }      This is a code example I'm having here. I generated the DataManagement tools assembly, but I could not find the file to sign it. This code just gives me an error. is it because of the signing? I've tried the other way too, using IVariantArray and calling from the tool name, without sucess. Is it just me or...? Can anyone point me a "nicer" solution? I need to run several process that are already built in ArcToolbox that I really don't want to duplicate. Any tips or pointers? Thanks guys!


title: unique value symbolize layers from a lookup table

text: I need to create a large number of Arcmap layer files (v9.3 or 10), perhaps a hundred, from lookup tables. I don't have to set the symbol patch itself, that can be random, but I do need the value and label attributes filled out. The lookup tables look like:               "FC","VALUE_FIELD","VALUE","LABEL"     "Wooded_area","CODE","1240012","Extracted"     "Wooded_area","CODE","1240022","Interpreted"     "Wooded_area","CODE","1240032","FCS-EOSD"     "Wooded_area","CODE","1240042","Land Cover, Circa 2000 - Vector"      And the result would resemble (except the symbols themselves): ![screen shot of layer properties](http://i.stack.imgur.com/fkerZ.png) I'd prefer a python based solution, via arcobjects+python if required, C# next, but anything that works, even partially or only conceptually, is welcome.


title: Quickly create contour polygons from cost points using OSS for web app

text: This is related to a previous question I've asked: Create filled contour map with GRASS I've created a point layer using OpenStreetMap data and pgRouting's driving_distance() function, where the points contain a cost attribute which contains the walking distance from a school. I've used QGIS and the contour plugin to create what I want: ![walking distance isolines](http://i.stack.imgur.com/ieyuV.jpg) **However I'd now like to do this programmatically for a web app, and in addition result in polygons rather than a raster.** I've got something working using GRASS, however it is quite slow, around 30 secs on my machine versus less than 1 sec using QGIS and the contour plugin. Here the short version (or see Python version):               # read file and set resolution and bounds     v.in.ogr dsn=cost_points.shp output=cost_points     g.region vect=cost_points rows=256 cols=256     # interpolate points and make bands of 0.4km increments     v.surf.rst input=cost_points elev=interp zcolumn=cost     r.mapcalc res='int(interp/0.4)'     # convert to polygons and write out file     r.to.vect -s input=result output=polys feature=area     v.out.ogr -c input=polys dsn=polys.shp type=area      **Are there other GRASS functions I should use?** r.surf.rst takes the majority of the time, so perhaps it is too robust for my needs, but even then the results aren't smooth enough for my liking. **Are there other open source libraries I should be considering for this?** If you'd like to know more about my overall goals, see this page.


title: OpenLayers + Geoserver: how to dynamically retrieve the "geometryName" of a feature layer in GeoServer

text: I am trying to get the list of feature layers on GeoServer and then user can add layers into map at anytime they want. my question is: when using `OpenLayers.Protocol.WFS` loading a layer, is "geometryName" required? since from the response of "WFSGetCapabilities" request to Geoserver, for each featureType, there are only "name, title, srsName, and featureNS". If “geometryName” is required, I have to retrive this from "DescribeFeatureType" request by enumerating the feature's properties . however, these two requests are asynchronous, I am not sure how to load this layer without user's interaction input. anybody helps? appreciated!!


title: Network Analyst: Location Allocation for Model Builder

text: I am trying to create a model to answer a question using a network analysis. Background:   * 10 possible fire station locations (shapefile)   * Calls for demand points (shapefile)   * Impedance is 5 minutes from facility to demand Product: I need to make a model to select 1-3 stations from the list of 10 for the best locations to meet as many calls as possible with a 5 minute drive from the station. Below is an image of what I have so far, but when I try and run it I get an Error saying:   * Error 000735: Input Network Analysis Layer: Value is required   * WARNING 030036: A layer with this name exists Not sure how to fix these two items. ![current model](http://i.stack.imgur.com/u9pqI.png)


title: Data Security for Geoserver - Users errantly seeing Workspaces not assigned to them

text: I'm having trouble with Data Security in my Geoserver 2.5 application. Let me preface my situation. I am using Geoexplorer composer which is making GeoServer API calls through ows and wms. I have multiple clients which should not be allowed to see each other's wms services. The implementation of my data security rules are scatter shot when ows and wms are used through GeoExplorer. CompanyAUser can see all layers in the CompanyAWS and 1 layer from CompanyB's ws. CompanyBUser can see 10 out of the 60 layers in CompanyA's ws. A further twist to this problem is that if I sign in as each company's user in Geoserver's webapp and I navigate to Layer Preview each user can only see the layers which I have specified that they are allowed to see in Data Security. I have made some edits to the Geoexplorer Composer code but nothing major. I don't think that this should effect the security in any way since security is implemented by the Geoserver application. I have included my security parameters below.               Groups:     group, role     CompanyA, COMPANY_A     CompanyB, COMPANY_B     CompanyC, COMPANY_C          Users:     user, group     CompanyAUser, CompanyA     CompanyBUser, CompanyB     CompanyCUser, CompanyC          Data rules:     Roles, rule path     Company_A, COMPANYAWS.*.r     Company_B, COMPANYBWS.*.r     Company_C, COMPANYCWS.*.r          Service rules:     *.*     ADMIN,COMPANY_A,COMPANY_B,COMPANY_C          Roles are independent and do not have parent/child relationships.      Are my roles, groups, data rules and service rules set correctly? EDIT1 I've noticed that the only layers that are errantly shared to users belonging to different company groups are 'Group Layers'. This leads me to my next question. What causes a regular layer to be hidden and a group layer to be shown to roles that do not have permission to do so?


title: Graphical Modeler to rotate points based on line layer azimuth

text: I am trying to create a model using the QGIS 2.0 Processing Toolbox graphical model to replicate the process described in this thread: How to horizontally align irregularly shaped point symbol to vector layer in QGIS? ![enter image description here](http://i.stack.imgur.com/webgJ.png) However, I get an error when trying to execute the Distance to Nearest Hub algorithm: Invalid name attribute: bearing. When I go through the process manually, using the same data, this doesn't happen. UPDATE: running the field calculator using the python console doesn't populate my vector layer's bearing field as I expected it to. A field called "bearing" is created in a new output layer (i've corrected to link this to distance to nearest hub instead of how it is currently), but the field is blank. I ran the script like so after import processing:               processing.runalg("qgis:fieldcalculator", resLayer, 'bearing', 0, '(atan((xat(-1)-xat(0))/(yat(-1)-yat(0)))) * 180/3.14159 + (180 *(((yat(-1)-yat(0)) < 0) + (((xat(-1)-xat(0)) < 0 AND (yat(-1) - yat(0)) >0)*2)))', 'bearingdata.shp')


title: ArcGIS JavaScript manually render a layer with supplied JSON

text: Given the following JSON which is retrieved outside the ArcGIS/dojo framework, how would I go about rendering this layer? Thanks in advance!                   {             "currentVersion": 10.01,             "serviceDescription": "",             "mapName": "Layers",             "description": "",             "copyrightText": "",             "layers": [                 {                     "id": 0,                     "name": "Transportation Stations-unexpired",                     "parentLayerId": -1,                     "defaultVisibility": true,                     "subLayerIds": null,                     "minScale": 400000,                     "maxScale": 0                 }             ],             "tables": [],             "spatialReference": {                 "wkid": 102100             },             "singleFusedMapCache": false,             "initialExtent": {                 "xmin": -10398329.13791846,                 "ymin": 2687052.938617364,                 "xmax": -7541451.81403804,                 "ymax": 5379667.868452799,                 "spatialReference": {                     "wkid": 102100                 }             },             "fullExtent": {                 "xmin": -10152568.449939176,                 "ymin": 2795685.894150028,                 "xmax": -7787215.334097848,                 "ymax": 5828927.975739868,                 "spatialReference": {                     "wkid": 102100                 }             },             "units": "esriMeters",             "supportedImageFormatTypes": "PNG24,PNG,JPG,DIB,TIFF,EMF,PS,PDF,GIF,SVG,SVGZ,AI,BMP",             "documentInfo": {                 "Title": "test.mxd",                 "Author": "test",                 "Comments": "",                 "Subject": "",                 "Category": "",                 "Keywords": "",                 "Credits": ""             },             "capabilities": "Map,Query,Data"         }


title: Is it possible to install 9.3 Desktop Administrator as a standalone program?

text: I have ArcServer 9.3 installed on a server that I use for running scheduled Python geoprocessing scripts. I have several running, and now I'm trying to schedule one to rec/post, delete and recreate our database versions. The issue that I'm having is this line of code               gp.SetProduct("ArcEditor")      Apparently to run the versioning tools, I need to have valid license of ArcEditor or higher availible to the server. I would really rather not install ArcGIS desktop on my server box, but I'm stuck for a solution. My main question is, is it possible to install the 9.3 Desktop Adminstrator as a standalone program, then point it at my license manager so that the script can find a license file? If not, that's fine, any ideas are appreciated.


title: Use Python lists in Field Calculator?

text: I couldn't find this question, although I am sure it has been asked. How could I use a python list to condense this script? This code block currently works, but is clearly not very efficient.               def Reclass(test):     if test == '2 7':         return 'Padmounted Single Phase'     elif test == '2 4':         return 'Padmounted Three Phase'     elif test == '2 2':         return 'Padmounted Three Phase'     elif test == '2 1':         return 'Padmounted Three Phase'     elif test == '1 4':         return 'OH Single Phase'     elif test == '1 2':         return 'OH Single Phase'     elif test == '1 1':         return 'OH Single Phase'     elif test == '1 5':         return 'OH Second Phase'     elif test == '1 3':         return 'OH Second Phase'     elif test == '1 7':         return 'OH Third Phase'      I was attempting to use (as an example of a list):               elif test == ['1 1', '1 2', '1 3']      but the code was being skipped over and not working.


title: Persisted member in an application extension does not save to stream

text: I have an ArcMap extension in which I'm trying to persist a string value. When the document is saved, the public member m_LayerName in the application extension class is not being persisted. Either it's not being written to the stream or it's not being loaded correctly once the document is reopened. I'm not sure what I'm doing wrong. The Save and Load events are definitely being executed appropriately, but no data is being saved. This is my first attempt at persisting data with an extension so I'm missing something. Below is the relevant code.               Public m_LayerName As String          ...          Public ReadOnly Property ID() As ESRI.ArcGIS.esriSystem.UID Implements ESRI.ArcGIS.esriSystem.IPersistVariant.ID         Get             Dim typeID As UID = New UIDClass()             typeID.Value = Me.GetType().GUID.ToString("B")             Return typeID         End Get     End Property          Public Sub Load(ByVal Stream As ESRI.ArcGIS.esriSystem.IVariantStream) Implements ESRI.ArcGIS.esriSystem.IPersistVariant.Load         m_LayerName = CStr(Stream.Read())         Marshal.ReleaseComObject(Stream)     End Sub          Public Sub Save(ByVal Stream As ESRI.ArcGIS.esriSystem.IVariantStream) Implements ESRI.ArcGIS.esriSystem.IPersistVariant.Save         Stream.Write(m_LayerName)         Marshal.ReleaseComObject(Stream)     End Sub


title: Identifying user triggering ArcGIS Server invalid username or password error

text: Can anyone tell me if there is a way to identify the culprit of an "Invalid Username or Password" error in ArcGIS Server logs? For source, it just says 'Admin' for source, and nothing for User. Is there a .txt log stored somewhere that would tell me which username is causing the bad login?


title: Database connection failed

text: I am trying to connect from Arc Cataloge (version 10.1) to postgresql 9.2. But every time I try to make "Database Connection" I am getting this error:               Failed to connect to the specified server.     Underlying DBMS error[could not connect to server: Connection timed out        Is the server running on host"x.y.z.w" and accepting TCP/IP connections on port 5432?      Postgresql is installed on Centos server (PORT 5432) and I can access it using pgadmin. Using `psql --h X.Y.Z.W -U user -d database` on another server I am getting the same error. here is the content of pg_hba.conf:                local   all             all                                md5      # IPv4 local connections:      host    all             all            X.Y.Z.W/32          md5      host    all             all             0.0.0.0/0          md5      host    all             all             127.0.0.1/32       md5      # IPv6 local connections:      host    all             all             ::1/128            md5      Also in postgresql.conf:               Listen_addresses = '*'     port=5432      in etc/sysconfig/iptables:                *filter       4 :INPUT ACCEPT [0:0]       5 :FORWARD ACCEPT [0:0]       6 :OUTPUT ACCEPT [0:0]       7 :RH-Firewall-1-INPUT - [0:0]       8 -A INPUT -j RH-Firewall-1-INPUT       9 -A FORWARD -j RH-Firewall-1-INPUT      10 -A RH-Firewall-1-INPUT -i lo -j ACCEPT      11 -A RH-Firewall-1-INPUT -p icmp --icmp-type any -j ACCEPT      12 -A RH-Firewall-1-INPUT -p 50 -j ACCEPT      13 -A RH-Firewall-1-INPUT -p 51 -j ACCEPT      14 -A RH-Firewall-1-INPUT -p udp --dport 5353 -d 224.0.0.251 -j ACCEPT      15 -A RH-Firewall-1-INPUT -p udp -m udp --dport 631 -j ACCEPT      16 -A RH-Firewall-1-INPUT -p tcp -m tcp --dport 631 -j ACCEPT      17 -A RH-Firewall-1-INPUT -m state --state ESTABLISHED,RELATED -j ACCEPT      18 -A RH-Firewall-1-INPUT -m state --state NEW -m tcp -p tcp --dport 22 -j ACCEPT      19 -A INPUT -m state --state NEW -m tcp -p tcp --dport 5432 -j ACCEPT      20 -A RH-Firewall-1-INPUT -m state --state NEW -m tcp -p tcp --dport 80 -j ACCEPT      21 -A RH-Firewall-1-INPUT -j REJECT --reject-with icmp-host-prohibited      22 -A RH-Firewall-1-INPUT -m state --state NEW -m tcp -p tcp --dport 5432 -j ACCEPT      23 COMMIT      I also put "postgresql-client-libraries" in /bin directory of ArcGIS for Desktop. I appreciate any help.


title: Display Leaflet Markers with UTM coordinates?

text:
is it possible to display markers with UTM coordinates (zone 32) with leaflet and OSM? How can I do this? If I use UTM Coordinates like latlng, the position is totally wrong:                       L.marker([5698015.694, 312400.093]).addTo(map)             .bindPopup("Hello world!
I am a popup.").openPopup();      Did I have to translate them first, or is it possible to display UTM coordinates directly with leaflet?


title: Merging more than two shapefiles in QGIS

text: I want to merge 9 shapefiles (all polygons) into one file. The problem is, when I try to merge more than two files at once, I seem to loose attributes on the way? Is there any other way of merging shapefiles than using the 'Merge shapefiles to one'-tool? In this case it's only nine files, so it's possible to merge whem 2 and 2, but I want to do the same thing with other files as well, and then merge about 20 files to one. I'm using QGIS 2.0.1. Thanks! EDIT: I have data of 4 different kinds of natural reserves (polygons). Each reserve has data seperated into three different shapefiles, where each shapefile is cooected to a powerline. The three powerlines are parallell at first and then turn in different directions. Each kind of natural reserve contains slightly different information, but they all contain the columns 'Powerline' (powerline identity) and 'Type' (type of natural reserve). I want my output to contains all of the polygons and all attribute data. What happened was that information like 'Powerline' was lost when I tried to merge all files at once.


title: How to convert a bil image into grid in ArcGIS?

text: GIS beginner here. I've downloaded **RFE data from FEWSNET**. It's a compressed file which includes .bil files. Now i want to convert this ".bil image" into a grid to use the precipitation data in an area. I'm using ArcGIS. How do i do it?


title: Maximum scale dependent visibility changes from 1:0 to -2147483648:1

text: QGIS 2.0.1 Dufour. When i open the layer properties and the scale dependent visibility is unchecked and when i check it, the maximum (inclusive) shows 1:0. I click Apply, OK, reopen layer properties and it says -2147483648:1. Everytime i change it to 1:0, apply and reopen the layer properties, the maximum (inclusive) gives me the value -2147483648:1. Any ideas why?


title: Multiprocessing issues with ArcPy

text: I am using ArcGIS 'watershed' routine in a script which makes use of multiprocessing. This works fine, but during execution I get the message:               Unable to remove directory.  Possible causes:     1- Not owner of the directory     2- Another person or application is accessing this directory      EDIT 27/02/14: I've incorporated the suggestions of using `arcpy.Exists` (with a very rudimentary way of checking if something didnt exist) and writing all results to disk:               import arcpy     from arcpy import sa     arcpy.CheckOutExtension("Spatial")          def watershed(pnts, flowdir, flowacc):         direc = tempfile.mkdtemp(dir = "C:\\temp")   #Create a separate directory for files to be written to.         arcpy.env.scratchWorkspace = direc         arcpy.env.workspace = direc              res = []         for i, p in enumerate(pnts)             pnt = arcpy.PointGeometry(arcpy.Point(p.x, p.y, ID = i)) # Convert the Shapely point to arcpy point             pourpt = sa.SnapPourPoint(pnt, flowacc, 10000) #Snap point to high flow acc cell             if arcpy.Exists(pourpt):                 ws = sa.Watershed(flowdir, pourpt) # Calculate watershed                 if arcpy.Exists(ws):                    out = os.path.join(direc, "poly_%i"%i)                    poly = arcpy.RasterToPolygon_conversion(ws, out) #Convert to polygon                    res.append(poly[0]) # Put the polygon in the results list for this set of points                 else polylist = "NoWS"             else polylist = "NoPourPt"              return res      The parallel function remains the same:               def watershed_pll(mod, proc=6):         """         Calculate the watershed for each station point using parallel processing          All relevant data is held in the mod object         """         pool = Pool(processes = proc)              # Iterate over each feature of the drainage path and submit the list of pour points to the worker function         jobs = []         for key, val in mod.stations.geometry.iteritems():             jobs.append((key, pool.apply_async(watershed, (val.points, key, mod.flowdir,                                                mod.flowacc, rdict))))              pool.close()         pool.join()         return jobs      Sometimes it runs perfectly, other times I get one of various errors such as the initial one above, or also:               FATAL ERROR (INFADI)     MISSING DIRECTORY      and:               ERROR 010088: Invalid input geodataset (Layer, Tin, etc.).]      and:                ERROR 010050: Cell size is not set.      Which is odd as I explicitly set the cell size within the function... When it works and when it doesnt seems to be hit and miss/random. Such errors lead me to believe that arcpy is still trying to delete stuff behind the scenes. Is there any way to explicitly prevent it from doing so? * * * I wonder if my problem could be related to one of two things:   1. Environments/folders getting mixed up. The input data is passed as complete filepaths to some raster datasets which are outside of the scratchWorkspace (which is set locally for each process) where intermediate/output data is created. However, I have noticed that Arc may make folders (typically and 'info' folder) in the directories of the input data. Why is that? Can it be prevented, or can I prevent Arc from then trying to delete it?   2. Are there any potential problems with accessing the input raster data sets at the same time? I.e each process will be attempting to open and read from the flow direction and flow accumulation rasters which are passed into the function.


title: "show only visible features" doesn't seem to work

text: I'm trying (unsuccesfully) to place the attribute table next to a map containing the same in print composer. Let's say for example my map in the primary window contains 6 districts and my map in the print composer only shows 1 district. If I unselect the "show only visible features" box then my table shows all entries (e.g. from 6 districts). However when I select the "show only visible features" it shows **only** the title row of the attribute table: **it does not show any of the data entries**. What am I doing wrong?


title: PostGIS: incorrect interpretation of a polygon that intersects the 180th meridian

text: I try to find out all geopoints which intersect the polygon set as a parameter. The problem is when I pass polygon that roughly covers an area of Bering Strait (nearby 180 longitude): ![enter image description here](http://i.stack.imgur.com/hdpk3.jpg) So I use the query:               SELECT ST_AsText(l.geo_point)     FROM "lightnings" "l"     WHERE (ST_Intersects(ST_GeomFromText('Polygon((132.0 40.0, -148.0 40.0, -148.0 -8.0, 132.0 -8.0, 132.0 40.0) )', 4326), geo_point));      As you can see, vertexes are set in a correct order, clockwise, from the North-West. But the result covers outside area and including whole other world. For example, in result:               POINT(75.5637 40.0434)      The problem doesn't touch the 0th meridian. Elementary test:               SELECT ST_Area(ST_GeomFromText('Polygon((0.0 60.0, 10.0 60.0, 10.0 40.0, 0.0 40.0, 0.0 60.0) )', 4326))     UNION ALL     SELECT ST_Area(ST_GeomFromText('Polygon((-5.0 60.0, 5.0 60.0, 5.0 40.0, -5.0 40.0, -5.0 60.0) )', 4326))     UNION ALL     SELECT ST_Area(ST_GeomFromText('Polygon((175.0 60.0, -175.0 60.0, -175.0 40.0, 175.0 40.0, 175.0 60.0) )', 4326))      gives the result:               200     200     7000      Is there any simple trick to force PostGIS understanding me? I don't like an idea to divide the polygon...


title: ArcGIS Python Error Raster Domain

text: I´m always getting an Error with my Python-Code. Even though it is almost the same as on the ArcGIS Help Website. The Error occurs in line 19: > Traceback (most recent call last): > > File "C:_geodata\check.py", line 19, in > > arcpy.RasterDomain_3d(rasterList,out_file,geometry) > > File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\ddd.py", line > 998, in RasterDomain > > raise e RuntimeError: Object: Error in executing tool > > Failed to execute (IterateRasterDomain). My Code looks like this:               import arcpy     from arcpy import env          arcpy.CheckOutExtension("3D")          env.workspace = arcpy.GetParameterAsText(0)     rasterList = arcpy.ListRasters()          if rasterList:         for raster in rasterList:             geometry = "POLYGON"             out_file = "domain_" + raster[:-4] + ".shp"             arcpy.RasterDomain_3d(rasterList,out_file,geometry)         print "Finished"     else:         print "No Data available"


title: Correcting gaps and overlapping polygon (parcels) geometry

text: I have come across a resource that we rely upon almost daily that is in terrible need of repair. It is not our own product and the supplier is generally not interested in fixing the problem. We are using a polygon shapefile which contains upwards of 2000 parcel boundaries. The parcels themselves are not parcels in the general sense of the term, but nevertheless, it describes the data quite effectively. No 'parcel' shall overlap or produce gaps between parcels. However, the data quality is terribly low and it almost seems like the parcels have been drawn by hand with snapping disabled! See picture below: ![enter image description here](http://i.stack.imgur.com/YvPQv.png) In some cases, it is much worse than in the photo, or the gaps/overlapping is so small that they are unidentifiable until after running spatial procedures (such as intersections). How would one go about fixing these kinds of issues? I have attached the shapefile which represents the 8 polygons in the photo. These 8 polygons show both my gap and overlapping issue. Can the process be automated to fix my 2000+ parcels? http://temp-share.com/show/dPf3mpihW Note: I have explored the "fix geometry" tool in ArcMap, but this only corrects the overlapping polygons. Albeit the overlapping geometry is the worse of the two problems, I'd like to fix both if possible!


title: How to label only some places instead of all?

text: How do I turn of some place names, or rather turn on some places names but not the rest.   How do I turn of or hide that?   [s] *  ~~I also get a red cross under the place names which I dont like to see. How do I turn of or hide that?~~ [/s]   The place names is in a line layer (using labeling), and there are also many numbers I dont need to see.   Some of them are Height above mean sea level (AMSL), but not all.   QGIS Lisboa versjon 1.8.0 ![enter image description here](http://i.stack.imgur.com/LTgXC.jpg)


title: Trying to use a QGIS plugin - getting classFactory error

text: I've had to create a GIS for the company I work for, choosing QGIS over ArcGIS due to price (company doesn't want to buy another copy of ArcMap). In any case, it's not a HUGE GIS but it would be good to be able to search for features quickly and easily - just by typing in. It seems QGIS doesn't support this as easily as the users want to use it, so I've tried downloading the Quickfinder plugin. In any case, after eventually finding out how to install it (as it will be offline when in use, so had to download the zip and extract to the correct place), QGIS has found it but can't install it, saying it's broken and can't call it's classFactory() method. The full code (with my project name removed due to confidentiality issues) is below:               Couldn't load plugin quickfinder due an error when calling its classFactory() method          Traceback (most recent call last):     File "C:/PROGRA~1/QGISDU~1/apps/qgis/./python\qgis\utils.py", line 204, in startPlugin     plugins[packageName] = package.classFactory(iface)     AttributeError: 'module' object has no attribute 'classFactory'          Python version:     2.7.4 (default, Apr  6 2013, 19:54:46) [MSC v.1500 32 bit (Intel)]               QGIS version:     2.0.1-Dufour Dufour, ebebdf3          Python path: ['C:/PROGRA~1/QGISDU~1/apps/qgis/./python/plugins\\processing',     'C:/PROGRA~1/QGISDU~1/apps/qgis/./python',     'C:\\Documents and Settings\\my.name/.qgis2/python',     'C:\\Documents and Settings\\my.name/.qgis2/python/plugins',     'C:/PROGRA~1/QGISDU~1/apps/qgis/./python/plugins',     'C:\\PROGRA~1\\QGISDU~1\\bin\\python27.zip',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27\\DLLs',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27\\lib',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27\\lib\\plat-win',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27\\lib\\lib-tk',     'C:\\PROGRA~1\\QGISDU~1\\bin',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27\\lib\\site-packages',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27\\lib\\site-packages\\PIL',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27\\lib\\site-packages\\win32',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27\\lib\\site-packages\\win32\\lib',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27\\lib\\site-packages\\Pythonwin',     'C:\\PROGRA~1\\QGISDU~1\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode',     'C:\\PROGRA~1\\QGISDU~1\\apps\\qgis\\python\\plugins\\fTools\\tools',     'C:/map stuff/MyProject']      Does anybody have any ideas?


title: ArcGIS 9.3 Export network dataset from geodatabase

text: I've got a set of feature classes with network dataset in a geodatabase *.gdb format and would like to export out to *.nds format, may I know how to do that on ArcGIS/ArcCatalog ?


title: Trace vector map from GeoTiff, Garmin

text: My goal is to create a vector map for my Garmin GPS by tracing over a topo raster. The vector data would be exported in ESRI SHP format and later converted for Garmin. So far, I can load the geo tiff layer, however am stumped on creating a vector layer. For instance, where do all the shapes come from to trace on the map? Does anyone have advice on this process? Any pointers are appreciated. Thanks, Nick


title: How to add title to mapbox <iframe> embed

text: I am embedding a map I made in Tilemill/Mapbox onto a Wix website. Items of the map such as legend, title, and description are not automatically included in the embed. I learned how to add the legend by simply adding ,legend, in the code below. However, I do not know how to add the title and description of the map as you see it here: https://a.tiles.mapbox.com/v3/yolivens.hgck7b4h/page.html?secure=1#4/40.85/-94.48 Here is my current code: ![enter image description here](http://i.stack.imgur.com/1jxu1.png) What would I need to do to add the title and description? I am above all concerned about the title.


title: Setting the output parameter in a toolbox to be a PDF

text: I need to allow a user to input a path and name where they would like to save their new map book PDF on their computer/drive. I have created a script where the user can input their mxd with data driven pages, input their title map, append the two PDF documents and save the appended PDF, but it's all hard coded and I'd like to create a script where they can select where to save their new appended PDF when they open it in the toolbox. That might be confusing, sorry. I need the script to be dynamic. I'm stumbling on setting up the parameters. I'm unsure what the Output PDF data type should be set as in the toolbox. None of the selections seem to fit. Aside from "Output" in the "Parameter Properties" should I change anything else? Thanks!


title: Table record counts change when using python tool

text: I have a python tool which investigates mxd files and clips/projects and exports all visible datasets. Recently I've found that the record counts (even on non-clipped data) is different. It's ok if I right-click save data. Not sure why this is. ![enter image description here](http://i.stack.imgur.com/ZFZ2y.gif) OutFile should be 716k records and not just 101k records as reported here. I am running a difference process in QGIS to see where the issue is in the file but I don't want to have to do this if the export works correctly.               ### Copies, clips and creates symbology of all visibile layers in mxd's in the current directory. Also creates a text file with metadata for use.     # Author: George Corea, Atherton Tablelands GIS     # georgec@atgis.com.au; info@atgis.com.au          # Licence:Creative Commons          import arcpy, string, datetime, shutil, os, glob     import arcpy.mapping as MAP     from arcpy import env          arcpy.env.workspace = os.getcwd()     arcpy.env.overwriteOutput = True          MXDList=string.split(arcpy.GetParameterAsText(0), ";")     ProjectPath=arcpy.GetParameterAsText(1)     CreateData=arcpy.GetParameterAsText(2)     clip_features=arcpy.GetParameterAsText(3) # polygon to clip data to AOI     AOI=arcpy.GetParameterAsText(4) #appended as suffix to dataset     arcpy.env.outputCoordinateSystem=arcpy.GetParameterAsText(5) #to maintain all datasets in the same projection     createdir=arcpy.GetParameterAsText(6)          rootPath=ProjectPath          #MXDList=glob.glob('*.mxd')     #ProjectPath=r'P:\2012\183_TownPlanning_Symbology\Working' # root output directory     #clip_features=r'L:\Vector_Data\Administrative\Boundaries\Local_Govt\TRC\trc_boundary_Polygon.shp' # polygon to clip data to AOI     #AOI='_trc' #appended as suffix to dataset     #arcpy.env.outputCoordinateSystem=r'L:\Vector_Data\Administrative\Boundaries\Local_Govt\TRC\trc_boundary_Polygon.prj' #to maintain all datasets in the same projection          xy_tolerance=1 #clip tolerance          # No edits should be required below this line.          VisibleLyrList = []     count=1               def layer_details(outPath, outFileN, lyrFile, type): #Generates the metadata              descLayer = arcpy.Describe(lyrFile.dataSource)         ReviewLog=outPath+'\\'+type+'_'+outFileN+'_log.txt'         f = open(ReviewLog, 'a')         try:             f.write(str(lyrFile.name)+': name{'+str(lyrFile.datasetName)+\         '}; query{'+str(lyrFile.definitionQuery)+\         '}; source{'+str(lyrFile.dataSource)+\         '}; description{'+str(lyrFile.description)+\         '}; symbology{'+ str(lyrFile.symbologyType)+\         '}; original projection{'+str(descLayer.spatialReference.name)+\         '}; extent(x,y){'+str(descLayer.extent.XMax)+','+str(descLayer.extent.XMin)+','+str(descLayer.extent.YMax)+','+str(descLayer.extent.YMin)+\         '}; format{'+str(descLayer.shapeType)+\         '}; size(bytes) ~{'+str(os.path.getsize(lyrFile.dataSource))+\         '} @{'+str(datetime.datetime.now())+'}'\         )             f.close()         except:             pass          print 'starting...'+str(MXDList)          for MXDFile in MXDList:              mxd=arcpy.mapping.MapDocument(MXDFile)         if createdir=="false":             outPath = ProjectPath         else:             outPath = ProjectPath+'\\'+mxd.filePath[mxd.filePath.rfind('\\')+1:mxd.filePath.rfind('.')]              arcpy.AddMessage ('Working on file #' + str(count) +' ...'+str(mxd.filePath))         try:             os.mkdir(outPath)         except:             pass         dfList = arcpy.mapping.ListDataFrames(mxd)         #print arcpy.GetMessages()         #print 'working on 2...'+str(dfList)+str(outPath)              for df in dfList:         #msd = outPath+'.msd'         #arcpy.mapping.ConvertToMSD(mxd, msd, df, "NORMAL", "NORMAL")         #print (str(count)+'...' +str(outPath)+'...'+'\n')             arcpy.AddMessage ('Working on dataframe ... ' +str(df.name))             lyrList=arcpy.mapping.ListLayers(mxd, "", df)             outPath = ProjectPath+'\\'+mxd.filePath[mxd.filePath.rfind('\\')+1:mxd.filePath.rfind('.')]+'\\'+str(df.name)             try:                 os.mkdir(outPath)             except:                 pass             print ('working on 3...in: '+str(outPath))             for lyrFile in lyrList:                 #print (str(lyrFile))                 if lyrFile.isFeatureLayer == True:                     if lyrFile.visible == True:                         if lyrFile.name not in VisibleLyrList:                             VisibleLyrList.append(lyrFile.name)                             arcpy.AddMessage (str(lyrFile)+' is visible PROCESSING...')                             outFileN=str(arcpy.ValidateTableName(lyrFile.longName[lyrFile.longName.rfind('\\')+1:]))                             try:                                 rows = arcpy.SearchCursor(lyrFile.dataSource)                                 row = rows.next()                                 if row:                                     if CreateData=="true":                                         arcpy.FeatureClassToFeatureClass_conversion(lyrFile.dataSource, outPath, outFileN)                                         #arcpy.Copy_management(lyrFile.dataSource, outPath+'//'+outFileN, "")                                         arcpy.Clip_analysis(outPath+'\\'+outFileN+'.shp', clip_features, outPath+'\\'+outFileN+AOI+'.shp', xy_tolerance)                                         #updateLayer = outPath+'\\'+outFileN+'.lyr'                                         sourceLayer = arcpy.mapping.Layer(outPath+'\\'+outFileN+'.shp')                                         sourceLayer_AOI = arcpy.mapping.Layer(outPath+'\\'+outFileN+AOI+'.shp')                                         #arcpy.mapping.UpdateLayer(df, updateLayer, sourceLayer, True)                                         arcpy.SaveToLayerFile_management(lyrFile,outPath+'\\'+outFileN+'_sym.lyr', "ABSOLUTE")                                         arcpy.SaveToLayerFile_management(sourceLayer_AOI,outPath+'\\'+outFileN+AOI+'.lyr', "ABSOLUTE")                                         arcpy.ApplySymbologyFromLayer_management (outPath+'\\'+outFileN+AOI+'.lyr', outPath+'\\'+outFileN+'_sym.lyr')                                         descLayer = arcpy.Describe(sourceLayer)                                         #QA Row Count                                         QA_rowCount_InFile = arcpy.GetCount_management(lyrFile)                                         outLyrFile=arcpy.mapping.Layer(outPath+'\\'+outFileN+AOI+'.shp')                                         QA_rowCount_OutFile = arcpy.GetCount_management(outLyrFile)                                         if int(str(QA_rowCount_InFile))==int(str(QA_rowCount_OutFile)):                                             arcpy.AddMessage("Count same - QA should not be Required..."+str(QA_rowCount_InFile)+" Out="+str(QA_rowCount_OutFile))                                         else:                                             arcpy.AddMessage("Count NOT same - QA IS Required...In="+str(QA_rowCount_InFile)+" Out="+str(QA_rowCount_OutFile))                                             layer_details(outPath, outFileN,lyrFile,"Row Count ISSUE --- Count NOT same - QA IS Required...In="+str(QA_rowCount_InFile)+" Out="+str(QA_rowCount_OutFile))                                             IssueLog=rootPath+'\\'+'Issue_log.txt'                                             f = open(IssueLog, 'a')                                             f.write(str(lyrFile.name)+': name{'+str(lyrFile.datasetName)+ '}; source{'+str(lyrFile.dataSource)+'}'+'\n')                                             f.close()                                         #QA Row Count process                                         layer_details(outPath, outFileN, lyrFile,"COMPLETED")                                         arcpy.Delete_management(outPath+'\\'+outFileN+'_sym.lyr')                                         arcpy.Delete_management(outPath+'\\'+outFileN+'.shp')                                     else:                                         layer_details(outPath, outFileN, lyrFile,"SIMULATED")                                     #arcpy.LayerToKML_conversion(outPath+'\\'+outFileN+'.lyr', outPath+'\\'+outFileN+'.kmz')                                 else:                                     arcpy.AddMessage ("!!!Datasource Issue!!!...continuing")                                     layer_details(outPath, outFileN,lyrFile,"ISSUE")                                     IssueLog=rootPath+'\\'+'Issue_log.txt'                                     f = open(IssueLog, 'a')                                     f.write(str(lyrFile.name)+': name{'+str(lyrFile.datasetName)+ '}; source{'+str(lyrFile.dataSource)+'}'+'\n')                                     f.close()                             except:                                 errorm=arcpy.GetMessages()                                 arcpy.AddMessage ('!!! ERROR !!!...'+str(errorm)+' ...continuing')                                 try:                                     layer_details(outPath, outFileN, lyrFile,"ERROR")                                 except:                                     pass                                 #break                                 ErrorLog=rootPath+'\\'+'Error_log.txt'                                 f = open(ErrorLog, 'a')                                 f.write(str(lyrFile.name)+': name{'+str(lyrFile.datasetName)+ '}; source{'+str(lyrFile.dataSource)+'}'+'\n')                                 f.close()                                       else:                                 pass                         else:                             arcpy.AddMessage (str(lyrFile)+' Is NOT Visible. Not Processing')                     else:                         arcpy.AddMessage (str(lyrFile)+' Is NOT Feature Layer. Not Processing')             count=count+1                 #print str(lyrFile)+' is not visible'              #Remove variable reference to file         del mxd, outPath, lyrFile      See code from rows = arcpy.SearchCursor(lyrFile.dataSource) onwards


i am using arcgis 10.1 i want to calculate the shortest distance between a raster and another raster. i have a bathymetric raster and an underlying fault plane raster. for each bathy raster cell, i need to know the shortest distance to the underlying fault plane. the bathymetry and the fault plane rasters are not planar nor horizontally level. i am not trained in python and hope there is an existing tool that i could use (though i could learn python if i need to).


title: How does the geoprocessing results window store the tasks?

text: Yesterday I was working in a single mxd all day, performing dozens of calculations and other analysis using multiple tools. A lot of it was trial and error, so when I found the correct combination I would label it as such in the `Results` window before documenting it in OneNote. ![Geoprocessing results](http://i.stack.imgur.com/gNELX.png) I opened and closed that mxd multiple times, and was able to rerun those tools I had labelled, as expected. When I opened that mxd this morning at work, I found that the tools I had renamed had been reset to their default names - as a result, I've had to trawl through each result until I found the ones I wanted. What I would like to know is, seeing as the processed tasks persist between sessions as it stored inside the mxd, so even if I open it on a different computer the results are still there, why does it "forget" that I renamed some of them?


title: EPSG 3857 or 4326 for GoogleMaps, OpenStreetMap and Leaflet

text: The discussion at What is the difference between WGS84 and EPSG4326? shows that _4326 is just the EPSG identifier of WGS84._. Wikipedia entries for Google Maps and OpenStreetMap shows that they both use WGS 84. http://wiki.openstreetmap.org/wiki/EPSG:3857 states that _EPSG:3857 is a Spherical Mercator projection coordinate system popularized by web services such as Google and later OpenStreetMap._ Leaflet's help states: _EPSG3857 The most common CRS for online maps, used by almost all free and commercial tile providers. Uses Spherical Mercator projection. Set in by default in Map's crs option.|_ _EPSG4326 A common CRS among GIS enthusiasts. Uses simple Equirectangular projection._ This is confusing - it seems that Google Maps and OpenStreetMap use EPSG3857 but they use WGS84 which 'is' EPSG4326. Something can't be right here, most likely my understanding. Could someone help me understand?


title: Convert ESRI Shapefile to Raster using R - biggest polygon area to raster cell value?

text: I want to create a raster from an ESRI shapefile. Multiple polygons may exist in one raster cell. If so, I want to assign the attribute variable which occupies the biggest area of the raster cell. If, for example, _Polygon1_ occupies 10%, _Polygon14_ 25% and _Polygon5_ 65% of the raster cell, the cell should contain the attribute variable of _Polygon5_. ![rasterize by polygon with maximum area](http://i.stack.imgur.com/0ILu4.png) Using the R-pacakge rasterize{raster}, values are transferred if the polygon covers the center of a raster cell. Of course, functions like min, max, mean, or character values like 'first', 'last', 'sum', 'count' are available. Default is 'last'. You can download the sample shapefile here. Up to now my code is this:               setwd("your_wd")          require(sp)     require(raster)     require(rgdal)          sp_df <- readOGR(getwd(), "stok_subs") # reads the shapefile into a large polygon                                            # data frame          extent(sp_df) # shows the extent of sp_df          summary(sp_df) # simple summary          # Object of class SpatialPolygonsDataFrame     # Coordinates:     #        min      max     # x 32539350 32544695     # y  5732607  5735439     # Is projected: TRUE      # proj4string :     # [+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=32500000 +y_0=0 +ellps=GRS80+units=m+no_defs]     # Data attributes:     # X_CENTR           Y_CENTR            ATTRIBUTE      # Min.   :3539714   Min.   :5734725   9.4-.2.3 :21       # 1st Qu.:3540436   1st Qu.:5735323   23.4-.2.3:19     # Median :3541226   Median :5735830   9.4.2.3  :19       # Mean   :3541263   Mean   :5735824   23.4.3.5 :18       # 3rd Qu.:3542031   3rd Qu.:5736358   23.4.2.3 :16       # Max.   :3543461   Max.   :5737024   19.4.3.5 :12       #                                     (Other)  :89            r <- raster(extent(sp_df)) # creates a raster with the extent of sp_df     projection(r) <- proj4string(sp_df) # uses the projection of the shapefile                                         # for the raster         res(r)=10.0 # sets the resolution of the raster to 10 m          r10 <- rasterize(sp_df, field="ATTRIBUTE", r) # converts sp_df to a raster                                                   # assigning the last attribute                                                   # variable in a cell (default)          writeRaster(r10, filename="stok_subs.tif", format="GTiff") # writes the raster                                                                # to a file          plot(r10) # plots raster      ![r10](http://i.stack.imgur.com/FIWv1.png) None of the provided functions seem to cover my needs. However, I'd rather operate with R, as operations with very large raster data (e.g. whole Lower Saxony @ a resolution of 10 m) may be easier to handle than with GIS solutions. Anyone out there who has a (or any) solution for my problem?


title: Installed PostGIS extension not listed for database

text: I have a PostgreSQL database with the PostGIS extension enabled. I know that this extension is enabled, since I am using PostGIS functions, and I can actually see the PostGIS functions. However, PostGIS is not listed in the installed extensions section, and thus I can't alter it nor uninstall it. If I try               CREATE EXTENSION postgis;      it tells me that PostGIS is already installed, and that I should uninstall it first. However               DROP EXTENSION postgis CASCADE;      tells me that no PostGIS extension is installed. What can I do? Any ideas? Thanks in advance.


title: Using qgis.core on redhat/linux

text: I just need the path to the api so I can import qgis.core, etc. I was able to import on my mac, run a standalone code, and test; however, my team requires it to be put on jenkins. qgis was installed via redhat on the linux-skybox. I found where qgis using "where qgis" and got "/usr/bin/qgis" to find qgis.core inside I tried the following: I believe they've installed version 1.8 of qgis cd the following /user/bin/qgis /user/bin/qgis/share /user/bin/qgis/share/qgis /user/bin/qgis/Contents /user/bin/qgis.app /user/bin/qgis.app/Contents And I got "no such directory exists" edit 1: I can import qgis.core on the main computer of linux-skybox, but I can't on another edit 2: I needed to use virtualenv jenkinsEnv --system-site-package instead of virtualenv jenkinsEnv --system-site-package


title: Converting LatLng coordinates to point shapefile

text: I have more than 50,000 LatLng coordinates defined under SWREFF99 1800 projection. I want to know, is it possible to convert these coordinates into a point shapefile ? if yes, then how can I do this ?


title: Why would programmatically setting the extent of a map on Page_PreRender be really slow?

text: We are using the ArcGIS 10 web adf with C# in a web application. We just recently upgraded everything from 9.3 to 10 and moved all the hosted stuff over to a new (more powerful) server. Still on the same network. During the Page_PreRender event, we get the full extent envelope of our map, which is fast, and then set our map's extent to it. Code in question...               Envelope fullExtent = AgsHelper.GetExtent(Map1, ExtentType.Map);     Map1.Extent = fullExtent;      The second line there takes a really long time (like 30-40 seconds). I can't claim to be an ArcGIS expert. I'm wondering if anyone has ideas of what I could look at or what it might be. I decompiled the code that does the setting of the extent and took a look at the code, but I can't really tell much without being able to debug into it and I can't do that since I don't have a reflector license and my trial is up. Pretty sure this is some sort of SDE or ArcGIS Server issue, I just don't know what it would be. Thanks!


title: Complete lookup table for German admin codes in Geonames data? (admin4_code)

text: I am working with both the postal and Geonames data for Germany. I noticed that while the Postal data may be used as a lookup for ADMIN_CODES 1-3 (ADM1,ADM2,ADM3), it does not mention ADM4 at all. Does anyone know where to find the lookup table for 'admin4_code' values? According to http://forum.geonames.org/gforum/posts/list/538.page, this admin code references a 'Gemeinde', an admin field I would like to include in my work, if possible. * * * The German dataset I'm working with: http://download.geonames.org/export/dump/DE.zip The German dataset I'm using as a lookup table for ADMIN CODES: http://download.geonames.org/export/zip/DE.zip Both zipped folders contain a 'readme.txt' file that outlines the fieldnames in the data. I am also tracking this question on the Geonames forum: http://forum.geonames.org/gforum/posts/list/2623.page


title: My woes rendering shape files on web-based mapping application

text: First, I'd like to let it be known that this is my first journey into developing mapping applications with such a large set of data. So, I am in way locked into the techniques or technologies listed below... The Overview What I am trying to accomplish is the ability to load my company's market information for the United States onto some form of web-based mapping software. Currently I am trying to accomplish this with Google Maps JS api and GeoJSON data, but am not against other alternatives. There are roughly 163 files (ranging from 6MB to 200KB) that I have exported from ArcGIS into GeoJSON files, and then loaded into a database (currently as GeoJSON strings). We have developed a UI that loads the data based on the current map bounds and Max/Min calculations in the corresponding records. The problem I'm running into is the render time on the map itself, which is annoying when switching between different regions, states, or zoom levels. The API calls to load the data are acceptable in regards to the size of the data being retrieved. My boss said it is great for a proof of concept, but would like to see it much, much faster. The Questions What suggestions could you offer to increase the render time? Are there better options (3rd party libs, techniques, etc) for what I'm trying to accomplish? It was suggested by a co-worker to export the map shapes and just use the images for overlaying the information based on the coords. Any input on this? **Edit #1** The shapes consist of polygons and multi-polygons. They essentially trace around all of our markets with distinctive opaque fill and border. The interactive ability would be nice, to be able to click on a given area and retrieve property, tract, and other data points, but is not essential to it's success. I could overlay the data if necessary.


title: Geoserver Datastore Property Update via REST API

text: I am trying to change the properties of an existing Geoserver data store via the REST API. I have created the data store from an existing folder of shapefiles. However, I need to change the "charset" and "Create spatial index if missing/outdated" properties but I cannot for the life of me figure out how to do it. I have tried sending the following XML as a PUT request to the data store but it doesn't work:               <DataStore>       <parameters>         <Parameter>           <name>create spatial index</name>           <value>true</value>         </Parameter>         <Parameter>           <name>charset</name>           <value>ISO-8859-1</value>         </Parameter>       </parameters>     </DataStore>      REST:               http://localhost/geoserver/rest/workspaces/uk/datastores/roads.xml     PUT /geoserver/rest/workspaces/uk/datastores/roads.xml      Anyone have any ideas? Thanks


title: Preparing Huff model with data to use in QGIS?

text: I'm trying to find out about the notion of potential of acessibility using the methods of Huff. I´m using QGIS, with point data (with the variables that I would like to use) and polygons for Europe NUTS and LAU-2. What I am looking for is the ability **to calculate areas of influence of the different points using the method of Huff.** With QGIS, I cannot find tutorials or videos about how to make it work with the python console. **So can there be any examples about to prepare a Huff model with data that we might use on QGIS?** _(I´ve already seen the posts relative to the Huff method on the forum but I do not know how to make it work properly with the python language)_ _PS: I have a correct background with python programming, but I have problems with the way of making it work with GIS methods for this kind of theme._


title: How to share styles between layers in QGIS

text: Is it possible to share one style between different layers in QGIS? Now if I apply the same style from *.qml or use the same style from Style Manager for different layers then in *.qgs project file each layer contains full copy of the same symbology description. In this case not possible to discriminate layers which use the same style from others. Maybe possible to set name for used style?


title: Rendering performance

text: I am using ArcGis Java SDK for creating a map in my JavaFX application. I have a requirement to plot up to 500,000 graphics on graphics layer. When I do 100,000 graphics its performance is really bad. It takes 10-15 seconds to render the Base map and Graphics layer. I cannot use any online map because this application is running on a machine that has not internet connection.                   jMap = new JMap();           ArcGISLocalTiledLayer tiledLayer = new ArcGISLocalTiledLayer(“Texas_UL3D_Test.tpk");           GraphicsLayer graphicsLayer = new GraphicsLayer();                jMap.getLayers().add(tiledLayer);           jMap.getLayers().add(graphicsLayer);           jMap.addMouseListener(new MouseClickListener());                final HitTestOverlay overlay = new HitTestOverlay(graphicsLayer);                overlay.addHitTestListener(new                     jMap.addMapOverlay(overlay);                        jMap.addMapOverlay(new MouseMoveOverlay());                Map<String,Object> attributes = new HashMap<String, Object>();                jMap.addMapEventListener(new MapEventListenerAdapter() {                   @Override               public void mapReady(MapEvent event) {                  SpatialReference mapSR = event.getMap().getSpatialReference();                  for(Data data : getAllData()){                     attributes.clear();                     Point point = GeometryEngine.project(data.getLongitude (),data.getLatitude (), mapSR);                     attributes.put("number", data.getNum());                     attributes.put("name", data.getName());                     Graphic pointGraphic = new Graphic(point, getSimpleMarkerSymbol(data.getStatus()),attributes);                   graphicsLayer.addGraphic(pointGraphic);                                                }               }           });                                                 chartSwingNode.setContent(jMap);       Please help me to improve its performance. Tell me what I am doing wrong. ![A map with 200k+ graphics](http://i.stack.imgur.com/05HUw.png) A map with 200k+ graphics


title: Measure Snapping with ArcGIS Javascript API

text: I would like to know if the ArcGIS Javascript API has built in snapping functionality for the measure tool. I am trying to create a measure tool which only measures distances along the edge of roads. Is there anything built-in or would I need to create this functionality by linking into ArcObjects?


title: comparing two TINs

text: i have the following two files: 1) LiDAR las file1 2) Let's give it a name LiDAR las file2, contains the same file1 data, where number of points randomly removed. So if file1 contains X number of points, file2 would contain Y number of points, where Y is less than X. Now, i am using ArcGIS "LAS to Multipoint" to create from both files above their shape counterparts. Next i am creating two empty TIN(triangulated irregular network) files using ArcGIS. Then i edit both TINs, let's say tin1 and tin2, with the appropriate input feature clasess, i.e. file1 to tin1 and file2 to tin2. Is there any tool/package/library or any software that would let me calculate the RMSE(root mean square error) or any difference information that would give me information on the accuracy of the second file, i.e. file2, with respect to file1? I tried to use the TIN difference but it gives system error and cannot obtain the result. I would really appreciate if anyone could show how to obtain accuracy information and interpretation of the results obtained... thanks a lot in advance!!!


title: Join does not return matching records after 1001 records have been successfully joined

text: I'm working on a process that joins a feature layer to another feature class. The code is as follows:               arcpy.MakeFeatureLayer_management ("New_FC", "newFC_layer")     arcpy.AddJoin_management("newFC_layer", "ID",Join_FC, "ID", "KEEP_ALL")     env.qualifiedFieldNames = False     arcpy.CopyFeatures_management("newFC_layer", newFC_Output)      The problem I am having is there are over 1000 records and the result of the join only contains 1001. I have tried this manually in Arc by first creating the layer from the original feature class, joining the layer to the desired feature class and using the copy features tools just as the scripting would and I get a desirable outcome. I do not have domains in the ID field. I've checked the values in the field to be certain there is no variation between the 2 ID fields. I'm not sure what else could be causing this. Your suggestions are much appreciated.


title: Multiply Rasters with Scalar Numbers as Input in QGIS Modeler

text: I built a little model in the QGIS modelbuilder which uses the raster calculator with the command > gt(a,threshold) I need to run the model with varying thresholds from 2 to 30 so it would be great to create a _Number_ input in the model for this threshold so I'm asked for its value when starting the model, but the raster calculator only accepts rasters as inputs. Maybe someone has an idea how to get around this problem? For example by generating a raster containing only the threshold value in every cell in a previous step if there is an easy way to do this?


title: How to MOVE and RESIZE Features within a Layer?

text: I am still new to ArcEngine 10, I currently am able to create a feature layer that is associated with a feature class. Its a polygon layer. So I can create polygon, give it attritbute table and put it on a layer using UNIQUEVALUERENDERER. My question is, is there a tool bar control, or any way I can MOVE or RESIZE my polygon easily. I current use `"esriControls.ControlsSelectFeaturesTool"` to select the tools. Is there an easy way to do this? Thanks and Regards


title: Source for UTM Zone File

text: I'm looking for a complete, freely available file showing UTM Grid zones. Preferably _not_ a simplified version, but one that shows the full grid zones and all 'exceptions' (e.g. in Northern Europe) - as per http://whatutmzoneamiin.blogspot.com/ Does anyone know where I can download such a file?


title: How to tell if an ArcGIS Service is secure?

text: I am creating a .NET REST SOE for ArcGIS Server 10, and I will need to determine if a map service (or image) is secure or not. Is there a property on a MapService that would tell me this in either ArcObjects or REST?


title: What server and API to use for displaying complex polygons with dynamic symbology?

text: **Background** I need to interface a map into an existing website which allows the user to query a database and which currently results in a table listing population statistics per county. For example, display the population under 15, for 2003 - the result is a table listing each county's population values. I need to display a map showing county polygons which are symbolised with those values. If the user modifies the query, the map needs to update (eg, if the year is changed to 2004, or the statistic is changed to population over 65). Performance will be a factor, so the map should update rapidly when the query is run/modified. The map should support mouse-over events on the counties, so that statistics can be displayed in a popup as the user's mouse moves. **My questions** What are the options for the GIS server (to serve the polygons) and the mapping API (to embed the map on the website)? **Options** I like the look of Leaflet, but it seems that the county polygons would need to be supplied as coordinates/rings. Would this be too slow when zoomed to national scale? Can Leaflet handle on-the-fly simplification (as with feature layers in ArcGIS Server)? Using OpenLayers with GeoServer could work too, but the performance is quite slow for hover events, and it appears that the symbology is hard-coded in the WFS service(?) I'd need the ability to update the symbology when the query was modified. ArcGIS Server 10.1 has advanced on-the-fly symbology options, but hasn't been released yet. The client will be wary of deploying beta software on a production site. Thanks for any advice.


title: Looking for large GML files containing vector polygon data

text: I have few vector data files with 500 to 700 MB GML files. But I am looking for files larger than these preferably 1 GB, 2 GB. It would be helpful if somebody can provide some pointers regarding the place to find such files. Note: I need to test polygon overlay over large data sets, so the two layers that I need should be intersecting.


title: How to get the second nearest neighbor between two point patterns in R?

text: Is there a way to get the distances for the second nearest neighbor between two point patterns in R? The spatstat package has a function called nncross but it only applies to the nearest neighbors between two patterns and I need the distances to the second nearest neighbors.


title: Does anyone know of a 'best practice' Seabed Survey Data Model (SSDM)?

text: Many of the offshore oil/gas companies are adopting an initiative from Shell for delivery of spatial data. This is called an SSDM, which is really an ESRI file geodatabase with a standardised structure of datasets, and featureclasses. An empty template (geodatabase) in the SSDM structure is freely available from the OGP website, but it would be really good if anyone knows of a 'best practice' example geodatabase which is fully verified so we can use it as a 'how to' reference guide. Does anyone know if such a thing exists?


title: PostGIS spacial data type to represent geo-located polygons

text: I have a PostgreSQL database with circle shaped objects called 'spots'. Each of those objects has latitude, longitude and a variable radius, information that I need to define this circles. Given an specific coordinate, I select all the objects that contain that given point. The thing is that I want to switch to PostGIS this stuff because of performance. I've started documenting about PostGIS but as a complete newbie I still have some doubts about data types. What is the best way to store this kind of information in PostGIS? I don't need precision, because the radius can only have a limited lenght of 100km. I've thought about creating a polygon using `ST_Buffer` and save this into a `geometry` field in the 'spots' table. I also do not need an exact circle, a decagon is enough. The test code I'm using to know if this is meeting my needs is the following               create table spots (name varchar, geom geometry);          -- Insert a 2500m radius polygon with center in given coordenates     insert into spots      values (         'Test',          ST_Buffer(             ST_SetSRID(ST_MakePoint(39.63467934415986, 2.6308822631835938), 4326),             2500,             4         )     );      I have several doubts about this code. Once the geometric figure is created with `ST_Buffer`, is it still preserving it's position? If I provide the SRID code with given lat-lon coordinates, do I need to use `ST_SetSRID` or it's completely unnecessary? The thing is that I've discovered that using `ST_Contains` is always returning **true** , with no matter about the lat-lon values I'm providing to the function.


title: Which Geodatabase supports SQL and is easily portable?

text: I want to choose the geodatabase type, which could be:   * permits direct sql queries.   * easily portable, because the app would installed in one computer. I know that the personal geodatabase is a good choice, but I would like to know if there's some limitations in personal geodatabase apart the size (the size won't be a problem). The geodatabase would be used for some geoprocessing operations and some sql queries.


title: How to get nearby tiles from a location?

text: I want to get the tile in which a particular location lie and also its nearby tiles. Is there any google map api which could return this information.


title: What is a shaved polygon?

text: Just trying to find a definition. I was asked this in an interview and I assumed it was a polygon created by incorrectly modifying a another polygon in a coverage. According to my best googling resulting in this code: http://www.falconview.org/svn/FalconView/trunk/public/fvw_core/Utility/GeometryUtils.cpp and my limited understanding of it. A shaved polygon occurs when you eliminate a polygon and your surface no longer forms a convex hull. I would assume this is important in a situation where keeping area constant (such as city limits).


title: How can I snap 2 lines in QGIS so the st_distance between them is 0?

text: Sometimes when I snap the ends of 2 LineStrings together and I check the st_distance between those endpoints I get a value in the order of 1.1E-10 while at the same time when I check for equality of those endpoints I get the value 'true'. How can I snap 2 lines so the st_distance between them is 0? I am using QGIS 1.9.0-Master and a PostgreSQL database 9.1.7 with PostGIS 1.5.3-1.


title: Join/Relate in ArcGIS 9.3.1 ArcEditor through a python script - is it possible?

text: I'm working on a ArcGIS 9.3.1 ArcEditor and I want to write a python script that will do a Join or Relate. Problem is, that join table\join field are restricted to ARCINFO (higher license), and relate isn't doable by code (or so it seems from the ARCGIS help). As upgrading the license isn't an option, Does anybody have a way I can do a join between two shapefiles in ArcGIS or Opensource tools?


title: How to create a shapefile in ArcGIS 10 using arcpy?

text: How can I create a shapefile using Python in ArcGIS 10? I have lat & long. From this I need Python code which will create a shapefile and display it in ArcGIS 10. Thanks.


title: ArcGIS online options?

text: There are lots of comments on different forums about different types on online GIS data/map sharing tools. What i need to know is what ArcGIS online CANT do compared to the functionality of the open source options? I have big geodatabases that contain dynamic data that is constantly being updated and shared across multiple offices. The geodatabases contain raster data with information that needs querying - Lidar/bathymetry data. We also have vector features that could be as many as 2000 features in any one feature layer. I have heard that it is a struggle to get raster layers into AGOL and they will just be an image with no data behind them I have also heard that there is a limit on the number of features AGOL can display at any one time? It seems a total minefield as to which one to use but i would prefer the free option if the functionality is the same or better in opensource.


title: Migrating to Pgrouting pgr_dijkstra?

text: I have a system that uses shortest_path function pgrouting 1.x. At the moment, I want to migrate to pgrouting 2. And I assume the equivalent function to shortest_path is pgr_dijkstra. (Correct?) When running this, (Hibernate Query Language)               SELECT SUM(t2.length) as totalLength FROM (         SELECT * FROM pgr_dijkstra(             'SELECT gid AS id, source::int4, target::int4,                      cost::double precision AS cost,                      reverse_cost::double precision as reverse_cost                  FROM openstreetmap',5374,5364, false, true         )) AS t1, openstreetmap AS t2 WHERE t1.edge_id = t2.gid      I get error ERROR: column t1.edge_id does not exist What is the edge_id in the output of pgr_dijkstra ?


title: What problems do special characters in file names cause?

text: When defining file names and paths, I have always used a simple rule - No special Characters including foreign letters Ä, Ü, Ö, ß (only _ and -) and no spaces. I used these rules to avoid problems with geoprocessing in older versions of GIS software products which used to complain about such paths and file names. I am now wondering what is the situation now? Are the rules still relevant to todays software products (not just GIS Software) and are there actually any guidelines or ISO standards regarding this?


title: What is the correct order for compressing, rebuilding indexes and updating statistics in ArcGIS?

text: I want to script out the compressing, rebuilding of indexes and updating statistics on our SDE databases, but I want to make sure I have the correct order for these to run. What would be the correct order? Should I run certain ones more than others? For the order, I would assume:               1) Compress     2) Rebuild indexes     3) Update statistics      But really I'm not positive. Thanks!


Is there a way to close a geoprocessing operation when the dialog box has been hidden/pushed to the back of the windows? I ran a python process as VB Script and when the error message appeared I misclicked. The error popup was hidden behind the attributes table. I finally terminated ArcMap. Is there an easier solution?


title: How to choose a GRASS kernel density function?

text: I'm starting to use the Grass module on QGIS, and I want to do a heatmap, so, I'm creating a new surface using the kernel function. But, GRASS come with 7 types of kernel functions, how do I know which one to choose? Regards


title: How can I create a multipoint geometry in QGIS?

text: I've got a Spatialite layer in QGIS with geometry type MULTIPOINT. Now I'd like to add some data in this table in the 'usual' way, so by doing multiple clicks in the map (e.g. left-click, left-click, right-click -> enter attribute values in the popup dialog for creating a record with a 3-points geometry. However, QGIS only allows me to add a single point by popping up the attribute dialog right after the first left-click. How to do this in the GUI? (I don't want to enter coordinates by hand in a query :-) )


title: How to export from PostGIS to KML

text: What is the best way to export from PostGIS to KML using open source tools? Can I do it directly, or do I need to first export to shape files?


title: How to redraw WMS layer on OpenLayers without flickering from geoserver?

text: I have a WMSlayer that is updated every 10sec with new data from SQL. Whenever I redraw the layer using redraw(true) , the map screen flickers. Is there a way to do it without flickering ?


title: GRASS GIS SDTS Manipulation

text: Preface: I'm a total GIS noob. Flame away. I'm learning GIS in my free-time and with free resources. I realize that is likely going to force me on the learn-fast-or-giveup track, but here it goes... My first project: Create a laser cut 9-layer topographic map of the 7.5 minute area where I work. My Plan: 1) Download SDTS DLG Hypsography data from USGS. 2) Unzip using Peazip. 3) Import into GRASS GIS. 4) Remove/edit unnecessary topo lines (I was able to do it with dlgv32 available at USGS website, but I can't export and don't want to pay $500 for a full version of global mapper). 5) Export to DXF or DWG 6) Open in Corel and create layers 7) Laser cut layers, glue them together. I've gotten through steps 1, 2, and 3 with minimal frustration, although peazip throws an error on extract that says "data after archive", which I hope has no effect on the usability of my data. I'm at a loss on step 4. That data set I'm working with is DLG STDS data in large format (7.5 minute quadrangles). Specifically, the data is from: http://dds.cr.usgs.gov/pub/data/DLG/LARGE_SCALE/I/indian_hills_CO/hypsography/1/ I can get the DDF files to import into GRASS and I can view the topomap in the display; however I'm having trouble manipulating what I'd like to see. Basically I want to display only the major topolines (for example, 6000', 6400', etc.). When I open this data set in Global Mapper, this is very easy to do. I imagine I can accomplish this with some kind of SQL select function in GRASS; however, I do not understand how the data is organized enough to do this. When I import the SDTS file, I get 10 database links. Most of the databases look 'empty', so I imagine there is just one file I need to be working with. Do any of you have any experience in working with Hypsographic data? Do you know how I can simply delete the entries that do not correspond to the lines I want?


title: Is there a 3rd party ArcGIS viewer (builder) for Javascript?

text: Esri provides an ArcGIS Viewer for Flex as well as an ArcGIS Viewer for Silverlight. Strangely though, there are no plans for a similar viewer for Javascript. **Are there any 3rd parties that have created a "viewer" for javascript?** Though they are referred to as "viewers", they are actually viewer-builders: they allow administrators to create web applications without writing any code. The javascript viewer blogged about here does not allow user to choose tools to include in the app, as can be done with the flex and silverlight versions. * * * **Update** I have developed some capabilities that have been packaged as a combination of an SOE and an Add-in for the silverlight viewer. This works great - when the admin builds a viewer application, he can pick my add-in tool from the list and configure it to use the URL of a server where my SOE is installed. However, from what I can tell there are no add-ins for the javascript viewer. Or am I missing something?


title: Can I use predefined list of values from Postgis when I draw lines in QGIS?

text: Lets say that we have 2 tables in Postgis:               Table 1 (PIPE):     -ID     -pipe_type_id     -geometry          Table 2 (PIPE_TYPE):     -ID     -type      foreign key `PIPE.pipe_type_id - PIPE_TYPE.ID` Table PIPE is empty. Table PIPE_TYPE is filled with values. I need to draw a line in QGIS directly to PIPE.geometry column. After the line is drawn I need a popup window to appear with an option to choose PIPE_TYPE.type value for that line. After submitting, correct value of PIPE_TYPE.ID should be stored in PIPE.pipe_type_id column. I did this with Bentley Map using their Map Administrator and Oracle DB, but I want to try this with QGIS + Postgis. Is it possible?


title: How to solve Microsoft .NET Framework error message.?

text: The below error message displaying when user leaves the arcmap desktop for more than one hour without any performance on the map, if this error comes then user should close the Arcmap and re-open it again, any suggestions how we can over come this problem. **Arcmap Version 9.3.1** ![enter image description here](http://i.stack.imgur.com/FBInI.png) I got this message in the details box. _**_** _*_ ** _**JIT Debugging **_** _*_ ****** To enable just-in-time (JIT) debugging, the .config file for this application or computer (machine.config) must have the jitDebugging value set in the system.windows.forms section. The application must also be compiled with debugging enabled. For example: When JIT debugging is enabled, any unhandled exception will be sent to the JIT debugger registered on the computer rather than be handled by this dialog box.


title: What is the difference between ArcSDE and ArcGIS Server?

text: I have been working with ArcGIS Desktop where I use ArcMap for doing GIS work. Can anyone please explain to me in simple words What is ArcSDE and ArcGIS server? These can be independent dedicated server machines? Any simple and full explanation on this is most welcome.


title: Creating getCapabilities document in more languages

text: I'm providing INSPIRE metadata and services and I want them in more languages than one. I expect that most users are locals, but some users are definitely from other European countries. Therefore I want to support national language and english. In metadata documents, I have used gmd:PT_FreeText_PropertyType, but getCapabilities file follows other XSD schemas (wms or wfs, according to the service). Do you know about some possibilities of adding additional language to getCapabilities document?