title: From a trajectory of GPS points to a multiline segment

text: Given a trajectory as a sequence of GPS points I would like to create a smooth linestring representing a path where a car drove. The data is stored in postgresql: the input trajectory in a table, where each row represents an individual GPS point, the output in a table, where each row represents the whole smooth tragectory. All points are stored in WGS84. My approach for this would be the following:   * implement everything in C++   * retreave the whole trajectory from the DB, convert it to UTM32N (as this is where the trajectories are comming from), flattern the points with the help of GDAL library   * compute a some sort of approximate spline   * linearize the approximate spline   * transform it back to WGS84   * send the data back to Postgres Any comments on the whole strategy and the implementation language? What would be the best tools/libraries to use for approximate splines and linearization? Java/C++/Postgis solutions are preferable, as currently we are using only these subsystems. I suppose this is a very common task people need to do in gis environment, however I haven't found close matches to similar questions on gis stackexchange. Any comments would be very helpful.. [Update:] An alternative would be to:   * collect all points into linestring   * ST_LineToCurve    * apply Bezier smoothening to the curve   * ST_CurveToLine Comments are very helpful..


title: Where are Print Composer Template files stored?

text: When I saved my Print Composer template, I cleverly navigated to a logical (to me) folder among my documents to save it. Now I find that it does not show up on the list of Composers under “File”. Can someone please tell me where it is supposed to be, so I can move it to there and open it?


title: Make Graph (Data Management), Changing item names in legend

text: I'm having python produce a multi-series graph. How do you change the series name in the legend, the default being Vertical Bar, Vertical Bar1, Vertical Bar2.... Thanks for the help               arcpy.MakeGraph_management(graphTemplatePath,'''SERIES=bar:vertical DATA={0} X=ParkUnit Y=Coverage LABEL=ParkUnit SORT=ASC;                                                  SERIES=bar:vertical DATA={1} X=ParkUnit Y=Coverage;                                                  SERIES=bar:vertical DATA={2} X=ParkUnit Y=Coverage;                                                  SERIES=bar:vertical DATA={3} X=ParkUnit Y=Coverage;                                                  SERIES=bar:vertical DATA={4} X=ParkUnit Y=Coverage;                                                  GRAPH=general TITLE={5}, Total Coverage by Unit;                                                  LEGEND=general TITLE=Year Total Coverage;                                                  '''.format(scratchReport + str(year), scratchReport + str(year-1), scratchReport + str(year-2),                                                             scratchReport + str(year-3), scratchReport + str(year-4), exoticPlant),"graphTemp")


title: Reading shapefile tolerance in python

text: Is there a way to get the cluster or xytolerance from a shapefile into python using free/open source software? I'm currently using the `pyshp` package and can't see methods for this anywhere. Also it doesn't seem to appear in the WKT of the `.prj` file.


title: Using search cursor to determine renaming of feature class

text: What I am trying to do is create a list of shapefiles based on field name values. From there I would like search on 1 field (each shapefile will only contain 1 row) and if the value in that field is greater than X have the code rename the shapefile. What I have so far is:               polyFCs = arcpy.ListFeatureClasses("*_Value.shp", "POLYGON")     for polyf in polyFCs:         field = "Area"         with arcpy.da.SearchCursor(polyf, (field)) as cursor:             for row in cursor:                 row.getValue(field)                 if field > X:                     arcpy.Rename_management(polyf, polyf.replace('.shp','') + "_NewValue")      When I run my code (and I'm pry missing or misunderstanding something) I get a runtime error stating 'A column was specified that does not exist'. I am in 10.1. Thanks for your help.


title: How to retain Field Name / Alias on Export of Table?

text: Is there a way I can modify the Field Name or to make him keep the Alias when a table is exported? Or is there a way to keep the Alias instead of the Field Name, when exporting for example to .dbf ? I want to be able to do this in ArcMap.


title: How can I change SQL query through a layer of an action (python)?

text: I have three layers from three diferent database tables. When I click in other layer I want to change (with an action) the query of that three layers. How can I do that? Thank you


title: Calculate Field tool to calculate on null fields

text: I'm using ArcGIS 10.1. I want to use the calculate field tool in ModelBuilder on a table where if field1 is null, it should take the value of field2. If field1 has a value then it should stay as is. I'm sure this is a pretty basic question but I don't have any experience writing expressions. Any input would be great! PolyGeo, your solution works like a charm. I also want to know what the expression would be for the exact same scenario above except instead of null, field1 is blank.


title: Mapping application basic for ruby on rails

text: I am new to build map application using ruby on rails, with rails 4.1.1. I want to create map using map tiles. i have to use PostgreSQL database with postgis. Can you please explain how this flow will work ? Form where i have to start for map creation and display. Thanks, Ghanshyam


title: How can I clip a shapefile in java using geotools?

text: Does anyone have any advice or sample code that shows how to clip a shapefile using geotools? I am looking to build a simple GUI to allow a user to input a bounding box and polyline shapefile, then clip the contents of the shapefile to the bounding box the user entered. The functionality I am trying to achieve is the same as ArcGIS's "Clip (Analysis)" Geoprocessor (http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//000800000004000000). UPDATED: I came up with the following code that seems to work for what I will be trying to do. I have no idea if it is the best way to go about accomplishing this functionality. I am brand new to geotools, but here it is for others.(use at your own risk ;D)               public static void main(String[] args) {         try {                  File file = JFileDataStoreChooser.showOpenFile("shp", null);             if (file == null) {                 return;             }                  FileDataStore store = FileDataStoreFinder.getDataStore(file);             SimpleFeatureSource featureSource = store.getFeatureSource();                  // Create a map content and add our shapefile to it             MapContent map = new MapContent();             map.setTitle("Quickstart");                  Style style = SLD.createSimpleStyle(featureSource.getSchema());             Style style2 = SLD.createSimpleStyle(featureSource.getSchema(),                     Color.RED);                  Filter filter = CQL.toFilter("BBOX(the_geom, -60, -50, 30, 40)");                  SimpleFeatureCollection featureCollection = featureSource                     .getFeatures(filter);                  SimpleFeatureType featureType = featureSource.getSchema();                  int geometryAttributeIndex = -1;             GeometryDescriptor geometryDescriptor = featureType                     .getGeometryDescriptor();             List<AttributeDescriptor> attributeDescriptors = featureType                     .getAttributeDescriptors();             for (int i = 0; i < attributeDescriptors.size(); i++) {                 AttributeDescriptor descriptor = attributeDescriptors.get(i);                 if (descriptor.equals(geometryDescriptor)) {                     geometryAttributeIndex = i;                     break;                 }             }             if (geometryAttributeIndex == -1)                 throw new Exception("Unable to find geometry attribute index.");                  ReferencedEnvelope clipEnvelope = new ReferencedEnvelope(-60, -50,                     30, 40, DefaultGeographicCRS.WGS84);             GeometryClipper clipper = new GeometryClipper(clipEnvelope);                  SimpleFeatureIterator iterator = featureCollection.features();                  List<SimpleFeature> clippedFeatures = new ArrayList<SimpleFeature>();                  while (iterator.hasNext()) {                 SimpleFeature feature = iterator.next();                 Object geometry = feature.getDefaultGeometry();                 Geometry clippedGeometry = clipper.clip((Geometry) geometry,                         false);                 if (clippedGeometry != null) {                     List<Object> attributes = feature.getAttributes();                     attributes.set(geometryAttributeIndex, clippedGeometry);                     SimpleFeature clippedFeature = SimpleFeatureBuilder.build(                             featureType, attributes, feature.getID());                     clippedFeatures.add(clippedFeature);                 }             }                  SimpleFeatureCollection clippedCollection = new ListFeatureCollection(                     featureType, clippedFeatures);             SimpleFeatureSource clippedSource = DataUtilities                     .source(clippedCollection);                  Layer layer = new FeatureLayer(featureSource, style);             Layer layer2 = new FeatureLayer(clippedSource, style2);             map.addLayer(layer);             map.addLayer(layer2);                  // Now display the map             JMapFrame.showMap(map);         } catch (Exception ex) {             JOptionPane.showMessageDialog(null, ex.toString());         }          }


title: ArcCatalog 10.0 complains that spatial references don't match between a feature dataset and a feature class that have the same coordinate system?

text: I have a feature dataset that is using the GCS_WGS_1984 as the geographic coordinate system. It includes a few feature classes. This feature dataset is itself in a file geodatabase that contains a few more feature classes in the root of the gdb. One of these is a feature classes called "points" that is also using the GCS_WGS_1984 as its coordinate system. I thought it would be trivial to use ArcCatalog to drag the feature class from the root of the gdb into the feature dataset since they have the same coordinate system. Unfortunately ArcGIS gives an error dialog that says:               Failed to paste points     The spatial references do not match      Why does this fail? If the GCS of the feature class and the feature dataset are the same, shouldn't it just work? I've verified this is the case via the dialogs as well as exporting the prj file for both the dataset and the class and using a diff tool to compare the two. They are identical. Is the spatial reference of a feature class different than the coordinate system/projection? Attemping to copy the feature class into the feature dataset by holding ctrl while dragging results in ArcCatalog crashing each and every time (shame on ESRI). I guess I could try alternate means of moving feature classes. Use CopyFeatures into the dataset. Project from the feature class into a new feature class within the feature dataset?


title: Modify 'Codepoint with Polygons' to add new fields

text: Sorry if this is a simple question - I'm just beginning in QGIS. I'm using OS Codepoint With Polygons data set which provides shape files of each Postcode area within the UK. Within each postcode area are boudaries of individual postcodes. I'm trying to merge these layers together to make another shape file of our business boundaries that are made up from groups of postcodes. From reading on here, it seems that i need to add a field to the Dbf file to name the business area that i want to 'dissolve' to and fill that field with the business area attributes. Am i barking up the right tree here? Sorry for the long post - Thanks in advance.


title: Transform RGB-Map to CMYK (qgis, opensource tools)

text: I want to print a map I made with qgis. For the print I need CMYK-colors but the print-composers pdf-output is RGB. How can I convert the RGB-map/-pdf to CMYK?


title: How to access selected feature attributes using Python Add-in button?

text: I'm trying to extract the attributes of a selected feature in ArcMap using an add-in button, but I don't know how to refer to this selection in the script. Later I need to copy these attributes into another feature that I want to select using another button and update that feature. Does anyone know how to refer in the script to the selected feature in a layer? Thanks.


title: Run a Python Script when opening an MXD?

text: Can anyone tell me how to accomplish this? I want a Python script to run on opening the MXD which would pull an excel worksheet into a GDB and refresh the map. I only want this to happen in this specific MXD. Thanks


title: Rendering error ArcGIS javascript API IE8

text: I have an issue when displaying features on a map in ArcGIS Javascript API V. 2.8. I want to know if this behavior is seen before, and what could cause it. In Internet Explorer 8 the map looks like this: ![Internet Explorer 8 issue](http://i.stack.imgur.com/Iuqwi.png) In Internet Explorer 9 and Chrome and any other browser I have tested, the map looks like this: ![Chrome screenshot](http://i.stack.imgur.com/bROKN.png) The red color we see is the color of the outline of the polygons. Below is a part of the code for drawing the polygons. As the application is quite complex, as will avoid posting the whole application here. Unfortunately I am not able to provide a direct link to the page, as it is not released yet.               var symbol = new esri.symbol.SimpleFillSymbol(         esri.symbol.SimpleFillSymbol.STYLE_SOLID,         new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_DASHDOT,             new dojo.Color([100, 0, 0, 0.5]), new dojo.Color([100, 100, 100]))     )          ...          dojo.forEach(features, function (feature) {         var graphic = new esri.Graphic(feature.geometry, symbol);         graphic.attributes = feature.attributes;         map.graphics.add(graphic);     });      The flow in the program is basicly this:   1. Page is loaded   2. Map is initialized   3. Polygons are loaded with query task   4. Polygons are drawn on map   5. Event handlers are attached ...


title: Lat Long line intersecting with circle

text: Can anyone help me what's the algorythm to calculate points of intersection whet is given one point p1 [lat,long] , bearing x (let say 137 degrees), arc center [lat,long] and radius? I woulde like to get p2[lat,long] and p3[lat,long]. For better imagination what I mean something like on this picture but on sphere with only one point and bearing: ![enter image description here](http://i.stack.imgur.com/WOFhC.jpg) Thanks in advance.


title: How to write a script to square polygon corners?

text: I'm looking for a script that will square the corners of polygons that are near 90 degrees but not exactly. Near could be dealt with via a variable, for example within 5 degrees of 90 (85 - 95 degrees). For a more detailed description and some rather formal pseudo code, please refer to the following ESRI thread (particularly William Huber's post) : http://forums.esri.com/Thread.asp?c=93&f=987&t=106494


title: In Ogr2OGR: what is SRS?

text: I am trying to understand this command:               ogr2ogr -f "MapInfo File" original.tab wgs84.shp -s_srs EPSG:27700 -t_srs EPSG:4326      I don't understand the s_srs and t_srs paramaters. The help file says that the commands "override the source srs." I am new to GIS. What is SRS? Is this a map making term? It is hard to google.


title: ArcGIS, Calibrating routes, interpolation error

text: I know that copy a question is not exactly great way to ask on this forum, but my situation is exactly the same as this unresolved issue on the (closed) ESRI forum: > I have a stream centreline (fairly curvy) which can be converted to a route. > I have a few points on that line which have stream bed elevation. The bed > elevation varies along the stream, and does not necessarily increase as you > go up the stream (e.g., goes up and down a little bit). > > I want to calibrate the route using the bed elevations, so that I can > interpolate an elevation value between the known points. > > When I attempt to calibrate the stream route with the known points, it > appears to work. But if I examine the measure values of individual vertices > along the route, they are incorrect. Vertices that fall exactly where there > is a known value (a calibration point) have the correct elevation, but > vertices between the known points have incorrect values. The incorrect > values are not 0, which makes me think it is calculating a value somehow. > But it is not doing as I expect. > > For example, between known points of 3.11 and 3.56, you'd expect values > between 3.11 and 3.56. Instead, I might get values of 3.78, 3.78, 3.78, > 3.79, etc. May be it can be resolved in here. Best regards Tomek


title: openlayers 3 creating popups on a click

text: I m very new to ol3. I have a osm map on which i display a geojson file (point). I'd like to create a popup on click for each point to show attributes. I cant find a sort of very basic step by step tuto to do that, even on openlayers-3 examples. Does anyone have some useful links ? thx a lot


title: Setting custom valid values to a MapInfo field?

text: I was wondering if it was possible in MapInfo to set valid values for a field, e.g a field called material and only valid values being BIT, CON or STO? I have a basic knowledge of MapBasic. Thanks, Alan


title: SQL server views in geoserver - problem with sql

text: I try to extract non-spatial data from mssql as spatial data to my Geoserver 2.1.0 installation from a sql statement looking like this:               SELECT  database.table1.id,         geometry::STPointFromText('POINT(' + CAST(database.table1.YCOORD AS VARCHAR(6)) + ' ' + CAST(database.table1.XCOORD AS VARCHAR(7)) + ')', 3008) AS geometry     FROM database.table1      The above given statement returns good spatial info in mssql management studio but I'm not sure how to get it working in geoserver (I have to CAST the data as it is stored as integer). Setting up mssql as a store works fine and I know I can access its tables. I then try to execute the statement as a sql view. That does not work :( The geometry attribute shows up as byte[] and I cant change that. I cant write anything in native SRS and filling in epsg:3008 in declared SRS don't seem to do anything when calculating the bounding box, giving the values MIN X & Y as -1 and the MAX X & Y as 0 in native bounds. * * * @jdeolive - not an option to try trunk at this stage, I have to convince my boss who is a diehard ESRI-fan and stability is a must. @unicoletti - I guess I have to read my statements better not to mix geometry and geography objects :-) thank you. It makes sense when calculating bonds, will try it tonight. Creating a view of it is what I have in mind later on. The only thing is that Im taking my information from another view but as it is only about 6500 points I dont think it will give any larger performance issues.


title: How to Enabling the option to convert layer transparency to color transparency for ArcGIS Server?

text: i have a trouble for publishing map service with transparent color, but when i make analyze map document, there error "warnings",told me (Enabling the option to convert layer transparency to color transparency for ArcGIS Server). ![here the error](http://s3.postimg.org/u8tesunhv/2014_09_06_16_40_46.png)


title: ArcObjects - Adding ArcGIS Server Layer to MXD Programtically

text: I am trying to add a layer from ArcGIS online to a mxd programatically with no luck. See code sample below. The error is happening at pAGSServerConnectionFactory.Open(props, 0); Anyone been able to do this?               private void AddArcGISOnlineLayer(IMap map)             {                 IAGSServerConnectionName servConnName = new AGSServerConnectionNameClass();                 IPropertySet props = new PropertySetClass();                 props.SetProperty("URL", "http://server.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer");                 servConnName.ConnectionProperties = props;                 IAGSServerConnectionFactory pAGSServerConnectionFactory = new AGSServerConnectionFactoryClass();                 IAGSServerConnection pAGSConnection = pAGSServerConnectionFactory.Open(props, 0);                 IAGSEnumServerObjectName pEnumSOName = pAGSConnection.ServerObjectNames;                 IAGSServerObjectName pSOName = pEnumSOName.Next();                      while (pSOName != null)                 {                     if (pSOName.Name == "Layers")                         break;                     pSOName = pEnumSOName.Next();                 }                      IName pName = (IName)pSOName;                 IMapServer mapServer = (IMapServer)pName.Open();                 IMapServerLayer msLyr = new MapServerLayerClass();                 msLyr.ServerConnect(pSOName, mapServer.DefaultMapName);                 IMapServerGroupLayer group = (IMapServerGroupLayer)msLyr;                 ILayer msLayer = (ILayer)msLyr;                 map.AddLayer(msLayer);                  }


title: Using local mbtiles in openlayers?

text: \--The problem-- I've successfully added mbtiles from Mapbox free online source to OpenLayers, but cannot figure out how to get the mbtiles I've created/stored locally to show up in OpenLayers. \--The process-- Style data in Tilemill, export mbtiles, add mbtiles to Tilestream running on localhost, then add to openlayers using the localhost:8888 url. \--The code-- This works:               // Add mapbox mbtiles from web     map.addLayer(new OpenLayers.Layer.TMS('geography-class',     'http://a.tiles.mapbox.com/mapbox/', {     maxResolution: 156543.0339,     type: 'png',     layername: 'geography-class',     isBaseLayer: true     }));      This doesn't work:               // Add tilestream mbtiles local data layer     map.addLayer(new OpenLayers.Layer.TMS('overlayWorld',     'http://localhost:8888/2.0.0/overlay_62ee43/{z}/{x}/{y}.png', {     maxResolution: 156543.0339,     type: 'png',     layername: 'overlayWorld',     isBaseLayer: false     }));      That localhost url comes from the 'TILE URL' field in Tilestream. I've tried multiple different edits to the localhost url but cannot get it right. Any help is appreciated.


title: ESRI Desktop Zoom to the feature

text: i find a feature by using IFeatureClass.Search then i want to show the feature on the screen.How can i do that ? is there a zoom interface or else ?


title: Is there a way to increase the search radius of nearest neighbors for the Nibble tool?

text: I am working with a raster map vegetation types. I have about 20 types with unique raster code values that can be sorted into three classes: 1 - vegetation types I want to expand using the nibble function (e.g. forest) 2 - vegetation types I consider static that I do not want to either nibble or be nibbled (e.g. water, developed areas) 3 - vegetation types that I want to eliminate through the nibble function (early seral stages, such as grass and chaparral) The first time I ran this I was under a deadline and didn't have time to figure out a better solution, and wasn't able to find one googling/searching. I'm hoping to find a way to make this a shorter geoprocessing slog. As inputs to the Nibble, I reclassified my veg layer so that the type 2 (static) values were nodata, since I didn't want them to nibble. That was my input layer. For the input MASK, I reclassified the original raster but made the type 3 values nodata, since they are the ones I want to eliminate through nibbling. The catch is that I tend to have large amounts of type 3 clustered together, and on top of that there are several places where the type 2 and type 3 values are clustered near each other. The upshot is that the nibble gave me an incomplete output. It ran correctly, as far as I could tell, but it would stop before nibbling away all of the masked values. I think the problem is that the nibble function has some predetermined search radius for the nearest neighbor, and if a cell is surrounded by type 2 and 3 cells then it can be so far away from a type 1 cell that the nibble doesn't reach it. However, looking into the tools, the whole Nibble tool is just one function that is called with those two inputs, so I can't modify the python code to change the search radius. I am not skilled enough in python to write a loop that would run the nibble over and over again. So, is there a way to get the Nibble to fill in more of the masked values at a time? Is there a place I can get the Nibble code to try and modify it? I would imagine at some point a distance threshold was specified. Thanks in advance! I hope my explanation was clear! Edit: More information on the data I am using               +-------------+-------------------+--+-------------+--------+-------------+--------+-------------+---------+     |    Name     | Integer in Raster |  |   Static    | Value  |  Dissolve   | Value  |  Nibblers   |  Value  |     +-------------+-------------------+--+-------------+--------+-------------+--------+-------------+---------+     | SMCM        |                 0 |  | URB         |      4 | MCH         |      3 | SMCM        |       0 |     | OCFW        |                 1 |  | BAR         |      8 | GRASS       |      7 | OCFW        |       1 |     | SMCX        |                 2 |  | WAT         |      9 | MCP         |     10 | SMCX        |       2 |     | MCH         |                 3 |  | AGR         |     11 | ADS         |     19 | MHW         |       5 |     | URB         |                 4 |  | MRIP        |     12 | ASPE        |     23 | BLU         |       6 |     | MHW         |                 5 |  | MED         |     13 | ASPW        |     24 | DFTO        |      14 |     | BLU         |                 6 |  | Pixel Count | 322591 | Pixel Count | 333625 | RF          |      15 |     | GRASS       |                 7 |  |             |        |             |        | LPN         |      16 |     | BAR         |                 8 |  |             |        |             |        | SCN         |      17 |     | WAT         |                 9 |  |             |        |             |        | WWP         |      18 |     | MCP         |                10 |  |             |        |             |        | YPN         |      20 |     | AGR         |                11 |  |             |        |             |        | SAGE        |      21 |     | MRIP        |                12 |  |             |        |             |        | CMM         |      22 |     | MED         |                13 |  |             |        |             |        | LSG         |      25 |     | DFTO        |                14 |  |             |        |             |        | Pixel Count | 3744816 |     | RF          |                15 |  |             |        |             |        |             |         |     | LPN         |                16 |  |             |        |             |        |             |         |     | SCN         |                17 |  |             |        |             |        |             |         |     | WWP         |                18 |  |             |        |             |        |             |         |     | ADS         |                19 |  |             |        |             |        |             |         |     | YPN         |                20 |  |             |        |             |        |             |         |     | SAGE        |                21 |  |             |        |             |        |             |         |     | CMM         |                22 |  |             |        |             |        |             |         |     | ASPE        |                23 |  |             |        |             |        |             |         |     | ASPW        |                24 |  |             |        |             |        |             |         |     | LSG         |                25 |  |             |        |             |        |             |         |     | Pixel Count |           4401032 |  |             |        |             |        |             |         |     +-------------+-------------------+--+-------------+--------+-------------+--------+-------------+---------+      I tried the first answer below and it works, but doesn't necessarily solve the problem. Some of my patches of veg that need to be dissolved are pretty large. Therefore, I need an enormous neighborhood to eliminate the whole thing. Given this, it would seem I'm looking at using the raster calculator in an iterative way as well. So, now my question becomes, which is better, and if both, in what order? For example, I could Nibble for a few iterations, then switch to the raster calculator and iterate that over some reasonable-sized neighborhood (not sure how big that would be - 10x10? 50x50) until I eliminated all the type 3 pixels. Or I could do things in the opposite way. Raster Calculator is a little faster, but does it do a better job with respect to landcover assessment if there is a large neighborhood?


title: How to insert a quotation mark through the Field Calculator of attribute table of a Shapefile?

text: I have a shapefile of polylines with information of diameter. But the problem is that i didn't enter the inch sign(") with the values of diameter. What would be the solution that qutation mark would place along with diameter values?


title: How do I force the label to display above the lines

text: How do I force the label to stand out above the red lines to make the label more readable ? Here is the sample of the image ![enter image description here](http://i.stack.imgur.com/C8ujG.jpg)


title: Which CRS fits "Lambert Conformal Conic Projection, Standard Parallels 46.40 and 49.20"?

text: I am trying to georeference a VFR chart in qgis and I can't figure out which CRS to use. Map projection - Lambert Conformal Conic Projection, Standard Parallels 46.40 and 49.20. Any ideas which one on the list will work?


title: Does anyone use both ArcGIS and Python(x,y)?

text: I'm trying to get ArcGIS 10.1 and **Python(x,y)** to peacefully co-exist with a unified installation of Python. Surely someone else is using both these tools... right? Pointers would be much appreciated. Here's where I've gotten to:   1. VERSION: Since ArcGIS 10.1 uses Python 2.7.2 (which contains NumPy 1.6.1), I went back to the pythonxy install mirrors and found version 2.7.2.0 of Python(x,y) which includes Python 2.7.2 and NumPy 1.6.1.   2. INSTALL: I crossed my fingers and just installed Python(x,y) after copying my original "C:\Python27\ArcGIS10.1\" folder to a safe location. It installs most of the files at "C:\Program Files (x86)\pythonxy\", and the Python files in "C:\Python27\". This is the same folder where my Arc install of Python lives, except Arc placed its version in the sub-folder "\ArcGIS10.1\". So to recap, I now have a "python.exe" and associated files in "C:\Python27\" as well as "C:\Python27\ArcGIS10.1\".   3. NEXT? Both Spyder and ArcGIS 10.1 seem to work so far but I'm guessing this won't last long. Two separate installs on the same machine seems like a bad move in principle. I wonder: Is it possible to change Arc's preferences to point to the version in the parent folder, and condense versions? I've gone through folder by folder and in almost every case the Python(x,y) install includes all the files that ArcGIS 10.1 Python does (and then some). The exception is that the Arc install sometimes includes '.pyc' and '.pyo' versions of some '.py' files. But I'm guessing this isn't a critical difference since both directories contain the actual '.py' file.)   4. Forgot to mention: I added "C:\Python27\ArcGIS10.0\" to the list of system paths, per the instructions in Method 2 at http://forums.arcgis.com/threads/34418-is-it-possible-to-install-other-python-distributions-on-top-of-the-arcgis. I really don't know what this does, just experimenting here... Please advise! If anyone can offer further guidance on this, it would be much appreciated!! (PS - I may post the same question at forums.arcgis.com, so if interested in the answer check there too) (PPS - Tried to add 'python(x,y)' as keyword but was unable to create a new keyword)


title: How to add QGIS Python API Intellisense support in Aptana Studio3/Pydev

text: I'm developing QGIS plugins using PyDev in Aptana Studio 3 and I would like to know if it's possible to make it aware of the QGIS API and add some "intellisense" functionnalities when programming. For example, I would like to see the available methods of an object when typing the dot "." after it... UPDATE: I also see this message: > module not in pythonpath When Hovering any QGIS object in code... Is it supposed to show me something more usefull? ![enter image description here](http://i.stack.imgur.com/SJRTz.png) Thanks a lot Etienne


title: Diffusion/Kernel Interpolation with barriers?

text: I need to interpolate water quality data in a rather small area with peninsulas, islands and developed coastline that will act as barriers between my testing/sampling points. I understand I have to use non euclidean method to do so, but does anyone know the difference between "Diffusion Interpolation with Barriers" and "Kernel interpolation with barriers"?


title: Using threading module with Arcpy: Manager/Worker architecture

text: I asked this question on Stack Overflow on how to setup a manager/worker scenario. But I am getting some error when trying to use threading on my Python install with ArcGIS 10 (just the python install from ArcGIS 10.0). Is anyone using threading successfully in arcpy. My end goal is to be able to poll an arcgisjobs directory on ArcGIS Server 10.0 and call worker processes for each new folder. Each new process should be called and set to go and return control back to the manager to call more processes before previous ones are finished.


title: What is the easy way for extracting a dted raster file?

text: i have some srtm data which was converted to dted file with globalmapper or arcmap ascii to raster. i want to extract raster from vector with easy way.. i have tried "extract by mask" method but it depends raster files for i/o. do you know any easy way extracting raster? A section is raster file B section is vector file... **B section is what i want to extract from A Raster File.** **Raster File:** ![enter image description here](http://i.stack.imgur.com/qbqPK.jpg) **Raster and Vector File:** ![enter image description here](http://i.stack.imgur.com/QM01z.jpg)


title: Promotion radius in miles around venues using PostGIS

text: Bear with me as I'm new to PostGIS and GIS in general. I also don't know the domain well enough to create a good title for this post! Here is what I'm trying to accomplish: I have Venues in my system. Each Venue has a promotion radius in miles. Given an arbitrary point (lat/lng), I want to find all the Venues which encompass the arbitrary point in their promotion radius. Here's what venues tables look like:               venues     id | name | latitude | longitude | promo_radius      Is this possible with PostGIS? Are there any resources that would help me further understand the query? Thanks!


title: Explode MULTILINESTRING into individual segments in PostGIS 1.5

text: I have a `MULTILINESTRING` in PostGIS consisting of several `LINESTRING`'s. Consider each `LINESTRING` to have a set of `POINT`'s > P1, P2, ... PN forming the line segments > L1, L2, ... LN-1 where LX is the line segment connecting PX to PX+1. **Using PostGIS 1.5 how can I explode the`MULTILINESTRING` into individual line segments of each `LINESTRING`?**


title: Geo and Population Coverage Calculation

text: We are calculating Population and geocoverage in mapinfo. We have calculated population coverage as an area overlap. **How do we calculate population coverage if we have population defined for each polygon?** ![coverage](http://i.stack.imgur.com/I4ju6.png) **1.** **Are there any particular method to calculate population coverage for each boundary?** **2\. Any documentation and detail for analysis (not only steps but reason) would be helpful.**


title: How to dynamically add all feature layers using ArcGIS JavaScript API?

text:
I have a mapserver > http:// mymapserver/services/rest/xxx/mapserver/ which contains 7 layers with /mapserver /1 /through mapserver/7 If a define an ArcGISDynamicMapServiceLayer and 7 FeatureLayers in OnDemand Mode, everything gets displayed and I can query every layer based on the infoTemplate. There is however one small glitch: the ArcGISDynamicMapServiceLayer has 1-x default visible layers which are not queryable and have a slight offset. Is there a way, to define the ArcGISDynamicMapServiceLayer, retrieve all included FeatureLayers and add them on the fly to a) the basemap and b) the Layerlist, while retrieving c) the Legend? Bonuspoints for keeping the hints borderline simple; I am a javascript novice and I hate dojo. Kind regards this is the actual code:               var visible = [];     var ml = [];          require(["dijit/layout/BorderContainer", "dijit/layout/ContentPane", "dijit/layout/AccordionContainer", "esri/dijit/Legend", "dojo/_base/array", "dijit/form/Button"]);     require(["dgrid/OnDemandGrid", "dgrid/Selection", "dojo/store/Memory", "esri/map", "esri/layers/FeatureLayer", "dojo/_base/declare", "dojo/number", "dojo/parser", "dojo/domReady!", "dojo/on", "dojo/dom", ],     function (Grid, Selection, Memory, Map, FeatureLayer, declare, dojoNum, parser, array, on, dom) {              parser.parse();         window.infoTemplate = new esri.InfoTemplate("Attributes", "${*}");         window.map = new esri.Map("map", {             basemap: "topo",             center: [12.8, 51],             zoom: 9         });              on(dom.byId("btnStreets"), "click", function () {             window.map.setBasemap("streets");         });         on(dom.byId("btnSatellite"), "click", function () {             window.map.setBasemap("satellite");         });         on(dom.byId("btnHybrid"), "click", function () {             window.map.setBasemap("hybrid");         });         on(dom.byId("btnTopo"), "click", function () {             window.map.setBasemap("topo");         });         on(dom.byId("btnGray"), "click", function () {             window.map.setBasemap("gray");         });         on(dom.byId("btnNatGeo"), "click", function () {             window.map.setBasemap("national-geographic");         });              window.MapServiceURL = "http://theserver/ArcGIS/rest/services/xxx/MapServer/";         window.HWSUrl = "http://theserver/ArcGIS/rest/services/xxx/MapServer/1";         window.outFields = ["OBJECTID", "HOCHWERT", "SchadenNr", "STATUS"];              var fl = new FeatureLayer(window.HWSUrl, {             id: "schaden",             mode: esri.layers.FeatureLayer.MODE_ONDEMAND,             outFields: window.outFields,             infoTemplate: infoTemplate         });              /** excluded FeatureLayers 2-7 **/              ml = new esri.layers.ArcGISDynamicMapServiceLayer(window.MapServiceURL);              dojo.connect(fl, "onLoad", function (fl) {             fl.maxScale = 0;             fl.setSelectionSymbol(new esri.symbol.SimpleFillSymbol().setOutline(null).setColor("#AEC7E3"));         });              dojo.connect(map, "onLoad", function (map) {             // show the border container now that the dijits              // are rendered and the map has loaded             dojo.style(dijit.byId("container").domNode, "visibility", "visible");         });              //add the legend         dojo.connect(map, 'onLayersAddResult', function (results) {             var layerInfo = dojo.map(results, function (layer, index) {                 return { layer: layer.layer, title: layer.layer.name };             });             if (layerInfo.length > 0) {                 var legendDijit = new esri.dijit.Legend({                     map: map,                     layerInfos: layerInfo                 }, "legendDiv");                 legendDijit.startup();             }         });              if (ml.loaded) {             buildLayerList(ml);         } else {             dojo.connect(ml, "onLoad", buildLayerList);         }         map.addLayers([fl]);                    });          /*building LayList function*/     function buildLayerList(layer) {         var items = dojo.map(layer.layerInfos, function (info, index) {             if (info.defaultVisibility) {                 visible.push(info.id);             }             return "<input type='checkbox' class='list_item' " + (info.defaultVisibility ? "checked=checked" : "") + "' id='" + info.id + "' onclick='updateLayerVisibility();' />" + info.name + "
";         });         dojo.byId("layer_list").innerHTML = items.join(' ');              layer.setVisibleLayers(visible);         map.addLayer(layer);     }          function updateLayerVisibility() {         var inputs = dojo.query(".list_item"), input;              visible = [];              dojo.forEach(inputs, function (input) {             if (input.checked) {                 visible.push(input.id);             }         });              if (visible.length === 0) {             visible.push(-1);         }         ml.setVisibleLayers(visible);     }


title: How to animate raster time series using open-source GIS?

text: I have a bunch of rasters taken about 10 minutes apart, and I have a nice QGIS project that produces pretty enough maps of them. I'm wondering if there's a way to animate this raster series, with a color map and overlaid with (static) vector data, using open-source GIS. A QGIS solution would be ideal since I already have the project set up there, but other suggestions are welcome as well.


title: GeoEyes stereo Image

text: We have 4 stereo GeoEyes images. We work with ErdasImagine 2014 when we generate a blk file after reading the images we see this error: "Unsupported Compression type". Please help us. Best Regards


title: Points appear at different location when using google map overlay

text: I have used vector point from database and also polygon country vectors from database.Points appear with in this country vectors but as I switch to google earth map as base layer the points at smaller zoom level would appear outside boundary where as at other zoom level it would appear inside. I have tried to add projection too but in vain.Do I need to reproject google earth's base layer and How?               map = new OpenLayers.Map('map', {                         controls: [                             new OpenLayers.Control.Navigation(),                             new OpenLayers.Control.PanZoomBar(),                             new OpenLayers.Control.LayerSwitcher({'ascending':false}),                                                    new OpenLayers.Control.ScaleLine(),                             new OpenLayers.Control.MousePosition(),                             new OpenLayers.Control.OverviewMap(),                             new OpenLayers.Control.KeyboardDefaults()                         ],                                   projection: new OpenLayers.Projection("EPSG:900913"),                  displayProjection: new OpenLayers.Projection("EPSG:4326"),                  numZoomLevels: 21,                  maxExtent: new OpenLayers.Bounds(-20037508, -20037508,20037508, 20037508.34)                          });                                     var border = new OpenLayers.Layer.WMS( "Border",                     url,                      {layers: 'GisData:BORDER'}                 );                          var dsBand = new OpenLayers.Layer.WMS( "District Band",                     url,                      {layers: 'GisData:BND'}                 );                                                 var Site1 = new OpenLayers.Layer.WMS( "Site 1",                     url,                      {layers: 'GisData:sites', transparent: "true",format: "image/png"}                              );                     dwnSite3G.setIsBaseLayer(false);                          var Site2 = new OpenLayers.Layer.WMS( "Site 2",                     url,                      {layers: 'GisData:sites2', transparent: "true",format: "image/png"}                              );                               var g_pmap = new OpenLayers.Layer.Google(                     "Google Physical",                     {type: G_PHYSICAL_MAP}                     );                          var gmap_gs = new OpenLayers.Layer.Google(                         "Google Streets" // the default                          );                          var gmap_ghy  = new OpenLayers.Layer.Google(                         "Google Hybrid",                         {type: G_HYBRID_MAP}                     );                          var gmap_gsat = new OpenLayers.Layer.Google(                         "Google Satellite",                         {type: G_SATELLITE_MAP}                     )


title: QGIS - Export PDF map with georeferenced data

text: I am used to working with ArcMap but have been playing around with QGIS. When I export a PDF map with ArcMap, the PDF is spatially referenced. This is a handy tool, I can put the PDF map on my iPhone and it tracks my location on the map (even when I'm out in the woods away from cell coverage). It's very handy for forestry work. Are the PDF maps exported with QGIS also spatially referenced? Is there a plugin that allows for this feature?


title: Using SearchCursor to access & export values in a raster table

text: I am desperately trying to export the values in a raster table to a .txt file. I've seen some other folk post about this, but for some reason can't get the code to work in my case. I have a raster, with table already built, called "shed". I simply want to export two fields, "Value" and "Count", to a .txt or .csv file (called "try" in this case). So far, I haven't even been able to get a script that even looks at the raster's table. So far I have:               outfile=open("E:/Weathering_GIS/try.txt",'w')     rows=arcpy.SearchCursor('shed',"","","Value;Count","")     for row in rows:         val=row.getValue('Value')         count=row.getValue('Count')          outfile.write(val)      But "val" and "count" always fail to become anything- they just remain empty values! Am I missing something? Maybe some sort of import I have overlooked?


title: SVG iconed features not rendered correctly during zoom

text: Im using Geoserver with PostGIS database to put points to map using GeoExplorer. The points are styled to use SVG icon with fixed size (uom metre, size) so they scale together with the map as user zooms in and out. For now, there are no filters or any other issues. Style code :               <sld:PointSymbolizer uom="http://www.opengeospatial.org/se/units/metre">             <sld:Graphic>               <sld:ExternalGraphic>                 <sld:OnlineResource xmlns:xlink="http://www.w3.org/1999/xlink" xlink:type="simple" xlink:href="load_icon.php?id=${param}"/>                 <sld:Format>image/svg+xml</sld:Format>               </sld:ExternalGraphic>               <sld:Size>80000</sld:Size>             </sld:Graphic>           </sld:PointSymbolizer>           <sld:TextSymbolizer/>      load_icon.php is returning the SVG code for given param :               <svg xmlns="http://www.w3.org/2000/svg" version="1.1" width="220" height="110">     <g transform="scale(0.5)">     <g transform="translate(90,50)">                     <rect x="1.5" y="1.5" height="103" width="167" fill="#fff" stroke="red" stroke-width="3"/>     <path d="M3,3 l164,100 m0,-100 l-164,100" stroke="red" stroke-width="3" fill="none"/>     </g>     <text x="173" y="49" text-anchor="middle" font-family="Verdana" font-size="25" fill="red">XXXXX< /text>     <text x="270" y="153" font-family="Verdana" font-size="25" fill="red">G< /text>     </g>     </svg>      I use the plain GeoExplorer to view the icon on the map, but same issue is when connecting the layer as WMS to Leaflet and other libraries. The icon is showing ok in some resolution : ![icon ok](http://i.stack.imgur.com/zsgpN.png) but gets cropped as user is zooming in : ![icon croped](http://i.stack.imgur.com/MooMk.png) and even more cropped with further zoom : ![icon more cropped](http://i.stack.imgur.com/zupKx.png) with tiles that do not contain the part of the icon. There is caching turned off for this layer, no help. What I am doing wrong here ?


title: accessing mapinfo data without using mapinfo api

text: I need to be able to access various mapinfo layers but without using the mapinfo api. I have found that I can convert the mapinfo layers to an access database to get the attribution and I can use sharpMap to load and view the layer. The problem is now I'm not sure how to get a location shown in sharpmap and then query the objects to get to the attributes, or if its even possible. It seems that mapinfo might keep this kind of link internaly perhaps


title: How to install fTools and Table manager

text: I've just make clean install of Qgis Lisboa under Ubuntu Precise. Latest Qgis version under Windows XP is very unstable (crashes every few minuter while editing shp file). After new installation there are some problems for me:   1. I can't see or add independent repositories. So I can't add them under "Repositories" tab. Can I add them anyway? What I need is Table Manager at least. Possibility of install some other plugins could be usefull as well.   2. I can't find "Join attributes" under fTools menu (it was the same problem after upgrading to Lisboa under Win XP, but deleting fTools directory and reinstalling this plugin solve problem but only under Win XP). So, I can't work under Windows due to crashes nor under Ubuntu :( Could you help me to make Qgis full functional under Ubuntu? TIA Best regards, Cyfranek


title: Python: shp file to GML using OGR

text: I am trying to convert a shapefile to GML 3.1.1 . I am using gdal 1.9.1 however the GML conversion it generates for a geometry feature is GML 2 and it also swaps the long/latitude coordinates. Are there simple workarounds to keep the coordinate positions the same, while changing the output to GML 3.1.1? Thanks.


title: How to create TINs in ArcGIS from triangular elements

text:
I am working with FEM models with triangular elements, and I'd like to view the results as TIN surfaces in ArcGIS (and in ArcScene). In my external source, I have **nodes** (with X,Y,Z coords) and **elements** (forming a triangle from three nodes), which are the fundamentals of a TIN. How can I bring these together in ArcGIS to convert my triangular elements into a new TIN surface? This should ideally be a 1:1 conversion, but the documentation isn't helping me. Here is an illustrated example of the data I have: ![Simple TIN](http://i.stack.imgur.com/TaS5A.png)               Nodes: n, X, Y, Z     1, 51.06, 84.98, -7.46     2, 30.45, 80.67, -7.47     3, 44.83, 72.70, 19.46     4, 64.99, 69.91, 61.48     5, 30.33, 63.40, 71.31     6, 46.95, 61.96, 48.44     7, 62.04, 51.97, 82.51     Elements: n1, n2, n3     6, 7, 4     4, 3, 6     3, 5, 6     5, 3, 2     1, 3, 4     2, 3, 1      I'm currently using ArcGIS 9.3 with the 3D Analyst Extension (among others), and would like some hints for a solution targeted for either Python or VBA. * * * _**Alternative approach_** : The raw data -> Esri TIN conversion doesn't look very easy, and I can't make heads or tails of the API for 3D Analyst, so I'll hold the Python/VBA solution off. I have the same data represented in a Shapefile where all features are triangular PolygonZM shapes and the vertices have the height-field. (If you ask, this Shapefile was made using OGR via Python from the original dataset structured similarly as above). Viewing this file in ArcScene is really really slow, as I have ~100k triangular polygons in the Shapefile. This PolygonZM Shapefile should ideally convert triangle-by-triangle to a TIN dataset. I have the data, so interpretation or interpolation is not wanted! So, how can I convert this PolygonZM Shapefile (pictured below) into a TIN? In ArcScene, in the "3D Analyst" menu, there is a "Create TIN From Features" tool, but I don't understand the language "soft clip" etc. The dialog provides no help. ![ArcScene: Create TIN From Features](http://i.stack.imgur.com/qaCeA.png) * * * # LandXML It turns out the best format is LandXML, thanks to the related answers below. For the example above, the file looks like this (note: the coordinate order is Y X Z):               <?xml version="1.0" ?>     <LandXML version="1.2" xmlns="http://www.landxml.org/schema/LandXML-1.2" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.landxml.org/schema/LandXML-1.2 http://www.landxml.org/schema/LandXML-1.2/LandXML-1.2.xsd">       <Surfaces>         <Surface name="My TIN">           <Definition surfType="TIN">             <Pnts>
84.98 51.06 -7.46

80.67 30.45 -7.47

72.70 44.83 19.46

69.91 64.99 61.48

63.40 30.33 71.31

61.96 46.95 48.44

51.97 62.04 82.51
             </Pnts>             <Faces>               <F id="1">6 7 4</F>               <F id="2">4 3 6</F>               <F id="3">3 5 6</F>               <F id="4">5 3 2</F>               <F id="5">1 3 4</F>               <F id="6">2 3 1</F>             </Faces>           </Definition>         </Surface>       </Surfaces>     </LandXML>      Then use the LandXMLToTin_3d tool to process the result to an Esri TIN file: ![Resulting TIN](http://i.stack.imgur.com/SmMke.png)


title: How can you make RasterToASCII_conversion observe environment setting output cell size?

text: I have two rasters (an ESRI grid and an IMAGINE image) that differ in projection, extent and cell size. I want to convert the IMAGINE to an ASCII, and want it to have the extent, cell size and projection of the ESRI grid. I've set global environment variables (extent, output CRS, cell size, snap raster) to match the values in the ESRI grid, and the ASCII that results from RasterToASCII_conversion has the correct projection and extent, but ignores the cell size setting. However, if I use the Raster and ApplyEnvironment tools (see code below) from the Python console, the output is as expected (correct cell size, extent, etc).               import arcpy     arcpy.env.extent = 'c:/path/to/grid'     arcpy.env.cellSize = 'c:/path/to/grid'     arcpy.env.outputCoordinateSystem = 'c:/path/to/grid'     arcpy.env.snapRaster = 'c:/path/to/grid'     r = arcpy.Raster('c:/path/to/imagine')     arcpy.CheckoutExtension('Spatial')     rclip = arcpy.sa.ApplyEnvironment(r)     rclip.save('c:/path/to/outfile')      While I could do the above and then convert the resulting grid to ASCII, I'd prefer to do the whole procedure without using Spatial Analyst. Is this possible? I'm using ArcGIS 10.


title: How to "lower the resolution" of a lat long csv

text: I have a csv with 12.5 million date/lat/long/value points, it's a point every 800 meters ( in a grid ) of the united states what's the best way to convert the points to every 1600 meters? or 3200? etc... Really need to reduce the file size


title: Spatially enabled Postgres database

text: I am trying to create a spatially enabled PostGIS database. I am following the PostGIS documentation, http://postgis.net/docs/manual-1.5/ch02.html#id2648455. In the short version, i executed,               createdb yourdatabase     createlang plpgsql yourdatabase      But I did not execute the last three commands,               psql -d yourdatabase -f postgis.sql     psql -d yourdatabase -f postgis_comments.sql     psql -d yourdatabase -f spatial_ref_sys.sql      Is my database still spatially enabled? If not, then can you tell me where these .sql files are located so I can execute those last commands, since I am getting the following error,               postgis.sql: No such file or directory      I am using Ubuntu 12.04


title: Qgis 1.8 Lisboa does not have Join attributes (only Join attributes by location)

text: I installed 1.8 Lisboa. I tried to join to the attribute table more data (more columns ), but Vector/Datamanagement Tools/ Join attributes by location there were also not "Join attributes" alternative. b) I installed 1.4 and there was Vector/Datamanagement Tools/ Join attributes 2) what to do, when attribute table has only Nulls, how to edit, the edit button does not work nor 1.4 and 1.8 operation system: Windows XP IBM ThinkPad T60


title: How to rotate a feature in openLayers

text: I have an svg image and I want to rotate it around itself so I do the following:               var origin = new OpenLayers.Geometry.Point(vertices[0].x,vertices[0].y);                 feature.geometry.rotate(90 ,origin);     feature.layer.drawFeature(feature);      but nothing happens. Does anyone know why? Thanks in advance.


title: Is there a workaround for 80+ layers in OpenLayers?

text: My Map web app reads client list from SQL server and for each client it draw points on a different layer so that the user could turn on and off that layer and see his entire customers. Problem is that once I have around 80+ layer than popups are appearing below the layers. Is there a workaround how I can show different features and turn them on / off without screwing with the popups ?


title: How to cache PostGIS views?

text: Does anyone knows if there exists a caching system for views or queries in postgis? I have a system that gives information through a webpage using a postgis view as backend. The problem is that, sometimes the view is big enough that takes 5 seconds per user to serve the info. I know there's this system named geocache that works very good for raster data using wms protocol. I was wondering if there's a similar tool for vector / table data.


title: How to convert polygon layer to txt files?

text: Is there any way to convert polygons that do not have coordinates in the attribute table to a txt file with coordinates. I will be using the txt file in a different software (not ArcGIS), so I need the txt file with polygon coordinates. Hope someone can help me. Thank you in advance. Lina


title: Retrieving routedescription geometry as GeoJSON

text: I'm building a small web-map and would wich to display a polyline representing a specific route via roads. With leafletjs it is easy to display the geometry on/over the map. So that's not the problem. I'm looking for a webservice where I can 'download' the route as a polyline/Multipoint. Does anyone know if this exists?


title: Can esri ipad app consume web mapping service?

text: Can the ESRI app for the iPad consume regular web mapping services, or do I need to publish services for this app as Mobile project services?


title: When a class in an ArcGIS extension should be an ArcGIS class

text: Visual Studio templates provided by the ArcGIS SDK offer a "ArcGIS Class" which is a simple C# class that has the attributes Guid, ClassInterface, and ProgId before the class declaration.   While adding a new class to an ArcGIS extension project, when should I use an ArcGIS class and when should I use a regular C# class?


Traceback (most recent call last):       File "/home/kasun/.qgis2/python/plugins/floogin/floogin.py", line 246, in findRoutes         path=self.astar.astarSearch(graph, start, end)       File "/home/kasun/.qgis2/python/plugins/floogin/astar.py", line 110, in astarSearch         for tile in graph[current]:     TypeError: 'QgsGraph' object does not support indexing      Above Error occured when I run the following code. I'm new to qgis and I don't know why the error occurs.I wanted to find path between two selected points.I would really appreciate if you could please help me out to find solution to the problem. Thank You.               def astarSearch(self, graph, current, end):         openSet = set()         openHeap = []         closedSet = set()              def retracePath(c):             path = [c]             while c.parent is not None:                     c = c.parent                     path.append(c)             path.reverse()             return path              openSet.add(current)         openHeap.append((0,current))         while openSet:             current = heapq.heappop(openHeap)[1]             if current == end:                     return retracePath(current)             openSet.remove(current)             closedSet.add(current)             for tile in graph[current]:                     if tile not in closedSet:                         tile.H = (abs(end.x-tile.x)+abs(end.y-tile.y))*10                          if tile not in openSet:                                 openSet.add(tile)                                 heapq.heappush(openHeap, (tile.H,tile))                         tile.parent = current         return []               def findRoutes(self):     vl=self.canvas.currentLayer()     director = QgsLineVectorLayerDirector( vl, -1, '', '', '', 3 )     properter = QgsDistanceArcProperter()     director.addProperter( properter )     crs = self.canvas.mapRenderer().destinationCrs()     builder = QgsGraphBuilder( crs )          global x1     global y1     global x2     global y2     pStart = QgsPoint( x1, y1 )     pStop = QgsPoint( x2, y2 )          tiedPoints = director.makeGraph( builder, [ pStart, pStop ] )     graph = builder.graph()     tStart = tiedPoints[ 0 ]     tStop = tiedPoints[ 1 ]     idStart = graph.findVertex( tStart )     idStop = graph.findVertex( tStop )          start, end = idStart, idStop     path=self.astarSearch(graph, start, end)     if path is None:             print "No path found"     else:         print "Path found:", path


title: QGIS is not working with Postgres Plus Advanced Server

text: We have a problem to integrate GIS software with our GIS database which is populated on Postgres Plus Advanced Server 9.2 from EnterpriseDB. We are using QGIS to maintain and analyze our GIS repository, the software is working just fine when the database in the postgres community version, but when the database in ppas; QGIS return the following error while trying to fetch the list of GIS layers: > Erroneous query: >   >   >      SELECT > l.f_table_name,l.f_table_schema,l.f_geometry_column,upper(l.type),l.srid,c.relkind >      FROM geometry_columns l,pg_class c,pg_namespace n >      WHERE c.relname=l.f_table_name >        AND l.f_table_schema=n.nspname >        AND n.oid=c.relnamespace >        AND has_schema_privilege(n.nspname,'usage') >        AND > has_table_privilege('"'||n.nspname||'"."'||c.relname||'"','select') >      ORDER BY n.nspname,c.relname,l.f_geometry_column >   > > returned 7 [ERROR: schema "dbms_profiler" does not exist] And the following is the answer from EnterpriseDB: > We looked at query and it seems, issue is the way that query built. By > default, PPAS stores package names in namespace table and so it is > retrieving the names of packages as name of schemas hence query is failing. > It works in PostgreSQL. Though, you are in postgres compatibility, we still > retains some packages. > > So we would recommend you to contact QGIS vendor on the same as it is not > tested with PPAS. Your help will be highly appreciated.


title: Flannery compensation in Leaflet JS

text: I'm creating a bubble map using Leaflet JS. Everything is working smoothly (LeafletJS is awesome!) but now I'd like to try flannery compensation on the bubbles. I can't find the math anywhere though. Best explainations I've been able to find say you have to multiply radius by either 1.4 or 1.8, which just seems a bit useless. What is the math used for Flannery compensation?


title: Snag loading Gdal 1.11 for Qgis 2.2, on a Mac

text: Can't load Q2.2 on a Mac without Gdal 1.11. but when I try to load 1.11 'Complete' version, it fails and gives me an error message to contact the developer. Anyone else had this kind of trouble? David A UK


title: csv to INSPIRE gml problem in FME Desktop 2014

text: I have .csv file with id, name, lat, long and other. I just cannot make gml writer to write geometry tags. I mean something like . How can I read this things from my .cvs? What am I missing? fme is new for me... Thank you.


title: How to extract vectors from map using QGIS?

text: I am trying to extract building footprints from a raster map (as shown in the example below). Please note that the maps are out of copyright, have been scanned and georectified by a library and have been provided to me for exactly this purpose. ![Scaled example - original geotif available on request](http://i.stack.imgur.com/fQotP.jpg) So far I have been using the NYPL Map-Vectorizer which works by first adjusting the brightness/contrast of the image, then setting a colour threshold to extract a black and white image, which in turn is processed using gdal_polygonize.py (followed by some simplification and removal of roads). The issue is that this gives polygons of the _inside_ of the building. Terrace buildings end up with a gap between them. **How can I post-process to remove the gap between the terrace buildings? Is there an altogether better method for extracting building footprints from this map?** Note: Solutions need to be able to run as a script. I'm able to test proposed techniques in QGIS.


title: Automating Intersect between two polygon layers

text: I have a spreadsheet that has county-level occurrences of around 1000 species within the US in the column headings (Species State County). Thus, for some counties there may be multiple species. I also have a county shapefile with the County field linking the shapefile to the spreadsheet. Finally I have a shapefile with multiple, separate polygons and I want to know which rows from the spreadsheet intersect or contain a polygon from the second polygon shapefile. I am looking for the best way to automate this process so that I create a new column in the spreadsheet that has a 1 if the county occurrence of the species intesects with the second polygon shapefile. For example, Species State County Intersect sp. A NC Durham 1 sp. A NC Orange 0 where species A county-level occurence intersects with the polygon in Durham County, but not in Orange County. I have not specified a platform, but I have worked with R, ArcGIS, and Python. Currently I am most comfortable working in R (but not for spatial data, there I am a newbie), but not sure how to automate this process. thanks for any information.


title: Weighting a data frame within another frame for map creation

text: I've currently got an inset data frame sitting in a corner of a larger data frame. Some labels are being half covered by the smaller data frame. I was wondering if there are any tips on how to weight it (either in Maplex or Standard) so that the labels that are half covered are either removed, or moved so that they aren't behind the smaller data frame? Thanks


title: Remove overlap in raster using ERASE command line?

text: I am trying to erase a raster and is there a erase command for the Raster side? I am using ArcGIS 9.3.1. I have 2 rasters that are in the Landcover and there is overlap in between and would like to remove the overlap and then merge them later. Thanks, Robert


title: Waiting for Vector & WFS loading

text: _user borrowed from friend_ since i use vector layer with wfs protocol, there is a little while before features loaded on the map. So i need to tell user to wait a moment until they load. what best way to do that ? i found this http://osgeo- org.1560.n6.nabble.com/Patch-quot-waiting-quot-mouse-cursor-on-WFS-loading- td3961983.html but not sure if it will work. any body tried it ?


title: Trying to make the private land white overlay hillshade but

text: I'd like to have anyone's feedback on this draft map I created. Tried to make the private ownership white but since it's RGB is not near to 255. I've seen someone created this map similiar to it and leaves it as light grey color. Thanks. http://www.cartotalk.com/index.php?showtopic=7125 ![Gray Image](http://i.stack.imgur.com/RqiR0.jpg) I am using ArcMap 9.3.1. I have Photoshop 7 on my computer but no training background on it. Have use couple of times in the past not related to GIS.


title: Is it possible to use a Multipatch to implement IFeatureRenderer.Draw in a custom renderer

text: I have a custom renderer which need to render modified geometry. The input features are all coming in via the `IFeatureCursor` as polygon elements. (This is similar to this example, though I'm not rendering points.) I can easily render Polygons, Lines, and Points via either passing the new geometry into `IFeatureDraw.Draw` or calling `IDisplay.DrawPolygon`/`IDisplay.DrawPoint`/etc. However, `IDisplay` does not provide a mechanism for drawing a multipatch, and passing this in as the geometry to `IFeatureDraw.Draw` results in nothing drawn to the screen. Is there a way to directly "draw" a multipatch element (preferrably with texturing), within this method?


title: Does a geotagged image contain information about its accuracy?

text: When i view my current location on Google maps on my android it shows me a circle of uncertainty. If i take a geotagged picture it will off course contain information about its Lat or Lon. Does it also contain information concerning accuracy of this geotagging ?


title: Toggle editing option in plugins

text: I am working on python plugins for QGIS.I developed the plugins same as RT SQL LAYER plugins. When we draw the layers using query builder of " **Add Postgis Layer** " option, **toggle editing** option is available for the same. But when i try to draw layers using **RT SQL LAYER** plugin ,it is not available for editing. Same happens with my plugins too. How do I make **toggle editing** option available for my plugins.?? My plugins are in **python**.


title: extract vertices from dbgeometry value

text: Who knows the elegant way to extract coordinates of polygon vertices from well known text representation of DbGeometry value? Here's my polygon for example:                {POLYGON ((52.808019 -1.345367, 53.095124 -1.345367, 53.095124 -0.864716, 52.808019 -0.864716, 52.808019 -1.345367))}      Thanks in advance


title: New Install Quantum GIS - .ind and .idm files not valid

text: I downloaded and installed the Windows version of Quantum GIS to view Pinellas County data. I can add the .shp, .shx, and .dbf layers, however do not see any data showing as far as road names, etc. I tried to add the .idm and .ind layers that come with the Pinellas County downloads and they show as not valid. Do I need those to see the data on the map? If so, how do I add them.


title: django + apache mod_tile

text: I'm not sure if I'm doing this right and I need your opinion on this current situation. I'm creating a Web Interface using django and postgres as the database. With regards to the slippy map, I went with using apache and mod_tile. The tiles where first designed in tilemill and then used the generated xml file with the generate_tiles.py. Am I doing this right or are there other alternatives? Just in case, the web interface will be used to draw lines on the map as routes and store them for querying.


title: How to install SciPy for QGIS Standalone

text: A number of QGIS plugins require the Python package SciPy to be installed. This package is available using the OSGeo4W package based installer or precompiled for Windows, but is not provided as part of the Windows standalone installer. The SciPy windows installer does not recognize the Python instance installed with the QGIS standalone install, and so will not install to it. Can anyone provide a simple solution to installing the SciPy package into the QGIS standalone python instance that is suitable for an inexperienced python user? Any help would be really appreciated.


title: How to generate StreetFrom and StreetTo data from street segments?

text: ## **Background:** I have a polyline feature class of about 4000 street centerline segments for which I would like to add data on what streets the segments intersect with. I would like to add some new fields to my dataset that capture information on what street the polyline segment begins at and what street the polyline segment ends at (essentially the intersections) ## **What I have so far:** I wrote a quick Python script to calculate these values by doing Select By Location -> Boundary Touches:               # Import arcpy and set path to data     import arcpy     arcpy.env.workspace = "c:\Users\me\Documents\ArcGIS\Default.gdb"          #create a master search layer to encompass all streets.      searchLayer = arcpy.MakeFeatureLayer_management('StreetsCP', 'searchLayer')      #set the fields that we are going to be updating     fields = ('FULL_NAME', 'OBJECTID', 'StreetFrom', 'StreetTo')     #create a cursor to encompass all rows in streets     with arcpy.da.UpdateCursor('StreetsCP', fields) as uCursor:             #loop over every row in streets             for uRow in uCursor:                 #set where clause to create a feature layer from the row                   whereClause = '"OBJECTID" = ' + str(uRow[1])                 #create feature layer consisting of just the row we are on                 inputLayer = arcpy.MakeFeatureLayer_management('StreetsCP', 'inputLayer', whereClause)                       #select the streets that touch the boundary of the street we are on                 arcpy.SelectLayerByLocation_management(searchLayer, "BOUNDARY_TOUCHES", inputLayer, "", "NEW_SELECTION")                      touchingStreets = []                 uniqueStreets = []                      #populate the unique streets list                 with arcpy.da.SearchCursor(searchLayer, 'FULL_NAME') as sCursor:                     for sRow in sCursor:                         if sRow[0] != uRow[0]:                             touchingStreets.append(sRow[0])                     #eliminate duplicate street values                     uniqueStreets = list(set(touchingStreets))                 if len(uniqueStreets) <= 2:                       i = 2                     for street in uniqueStreets:                             uRow[i] = street                          i += 1                       uCursor.updateRow(uRow)                      arcpy.Delete_management("inputLayer")      ## **The Problem:** The issue I am running into is the directionality of the street. In my street data I have address directionality using Left From/Left To/Right From/Right To, like in the below table: ![Example Table](http://i.stack.imgur.com/qpsg7.jpg) The Street From and Street To should reflect this same directionality. **What are some ideas for capturing the directionality of the street?** ## **Some ideas I've been exploring:**   * Is there something exposed through ArcObjects that stores the endpoint and beginning point of a polyline using these metrics that I can use to select by location?   * Can I somehow determine the directionality just in ArcPy?   * As a last resort, Network Analyst might be able to offer me this functionality, but due to not having used Network Analyst before the ramp up time may be too much for my purposes.   * If there is functionality that may help me in GRASS/QGIS I'd be willing to explore that as well.   * I also have an "intersections" point file that could be used.


title: Rendering a QGIS WMS layer in Firefox hogs the machine

text: I am less than a novice beginner on GIS Systems and am groping for answers. I followed instructions for WIN 7 Install at both http://underdark.wordpress.com/2012/04/06/qgis-server-on-windows7-step-by- step/ and http://hub.qgis.org/projects/quantum- gis/wiki/QGIS_Server_Tutorial#Windows All appears to work well. I can retrieve a capabilities-XML. Next, I add an WMS Layer to a shapefile in QGIS. Again, no issues. Then I attempt to display via Firefox 18.0.2 with address call > > [http://localhost]/qgis/qgis_mapserv.fcgi.exe?SERVICE=WMS&VERSION=1.3.0&SRS=EPSG:4326&REQUEST=GetMap&map=C:/OSGeo4W/apps/qgis/bin/VA_VOTES_BUSINESS.qgs&BBOX=-36.215710,6.817432,-9.700450,58.551008&WIDTH=500&HEIGHT=1000&LAYERS=VA_VOTES_BUSINESS&FORMAT=image/png The entire operating system becomes unstable at the page call. Firefox displays `Connected to LocalHost`; there is 100% CPU usage (Pentium Dual Core 2.2 GHz); and 80% memory usage (4 GB). I have to power off to regain control. I see nothing but errors in the Apache error log but of course I'm not familiar with the references. Spent 10 hours+ trying to troubleshoot and this is my last resort. Any direction would be most helpful and appreciated. QGIS 1.8.0 Lisboa Desktop


title: Unable to publish geocoding service [ERROR: 001369]

text: I manage to create a geo locator on arcmap and tested it. but when I publish to arcgis server using arcmap or arcgis server manager. I encounter the error 1) Fail to publish service. The error is very generic on logs. may I know is there anyway to to troubeshoot this issue?


title: What might slow Inkscape down when working with QGIS SVG or PDF output?

text: When I import an SVG or PDF exported from QGIS print composer (2.2.0) into Inkscape 0.48 the vector graphic has large black lines and squares. These change size when I zoom in and out, and the graphic is very difficult to manage in Inkscape (slow reaction time when moving around the page and during difficult to group and ungroup vector layers). The map is not large (800 KB) and because it works perfectly other graphic programs (e.g., CorelDraw X5 and LibreOffice Draw 4.1) I believe it is a direct problem with Inkscape. Any suggestions?


title: Is there a QGIS bezier curve tool?

text: I'm a recent convert to QGIS but have been confused about the lack of support for drawing curved shapefiles, especially when Arc has almost too many. A lot of the work that our consultancy do is based around natural features such as lakes/rivers so there aren't many straight lines and being able to accurately represent or work with these features requires smooth curves as opposed to the tools i've been able to find on QGIS so far. Other threads i've found relating to this have stated that being able to draw curves would require a fundamental redesign of QGIS, how is this though? Other threads also recommend to use a rough shapfile first then simplify/smooth over to produce your curve but i would be worried about using this simply because you can't see the finished lay of the line when you draw it, which could be frustrating if your boundary moves after smoothing. I really love QGis and dont want to have to go back to Arc, but being able to draw curves is a big sticking point sadly. Any input/idea's would be greatly appreciated. Thanks, Albert


title: rasterlite build issues

text: i am able to build rasterlite library but when it is building .exe for tools(rasterlite_load) it gives linking error for external symbole _SPLite3_ like this rasterlite_load.obj : error LNK2019: unresolved external symbol _SPLite3_last_in sert_rowid referenced in function _export_tile there are many like this but all are due to splite I don't know what is missing. Searching google gives sqlite reference. Development Environment - windows and vc libraster version -1.1c Thanks


title: ArcGIS 10 -- Polygon Centroid Options and Appropriateness (FeatureToPoint, Labels, etc)

text: I've found a couple of other Stack Exchange questions related to this, but just like the ArcGIS online help, they were not satisfying. For instance, How does ArcGIS calculate the "inside" point for "Feature to Point" Doesn't tell me how the "Inside" point is determined. While the ArcObjects may refer to centroids and label points, I have at least 4 outcomes in my situation that vary differently, wildly in some cases. I have unidentifiable points inside zip codes, so I'm assigning each one to the centroid of that zip code. Then I will calculate straight line distances from these centroids to various locations of interest. Therefore, my choice of centroid is important. Some zip codes are also of odd shapes and may possibly be discontinuous. At least that is what I assume is generating some of the odd behaviors. Now, I have found four different ways to identify possible centrality of these polygons.   1. The polygon geometry contains a centroid information that can be directly extracted (and supposedly is used for labeling?).   2. The labeling location of the polygon (which I found does not always match where the above point is located--hence, the confusion).   3. Feature to point method with 'centroid' option.   4. Feature to point method with 'inside' option. Some of these or most of these may overlap. In some zip codes, they are all completely different. I would have expected 1. and 3. to be similar and maybe even 1. and 2. to be identical, depending on how labeling is set up (I'm using default). My question is not to ask "so how are these calculated?" That's neither here nor there since their algorithms are proprietary, and there's no point in trying to figure it out. The real question is practical: **how should one go about determining the appropriate centroid to use?** My conclusion is the feature to point method with the 'inside' option will be the most appropriate because it is guaranteed to reside within the polygon boundary. I also assume it will be relatively central. But the descriptions about this are nebulous. Hell, consider a crescent shaped polygon. The centroid would be outside the crescent but in the middle of the encapsulating envelope. Would an "inside" option literally just "push" it over into the crescent? We may not know the algorithm, but does the "inside" option actually identify centrality? Is it better than 1.? Maybe I'm missing something in my approach to this problem? Can the projection impact these comparisons? Is there a post-processing review of the polygons I should consider? Note, I'm completely open to using other tools (Python, R, QGIS, etc.). Your answers and advice are appreciated.


title: PostGIS query vs QGIS filter

text: I have a PostGIS dataset with around 3 million shapes. I only need to view about 1,500 of these shapes at any given time, but the shapes I need to see depends on which feature is selected in a different dataset. My question is: what is the speed difference between filtering my large dataset within PostGIS and then serving it up to QGIS, vs loading the entire dataset into QGIS and then using the baked-in layer filtering? I'm running PostGIS 2 and QGIS 1.8. I assume PostGIS is faster, but how much faster? I'd like a second opinion before I spend my time coding for one solution over the other. Thanks!


title: How to integrate CAD and GIS layers in ArcGIS?

text: I am supposed to integrate CAD files and different GIS layers in a project (as a job interview) in ArcGis. I am writing to know whether there will be any special consideration regarding to this procedure which should be taken into an account that effect the accuracy.


title: What happened to Live GPS Tracking in 1.8.0-Lisboa?

text: Dang, I was just getting used to Live GPS Tracking in 1.7.4-Wroclaw. I downloaded the new version 1.8.0-Lisboa on another computer. The Live GPS Tracking thing was in the View pull-down menu in 1.7.4. But, in 1.8.0-Lisboa, it's not there!!! Was Live GPS Tracking taken out of 1.8.0-Lisboa? or maybe Live GPS Tracking is hidden somewhere?


title: How to rotate the screen in ArcGIS-Online?

text: I need to print out my map in ArcGISOnline, however I'm using a grid system and need to rotate the screen so that I can print each map individually and have them line up. Any thoughts on how to do this? I will attach a screenshot to help better understand what I mean. ![enter image description here](http://i.stack.imgur.com/o0Hjb.png)


title: Can ArcPy 10.1 be used to generate true curve elliptical polygons in a file geodatabase?

text: As background, this question arose from trying to answer a question on **How to generate overlapping polygons from lines output from Table To Ellipse tool?** Using the ArcMap GUI it is very easy to digitize true curve polygons into a file geodatabase feature class using the **Ellipse** Construction Tool but ... Is it possible to write true curve elliptical polygons while reading rows (arcpy.da.SearchCursor) from a table containing a centre point, a major axis, a minor axis and an azimuth for each? I had hoped that there might be a token available with **arcpy.da.InsertCursor** to do this, but `SHAPE@` seems to be limited by what the **Geometry** object supports, and that does not appear to include true curves.


title: How do I edit a shapefile DBF to repair an extra row?

text: Something got screwed up and now the rows of my DBF beyond #1210 are offset by one, meaning that each shape has the attributes that should be assigned to the next shape. If I were in a spreadsheet program I would be able to select rows 1210+, copy, then select 1211+ and paste. But I can't do this in the QGIS attribute table editor, and I can't find an editor that will save shapefile DBFs without damage. * * * Apparently I have too little mojo to answer my own question this early. * * * Never mind - I fixed it in Excel. I had previously tried swapping columns in Excel and it didn't work, but I guess rows can be edited. The only minor issue is that the column titles are now all caps. But all the expressions that call them still work, so I guess that's not a problem.


title: gdalcompare.py not found in GDAL 1.10.1

text: I cannot find the script gdalcompare.py within the osGeo distribution of GDAL by http://www.osgeo.org/. It is the 1.10.1 though. Is there any way of downloading that script separately from any other source?


title: Solid Laptop Specs for running ArcGIS

text: I'm getting a Lenovo thinkpad w520 with 8gb of ram/2000m quadro GPU/2760qm Corei7 to use as my secondary work laptop. I mostly deal with large amounts of ortho photos and some GPS software. Do you guys think this is overkill or any of these specs useless? I'm just hoping this was a good purchase, I didn't think a mid consumer laptop would be able to run GIS very well (without crashing)....given that I do not do much 3D work, just alot of different GIS tasks, is the Quadro 2000m going to benefit me more than say a Nvidia 540m mid level card? Thanks for your help!!!


title: Alternatives to traditonal map server tile caching

text: Are there any alternatives to map server tile caching? It seems like there have been many improvements to computers since tile caching has become main stream. With SSD storage and large amounts of memory available on servers, is anyone forgoing all hassle of pre-rending tiles (imagery in particular) and does anyone have any benchmarks comparing dynamic map rendering vs cached?


title: Reprojection problem with gdalwarp

text: I have a GeoTiff file with WGS84 datum projection and I'm trying to reproject it using gdalwarp. Using this command               gdalwarp -t_srs '+proj=tpers +azi=19 +tilt=30 +lon_0=-109.5 +lat_0=23 +h=1274199.4 +wktext' atnight.tif atnighttilted.tif      returns this               ERROR 1: Too many points (10000 out of 10000) failed to transform,     unable to compute output bounds.      I'm using GDAL version 1.10.0. This is the contents of the .tfw file:               0.0066667901     0.0000000000     0.0000000000     -0.0066669136     -32.8672753199     72.8560317049      I've just started using gdalwarp. Thanks in advance, Tom The output of `atnight.tif` using `gdalinfo` is:               Driver: GTiff/GeoTIFF     Files: atnight.tif            atnight.tfw     Size is 9123, 4957     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (-32.870608714976186,72.859365161672642)     Pixel Size = (0.006666790125743,-0.006666913589395)     Metadata:       AREA_OR_POINT=Area       TIFFTAG_DATETIME=2012:11:26 15:13:57       TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)       TIFFTAG_SOFTWARE=Adobe Photoshop CS6 (Macintosh)       TIFFTAG_XRESOLUTION=72       TIFFTAG_YRESOLUTION=72     Image Structure Metadata:       COMPRESSION=LZW       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  ( -32.8706087,  72.8593652) ( 32d52'14.19"W, 72d51'33.71"N)     Lower Left  ( -32.8706087,  39.8114745) ( 32d52'14.19"W, 39d48'41.31"N)     Upper Right (  27.9505176,  72.8593652) ( 27d57' 1.86"E, 72d51'33.71"N)     Lower Right (  27.9505176,  39.8114745) ( 27d57' 1.86"E, 39d48'41.31"N)     Center      (  -2.4600456,  56.3354198) (  2d27'36.16"W, 56d20' 7.51"N)     Band 1 Block=9123x115 Type=Byte, ColorInterp=Red     Band 2 Block=9123x115 Type=Byte, ColorInterp=Green     Band 3 Block=9123x115 Type=Byte, ColorInterp=Blue


title: Spatial Join Get ZERO Value in ArcMap

text: The dataset contains about 70,000 point features. And I want to use spatial join tool in arctool box, to join a polygon feature. But it went that the values in the joined column are all zero, values in the 'join_count' are also zero . And the result was right when the point dataset is not that large(actually 40,000 is OK). What's the problem? Do I have any solutions? p.s. I don't want to use the join function in the right-click join. Although it went right, but I want to do spatial join using Geoprocessing.


title: How to export multiple Shapefiles to one KML file in QGIS?

text: As the title says, is it possible to export multiple Shapefiles to one single KML in Quantum GIS (Similar to "Map to KML" in ArcMAP)?


title: How can I reclassify a raster to floating point values?

text: I have a raster of landuse values that I need to reclassify into pollution concentration values. I have the pollution values in a separate table. They are decimal values (ex. 0.0056), and my landuse values are integers 1 through 6. I have a table mapping the landuse values to the pollution values. I tried using the Reclass By Table tool, setting my pollution values as the "Output value field", but the tool description says that this field must be an integer field. **How can I reclassify the values into decimal values?** I feel like I must be missing a much easier, alternate solution. My goal is to have a pollution concentration raster that I can use later in further calculations.


title: How to query nearest date in ArcGIS Server REST?

text: How to query nearest date (today) in ArcGIS Server REST? Table:               ID    NAME            DATA     1     test           2013-06-01 13:00:00     2     test1          2013-06-01 16:00:00     3     test2          2013-07-01 18:00:00      Expected result : query ID = 3


title: GIS Layer Versioning In GeoServer

text: i have developed the application for editing the map features online but i want the features to be saved as version of the layer not in original layer.this work is done using GeoExt +GeoServer+OpenLayers.


title: Is it possible to use commands from ArcMap's main menu in Model Builder

text: Is it possible to use the functions from ArcMap's main menu (for example Add Data, Start Editing) in my model? I can not find a way to add any other functions except Toolboxes.


title: Why is there no Map.SnapToLevels property in the ESRI JS API?

text: The ESRI Silverlight and WPF API's both have a **Map.SnapToLevels** property that can be set on or off. It "gets or sets whether zooming will snap to scale levels if the map contains cached layers". But with the JS API, you always snap to levels and there appears to be no way to turn that behavior off. Is there a way to disable the snap to levels behavior in the JS API? Or can someone provide the rationale behind ESRI's decision to not include that property in the JS API? (I'd like to be able to explain _why_ I cannot provide the same type of zooming in JS API as with WPF API.)


title: Checking "connectedness'' of line shapefile in ArcMap

text: I inherited a line shapefile in ArcMap that is used to represent a road network. The problem is that the road network is pretty large and it is impossible for me to: a. Tell if the end of each line (road segments) is "connected" to either an edge, vertex, or the end of another line (road segment); and b. Tell if the entire road network is "connected" - i.e., can I trace a route from any point on the road network to any other point on the road network? Is there a way to accomplish this in ArcMap? That is, is there a function that will allow me to highlight line segments that do no connect to other line segments, or something similar that would allow me to answer the two questions above?


title: How to connect lines in a network?

text: I have a shapefile of lines and another of polygons. When two or more lines intersect a polygon they should be connected (snapped), but sometimes during manual digitization this doesn't occur. I need to find and correct all these errors. I used `v.clean` from GRASS GIS and it automatically snapped a lot of vertices, which was good, but some of those disconnected lines are further from the tolerance I gave and if I increase the tolerance too much it begins to generate other types of errors. ![enter image description here](http://i.stack.imgur.com/cUNFp.png) Any ideas?


title: What arguments does MapBasic's ColumnInfo() take?

text: MapBasic's `tableInfo()` function has a KnowledgeBase docs page here listing the arguments and all the attributes it takes, and a usage example in the reference document. `columnInfo()` has neither of these (that I can find). I can't find any docs for it anywhere, other than a passing mention in the reference docs confirming that it exists: > The functions **TableInfo( )** , **ColumnInfo( )** and **NumTables( )** let > you determine information... > >   * **ColumnInfo( )** returns information about a column in a table, such as > the column’s name, the column’s data type, and whether the column is > indexed. > That's it. What arguments does it take (e.g. is it `columnInfo( nTabNum, nColNum, Attr )`? Or maybe `columnInfo( objColumn, Attr )`?) and what are the attribute names (e.g. are they named COLUMN_NAME, COLUMN_INFO_NAME, TAB_INFO_COL_NAME... etc?)


title: How can I use a sub-query in where statement of a query in ArcGIS 10.2 Map service?

text: I'm migrating my **ArcGIS** map services from version 10.0 to 10.2. It uses ArcSDE to manage spatial data in a **SQL Server** database. I was able to create the map service and everything seems OK. However if I add a sub-query in the WHERE statement of a query form in any map service, I receive "Unable to complete operation." error as the following image: ![enter image description here](http://i.stack.imgur.com/qb6vx.png) The same query was working in ArcGIS server 10.0. Also I can run the same query to select data, when I use the map document directly from ArcMap 10.2. The ArcGIS Server 10.2, ArcSDE, and SQL Server are installed on a server which is internally accessible from the web server. The ArcGIS Web Adaptor is installed on the web server. Both servers have windows Server 2008 and IIS. I appreciate your advice. Thank you


title: ArcObjects - cheapest way to check if license has been checked out

text: So I noticed a big performance hit when I add the following line to my code:               ESRI.ArcGIS.RuntimeManager.BindLicense(ESRI.ArcGIS.ProductCode.EngineOrDesktop);      I have two separate classes that call the BindLicense method. After the first time it is called, it isn't necessary to call again, but neither class can assume that has happened. What I want to know is if there is a cheaper way of checking if a license exists before attempting to bind? I feel like I've seen code to do this in Python, but don't know if there is the equivalent in ArcObjects...


title: QGIS API - Print Composer Template produces empty map?

text: I have a **print composer template** which I'm trying to populate with different subset queries from a CSV-derived vector layer using the **Python** API. However, all the output PDFs contain everything **except the layers on the map**! I'm doing this with a **standalone app** , and it functions fine in that it does produce PDFs and the data is all loaded and interrogated just fine (I can subset query, count features etc.) Here's the method I'm using to render the composition, as you can see I've tried to **set the layers** on the map in different places (actually there are 2 maps - one shows the extent rectangle of the other, but I can't get either to show anything) - this is my first QGIS scripting experience so I'm sure I'm missing something very obvious...               def compose(layers, species_binomial):         # vlayer is the CSV vector layer, outline is the ADM0 country vector.         vlayer, outline = layers              # construct a QStringList for the renderer setLayerSet         layers = QStringList(vlayer.id())         layers.append(outline.id())              # 1 maybe set the registry layers???         QgsMapLayerRegistry.instance().addMapLayer(vlayer)         QgsMapLayerRegistry.instance().addMapLayer(outline)              # 2 maybe set the canvas layer???         canvas = QgsMapCanvas()         canvas.setExtent(outline.extent())         canvas.setLayerSet([QgsMapCanvasLayer(outline), QgsMapCanvasLayer(vlayer)])              # 3 set the renderer layers???         renderer = canvas.mapRenderer()         renderer.setLayerSet(layers)         renderer.setExtent(outline.extent())              composition = QgsComposition(renderer)         with open(os.path.join(ROOT, 'species_composer.qpt'), 'rt') as fp:             composer_template = fp.read()         doc = QDomDocument()         doc.setContent(composer_template)         composition.loadFromTemplate(doc)              # 4 maybe set the map layers of one of the maps???         map1 = composition.composerMapItems()[1]         map1.setLayerSet(layers)              title = composition.getComposerItemById('title')         title.setText(" ".join(species_binomial))  # this works fiiiine              out_path = os.path.join(ROOT, 'out',  "%s_%s.pdf" % species_binomial)              # make a lovely blank PDF with a title and a blue map rect         ...      What am I doing wrong? Here's an example output... ![PDF output from print composer template](http://i.stack.imgur.com/5dSHY.png)