title: QGIS layer in .svg or with a better resolution?

text: I have a georeferenced layer, built from a nice map drawn with **Inkscape**. But when I convert it to *.png, *.jpg or *.tif, there is a pixelisation of the map. After georeferencing, it gives a pixelised map. I tried different definitions, sizes, etc but it always has much pixelisation. How could I import the same map, but with nice contours, like in a *.svg file? I have thought about converting it into a shapefile, but I haven't found out if it is possible?


title: How to clip by inidividual polygon shapes rather than by an aggregated shape consisting of all polygons in the shapefile.

text: I have a shapefile with overlapping polygons - network buffers around people's homes. I am trying to clip features such as tree canopy polygons and street polylines so that I can join the set of clipped features contained within each buffer to each network buffer polygon. Unfortunately I cannot get the Clip tool to clip the input feature polygons at the edge of the clip feature polygons if the latter overlaps with another clip feature polygon. How can I get the Clip tool to treat each of my clip features separately? Or is there another approach? I'm working in ArcGIS 10.


title: How to use QGIS Atlas with multiple layers/files per coverage?

text: I've got 30 contour datasets that all sit over the same extent. Each contour set is in its own file. I also have 30 extent polygons for each contour set as a separate file. Can I use Atlas to generate a series of maps - each displaying a 'pair' of datasets - ie. contours and extent for each pair over the same extent? Everything I've read indicates that you can loop through features in a layer with Atlas, but not through multiple layers. Could I work around this by adding a 'layer display list' as an attribute of the coverage layer? Maybe with a python script?


title: Export Map to GeoPDF

text: I am using GeoServer/OpenLayers to build a map website. One of things I want to do is to export a map to GeoPDF. It's not an export from desktop software. More specifically, map tiles are added as images and GML features are added as vector data. The only thing I have found that could allow me to do that is TerraGo SDK. Has anybody tried to do the same thing? Could anybody shed some lights? Thanks,


title: Labels appear and disappear for no apparent reason using QGIS?

text: I'm going through the Alaska sample data tutorial for QGIS 2.2.0. Some of the lake labels appear and disappear at different scales when zooming. At 1:314,460 Becharof is visible, but Lower Ugashik doesn't appear until 1:157,230. It seems like others have had similar problems. I have played with the rendering and placement dialogs to no avail. Even if I don't have the "magic combination" of parameters set, shouldn't they all behave the same? For my needs, this may well render (no pun intended) the program useless. I need it to print labels when I need them!


title: How to use Javascript Esri map as GPS Tracking?

text: I have an android application that get the gps coordinates of the device and stores it in a remote mysql database. I want to use ajax to display the coordinates in a javascript api using esri. I just need an example on how to refresh the function to keep calling for gps and plot in real time where that device is at. Any help is appreciated any attempts at such project would be great.I am using ArcGIS API for Javascript thats what I am using.


title: Looking for material for learning geospatial development in Java

text: I want to begin using Java for geospatial development. For the purpose of learning I am interested in suggestions on some books, web sites, etc. (i.e. reference materials) that would help me out. **Note:** I am using JDk and am also looking for suggestions on what geospatial libraries are available in Java.


title: Can I use a QGIS geometry in a Python PostGIS query?

text: I have a "Python Init function" for the attribute form in QGIS, which is called when I add a new geometry. I can access the geometry of the newly created feature like this               geom = feature.geometry()      Now, can I use this geometry in a PostGIS query, like for ST_Intersects? Regards, tkr


title: Plugins are not installed in qgis 1.8.0 Lisboa

text: I am trying to install plugins in qgis1.8.0 Lisboa but error message is appear as follows, > Failed to unzip the plugin package. Probably it's broken or missing from the > repository. You may also want to make sure that you have write permission to > the plugin directory: /home/hebert/.qgis//python/plugins Please tell me solution for it.


title: Filtering a Feature Class in using ArcGISDynamicMapServiceLayer

text: I'm using the following code in Flex to display all the features in a Feature Class defined in the "Data Layer".               var dataLayer:ArcGISDynamicMapServiceLayer= map.getLayer("Data Layers") as ArcGISDynamicMapServiceLayer;     var dataLayerId:int;          dataLayerId = 5;          dataLayer.visibleLayers.removeAll();     dataLayer.visibleLayers.addItem(dataLayerId);     dataLayer.refresh();      But I'd like to apply a filter like "WHERE MGD = 1" (MGD is an attribute defined in the feature class) to display only a subset of that feature class. Any help or autosuggestion would be much appreciated.


title: Default settings in Print Composer?

text: Is it possible to change default settings for things like font for labels in Print Composer? -TVZ


title: One to Many / Many to Many

text: I have a sanitary sewer line feature class where each line has a unique identifier. I have a table that lists the date a line was cleaned. It may have been cleaned multiple times so I can have multiple dates for each unique identifier. I also have another table that lists the date a line was inspected. It may have been inspected multiple times so I also have multiple dates for each unique identifier. I need to select each cleaned date and find an inspection date that lies within a 30 day range from when it was cleaned. For Example, From _Table_Clean_ :               Line A is cleaned 1/1/2014     Line A is cleaned 4/1/2014     Line A is cleaned 8/1/2014      From _Table_Inspected_ :               Line A is inspected 1/2/2014     Line A is inspected 4/2/2014     Line A is inspected 11/2/2014      I need to list that the first two were cleaned and inspected within 30 days. The third was cleaned but not inspected within 30 days. Is this many to many relationship possible? And does anyone have any suggestions on how to symbolize this on a map? Obviously the same exact line feature may be represented twice being cleaned and inspected on different dates.


title: Updating 30 Shapefiles, adding attribute and populating based on data

text: I need to add an attribute to each shapefile attribute table in ArcMap, I was thinking to do it in SQL, but I'm not sure how to go about doing it since I'm very new to ArcMap. Can someone explain to me how I would go about doing this? The added attribute would need each value populated based on another attribute in the table. Thanks.


title: GIS software with Bayesian Network capabilities

text: I am looking for any GIS software that has Bayesian Network capabilities. Ideally I'd be able to define grid layers as discrete or continuous inputs in to the network and then be able to produce layers based on the predictions of given nodes. I note there are several academic papers that seem to use a variety of different GIS and Bayesian tools but there doesn't seem to be any mention of an integrated solution. I'm starting to think that a solution doesn't exist and I'm going to hack my own solution together.


title: function ST_Line_Interpolate_Point with multilinestring

text: Is there a function like ST_Line_Interpolate_Point but using a multilinestring?


title: Error building pyramids/statistics on raster with arcpy

text: I am building a script to process rasters with the option of adding the extra step of building pyramids and statistics at the end. For some reason, it will not recognize inputs as legitimate raster datasets, despite them all being .tifs. The script has no problem processing the rasters prior to this point using the same conventions. I've tried this code, the "Build Pyramids and Statistics" tool:               outWorkspace = arcpy.GetParameterAsText(2)     #...lots of stuff happens here     #Build pyramids?     try:         if pyrs == 'true':             arcpy.AddMessage("Building pyramids...")             arcpy.BuildPyramidsAndStatistics_management(outWorkspace)         else:             arcpy.AddMessage("No pyramids will be built.")     except:         import traceback         arcpy.AddError("An error occurred building pyramids.  Here's the ArcGIS message and a traceback:")         arcpy.AddError(arcpy.GetMessages(2))         arcpy.AddMessage(traceback.print_exc())      This returns an error that is undefined. I've also tried this code, the "Build Pyramids" and "Calculate Statistics" tool:                   #Build pyramids?     try:         if pyrs == 'true':             arcpy.env.workspace = outWorkspace             arcpy.AddMessage("Now building pyramids on output images...")             rasters = arcpy.ListRasters("", "TIF")             for raster in rasters:                 arcpy.BuildPyramids_management(raster)                 arcpy.CalculateStatistics_management(raster)         else:             arcpy.AddMessage("No pyramids will be built.")     except:         import traceback         arcpy.AddError("An error occurred building pyramids.  Here's the ArcGIS message and a traceback:")         arcpy.AddError(arcpy.GetMessages(2))         arcpy.AddMessage(traceback.print_exc())      Thanks in advance for the help, as always.


title: How to change class break schemes in QGIS

text: When using Graduated Symbology in QGIS, by default the bounds seem to be overlapping. For example, class 1 shows 0-100.000; class 2 shows 100.000-200.000; class 3 shows 200.000-300.000 where it appears as though the value 100 is captured in class 1 & 2; 200 is captured in class 2 & 3. Is there a way to change the default class breaks so they are not overlapping (e.g. Class 1 = 0-100; class 2 = 101-200; class 3 = 201-300). I know in ArcGIS you can adjust the number of decimal points in the label - is this also a capability of QGIS? Thanks,


title: Scripting a Loop of Clip (Analysis)

text: I wish to loop function that clips a shapefile by a feature class. I want to make the clip the same as the input name. The names of the files are WVHGT_201403281200_f180_201404050000.shp and I just want f180 part as the output feature class. The clip feature will remain consistent throughout. I have the syntax down from here. I am just not sure how the write the for loop to do this.


title: Is the SAGA tool Basic terrain analysis included in Processing?

text: I can't find the saga tool Basic terrain analysis in the Sextante/processing toolbox. Is this not included or is there a chance to install it?


title: How does Arc projecting data without spatial reference?

text: I have a raster file downloaded from here. When I load it into ArcMap, it says the data is missing spatial reference, and that it can be drawn in ArcMap but not projected. Nevertheless, using the Define Projection tool, I can assign a coordinate system that seems to work pretty well. How does Arc do this if there was no spatial reference to begin with? Is this a valid approach to data missing spatial reference?


title: Openlayer not support touchstart for map object?

text: For a while, our website only supports desktop environment so we didn't set up any touch handlers. Openlayer working fine so far. But today we got a request to add support for touch devices too. During our tests, OpenLayers only triggers touchmove and touchend for map, but no touchstart. Openlayer đi trigger touchstart for overlay object though. So does OpenLayers support touchstart for map? If not, is there any event like touchstart which OpenLayers would trigger? We need touch start to determine the type of touch: Drag / Click / Double click. We're using jQuery to catch touchstart at the moment to overcome this problem ran into few problem.


title: Third party WMS getcapability URL

text: i need some third party WMS Getcapabilities URLs. From where i can get those URLs. I want to test those URLs in my Geoserver.


title: How to convert a set of points to a solid object?

text: I have a data generated as the output of some numerical processes. I need them to be converted into 3D solid objects (as CAD solid format) to make them importable for some particular numerical analysis software such as FEM. One initial stage is to do some Boolean operations on them. If they were is solid format it is OK, but the data is only coordinates and that's all! The data set consist of points in which a set of points makes a particular shape say a cube. Each object has own associated points. If we draw lines between points of each object a proper object appears. These data need to be exported in 3D format consisting of only solid objects. Then a series of Boolean operations will apply. So the questions are: **How to convert such a data set to for example a _dxf_ format consisting only _solid_ objects?** In particular purpose, **what is the structure** (format or so) **of a solid object in CAD** so we may write some simple code to convert our data set to solids? Any other idea is more than welcome. _(Note that the data set is large so we need finally to do some programming in any language required.)_ Still waiting for your responses.


title: Allowing cross-domain access of third-party WMS server from OpenLayers

text: I have a web site with an OpenLayers client. I also have a GeoServer running on the same domain. My OpenLayers uses services like WMS, GML, etc. from this GeoServer installation and everything runs fine. Now, I want to add another WMS layer to my existing set of layers. This layer is provided by a third-party WMS (not owned by me). How can I add this layer to my OpenLayers by bypassing the cross-domain issue. PS: I am aware that GeoServer allows adding third-party WMS as a new data store and lets me publish the available layers, but I am facing some projection issues there which I will be asking in a different question.


title: Using Tilestache and Leaflet To Display Rasters In Native Projection

text: I have a collection of old Idaho state highway maps in GeoTiff format. The projection that these maps have been georeferenced in is a transverse mercator projection, called IDTM83 (Idaho Transverse Mercator - NAD 83). The projection is defined as: > +proj=tmerc +lat_0=42 +lon_0=-114 +k=0.9996 +x_0=2500000 +y_0=1200000 > +ellps=GRS80 +units=m +no_defs The projection is basically a modified UTM projection that compensates for that fact that the state of Idaho is split by UTM zones 11 and 12. My goal is to serve the maps in the native IDTM83 projection using Tilestache and Leaflet. I am able to serve these same rasters that have already been projected into world mercator, but the the projection distorts much of the line work from the original map. There have been some helpful resources on the internet on how to serve projected data using Tilestache and Leaflet, here:   1. https://gist.github.com/skogsmaskin/6452343   2. http://blog.thematicmapping.org/2012/07/using-custom-projections-with-tilecache.html.  The mapnik configuration for one of my raster maps looks like:               <Map background-color="steelblue"      srs="+proj=tmerc +lat_0=42 +lon_0=-114 +k=0.9996 +x_0=2500000 +y_0=1200000 +ellps=GRS80 +units=m +no_defs"     maximum-extent="2205517,1134032,2834671,2040323">     <Style name="raster-style">       <Rule>           <RasterSymbolizer/>       </Rule>     </Style>          <Layer name="idaho-1947" status="on" srs="+proj=tmerc +lat_0=42 +lon_0=-114 +k=0.9996 +x_0=2500000 +y_0=1200000 +ellps=GRS80 +units=m +no_defs">     <StyleName>raster-style</StyleName>     <Datasource>         <Parameter name="type">gdal</Parameter>         <Parameter name="file">GIS/IdahoMaps/1947_geo_idtm.tif</Parameter>         <Parameter name="format">tiff</Parameter>     </Datasource>     </Layer>     </Map>      Using Python, I am able to generate maps using bounding boxes from the IDTM coordinate system with the Mapnik API. So I know that Mapnik is not the problem. My Tilestache config is pretty basic:                   {       "cache": {"name": "Test","verbose":"True"},       "layers": {         "idaho-1947": {             "provider": {"name": "mapnik", "mapfile": "./mapnik.xml"},             "projection": "TileStache.Goodies.Proj4Projection:Proj4Projection('+proj=tmerc +lat_0=42 +lon_0=-114 +k=0.9996 +x_0=2500000 +y_0=1200000 +ellps=GRS80 +units=m +no_defs',[5545984, 2772992, 1386496, 693248, 346624, 173312, 86656, 43328, 21664, 10832, 5416, 2708, 1354, 677, 338.5, 169.25, 84.625, 42.3125, 21.15625, 10.578125, 5.2890625, 1])"         }}      My leaflet viewer is very basic, set up like on Link #2 above:                 var crs = L.CRS.proj4js('SR-ORG:6677', '+proj=tmerc +lat_0=42 +lon_0=-114 +k=0.9996 +x_0=2500000 +y_0=1200000 +ellps=GRS80 +units=m +no_defs', new L.Transformation(1, -432000, -1, 2040323));       crs.scale = function(zoom) {           return 1 / (234.375 / Math.pow(2, zoom));       };            var map = new L.Map('map', {         crs: crs,         layers: [           new L.TileLayer('http://localhost:8080/idaho-1947/{z}/{x}/{y}.png', {             minZoom: 0,             maxZoom: 4,             continuousWorld: true           })         ],         center: new L.LatLng(45, -116),         zoom: 3,         continuousWorld: true       });      I intend to refine this once I can get the maps to display, but the example is using UTM data and the scales should be pretty similar. I am stuck here. Using the JavaScript examples provided in both the links above, I cannot display any of the maps - only the lightblue background. I know what I am trying to do is possible in ArcGIS Server, but I would really like to leverage some of the OpenSource tools.


title: How to create custom markers for QGIS?

text: Is it possible to define a custom marker as symbolizer for a point in QGIS? I'm looking into using a variation on the "Vector Field Marker" with a couple of arrows pointing at dynamic angles (the angles being read as point attributes). I imagine it's doable in C++ (within the QGIS code) but I'm looking for a scripting solution (python). I'm adding a representation of what I'm trying to achieve - for two points A and B with attributes x,y,z. ![custom markers](http://i.stack.imgur.com/Af3Al.png)


title: geodatabase and relationship

text: I have two tables in a personal geodatabase and one of them has 13 records and 6 fields, the other has 118 records and 125 fields. There is a relationship class between the two tables where the first table is the origin and the other (bigger) table is the destination. What I need is to take the records satisfying a condition from the second table and sum those records and put into a new field in the origin table


title: what happens to attribute table while performing union, intersection, difference of gml and shp file in QGIS

text: I am trying to do an analysis of the data conversion from shapefile to GML file using QGIS. For the analysis I have taken a point, line and polygon file. I first took a shapefile and converted it into a GML file; after converting it I compared the feature count of both files, which were same. I then performed geoprocessing analyses (difference, intersection, union, etc) on both of these files. I am doing a geoprocessing analysis because if the shapefile has been successfully converted into GML without errors, then a difference of both these files will give me a null set. Similarly for others: http://infogeoblog.wordpress.com/2013/01/08/geo- processing-in-qgis/#comments I found that for union, it is not showing the number of tuples in the resultant attribute table as sum of tuples in both GML file attribute table and shapefile attribute table. Similarly, for intersection and difference it is not showing the results as expected (by expected I mean, as we do in set theory). For some files it also shows null values. I want to know what happens to tuples in the attribute table while performing union (intersection, difference) of GML and shapefiles. Do they get added up or is some join performed?


title: PostGIS | Get degrees clockwise from north of a linestring's beginning and end

text: I have a table full of linestrings in PostGIS. I would like to know the degrees clockwise from north of the starting point of the linestring and also of the end point. I was thinking that maybe creating a circle with both points could help to get this two values (example: X=32º, Y=209º) but i have no clue how to achieve this.... ![enter image description here](http://i.stack.imgur.com/kzFrn.png) Regards


title: How do you assign classes to percentage data in a data layer?

text: As a learning project I need to coneptualize a GIS output for assigning classes to data in a raster layer. In this case i have the unemployment rates for relevant magisterial districts in my database. So I need to create a layer where I can assign classes to the unemployment rate (percentage of people unemployed) as low (0-10%), medium (10-30%), or high (40&+) within each district. I've thought that a raster layer containing the unemployment percentages as cell values could work, but how would I go about outputting this through a GIS? I'm sorry if this question is a bit non-specific, but any tips will be greatly appreciated! Thank you


title: How to check if REST plugin is installed for GeoServer

text: Hi I am having trouble uploading vector files onto GeoServer 2.4 using 'Import Vector' in the QGIS Processing toolbox. After reading around about the 'REST' plus in I went to: http:// Also when I run the tool I find an error code 404 in the log. GeoServer is running on a company server along witch apache 6.0 Could someone please tell me how to check the REST plug in, and if it is not installed how to get the REST plug in for GeoServer 2.4.0. Thank you, Sean0


title: How to Change the Default Displaying Order of Esri Basemaps in basemapGallery?

text: I just started using ArcGIS JS APIs. In my app, I want to combine a few ESRI basemaps with my own basemap service. I am using basemapGallery dijit to change basemaps. By default, basemapGallery shows the imagery map. I want to show OpenStreetMap as the default. How can I make this happen? The code in my app is simple and directly borrowed from ArcGIS samples. So I did not put it here. Thanks!


title: Don't move Layout

text: I'm doing a vb.net button to run in ARCGIS. the button is working, the problem is as I'm working in the layot when I press the button, change the layout coordinates and the map coordinates. Is a way that the layout view don't move and keep the Coordinates (0,0)               I leave the code:          MxDoc = My.ArcMap.Document          ActiveView = MxDoc.ActiveView          Dim x As Double = txtEast.Text          Dim y As Double = txtNorth.Text          Map.MapScale = 4000 'zoom to this scale          Envelope = ActiveView.Extent          ,Create and Add the point coordinates     Point = New ESRI.ArcGIS.Geometry.Point     Point.PutCoords(x, y)          'Envelope Position go to the middle of the point     Envelope.CenterAt(Point)               'Refresh View      ActiveView.Extent = Envelope      MxDoc.ActiveView.Refresh()Map = MxDoc.FocusMap      Thank You!


title: Openlayers map save feature SRID error into PostGis

text: OpenLayers.ProxyHost = "proxy.ashx?url=";                  //set up the modification tools             var DeleteFeature = OpenLayers.Class(OpenLayers.Control, {                 initialize: function(layer, options) {                 OpenLayers.Control.prototype.initialize.apply(this, [options]);                 this.layer = layer;                 this.handler = new OpenLayers.Handler.Feature(                     this, layer, {click: this.clickFeature}                 );                 },                 clickFeature: function(feature) {                 // if feature doesn't have a fid, destroy it                 if(feature.fid == undefined) {                     this.layer.destroyFeatures([feature]);                 } else {                     feature.state = OpenLayers.State.DELETE;                     this.layer.events.triggerEvent("afterfeaturemodified",                                                 {feature: feature});                     feature.renderIntent = "select";                     this.layer.drawFeature(feature);                 }                 },                 setMap: function(map) {                 this.handler.setMap(map);                 OpenLayers.Control.prototype.setMap.apply(this, arguments);                 },                 CLASS_NAME: "OpenLayers.Control.DeleteFeature"             });                  function showMsg(szMessage) {                 document.getElementById("nodelist").innerHTML = szMessage;                 setTimeout(                 "document.getElementById('nodelist').innerHTML = ''",2000);             }                  function showSuccessMsg(){             alert("Succes");                 showMsg("Transaction successfully completed");             };                  function showFailureMsg(){             alert("NoneSucces");                 showMsg("An error occured while operating the transaction");             };                  $(document).on("pageshow", "#mappage", function( event ) {                          haritacontrol=1;                      //set up a save strategy                 var saveStrategy = new OpenLayers.Strategy.Save();                     saveStrategy.events.register("success", '', showSuccessMsg);                     saveStrategy.events.register("failure", '', showFailureMsg);                        var bounds = new OpenLayers.Bounds(417138.6205791628, 4243634.852681775,486211.1950145211, 4331542.555877601                     );                     var options = {                         controls:[                         //allows the user pan ability                         new OpenLayers.Control.Navigation(),                         //displays the pan/zoom tools                                            new OpenLayers.Control.PanZoom(),                         //displays a layer switcher                         new OpenLayers.Control.LayerSwitcher(),                          //displays the mouse position's coordinates in a                          //<div> html element with id="coordinates"                         new OpenLayers.Control.MousePosition(                         {div:document.getElementById("coordinates")})                        ],                         maxExtent: bounds,                         maxResolution: 343.3894656086959,                         projection: "EPSG:2322",                         units: 'm'                     };                     map = new OpenLayers.Map('map', options);                          // setup tiled layer                     tiled = new OpenLayers.Layer.WMS(                         "Geoserver layers - Tiled", "http://192.168.1.236:8080/geoserver/Kaski/wms",                         {                             LAYERS: 'KaskiGrup',                                       tiled: true,                             tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                         },                         {                             buffer: 0,                             displayOutsideMaxExtent: true,                             isBaseLayer: true,                             yx : {'EPSG:2322' : true}                         }                      );                      //wfs-t editable overlay                 var wfs_layer = new OpenLayers.Layer.Vector("Editable Features", {                     strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],                          protocol: new OpenLayers.Protocol.WFS({                         version: "1.0.0",                         // loading data through localhost url path                         url: "http://192.168.1.236:8080/geoserver/wfs",                         featureNS :  "http://192.168.1.236:8080/KASKIDB",                         maxExtent: bounds,                         // layer name                         featureType: "DBadres_mahalle",                         // geometry column name                              geometryName: "geom",                         schema: "http://192.168.1.236:8080/geoserver/wfs/DescribeFeatureType?version=1.0.0&typename=KASKIDB:DBadres_mahalle"                     })                     });                      map.addLayers([tiled, wfs_layer]);                 map.zoomToExtent(bounds);                      // add the custom editing toolbar                 var panel = new OpenLayers.Control.Panel(                     {'displayClass': 'customEditingToolbar'}                     );                      var navigate = new OpenLayers.Control.Navigation({                     title: "Pan Map"                     });                      var draw = new OpenLayers.Control.DrawFeature(                     wfs_layer, OpenLayers.Handler.Polygon,                     {                         title: "Draw Feature",                         displayClass: "olControlDrawFeaturePolygon",                         multi: true                     }                     );                                var edit = new OpenLayers.Control.ModifyFeature(wfs_layer, {                     title: "Modify Feature",                     displayClass: "olControlModifyFeature"                     });                      var del = new DeleteFeature(wfs_layer, {title: "Delete Feature"});                      var save = new OpenLayers.Control.Button({                     title: "Save Changes",                     trigger: function() {                         if(edit.feature) {                         edit.selectControl.unselectAll();                         }                         saveStrategy.save();                              // alert('saved');                     },                     displayClass: "olControlSaveFeatures"                     });                      panel.addControls([navigate, save, del, edit, draw]);                 panel.defaultControl = navigate;                 map.addControl(panel);                               });      My code is like that, but when I save feature, it is added but there is no information on table. For example fid columns nothing. Also when I click saved feature, i got this error :               Internal error occurred java.io.IOException null ERROR: Operation on mixed SRID geometries      How can I give columns parameters. Also I found another example about insert but it gives success WFS transaction but nothing change. That code :                   var map;         OpenLayers.ProxyHost = "proxy.ashx?url=";                    $(document).on("pageshow", "#mappage", function( event ) {           var bounds = new OpenLayers.Bounds(417138.6205791628, 4243634.852681775,486211.1950145211, 4331542.555877601             );             var options = {                 controls: [ new OpenLayers.Control.Zoom(),                     new OpenLayers.Control.Navigation()],                 maxExtent: bounds,                 maxResolution: 343.3894656086959,                 projection: "EPSG:2322",                 units: 'm'             };             map = new OpenLayers.Map('map', options);                  // setup tiled layer             tiled = new OpenLayers.Layer.WMS(                 "Geoserver layers - Tiled", "http://192.168.1.236:8080/geoserver/Kaski/wms",                 {                     LAYERS: 'KaskiGrup',                               tiled: true,                     tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                 },                 {                     buffer: 0,                     displayOutsideMaxExtent: true,                     isBaseLayer: true,                     yx : {'EPSG:2322' : true}                 }              );         var roads = new OpenLayers.Layer.WFS(             "Roads",              "http://192.168.1.236:8080/geoserver/wfs",             {typename: 'KASKIDB:DBadres_mahalle'},             {                 typename: 'DBadres_mahalle',                  featureNS: 'http://192.168.1.236:8080/KASKIDB',                  geometryName: "geom",                 extractAttributes: false             }         );         roads.style = OpenLayers.Util.applyDefaults({strokeColor: "#ff0000"},                  OpenLayers.Feature.Vector.style["default"]);                        map.addLayers([tiled,roads]);              var panel = new OpenLayers.Control.Panel(             {displayClass: 'olControlEditingToolbar'}         );              var drawPol = new OpenLayers.Control.DrawFeature(             roads, OpenLayers.Handler.Polygon,             {displayClass: 'olControlDrawFeaturePath'}         );         drawPol.featureAdded = function(feature) {             feature.layer.eraseFeatures([feature]);             // cast to multilinestring             feature.geometry = new OpenLayers.Geometry.MultiPolygon(                 feature.geometry             );             feature.style.strokeColor = "#ff0000";             feature.state = OpenLayers.State.INSERT;             feature.layer.drawFeature(feature);         }                   var drawPoint = new OpenLayers.Control.DrawFeature(             roads, OpenLayers.Handler.Point,             {displayClass: 'olControlDrawFeaturePoint'}         );         drawPoint.featureAdded = function(feature) {             feature.layer.eraseFeatures([feature]);             // cast to multipoint             feature.geometry = new OpenLayers.Geometry.MultiPoint(                 feature.geometry             );             feature.style.strokeColor = "#0000ff";             feature.state = OpenLayers.State.INSERT;             feature.layer.drawFeature(feature);         }                   panel.addControls(             [new OpenLayers.Control.Navigation(), drawPol]         );              map.addControl(panel);         map.addControl(new OpenLayers.Control.LayerSwitcher());              map.zoomToExtent(bounds);                                                                                });


title: QGIS 2.0 WMTS Issue

text: When using a WMTS within QGIS 2.0 (using the "?SERVICE=WMTS&REQUEST=GetCapabilities" addition) only one tile per layer is loading. The WMTS works in other GIS software and this problem would appear to be unique to QGIS 2.0. Does anyone know why this would he happening and if there is a workaround to this issue?


title: (How) Can I define my own placement settings for labels?

text: Is it possible to define my own Placement settings for my labels? ![Placement properties](http://i.stack.imgur.com/kU8da.png) The predefined ones are not bad, but I would need a few more of my own design. How can I accomplish that?


title: Using ECQL filter on user form using GeoServer and OpenLayers

text:
I know this is more of a JS question (skills of which I am inadequately inept) but it is using OpenLayers and GeoServer. I am trying to apply a filter with a user form on a Layer, but am not sure why it is not working. If I manually write the code into the CQL_FILTER, it works, but it won't work with the UpdateFilter button.               function init_map() {              var layerExtent = new OpenLayers.Bounds( -9816656, 5296800, -9772016, 5345759);              map = new OpenLayers.Map('map',                  {'restrictedExtent': layerExtent});         basemap = new OpenLayers.Layer.OSM("Simple OSM Map");         map.addLayer(basemap);         markers = new OpenLayers.Layer.Markers("Markers");         map.addLayer(markers);         map.setCenter(             new OpenLayers.LonLat(-71.147, 42.472).transform(                 new OpenLayers.Projection("EPSG:4326"),                 map.getProjectionObject()             ), 12         );              var tiled = new OpenLayers.Layer.WMS(                     "Verizon Cell Towers", "http://mpdflash:7070/geoserver/Milwaukee/wms",                     {                         LAYERS: 'Milwaukee:Verizon Cell Towers',                         tiled: true,                         transparent: true                     },                     {                         displayOutsideMaxExtent: true,                         isBaseLayer: false,                         yx : {'EPSG:6948' : false}                     }                  );                map.addLayer(tiled);              tiled.mergeNewParams({'CQL_FILTER': "filter"})     }          function updateFilter(){             var filter = document.getElementById('filter').value;         }          function resetFilter() {             document.getElementById('filter').value = "";             updateFilter();         }         and
                 <input type="text" size="80" id="filter"/>


title: Where to get Annual Lightning Data (Zipcode) to use with ArcGIS for Desktop and AnyLogic?

text: I am building a Dynamic Agent Based Model for School and was looking for Annual Lighting Strike Data for zip codes in Ohio. I can get the state count from NOAA but want to get a bit smaller. I was wondering if anyone knew of a source? (I am using the data to analysis the risk of lighting in and around the Cincinnati Area.) I am using ArcGIS for Desktop to do the maps and **AnyLogic** for my Model.


title: MapInfo 10.5.1 Update and MS Access Database

text: I am having problems with a database that I have developed. The data uses MapInfo for any of the spatial parts of it. I installed the update to MapInfo 10.5.1 and I now get an error when I am trying to access an ADODB.recordset. object which previously worked. I get the error: Method 'Connection' of object '_CurrentProject' failed when the database tries to run the code:               rs5.ActiveConnection = CurrentProject.Connection       (where rs5 is a ADODB.Recordset): Has anyone else seen anything like this? Is this related to the updated drivers for Office 2010 that were installed with the update? I am using Access 2007. The database still works on machines that are 10.5.


title: Configuring & debugging a new mapnik / tirex installation

text: I've a stack OpenLayers > Apache mod_tile > tirex working with the tirex test backend on a Ubuntu 10.04 virtual machine. However, I now have troubles with the mapnik-example backend included with tirex. OpenLayers works fine, mod_tile does not complain about the newly added "AddTileConfig /tiles/example/ example" to my Apache httpd.conf. Also, the tirex-master receives a request from mod_tile. I suspect this problem to lie with the configuration of the tirex Mapnik rendering backend. Using tirex-status command I can see tirex-master currently rendering the tiles for the example map. The rendering process times out after about 300 seconds and no (meta)tiles are being written to the tile directory. I suspect a wrong configuration or problems with user/write access. Does anyone know where I can find debug information from tirex or mapnik that can help me solve this problem? Thanks in advance. * * * EDIT: I am looking for clues about what is going on in some logfiles. I suspect the problem to lie with the tirex-backend-mapnik process, which seems to be managed by the tirex-backend-manager process. In the [OSM Tirex Logging wiki page][1] it says that the tirex-master and tirex-backend-manager will log to the syslog by default, which should be located under /var/log/messages. This logfile has the following entry at the end ` Dec 23 12:50:13 ubuntumapserver kernel: [10967.139797] mapnik[2079]: segfault at 0 ip 000000000041314b sp 00007fffacfa9530 error 4 in mapnik[400000+1f000]   Dec 23 12:50:14 ubuntumapserver kernel: [10968.391401] mapnik[2077]: segfault at 0 ip 000000000041314b sp 00007fffe23608a0 error 4 in mapnik[400000+1f000] ` Does anyone know have any clues about what can be causing this error? * * * EDIT 2: Hi Dane, Thanks for your reply. I installed Mapnik from source using the [Mapnik wiki tutorial][2] (Installation on a VM guest running Ubuntu 10.04 Server x64) ` python scons/scons.py configure DEBUG=True INPUT_PLUGINS=all OPTIMIZATION=3 SYSTEM_FONTS=/usr/share/fonts/ python scons/scons.py sudo python scons/scons.py install ` I have successfully used Mapnik standalone via the Python commmand line to render images. Furthermore I did an strace during a request from OpenLayers. (I saw you recommending this elsewhere). I can see mod_tile.so being called. However, there is no data about tirex or mapnik. [strace][3]   [config.log from Mapnik installation][4]   Therefore, I guess you might be correct about tirex **not working with a compiled version of Mapnik**. I will try to find out if this is indeed the problem. 1_http://wiki.openstreetmap.org/wiki/Tirex/Logging (About Tirex Logging)   2_http://trac.mapnik.org/wiki/UbuntuInstallationOld (Mapnik installation)   3_http://pastebin.com/L8VBzvgM (strace of a tile request)   4_http://pastebin.com/yvLQDr39 (Mapnik scons config.log)


title: How to create regular point grid inside a polygon in Postgis?

text: How to create, inside a polygon, a regular grid of point spaced x,y in postgis? Like the example: ![alt text](http://i.stack.imgur.com/ynKBu.png)


title: Transform wkb to gml geometry using FME

text: Maybe (and probably) this is stupid question, but I am totaly new to FME and geometry conversion. I need to transform WKB field from SQLServer to GML geometry using FME software. How can I do this?


title: Count points in a polygon that changes form ("on the fly calculation")

text: I have tried to find a plugin for this, with no luck. Perhaps someone knows of a plugin that can do the following: In this project I create a lot of polygons that contain points, and the size of the polygons is dependent of how many points there are inside the polygon. Therefore I need to manually change the size of the polygons to get a certain number of points inside it. It is tedious work to run a new "Points in polygon" command as soon as I have changed a polygon and want to know the number of points. Is there any way to get this calculations directly on screen some how (a sort of "on the fly" point calculation)?


title: Tool in ArcGIS for stream flow analysis?

text: I'm just wondering if there is any tool in ArcGIS that can be used to analyse stream flow measurements. For example, estimating runoff, baseflow and baseflow index based on stream flow measurements.


title: How to store esri geometry types in sql server

text: I have a txt file containing some data of polylines. Actually it is saved json objects. The features in it are esriPolyLine objects. Now i want to store it in my sql server db as geography data types through EF. Is there anywhere some great examples for doing this. How should the string value loop like to be able to enter the values using DbGeography.From(polylinestring, 4326) Or just a simple example of storing esriGeographyPolyline objects in sqlserver


title: Using Arcgis geometry in OpenLayers

text: I am sending request to Arcgis Rest service and getting result as JSON object. ![enter image description here](http://i.stack.imgur.com/rbomz.png) I want to zoom to this Ring array.               var extent = new OpenLayers.Bounds(extent.geometry.rings[0]);      I tried but did not work. I am new at OpenLayers.


title: Open source alternative to ArcGIS geoprocessing service concept?

text: I am searching for any examples of implementing a logic that Esri has built with geoprocessing services. You author a custom tool or it can be a system tool >> you publish it exposing this functionality as a web service >> you consume the web service with a client. A classical simple example: user can click on the map to create a point > the coordinates are sent to the server > buffering operation is being executed at the back-end > the buffer zone feature (or just an image) is sent back to client and shown on the map canvas. I've been developing Esri GP services for last 3 years and it has been a nice experience specifically since version 10.1. **Are there any examples of the open source system/solution that one could use to build a similar logic?** I've googled and found sextante gis but there was nothing specific on how this works and how much of logic that is available through Esri is available. 52North also seem to be working on that, but again I am missing the clean explanation of the concept.


title: Concatenate Fields in ArcMap 10

text: I added a field to an address layer in ArcMap 10 to concatenate separate address fields. I used the field calculator and a VB Script to accomplish this. So, for example, if the separate fields are "10" "Main" and "Street", then the new field would be "10 Main Street". The problem is that the concatenated field needs to update dynamically. So if the user changes "Main" in the database to "East Main" then that field will update in ArcMap. However, the concatenated field will not update to "10 East Main Street" unless I rerun the Field Calculator script. I need a way for this new field to always reflect the current data, but still do the concatenation. Does anyone have any suggestions? Thanks, Seth


title: Centralizing shapefile data into database

text: I've got hundreds of shapefiles from various different GIS projects that I want to start consolidating into a single database platform, currently attempting this with Postgres / PostGIS. Hardly any of the data is standardized - meaning that it's a lot of the same data _types_ , but the particular attribute names / types don't match. Where should I start tackling this? Should I develop a standard model to migrate each shapefile into first (e.g. Hydro_line, transport_line, Hydro_poly standards, etc)? An alternative is to just import each shapefile into Postgres individually, so each shp becomes a table in the database, but I'm not sure about this in terms of performance and organization. Feels sort of like delaying the inevitable... Any advice on dealing with this daunting task?


title: Autodesk MapGuide Enterprise MrSid on Linux

text: I am running MapGuide Enterprise 2011 on CentOS 5 linux. Using the Autodesk Raster FDO, i loaded a mrsid file (through the "Load file-based data" method). Whenever i try to view the layer with the mrsid file, the server crashes. most of the time it says "Segmentation Fault", but sometimes it just dies without saying anything. Here's what the error log says when i try to view the layer:               <2011-01-14T10:39:25>   53263248        Autodesk MapGuide Studio v2.2.0.5305       198.245.206.237 Administrator     Error: An exception occurred in FDO component.         Only MOSAIC(), CLIP() and RESAMPLE() are supported.     StackTrace:     - MgFeatureServiceHandler.ProcessOperation() line 83 file FeatureServiceHandler.cpp     - MgOpSelectFeaturesSpatial.Execute() line 110 file OpSelectFeaturesSpatial.cpp     - MgServerSelectFeatures.SelectFeatures() line 266 file ServerSelectFeatures.cpp      I've tried a few other mrsid files with the same problem. TIFF files load fine. Anyone has any insight on this? Thank you!


title: Problem when converting sosi to shapefile with ogr2ogr on Linux

text: I am trying to convert a SOSI file (Norwegian geodata format) to a shapefile. I have set up gdal, but apparantly there is a lot of things that I dont have support for. When I do a ./configure, this is the result:                 LIBTOOL support:           yes            LIBZ support:              external       LIBLZMA support:           no       GRASS support:             no       CFITSIO support:           no       PCRaster support:          internal       LIBPNG support:            internal       DDS support:               no       GTA support:               no       LIBTIFF support:           internal (BigTIFF=yes)       LIBGEOTIFF support:        internal       LIBJPEG support:           internal       12 bit JPEG:               yes       12 bit JPEG-in-TIFF:       yes       LIBGIF support:            internal       OGDI support:              no       HDF4 support:              no       HDF5 support:              no       NetCDF support:            no       Kakadu support:            no       JasPer support:            no       OpenJPEG support:          no       ECW support:               no       MrSID support:             no       MrSID/MG4 Lidar support:   no       MSG support:               no       GRIB support:              yes       EPSILON support:           no       WebP support:              no       cURL support (wms/wcs/...):no       PostgreSQL support:        no       MySQL support:             no       Ingres support:            no       Xerces-C support:          no       NAS support:               no       Expat support:             no       libxml2 support:           yes       Google libkml support:     no       ODBC support:              no       PGeo support:              no       FGDB support:              no       MDB support:               no       PCIDSK support:            internal       OCI support:               no       GEORASTER support:         no       SDE support:               no       Rasdaman support:          no       DODS support:              no       SQLite support:            no       PCRE support:              yes       SpatiaLite support:        no       DWGdirect support          no       INFORMIX DataBlade support:no       GEOS support:              no       Poppler support:           no       Podofo support:            no       OpenCL support:            no       Armadillo support:         no       FreeXL support:            no       SOSI support:              no                 SWIG Bindings:             no            Statically link PROJ.4:    no       enable OGR building:       yes       enable pthread support:    yes       enable POSIX iconv support:yes       hide internal symbols:     no      And doing the command results in this: ~/gdal$ ogr2ogr -f "ESRI Shapefile" /media//Geodata/testing/sosi.sos /Documents/test.shp FAILURE: Unable to open datasource `/Documents/test.shp' with the following drivers. I am fairly inexperienced with this, so I am trying to learn. What is going on here, what am I doing wrong?


title: Problems with ESRI JavaScript API

text: I'm having some problems to get started using the ESRI JavaSript API. I have a test map service published on my notebook but I can't get it to display. I don't have problems with esri hosted services, just the ones I'm trying to host. Here is some code:               <script src="http://serverapi.arcgisonline.com/jsapi/arcgis/?v=2.1" type="text/javascript"></script>     <link href="http://serverapi.arcgisonline.com/jsapi/arcgis/2.1/js/dojo/dijit/themes/claro/claro.css" rel="stylesheet" type="text/css" >     <script type="text/javascript">         dojo.require("esri.map");         dojo.require("esri.geometry");         dojo.require("esri.layers");         var map;         function init() {             esri.config.defaults.io.proxyUrl = "/proxy/proxy.ashx";             map = new esri.Map("mapDiv");             // var basemapURL = "http://server.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer";             var myURL = "http://brpnb010/arcgis/rest/services/teste/MapServer";                  // var basemap = new esri.layers.ArcGISTiledMapServiceLayer(basemapURL);             var myLayer = new esri.layers.ArcGISTiledMapServiceLayer(myURL);                  // map.addLayer(basemap);             map.addLayer(myLayer);                  dojo.connect(map, "onExtentChange", showExtent);         }              function showExtent(extent) {             var s = "";             s = "XMin: " + extent.xmin +             " YMin: " + extent.ymin +             " XMax: " + extent.xmax +             " YMax: " + extent.ymax;                  dojo.byId("onExtentChangeInfo").innerHTML = s;         }              dojo.addOnLoad(init);     </script>      I'm using Visual Studio to develop, and when the page loads, I get a "Could not load cross domain resources" error. **edit** I'm using version 2.1 of the Javascript API. And I'm developing on localhost, which has the brpnb010 name. Any ideas? Thanks, George


title: What is the difference between ESRI's SDE and Geospatial dbs?

text: When would you want to use ESRI SDE (I think it is called Server but I forget the level - please feel free to edit) versus a spatially enabled database. What are the trade-offs on either side? What are the benefits on either side?


title: What are the differences between spatial databases?

text: What is the difference between Oracle Spatial, PostGIS, SQLServer2008, SpatiaLite, and MySQL in terms of their spatial support. Are there specific use cases where each shines. What might you want to avoid with any of these databases.


title: Good conferences for Spatial/GIS developers

text: What are good conferences for GIS and spatial developers to attend. For ESRI desktop users I think the User Conference is a must attend but if you are an ESRI developer I think the ESRI developer conference is a better event. If you can't make the big shows then the regional events are also a really good deal. FOSS4G is great when it is in your neighborhood but that is a big if.


title: How can I print a bunch of Arcmap mxd's all at once?

text: I have a folder with a dozen Arcmap compositions (.mxd files) I need to print. They're all ready to go, I don't need to do any editing. Do I really have to open up each one, wait to refresh, and press print?


title: Can ArcGIS and OSGeo4W share the same python install?

text: I don't like unnecessary duplication. Is it possible to install python only once and have both ArcGIS and OSGeo4W use it?


title: Allow multiple users to share personal geodatabases safely?

text: We're using ArcGIS 9.3 but for reasons beyond the scope of this question we have to stick with personal geodatabases (MS Access .mdb's) and multiple users. How can we set this up so only a limited few have write access and still allow the bulk of the users to use it?


title: Calculating travel time from residential points to stores via public transit using ArcGIS for Desktop?

text: I don't know how to summarize it in just one title, but here is the problem: My goal is to calculate the travel time from residential points to stores via public transit. Basically this period of time consists three parts:   1. The walking time from a residential point to the nearest bus stop "t1", which can be calculated by the Nearest Facility tool from Network Analyst. The resulting table is called "Resident_Stop", of which the number of records is m;   2. The bus-riding time "t2" from this bus stop to the nearest bus stop of the target store via bus route, which is a record from Network Analyst's OD Time Matrix table, namely, "Stop_Stop". The origins and destinations are all bus stops on one single bus route so that the number of records is s*s;   3. The walking time from the ending bus stop to the store "t3", which is basically the same with the first step. The table is "Stop_Store", of which the number of records is n. So the total travel time "t" equals to (t1+t2+t3) and here is where problem occurs: I haven't figured out a way to attach the three tables reasonably to add the times on using a join. Currently I have an OD Time Matrix table of the walking time directly from residential points to stores via road network (no transit involved) and its number of records is m*n. I also separatly joined the table "Resident_Stop" in step 1 and the table "Stop_Store" in step 3 to the m*n-sized table based on residential point IDs and store IDs. So the remaining thing now is to assign the bus-riding time in table 2 based on table 1's stop IDs and table 3's stop IDs to each of the m*n records. Then I can get the total time. Should you get any ideas, I'm willing to discuss with you. I'm now thinking about using Python to handle this but I'm a beginner on that. If you have any books or tutorials to recommand please tell me.


title: provide an alternative "label" when styling vector layer with openlayers

text: I'm using openlayers 2.12 and I am styling a vector layer using your typical `new OpenLayers.Style({` Now, this layer does two things. First , renders polygons from Geoserver, so in the `style` I have `'label': '${p_name}'`. Works fine. The second thing the layer does is providing tools to the user to edit an already existing polygon. So when a user chooses to edit a polygon, (resize it for example) the familiar little points appear , as handlers. Now, these points/handlers render "undifined", as name because the `'label': '${p_name}'` is actually undifined. Did not came from Geoserver. Is there a workaround? How can I set an alternative name for the handlers? Thanks in advance


title: Determining projection of Texas county shapefile?

text: I have a Shapefile of San Jacinto County in Texas. It thinks it is NAD83 UTM Zone 15N, but the dataset isn't jiving with any other dataset. Max X: 1,073,721 Max Y: 11,223,970 Min X: 895,728 Min Y: 11,013,769 I don't usually work with datasets from this part of the world, so I don't know which projections are in common use. Anyone who is familiar with the area know what projection these extents might fit? This is the county I'm talking about: ![enter image description here](http://i.stack.imgur.com/E5iIh.png)


title: Setting default browser for QGIS actions

text: I am trying to open a pdf using QGIS action tool with Adobe Reader. Currently, the action defaults to Internet Explorer which doesn't succeed in opening my document. I'm currently able to open AutoCAD files and photos. Does anyone know how to set the default software used to by QGIS actions? I'm using QGIS version 2.0. my pathname in the attribute field is: explorer C:\Documents and Settings\user\Bureau\file.pdf Action Type: Windows. Name: Open file. Action: [%"pathname"%]. I didn't click the Capture checkbox. Thanks in advance!


title: How do I create a static popup in OpenLayers that contains HTML content?

text:
I would like to create a static popup using openLayers 2.11 which is always in the lower-right corner of the screen and appears on map load. The idea behind this is to have a logo for my school always appear. I already have a marker layer setup. Maybe I am approaching this incorrectly. Any help would be appreciated. I have tried using some of the examples on the OpenLayers website regarding anchored popups with no luck. Thank you in advance from a poor, starving student. :) The contents of the popup are as follows:
         <link rel="stylesheet" type="text/css" href="css/cube.css">

front

back

top

bottom

left

right


title: How does a GPS/GSM/GPRS Tracker work without a monthly fee?

text: GPS newbie here. Looking at this product: http://dx.com/p/gps-gsm-gprs-tracker-for-personal- remote-positioning-29284 How does a GPS tracker work without a monthly plan? Thank You.


title: Does OSM2PO take into consideration turn restrictions?

text: It is supposed that for using Shooting star function in pgrouting, we need to have RULE and TO_COST fields in our table, but OSM2PO hasn't created those fields. what must i do for OSM2PO takes into consideration the turn restrictions? Thank you! (sorry, my english is so poor...)


title: Info from the getResultImageLayer missing

text: I was finally able to get the result layer overlaid on my application but now I am facing another problem. The output of the layer url does not show any field or attribute details. I am not sure why this is happening? The output looks like this, ![enter image description here](http://i.stack.imgur.com/thvkY.png) As you can see that it is only describing the Extent and other Map properties but the attribute details, which are visible in the jobs folder of the newly created final_shp, is missing on the url and so I am not able to fetch anything onto the map. Any Idea how I can resolve this?


title: Creating tiles in Mapnik

text: I have installed mapnik and successfully tested it. Now I am facing problem to create my own customized tiles. I have google it but everywhere it was telling about osm tile generation. My vision is simple. I have a huge map data of my target area. I want to create my own customize map in Mapnik/tilemill. Then create tiles and upload them into server. A step by step instruction will be helpful to me as I am newbie in Mapnik. This my main problem. How will I create tiles server of my own data??? Thanks in advance


title: How do i get x, y coordinate values when click on map?

text: I am very preliminary in arc gis and in fact i am a web programmer. I created a local serviece and use this to create a web applcation in ArcGIS > ArcGIS Server for the Microsoft .NET Framework > ArcGIS Server Manager.exe . Then i copy website from "c:/inetpup/wwwroot/projectName" and open this to VisualStudio. Now i want when i click on map, display an pushpin on click location and send x, y to server.


title: How to add point to a polyline?

text: I have a `polyline`, and i want to add `point (vertex)` on a specific distance, like two meters for example. How can i do that?


title: How can I make the layout size 36" x 72" in ArcGIS?

text: I need to make a poster to present at a conference and I'm wondering if there's anyway to do this all in ArcGIS 10.1. Is it possible to change the layout/page size to 36x72 inches? Is it practical or even feasible to do this? Thanks! Andy


title: Why is arcpy.mapping.ConvertWebMapToMapDocument slow?

text: We use arcpy.mapping for print tasks. The slowest part in workflow is arcpy.mapping function **ConvertWebMapToMapDocument** ; it takes 10-12 seconds each time, doesn't matter if map is big or not. We use ConvertWebMapToMapDocument(mapJson, templateMxd) form. templateMxd is empty A3 document with empty data frame and north arrow. Function ExportToPNG takes only 2 seconds while ConvertWebMapToMapDocument eats most of the total printing time. What can be tips to improve performance here? **Edit: ------------------**   Looks it's not possible to solve it for myself as ConvertWebMapToMapDocument implementation spends a lot of time in unnecessary web requests.   Simplest example:               import time     import arcpy          templateMxd = 'C:\\Program Files\\ArcGIS\\Server\\Templates\\ExportWebMapTemplates\\A3 Portrait.mxd'     mapJSON = '{"layoutOptions": {"titleText": "Simple WebMap JSON example"},"operationalLayers": [{"url": "http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer","visibility": true}],"exportOptions": {"outputSize": [1500,1500]},"mapOptions": {"extent": {"xmin": -13077000,"ymin": 4031000,"xmax": -13023000,"ymax": 4053000}},"version": "1.4"}'          timer = time.time()     aTuple = arcpy.mapping.ConvertWebMapToMapDocument(mapJSON, templateMxd)     print('ConvertWebMapToMapDocument took %.2f s' % (time.time() - timer))           Output:               ConvertWebMapToMapDocument took 2.10 s      Now where all this time went? I used Fidler2 to track web requests and got following:               http://maps1.arcgisonline.com/ArcGIS/rest/info?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer/layers?&f;=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer/legend?&f;=json     http://maps1.arcgisonline.com/ArcGIS/rest/info?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/info?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer/layers?&f;=json     http://maps1.arcgisonline.com/ArcGIS/rest/info?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/info?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer?f=json     http://maps1.arcgisonline.com/ArcGIS/rest/services/USA_Federal_Lands/MapServer/layers?&f;=json          Aggregate Session duration: 00:00:01.687      Only first 4 requests are actually needed other requests are duplicates. For these 4 requests               Aggregate Session duration:    00:00:00.683 Such optimization would save 1 second of 2.10 s spent in ConvertWebMapToMapDocument function.


title: Is there a way to clip a raster using a vector polygon without having Spatial Analyst?

text: I have ArcView 10 and I'm trying to clip a DEM (TIFF) to a wilderness boundary (polygon). I know ArcView won't let me do this without Spatial Analyst, but I don't have $2500. I tried exporting my raster and vector in QGIS, but for the life of me, I can't get the the two to line up (different coord systems). I've tried to "DEFINE COORD SYSTEM" in Arc Toolbox, but something just isn't right... I'm about to go insane after 5 hours of failures... Can anyone alleviate this pain??? PLEASE


title: Editing a KML file or converting KML to Editable format with attributes

text: I'm using QGIS. I'm trying to either edit a .kml file someone else created in Google Earth or convert it to another editable format (like a .shp file). I've read the following threads without really being able to figure out what the solution is: Any luck with qgis2kml? (seems like this plugin doesn't exist? couldn't find it) How to convert KML to shape file without losing attributes? (I tried the zonums converter mentioned in the answers section, but I lose most of the data in the original .kml file) Converting kml to shape file via ogr2ogr including all attributes (I'm enough of a newbie that I didn't really understand this thread) I've also tried the kml to csv converter here http://www.gps-data- team.com/convert.php. Again, it strips most of the data, giving me only the points and names. Edited to add: saving my .kml file as a .shp file does not work. I get the following error: Export to vector file failed. Error: creation of data source failed (OGR error:Failed to open shapefile .shp/.shp. It may be corrupt or read-only file accessed in update mode. ) Here is a snippet from the beginning of the file. I'm guessing by the comments below, it's just really poorly organized. You have these pieces with the pushpin markers:                   <?xml version="1.0" encoding="UTF-8"?>     <kml xmlns="http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2" xmlns:kml="http://www.opengis.net/kml/2.2" xmlns:atom="http://www.w3.org/2005/Atom">     <Document>     <name>Libanius' Letters, Cribiore.kml</name>     <open>1</open>     <Style id="sn_ylw-pushpin00">     </Style>     <Style id="sh_ylw-pushpin2">         <IconStyle>             <scale>1.2</scale>         </IconStyle>     </Style>     <StyleMap id="msn_grn-pushpin">         <Pair>             <key>normal</key>             <styleUrl>#sn_grn-pushpin0</styleUrl>         </Pair>         <Pair>             <key>highlight</key>             <styleUrl>#sh_grn-pushpin</styleUrl>         </Pair>     </StyleMap>     <Style id="sn_wht-pushpin">         <IconStyle>             <scale>1.1</scale>             <Icon>                 <href>http://maps.google.com/mapfiles/kml/pushpin/wht-pushpin.png</href>             </Icon>             <hotSpot x="20" y="2" xunits="pixels" yunits="pixels"/>         </IconStyle>     </Style>     <Style id="sn_grn-pushpin0">         <IconStyle>             <scale>1.1</scale>             <Icon>                 <href>http://maps.google.com/mapfiles/kml/pushpin/grn-pushpin.png</href>             </Icon>             <hotSpot x="20" y="2" xunits="pixels" yunits="pixels"/>         </IconStyle>     </Style>      .... Then about halfway through, there's a folder tag that contains the place marks :                   <Folder>         <name>Libanius' Letters, Cribiore</name>         <description>Exported on Mon Oct 01 2012 16:11:04 GMT-0700 (PDT)</description>         <Placemark id="18">             <name>Antioch</name>             <description>Home</description>             <styleUrl>#msn_star</styleUrl>             <ExtendedData>                 <Data name="searchrename">                     <value>Antiochia ad Orontem</value>                 </Data>                 <Data name="en_name">                     <value>Theoupolis</value>                 </Data>                 <Data name="featuretyp">                     <value>settlement</value>                 </Data>                 <Data name="pid">                     <value>658381</value>                 </Data>                 <Data name="gr_name">                     <value>Ἀντιόχεια ἡ Μεγάλη</value>                 </Data>                 <Data name="la_name">                     <value>Theoupolis</value>                 </Data>                 <Data name="timeperiod">                     <value>AHRL</value>                 </Data>                 <Data name="perseus_li">                     <value>http://www.perseus.tufts.edu/hopper/text?doc=Perseus%3Atext%3A1999.04.0006%3Aalp</value>                 </Data>                 <Data name="wiki_link">                     <value>http://en.wikipedia.org/wiki/Antioch</value>                 </Data>                 <Data name="path">                     <value>/places/658381</value>                 </Data>                 <Data name="x1">                     <value>36.15818</value>                 </Data>                 <Data name="y1">                     <value>36.202097</value>                 </Data>             </ExtendedData>             <Point>                 <coordinates>36.15818,36.202097,0</coordinates>             </Point>      ..... I tried deleting all the pushpin crap in Oxygen and then reopening the layer in QGIS. The new kml file opens, and the points are there with the data I want. Still not editable, and when I try to save as shapefile I get: "Export to vector file failed. Error: creation of layer failed (OGR error:Failed to open Shapefile `LibaniusShapefile.shp'. ) I think I'm not sure what a well-formed KML file should look like. Thanks.


title: Why does Kernel Density fail in model with ERROR 000865?

text: I am attempting to build a model using variables for automation. (see below) ![Model](http://i.stack.imgur.com/HkYDk.jpg) The problem I am having is when the model attempts to run the Kernel Density. It seems as if the Kernel Density will not recognize %FilePreFix%_ts_cd even though a shapefile is created. Below is the error I get. ![Error](http://i.stack.imgur.com/R7Lcg.jpg)


title: TileMill raster-colorizer

text: I am trying to use the raster-colorizer in Tilemill (v0.10.1) but it is giving me an error. I did set in the advanced section the band to `band="1"`. What else do I have to do in order to use raster-colorizer? ![enter image description here](http://i.stack.imgur.com/gvy3G.png)


title: ArcMap lagging when I use the measure tool

text: ArcMap is lagging when I click the measure tool, my data isn't particularly big (biggest .shp is ~150mb) and my PC is brand new. Any idea what may be causing this?


title: Raster Analysis - Identifying features from Raster

text: I have an IDW, tension spline and hill-shade rasters that show an archaeological landscape and from these rasters I want to identify new features which are seen in the landscape as depressions. And from this raster I want to calculate the size and shape of these features. I have not yet come across a form of analysis that would do this, or I am just not understanding the different types of analysis I have looked at. I have been using GRASS, QGis, ArcMap/Secne and surfer, but I'm willing to give anything a try. Thanks Christina


title: How can I select which interpolation technique will be used in a gdal command? (gdal_fillnodata.py)

text: I work on Landsat 7 slc-off images and I would like to use the gdal_fillnodata.py command to fill-in the missing pixels. The problem that I have is that the gdal_fillnodata.py uses the "inverse distance weighting" (IDW) interpolation technique and I want to use the kriging geostatistical technique. Do you have any idea how can I do it? Thank you!


title: How to use IDatabaseCompact.Compact() on an used IWorkspace, without permission denied

text: I am writting an ArcMap-AddIn with vb.net. I get an error when trying to `compact` my IWorkspace(mdb), the error is: > COMException > > You tried to open a database, which was already opened by the user 'ADMIN' > on Computer 'XXXXXX'. Try it again when the database is available. > > On ESRI.ArcGIS.Geodatabase.IDatabaseCompact.Compact() on > MyProject.MyClass.CompactGDB(IWorkspace pWS) How can i access that workspace, which i used before, to compact it? The Code:               ' CompactGDB         Public Sub CompactGDB(ByVal pWS As IWorkspace)                  Dim pDatabaseCompact As IDatabaseCompact             If (TypeOf pWS Is IDatabaseCompact) Then                 pDatabaseCompact = CType(pWS, IDatabaseCompact)                 If (pDatabaseCompact.CanCompact) Then                     Try                         pDatabaseCompact.Compact()                     Catch ex As Exception                         MessageBox.Show(ex.type & ex.Message & ex.StackTrace, "Error", MessageBoxButtons.OK, MessageBoxIcon.Error)                     End Try                      End If             End If              End Sub      I've got 8 other functions, which work with my workspace. Any suggestions?


title: Localized Linear Histogram Match using python

text: I would like to use the local linear histogram match for filling the missing pixels of a Landsat 7 image with another older image (SLC Gap-Filled Products). Do you have any idea how can I do it by using the gdal commands? Thank you!


title: What kind of different results or what is the difference between two URL based WMS requests

text: I am trying to retrieve some WMS layers from geoserver so that I can display them in leaflet and I need to populate a URL field and I noticed that many of the layers that I opened up in openlayers had "GetMap&layers" vs examples that were successful with "GetFeatureInfo&Layers" when being served via URL.


title: What is this PostGIS query doing to show great circle connections?

text: I'm trying to draw a flight map in CartoDB. My starting point is a single lat long point. We'll pretend it is this:               ST_GeomFromText('POINT(-71.060316 48.432044)', 4326)      And I want draw lines between that point and the points in my CartoDB table. I had a hard time even getting started, but someone pointed me to this map: http://bl.ocks.org/andrewxhill/raw/8406976/ Which is fancier than what I need (I don't need to redraw everytime you click someplace) but includes the basic query:               SELECT ST_Transform(ST_Segmentize(ST_MakeLine(        ST_Transform(the_geom, 953027),        ST_Transform(               (SELECT ST_PointOnSurface(the_geom)                FROM lower_48_zips                WHERE cartodb_id = "+data.cartodb_id+"), 953027)), 100000), 4326)         the_geom FROM us_ski_areas WHERE mindist < 1000      Someone also pointed me towards an example which uses a very similar query:               UPDATE experimental.airroutes     SET the_geom =     (SELECT ST_Transform(ST_Segmentize(ST_MakeLine(            ST_Transform(a.the_geom, 953027),            ST_Transform(b.the_geom, 953027)          ), 100000 ), 4326 )      FROM experimental.airports a, experimental.airports b      WHERE a.id = airroutes.source_id         AND b.id = airroutes.dest_id     );      But I'm having trouble pulling either of these apart. I'm decently conversant in MySQL but totally new to PostGIS and I would love help taking apart what is happening here and replacing source_id with a specific point. I thought I could do something like               UPDATE my_table     SET the_geom =     (SELECT ST_Transform(ST_Segmentize(ST_MakeLine(            ST_Transform(my_table.the_geom, 953027),            ST_Transform(ST_GeomFromText('POINT(-71.060316 48.432044)', 4326), 953027)          ), 100000 ), 4326 )      FROM my_table     );      But I know I need a "where" clause in there and I'm not even sure what it should be.


title: Stitching MapProxy tiles to receive one georeferenced image?

text: In my struggle to get a decent aerial image for playing around with as a layer in Biosphere3D (and for other purposes as well), I managed to install Mapproxy and then cached the area in my desired zoom level from a 2m-Orthophoto WMS through the proxy within QGIS.. Well, now I am in possession of tiles worth 60MB, but i see no way how to stitch them together to one image. The tiles are stored in a folder-structure like this: ![enter image description here](http://i.stack.imgur.com/eRGsK.png) The following sets of tiles are in 32/000/781/000/25, 32/000/781/000/25 etc. containing 408.png to 431.png I have no clue how to interpret this structure. I wonder if there is any way (script) to stitch those tiles to receive one georeferenced image. The tiles were captured in WGS84. Grateful for any hint on how to proceed.


title: How can I determine the location of a visitor to my website?

text: If I want to provide a website with location appropriate content, what ways can I use to find the visitor's rough location using their IP or header information?


title: pgrouting download link not working

text: I am trying to download pgrouting but the download links are not working. I have tried many links but no link is working. can you please help me in this regard. will be thankful


title: Create Materialized View to Flatten ArcGIS JavaScript API Database

text: I'm doing a project using ArcGIS JavaScript API v3.5 hosted on ArcServer v10.1. While reading the question below, remember that I **do not** have control over the datasource and cannot just create a materialized view in the database to solve the issue. **Problem:** The issue I'm having is that parts of the database I'm querying are very very normalized -- one table has 11 Foreign Keys that point to parent tables that are even more normalized. You can imagine that to join all the information for a single record for all these relationships is cumbersome. It requires many requests and lots of looping and checking which is slow and error prone. **Questions:** I have 2 questions:   1. What other options are there in the JS API to make querying easier from the client? I see layer types like QueryDataSource in the API reference that seem to be what i want. But there are no examples of how that works and no one online seems to be using that solution. Has anyone here used QueryDataSource?    2. I'm open to options on how to request that the datasources be modified. The most simple solution is to make them create a flattened table joined with all the data. What other options are there? Since I don't know the ArcGIS Server aspect as well I'm a little blind in what to request. For example, is there a way in ArcMap to create materialized-like views for a normalized datasource? Another possible solution is that the API should gather all related records up for me and spew out one big JSON object when i request a child table with many Foreign Keys. This is how I create API's in Django. All related records become nested hashes of a record. Is there someway to construct your API like this with ArcServer? I can't wait to hear what thoughts/solutions others have for this problem


title: Displaying DEM raster in QGIS

text: Using Raster>Analysis>DEM to display hillshade produces the message "gdaldem.exe had stopped working" so I asked the community. I'm new to QGIS, and I now realise after a little experimentation that I should have used Raster>Terrain Analysis>Hillshade which now works OK - just in case another newbie has encountered the same problem. Why the former action should crash gdaldem.exe remains a slight concern though.


title: How to convert Shapefile to TopoJSON with grouping neighbor shapes by some criteria

text: I have a shapefile with countries, which I convert to GeoJSON using `ogr2ogr` and then converting that json to another TopoJSON, using `topojson` command- line tool. Everything works like charm, but now I want to learn to concatenate/merge/join/cluster (just don't know what is the most appropriate term in GIS context) different neighbor shapes to a bigger ones - say, all African countries to one big Africa shape. Let's assume that all Hdata needed (like continent name) is present in the original shapefile. So, what is the simplest way to merge neighbour shapes while converting shapefile to topojson format? Sorry if this question is about something trivial, it's just that this is a new topic to me, though I'm very glad I should deal with such issues.


title: Separate polygons based on intersection using PostGIS

text: I have a PostGIS table of polygons where some intersect with one another. This is what I'm trying to do:   * For a given polygon selected by id, give me all of the polygons that intersect. Basically, `select the_geom from the_table where ST_Intersects(the_geom, (select the_geom from the_table where source_id = '123'))`   * From these polygons, I need to create a new polygons such that intersection becomes a new polygon. So if polygon A intersects with polygon B, I will get 3 new polygons: A minus AB, AB, and B minus AB. Any ideas?


title: The mysteries of the ArcGIS Server Identity - can you help?

text: Just getting back into ArcGIS for .NET again after about 2 years of doing the training course. I have a question about the "Add ArcGIS Identity" that you choose in Visual Studio - it doesn't seem to matter what combination of username/password I think of, it always says I have the an invalid username/password (by the way I am using accounts that are in the agsadmin group) on the AGS server. Yet if I copy a ciphered identity tag from the web.config of some random AGS application I made via ArcGIS Manager - it works!!, i.e. the web app loads up the map services and I can see a map in my application. Now, I can't decipher this identity tag as you can probably guess - but does anyone know what I might be doing wrong or what I need to look, or check for in order to get the correct credentials right. I just can't rest until I have this figured out - instead of doing this hack. Cheers **UPDATE - Answer** I made sure my dev PC had the ArcGISWebServices user set with the same password as the AGS physical server - so the "impersonation" was consistant between the two machines and accounts. I replaced out the ciphered details in the web.config and replaced with the ArcGISWebServices one and it works fine. I also tried doing it via the "Add ArcGIS Identity" and now it works but make sure you put the machine name (i.e. the name of your dev pc) in the "Domain or Machine Name" textbox.**


title: how to retrieve ATMs List from OSM

text: I am using this command to get the ATMs list `osmosis --read-xml my_map.osm --tf accept-nodes amenity=atm --write-xml output-nodes.osm` But it is giving the list of nodes and ways also...I want only nodes and is there any relation for these nodes and ways.


title: What is the best way to protect data online?

text: What techniques do people use when protecting data ( _e.g._ , preventing a user downloading the whole dataset) they are overlaying on top of Google/Yahoo/Bing etc.?


title: PostGIS importing shapefile error w/ official manager

text: I'm having some trouble when trying to import an shapefile of the Sao Paulo Brazil's state (I downloaded it from the web and I'm using the own PostGis manager to import it). Here's an image! ![enter image description here](http://i.stack.imgur.com/vGZQv.jpg) The connection is already on. In the "Options" menu, i had to change from "UTF8" to LATIN1, cause if I didn't, I'd get this error message: "Unable to convert field name to UTF-8 (iconv reports "Illegal byte sequence"). Current encoding is "UTF-8". Try "LATIN1" (Western European), or one of the values described at http://www.gnu.org/software/libiconv/." I also wanted to now what is SRID and wich value should i put there. When I changed to LATIN1, I get this other error message: "Failed in pgui_exec(): ERRO: sequência de bytes é inválida para codificação "UTF8": 0xe3 0xa7 0xe3". In english, it says the byte sequence is not valid to UTF8 codification. I'm a bit desperate cause I've never used GIS before and I'll need to deliver a project using it on the end of the year. Thank you guys since now.


title: Bing maps get Distance between two places

text: I have referred the sample http://help.arcgis.com/en/webapi/silverlight/samples/start.htm#BingRouting and developed a Bing maps routing tool on my Silverlight GIS website.. now my client wants me to also show the distance between two places whos direction is shown.. I took a look at the Routing class but failed to find anything to get the distance.. I used the ‘haversine’ formula, but it does not give me the exact distance along the route.. or Can I calculate the length of the geometry/polyline drawn on the map.. Will this result be accurate? Please guide me so that I can get the distance between my start and destination..


title: Comparison of SDI frameworks

text: What are the main differences between these SDI, especially in the level of geospatial features and ease of integration with existing CMS?   * Cartaro   * Geonode   * GeOrchestra   * EasySDI I wonder for the reasons why SDI software is not so popular?


title: How to restrict portions of road network from the analysis in ArcGIS 10 Network Analyst?

text: I'm trying to sort out network analysis in ArcGIS 10. Toy example of finding shortest distance between two points looks something like that : ![closest facility output](http://i.stack.imgur.com/xpDPM.jpg) Where squares are two locations, and the green line is an output from Closest Facility tool. What I need to do, is to restrict motorways and slip roads from the analysis. So the new route should avouid the red segments: ![restricted](http://i.stack.imgur.com/gRUas.jpg) What is the best strategy to proceed in such case? Should I edit my network and delete all the segments I do not want to have in analysis? Or perhaps use linear Barriers? Thanks for help!


title: Is there a st_buffer inverse function, that returns a width estimation?

text: When you supply a geometry `g1` and a width `w` to the standard OGC function ST_Buffer, it returns a new geometry `g2`,                 g2 = ST_Buffer(g1,w)      This question is about the inverse: a function where you can supply two related geometries -- `g1` and `g2` satisfies ST_Within(g1,g2) -- and returns the width parameter,                 w = Unbuffer(g1,g2)      Now the interpretation of `w` is "average width". So, **where is this function _Unbuffer_ defined? Are there references for it? Is there a standard library for it?** * * * Complementing after bounty. Examples of applications (estimation of): street width, sidewalk width, river width, canal width, internal polygon semelhance, riparian vegetation buffer width, etc. Example of proposals of credible sources: Pavelsky & Smith 2008 ; ArcMap Script Procedure. They are **not** solutions here because are (very) CPU consumer, "river-specialized", and _can not be implemented by OGC standards_. They calculate average of _river widthes_ by each of many perpendicular (of a central line) lines; so, they are "brute force" algorithms. Example of proposal _without_ credible and/or official sources: troubleshooter here. In this case, you can get the bounty finding sources that makes credible the posted solution.


title: Why can't I install python-qgis on Ubuntu 12?

text: I am with ubuntu 12 and Qgis 1.8. I have installed QGIS and grass plugin but I can't install python plugin. How to install python plugin with dependeny in Ubuntu?               ~$ sudo apt-get install python-qgis     [sudo] password for sugandha:      Reading package lists... Done     Building dependency tree            Reading state information... Done     Some packages could not be installed. This may mean that you have     requested an impossible situation or if you are using the unstable     distribution that some required packages have not yet been created     or been moved out of Incoming.     The following information may help to resolve the situation:          The following packages have unmet dependencies:      python-qgis : Depends: python-qt4 (>= 4.1.0) but it is not going to be installed     E: Unable to correct problems, you have held broken packages.


title: tin model in Qgis

text: I am having DSM and DTM data .My primary objective is to develope a CHM(Canopy Height Model) in Qgis.which i have done already but in my company they are stressing that CHM model which i already achieved is in Raster but they want the result in TIN model how to do this in Qgis?


title: "XY to Line" not converting all 58,000 records

text: I have a dbf (and xls) that contains ~58,000 records of latitude and longitude for both start and end points. I would like to connect each XY coordinate pair using the XY to Line tool found inside the Data Management toolbox. I will be using the Great Circle method fyi. The resultant shapefile contains ~16,000 records/polylines. As far as I can tell, there is no difference between processing the table inside or outside a file geodatabase. Anyone have any thoughts as to why it may not be processing the entire table? This issue persists with both the dBaseIV and XLS file. Cheers, Michael


title: Getting Field value for a given coordinate out of ESRI Shapefile using GDAL/OGR

text: I have a bunch of ESRI shapefiles that contain Land Cover data across Canada (provided by Geobase). I need to query the files and obtain the value of a certain field (called COVTYPE) in a given location specified either by lat/lon or UTM coordinates. I can open a file in Quantum GIS and find the needed values manually, but this process needs to be done automatically for a large number of values so I'm using C++ and GDAL. I'm very new to the whole GIS area and I'm just starting to learn GDAL. So far I've learned how to loop over layers, features, fields etc, but I can't figure out how to relate coordinates to features. I went through the official GDAL/OGR tutorials, but I'm still stuck. Any pointers would be appreciated! **Update:** Just to clarify. My input data would be coordinates (lat/lon or UTM). So the way I imagine this should work is:   1. I find which feature the given point is located within   2. read the field value in that feature I know how to do step 2, but I'm stuck at step 1. **Final solution** Technically, the update to Ragi's answer comes very close to what I need to do, but I'll post here a quick summary of how I ended up using it. My way is easier for constructing a point geometry object. So leaving all the boilerplate code out, the essential part is:               OGRPoint* point = new OGRPoint(some_longitude, some_latitude);     layer->SetSpatialFilter(point); //getting only the feature intercepting the point     OGRFeature* feature = layer->GetNextFeature();     int value = feature->GetFieldAsInteger(0); //getting the value - DONE!     //in my case it's 1st field, hence the 0th index


My OpenLayers map was working fine before, showing a KML overlay. Suddenly it broke. Now I see this error in my browser error console. > TypeError: 'undefined' is not a constructor (evaluating 'new > OpenLayers.Layer.GML') My code looks like this:               var kmllayer = new OpenLayers.Layer.GML("KML", "./my-kml-file.kml",         {             format: OpenLayers.Format.KML,          formatOptions: {           extractStyles: true,            extractAttributes: true,           maxDepth: 2         }        });     map.addLayer(kmllayer);


title: Can I split up v.in.ogr into multiple v.clean commands?

text: I'm trying to load a large set of polygons into GRASS with " _v.in.ogr_ " and am running out of memory at the "Break Boundaries" stage. I'm wondering if I can break into steps all of the different things that _v.in.ogr_ does to a data set when it is run without the " _-c_ " flag? For example, I know that it must run the equivalent of _"v.clean tool=bpol_ " and " _v.clean tool=rmdupl_ ", but I'm not exactly sure what to run after those 2 to complete the polygon cleaning that normally happens. I figure that I might be able to save some memory if I don't have to run all the operations in one go. Does this makes sense? Thanks. * * * Ok, so it was like digging for gold to find this, but in a presentation given by some people who seem to know GRASS pretty well, they said that the following steps were needed to clean polygon geometry after import: **Cleaning procedures for area import**   * break polygons   * remove duplicates   * break boundaries   * remove duplicates   * clean boundaries at nodes   * change dangles to lines   * remove bridges Link to presentation here: http://geostat- course.org/system/files/GRASS_geostat_landau_2011_intro.pdf * * * Wait, no... there's more to this apparently. I tried those steps listed above (or attempted to at least). I think I'm worse off now than when I started.               v.in.ogr -c     Importing map 52 features...     Number of boundaries: 52     Number of centroids: 52     Number of areas: 42     Number of isles: 42     Number of incorrect boundaries: 10     Number of centroids outside area: 10          1. break polygons       v.clean type=boundary tool=bpol     Number of boundaries: 270     Number of centroids: 52     Number of areas: 2     Number of isles: 4     Number of incorrect boundaries: 268     Number of centroids outside area: 50          2. remove duplicates       v.clean type=boundary tool=rmdupl     Number of boundaries: 190     Number of centroids: 52     Number of areas: 40     Number of isles: 4     Number of incorrect boundaries: 126     Number of centroids outside area: 12          3. break boundaries       v.clean type=boundary tool=break     Number of boundaries: 5147     Number of centroids: 52     Number of areas: 151     Number of isles: 0     Number of incorrect boundaries: 5132     Number of centroids outside area: 40     Number of areas without centroid: 139          4. remove duplicates       v.clean type=boundary tool=rmdupl      Number of boundaries: 4803     Number of centroids: 52     Number of areas: 153     Number of isles: 0     Number of incorrect boundaries: 4781     Number of centroids outside area: 38     Number of areas without centroid: 139          5. clean boundaries at nodes       (No idea what-the-flip this means.  Clean with what?  Do it by hand?)          6. change dangles to lines        v.clean type=boundary tool=chdangle thresh=-1     Number of boundaries: 4803     Number of centroids: 52     Number of areas: 153     Number of isles: 0     Number of incorrect boundaries: 4781     Number of centroids outside area: 38     Number of areas without centroid: 139          7. remove bridges       v.clean type=boundary tool=rmbridge in=eq1_6 out=eq1_7     Number of boundaries: 4420     Number of centroids: 52     Number of areas: 153     Number of isles: 22     Number of incorrect boundaries: 4352     Number of centroids outside area: 38     Number of areas without centroid: 139      Guess I'm off to read the source for _v.in.ogr_ now.


title: ArcSDE 9.3.1 -> Performance

text: We have a number of mapservices in ArcGIS Server 10.0 that reads data from FGDBs. We are now moving that data in to ArcSDE and the plan is to use this data from AGS. Is it better to connect to SDE using Direct Connect or going through the ArcSDE services performance wise?


title: Saga Reclassify grid problem

text: I am using Qgis 2.1 ver. 5f5cd4c (64 bit version on win7 64) OSGeo installed. I was trying to reclassify a DEM in EPSG:32614. I used the SAGA's Reclassify algorithm in the processing toolbox. I defined that I wanted to use a simple table, defined the table and that I wanted the reclassified grid saved to a tif file. I also wanted the resulted reclassified DEM to be opened after the algorithm finished the calculations. This is what I am repeatedly getting: ![enter image description here](http://i.stack.imgur.com/nulBb.png) I tried to open a new project and load the resulting raster manually, but Qgis does not recognize it as a valid raster data source: ![enter image description here](http://i.stack.imgur.com/EWNl0.png) I was able to open the file and convert it to GeoTIFF in SAGA. I wonder if I am missing something in the reclassifying procedure, or something related with the format of the resulting raster. Thanks in advance for any advice. PS. I was able to use r.reclass to do what I wanted. I also tested the saga reclasify algorithm in other two installations and I have the same problem. Gerardo