title: Does GeoServer's WorldImage plugin use UTM or WGS84 (EPSG:4326)?

text: I have some PNG images that I want to put on an OpenLayers map using GeoServer as a WMS, and it seems as though the easiest way to do that would be to create an accompanying `.pnw` file (ESRI world file) for each `.png` and load them using GeoServer's WorldImage plugin. I also have the longitude, latitude coordinates of nine points on each rectangular image: {top, middle, bottom} x {left, middle, right}, and I can use five of them to convert pixel x,y into longitude, latitude x',y': ![x' = \\left\(\\frac{mr_{lng} - ml_{lng}}{width}\\right\) x + \\left\(\\frac{tm_{lng} - bm_{lng}}{height}\\right\) y + \\left\(mm_{lng} + \\frac{mr_{lng} - ml_{lng}}{2} + \\frac{tm_{lng} - bm_{lng}}{2}\\right\)](https://chart.googleapis.com/chart?cht=tx&chl=x%27%20%3D%20%5Cleft\(%5Cfrac%7Bmr_%7Blng%7D%20-%20ml_%7Blng%7D%7D%7Bwidth%7D%5Cright\)%20x%20%2B%20%5Cleft\(%5Cfrac%7Btm_%7Blng%7D%20-%20bm_%7Blng%7D%7D%7Bheight%7D%5Cright\)%20y%20%2B%20%5Cleft\(mm_%7Blng%7D%20%2B%20%5Cfrac%7Bmr_%7Blng%7D%20-%20ml_%7Blng%7D%7D%7B2%7D%20%2B%20%5Cfrac%7Btm_%7Blng%7D%20-%20bm_%7Blng%7D%7D%7B2%7D%5Cright\)) ![y' = \\left\(\\frac{mr_{lat} - ml_{lat}}{width}\\right\) x + \\left\(\\frac{tm_{lat} - bm_{lat}}{height}\\right\) y + \\left\(mm_{lat} + \\frac{mr_{lat} - ml_{lat}}{2} + \\frac{tm_{lat} - bm_{lat}}{2}\\right\)](https://chart.googleapis.com/chart?cht=tx&chl=y%27%20%3D%20%5Cleft\(%5Cfrac%7Bmr_%7Blat%7D%20-%20ml_%7Blat%7D%7D%7Bwidth%7D%5Cright\)%20x%20%2B%20%5Cleft\(%5Cfrac%7Btm_%7Blat%7D%20-%20bm_%7Blat%7D%7D%7Bheight%7D%5Cright\)%20y%20%2B%20%5Cleft\(mm_%7Blat%7D%20%2B%20%5Cfrac%7Bmr_%7Blat%7D%20-%20ml_%7Blat%7D%7D%7B2%7D%20%2B%20%5Cfrac%7Btm_%7Blat%7D%20-%20bm_%7Blat%7D%7D%7B2%7D%5Cright\)) where (mm_lng, mm_lat) is the position of the middle-middle point, (mr_lng, mr_lat) is middle-right, (tm_lng, tm_lat) is top-middle, etc. for the middles of the four sides. This is a linear approximation to the full coordinate transformation, minimizing error at the center of the image. The World Image format describes an affine transformation in the same way: the six parenthesized terms above can be interpreted as the six lines of the `.pnw` file (in A, B, C, D, E, F order) if the map units are longitude, latitude coordinates (that is, WGS84 a.k.a. EPSG:4326). The World file Wikipedia page says "World files describing a map on the Universal Transverse Mercator coordinate system (UTM) use these conventions..." and goes on to describe the UTM application, though it also says, "World files do not specify a coordinate system; this information is generally stored somewhere else." If the result of the `.pnw` conversion is interpreted as UTM, then I would need to do additional work to transform from WGS84 into UTM, and GeoServer would also need more information to specify which of the 60 UTM zones that the `.pnw` converts into. I don't see how to provide that information to GeoServer. Am I right in assuming that GeoServer interprets the `.pnw`-converted coordinates as WGS84?


title: How to render a static map from CartoDB?

text: CartoDB is working great for dynamic maps, but I have a use case to render a totally static map, just an image online. Google offers its static map api. Afaik there is no similar API for CartoDB. But does anyone have a work around that would let me display a single, non-zooming non-panning image? It'd be ideal if I could just put it in an html img tag, but something with javascript could potentially work. It'd be ideal if it's lightweight, as many of these users have very slow connections.


title: CartoDB data not displaying on Google basemap

text: I'm new to CartoDB but hoping to use it to display various types of real estate data on my website. That said I've been playing around with some of the sample code (Google Maps V3 integration example - https://github.com/CartoDB/cartodb.js/blob/develop/examples/gmaps.html) and cannot get this data to display within my website (see test page - https://metropolisiq.com/MapTest). I've tried the same code within Dreamweaver and it works fine. Not sure what the issue is, but perhaps it relates to the differences between Dreamweaver and my website:   1. My site is encrypted (HTTPS)   2. My site was built using the Joomla CMS. Any thoughts on how to solve are greatly appreciated.


title: ArcGISDynamicMapServiceLayer is not refreshing properly

text: I have a DynamicMapServiceLayer that I am refreshing every 30 seconds. I want the data in the map to refresh every 30 seconds too but it doesn't appear to be working properly. The map refreshes but then the service goes away and I just see a blank map. The service will stay off for about 20 seconds then come back on again. My code is here too...Is there a way to keep the service on all the time and still refresh it or make it so it doesn't go away? Thanks in advance!               var map;         require(["esri/map", "dojo/domReady!", "esri/layers/ArcGISDynamicMapServiceLayer"], function (Map, domReady, ArcGISDynamicMapServiceLayer) {             map = new Map("mapDiv", {                 center: [-56.049, 38.485],                 zoom: 2,                 basemap: "streets"             });             var liveLayer = new ArcGISDynamicMapServiceLayer("http://server/MapServer",              {                  refreshInterval: .5,                 disableClientCaching: true             });             console.log(liveLayer);                  map.addLayers([liveLayer]);         });


title: Selecting features through a many-to-one relationship

text: I want to select multiple counties based on the occurrence of a species present within them. I have a County polygon layer with CountyName as an attribute, a flat table called tblSpCnty with the attributes SpeciesName and CountyName (multiple CountyNames for each SpeciesName, representing all the counties in which each species occurs), and a flat table called tblSpeciesList with the list of species (SpeciesName, one record for each species). What I'd like to do is select a single species from tblSpeciesList, resulting in the selection of all the county polygons in which the species occurs.   I used to do this by linking tblSpeciesList through tblSpCnty to the County layer in ArcView, but I haven't been able to replicate the functionality in ArcMap. Using relates doesn't seem to achieve the same result. I'm not tied to ESRI, in fact if there's an Open Source solution all the better.


title: What are lazy geometries

text: What are "Lazy Geometries"? How are they different from other types of geometries? I'm an experienced Python/Django developer but am very new to the world of GIS, and understand that this is a newbie question. I've been unable to find the answer on google (if anyone can direct me towards a good tutorial site, I'd be grateful). I encountered the term when reading up on GeoDjango (and it's use of the GEOS library). Thanks in advance.


title: How can I speed up my arcpy script?

text: I have written an application that does a lot of geoprocessing using arcpy. Once started, the application itself runs at an acceptable speed, but it is very slow to boot. It can take 30 seconds or so to start up the GUI, so I'd like to find ways to make it start faster. As far as I can tell, the main source of delay is when I set the dataframe as follows:               self.df = arcpy.mapping.ListDataFrames(self.mxd, "")[0]      Does anyone have any tips / tricks to make either the whole application or the ListDataFrames run quicker when the application is started? I use ArcGIS 10.1, python 2.7 and tkinter for the GUI. Feel free to ask if I've forgotten to mention anything. EDIT: As requested, here is the (very!) stripped down version of the start my code:               import arcpy     import sys     from Tkinter import *     import ttk     import tkMessageBox     import tkFileDialog     import webbrowser     import time     import csv     import os     import textwrap          arcpy.env.workspace = r".\..\Append\temp.gdb"          Title = "Automated Map Generator"          class Application(Frame):         def __init__(self, master=None):             Frame.__init__(self, master, background = "white")             self.grid()                  self.createWidgets(master)              def createWidgets(self, master):                  # setup map document             self.mxd = arcpy.mapping.MapDocument(r".\autolim_mapping.mxd")             arcpy.env.overwriteOutput = True             self.df = arcpy.mapping.ListDataFrames(self.mxd, "")[0]          #[Lots of missing code here]          root = Tk()     root.tkraise()     app = Application(master=root)     app.master.title(Title)     app.mainloop()      How do I know what it's the ListDataFrames that takes so long? I put a whole heap of print statements through the startup code and noted how long they took to come out... I bet there are better ways to do this, but I'm pretty new to python. EDIT 2: @dassouki - thanks! Here are the output times from after 'import time':               Start Class Application (Frame):  1.53200006485     Start def Create Widgets:         1.59400010109     Start  mapping.MapDocument:       1.6099998951     Start env.overwrite:              1.64100003242     Start mapping.ListDataFrames:     1.64100003242     Finish mapping.ListDataFrames:    12.2660000324


title: interpolate point values to line

text: This might be a stupid question and I apologize if so. I have count data for different trails. They are points and I would like to make an average of those counts for each trail so that each trail has a count value. I don't know if I am being clear... Later on I would like to have a probability of encoutering the different species that were counted along those trails depending on the average counts. I am fairly new to GIS and I have no idea which tool I could use to do that... I was thinking of interpolation but I want to model the trails only, not the whole layer. I hope I have been clear enough and you will not mind my ignorance, Thank you,


title: Using Value from Iterate Feature Selection to Calculate Field in ArcMap 10.1 ModelBuilder

text: I have two polygon feature classes, an alignment layer and a parcels layer. I am trying to use ModelBuilder to create a tool that will do the follow:   1. iterate through groups in a single field of a polygon alignment layer (ie select and iterate through certain attributes in a field with same names such as "tunnel", "portal", etc)   2. Clip a parcel layer by each group    3. Populate a new field in the parcel layer called "EsmtType" with the name of each group (group is called "Value" in model)   4. Merge the clipped groups into one feature class I have gotten steps 1 and 2 to function but I am stuck on step 3. I am not sure how I need to configure the "Calculate Field" tool in order to recognize the "Value" variant of the "Iterate Feature Selection" tool. I tried using the advice from this post Using a value from "iterate feature selection' to calculate a field but I couldn't figure out how to implement it for my situation. Here is the model I have thus far: http://imgur.com/ji2nAJJ I would also like advice on how to implement step 4, I am not sure if it would require integrating a model within a model. Thank you so much!!! -kkaszas


title: Part of text missing in legend and text boxes

text: I have had issues lately when exporting to pdf from ArcMap 10.2. Parts of the text sometimes fall off and goes missing. Mostly it's been some lines in legends, but now also text boxes. The text is also missing in legends displayed in Arcmap, but not always. ![enter image description here](http://i.stack.imgur.com/QkEUo.png) This is not how the text in the textbox supposed to look, obviously. ![enter image description here](http://i.stack.imgur.com/Vl9p7.png) This screenshot shows a legend, but some of the descriptions are missing. In Arcmap the legend sometimes looks like this: ![enter image description here](http://i.stack.imgur.com/MDoJD.png) When it's supposed to look like this: ![enter image description here](http://i.stack.imgur.com/axc9X.png) It happens that this gets solved when exporting a second (or third) time. Which is good, of course, but does not help me in trying to figure out why it's happening in the first place. Has anyone else encountered this problem? Is there a solution? This is really starting to get on my nerves, and not even creating a new legend helps every time (which is very tideous in the first place).


title: What if format for Getcapabilities document

text: I want to add capabilities URL for cascading a layer from remote WMS in Geoserver. What will be the exact format for Getcapabilities URL ?


title: Pick object by properties (and centre on it)

text: I have a Leaflet map with a Cloudmade base layer and a single layer of polygons. Each of these polygons has a reference number attached as a property. I would like to select a reference number from a list (or whatever) and then pick out the object by its reference number and zoom to the extents of that object. How can I pick an object by its properties? (Presumably once I have done that zooming to extents shouldn't be too hard)


title: reproject wms in geoserver

text: > **Possible Duplicate:**   >  Re-projection in Geoserver I've got a WMS service that publishes in WGS84 (EPSG:4326). I have several AutoDesk users that need to use the imagery in NAD83 (EPSG:2227). I have geoserver 2.1 setup to cascade the wms and reproject it. Problem is that as soon as I turn on SRS handling to reproject to declared SRS the reprojected wms fails. What am I doing wrong? Is Geoserver capable of doing the reprojection.


title: Attribute Table Decimal Places

text: As we all know, the attribute table field properties has number format to display number of decimal places to up to maximum 15. If we have many fields/columns in attribute table, how can we change all fields/columns decimal places roundoff/display at a time ( I mean in single operation. Not changing each and every field properties). Ex:               ColA:12.987654 ! ColB:34.20987766  ! ColC:89.654309877      In one operation my output should be (Rounding/display to 2 digits) For Rounding:               ColA:12.99 ! ColB:34.201  ! ColC:89.66      For Display:               ColA:12.98 ! ColB:34.20  ! ColC:89.65      is there any method/tool to do this in Arcgis?


title: Fix topology error tool- can't automatically fix multiple errors?

text: I am a ArcGIS.10. user. I'm want to solve multiple errors using the fix topology error. So I select all my features I click right but I don't have the correcting options (Snap, trim, extent), I only have Delete. Does anyone know why? Thank you very much,


title: Can someone help me to convert this Python to VBscript for labellling

text: I have two Python script snippets that need to convert to VBscript. Can anybody help me because I am really not familiar and know Python. the first script is               def FindLabel([BA2]):         name = str([BA2])         if "-" in name:             first, second = name.split("-",1)             for letter in first:                 if not letter.isalnum():                     first = first.strip(letter)             return first         else:         return name      the second script is:               def FindLabel([BA2]):         name = str([BA2])         if "-" in name:             first, second = name.split("-",1)             for letter in second:                 if not letter.isalnum():                     second = second.strip(letter)             return second            else:             return name


title: Problem with snap layer when in name is space

text: I'm trying to make a option, where user can choose layer to snap. All is okey when layers names are without space, but if in name is space (like 'something layer name') snapping does not work. How can I solve it? thank for all posts.


title: Pause ArcPy script until feature class lock goes away

text: My script needs to copy data into a feature class to an SDE geodatabase that might have a lock on it from another application. I shouldn't break the lock (that app has priority over a data update), but I also don't want the script to quit just because it tried at the wrong moment. How can I set it to keep trying until there is no lock? I'm looking for something like a try/except, but more like try/keep-trying-every-15-seconds. I considered making a separate function and adding a pause:               try:         funcCopyData()     except:         pause 10 seconds         funcCopyData()      but wouldn't this still be limit me to only two attempts?


title: QGIS : launch Atlas composer generation from command line

text: I am using QGIS 2.2 and have configured a composer with atlas generation. I would like to launch this atlas generation every day. The idea is to write a batch and launch it with a cron. Is it possible to launch atlas generation from command line ? Thanks


title: Export Stretch Function values - Arcpy

text: I'm looking for a way to export Stretch Function values from images that’s stored in a Mosaic dataset. Not values from the whole set, but specific stretch functions per image. The set is stored in either a filegeodatabase or SDE. The values I’m looking for is eg. Band[] Min, Max, Mean values and Gamma values.... EditRasterFunction_management() allows me to insert, update or delete which I’ve used to initially create the Stretch Function. But after some manual editing for the raster images I need a way to export the new values. I can’t manage to get access to this information using cursors since it’s a function? Any ideas how to maintain this function values? Preferably using Arcpy/python. Thanks!


title: arcgis silverlight - how to delete/cut a part of a multipart polygon

text: Say I cloned a single-part lake polygon that had a hole for an island. Then I used the Editor Reshape command to cut a piece off. That incidentally caused a second part to be created where the island was, as seen here. ![enter image description here](http://i.stack.imgur.com/v1J0v.png) Now I want to get rid of that bogus island. I've tried drawing a circle around it with the Cut command, tried the Reshape command again, no joy. Help a dummy out? (yeah I know my links are to ArcGIS Desktop editing tools, the better descriptions of the apparently more or less equivalent Silverlight API Editor commands) ;-)


title: Missing tiles when django application is hosted on an apache server

text: I'm building a GIS application using Django as my web framework and OpenLayers to implement a slippy map interface. I was initially using the django runserver for development and everything was working fine. However, I recently ported my app on an Apache server and I started seeing some missing tiles (see here): It's not always the same tile that is missing and usually there are two or three missing tiles. Sometimes, after refreshing the page several times the tiles appear properly. There are no missing tiles when I zoom in or out. Note that the page keeps on loading for an indefinite amount of time when I experience this issue. I'm using a TMS layer to display the map using data from the GSHHS database.The map is being rendered via a mapnik XML file which uses the GSHHS shapefile. Unfortunately, at this point I'm unable to paste any code over here since I'm unable to login to SE from work. I'll try to post some code later. I have no idea as to what might be causing this problem. There is nothing helpful in the Apache and Django log files. Any help would be greatly appreciated. Thanks in advance.


title: Problem with MIF/MID allignment

text: After creating a rather large map in MapInfo using US Census Tiger data and exporting it to the MIF/MID format, it appears the files are no longer properly aligned. When I try to open the MIF/MID in some other software it gives me an error. "Got Pline in the middle of reading other Pline" I assume this means that the MIF formatting is incorrect. I have recreated the MIF/MID files and the problem continues. Maybe a problem in the US census data? I am quite experienced with doing this process and I have never had this problem before. Help!


title: Why Add Join does not work in ArcMap and Arcpy?

text: I have problem with joining feature class with table from the same location (GDB) using Python script. When I use join Join and Relates -> Join... in ArcMap TOC it does not work because I do not get any joined records despite the fact there are equal values in fields used as a keys . The same problem occurs using arcpy. AddJoin_management(). When I use arcpy.JoinField_management () or the same tool from Toolbox in ArcMap I get the results. What is the reason that JoinField works and it gives me a good results, but in the case of AddJoin I do not get any joins? The field type in the Feature Class is the same as in the table to join. Any suggestions?


title: Find tunnel 'center line'?

text: I have some map files consisting of 'polylines' (each line is just a list of vertices) representing tunnels, and I want to try and find the tunnel 'center line' (shown, roughly, in red below). ![alt text](http://i.stack.imgur.com/jZdYb.png) I've had some success in the past using Delaunay triangulation but I'd like to avoid that method as it does not (in general) allow for easy/frequent modification of my map data. Any ideas on how I might be able to do this? Thanks. P.S. I'm working in fairly raw C++, and was referred here from Stack Overflow.


title: What are the most tangible benefits to getting your GISP or ESRI Certs?

text: > **Possible Duplicate:**   >  GISP - Value of being certified? Looks like I'm only a few presentations or other "contributions to the Profession" points away from being able to apply for my GISP. Has anybody found any actual benefits to getting the certification? Seems like it might be somewhat arbitrary. Of course it lets people know you have a high level understanding of GIS but beyond that are there any benefits? Salary bumps? Other? **ADDED** Additionally, what, if any are the benefits to getting the new ESRI Certifications? ESRI is obviously the most well known GIS software, but there are others... so if one is looking to move forward in their GIS career, what is the best way to do it as far as certs/continuing education go?


title: Point in polygon algorithm for multiple polygons

text: I've a Google map with bunch of polygons on it. Here is a problem I'm interested in: Given a lat,lng point, what's the best way to determine all the polygons that this point lies in? The obvious way is to run a "point in polygon" algorithm iteratively for each polygon, but I was wondering if there an efficient algorithm to answer such queries especially if you have thousands of polygons.


title: Is there a way to add turn restrictions in A* and Dijkstra?

text: We are currently working with pgRouting and found that there's no way of implementing turn restrictions (e.g. banned right or left turns). While it is possible to assign "to_cost" and "rule" in Shooting * algorithm... I couldn't find a way of implementing these restrictions using A star and Dijktra algorithms. Is there a way of implementing special turn restrictions in A star and Dijkstra methods ?


title: Command r.sunmask in GRASS

text: I have hypsometric model of relief, and I need to visualize shadows. I was suggested to use GRASS command r.sunmask, so I tried to set time values and run command, but it gives me an error               "ERROR: Neither sun position nor east/north, date/time/timezone definition are complete"      So then I tried to put values in Position tab, which I founded by SOLPOS Calculator (http://www.nrel.gov/midc/solpos/solpos.html)               r.sunmask elevation=model@PLANINA output=model_sunmask year=2011 month=1 day=15 hour=9 timezone=1 east=45.909361 north=15.965538      But it didn't work out, too, it says "Calculating shadows from DEM..." but nothing's going on for hours of waiting. So I ask if that command works at all in GRASS, or I'm doing something wrong? You can suggest some other application for same thing, except ArcGIS applications. Thank you.


title: Save spatial reference to application settings?

text: Having created a tool as a macro for several geoprocessing tasks, I'm saving user input to application settings so the user doesn't have to enter all required data on each use. I'm stuck on saving the spatial reference. ESRI code shows using the FactoryCode and an enumeration to create either a projected or a geometry spatial reference. What I'm trying to figure out is how to get the enumeration, and how to identify whether the factorycode/enumeration points to a projected or geometry spatial reference. I will greatly appreciate any answers/directions provided.


title: How to set the value of a field in all selected rows of a table using ArcGIS 10

text: I have an attribute table in ArcGIS (it's large - around 300,000 elements). I've selected some rows of the table using 'Select by attributes' and I want to set the value of a field for all of those rows. How should I go about doing this? I'm hoping there will be a nice, easy, GUI- based way - as I'm teaching someone with very little GIS experience.


title: What is future of Silverlight with ArcGIS Platform?

text: At the September 2011 Build Conference, Microsoft said a lot about Windows 8 and the Future of JavaScript over Silverlight. No Plugins wanted!! JavaScript with html5 ftw! What is the future of Silverlight with the ArcGIS Platform (ArcGIS for Server, ArcGIS Online, etc) in 2014 and beyond?


title: cad and geoserver

text: Hi Guys i have geoserver and openalyers running quite nicely and now need to overlay a cad drawing. I have looked at the geographic location part of cad and set the long / lat to the appropriate values. I saved the drawing and then converted to a shapefile and then imported in to postgis. When i add to geoserver the nativebounds and lon/lat generated in geoserver are no where near where i would expect them to be. Can anyone point out what i am doing wrong. When i set the long/lat in my drawing where does this point reference. is it the center, top and left, bottom and left of the drawing? Also how can i find out what projection the drawing is in. It doesnt seem that you can set this in cad. Ive asked my cad guys and no one seems to know. Thanks in advance


title: Is there any feature limit for Spatial data migration to Postgre SQL9.2...?

text: I am trying to migrate spatial data to Postgre SQL9.2 using 'Postgis Shapefile & DBF Loader 2.0'....Each time it returns a "Runtime error'. Primary data extension - .shp ; Data type - Multipolygon ; Feature Count - 1116930 . Is this a problem due to the lage amount of data involved....!!! Thanks in advance..


In the definition of **surface plots** in my book, (An introduction to data mining, chapter 3) it is said that: > Surface plots are often used to describe mathematical functions or physical > surfaces that vary in a relatively smooth manner. I don't understand what exactly is meant by the above sentence. What special property do mathematical functions and physical surfaces have that is appropriate for the case of surface plot?


title: maxScale or minScale dont working in vector (from wfs) layer

text: var options =          {             projection: "EPSG:900913",             displayProjection: "EPSG:4326",             units: 'm'         };         map = new OpenLayers.Map('map', options);                    ilceSinir = new OpenLayers.Layer.Vector("İlçe Sınırı",         {             minScale: 150000,             strategies: [new OpenLayers.Strategy.Fixed()],              projection: new OpenLayers.Projection("EPSG:4326"),             protocol: new OpenLayers.Protocol.WFS(             {                 version: "1.1.0",                 url: "http://localhost:8080/geoserver/wfs",                 featurePrefix: 'topp',                 featureType: "g_ilce_siniri",                  featureNS: "http://www.openplans.org/topp",                 geometryName: "bounds"              })         });      scale isn't working. layer is visible everyzoom.


title: Heat Map trouble: Is this a data formating problem?

text: I'm playing with the heat map plugin with QGIS. I've run through QGIS tutorial: http://docs.qgis.org/html/en/docs/user_manual/plugins/plugins_heatmap.html,   using the airports data set as explained and I reproduced tutorial results. I have a 600,000+ data set (five years worth) of police calls from Philadelphia, which I've reduces to 185 points (a day) for testing. The Meta data for both sets is below although I've used both the dbf file and saved shp file as input. For larger Radius settings (above about 0.1) I get a rainbow bulls Eye centered on but many time lager than my data set range. For smaller Radius settings I get a uniform single color covering the data set range and transparent holes opening over (but not centered on) each data point. I get just what I should with the Alaska data set. General:               Storage type of this layer: ESRI Shapefile     Source for this layer: C:\Users\Gallery of Tools\Documents\GIS DataBase\Alaska\shapefiles\airports.shp     Geometry type of the features in this layer: Point     The number of features in this layer: 76     Editing capabilities of this layer: Add Features, Delete Features, Change Attribute Values, Add Attributes, Delete Attributes, Create Spatial Index, Fast Access to Features at ID, Change Geometries     Extents:     In layer spatial reference system units : xMin,yMin -4480198.52,1433525.80 : xMax,yMax 4615124.98,6502586.83     Layer Spatial Reference System:     +proj=lcc +lat_1=40.96666666666667 +lat_2=39.93333333333333 +lat_0=39.33333333333334 +lon_0=-77.75 +x_0=600000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs               General:     Storage type of this layer: ESRI Shapefile     Source for this layer: C:\Users\Gallery of Tools\Recovered from the Dell600GB\Image Catalog\EngravingsNMaps\Philadelphia\Police Incident\PoliceIncJune302012.dbf     Geometry type of the features in this layer: Point     The number of features in this layer: 185     Editing capabilities of this layer: Add Features, Delete Features, Change Attribute Values, Add Attributes, Delete Attributes, Create Spatial Index, Fast Access to Features at ID, Change Geometries     Extents:     In layer spatial reference system units : xMin,yMin -75.2515,39.9047 : xMax,yMax -74.9594,40.1226     Layer Spatial Reference System:     +proj=lcc +lat_1=40.96666666666667 +lat_2=39.93333333333333 +lat_0=39.33333333333334 +lon_0=-77.75 +x_0=600000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs      P.S. And another thing. When I change the Radius setting, the cell size changes automatically. Is this supposed to happen?


title: How to intersect 2 polygon layers with a given tolerance for features in different layers only?

text: I can intersect 2 layers with polygons using ArcMap. However I can specify only overall tolerance no matter what layer a feature is in (even same layer). So I either end up with lots of slivers if I don't tolerate snapping, or I end up collapsing small features from the first layer that I'd want to keep. The problem is that I would like for features let's say less than 10 m apart in the first layer to stay as they are if they happen to be contained within some polygon of the second layer. However I'd like for edges of polygons of the 2nd layer to snap to those in the first layer if they are within 10 m. I'm open to a solution using any tool be it ArcMap or some python library/script. Here is the rough idea of what I want. ![enter image description here](http://i.stack.imgur.com/0XHhb.png)


title: Tools to calculate routing distance

text: I'm working to a web application, and I need to calculate the routing distance (not the distance "as the crow flies") between a starting address and several final addresses. I don't need to render the routes on the map. The application will have a lot of queries so I can't use Google APIs and probably the best solution is to perform the task not with an external service but on a local server. I think the only way is to use Open Street Maps (or not?). Which tool to perform the calculation do you recommend? Thank you in advance! * _*_ edit:   * At the moment geocoding (700 conversions daily) is performed by GoogleMaps and Nokia's HERE APIs. Do you recommend me some other tools also for this?   * The number of daily point to point routing calculations are approximately 100.000


title: ArcToolbox 9.3 converted to ArcToolbox 10 and now won't work in ArcMap 9.3

text: I added a custom toolbox created for 9.3.1 to a ArcGIS Desktop 10 ArcMap. Now it seemed to have disappeared. I can't add it to ArcMap 9.3.1. Is there a way to convert it back if that is even the issue? Or is there something else going on? Thanks.


title: grass gis- ERROR: nwalk (7000001) > maxw (7000000)!

text: i want to make r.sim.water (hydrological modelling) in Grass Gis, and this ERROR: nwalk (7000001) > maxw (7000000)! always appears, any ideas how to solve it? number of cells of dem is below 7000000.


title: Geoserver administration - defining different privileges

text: I need to create separate admin accounts for Geoserver. I'd like to have one superadmin account that would be able to read/edit/delete any style, layergroup, workspace, store or layer - that is what the default admin is already able to. I also need to have several accounts that would be able to create workspaces, layers and stores, would not be able to edit styles and would be able to edit and delete only the workspaces, layers and stores they created (e.g. if superadmin creates a workspace, he is the only one to delete it). Is this possible to achieve with Geoserver?


title: Convert Shape Polyline-M to Polyline

text: I have some pretty old software that I am forced to use (no option). It's old enough that it doesn't recognize shape files of type Polyline-M, it only recognizes shape files of Polyline. I have a shape file that is completely Polyline-M that I need to convert to plain old Polyline. I'm struggling to find a way to do this in ArcMap 10. Does anyone know how to do this? A helpful link, or even just a point in the right direction to where I could read how to do it would be really handy. I'd rather not have to create a tool to do it, surely something like this exists out there already! :) Thank you


title: Getting what seems to be incorrect data from GeoNames

text: I have a hiking website. I am trying to query for parks big and small around the world where people could potentially go walking or hiking. I was trying to use the GeoNames data and was using query arguments from this table: http://www.geonames.org/export/codes.html Here is an example query I made trying to get parks:               http://api.geonames.org/search?maxRows=10&username=demo&q=PRK      and it gives me some things which are I have no idea what. Here is also another thing where I try to query for arctic spaces:               http://api.geonames.org/search?maxRows=10&username=demo&q=LAND       And that gave places in the Carribean :) So I must be doing something totally wrong. My goal was to look up all parks in the world, and also cities (with state/country info) so that I can make that connection on my site for my users. Is my goal possible to achieve? What was I doing wrong in my queries? Thank you!


title: Is there a tool for automated quality assurance for vector data captured on raster data?

text: I wonder whether is any tool(or script) in ArcGIS or QGIS which allows to make an automated quality assurance for vector data captured on raster data? vector data on binary raster: ![enter image description here](http://i.stack.imgur.com/ewmav.jpg)


title: Import KML attributes into the correct schema

text: I have a kml file containing over 4000 polygons with attributes which must be retained (in 3 different fields). After using the standard `KML to Layer` tool, all the attributes are now in one `PopupInfo` field, along with the html code. I have the `Data Interoperability` extension, so after looking at the questions here and here, I came across a post on the ArcGIS forums: > The Data Interoperability extension will read KML (including folders and > subfolders) and use any schema object defining the attributes. Thus what it > writes out will have the correct field definitions. The same goes for > writing KML with a schema object, given geodatabase inputs. What I understand from that is that using the `Quick Import` tool will read the KML and export it to a gdb based on a schema which I defined. This seems to be exactly what I want, except I cannot find the option to specify the schema. I exposed the `kml_name` attribute as mentioned in the other question, but the field remained empty. I tried downloading the KML2shapefile plugin with MapWindow, but it is only for 32 bit and I am running 64 bit. I am unable to download the FME trial at work.


title: Background (sky/terrain) of a viewshed

text: Is it possible to determine whether a point feature of a particular height would be backed by sky or by terrain? I've done the viewsheed analysis of a wind turbine. Now I would like to make a raster map representing the background (sky/terrain) from each cell of the viewshed.


title: Nodata issues in map algebra

text: A dummy question. I need to change all values in my rasters into ones. I've been able to do this fine in the Rater Calculator just by multiplying the layer by zero and adding one (my_layer*0+1), but when I'm doing the same in arcpy, for some reason it works for a part of my raster but for some it changes all the nodata cells to ones as well. The nodata values in different rasters are different because my rasters have been generated from different source data. Does it not recognize some nodata values? If so, this is odd as it displays them correctly when I just import them on the mxd. Any ideas on why this works sometimes, but sometimes not? Here's the exact script:               bin_test = arcpy.Raster("F:/datasets/map1.img")*0+1


title: can't generate/display signed float values using python/gdal

text: I'm trying to create a raster with signed float values as the pixel values and display the results in MapWindow on windows. What is the proper gdal band type to use? I've tried both **gdal.GDT_Byte** and **gdal.GDT_Float32** and using both methods I can't seem to get the pixels displayed. If I change the pixel value to the abs() value they will display.               dataset = driver.Create(filename, cols, rows, number_of_bands, gdal.GDT_Byte)      **gdal.GDT_CFloat32** allows for negative numbers but doesn't properly display the floating point value past the point.


title: Integrate QGIS with ASP.NET MVC3

text: I'm completely new to GIS and i'm trying to work my way into this. My goal is to integrate a simple map created using QGIS with my asp.net web mvc3 application. The main objectives would be to show the map on the application as well as performing some operations with the different points and lines (points of interest and roads respectively). Can anyone point me in the right direction on how to achieve this?


title: Does the esri javascript api Identify Task Return a FeatureSet?

text: I am trying to get the results of an identify task using ArcGIS JS API 2.7 into a featureSet so I can utilize this sample http://arcscripts.esri.com/details.asp?dbid=16528 that takes the results of a datagrid store and exports them to an excel file. I have it working where one feature is being exported, I just need to have all selected features export. Does an Identify Task produce a FeatureSet? Thanks!                function executeIdentifyTask(geom) {       //clear the graphics layer        map.graphics.clear();            identifyParams.geometry = geom;       identifyParams.mapExtent = map.extent;       identifyTask.execute(identifyParams,function(response){           var polygonSymbol = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_SOLID, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_DOT, new dojo.Color([151, 249,0,.80]), 3), new dojo.Color([151, 249, 0, 0.45]));          var markerSymbol = new esri.symbol.SimpleMarkerSymbol(esri.symbol.SimpleMarkerSymbol.STYLE_CIRCLE, 20, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID, new dojo.Color([0, 0, 0]), 1), new dojo.Color([25,50,225,0.3]));               var controlItems = [];          var surveyItems = [];               dojo.forEach(response,function(result){                 var feature = result.feature;                 if (result.layerName =="surveys"){         showSurveysNameGrid();             feature.setSymbol(polygonSymbol);             searchType="selSurvey";         surveyItems.push(attributes);                   var attributes = feature.attributes;         var csvTextSurv = "DOCUMENT_N;TOWNSHIP_RANGE;SECTION\n";                     csvTextSurv += attributes["DOCUMENT_N"] + ";" +                                attributes["TOWNSHIP_RANGE"] + ";" +                                attributes["SECTION"] + "\n";                  document.getElementById("csvSurv").value = csvTextSurv;              map.graphics.add(feature);        }else{          showPointNameGrid();         feature.setSymbol(markerSymbol);         var attributes = feature.attributes;         controlItems.push(feature.attributes);          searchType="selControl";          var csvTextMons = "POINT_NAME;SECTION\n";                csvTextMons += attributes["POINT_NAME"] + ";" +                                attributes["TOWNSHIP_RANGE"] + ";" +                                attributes["SECTION"] + "\n";        }               document.getElementById("csvMons").value = csvTextMons;         //add selected feature to graphics layer                   map.graphics.add(feature);      });               if(surveyItems.length >0){          showSurveysNameGrid();          var surveysStore = new dojo.data.ItemFileReadStore({data:   {identifier:'DOCUMENT_N',items:surveyItems}});             var grid = dijit.byId('grid5');             grid.setStore(surveysStore);          }      if(controlItems.length >0){         showPointNameGrid();         var controlStore = new dojo.data.ItemFileReadStore({data:{identifier:'POINT_NAME',items:controlItems}});         var grid = dijit.byId('grid4');             grid.setStore(controlStore);                }                         });      }


title: Bearing/Course/Heading in WGS84 ellipsoid model

text: There is a common formula for computing the initial bearing for a course from `lat1, lon1` to `lat2, lon2`. (And yes, I'm aware that the course will change when traveling).               atan2( sin(lon2-lon1) * cos(lat2), cos(lat1)*sin(lat2)-sin(lat1)*cos(lat2)*cos(lon2-lon1))      The formulary where I took this from mostly uses a spherical model for simplicity. My questions about this formula are:   1. What is the appropriate formula for an ellipsoid model such as WGS84?   2. If it is the same formula, does this hold generally? I.e. is it the same formula for any ellipsoid model because we deformed the coordinate system, not the coordinates? Somewhat like Vincentry's formula, which uses a reduced latitude `u = atan((1 - f) * tan(lat))`. Also note that I'm looking for a global solution, including for the poles.


title: Confusing offsets in shapefile written by QGIS/OGR

text: I think I'm stuck here. A shapefile has been updated in QGIS. Records have been deleted and those records have been marked with an asterix in the first byte of the record in the dbf file. I open the shapefile in a C program to convert it to postGIS. The deleted records are skipped. Things happen to the next record after a deleted record, not after the first deleted record though. The shx file says that the offset for the next record after the deleted record is 30948 and the accumulated offset in the shp file says 28740. Neither is right, or rather when you start reading in the shp file at the two offsets you land up in the middle of a point stream and not at the start of a record header. The offset 28740 should be right because that is what you get when you add up the offset of the deleted record and the contentlength of the deleted record. The conclusion is that the shape is corrupt. But what is driving me nuts is that you can load the shape file into QGIS and it is shown correctly. What have I missed.How can QGIS handle this.


title: Google Maps API: how to get the street name in multiple languages

text: I'm trying to get the street name in multiple languages from the Google Maps API V3. For example in Belgium there are multilingual street names: ![enter image description here](http://i.stack.imgur.com/HU0JX.png) If i set the `language=FR` i will get "Rue Brogniez" and `language=NL` results "Brogniezstraat". How can i get these names together in one `geocoder.geocode` call?


title: How do I obtain layout coordinates of features using arcpy?

text: I just read the topic about Convert point XY to page units XY using arcpy? but I didn't find a solution to my problem. I have a point layer file, I know coordinates in map view, I'd like to know coordinates in layout view. I suppose the main difference with the case of a dataframe is that the point is not a graphic element itself in layout view, so the referencing of the point is more tricky. Is there anyone who solved this (or similar) issue? ps. the dataframe where I'm working is rotated, the solution proposed in the Convert point XY to page units XY using arcpy? post doesn't work for rotated dataframe. Does anybody knows what is the base point for rotation of dataframes in ArcGIS?


title: Oracle spatial. Cannot load data

text: Have just installed QGIS 1.9 with Oracle spatial abilities. We maintain geodata in an Oracle 9.2. spatial database. I get connected to the database but I can not load the data. The table shows up with no geometry and contains "nothing". Is there anyone else who has tried this ? Any suggestions ? Svein Inge Gausemel Lillehammer


title: c# How can I get the spatial reference of a ArcGIS layer?

text: I m using the following c# code to add my layer dynamically on a arcgis silverlight application:               ArcGISDynamicMapServiceLayer mylayer = new ArcGISDynamicMapServiceLayer()     {        ID = layerName,        Url = mapServiceUrl,      };      mylayer.Initialize();      mylayer.Initialized += mylayer _Initialized;      mylayer.InitializationFailed += mylayer _Failed;      myMap.Layers.Add(mylayer );      But when I try to get the spatial reference of the layer using:               SpatialReference sp = mylayer.SpatialReference;     int testww = sp.WKID;      I get an error "Object reference not set to an instance of an object." because mylayer.SpatialReference is Null. Where m I wrong? Thanks in advance!


title: Display GeoJSON layer on top a SVG image instead of a map.

text: I'm very new to GIS, so hope this is the correct place to ask. Not sure if this will even be possible. **This is what I want to achieve:** In short, I want to display a set of polygons, lines etc. in a browser, and then display an overlay which marks specific polygons and lines etc. **This is the way I think it should be approached:**   * Take a SVG which represents a set of lines and polygons.    * Display this SVG using something like leaflet in a browser, but without displaying a world map which seems to be the basis of leaflet, open layers etc.   * Overlay a GeoJSON layer on top of the SVG.   * If possible, I want the coordinates in the GeoJSON file to be relative to a reference point in the SVG. Is this possible at all? Also, is the last point possible, that is: specifying what the coordinates in the GeoJSON file should be relative to, or are they always interpreted as GIS coordinates? Thanks for any inputs, Jaco


title: Configure QGIS to show only the toolbars wanted.

text: I friend of mine who are using QGIS in their organization for non GIS- experienced people asked me a question the other day that I didn't knew the answer to. He would like to configure QGIS to only show the tools/toolbar that those non-experienced people need to see. A common complaint is that QGIS looks too so difficault with all toolbars. QGIS is installed to this users on a server and he would like to control this so that people don't need to change the toolbars on every computer that are using QGIS. Has anyone tried to solve this kind of problem (or that understand what the hell I'm talking about ;))


title: How to import a multiband image into IDL?

text: I want to run mathematical calculation on monthly NDVI images from april to october. I stacked monthly NDVI images through april to october into an image (with 7 band) using ENVI. Now I wanna import this image into IDL. How can I do this? What kind of function or procedure should be written? Please help me?


title: Which tool is recommended for a PostGIS design?

text: I need to create a design for a spatial database. I've found this list http://wiki.postgresql.org/wiki/GUI_Database_Design_Tools, but it seems to me that this software is orientated towards a regular (not spatial) DB design.


title: Creation raster image from interpolated image

text: I have interpolated temperature data that I have. This image's max and min values are 11.9 and -1.7, respectively. And then I exported this image to raster using 'export to raster'. But the max and min values of the output raster image are different from that of the interpolated temperature image. I'm wondering why values of both images (temperature and raster) are not same. They need to be same. Is there any way to make the values of the raster image be the same as the interpolated image?


title: Changing Buffer table geometry dynamically in PostGIS using plain SQL

text: I have two tables in PostGIS, a point and a point_buffer table. The point table has a buffer_distance field, with a default value say 200. Now I want to change the buffer table geometry whenever I change the buffer_distance value in my point table. I can do this for a single row in my point_buffer table using the following:               UPDATE point_buffer     SET the_geom = (SELECT ST_Buffer(the_geom,500) FROM point WHERE gid = 1)     FROM point     WHERE point.gid = point_buffer.gid      But whenever I am trying to change the whole of point_buffer table (droping WHERE clause in my sub-query) I am getting error message: > 'ERROR: more than one row returned by a subquery used as an expression'. My question is, can I change the entire point_buffer table in one go? I know one option is using a for lopp, with the upper bound being the count value of my point table and incrementing the point.gid value . But I want to do this in **plain SQL**.


title: Process to load elevation data into Oracle 11G

text: We have the National Elevation data set from USGS and considering loading it into Oracle (11G) database.   Can someone help with the process/steps?


title: How to store a password during one QGIS work session?

text: I'd like to get into qgis-plugin development, but am rather new to both python and QGIS development. I intend to write a plugin which connects data from a mysql database to shapefiles. As the data is related to persons, I need a kind of login before accessing the data. To avoid the user from typing his password every time he invokes the script, I wonder if it is possible to retain values between calls of the script. In the python cookbook, I found global project settings - but that is not quite what I'm looking for as the values should be forgotten on closing QGIS. I hope I made my intentions clear in my clumpy English and somebody could give me a hint with which way to proceed.


title: Updating and exporting XY positions in UTM and WGS84 to spreadsheet using ArcGIS ModelBuilder?

text: On a daily basis I'm editing the geographic positions of various points of a shapefile (corner_XY.shp) by dragging them on my map. When this editing task is done I'm updating the points coordinates in their attribute table column by column by right clicking on the column and calculate geometry. I have 4 columns: Latitude - Longitude - X_UTM - Y_UTM Latitude/longitude in WGS84 (format= DD MM SS,SSS) and X_UTM - Y_UTM in WGS84 UTM Zone 15 N. I'm then exporting the attribute table with the updated fields (Latitude - Longitude - X_UTM - Y_UTM) to an excel spreadsheet to share it with my team. Exemple of output:               X_UTM        Y_UTM      Latitude            Longitude     596697.989  2276072.78  20° 33' 58.529" N   91° 9' 16.333" W     599675.251  2123364.47  20° 32' 28.838" N   91° 12' 18.624" W     575317.464  2263171.52  20° 23' 51.409" N   91° 6' 35.489" W     ....        .....       ......              ......      My project is a projected one in WGS84 UTM Zone 15 N and my shapefile containing the points that I'm moving on a daily basis (corner_XY.shp) is in WGS84 not projected. I'm using ArcGis 10.1. GOAL: I'm trying to replicate those tasks using the model builder in order to automatically update and export the new position of my points in WGS84 and UTM to an excel spreadsheet.


title: How to simplify vector lines into an overlaid grid in QGIS?

text: I'm a relatively new user with GIS software (primarily QGIS) and I'm sure this is probably quite simple but I've been trying to find an answer for ages! I have transportation network data for an area (as vector lines) and I'm wanting to simplify these to then attempt to find a spatial correlation with another set of points that I'm interested in. I was told the best way to do this was to overlay a grid over the vector lines, so that each grid cell could be classified as to whether or not it contains one of the lines or not.


title: Shapefile to Rectangle Overlay/Intersect Test?

text: I have a shapefile and a repository of images that have corner Lat/Long points. I want to search that images who overlaps/intersects to shapefile. Here is the scenario that I am doing, creating a polygon from each image corner lat/long points and placing it in temporary memory. Now how can I check whether this created polygon intersects/overlaps the shapefile? Any solution recommended for this issue?


title: Digitizing a city map for thesis project; Transportaion

text: So, basically I need a digitized map of a city that can be used for analysis. My thesis will be focusing on land-use destination, so I'll be needing the fully digitized map that has all kind of land-use categories plus the road networks, precisely the measurement of accessibility. The biggest hurdle that I've came across is, as I belong to the third world country i.e. Pakistan, the nonavailability of data for GIS. Besides, I am not an expert in the GIS. However, I am major in urban planning. Things that I would like to ask are:   * How can I extract image for the study area?   * How can I digitize the map in the fastest way?


title: Why do ArcGIS Area Calculations not match Official Statistics?

text: I am trying to calculate area of polygons using a shapefile from the world bank where countries are defined. When however I am trying to calcuate the area, the numbers that emerge are different from the official statistics on area. Morever, whenever I use a different projection (e.g. equal area) and recalculate the area, I always get the same results. Can someone tell me where the mistake could come from? Thanks!


title: Precision for an FGDB raster dataset?

text: I would like to know the precision for an ArcGIS raster stored within a file geodatabase (FGDBR format); I can't seem to find this information anywhere in their documentation (apologies if I have missed the obvious...). By precision, I mean the number of decimal places that are used to store the value for each grid-cell (not the coordinate values for spatial referencing). I know this is 6 for the GRID format, but can't find the equivalent info for the FGDBR format. The reason for asking is that I get different results for exactly the same processing when I use FGDBR vs. file stored GRID raster datasets. I am interested in the cumulative sum over the entire raster grid; there is a difference of approx. 0.5% between the two methods (i.e. a result of 5301.316 when the analysis is carried out in a GDB, compared with 5323.825 in a file storage system). EDIT: In addition, students have reported getting vastly different answers when using 64 vs 32 bit machines; I don't have the values for these though. I will add when I can repeat the results myself.


title: Print Google and WMS layers with MapFishPrint

text: My application print maps using mapfish print with geoserver. But an problem, when print layer google and wms the two images are not properly overlaying, in other words my layer wms of Brazil not not overlaying territory of Brazil in Google Maps Street. app/pdf/print.pdf?spec= {"title": "Printing Demo","layout":"A4 portrait","srs":"EPSG:900913","units":"m","dpi":254, "pages": [ {                   "bbox": [-6227217.67059, -3780794.2499401,  -691301.11508219,  586512.2832727],         "scale": 28000000,         "rotation": 0,         "comment":"",         "mapTitle":"",          "geodetic": true     }], "layers":[      {          "baseURL": "http://maps.google.com/maps/api/staticmap",          "type" : "Google",         "maxExtent": [-20037508.3392,-20037508.3392,20037508.3392,20037508.3392],         "resolutions": [156543.03390625, 78271.516953125, 39135.7584765625, 19567.87923828125, 9783.939619140625, 4891.9698095703125, 2445.9849047851562, 1222.9924523925781, 611.4962261962891, 305.74811309814453, 152.87405654907226, 76.43702827453613, 38.218514137268066, 19.109257068634033, 9.554628534317017, 4.777314267158508, 2.388657133579254, 1.194328566789627, 0.5971642833948135, 0.29858214169740677],          "extension" : "png",          "format" : "png32",          "maptype": "street",          "sensor": "false",          "tileSize": [512,512],              },     {         "baseURL":"http://mygeoserver/geoserver/datageo/wms",         "opacity":0.5,         "singleTile":false,         "type":"WMS",         "layers":["BrasilUF"],         "format":"image/png",         "styles":[],         "overview":false     }         ]}


title: Arcmap [ArcView License only] 10 labels at different scales

text: Is it possible to have different labels at different scales (without duplicating layers)? For e.g. I have loads of point shapefiles and would like to label them with ID, Name and postcode. My idea was to give more information as we zoom in - meaning 1) at full extent - no labels 2) at 1: 1,000,000 - ID only 3) at 1: 500,000 - ID and Name 4) at 1: 100,000 - ID, Name and Postcode Any ideas or thoughts?


title: How can I get a single part dissolve using ST_Union or ST_Aggr_Union?

text: I am trying to get dissolve functionality using `ST_Union` using this query :               SELECT c.fid, ST_Union(c.boundaryshape) FROM c GROUP BY c.fid;       Than I found instead of using `ST_Union` I should use `ST_Aggr_Union` :               SELECT c.fid, ST_Aggr_Union(c.boundaryshape) FROM c GROUP BY c.fid;      I am getting results like dissolve but it looks like this is multi-part. Can somebody let me know how can I get the single-part dissolve?


title: ArcGIS Server Manager (10.2.2) - permissions applied to folder not "sticking"

text: I was able to do this before updating to 10.2.2, because I have several other folders with specific permissions applied via local role store. User store is Active Directory. Data store is flat files via UNC to shared drive. Each time I edit the permissions of a newly-created folder, I choose private, un-check allow access to all users, & add the role to allowed roles. Saving says that, "permissions applied successfully," but the lock icon remains open, & checking the permissions shows that the folder is public. I have tried deleting & re-creating the folder in the server manager.


title: Understanding URL that IdentityManager.findServerInfo() method expects?

text: Does anyone have experience working with the IdentityManager object in the JS API 3.9? I'm trying to understand what URL the findServerInfo() method expects. I've tried the base REST url, and individual service URLs, but so far no luck.


title: OpenLayers code for not showing intended area from Google Streets

text:
Would anyone be able to help me with the following code? I have downloaded googleclient.htm world map and want to change it to display Australia only. So I have changed the max.Extent:new OpenLayers Bounds(-20037508, -20037508, 20037508, 20037508) from displaying the whole world to the coordinates I received when I zoomed in to Australia (12249491, 15302081, 17337140, 14910723). For some reason, I'm not getting Australia, and I don't understand what I'm doing wrong.                                  <title>Webmaps</title>

This is my map:


title: Is there a Targer-like tool for QGis?

text: I'm trying to make subsurface drillhole and I heard that Target (an AutoCad tool) can help me with this task but the issue is that I work with QGis.


title: QGis & GeoServer WPS error

text: When I try to use QGis WPS with GeoServer, I get the following error in QGis :               java.lang.IllegalArgumentException: argument type mismatch      Here are the GeoServer logs :               2013-09-19 15:18:50,934 WARN [geoserver.ows] - Could not get a ServiceInfo for service wps thus could not check if the service is enabled         2013-09-19 15:18:52,809 INFO [geoserver.wps] -          Request: getCapabilities         service = wps         baseUrl = http://localhost:8080/geoserver/     2013-09-19 15:18:58,121 WARN [geoserver.ows] - Could not get a ServiceInfo for service wps thus could not check if the service is enabled     2013-09-19 15:18:58,121 INFO [geoserver.wps] -      Request: describeProcess         service = wps         version = 1.0.0         baseUrl = http://localhost:8080/geoserver/         identifier[0]:             value = gs:HelloWPS     2013-09-19 15:19:00,824 WARN [geoserver.ows] - Could not get a ServiceInfo for service wps thus could not check if the service is enabled     2013-09-19 15:19:00,918 ERROR [geoserver.ows] -      java.lang.IllegalArgumentException: argument type mismatch      I do have a wps.xml file in geoserver2.4-RC2\data_dir :               <wps>       <id>WPSInfoImpl-15378ab6:140ed6a544e:-8000</id>       <enabled>true</enabled>       <versions>         <org.geotools.util.Version>           <version>1.0.0</version>         </org.geotools.util.Version>       </versions>       <citeCompliant>false</citeCompliant>       <schemaBaseURL>http://schemas.opengis.net</schemaBaseURL>       <verbose>false</verbose>       <connectionTimeout>30.0</connectionTimeout>       <resourceExpirationTimeout>300</resourceExpirationTimeout>       <maxSynchronousProcesses>2</maxSynchronousProcesses>       <maxAsynchronousProcesses>2</maxAsynchronousProcesses>       <processGroups>         <processGroup>           <factoryClass>org.geoserver.wps.DeprecatedProcessFactory</factoryClass>           <enabled>true</enabled>         </processGroup>         <processGroup>           <factoryClass>org.geoserver.wps.jts.SpringBeanProcessFactory</factoryClass>           <enabled>true</enabled>         </processGroup>         <processGroup>           <factoryClass>org.geotools.process.geometry.GeometryProcessFactory</factoryClass>           <enabled>true</enabled>         </processGroup>         <processGroup>           <factoryClass>org.geotools.process.raster.RasterProcessFactory</factoryClass>           <enabled>true</enabled>         </processGroup>         <processGroup>           <factoryClass>org.geotools.process.vector.VectorProcessFactory</factoryClass>           <enabled>true</enabled>         </processGroup>       </processGroups>     </wps>      Anyone knows where the problem comes from ? Thanks * * * Running on Win XP


title: How can I import data from Excel to SAGA?

text: I have just started using SAGA and I want to plot some geographic coordinates (latitudes and longitudes) and I want to import them from Excel. Please tell me what to do from Excel to when the points are plotted. Thank you.


title: How to Extract Data by Attributes (Clip/Zip/Ship) using Attribute Where Clause - Esri Server Tools

text: I want to use the out-of-the-box Esri geoprocessing task to extract data via ArcGIS Server (Server Tools --> Extract Data). Its pretty straight-forward, but it doesn't seem to support use of a "where clause" as an input parameter for the service. I could add this in using python, but would love if it just supported it without custom development. Has anybody seen a way to use ArcGIS Server to extract data using an attribute where clause WITHOUT custom development. Thanks ahead of time.


title: trsp route between two points ( not edge to edge )

text: I'm trying to solve common problem, how to generate route/navigate between two GPS coordinates and not between two edges using trsp algorithm. I did read many articles here but i'm not sure how should i create temporary edge on huge psql database where even now routing taking several seconds and insert/update can take forewer. Solution where i will navigate to nearest edge and then calculate distance between edge coordinates and destination coordinates is in my case not good b/c i'm trying to avoid one way streets and this solution is ignoring it. Another solution i have is to find edge where both start(x1,x2) and end(x2,y2) coordinates are nearest to the destination GPS coordinates ( which does not need to be nearest edge ) and navigate there but found this edge is also very slow.


title: Which bins should be used for demographic data for a choropleth map?

text: I'm currently working on a project which produces a choropleth map similar to a map like this: http://dl.dropboxusercontent.com/u/54512470/Telefondaten_VIZ/index.html That map has 2x5 classifications (5 for positive migration numbers, 5 for negative ones.) I have different underlying numbers than the one I linked. How do I decide the following things:   * Linear or Exponential increase of the break points (10, 20, 30, 40, 50 vs. 10, 20, 40, 80, 160)? Or even square-root-shaped?   * related: If linear(ish), when is it ok to make some bins larger than others?   * 'final break point', i.e. what would be the equivalent of -500 and 500 in the linked example Is there generic advice on these things which can be given without looking at the data? **Edit:** The reason I don't want to post the data here is that I'd like to use some kind of algorithm to automatically create bins and resulting choropleth maps. Unlike the linked project, I don't have the migration data for a fixed timespan, but I have data for 2007, 2008, and so on up to 2011. I'd like a radio button to switch between the years (can combine several years). I don't want the colors to be consitent between different different settings (i.e. + 5000 in year 2007 may have a different color than for 2010 + 2011), but to be informative for a fixed year. One approach I was thinking of is: Take the largest absolute migration number (e.g. 1200) and round down to the nearest 50, 100, 500, 1000 and so on. This is the lower border of the dark orange, let's call this threshold number T (i.e. 1000 in this example). The threshold for the dark blue is everything below -T (i.e. -1000). The binning points of the linked project are T/10, T/5 and T/2 (in the example 100, 200, 500). I'm unsure about these particular break points. Why not T/10, T/4, T/2 (100, 250, 500)?


title: How to convert 'RasterBrick' object to 'im' object in R?

text: During script execution I need to convert 'RasterBrick' to 'im' object. Package 'geostatsp' is supposed to do this job:               library(geostatsp)     myraster = raster(matrix(1:100,10,10), xmn=0,xmx=10,ymn=0,ymx=10)     class(myraster)     >[1] "RasterLayer"     >attr(,"package")     >[1] "raster"          myIm = as.im(myraster)     class(myIm)     >[1] "im"      But when I try it with real data I catch an error:               library(geostatsp)          b <- brick(some_file)     r <- raster(b)     class(r)     >[1] "RasterLayer"     >attr(,"package")     >[1] "raster"          bim <- as.im(r)     >Error in .local(x, ...) : 'x' has no values


title: Get current route segment from Route type using ArcGIS Runtime LocationListener?

text: I am trying to get the route segment on which I am currently navigating. i.e. I want to get the current route segment when `onLocationChanged` method of my `LocationListener` is called. A route is constituted from small segments. And we have the routing directions related to a segment but I do not know how to get the current segment in my `LocationListener`.


title: Cleaning up overlapping polylines

text: I am working on creating a project boundary for a specific area. To do this I merged several data sets including roads,rivers & railroads (polylines) and lakes, civil townships & congressional townships, (polygons converted to polylines) into a single polyline feature class. I then split the polylines at every intersect in the data. I then created a subset of the data using only the nearest features to the project area to create a boundary. My problem is my data has some issues with overlap. In some areas the civil townships are coincident with bodies of water or roads and the polyline for both types of feature are overlapping (or close to overlapping) in my subset of the data. I would like to clean up the data so there is only 1 line for each section of the boundary and duplicate boundary lines are deleted. after I get a clean polyline feature I was planning and using the feature to polygon tool to create the the final boundary. Is there a simple method of cleaning up this data? I already have used topology rules to clean up a few thing such as dangles, but I couldn't really find a rule that applied to my situation. any advice would be appreciated.


title: Using update cursor in a while loop to perform analysis on individual polygons gets stuck on first row

text: I am working through an open courseware project from PSU found at https://www.e-education.psu.edu/geog485/node/145 I am trying to use an update cursor in a while loop to perform an analysis on individual polygons within a multipart feature. However, the output is for the first row ad infinitum. Here is what i have so far:               import arcpy     arcpy.env.workspace = "c:/users/documents/arcgis/arcgis tutorials/psu tutorials/project3/policedata.gdb"     arcpy.env.overwriteOutput = True          # Declare variables for files and fields used in analysis          incidents = "c:/users/documents/arcgis/arcgis tutorials/psu tutorials/project3/policedata.gdb/GraffitiIncidents"     patrol_zones = "c:/users/documents/arcgis/arcgis tutorials/psu tutorials/project3/policedata.gdb/PatrolZones"     incidents_layer = "IncidentsLayer"     patrol_field = "NAME"          cursor = arcpy.UpdateCursor(patrol_zones)     zone = cursor.next()          # Loop through each patrol zone          while zone:         try:             arcpy.MakeFeatureLayer_management(incidents,"IncidentsLayer" )             zone_name = zone.getValue(patrol_field)             print zone_name                  sql = patrol_field + " = " + "'" + zone_name + "'"              # Make a feature layer for each zone             arcpy.MakeFeatureLayer_management(patrol_zones, "CurrentZone", sql)                   # Make selection             arcpy.SelectLayerByLocation_management(incidents_layer, "CONTAINED_BY", "CurrentZone")             num_incidents = arcpy.GetCount_management(incidents_layer)                       print "Graffitti incidents = " + str(num_incidents)     ##        arcpy.Delete_management("CurrentZone")     ##        zone         except:             print "Could not perform analysis."               del cursor      What in my code is causing the script to produce results only for the first row?


title: How to fix QGIS crashes when using GRASS toolbox?

text: I've installed QGIS 2.0.1 on Ubuntu 12.04. exactly as explained on qgis.org . Unfortunately QGIS crashes when trying to use any GRASS model (r.aspect as example): just opening QGIS, click on toolbox expand toolbox, click on GRASS function r.aspect (as example) and it hangs... no error message, no possibility to define data for input! It's probably a local problem but any help is appreciated. I've had QGIS 1.7.4 installed prior. Additionally I tried on a fresh ubuntu build (12.04) in a virtualbox. same results: when trying to us a function from Grass-toolbox it crashes QGIS :-(


title: How to repair a data source with joins to Excel?

text: I have a set of MXD files, that do not use relative paths, and related data that have been copied to a new computer. Normally it's not a big problem to repair the data sources in the project, but in this case some of the layers are joined to an Excel file. When I try to repair the data sources any layer that is joined to an Excel file looses it's symbolization. Has anyone run into this problem and figured out a fix?


title: How do I get the current Undo and Redo item in ArcMap Operation Stack?

text: I just posted this on ESRI forum and found this link so I am trying it here too: Research ArcGIS Resource Center document - How to work with the operation stack: http://help.arcgis.com/en/sdk/10.0/a...00000288000000 I found out the following: "The Undo and Redo methods do not remove operations from the operation stack, but rather move an internal pointer up and down the stack." Question: So how do I get access to the operation that ArcMap is getting ready to Undo or Redo? History and sample code: I have been working on trying to manage ArcMap operation stack from a custom form. In the custom form we are actively editing data using ArcMap controls (add a sketch, modify feature, etc..) and we are also calling custom edits which are enclosed within StartOperation and StopOperation("EditOperationName"). In my code, I would like to evaluate the actual operation that ArcMap is getting ready to Undo or Redo and then perform a specified action if the Operation is a specified custom edit operation that was added within the edit session. I am able to return the name of all the menu string edit operations in the IOperationStack, I am able to call the mxDoc.OperationStack.Undo and Redo methods but I have not found a way to evaluate the actual Undo or Redo operation that ArcMap is pointing to. Refer to the code sample below for the Undo button click:               Private Sub btnUndoLastTest_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles btnUndoLast.Click     'Dim pOperationStack As IOperationStack     'pOperationStack = _mxDocument.OperationStack          If _mxDocument.OperationStack.Count > 0 Then     Debug.Print("Before Undo is Selected")     For i = 0 To _mxDocument.OperationStack.Count - 1     Dim pOper As IOperation = _mxDocument.OperationStack.Item(i)     Debug.Print("operation number: " + i.ToString() + " Name: " + pOper.MenuString)          If pOper.CanUndo Then     Debug.Print("CanUndo Operation")     Else     Debug.Print("Can NOT Undo Operation")     End If     If pOper.CanRedo Then     Debug.Print("CanRedo Operation")     Else     Debug.Print("Can NOT Redo Operation")     End If     Next          '_editor.UndoOperation()     Dim pOperation As IOperation = _mxDocument.OperationStack.Item(_mxDocument.Operat ionStack.Count - 1)     Debug.Print("Operation Name = " + pOperation.MenuString.ToString())          If pOperation.MenuString = "Apply Rules" Then     'reset the IsApplyRules variable     _geoeditWorflow.IsApplyRulesSet = False     End If          pOperation = Nothing          _mxDocument.OperationStack.Undo()          Debug.Print("After Undo is Selected")     For i = 0 To _mxDocument.OperationStack.Count - 1     Dim pOper As IOperation = _mxDocument.OperationStack.Item(i)     Debug.Print("operation number: " + i.ToString() + " Name: " + pOper.MenuString)          If pOper.CanUndo Then     Debug.Print("CanUndo Operation")     Else     Debug.Print("Can NOT Undo Operation")     End If     If pOper.CanRedo Then     Debug.Print("CanRedo Operation")     Else     Debug.Print("Can NOT Redo Operation")     End If     Next     End If     End Sub      NOTE: 1) It is interesting to see that ArcMap acutually sets ALL Operations in its stack to either true or false for CanUndo or CanRedo properties regardless of wether this was the last operation that you applied undo or redo. 2) ArcMap keeps an Internal pointer to the stack and moves this pointer up and down. So, this means that the OperationStack is NOT reorganized (or even a true STACK.POP) based on an Undo and Redo. Therefore it doesn't appear like there is an easy way to get the operation Item that ArcMap is getting ready to undo or redo. You can see in ArcMap tools what the undo or redo operation name is by either hovering over the undo and redo buttons or in the standard Menu toolbar --> Edit -> Undo "Name" or Redo "Name" Any Ideas on how I can get the name of the NEXT Undo or Redo operation? I am trying to avoid having to write my own UndoStack and RedoStack if at all possible. Thanks for any help you can provide.


title: Upload Shape File using GeoServer REST API and C#

text: I'd like to upload shape files using GeoServer REST API and C#. So far, I can create workspaces using the REST API and C#. From GeoServer documentation (specific to cURL), the difference between creating a workspace and uploading a shape file is in the HTTP Method (workspace creation uses `POST`, shape file uploading uses `PUT`), ContentType (workspace creation uses `application/xml`, shape file uploading uses `application/zip`), and finally there is the `--data-binary` parameter in shape file uploading. The code I use for creating workspace is shown below. Can anyone kindly help me refactor this code so that I can use it to upload a shape file. Thanks in advance.                       String gUrl = "http://localhost:8081/geoserver/rest/workspaces";             WebRequest request = WebRequest.Create(gUrl);                  request.ContentType = "text/xml";             request.Method = "POST";             request.Credentials = new NetworkCredential("geoServer-username", "password");                  byte[] buffer = Encoding.GetEncoding("UTF-8").GetBytes("<workspace><name>new_workspace_csharp_2</name></workspace>");             Stream requestStream = request.GetRequestStream();             requestStream.Write(buffer, 0, buffer.Length);             requestStream.Close();                  WebResponse response = request.GetResponse();             Console.Write("Response from GeoServer: " + response);      EDITED: I've managed to make some headway, and I can now upload a shape file using C# and the GeoServer REST API. My refactoring effort, which now works, is shown below:                   public bool UploadShapeFile(string workspace, string dsName, Uri zipUri)         {             String fileUri = zipUri.AbsolutePath;             Console.Write(fileUri);                  byte[] localShapeFile = readLocalShapeFile(fileUri);                  String sUrl = "http://localhost:8081/geoserver/rest/workspaces/cs_ws1/datastores/cs_ds1/file.shp";             WebRequest request = WebRequest.Create(sUrl);                  request.ContentType = "application/zip";             request.Method = "PUT";             request.Credentials = new NetworkCredential("gs-username", "password");                  // byte[] buffer = Encoding.GetEncoding("UTF-8").GetBytes(@fileUri);             Stream requestStream = request.GetRequestStream();             requestStream.Write(localShapeFile, 0, localShapeFile.Length);             requestStream.Close();                  WebResponse response = request.GetResponse();             Console.Write("Response from GeoServer: " + response);                       return false;         }              private byte[] readLocalShapeFile(string filePath)         {             byte[] buffer;             FileStream fStream = new FileStream(filePath, FileMode.Open, FileAccess.Read);             try {                 int length = (int)fStream.Length;                 buffer = new byte[length];                  int count;                  int sum = 0;                       // Read until Read method returns 0 - End of stream reached                 while ((count = fStream.Read(buffer, sum, length - sum)) > 0)                     sum += count;             }             finally {                 fStream.Close();             }                  return buffer;         }      However, I still believe there can be a more elegant way of doing this. Any suggestions will be highly appreciated.


title: Problem displaying OpenLayers KML layer in IE

text:
I am serving the following files from a the same directory on localhost:               map.html     OpenLayers.js     output.kml      _map.html_ simply creates an OpenLayers map and displays the KML file as per the OpenLayers example. I access it via _http://localhost:8080/map.html_ and it all works fine in Firefox and Chrome but not in IE. In IE I get the navigation buttons but no map. The OpenLayers example does work in IE... The contents of _map.html_ are shown below:                         <title>Map</title>

                     Any ideas?


title: GeoServer container need to test locally?

text: GeoServer container need to test locally I did some tests with geoserver and openlayers but it did not work with the layers I have on my machine. I did not need the container to use locally calling my WMS layers using openlayers. My javascript pages are in the openlayers folder on my C :/ are not in JSF project. My question is I need the container to use locally in my tests?                   <!DOCTYPE html>     <html lang='en'>     <head>      <meta charset='utf-8' />      <title>My OpenLayers Map</title>      <script type='text/javascript' src='OpenLayers.js'></script>           <script src="http://maps.google.com/maps/api/js?sensor=false&v=3.2"></script>           <script type='text/javascript'>                var map;      var mercatorProjection = new OpenLayers.Projection('EPSG:900913');      var latLongProjection = new OpenLayers.Projection('EPSG:4326');      var bounds = new OpenLayers.Bounds(                     -51.2809219360352, -30.2445888519287,                     -51.0207977294922, -29.9661273956299                 );           function init() {            map = new OpenLayers.Map('map_element',{         maxExtent: bounds,         maxResolution:  0.0010877400636672,         units: 'm',         //allOverlays: true,         projection: mercatorProjection,         displayProjection: latLongProjection      });                var google_streets = new OpenLayers.Layer.Google(         "Ruas",         {numZoomLevels: 20}         );          var wms_layer = new OpenLayers.Layer.WMS(     'Linhas Ônibus',     'http://localhost:8080/geoserver/Teste/wms',     {layers: 'Teste:onibus', transparent: true},     {isBaseLayer: false,         opacity: 0.7}     );               //Adiciona as camadas ao mapa     map.addLayers([google_streets, wms_layer]);               var point = new OpenLayers.LonLat(-51.22,-30.08);      point.transform(new OpenLayers.Projection("EPSG:4326"),          map.getProjectionObject());      map.setCenter(point, 10);       //Camada de controle que vai mostrar as camadas no mapa      map.addControl(new OpenLayers.Control.LayerSwitcher({}));           //Verifica se o mapa tem um ponto central e o extende a sua extensão máxima      if(!map.getCenter()){      map.zoomToMaxExtent();           }      }           </script>     </head>          <body onload='init();'> <!-- Chama a função js init() -->           <!--Elemento HTML onde o mapa é exibido -->      <div id='map_element' style='width: 800px; height: 800px;'>      </div>          </body>     </html>


title: OpenLayers PanZoomBar css style

text: I need to use PanZoomBar in my application.               var panZoomBar = new OpenLayers.Control.PanZoomBar();     map.addControl(panZoomBar);      Also I need to overwrite the PanZoomBar css style. I overwrite the whole bar.               .olControlPanZoomBar {         top: 38px !important;         left: 18px !important;         width: 44px;         height: 214px;         background: transparent url(../images/zoombar_bg.png) no-repeat 50% 100%;     }      But how I can get access and overwrite the default buttons? I mean:               .whatName(Id)ShowIWriteHere {         some style     }


title: Check validity of file geodatabase using arcpy

text: AS the title suggests, how can I whether a GDB file is valid using arcpy? One approach is, only checking the whether the name ends with `.gdb`, but that's not a robust approach. Because a normal file can be ended with `.gdb`, where GDBs are essentially folders to OS. So how can I do that?


title: Extract information of a raster overlayed on another raster with a lower resolution directly in python

text: I have two aligned rasters with the same extent, one has a resolution of 5 meters and the other of 125 meters. I would like to place the 5 meter one over the 125 meter one and assign to these smaller pixels the value of the large pixel on which they land on. How can I do this in pure python code? The raster I/O is not a problem.


title: Is there mathmatical function/equation available for ArcGIS Network Analyst's Closest Facility Method?

text: I would like to use this function in other statistical programming, ideally SAS where td = (whatever the function is to calculate total distance from centroid of county to facility). The goal is to leverage other demographic weights in drive time and distance.


title: How does ESRI ArcGIS 10 recognize a WMS CRS?

text: I'm working on a WMS service with Mars data. Currently the GetCapabilities shows `<SRS>EPSG:4326</SRS>`. I assume ESRI ArcGIS 10 reads this and therefore the WMS has GCS_WGS_1984 as a Geographic Coordinate System? I would like to use a custom URL definition because I am working with data from Mars. Is it for example possible to add a URL with GML such as http://spatialreference.org/ref/iau2000/49900/gml/ ? Thanks in advance, Jelmer Oosthoek


title: Using gdal to convert from UTM to equidistant lat/lon

text: I am getting rather confused with converting a UTM gridded dataset into a equidistant/rectangular latlon gridded dataset. Can anybody help me? I have a UTM set with the following projection:               PROJCS["WGS 84 / UTM zone 18N",     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"]],     PROJECTION["Transverse_Mercator"],     PARAMETER["latitude_of_origin",0],     PARAMETER["central_meridian",-75],     PARAMETER["scale_factor",0.9996],     PARAMETER["false_easting",500000],     PARAMETER["false_northing",0],     UNIT["metre",1,         AUTHORITY["EPSG","9001"]],     AUTHORITY["EPSG","32618"]]     Origin = (556950.000000000000000,5573990.000000000000000)     Pixel Size = (90.000000000000000,-90.000000000000000)     Metadata:      AREA_OR_POINT=Area     Image Structure Metadata:      INTERLEAVE=BAND     Corner Coordinates:     Upper Left  (  556950.000, 5573990.000) ( 74d12' 0.38"W, 50d18'54.95"N)     Lower Left  (  556950.000, 4905470.000) ( 74d17' 9.82"W, 44d18' 0.91"N)     Upper Right (  655950.000, 5573990.000) ( 72d48'37.20"W, 50d17'50.61"N)     Lower Right (  655950.000, 4905470.000) ( 73d 2'43.67"W, 44d17' 8.74"N)     Center      (  606450.000, 5239730.000) ( 73d35'30.62"W, 47d18' 8.60"N)     Band 1 Block=1100x1 Type=Float32, ColorInterp=Gray      And want to convert this to a rectangular latlon set (with corners just as they are provided in the corner coordinates. What to do? I am sure it should be easy...The biggest problem is that if I directly convert this into latlon, it won't be equidistant... Cheers, help is appreciated.


title: Looking for best practice to connect multiple road vector line layers for routing (pgRouting) purposes

text: **The Data I have:**   * `Bridges` (line layer)   * `Municipal_roads` (line layer and does not include information in `Bridges`)   * `State_roads` (line layer that doesn't include `Municipal_roads` or `Bridges`) **Question:** Since these layers are continuously getting updated/upgraded especially in construction season; what is the _'best of practice'_ way to connect these layers to produce one layer for routing purposes? **Edit:** I would like to, if possible avoid merging the layers.


title: Delete identical by size?

text: I got a Feature class with some identicial fields! Now i only need one polygon which had the same ID and this is the one with the largest size? Is there a posibility to delete the Polygons which are identical by size so that only the largest Polygons(from the ones with the same ID) are left? For example:               ID   shape area          1    5     2    4     2    3     3    11     4    7     5    8     5    10      and i Need only the Polygons with the largest size that means the frist, second, fourth, fifth and the last!


title: How to interpolate GPS data correctly?

text: I'm working on a project that involves using linear interpolation to fit lines to lat/lng points. I'm pretty new at scientific programming related stuff, and I was wondering if people could steer me in the right direction.What is a good tutorial for linear interpolation in Python, preferably using GIS data?


title: Search cursor stops selecting after two iterations

text: I have a script which selects an extent from a grid and exports the data frame.               import os, arcpy          mxd = arcpy.mapping.MapDocument('CURRENT')     df = arcpy.mapping.ListDataFrames(mxd,"Layers")[0]     grid_lyr = arcpy.mapping.ListLayers(mxd,"grid")[0]     scur = arcpy.SearchCursor(grid_lyr)          for row in scur:         sqlExp = """{0} = '{1}'""".format(arcpy.AddFieldDelimiters(grid_lyr,"PageName"),row.getValue("PageName"))         arcpy.SelectLayerByAttribute_management(grid_lyr,"NEW_SELECTION",sqlExp)         print "Selecting {}".format(row.getValue("PageName"))         df.extent = grid_lyr.getSelectedExtent(False)         df.zoomToSelectedFeatures()         df.scale = df.scale*1.05         arcpy.mapping.ExportToTIFF(mxd,os.path.join(r'C:\temp\109147',row.getValue("PageName")+".tif"),df,817, 534,world_file=True)              arcpy.RefreshActiveView()     del scur, mxd      I ran this script last week, and it crashed ArcMap after exporting the first tiff. I then ran it in a blank ArcMap document, changed `CURRENT` to the full path and it worked (all 414 tiff files exported). I am running the same script again today, on the same data, from a blank document. It now successfully exports the first two tiffs, and then stops executing as if it is finished, without any errors. When I try to run it again, ArcMap crashes again. I am not sure what the problem is, last week when it failed I just opened a blank document, as that normally solves the issues half the time. I updated the code to include print statements for debugging:               for row in scur:         print "Row count: " +str(arcpy.GetCount_management(grid_lyr))         sqlExp = """{0} = '{1}'""".format(arcpy.AddFieldDelimiters(grid_lyr,"PageName"),row.getValue("PageName"))         arcpy.SelectLayerByAttribute_management(grid_lyr,"NEW_SELECTION",sqlExp)         print "SQL is: {}".format(sqlExp)         df.extent = grid_lyr.getSelectedExtent(False)         df.zoomToSelectedFeatures()         df.scale = df.scale*1.05         print "Data frame {0} has a scale of {1}".format(df.name, df.scale)         name = os.path.join(r'C:\temp\109147',row.getValue("PageName")+".tif")         arcpy.mapping.ExportToTIFF(mxd,name,df,817, 534,world_file=True)         print "Exported {}".format(name)          arcpy.RefreshActiveView()          del scur, mxd, row      The output was:               Row count: 500     SQL is: "PageName" = 'E3'     Data frame Layers has a scale of 2276.80888604     Exported C:\temp\109147\E3.tif     Row count: 1     SQL is: "PageName" = 'E4'     Data frame Layers has a scale of 2276.80888604     Exported C:\temp\109147\E4.tif      My grid layer only has 414 features, so the row count numbers are a bit weird. Also, if I try to define the cursor variable again in the same session, ArcMap crashes.


title: how to generate a NDVI in QGIS

text: i have tried using raster calculator to generate NDVI and tried searching for but couldn't find it. is there any other way to extract NDVI and what are the procedures to generate NDVI in QGIS.


title: Crop and reproject shapefile based on GeoTIFF sample

text: I have GeoTIFF with these attributes:               Size is 495, 309     Coordinate System is:     PROJCS["unnamed",         GEOGCS["unnamed ellipse",             DATUM["unknown",                 SPHEROID["unnamed",6370000,0]],             PRIMEM["Greenwich",0],             UNIT["degree",0.0174532925199433]],         PROJECTION["Lambert_Conformal_Conic_1SP"],         PARAMETER["latitude_of_origin",47.5],         PARAMETER["central_meridian",4],         PARAMETER["scale_factor",1],         PARAMETER["false_easting",0],         PARAMETER["false_northing",0],         UNIT["metre",1,             AUTHORITY["EPSG","9001"]]]     PROJ.4 string is:     '+proj=lcc +lat_1=47.5 +lat_0=47.5 +lon_0=4 +k_0=1 +x_0=0 +y_0=0 +a=6370000 +b=6370000 +units=m +no_defs '     Origin = (-2963997.870569999800000,-1848004.200800000000000)     Pixel Size = (11975.754948303031000,11961.191901488673000)     Image Structure Metadata:       COMPRESSION=LZW       INTERLEAVE=BAND     Corner Coordinates:     Upper Left  (-2963997.871,-1848004.201) ( 24d36'23.04"W, 26d22' 5.88"N)     Lower Left  (-2963997.871, 1848004.097) ( 45d39'38.75"W, 55d16'35.82"N)     Upper Right ( 2964000.829,-1848004.201) ( 32d36'23.13"E, 26d22' 5.85"N)     Lower Right ( 2964000.829, 1848004.097) ( 53d39'38.88"E, 55d16'35.76"N)     Center      (   1.4791350,  -0.0520200) (  4d 0' 0.07"E, 47d30' 0.00"N)     Band 1 Block=495x4 Type=Float32, ColorInterp=Gray       NoData Value=9.969209968386869e+036      and I want to overlay Natural Earth shapefile. So I have to crop and reproject global Natural Earth shapefile to this GeoTIFF sample, however I can't understand how to do this with such projection and was wondering if I can instruct somehow `ogr2ogr` to use my GeoTIFF as reference to crop and reproject the shapefile? Also I'm interested in other approaches if `ogr2ogr` does not have this feature.


title: Leaflet HTTP.Get bad request; not loading tiles

text: I am trying to walk through the Leaflet quick start tutorial and am having a really simple issue: None of my tiles are loading. I get this error at the image url and an HTTP Bad Request in the inspector: > {"message":"id must be a string with the format account.handle"} My code is the same as the tutorial:               <link rel="stylesheet" href="//cdn.leafletjs.com/leaflet-0.7.3/leaflet.css" />


title: Deleting unconnected lines with QGIS

text: I am using a road map in shapefile format in order to calculate a route between two points. The map is downloaded from OpenStreetMap and it has small unconnected lines that are messing up my routing. Is there any way that I could get rid of those without having to manually remove them one by one? The lines look like this: http://screencast.com/t/4jnQFx57


title: Optimum grid/pixel size for interpolating rainfall and temperature

text: I have a set monthly rainfall and temperature data within an approximately 100 x 100 km study area. The sampling points are unevenly distributed within the study area. Based on a paper by Janis et. al, the number of samples is more than enough. I plan to create a raster layer by interpolating the data using either Inverse Distance Weighting or Spline techniques. I can't decide on the best grid/pixel size to use. Is there a certain tull of thumb given an `x` number of samples and a `y` size of the study area to determine the best grid size to use?


title: How to send data from my application to QGIS?

text: I have an application that makes analysis on maps . After making analysis I have new map in my application. My application doesnt use QGis API. How can I send results back? Should I construct object QgsVectorLayer from my application and send new maplayer to QGis? What approach is better to use? I want to create a plugin for QGis to work with my application and use TCP IP protocol. My application will be SERVER and QGis plugin CLIENT. My application works with graphs. It accepts .mif files. Could I construct QgsVectorMaplayer from points, lines and polygons?


title: Export and edit annotation featureclass using arcpy

text: I would like to convert labels to annotations and make some advanced editing for these labels using ArcPy. How would I do this? I'm using a file geodatabase and ArcMap 10. Any help, please.


title: My house address isn't listed on Google Maps

text: When I search for my address on Google Maps it cannot be found which means I can't order takeaway etc. Is there a way of adding my address to Google Maps or even add a street address? I want the address to be publicly available so others can search it.