title: Stratified Random Sampling by ID in QGIS?

text: I've been playing around with the Random Points tool in QGIS to realize stratified random sampling, but I'm running into a problem where I want my strata to be based on a polygon ID, not each individual polygon. In other words: I may have multiple non-contiguous polygons that all represent a single strata, where I want a TOTAL of (e.g.) 10 points randomly distributed among the disconnected polygons, not 10 per polygon. Is there any way to do this? Would a multi-part polygon work, for this?


title: How to explode multipart polygons within shapefile in QGIS 1.8

text: I've got a shapefile with polygons, some of which are multipart polygons. I need to create individual polygons out of the multipart polygons. Is there a tool or plugin that can do this in QGIS 1.8?


title: Its possible to use a c# webmethod with a openLayers.Protocol.HTTP?

text: Hi Im trying to use c# as a backend of a layer:               var layerVehicles = new OpenLayers.Layer.Vector("Dispositivos", {             strategies: [new OpenLayers.Strategy.Fixed(), new OpenLayers.Strategy.Refresh({ interval: 30000, force: true })],             protocol: new OpenLayers.Protocol.HTTP({                 url: "mapa.aspx/obtenerRastreados",                 params: "{ 'idUsuario': '1' }",                 //headers: "application/json; charset=utf-8",                 format: new OpenLayers.Format.GeoJSON()             }),             styleMap: styles         });      and the c# code is:                   [WebMethod(EnableSession = true)]         [ScriptMethod(ResponseFormat = ResponseFormat.Json)]         public static string obtenerRastreados(string idUsuario)         {             int idUsuarioInt = 0;             int.TryParse(idUsuario, out idUsuarioInt);             List<dispositivo> dispositivosMapa = Negocio.NegocioRastreos.ratreoMapa(idUsuarioInt);                  crspropiedades urn = new crspropiedades();             urn.urn = "urn:ogc:def:crs:OGC:1.3:CRS84";                  crs_ crs = new crs_();             crs.type = "OGC";             crs.properties = urn;                  coordenadas coordenadas_ = new coordenadas();             coordenadas_.type = "Point";             decimal[] latlon = { 1.76m, 0.87m };             coordenadas_.coordinates = latlon;                  propiedades propiedades_ = new propiedades();             propiedades_.alias = "AUTO";                  dispositivosMostrar mostrar = new dispositivosMostrar();             mostrar.type = "Feature";             mostrar.id = "auto123";             mostrar.properties = propiedades_;             mostrar.geometry = coordenadas_;             mostrar.crs = crs;                  dispositivosMapa rastreos = new dispositivosMapa();             rastreos.type = "FeatureCollection";             List<dispositivosMostrar> listaMostrar =  new List<dispositivosMostrar>();             listaMostrar.Add(mostrar);             rastreos.features = listaMostrar;                  System.Web.Script.Serialization.JavaScriptSerializer serializer = new                  System.Web.Script.Serialization.JavaScriptSerializer();                  return serializer.Serialize(rastreos);         }      but thelayer never consume the webmethod, I have the followig error: Uncaught SyntaxError: Unexpected token <


title: How to write a signal for 'QgsMapCanvas::xyCoordinates (const QgsPoint & p)'

text: I am very new to QGIS and python. I am writing some python codes and I want to capture the x,y coordinates when the cursor moves on the map layer. (same as the coordinates show in the QGIS status bar) I found that there is a signal which emits the current mouse position. http://qgis.org/api/classQgsMapCanvas.html#abf90fbd211ea419ded7c934fd289f0ab How should I write the code for that signal? I was able to write(found) same like code for the canvas clicked event,               def __init__(self, iface, layer=None):         # Save reference to the QGIS interface         self.iface = iface     # a reference to our map canvas         self.canvas = self.iface.mapCanvas() #CHANGE         # this QGIS tool emits as QgsPoint after each click on the map canvas         self.clickTool = QgsMapToolEmitPoint(self.canvas)          def initGui(self):         result = QObject.connect(self.clickTool, SIGNAL("canvasClicked(const QgsPoint &, Qt::MouseButton)"), self.handleMouseDown)          def run(self):         # set the current layer immediately if it exists, otherwise it will be set on user selection         self.cLayer = self.iface.mapCanvas().currentLayer()         #if self.cLayer: self.provider = cLayer.dataProvider()         # make our clickTool the tool that we'll use for now         self.canvas.setMapTool(self.clickTool)         # show the dialog         self.dlg.show()         # Run the dialog event loop         result = self.dlg.exec_()         # See if OK was pressed         if result == 1:             # do something useful (delete the line containing pass and             # substitute with your code)             pass      Can someone please help me to write same like this code for 'xyCoordinates' signal. I tried as below but not worked,(no errors but not worked)                def __init__(self, iface, layer=None):         # Save reference to the QGIS interface         self.iface = iface     # a reference to our map canvas         self.canvas = self.iface.mapCanvas() #CHANGE         # this QGIS tool emits as QgsPoint after each click on the map canvas         #self.clickTool = QgsMapToolEmitPoint(self.canvas)     self.curMove = QgsMapCanvas(self.canvas)          def initGui(self):         result = QObject.connect(self.clickTool, SIGNAL("canvasClicked(const QgsPoint &, Qt::MouseButton)"), self.handleMouseDown)         QObject.connect(self.curMove, SIGNAL("xyCoordinates(const QgsPoint & p)"), self.cursorMove)          def cursorMove(self, point):         #print "inside function"         print str(point)          def run(self):         # set the current layer immediately if it exists, otherwise it will be set on user selection         self.cLayer = self.iface.mapCanvas().currentLayer()         #if self.cLayer: self.provider = cLayer.dataProvider()         # make our clickTool the tool that we'll use for now         #self.canvas.setMapTool(self.clickTool)         self.canvas.setMapTool(self.curMove)         # show the dialog         self.dlg.show()         # Run the dialog event loop         result = self.dlg.exec_()         # See if OK was pressed         if result == 1:             # do something useful (delete the line containing pass and             # substitute with your code)             pass


title: Making my GIS application OGC compliant

text: We have a stack of open source software's which we use to build a GIS. We would like to align our application to OGC as much as possible. I visited the site and read about the standards. Now my question is how to start with making the application OGC compliant. There are lot of documents which talk about standards but I couldn't find much on how to get started. Specifically given that am building on top of Open Source software's and Implementing is very subjective word. Can you guide me to some resources which can get me started? Ideally we would like address Web Mobile and Desktop, We have a existing solution for Desktop but are planning to go ahead with Web. We are aiming for a GIS product. Here is the list of Open Source Products: GDAL/OGR GeoServer OpenLayer Cesium


title: How to create a point along a line given distance

text: I have been given a non-spatial database which contains the names of lines which I have a shapefile for. The database also has names of points and the distance of each point from the origin along the route of these lines. How would I create a point shapefile which locates each of the points? Of course I would like to use QGIS, or GRASS and as a last resort I can get access to ArcGIS. I expect I will have to carry out a step to specify which end of the line is the origin.


title: How to create equidistant points in QGIS?

text: I am trying to create points (new layer) at specific distance along the road (existing layer) in QGIS. This gives solution for ArcGIS. How to achieve this in QGIS? This one  explains how to create points but does not do anything about the distance. **EDIT** (I applied the proposed solutions with different measures of lengths for I did not know the conversion) @Nathans's solution worked to some extent, I got... ![enter image description here](http://i.stack.imgur.com/AJLAe.png). Here, the projection of these equidistant points is different from the original line. With @underdark's suggestion, I got ![this image](http://i.stack.imgur.com/JdcsQ.png) where the points does not seem to be equidistant. I guess there is some projection issue with both of these which I am not understanding. Thank you very much.


title: OpenLayers spatial filter polygon problem

text: I've got a dataset of English Counties (currently cut down for testing) that I'm serving from SQL 2008 using MapServer and OpenLayers. The features in my dataset are appearing in the right location on my base map. I'm now trying to carry out a selection \ filter using a polygon that the user draws on screen. Set up is as follows: County Layer                var CountyLayerDay2 = new OpenLayers.Layer.Vector                 ("CountyDay2", {                     isBaseLayer: false,                     projection: new OpenLayers.Projection("EPSG:27700"),                     displayProjection: new OpenLayers.Projection("EPSG:27700"),                     strategies: [new OpenLayers.Strategy.Fixed()],                     protocol: new OpenLayers.Protocol.WFS({                         url: "http://localhost/cgi-bin/mapserv.exe?map=C:/DevProjects/fgsv2/mapfile.map",                         version: "1.0.0",                         featureType: "CountyDay2",                         featurePrefix: "ms",                         featureNS: "http://mapserver.gis.umn.edu/mapserver"      My polygon drawing layer is then added as follows:                var polyLayer = new OpenLayers.Layer.Vector                 ("Polygon",{projection: new OpenLayers.Projection("EPSG:27700"),                             displayProjection: new OpenLayers.Projection("EPSG:27700")                            });            When the user draws a polygon, it is being reported as being in the right place with regards to the coordinates reported in the "Post" tab of Firebug (user has drawn round the West Midlands):               <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.0.0" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.0.0/WFS-transaction.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">     <wfs:Query typeName="ms:CountyDay2" xmlns:ms="http://mapserver.gis.umn.edu/mapserver">     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">     <ogc:Intersects>     <ogc:PropertyName>the_geom</ogc:PropertyName>     <gml:Polygon xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:4326">     <gml:outerBoundaryIs>     <gml:LinearRing>     <gml:coordinates decimal="." cs="," ts=" ">-2.4628747596580483,52.856622995202905 -1.4081872598048397,52.77694589865457 -1.2983239785701848,52.01271002627193 -2.3749841346702887,52.32267900571302 -2.4628747596580483,52.856622995202905</gml:coordinates>     </gml:LinearRing>     </gml:outerBoundaryIs>     </gml:Polygon>     </ogc:Intersects>     </ogc:Filter>     </wfs:Query>     </wfs:GetFeature>      When I look on the "Response" tab I get some strange results. The initial gml:boundedby tag looks ok as the coordinates sit in the right place in terms of British National Grid coordinates as they correspond to the bounding box of the West Midlands county. However, when the response get to the gml:boundedby in the features, the coordinates hav expanded and I end up with a multi part polygon being returned as other counties are also selected:                xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"        xsi:schemaLocation="http://www.opengis.net/wfs     http://schemas.opengis.net/wfs/1.0.0/WFS-basic.xsd                         http://mapserver.gis.umn.edu/mapserver http://localhost/cgi-bin/mapserv.exe?map=C:/DevProjects/fgsv2/mapfile.map&SERVICE=WFS&VERSION=1.0.0&REQUEST=DescribeFeatureType&TYPENAME=ms:CountyDay2&OUTPUTFORMAT=XMLSCHEMA">           <gml:boundedBy>             <gml:Box srsName="EPSG:27700">             <gml:coordinates>386080.398996,272232.103705 439261.802108,307214.500253</gml:coordinates>         </gml:Box>       </gml:boundedBy>     <gml:featureMember>       <ms:CountyDay2 fid="CountyDay2.4">         <gml:boundedBy>             <gml:Box srsName="EPSG:27700">                 <gml:coordinates>386080.398996,272232.103705 475119.498244,576159.994876</gml:coordinates>             </gml:Box>         </gml:boundedBy>         <ms:msGeometry>         <gml:MultiPolygon srsName="EPSG:27700">         <gml:polygonMember>           <gml:Polygon>             <gml:outerBoundaryIs>      Can anyone spot what I've done wrong? My first concern was the conversion from lat long to British National Grid between the "Post" and "Response" tabs but it looks to work initially. I'm new to OpenLayers etc so may have set my layers up wrong. I've searched this site and been through the OpenLayers examples and can't find anything similar to what I need to do. Any pointers would be appreciated.... **EDIT** I'm thinking that it might be a coordinate conversion problem from when the orginal user drawn polygon is created using SRSName:4326 and is then converted to SRSName:27700. I've tweaked a couple of things in my map setup to make sure everything is set up in SRSName: 27700, so the polygon drawn by the user is now pretty much the only thing that references SRSName:4326.


title: Why does Mapserver (ms4w) not generate tmp images?

text: I am having trouble with my web-based application, where my tmp images are not being stored for later retrieval. The images are generated each time they have to be loaded, so zooming out and then back in to the same location requires mapserv.exe to re-generate the tiles for that view. The application lives here: http://rinnerweb.arts.ryerson.ca/moot2 A sample map file snippet is below:                   OUTPUTFORMAT                                        #-> BEGIN OUTPUT FORMATTING DEFINITION             NAME "png"                                      # Name of format             DRIVER AGG/PNG                                  # Image driver             MIMETYPE "image/png"                            # Type of image             IMAGEMODE RGBA                                  # Image mode             EXTENSION "png"                                 # Image filetype             TRANSPARENT ON                                  # Transparency         END                                                 #-> END OUTPUT FORMATTING DEFINITION              WEB                                                 #-> BEGIN WEB DEFINITION             IMAGEPATH "/ms4w/tmp/ms_tmp/"             IMAGEURL "/ms_tmp/"         END                                                 #-> END WEB DEFINITION      The web page's source will show the OpenLayers init code, but I've pasted some sample layer init code below:               hii_1 = new OpenLayers.Layer.MapServer("Land Cover","../../cgi-bin/mapserv.exe",{         map:'C:/ms4w/Apache/htdocs/hii/hii_landcover.map'},{         isBaseLayer:false,         transparent:true,         format:"image/png",         alpha:true,         displayInLayerSwitcher:false         });      Any ideas?


title: How to derive vehicle paths from timestamped location data?

text: I've got a bunch of location data for snow plows. There's a unique id for every plow, a timestamp, and x & y coordinates. Instead of knowing where hundreds street plows were at one minute intervals, I'd like to know which streets have been plowed. I know there are ways of turning points to polylines, but does anybody know of a method to snap the points to the street grid? In other words, of deriving vehicle paths that hug the streets? The reason I think that turning points to polylines probably isn't enough is because I only have location data for every plow every 60 seconds. This means there aren't coordinate points for when plows were turning at certain intersection, so if you just drew lines between points some of them would cut through blocks.


title: Remove duplicate nodes in QGIS

text: I am trying to create a polygon feature from an enclosed line feature using the Geometry Tool: Lines to Polygons. When I do this, I get several spikes with the newly created polygon. I then run the "Check Geometry Validity" tool and it notes that I have about seven duplicate nodes. I am trying to delete the duplicate node using the Edit tools: Node Tool, select the node, but have not been successful with deleting the extra (or duplicate) node. I was told to use Fn + Delete on the keyboard (Mac OS), but nothing seems to happen. Please tell me the best way to get rid of these duplicate nodes. Thx, DP at LMN


title: How do I publish a custom .NET GP Tool as a GP Service on ArcGIS Server 10.1?

text: I have a GP Tool that I am able to publish as a GP service at 10.0. I recompiled the DLL for 10.1 and have followed the steps in the Deploying custom .NET and C++ tools page of the 10.1 help doc. The tool works correctly in arcmap. However, when I get to the step where I create the geoprocessing package, or when I try to publish directly as a service, I get error 00151 "Tool is unsupported on server". _Has anyone been able to take a custom GP Tool developed with C# and publish it to ArcGIS server 10.1 as a GP service withoutwrapping it with python?_ (I do see this thread about difficulties with an executable - which I presume means an exe. I would like to avoid having to wrap all my tools in python, as suggested in that thread.) **Update** I probably should mention I also have desktop installed on the same machine with server (along with Visual Studio). It is interesting that even though I get the 00151 error, when I try to recompile the dll, Visual Studio fails, saying it can't delete the dll because it's locked by another process. It turns out that other process is arcgis server. If I stop the arcgis server instance I can recompile, implying that even though I get the 00151 error, arcgis server is still loading the dll into memory. **Update 2** Esri was able to provide me a simplified version of my gp tool. See the Deploying a Custom C# GP-Tool on 10.1 Server thread in the Esri forum. The project team will be deciding whether or not to port to python before I spend more time on trying to figure this out. Will update again once I get access to the test server.


title: How can I debug a geoprocessing service on the server?

text: I've a geoprocessing script that I'm trying to get working. I've turned on Show Messages, so I can see how far in the script I've gotten before the submitted job fails, but it's not helping, as I'm getting errors that I don't understand exactly, because I don't know where they occurred. Is there any way to debug server geoprocessing scripts on the server? EDIT: Error is occurring in this code, which worked fine locally. I'm not getting either of the messages returned, so something seems to be going wrong with the layers.               for layer in mxdLayers:         try:             arcpy.AddMessage("Layer name is %s" % (layer.name))             layer.visible = visibleLayers[layer.name]         except:             arcpy.AddMessage("Layer not set: " + layer.name)      Error message is Layer: Unexpected error Edit 2: The error is not in the code above. The following code takes the parameter and stores the passed in features in a scratch gdb. Currently this is passed as an empty string, as I haven't gotten to testing this part yet.               # Get the passed point graphics and save to scratch workspace     inPointFset = arcpy.GetParameterAsText(9)     if inPointFset == '#' or not inPointFset:         inPointFset = None     else:         arcpy.CopyFeatures_management(inPointFset, os.path.join(arcpy.env.scratchWorkspace, r"scratch.gdb\Points"))      This code then changes the data source of the first mxd layer (which is why it's blowing up; the first layer doesn't reference the points feature class in the scratch gdb), but since the parameter was an empty string, I'm not sure why this is evaluating to true.               mxdLayers = arcpy.mapping.ListLayers(mapDoc)     if inPointFset:         arcpy.AddMessage("Points")         mxdLayers[0].replaceDataSource(os.path.join(arcpy.env.scratchWorkspace, r"scratch.gdb"), "FILEGDB_WORKSPACE", "Points")


title: Python dataSource Path

text: I have the following code: When I run it:               conType = arcpy.mapping.ListLayers(imxd, "Construction_pt", df)[0]     print conType.dataSource     prints... "Database Connections\Edit - Editor.sde\Edit.SDE.Construction_pt"     and rest of program runs correctly      When another user runs it:               print conType.dataSource     prints... "\Edit.SDE.Construction_pt"     program crashes because cannot find path      My question is why is it only displaying the full path for me and not other users? ArcGIS 10, Windows 7


title: Cause of expected an indented block error in Python script?

text: I have been trying to list all the field names of a shapefile but it is showing 'expected an indented block'. What should I do?               import arcpy          fieldList = arcpy.ListFields("C:\Users\Mridul\Documents\GIS\Spring 2011\650\streets_dev.shp")           for field in fieldList:      print field.name + " is a type of " + field.type + " with a legth of " + field.length          rows = arcpy.UpdateCursor("C:\Users\Mridul\Documents\GIS\Spring 2011\650\streets_dev.shp", "SHAPE_len>400")            cnt = 0          for row in rows:          row.STREETNAME = "OLD " + row.STREETNAME          rows.updateRow(row)         cnt +=  1         print cnt          del row      del rows


title: Remove Polygon Islands

text: I'm really stumped with the problem. I have a polygon layer with no intersects or holes and I need to remove some 'islands', as defined by attributes of a certain size within the layer. It's a hard thing to explain, so refer to the picture below. The different colors represent different attributes within my layer. What I need to do is remove the orange islands of certain size (less than 5000m²) that do not share border with an Orange, Green or Neon polygon. ![Polygon Island Issue](http://i.stack.imgur.com/02JAF.png) I cannot simply select all of the orange polygons under a certain size because that would remove some polygons that make up large areas of polygons. I can't dissolve the orange than remove polygons under a certain size because I will lose the individual polygon shapes which are extremely important. I am really stumped here guys, I'll take a simple solution or any advice. I am pretty handy with python so if you think a programming solution is best let me know.


title: Saving VBA customization to a Normal.mxt problem

text: I have experienced problems with saving VBA customization to a Normal.mxt document within ArcGIS 10. I have a set of VBA macros developed originally within 9.3 version of ArcGIS that reside in a custom toolbar. These VBA forms, modules and customization have been saved within the Normal.mxt in ArcGIS version 9.3 and open and work fine while opening this Normal.mxt in ArcGIS 10. Recently I developed an additional VBA macro (within ArcGIS 10) which works and adds fine to the customized toolbar, however I am unable to save customization permanently to this Normal.mxt in ArcGIS 10. Please note that Customize>Options “Save all customizations to the document” button is unchecked. Despite saving changes to this Normal.mxt in Visual Basic Editor with a new macro/tool added to the customized toolbar, after closing and opening ArcGIS 10 again newly added macro/tool disappears from the customized toolbar. However, the new tool is never saved permanently within the Normal.mxt. It appears that ArcGIS 10 has some issues with VBA customisations. Please could you advice what could be the cause of it and how it could be resolved?


title: Image Pyramid + Geoserver + Google Earth slow with errors

text: I have hundreds of GB of somewhat small tif files. Here is gdalinfo of one of them:               F:\gis\geotiffs>gdalinfo finaloutputSSAM211Jun181306hfwaterfall0.tif     Driver: GTiff/GeoTIFF     Files: finaloutputSSAM211Jun181306hfwaterfall0.tif     Size is 2281, 4301     Coordinate System is:     GEOGCS["WGS 84",     DATUM["WGS_1984",         SPHEROID["WGS 84",6378137,298.257223563,             AUTHORITY["EPSG","7030"]],         AUTHORITY["EPSG","6326"]],     PRIMEM["Greenwich",0],     UNIT["degree",0.0174532925199433],     AUTHORITY["EPSG","4326"]]     Origin = (-85.749319836621979,30.000952072441574)     Pixel Size = (0.000000303816798,-0.000000303816798)     Metadata:      AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=PIXEL      Corner Coordinates:     Upper Left  ( -85.7493198,  30.0009521) ( 85d44'57.55"W, 30d 0' 3.43"N)     Lower Left  ( -85.7493198,  29.9996454) ( 85d44'57.55"W, 29d59'58.72"N)     Upper Right ( -85.7486268,  30.0009521) ( 85d44'55.06"W, 30d 0' 3.43"N)     Lower Right ( -85.7486268,  29.9996454) ( 85d44'55.06"W, 29d59'58.72"N)     Center      ( -85.7489733,  30.0002987) ( 85d44'56.30"W, 30d 0' 1.08"N)     Band 1 Block=2281x1 Type=Byte, ColorInterp=Red       Mask Flags: PER_DATASET ALPHA     Band 2 Block=2281x1 Type=Byte, ColorInterp=Green       Mask Flags: PER_DATASET ALPHA     Band 3 Block=2281x1 Type=Byte, ColorInterp=Blue       Mask Flags: PER_DATASET ALPHA     Band 4 Block=2281x1 Type=Byte, ColorInterp=Alpha      I used gdal_retile.py with the following command line:               gdal_retile.py -v -r bilinear -levels 4 -ps 2048 2048 -co "tiled=YES" -targetDir pyramid --optfile files.txt      I then create an Image Pyramid store in Geoserver and set multithreading to true, and USE_JAI_IMAGEREAD in the layer. When I go to Layer Preview and click on KML it opens in Google Earth. When I zoom in, it takes a very long time to render, and sometimes it times out and gives me various errors. Currently I built the Image Pyramids with only 30 or 145 of those tifs to see if maybe it was because my data set was so large, however I get the same performance with those smaller sets. What am I doing wrong? Here is a picture of one error I am getting: http://i.imgur.com/HHjzu.jpg I also get an error that just says "Error when parsing file" Thanks!


title: How to automate downloading "slices" of data from Web Coverage Service (WCS)?

text: I'd like to automate downloading "slices" of data from a WCS service. This is a WCS 1.1.2 service, supporting netCDF. The data are updated daily - forecasts from an atmospheric regional model - and I'd like to automatically download netCDF "slices" of the data for a fixed geographical bounding box and for a couple of days out into the future. I am quite happy processing the netCDF once I've downloaded it. The obvious brute-force approach is to read the OGC standard and devise the HTTP request myself, and code it up in some scripting language. I suspect I'll end up doing that. But it just feels _wrong_. Can anyone recommend a WCS client library that I can use from code? Or a WCS client that I can run from a command-line? I know I could probably script QGIS or Grass or ArcGIS or something, but that seems too heavy a requirement for this task. Can GDAL do it? Google has turned up the BADC Python WCS library (PDF) (now in OWSlib?) which I will take a look at. Java and GeoTools might be a possibility. And there's the OGC list of WCS clients which I could trudge through. This is on the Windows platform. A .Net library would be ideal :) I haven't linked to the actual WCS service as it's not yet publicly operational. But here is an edited `GetCapabilities` request below, in case anyone's interested in the details               <Capabilities xmlns:xlink="http://www.w3.org/1999/xlink"         xmlns:ows="http://www.opengis.net/ows/1.1" xmlns="http://www.opengis.net/wcs/1.1"        version="1.1.2">       <Contents>       <CoverageSummary>       <ows:WGS84BoundingBox crs="urn:ogc:def:crs:OGC:2:84">       <ows:LowerCorner>+0.0000 -90.0000</ows:LowerCorner>       <ows:UpperCorner>+358.8750 +90.0000</ows:UpperCorner>       </ows:WGS84BoundingBox>       <SupportedCRS>urn:ogc:def:crs:EPSG::4326</SupportedCRS>        <SupportedCRS>urn:ogc:def:crs:OGC:2:84</SupportedCRS>       <SupportedFormat>image/netcdf</SupportedFormat>       <SupportedFormat>application/x-netcdf</SupportedFormat>       <Identifier>------------</Identifier>      </CoverageSummary>


title: TIN to Raster Tool, output elevations are 3.40282e+038 and -3.40282e+038, tried calculate statistics didn't work

text: I have an input TIN file that I created using Create TIN. The elevations are all correct. When I convert the TIN to Raster using the tool, my output elevations are 3.40282e+038 and -3.40282e+038 I've checked projections (and probably do it again) and tried calculating statistics, but so far nothing has changed. Someone else used the same file I was using, output it to the same folder I have been, and they got it to work. EDIT: I edited out all but the end to file names. My S drive is a server drive.               arcpy.CreateTin_3d("C:/...Default.gdb/water_level","PROJCS['NAD_1983_HARN_Adj_MN_Washington_Feet',GEOGCS['GCS_North_American_1983_HARN',DATUM['D_North_American_1983_HARN',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500019.9183995],PARAMETER['False_Northing',100003.9836799],PARAMETER['Central_Meridian',-92.83333333333333],PARAMETER['Scale_Factor',1.000039836799],PARAMETER['Latitude_Of_Origin',44.74583333333334],UNIT['Foot_US',0.3048006096012192]]","B_XS100 Max_Water Hard_Line <None>;Boundary <None> Hard_Clip <None>","DELAUNAY")          # Replace a layer/table view name with a path to a dataset (which can be a layer file) or create the layer/table view within the script     # The following inputs are layers or table views: "water_level"     arcpy.TinRaster_3d("water_level","S:/.../water_level_3.tif","FLOAT","LINEAR","CELLSIZE 100.485944","1")


title: How to calculate a standard deviation across multiple raster layers?

text: I have 15 raster layers. Each one shows accumulated precipitation (rainfalls) for the month of April for the same US state during a given year. So, I cover the period from 1997 to 2012. I would like to have an idea of the variability of precipitation during the last 15 years. I would like for example to calculate the standard deviation at each pixel across the 15 layers. I am guessing I can do it through the raster calculator using sums, subtractions and divisions. But the final formula would pretty long. Is there a more efficient way to achieve my goal? (So far, I am only a QGis user and I don't know Python...but I can start learning if necessary).


title: How to register a Stored Procedure in ArcSDE?

text: I just registered a SQL Server view with ArcSDE through               sdetable -o register [...].       I wonder, how to register a SQL Server stored procedure with ArcSDE? Is it the same process as the view? **EDIT** It seems that to register a stored procedure, it is required to have a "stable" database object: the definition of the stored procedure by itself it is just a piece of code that has no sense until it is executed with some input parameters. Ok, so, if we consider the stored procedure as a "data collection generator", how would you execute this code with the input parameters, catch the tabular data output and register it to ArcSDE? What would you use? Some kind of ArcGIS processing? Please, could you provide any example or guidelines for this process? Thanks!


title: Need to use xlsx data to clip a shapefile

text: I have an Excel table with just one column, with 1150 records, each of them a species name. I also have a shapefile with 24,000 records, each of them a species name and its geographic range. I need to clip out of the shapefile the 1150 records in my Excel table. I thought I could just use model builder's convenient row iterator to do this, but while the row iterator works fine the clip tool will not accept the result as the clip feature input. I am now trying to convert my Excel xlsx file into a shapefile, but since it doesn't have any geographic locations, it's just a list of names, this seems to be impossible. Any advice would be greatly appreciated!


title: How can I get the ellipsoidal values for length and area using spatialite?

text: I have spatialite geometry data stored in a spatialite database in WGS 84 (EPSG 4328) format:               SELECT Srid(geometry) FROM shape_layer;     > 4326      When geometry is being returned as text, it is returned in meters (when the data is being worked on, the projection is more commonly EPSG 3857 to make use of Google Maps/OpenStreetMap). By example:               SELECT AsText(geometry) FROM shape_layer;     > LINESTRING(-381578.202 6574042, -378511.828 6574794.956)      All data works will display correctly (e.g. if using QGIS), but this is giving me issues when using spatial queries. GLength returns a non-projected length (3.16km) rather than the actual (Ellipsoidal) length (2km). I'm assuming I'm going to have similar issues with area. QGIS can do this simply by checking the Ellipsoid option for measuring length/area so I know its possible with the data I have. But how can I get the actual values for length and area from this data using spatialite? I'm hoping this is a pretty simple thing to resolve. Thanks in advance for any help you can provide!


title: Free Tide prediction service - World (a freeware program would work as an alternative)

text: Is there a free service where I could input a GPS location and it would return me the current tide situation at the position? As an alternative - some library / program that does the exact thing one could integrate? For my app I am currently using the public-domain, accurate NOAA predictions, but these are limited to the US. :( Several apps provide that data and I am curious what their source is.


title: If/then Field Calculator

text: I have two fields in my model, the first is a spatial join return for the township, the second is a spatical join to a city. If there is not city, then the field is left blank. The column that has the city name, I wanted to use a if/then statment to fill in the blank columns with the township name, does anyone have any idea how I would write this or if this is even possible?


title: Iterating through z-values in ModelBuilder

text: I have a series of fields that I am trying to run the IDW tool against. I'm looking for a way iterate through the feature class, using a different field as the z value each time. This is my data: http://imageshack.com/a/img849/115/zqkc.png I'm looking to use each field following Village as a z value. Typically I would just do this all manually but I have a ton a data that I need to run the IDW for so I need to automate it. I've looked through the forums already and I haven't been able to find any responses that include all the necessary details for me to make this happen. Any insight would be incredibly helpful.


title: How to automatically add identical new polygons to a layer?

text: I am seeking a way, using QGIS, to automatically add polygon features identical to others on the layer. In this case, they are water valves on a water network. Anyone know of a user-friendly way to add new valves that visually line up with the pipes without needing to manually draw the polygons through editor mode?


title: Export to PDF in python with parameters

text: I need to export to PDF with some parameters, and I went through the help and it only explains the data frame portion. When I need the entire map or layout view. Here is what I need: 1.400 dpi resolution 2.Embed all document fonts 3.The maps are 36x48 if that matters.               import arcpy     mxd = arcpy.mapping.MapDocument(r"C:\Project\Project1.mxd")     df = arcpy.mapping.ListDataFrames(mxd, "Entire Map?",)[0] #What is the zero for, frames?     arcpy.mapping.ExportToPDF(mxd, r"C:\Project\Output\Project1.pdf", df)          del mxd


title: Google Maps V3 + multiple Fusion Table layers = inconsistent markers?

text: One of my colleagues is working on a map for a news story and is having some difficulty getting the markers to consistently show up. Changing the order in which the layers instantiate is causing significant variance in which markers are actually visible: Compare: http://jsfiddle.net/aendrew/gTgHm/3/ ...With: http://jsfiddle.net/aendrew/hLVzM/1/ ![Version 1](http://i.stack.imgur.com/8YICH.png) ![Version 2](http://i.stack.imgur.com/RKpQn.png) For an obvious example, look to the right of the label for Glasgow — there are two green markers in the first version, three green in the second. **Why is this?** Is there some query limit I'm hitting with the Fusion Tables? Note that each layer has **less than 100 markers.** Any help would be most appreciated!


title: How to add other column into pgrouting query result in pgAdmin output pane?

text: I'm working on my final research project now and I'm a new user in pgrouting. Is there any chances to add another column name to display in the output pane? I'd like to add road name column so that I can easily know the road, for the final result I'll build a web app![here the screenshot of my pgAdmin window query result using A*](http://i.stack.imgur.com/Xw9Mv.jpg)


title: How to Print from free version of CartoDB?

text: I am a complete novice to CartoDB and slowly getting to grips with the free version. I am trying to export the maps into a MS Word document and cannot find a way of doing it. Currently I am having to screen print to Photoshop and then cut and paste into a document. The problem is the narrow viewable screen area limits what can be seen at once and the resolution of this method is poor. Is there another way? As I said I am a complete novice so if you can offer me advice I would be most grateful but can you do it in simple terms please :-)


title: Geodjango render_to_kml Multipolygons - Can be fast or faster as WMS?

text: I am learning now geodjango. I made basic web app with postgis database. With openlayers iam displaying basic wms map and I want to add other layers from my postgis database. I use render_to_kml function and then I visualize this kml with openlayers. The problem is that I can only display point layers instead of multipolygons...I tried to make another kml template bud I misunderstand it and there is no good documentation... I paste here my files ( view, template, urls) and pls give me advice if you can,,,where I have mistake ( I think it will be something with KML templates but i dont know.) ![this is my views.py](http://i.stack.imgur.com/7TY46.png) ![this is my html template](http://i.stack.imgur.com/WNIzO.png) ![this is my urls.py](http://i.stack.imgur.com/LNeFX.png)   1. this is my views.py   2. this is my html template   3. this is my urls.py This results into application with 3 layers (basic map and 2layers from my postgis database which are equivallent (both are point feature....Polygons feature is only copy of Point feature...Thats my problem))... Anyone have pls solution for me? Thank, Mchal Beniak


title: Getting error 'exceptions.NameError' in Python script?

text: I'm getting `<type 'exceptions.NameError'>: name 'B' is not defined` when I run this script. Why is my 'B' not defined and is it a syntax error? For clarification, I am trying to get the code to loop through specific shapefiles in the folder "\\My Documents\ArcGIS\..." and copy them into "\\My Documents\NewFolder\..." The A variable represents a list of the specific shapefiles I want to copy from a folder and into a new one; the issue I am having is defining B, where the individual variable name from A is automatically carried over/read in any part of the code. I am trying to tell Arc that every time it sees "+STATE+" (which should be seen as B: main issue) it must insert the specific variable name listed in A, process it (in this case copy), and then loop back to the next specified variable, until it reaches the end of the list. I would like to copy the 'FL_RechargeArea' _through_ 'GA_AreaofContribution' shapefiles from the "\\My Documents\ArcGIS\..." folder into the "\\My Documents\NewFolder\..." and get it to save as '(specific-varialbe- name)2.shp'. It's much less about copying and more about being able to loop through and process specific shapefile names in a folder, rather than looping through and processing perhaps the entire list within a folder.               # Import arcpy module     import arcpy          for A in B:              # Local variables:         A = ['FL_RechargeArea','FL_AquiferProtectionArea','FL_AreaofContribution','GA_RechargeArea','GA_AquiferProtectionArea','GA_AreaofContribution']         B = "STATE"         New_Shapefile = "\\\\My Documents\\ArcGIS\\"+STATE+".shp"         Output_from_New_Shapefile = "\\\\My Documents\\NewFolder\\"+STATE+"2.shp"              # Process: Copy Features         arcpy.CopyFeatures_management(New_Shapefile, Output_from_New_Shapefile, "", "0", "0", "0")


title: Field Calculator syntax

text: I can't get the syntax correct for a fairly simple if statement. Right now it looks like this:               dim Output as Integer     If [HECTARES] >= 30 AND [HEIGHT] >= 17 AND [CROWN] >= 59      THEN Output = 1      ELSE Output = 0     END IF      I have no idea whether that's the correct use of "AND" but the expression wouldn't evaluate with only one expression (i.e. [HECTARES] > 30) either. Thoughts? (and thanks) Phil


title: ArcGIS 10.1 License Issue via Parallels 9 on MacBook Pro

text: I've been using ArcGIS 10.1 on my MacBook Pro (2011 model) via Boot Camp for a few months. A few days ago, I just installed Parallels 9 in order to use ArcGIS and my Mac programs at the same time and I've run into an issue with the license needing to be repaired. When I try to open ArcMap via Parallels, I get an error message saying something to the effect of "the license is in need of repair". When I installed Parallels 9, I chose to import my Boot Camp so I didn't have to reinstall anything. My guess is that ArcGIS thinks Parallels is a different computer and that's why the license is failing? Any ideas on how to resolve this issue? I can get a hold of another ArcGIS license so if it's a matter of using another license, I should be ok. Thanks.


title: error message when clipping raster with vector layer

text: When I try to crop a raster using QGIS Clipper tool, I get the following error message: `ERROR 1: Attempt to create 0x0 dataset is illegal,sizes must be larger than zero.`. Why? The code QGIS generated to perform the clip was as follows: `$ gdalwarp -q -cutline <clip_mask>.shp -crop_to_cutline -of GTiff <input_raster>.tif <output_raster>.tif` The same error message is generated if I try to run the gdalwarp command by itself from the commandline.


title: Unit of $length in QGis Field Calculator

text: I'm calculating field of various features in my shapefile using QGis. I'm getting some sample length as `0.134914579766986`. I want to know in what unit (eg: meters, km etc) is the above value ?


title: OSM Address suggestion service

text: For geocoding we are using the Nominatim from the below URL http://nominatim.openstreetmap.org/ But i need to implement suggestion on the address, do OSM has this service like Google and ESRI is their any other method to implement?


title: how to extract a rectangle corner coordinates in python script

text: i want to extract a rectangle corner coordinates of a shape file in python script..please suggest me a solution. Thanks in advance, Bhavya


title: OSM DataImported in to Postgres: What is the significance of "_gen0" and "_gen1" suffixed tables

text: Imported .OSM file to postgres/postgis (IMPOSM). Along with the table "osm_new_motorways" I see two more tables "osm_new_motorways_gen0" and "osm_new_motorways_gen1". Same is there for few other tables also. What is its significance?


title: Using ArcGIS in mixed version environment?

text: What are some best practices for ArcGIS 9.3 and 10.0 in the same shop? We have a couple of dozen users on 9.3. The bulk of our data is in file geodatabases and geotiffs, with a respectable splash of shapefiles and occasional jaunts to ArcSDE (9.3.1) hosted by another department. We will be introducing 10.0 to the power users first and for an unknown period of time, likely 6 months, the two will be operative at the same time. All of the data as well as map compositions, layer files, symbol styles, etc. will be used by both groups. How do you structure things so these two sets can peacefully co-exist? What are the things to watch for?


title: Is there a way to code-complete arcpy scripts using an external editor?

text: I have been scripting in Python for ArcGIS 10, and so far I have just been using Notepad++. I like how the Python command line interface code-completes arcpy functionality for you, showing you your options, and I was wondering if there is a way to add this feature to an external editor.


title: Creating multiple feature classes from multiple text files using python

text: I have a folder with multiple text files. I am writing a code that reads all of the text files, and creates a point feature class (i.e. a shape file) for each one. The text files are all tab delineated with the same set up "Id, name, type, lat, long" I have posted the code I have below. A shape file is created for each text file, and the fields are created for the attribute table, however, the points are not created, and the rows are not being inserted into the attribute table. I had it working for an individual text file but I would like to do it for multiple text files. Any help would be appreciated. # importing arcpy and os import arcpy, os from arcpy import env arcpy.env.overwriteOutput = True Location = r"C:\temp" # Set environment settings env.workspace = Location for file in os.listdir(Location):               if file.endswith(".txt"):         try:                  fcName = os.path.splitext(file)                 outName = fcName[0]  + '.shp'                    outFolder = r'c:/temp/Data'                           #Creating the feature class             datasetCreate = arcpy.CreateFeatureclass_management(outFolder, outName, "Point", "", "", "")                  cursor = arcpy.InsertCursor(datasetCreate)                                #Creating the fields for the attribute table             arcpy.AddField_management(datasetCreate, "Name", "TEXT", "", "", 40)             arcpy.AddField_management(datasetCreate, "Type", "TEXT", "", "", 20)             arcpy.AddField_management(datasetCreate, "Lat", "DOUBLE", "", "", 20)             arcpy.AddField_management(datasetCreate, "Long", "DOUBLE", "", "", 20)                                    #Open the text file and read it              textfile = open(file, "r")                              #Create a variable for the headerline             headerline = textfile.readline()                                   #Get rid of the headerline, and split the text file into columns             valueList = headerline.strip().split("\t")             print valueList             Lat = valueList.index("Lat")             Long = valueList.index("Long")             Name = valueList.index("Name")             Type = valueList.index("Type")             ID = valueList.index("ID")                              #Insert a cursor using arcpy                              for point in textfile.readlines():                 segmentedPoint = point.split("\t")                 latvalue = segmentedPoint[Lat]                 longvalue = segmentedPoint[Long]                 namevalue = segmentedPoint[Name]                 typevalue = segmentedPoint[Type]                 idvalue = segmentedPoint[ID]                              #Create a point object and set the lat and long fields to the x and y coordinates                     vertex = arcpy.CreateObject("Point")                 vertex.X = latvalue                 vertex.Y = longvalue                                        #Use the cursor to insert each new row for the attribute table                   feature = cursor.newRow()                       feature.shape = vertex                 feature.setValue("ID", idvalue)                 feature.setValue("Name", namevalue)                 feature.setValue("Type", typevalue)                 feature.setValue("Lat", latvalue)                 feature.setValue("Long", longvalue)                      cursor.insertRow(feature)                                  # If an error occurred print the message to the screen              except:                 # If an error occurred print the message to the screen             print arcpy.GetMessages()                                          #delete the cursor             del cursor


title: GeoDatabase user access based on administration boundry

text: I have silverlight application consuming Map service. Now I want to create user login. The scenarion is somethinkg like when a user logs he should only be able to query on data belongs to his area. Is there any possibility to achieve it using single service?


title: What are the steps to go from raw address data to a locator report?

text: I'm a seasoned programmer, but I'm new to GIS in general. I've read a decent amount of ArcPY documentation but haven't seen an example that takes me from raw data to a locator report, from start to finish. I have all my addresses in a flat file and I want to programatically load them into a DBF, after cross checking their validity, I think? Then I need to add some feature classes? Then I can run a locator report? All of this through ArcPY. Is that thinking correct and is there a complete example anywhere?


title: Creating StreetView/Bing Bird's Eye Add-In?

text: I'm trying to create an Add-In for this tool: Link to Google Maps and Bing from MXD (Streetview and Bird's Eye supported) I'm not receiving an error however my browser is not opening after I select either the Yes or No option in the messagebox. Here is the code. I'm using Windows7 and IE8.               Public Class StreetViewBirdsEye         Inherits ESRI.ArcGIS.Desktop.AddIns.Tool              Const SW_SHOWMAXIMIZED = 3         Const SW_SHOWMINIMIZED = 2         Const SW_SHOWDEFAULT = 10         Const SW_SHOWMINNOACTIVE = 7         Const SW_SHOWNORMAL = 1              Private Declare Function ShellExecute Lib "shell32.dll" _     Alias "ShellExecuteA" (ByVal hWnd As Long, ByVal lpOperation As String, ByVal lpFile As String, ByVal lpParameters As String, ByVal lpDirectory As String, ByVal nShowCmd As Long) As Long              Public Sub New()              End Sub              Protected Overrides Sub OnUpdate()             Enabled = My.ArcMap.Application IsNot Nothing         End Sub              Protected Overrides Sub OnMouseDown(ByVal arg As ESRI.ArcGIS.Desktop.AddIns.Tool.MouseEventArgs)                       MyBase.OnMouseDown(arg)                  Try                      My.ArcMap.Application.CurrentTool = Nothing                      Dim pMxDoc As IMxDocument                 Dim pMap As IMap                 Dim pPoint As IPoint                 Dim pSpatialReferenceFactory As ISpatialReferenceFactory                 Dim pSpatialReference As ISpatialReference                      Dim pApp As IMxApplication                      pApp = CType(Hook, IApplication)                 pMxDoc = CType(pApp.Document, IMxDocument)                      pMap = pMxDoc.FocusMap                 pApp = My.ArcMap.Application                      pSpatialReferenceFactory = New SpatialReferenceEnvironment                 pSpatialReference = pSpatialReferenceFactory. _                     CreateGeographicCoordinateSystem(esriSRGeoCSType.esriSRGeoCS_WGS1984)                      ' convert mouse click to map units        pPoint = pApp.Display.DisplayTransformation.ToMapPoint(arg.X, arg.Y)                 pPoint.SpatialReference = pMap.SpatialReference                 If pPoint.SpatialReference.Name <> pSpatialReference.Name Then                     pPoint.Project(pSpatialReference)                      End If                      Dim URLstr As String = Nothing                 Dim returnLong As Long                      'more info on google map URL query string request setting parameters:                 'http://mapki.com/index.php?title=Google_Map_Parameters                      Dim webmapchoice                      webmapchoice = MsgBox("View this location in a browser?" & vbCrLf & vbCrLf _                 & "Click YES for BING Bird's Eye" _                 & vbCrLf & "Click NO for GOOGLE Street View" _                 & vbCrLf & vbCrLf & "or CANCEL to exit.", vbYesNoCancel, "Web Mapping Analysis")                      Select Case webmapchoice                          Case 6                         URLstr = "http://www.bing.com/maps/default.aspx?v=2&cp=" & _                         pPoint.Y & "~" & pPoint.X & "&lvl=/8&style=o&dir=0"                         ' goto http://www.bing.com/community/blogs/maps/archive/2008/04/10/live-search-maps-api.aspx                         ' for additional parameters to enter for Bing Maps                          Case 7                         URLstr = "http://www.google.com/maps? ie=UTF8&layer=c&cbll=" & _                         pPoint.Y & "," & pPoint.X & "&cbp=11,0,,0,5&ll=" & pPoint.Y _                         & "," & pPoint.X & "&z=16"                          Case 2                         Exit Sub                      End Select                      'Use one of the constants as the window state parameter                 returnLong = OpenLocation(URLstr, SW_SHOWNORMAL)                  Catch ex As Exception                 MessageBox.Show("Caught an unspecified error in the calling code: " & vbCrLf & ex.ToString)             End Try              End Sub              Private Function OpenLocation(ByVal URL As String, ByVal WinState As Long) As Long                  Try                 'PURPOSE: Opens default browser to display URL                      'RETURNS: module handle to executed application or                 'Error Code ( < 32) if there is an error                      'can also be used to open any document associated with                 'an application on the system (e.g., passing the name                 'of a file with a .doc extension will open that file in Word)                      Dim lHWnd As Long                 Dim lAns As Long                      lAns = ShellExecute(lHWnd, "open", URL, vbNullString, _                 vbNullString, WinState)                      OpenLocation = lAns                      'ALTERNATIVE: if not interested in module handle or error                 'code change return value to boolean; then the above line                 'becomes:                      'OpenLocation = (lAns > 32)                  Catch ex As Exception                 MessageBox.Show("Caught an unspecified error in the calling code: " & vbCrLf & ex.ToString)             End Try              End Function          End Class      Maybe this is not coded to open IE8? Any ideas of what needs updated? Thank you!


title: FOSS method to split raster into smaller chunks

text: I have a raster (USGS DEM actually) and I need to split it up into smaller chunks like the image below shows. That was accomplished in ArcGIS 10.0 using the Split Raster tool. I would like a FOSS method to do this. I've looked at GDAL, thinking surely it would do it (somehow with gdal_translate), but can't find anything. Ultimately, I'd like to be able to take the raster and say how large (4KM by 4KM chunks) I would like it split up into. ![enter image description here](http://i.stack.imgur.com/AKPJ8.png)


title: How to setup PyCharm to reference ArcPy?

text: How would I setup PyCharm to recognize import ArcPy? I'd like to use PyCharm to develop and debug geoprocessing tool.


title: How to Clip a Shapefile of Circle features using a Raster using ArcGIS 10

text: I want to keep all the features in the shapefile that lie in the raster region and reject others. How do i do this? ------Adding More Details-------- Data to be Clipped: Is A Shapefile of 1 mile circles. Data I want it clipped with : is a Raster (GeoTiff or ADF). Now I want to find out all the Circles from shapefile that overlap on the Raster. I dont want to modify the raster here. I want those features (circles) from the Shapefile removed which donot overlap with the raster. * * * Problem is I have tried the tool mentioned by Dano, but tools expects a feature input in the second field, and not a raster directly. Thanks Shaunak


title: Why does "Add Database Server" in Catalog 10 allow only ArcSDE Personal or Workgroup databases?

text: Arc Catalog 10 - ArcInfo - Database Severs - Add Database Server: "Add an ArcSDE Personal or Workgroup database server by typing in the SQL Server Express instance you wish to connect to." I want to connect to an MS SQL Server 2008 R2 Developer (like Enterprise) server on another machine, but it sounds like I cannot do that from here. If I try it anyway I get the message, "Connections to this edition of the database server are not supported." My keycodes file for licensing has "svrenterprise" in it, so I don't know why I cannot connect to this database server. I want to make a direct connection to the database server, through the open port 1433, not 5151 SDE port.


title: WFS Filtered : The filters that worked for me in 1.8 - gives no results in 2.0

text: Using Quantum GIS 1.8.0 and 2.0.1 desktop on windows XP 32bit Server: http://vmgeoserver.vd.dk/geoserver/VD/ows?version=1.1.0 Setting filter for VD:vma_admdata:               "BESTYRER"  = '230' AND  "KODE_ADMIN_VINTERVEJSERVI"  = '230_V1' AND  "KODE_ADMIN_VEDLIGEHOLD_NI" = '230_N'      In 1.8 this returns the results I'm looking for. Metadata for the resulting layer states _Source for this layer:_               http://vmgeoserver.vd.dk/geoserver/VD/ows?version=1.1.0&SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=VD:vma_admdata&SRSNAME=EPSG:25832&FILTER= BESTYRER 230 KODE_ADMIN_VINTERVEJSERVI 230_V1 KODE_ADMIN_VEDLIGEHOLD_NI 230_N       Doing the same in Quantum GIS 2.0.1 gives no error messages but the resulting layer is empty. Checking metadata the _Source for this layer_ is:               http://vmgeoserver.vd.dk/geoserver/VD/ows?version=1.1.0&SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=VD:vma_admdata&SRSNAME=EPSG:25832&FILTER= BESTYRER 230 KODE_ADMIN_VINTERVEJSERVI 230_V1 KODE_ADMIN_VEDLIGEHOLD_NI 230_N       ie. identical. What am I doing wrong?


title: Distributed WMS (and TMS, optionally) server architecture, GPU rendering

text: What if we could use the power of dozens of servers to process a lot of incoming WMS requests? And more, using the power of GPU's to do the job of styling the maps for scalability and speed! The theme of my senior research project is exactly this, distributing the WMS service (and probably, TMS) to reach a maximum horizontal processing of GIS data. I'm very excited about this, and I would like to hear from you, the professionals in this field, ideas and criticisms regarding this approach.


title: Help with null values in rasters

text: I need to add different boolean rasters of the distribution of some species, since we need to obtain a single file that shows the locations where there are more than one species occurring. To do that, I rasterized all of the shape files showing the distribution and then reclassified them to obtain the boolean rasters. The problem is, now that I'm trying to add them all together using rastercalc, there are files with different extends even though I clipped them with the same shape file (the final single file -the one that I obtain by adding all of the booleans, is all cut and incomplete with null values)... how can I change them to the same extent in order to obtain the complete map... ? the first picture is the shape file I used to clip the other shape files (this is the polygon of the state I live in). The second picture shows an example of one boolean images with the raster of the shape I clipped them (the state's shape) The third one shows the last raster with all of the boolean rasters together... if you compare images 2 and 3, you will see there is lots of missing data ![enter image description here](http://i.stack.imgur.com/aNZrY.png) ![enter image description here](http://i.stack.imgur.com/eYUsU.png) ![enter image description here](http://i.stack.imgur.com/gxHVt.png)


title: Postgis 2.0 - How to get a Point from given LinestringM and a measure?

text: If I have a cable route linear referenced and stored in a LinestringM with M value is the length (in meters) from the start of the route - for example:               LINESTRING M (0 0 0, 10 0 20, 12 0 40, 20 0 50, 21 0 70)      Now I know the cable route is broken at _45 meters from start point_ , that would be in the middle of two points `(12 0 40)` and `(20 0 50)` (see the M value). With Postgis, how can I get the broken point coordinate - which should be `POINT (16 0)` in this case?


title: How do i draw a route on Quantum GIS between nodes?

text: I have a requirement to trace path on the points present in the shapefile which is loaded in Quantum GIS. This should enable the user to specify their own routes within the present network. **There will be existing points and lines present in the shape file. Whenever the user clicks on the points there should be a line drawn between the points selected and the route should be saved as a shape file.** **I have searched and found in many places of options like pgRouting and shortest path routing and that is not what i require.My need is simple and due to these options being available i m not able to find a simple solution.Its very much similar to this question :How to route railcars along railroads?. which does not yeild a clear solution.** I am a Newbie in GIS. Any guidance will be very much appreciated. Thanks in advance.


title: How to set field values as labels in QGIS2 Print Composer

text: I am mapping the distribution of 325 bird species, displaying circles in the grid cells each species occurs in. The data contains as many duplicates of the NAME as there are cells in which it occurs. I produce a map by filtering for a NAME then want to print the maps using the NAME field in a label as a caption below the map. I would like to have it generated automatically rather than write it in manually - in anticipation of being able to generate all species maps as an Atlas (that's another issue). So far, all I get by using `[% name %]` is `[name]`, not the species name. Any ideas of what I might be missing? Nothing I've found by searching has helped so far. Many thanks Niven


title: ArcGIS Android MapView with custom Spatial Reference SVY21

text: I have a MapView that is using a custom spatial reference of SVY21 from the server http://e1.onemap.sg/arcgis/rest/services/SM128/MapServer. The problem is the map's location service's accuracy circle due to the SpatialReference.getUnit causing an exception shown below               Could not update the active symbol     java.lang.RuntimeException: Invalid SpatialReference.          at com.esri.core.geometry.SpatialReference.nativeGetUnit(Native Method)     at com.esri.core.geometry.SpatialReference.getUnit(Unknown Source)     at com.esri.android.map.LocationService$a.k(Unknown Source)     at com.esri.android.map.LocationService$a$a.onLocationChanged(Unknown Source)      I really need the accuracy circle. Is there any solution to this?


title: Is it possible show the feature attributes dialog when I add a feature in code?

text: I want to know if it is possible to make the feature attributes dialog in Qgis show up when I add a feature by code using: `QgsVectorLayer::addFeature( tFeature )` This is available when I add features using Qgis' digitizing tools but not when I add features by code!


title: Image classification for a study site falling over two adjacent Landsat 8 scenes (qGIS or ArcGIS)

text: This question concerns work-flow rather than specific methodological problems. I will be using GRASS to classify Landsat 8 images to create a vegetation map for a study site that happens to fall over 2 separate Landsat scenes. I would really appreciate some advice about which route to take and I'm not able to find any other information about this elsewhere: 1.) Should I **FIRST** find a way to combine the 2 scenes together (which presumably would require combining the 2 scenes of each channel separately - would that be right?), and then cut out the small area covered by my study site and THEN start the classification process. OR 2.) Should I conduct image classification on the 2 adjacent Landsat scenes separately and then combine right at the end? (Although I'm worried this will take much longer, not to mention use up more memory - the study site I am working with is only 1400 km2) N.B. I am using a Mac (Mavericks OSX) and have access to both QGIS 1.8.0 and ArcView 9.2 to conduct any manipulation on this Geotiff files. Thank you very much for your time and advice!


title: How to use Plus tool in ModelBuilder?

text: I want to have the **sum** of at least 2000 raster combined to one raster layer using ModelBuilder ( **Plus tool** ). I think I'm missing something here. This model does not work. ![sample model builder flowchart](http://i.stack.imgur.com/ig9gd.jpg)


title: Memory layer or QgsVectorFileWriter Error?

text: Currently this excerpt of code will not save to file (as it is supposed to). The memory layer is visable and editable. Yet even when manually saving the layer, any priorly added features (geometries) are non-present. The variables 'FileDirectory._destination' and 'GetMaps._TOPOcrs' are proper as confirmed via the print command. It should be noted, that the error prompt for the FileWriter indeed reads an error with the value '1'. Any clues?                   try:             self.estabLocal()         except ValueError:             iface.messageBar().pushMessage("Error", "Please select a central folder location to save your model to!", level=QgsMessageBar.CRITICAL, duration=3)             return         try:             self.estabCRS()         except ValueError:             iface.messageBar().pushMessage("Error", "Please select a topographic map for your simulation!", level=QgsMessageBar.CRITICAL, duration=3)             return         pointstress = QgsVectorLayer("Point?crs=" + GetMaps._TOPOcrs + "&field=Well:string(3)&field=Discharge:string(6)&index=yes", "Points" , "memory")         pr = pointstress.dataProvider()         QgsMapLayerRegistry.instance().addMapLayers([pointstress])         PointStressFeat = FindDirectory._destination+"Points.shp"         writer = QgsVectorFileWriter(PointStressFeat, "CP1250", pr.fields(), QGis.WKBPoint, pr.crs(), "ESRI Shapfile" )         if writer.hasError() != QgsVectorFileWriter.NoError:             print "Error when creating shapefile: ", writer.hasError()         iface.setActiveLayer(pointstress)         pointstress.startEditing()         iface.actionAddFeature().trigger()


title: Why do I lose my printer settings when saving a copy of mxd from ArcPy?

text: Judging from the answer to this question, one cannot change the page size etc from ArcPy, so I set up my "template" mxds so my script can save a new copy of it, change data sources etc. This methods works perfectly for my A4 maps. However, my A0 mxds are losing the printer and page settings from the template. My `A0_template.mxd` is set to print from the plotter on A0 paper using the default printer settings. After the copy of the mxd is produced and I open it, the printer has somehow changed to the normal small printer with the paper size matching the default tray. I am simply using the `saveACopy(filename)` method which is part of the MapDocument class. Why would the correct settings for the A4 mxd be saved, but not for the A0?


title: Obtain address from PostGIS OSM data like Nominatim does?

text: We currently use the Nominatim service to do address lookups. We want to store the OSM data in our own servers as our usage gets higher. How do I calculate the address of a given coordinate, like how Nominatim does?


title: How can I split an osm file into separated files with poly files using osmosis?

text: I have got a country.osm file and I want to split it alongside of admin_level=8. I extracted previously ~3.000 poly files into a separated folder. I have batch script, that split them, if I give the input and the poly:               osmosis --read-xml %1 enableDateParsing=no --bounding-polygon file=%2 completeWays=no completeRelations=no clipIncompleteEntities=true --write-xml file=clc_%2.osm      **The question:** Are there any way to do this parallel and the input file read just once?


title: Ways to Speed Up Python Scripts Running As ArcGIS Tools

text: This is a pretty general question. I just wondering what tips and tricks GIS programmers have used to speed up arcpy scripts that you import into the toolbox and run. I work most everyday writing little scripts to help non-GIS users at my office process GIS data. I have found that ArcGIS 10.0 processing in general is slower than 9.3.1 and sometimes it gets even slower when running a python script. I'm going to list a particular example of a script that takes over 24 hours to run. It's a loop that tabulates the area of a raster in a buffer for each shape in the buffer. The buffer has about 7000 shapes. I don't believe it should run this long. A               while x <= layerRecords:              arcpy.SetProgressorLabel("Tabulating Row: " + str(x) + " of " + str(ELClayerRecords))         arcpy.SelectLayerByAttribute_management(Buff,"NEW_SELECTION", "Recno = " + str(x))                                  # Selecting the record         TabulateArea(Buff, "Recno", MatGRID, "VALUE", ScratchWS + "/tab" + str(z) +".dbf", nMatGRIDc)                          # Tabulate the area of the single row              arcpy.AddMessage ("          - Row: " + str(x) + " completed")         x = x + 1         z = z + 1      Before anyone says it, I have run tabulate area on the entire buffer, but it produces errors if run on more then 1 record. It's a flawed tool, but I have to use it. Anyways, if anyone has any ideas on how to optimize, or speed up this script, it would be most appreciated. Otherwise, if you have any speed up tricks for python, when used in ArcGIS, I'd love to hear them. Thanks in advance for any replies!


title: How to test whether two countries are adjacent?

text: I have an interesting question:   Is there a dataset by which I can use to test whether two countries are adjacent? For example, if inputs are Mexico and U.S., they are adjacent; if inputs are England and Russia, they are not. **Is it easy to produce an adjacency matrix for all the countries?** **Edited:** This question might be extended as, if we get a political boundary map, is there any ways to produce an adjacency matrix for all the political entities (e.g. Given a map of Florida, can a counties matrix be built to indicate the adjacency of all the counties)


title: How to draw curves in Flex Viewer?

text: I want to draw curves from a given point to another in Sample Flex Viewer like shown below. Any help? ![Routes](http://i.stack.imgur.com/Nk20s.png) ![enter image description here](http://i.stack.imgur.com/vt6Gm.png)


title: How to fix symbols/labels in a tiled WMS?

text: I have a postgis stored dataset that is consumed in a qgis project and served on qgis server. I noticed however that symbols and labels of point features that overlap 2 tiles of the map are only rendered in the tile that actually contains the feature. Is there some kind of buffering option or a configuration option for rendering to fix this? ![enter image description here](http://i.stack.imgur.com/GLGsn.png)


title: Connect JSON to GeoExt Application

text: I am returning a JSON to a GeoExt app and this is the code that is collecting the return,               var control = new OpenLayers.Control.Click({                          trigger: function(evt) {                           var lonLat = evt.xy;                           if (popup) popup.close();                           var  request = Ext.Ajax.request({                           url:"http://128.196.142.94:9000/info?lat="+ panel.map.getLonLatFromViewPortPx(lonLat).lat +"&lon=" + panel.map.getLonLatFromViewPortPx(lonLat).lon,                           disableCaching: false,                           success: function(response) {                                   var obj = Ext.decode(response.responseText);                                   console.log(obj);                                   console.dir(obj);                                   addToPopup(lonLat, response.responseText);                               },                           failure: function(response) {                                   console.log('server-side failure with status code ' + response.status);                                   addToPopup(lonLat, "Server side Failure");                               }                                    //callback: myhandler                      });     }});      I wish to obtain the returned JSON and then assign it to the popUp but whenever I try to connect to the App it is returning the following error,               OPTIONS http://128.196.142.94:9000/info?lat=55.44873046875&lon=-118.17626953125 405 (Method Not Allowed)     XMLHttpRequest cannot load http://128.196.142.94:9000/info?lat=55.44873046875&lon=-118.17626953125. Origin null is not allowed by Access-Control-Allow-Origin.     test.html:222server-side failure with status code 0      If you simply pass the url `http://128.196.142.94:9000/info?lat=55.44873046875&lon=-118.17626953125` you will obtain the JSON but for some reason there is failure that keeps occurring in the application. Is there CORS that can be used in the OpenLayers or GeoExt. I have tried using ExtJS JSON-P but I am not sure what should I enter in the 'callback' of,               <script src="http://domainB.com/users?callback=someCallback"></script>


title: Network Analyst directions window options

text: I have created a vehicle route problem with multiple route zones with one depot in each zone. Each zone has a pickup route that delivers to the depot in the zone. Some addresses have special instructions for where to pickup the item from their property as well as a phone number. My question is whether or not there is a way to make any special instructions or any other information like a phone number appear in the directions window to be printed out for the driver? I'm using Network Analyst in ArcGIS 9.3.1.


title: How to use ArcView GIS extension in QGIS or how to convert one to do so?

text: I really would like to know whether is a way to can work with ArcView GIS extensions at QGIS? or at least a way to be converted to plugins at QGIS?


title: Force a new AddDataDialog to open at a defined directory path

text: Is there a way to force the AddDataDialog to open to a specific path? I created a tool that when clicked on a specific feature will open the Add Data dialog, but I need it to open starting at a specified location. AddDataDialog doesn't have a property analogous to OpenFileDialog.InitialDirectory.


title: Editing with ArcGIS Silveright

text: I am using the ArcGIS Silverlight API to edit data in an ArcSDE SQL Server geodatabase. I have been told that any editable feature classes must be registered as versioned in ArcCatalog before they can be edited via the Silverlight API. Does anyone know why this is the case? From what I understand, you can not use the Silverlight API to edit versioned data in a feature class, so I'm not clear on why you would have to register the feature classes as versioned. Thanks, Seth


title: Specify a starting location with IgxDialog

text: I would like to know what is the easiest way to set the starting location of a data dialog to a IWorkspace? My problem is that basically I do not want a connection (.sde file) for the users. I want to set the starting location to a workspace connected via IPropertySet. Is it possible? This is related: Force a new AddDataDialog to open at a defined directory path but does not help. EDIT: following Kirk's suggestion, here's a code snippet. This does not work.                   public static IDataset OpenDatasetDialog(IWorkspace workspace,int handle)         {             var workspaceName = (IWorkspaceName)((IDataset)workspace).FullName);             var gxDatabase = new GxdatabaseClass();             gxDatabase.WorkspaceName = workspaceName;                  // we need to set a different variable as .NET runtime does not allow ref parameters with casts             object gxObj = (object)gxDatabase;                  IGxDialog openDialog = new GxDialogClass();             openDialog.RememberLocation = false;             openDialog.set_StartingLocation(ref gxObj);             var filterColletion = (IGxObjectFilterCollection)openDialog;             filterColletion.AddFilter(new GxFilterTablesAndFeatureClassesClass(), true);             openDialog.AllowMultiSelect = false;             IEnumGxObject datasets = null;             IDataset dataset = null;                  if (openDialog.DoModalOpen(0, out datasets) && datasets != null)             {                 IGxObject gxObject = datasets.Next();                 if (gxObject is IGxDataset)                 {                     IGxDataset gxDataset = (IGxDataset)gxObject;                     dataset = gxDataset.Dataset;                 }                 else                     throw new ArgumentException();             }             openDialog.InternalCatalog.Close();             return dataset;         }


title: Flex popup from outside resource weather warnings

text: I am using the USGS service of NHSS NOAA weather warnings in an application. I have the service displaying in my map using   http://rmgsc.cr.usgs.gov/ArcGIS/rest/services/nhss_weat/MapServer   The service shows that identify is supported.   The question is how do I get a popup from this service with the warning type and other info?   If I try to use the same service I get no popup, If I point it at my own server I get a popup with no information.   I have also tried using the atom described here but not sure of the exact usage in this context. Any help would be appreciated.


title: How to put GeoServer layers behind authentication?

text: I want to publish map layers using OpenLayers and GeoServer. While this is easy enough to set up, I'm a bit lost when it comes to authentication possibilities. What I want is to create different users that will, by logging in with their own user name and password, be directed to their own map window. How can this be securely done so that none of my users will be able to access each other's layers? I have understood that it is quite easy to "hack" the setup by, for example, making queries to the GeoServer where one simply asks it to list all the available layers. What are my options? I suppose apache could be of help when setting it up in my ubuntu server? Any tutorials on this subject?


title: Rotate while in Layout View

text: I am using ArcGIS 9.3 without any extensions. How do you get your map to rotate? I've done it before but it has been awhile and now the "rotate" does not seem to function now. I've tried using while in edit and non-edit. Any suggestions?


title: How to find out Why a QGIS plugin is rejected for Version 2

text: I have a selection of python plugins (one rather large one) that I have developed under Qgis 1.8 and that I am migrating to qgis 2.0. The error that the Plugin Manager gives me when I load these plugins (This plugin was developed under 1.x) is rather unhelpful. It also arises in situations that have nothing to do with the API Version. For instance, a Network error yesterday introduced a stream of random characters into a few of my code files (Fun!) and this led to the Plugin Manager telling me my plugin was old. Eventually, after pain, I found out a tab and a dollar sign had been inserted in front of a comment in my metadata.txt, and that was the Problem. Screwed up, yes. Old, no. Is there any way to see a more detailed error when QGIS 2 rejects a plugin as 1.x? If not a proper trace, then at least what class, or even what file was rejected? Anything? The instructions regards API changes at the Qgis Hub Wiki are lovely and I sleep with them under my pillow, but they are not as lovely as a simple error message. Qgis 2, by the way, is beautiful. Thank you.


title: How to create .dem file in QGIS?

text: I have an elevation data in .asc file, but need to convert it to .dem file to later use in COMSOL. However, in Raster--> Analyzes-->DEM I got just error message: ERROR 1: DEM Default Setup: Importing source dataset projection failed What is wrong? regards, L.


title: Possible to run "Merge shapefiles to one" algorithm from Python script?

text: Developing a python script in QGIS 2.2, I can run algorithms from the **Processing** plugin fine such as those from _QGIS geoalgorithms_ and _SAGA_. Can you do the same for the algorithm " _Merge shapefiles to one_ " from the **Vector** menu and how can you call it from the Python console? EDIT: The " _Merge shapefiles to one_ " algorithm is from the Ftools plugin. Found a link which might be useful: Is there a way to access QGIS plugins in Python?


title: Drawing text and line-dasharray

text: Is it possible with TileMills and CartoCSS? ![enter image description here](http://i.stack.imgur.com/WmKWA.jpg)


title: ArcGis Flex Viewer Helloworld widget doesn't display

text: I'm trying to create a new widget, but it doesn't show. If I comment out the other widgets in the widgetcontainer in the config.xml file, they don't go away. I don't understand why the changes I make don't show... config.xml file:               <widgetcontainer layout="float">             <widget label="Bookmarks"         left="430" top="90"                     icon="assets/images/i_bookmark.png"                     config="widgets/Bookmark/BookmarkWidget.xml"                     url="widgets/Bookmark/BookmarkWidget.swf"/>             <widget label="Louisville Police" left="590" top="280"                     icon="assets/images/i_police.png"                     config="widgets/Query/QueryWidget_Louisville_PoliceStations.xml"                     url="widgets/Query/QueryWidget.swf"/>             <widget label="Search" left="80" top="280"                     icon="assets/images/i_search.png"                     config="widgets/Search/SearchWidget_Louisville.xml"                     url="widgets/Search/SearchWidget.swf"/>             <widget label="Earthquakes (GeoRSS)" left="410" top="280"                     icon="assets/images/i_rss.png"                     config="widgets/GeoRSS/GeoRSSWidget.xml"                     url="widgets/GeoRSS/GeoRSSWidget.swf"/>             <widget label="Draw" left="60" top="400"                     icon="assets/images/i_draw2.png"                     config="widgets/Draw/DrawWidget.xml"                     url="widgets/Draw/DrawWidget.swf"/>             <widget label="Print" left="390" top="400"                     icon="assets/images/i_print.png"                     config="widgets/Print/PrintWidget.xml"                     url="widgets/Print/PrintWidget.swf"/>             <widget label="My first widget"                     icon="assets/images/i_widget.png"                     config="widgets/Samples/HelloWorld/HelloWorldWidget.xml"                     url="widgets/Samples/HelloWorld/HelloWorldWidget.swf"/>             <!--             <widget label="My first widget"                     icon="assets/images/i_widget.png"                     config="widgets/Samples/HelloWorld/HelloWorldWidget.xml"                     url="widgets/Samples/HelloWorld/HelloWorldWidget.swf"/>             -->         </widgetcontainer>      I click the "Run" button, the application appears in my browser, but that widget isn't there. I've confirmed that the widget is in that directory, I've added the module. I'm following the directions here: http://video.arcgis.com/watch/539/arcgis-viewer-for-flex-creating-a-custom- widget http://resources.arcgis.com/en/help/flex- viewer/concepts/index.html#/Creating_widgets/01m300000010000000/ ![enter image description here](http://i.stack.imgur.com/GzNX1.png) ![enter image description here](http://i.stack.imgur.com/GI1US.png) ![enter image description here](http://i.stack.imgur.com/yRlVw.png)


title: How do I use wildcards to work around the GDAL character limit?

text: I am attempting to use the 'Merge' & 'Create VRT' tools in QGIS 2.0 on a folder with 3000+ raster tiles in. However whenever I try to perform this operation the following error appears: **'The process failed to start. Either the invoked program is missing, or you may have insufficient permissions to invoke the program.'** The process works when I select only a few tiles, so I can only assume it's due to a character restriction due to the number of files? Upon searching the website for a solution I came across the concept of using 'Wildcards' to get around this problem: Having problems creating VRT in QGIS from png files However I cannot find any clear explanation on how to use them, and I don't have much experience with using OSGeo4W outside of the QGIS GUI. I would appreciate any help or advice.


title: what is Network analysis and how to find shortest path

text: I don't know about Network Analysis.Please any one help to understand the Network Analysis with example. General Process flow and procedure is enough to learn my self Thanks in Advance.


title: I am having Problem with my geoserver using URL

text: I am having problem when I am Requesting following URL in GeoServer /geoserver/wfs?REQUEST=GetFeature&VERSION=1.1.0&typename=Final&propertyname=name I am getting following error http:///geoserver/schemas/ows/1.0.0/owsExceptionReport.xsd"> java.lang.RuntimeException: Unable to obtain connection: Connection refused. Check that the hostname and port are correct and that the postmaster is accepting TCP/IP connections. Unable to obtain connection: Connection refused. Check that the hostname and port are correct and that the postmaster is accepting TCP/IP connections. Connection refused. Check that the hostname and port are correct and that the postmaster is accepting TCP/IP connections. Connection refused: connect Can any one help in Resolving this issue


title: how to make a line from a set of gps points?

text: As a part of my research I need to create the trajectory of respondents according to the their collected GPS data (including X and Y). I am familiar with python programming, but I am a beginner in ArcGIS. I want to write a script to automate this process I really appreciate your help and comments.


title: Where does ArcMap GUI get font size and docking guide positions from?

text: In the screenshot below, ArcMap 10.1 appears to be making some arbitrary choices about the font to use for:   1. The Layers data frame name in the Table of Contents (larger font than expected) - **NOW FIXED**   2. The Main Menu pulldowns (smaller font than expected) - **_NOW RESOLVED_** Can anyone tell me where to find a setting that will enable me to make these the same size as the font used for the Main Menu items, please? ![enter image description here](http://i.stack.imgur.com/e6AIc.png) UPDATES: I finally resolved this by visiting **Control Panel\Appearance and Personalization\Display** and **resetting the font size I had set there from 150% to 125%**. This resolved all. I was also pursuing the mismatch in font size between Menu items and pulldowns via the ArcGIS Discussion Forums. The problem also affected the docking guides used to help resize windows (see graphic below). ![enter image description here](http://i.stack.imgur.com/FmMwX.png)


title: Why does QGIS scalebar show a wrong scale?

text: I'm using QGIS 1.7.1 and have added a scalebar. The scalebar usually displays in degrees, unless I switch the map's projection to meters, in which case the units are display in _cm_. This is not correct. Any thoughts on how I can get this in reasonable units? Thanks!


title: Leaflet Event Craziness: jQuery click events don't propagate to L.DomEvents

text: I noticed that calling `.click()` on jQuery selected DOM elements in Leaflet ( such as anchor icons ) will not do anything -- the click event isn't propagated to the events created with L.DomEvents. For example, go over to bboxfinder and try to trigger the zoom out button's click event:               $('a[title="Zoom out"]').click();      Nothing will happen, the map won't zoom out. I jumped down the rabbit hole to figure out the why(s) and ended up with more questions than answers. My two favorite questions are listed below. I'm hoping the answers for each question address two parts:   1. Why is this behavior happening? Be as technical as you can so people can learn!   2. What are the benefits of this behavior? This is a well-written, abstracted library. I assume they made these choices for reasons that are not easy for someone like me to see. Any guesses what those benefits are? For example, maybe these choices work well with mobile devices. **Questions** **1.** Why are the jQuery events not propagating to DOM elements? I can call click on the unwrapped jQuery DOM element ( see below ) and fire the event, but i assume jQuery would pass this event down for me:               // this works fine, notice we're grabbing the raw element     $('a[title="Zoom out"]')[0].click();      **2.** When i was digging around in the Leaflet source I noticed button events were created/bound with this code workflow below ( it also happens in 3rd party packages like L.Draw ). The question is why are we using all these events on a single button and how are they working together? And really, why is it like this? I've pulled the code below from this part of the source:                   L.DomEvent          .on(link, 'click', stop)          .on(link, 'mousedown', stop)          .on(link, 'dblclick', stop)          .on(link, 'click', L.DomEvent.preventDefault)          .on(link, 'click', fn, context)          .on(link, 'click', this._refocusOnMap, context);


title: select by attribute on shapefiles with open-source python modules

text: I usually use ArcGis 9.3.1, but when running a python script that includes Select_by_attribute (from gp), the script runs very very slowly. Is there a way I can perform a select-by-attribute on a opensource-python-module (such as shapely?)? Will it be faster?


title: Find Nearest features from an input location

text: I wanted to find the nearest features (it can be point, line and polygon layers) from an user input or search location and identify 5 nearest locations along with their distances. Have anybody worked on this link below. If so, can anyone point me to the right sample application, please. https://developers.arcgis.com/javasc...arest-amd.html https://developers.arcgis.com/rest/a...nd-nearest.htm On the above link, there is no sample application or code to show how this works and also it can only query against two feature layers. I wanted to query against multiple layers. Secondly, there should be a way to do without using FindNearest class, as this class has only been introduced with V3.7. If this is the case, then how GIS developers across the world would have created application without this class. I think, some one should really have an answer to this.


title: Changes in esri config in VS2010 has no effect

text: I've created a DockableWindow in my Add-In (vs 2010, c#, arcobjects 10 sp2) and the esri config was updated with the following:               <DockableWindows>         <DockableWindow id="CTArcMapAddin_AreasDockWindow" class="AreasDockWindow+AreasDockWindowImpl" caption="Areas" image="Images\AreasDockWindow.png">           <InitialPlacement height="300" width="250" state="unpinned" position="left" />         </DockableWindow>       </DockableWindows>      I run the Add-In and the dockable window shows up on the left side... Then when changing the "position" value to "right" and running the addin again it's still positioned to the left...why does that change not show through...a bug maybe?


title: Calculate ANOVA using GIS

text: It's possible to calculate a statistic like ANOVA using GIS. I'm using ArcGis, but there's any GIS software to calculate it? Thanks,


title: Run a model from python problem

text: I am using arcGIS 10.0 SP5 I have the following code :               import arcpy          arcpy.ImportToolbox(r"C:\Documents and Settings\Derek\FRAME\FORTEST.tbx", "Derek")          arcpy.TEST_Derek()      The model is names TEST The toolbox alias is Derek When I try and run this I get an error               AttributeError: 'module' object has no attribute 'TEST_Derek"      What I want to do is to run the model TEST from the toolbox FORTEST


title: Defining SQL 2008 store within GeoServer 2.3

text: I want to define a SQL server 2008 store within Geoserver 2.3. I installed the sql Extension. After defining the different attributes, I get the following error message :               Error creating data store, check the parameters. Error message: Unable to obtain connection: Cannot create PoolableConnectionFactory (The TCP/IP connection to the host localhost, port 1433 has failed. Error: "Connection refused: connect. Verify the connection properties, check that an instance of SQL Server is running on the host and accepting TCP/IP connections at the port, and that no firewall is blocking TCP connections to the port.".)      The data I entered : host -localhost/sqlr2 port - 1433 database - the database name within SQL scheme - dbo user & password = A geoserver user i defined within sql Thanks in advance Lauren


title: QGIS Quickly switching overlapping raster images

text: I'm new to QGIS, and haven't used arcGIS in about 5 years, however I've just seen someone using arcGIS and was able to quickly switch between overlapping rasters when checking registration. The top half of the screen was one raster, and the bottom was another, and the mouse was used to move up and down in real time switching between the two (I'm probably not explaining this very well). The images I'm working with are very large (2x1GB geotiffs) and of the same area. There is a good 5second lag between switching layers which makes it difficult to see if there is good alignment. Does this functionality exist in QGIS, to quickly switch rasters in real time?


title: Watershed layers with names for each one for Europe or Sweden

text: I am working on freshwater species and I would like, make a list of each watersheds used for each species. So I am looking for a watershed maps which contain name of the watershed. I found some good layers, but, they don't get names for each basins, only for the bigger ones. I would like to know if someone can give me advices, links or good sources? I need mainly the cover for Sweden at least, or Europe. I searched this site, but I found most of the elements to create watersheds. Regards


title: Is it possible to reproject KML using ESRI's Silverlight API?

text: My question is similar to this one, but for Silverlight. I need to reproject an entire KML file to NAD83, add it to the map with all the layers, symbology, etc. in tact. Is there anyway to do this with the Silverlight API? thanks


title: Does anyone display OpenStreetMap through MapServer, ArcGIS for Server or Geoserver?

text: I use Windows OS. The data (OpenStreetMap) are stored in a postgres/postgis database. How do I display this data correctly? **Update:** Now I use: geoserver-2.3.5.exe, postgresql-9.1.9-1-windows-x64.exe, postgis-pg91x64-setup-2.0.3-2.exe, osm2pgsql **Update2** I lookd at the source for "osminabox". https://git.hsr.ch/index.html#repositoryBrowser;26067c88-f4a8-46e8-be0d-95c1a5091c94;null;Code/osm2gis/geoserver/geoserver_data/styles   1. But I find many sld styles (I have only 4 layers: planet_osm_polygon, planet_osm_line, planet_osm_point, planet_osm_roads).   2. There are a lot of styles for points: osm_poi.sld, osm_building.sld, osm_pow.sld. I don't have this layer.    3. Also I find strings such: "img/recycling.png"  So, If I use intranet system I can't work? **Update 2 continue:** I try to use https://github.com/fotographs/OSM-Styles- for-GeoServer/tree/master/sld. But I cannot find sld styles for the points. https://help.openstreetmap.org/questions/7521/geoserver-style-files - Would it help?


title: QGIS 2.01 - raster transparency not working properly

text: ![Trasnsparent raster image with underlying vector layer](http://i.stack.imgur.com/TlqM4.jpg)I have loaded a color raster image (aerial orthoimagery) into a QGIS project. The layer is visible but very transparent (in the enclosed image the vector polygons are below the imagery), but when I check Properties/Transparency the transparency is set to 0%. The imagery is apparently color infrared. In QGIS 1.8 I simply went to the Color Map dialog and set the R-G-B color map to 4-1-2 respectively and the display was in perfect true color. I cannot figure out how to do the same with 2.01. I have looked at the options under multiband and singleband pseudocolor and neither seems to give me the option to reset the color map as I was able to do with 1.8. You can see that the image is displaying in true color, but highly transparent. I did not change any color map settings. When I view an older project that created with 1.8 with 2.01 imagery from this same data set displays perfectly, so it is not a problem with the source data. However, when I look at the Style settings for the older project with 2.01 I cannot find any differences between that project and the Style settings for the new project. I don't konw if this is a color map issue or a trasparency issue, but I would simply like to get the image ot display properly.