title: Finding points on an odd map projection using interpolation between known points

text: I have an image of a somewhat idiosyncratic projection that is a frame from a film: ![My image](http://blog.nuclearsecrecy.com/wp-content/uploads/2013/12/WAR- DANCE.png) A geographer friend of mine said it was probably some kind of Gnomic projection, but that from the resolution it wasn't possible to know for sure. I wouldn't know. You can see North America at top left, and the Soviet Union and China (!) in red with a black border around it in the middle-right of the frame. The dark vertical and lighter horizontal lines visible seem to go through the North Pole. As you can guess this is a declassified target map from the Cold War; those branched lines are bomber routes (some from bases, some from refueling points). I am interested in finding the rough lat/lngs of the ends of the branched lines (the "targets"), knowing, of course, there is going to be a large margin of error because of the resolution,. I could hunt and peck these by hand, working first to identify targets that ought to be on there (e.g. Moscow), but it seems like there ought to be an easier way. There are plenty of identifiable geographical locations on there, such as James Bay in Canada, the tip of Gujarat in India, the Caspian Sea, and so on. My first thought was to try line between a number of these by hand in Adobe Illustrator, then draw similar lines in Google Earth, and use them as a rough grid to triangulate the missing points by hand. This is of course difficult. I then thought maybe there was some way to do this algorithmically — if I know the lat/lng of fixed pixel locations, can I move from there to figuring out the missing lat/lngs of other pixel locations? I suspect this might be harder than it sounds, given that I don't really know the projection that produced those pixel locations. Is there a better way? Any thoughts would be appreciated.


title: Creating Random Points in ArcGIS 10.1

text: I would like to create random points across my study area (which is a DEM raster) using the Feature Class> Create Random Points. However, I need to select/copy/move/delete these points and re-arrange them in a line, with each point 1km away from each other. I was wondering whether there is a specific tool that I can use in ArcGIS 10.1 to achieve this -I mean a tool that can tell me how far each point is from another point or is there way to measure the distance between two points to ensure I am placing these points exactly 1km away. Any info would be appreciated!


title: How to use Maperitive with the two layers HikeBikeMap (Hillshading) and OSM CycleMap without the red and blue routes?

text: I want to use Maperitive for creating maps of National Parks in Europe. I choose HikeBikeMap (Hillshading) and OSM CycleMap, because it shows what I want regarding the levels of hights, but there are already made many routes in red and blue, but I don't want these, but to mark my own routes/trails. Are there any possibility for getting rid of the red and blue routes? Thanks in advance Lone


title: How to create fixed length (and width) buffer around a polyline in ArcMap 10?

text: Could someone please help me? I would need to create 100 buffers (500m x 1km) along a stream polylines. I know how to do it by using shortest distance but the buffers should follow the curves of the streams. Each buffer should start from certain point (sampling site) and go to upstream for a 500m. Some of my stydy-streams have multiple tributaries on upstream areas and in those cases the buffer should follow each branch. So far I have a polyline-layer with 100 streams and point-layer with 100 sampling points snapped together with the end of study-streams. How to proceed from here? With some network tools? I would really appreciate some help! ![Schematic figure describing my problem](http://i.stack.imgur.com/dLR8T.jpg)


title: How to create a QGIS road style for OSM data?

text: I am trying to use qgis to visualize the shortest path algorithm's result. I managed to import the osm map into qgis ![enter image description here](http://i.stack.imgur.com/ISt8F.png) But I think it looks kind of different. I was referring to this link and she has this nice map like this ![enter image description here](http://i.stack.imgur.com/fFqul.png) So I tried to import the database creating using osm2po. And this is what I got ![enter image description here](http://i.stack.imgur.com/pOWTE.png) This is still different from the one seen at link. I know the region is different, but still I followed the same step as defined in link link. I also used the same file for austria and imported it into my database. I changed the style as well and this is what I got. ![enter image description here](http://i.stack.imgur.com/O8xCf.png) Now the roads and highways in the link link look nice and fine. However, in my case they look kind of overlapped. Why is it so?


title: Generate a global view Image of where a Lat-Long location is

text: I would like to automatically generate an image which highlights a location from a global view. Here is an example of what I would like to achieve: ![enter image description here](http://i.stack.imgur.com/OOdKv.png) I am debating using python utilities such as mpl_toolkits.basemap,which examples can be seen here. But I was hoping to find an easier way from you GIS professionals. This is for a small webpage, and it would be nice to have it automatically generated. Thank you for your time, --Beau


title: ArcGIS exporting files to Access 2007

text: Does anyone know how to export a shape file out of ArcGIS into an Access 2007 database?


title: Adding multiple layers from an array in OpenLayers?

text: I have a lot of OpenLayers layers that change every day so instead of having to create, and change a long list of these:               var layerone = new OpenLayers.Layer.XYZ("mylayer1", "mylayer1/${z}/${x}/${y}.png", { });          var layercool = new OpenLayers.Layer.XYZ("mylayer1", "mylayer1/${z}/${x}/${y}.png", { });      etc I'd like to first make an array of layer names:               var mylayers = new Array("layerone","layercool","layeryo","layerhello");      and then create all these layers in a for loop, something like this:               var length = mylayers.length;     for (var i = 0; i < length; i++) {          mylayers[i] = new OpenLayers.Layer.XYZ("mylayers[i]", "mylayers[i]/${z}/${x}/${y}.png", { });           }      How do I first of all create unique var names for each layer in the foor loop and then do the loop? I read somewhere to use `window[layers[i]] = new OpenLayers.Layer...` or maybe `eval`? but I'm not sure how to access these later on when doing the `map.addLayers();`


title: Google Map API Drag event..?

text:
The below code is working perfectly on click event to retrieve Lat/Long and Current Zoom level. How do I get the same when I hover mouse/drag on map Window (not on submit button)?                   <!DOCTYPE html>                                                  <input type= "text" ID=latlng>     <input type="text" ID=Zoom>
     <button onclick="copyText()">Submit</button>


title: Does Leaflet has option for scale dependent layers?

text: Is it possible to show layers with Leaflet depending on scaling? For example not show counties until zoom level reaches 10. Thanks.


title: ListFeatureClasses brings back more than expected

text: Why, when I issue the function:               fcList = arcpy.ListFeatureClasses(feature_type='point')      does it bring back a list of grid/rasters as well as all my point featureclasses from ArcSDE. All I want is my **point** featureclasses, this doesn't make any sense! Anyone had this happen to them? Plus it is incredibly slow! Any help much appreciated :)


title: Bing map SetCenter doesn't work

text: Hi everybody am trying to use bing map on my html page using openlayers and geoEXT. the problem is when i change the center of the bing map nothing happens, also when i try to zoom on the map on certain point like this:                      Bing.setCenter(new OpenLayers.LonLat(rec.get('x'),rec.get('y')).transform('EPSG:4326','EPSG:26191'), 9);      it doesn't zoom on this point with the x and y , all i get is the zoom with degree 9 here's the code:                 var Bing;     var apiKey = "xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx";      Bing = new OpenLayers.Map("map");                  map.addControl(new OpenLayers.Control.LayerSwitcher());                  var road = new OpenLayers.Layer.Bing({                 name: "Road",                 key: apiKey,                 type: "Road"             });             var hybrid = new OpenLayers.Layer.Bing({                 name: "Hybrid",                 key: apiKey,                 type: "AerialWithLabels"             });             var aerial = new OpenLayers.Layer.Bing({                 name: "Aerial",                 key: apiKey,                 type: "Aerial"             });                  Bing.addLayers([road, hybrid, aerial]);                  var center = new OpenLayers.LonLat(-93.27, 44.98).transform(                 'EPSG:4326', 'EPSG:3857'             );                  Bing.setCenter(new OpenLayers.LonLat(-110, 45), 3);           bingpanel =  new GeoExt.MapPanel({               region:"center",               height: 400,               width: 600,               map: Bing                     });


title: What is needed to develop an ArcGIS application? (C#?)

text: I am trying to fully understand the genius that is ArcGIS. I want to develop a piece of software that uses the engine. This application needs to import user data (GeoTIFF, etc) and allow users to put points and polygons on a map. I would think this would be VERY straightforward. I am unable to find sufficient detail on how to begin this process and was hoping someone here could help me. There are APIs, there is the ArcEngine, there is ArcObjects, also SDKS. I don't know which I need to develop my application. I would love to make a Windows Form application in C#. What would I need? ArcEngine? I am working on computers that have ArcGIS desktop installed. How can I know what else is installed (Engine, Objects, etc...)? How can I get ArcEngine if that is what I need? Thank you.


title: I Need to add other data to points in ArcMap 10.1 document

text: I'm a user (not producer) of ArcGIS documents in my work and I need some help asking experts with the following: I have some Microsoft Word documents and some pictures that describe characteristics of buildings that I have represented in ArcMap 10.1 document as a points in a shape. I need to have access to the data and pictures when I select any of these points (buildings) clicking on ArcMap project and read this info when needed. Can I add this both type of data in ArcMap or should I create a separated database and link ArcMap and this DB in any way? Should I create a PDF file joining MS Word and images to have only one file and access them? I greatly appreciate your help and suggestions Thanks in advance Marcelo


title: Can OpenLayers layers in QGIS have transparency?

text: I often would like to set the basemap back by lowering the transparency. I would love to have this option in QGIS for OpenLayers layers, such as the stamen layers, google, etc. Is this possible? I don't see the ability in the layer properties or within the raster transparency plugin.


title: How can I snap lines in PostGIS?

text: I have a few linestrings that together form a long line. The start- and endpoints of the linestrings are equal with their direct neighbour, but when tested the ST_Distance between any 2 lines is like 1.0x10^-11. This seems to be the known issue with encoding floats. They are the same but not the same. However when I update 2 neighbouring lines with their existing Astext coordinates they are encoding by the same algorithm and their distance is 0. But then they have a small distance with their outside neighbours (that are not updated this way). What I want to do is snap a line left and right with its neighbours so the ST_Distance with both neighbours is 0. I think ST_Snap only moves a geometry to another geometry. That is not what I want. I only want the start- and endpoint of a line to snap to the neighbour. The rest of the line should stay where it is. How do I do that? With the given answer I fixed both ends of the conflicting line like this:               update network set the_geom = ST_RemovePoint(the_geom, (select st_npoints(the_geom)-1 from network where gid=6909)) where gid=6909;     update network set the_geom = ST_AddPoint(the_geom, (select st_endpoint(the_geom) from network where gid=6123), (select st_npoints(the_geom) from network where gid=6909)) where gid=6909;     update network set the_geom = ST_RemovePoint(the_geom, 0) where gid=6055;     update network set the_geom = ST_AddPoint(the_geom, (select st_startpoint(the_geom) from network where gid=6909), 0) where gid=6055;


title: Displaying graph with value scale on Y axis using ArcGIS for Desktop?

text: I have to symbolize my vector layer with Bar Graph in ArcGIS for Desktop. Along with the graph i have to display value scale on y Axis. Is this possible? I am using ArcGIS 10.2 for Desktop but could not find this option any where.


title: How can I split a multi-part polygon into single polygons?

text: i want to split a multipart polygon which consists of several polygons into single polygons so I can delete each polygon without deleting them all at one time. I am using ArcGIS 10. ![screenshot](http://i.stack.imgur.com/Xwwfs.png)


title: ArcGIS script tool: Can you have an input parameter set to optional and have the tool work when an input is not provided?

text: I'm building a script tool that processes shapefiles ( GPS data collected on weekly basis). The tool is set up to take an input for each shapefile ( in this case Sewer, Storm, Sign). On some days, data for one of these assests may not be collected, so there will be no dataset to provide to the tool. I thought O.K., i'll just set the input parameter to optional, and if no dataset is provided, it'll just continue on to the next dataset. Wrong answer, play again. Is it possible to make this work? Here's a snippet of the code.               #SEWER#     sewer = arcpy.GetParameterAsText(0)     arcpy.AddField_management(sewer, "UID", "TEXT", 10, "", "", "", "NON_NULLABLE", "", "")     arcpy.AddField_management(sewer, "Phot1_Link", "TEXT", 240, "", "", "", "NON_NULLABLE", "", "")          #Calculate UID field     fieldName = "UID"     expression = "id_gen()"     codeblock = """def id_gen(size = 9, chars = string.ascii_uppercase + string.digits):     return ''.join(random.choice(chars) for x in range (size))"""     arcpy.CalculateField_management(sewer, fieldName, expression, "PYTHON", codeblock)      ![enter image description here](http://i.stack.imgur.com/H84UD.png)


title: Order members of OSM multipolygon relations

text: Im writing an importer for OSM files. Now i'm importing the relations with multipolygoins. The documentation says               The order of the relation members does not matter (but properly sorted member lists can help human editors to verify completeness).       Is there a pre defined algorithm to find the ways which form one closed polygon? Preferably in java


title: How can I reclassify an integer raster to a float raster?

text: I have been using various arcGIS tools to try and do this. Unfortunately, the reclassify tool in arcmap can only output integers. A cool workaround that I found is using the raster calculator in arcmap, as suggested by Jeff Evans in this post: http://forums.arcgis.com/threads/15396-Reclass-from-Integer-to- Floating-point * **"** _Here is an example using a 4 value integer (21, 41, 42, 43) raster that results in a float with the values of 1.0, 2.0, 3.0, 4.0: Float(Con("nlcd" == 21, 1.0, Con("nlcd" == 41, 2.0, Con("nlcd" == 42, 3.0, Con("nlcd" == 43, 4.0, 0.0))))) **_"*** This works like a charm and gives me the desired output... HOWEVER, I cannot use this process in arcpy... Because I need to do many loops of this process, I need it to be its own python script and “The Raster Calculator tool is intended for use in the ArcGIS for Desktop application only as a GP tool dialog box or in ModelBuilder. It is not intended for use in scripting and is not available in the ArcPy Spatial Analyst module”. To be honest, I am starting to give up on using arcpy for this, does anyone out there know how to do this efficiently using arcpy? I also tried some R code but it took way too long (my grids are around 1500x2400 cells). I also tried a python reclass script my colleague he built scratch, and it also takes too much time. Any ideas/suggestions? Thanks for your time!


title: How to sort features in a layer using OGR/Python?

text: How can I sort the features of an ogr layer by the values of a field? I tried naive variations of               def sortlayer(ds, layer, field):         lname = layer.GetName()         lsort = ds.ExecuteSQL(b'select * from "{}" order by {}'.format(lname, field))         ds.DeleteLayer(lname)         ds.CopyLayer(lsort, lname)      but this either crashed or freezed. Bonus points for a solution that doesn't change the order of layers in the datasource and/or that doesn't rely on unique layer names. Actually the datasources are in-memory (memory-driver), and I may have up to 10 layers with up to 10.000 features over all. Maybe it would be better not to use SQL but create an index as list, and then swap features around in the layer with Layer.GetFeature()/.SetFeature()? Thanks, Redoute


title: How to de-noise a DSM

text: I have created a DEM, DSM, and nDSM from a LiDAR Point Cloud using SAGA-GIS. I extracted the appropriate returns from the point cloud, then converted to grid, and then closed gaps to fill in any nodata regions. I'm finding that my resultant rasters are still pretty noisy. For example, a region of pixels in the DSM that is obviously a tree or cluster of trees has several pixels with a very similar value to ground in stead of the surrounding higher value pixels. I've tried doing a majority filter on the DSM, but it didn't help. Any suggestions for other filters that might help create more contiguous regions of pixels?


title: Find GPS midpoint and SD of GPS coordinate along a certain heading in QGIS

text: I have a set of GPS coordinates (x, y) along with heading information. The heading for these coordinates is more or less the same and the coordinates are within a few meters of each other. How can I use QGIS (or any other tool) to find the mean and standard deviation of these GPS coordinates along the heading, i.e. looking to rotate the X, Y plane along the heading direction and then find the mean and SD. Thanks!


title: How to quantify unevenness of polygon boundaries?

text: I have two polygons: Polygon 1 and Polygon 2. Using two metrics, area and perimeter length, I want to express quantitatively that Polygon 1 has a more uneven/jagged/irregular perimeter than Polygon 2. ![enter image description here](http://i.stack.imgur.com/w5mhx.jpg) Each polygon has the same perimeter length but each covers quite different areas. To quantify the unevenness/jaggedness/irregularity of each polygon, should the calculation be: area/perimeter or perimeter/area I thought perimeter/area, but then I found this blog post which uses area/perimeter : http://www.r-bloggers.com/measuring-the-gerrymander-with- spatstat/


title: how to synchronously execute a queryTask

text: How can a queryTask be executed synchronously? I'd like a feature set to be returned instead of deferred object. Thanks,


title: Stuck installing gisgraphy: org.springframework.beans.factory.BeanDefinitionStoreException

text: I'm been trying to get gisgraphy to run for a few hours but seem to get stuck at the exception being thrown below. I thought it might have been the gisgraphy-3.0-beta2 version of gisgraphy but I now get the same with gisgraphy-2.1.1-final. I've followed the instructions pretty closely. I'm running on Windows 7 and Java version 1.8. Googling has lead nowhere -- I wonder if I have a configuration mismatch.               C:\Users\bob\apps\gisgraphy>java -Dfile.encoding=UTF-8 -Xmx1024m -Xms512m -jar start.jar      ERROR [main] ContextLoader | Context initialization failed     org.springframework.beans.factory.BeanDefinitionStoreException: Unexpected exception parsing XML document from URL [file:/C:/Users/bob/apps/gisgraphy/webapps/ROOT/WEB-INF/classes/applicationContext.xml]; nested exception is java.lang.IllegalStateException: AnnotationTransactionAttributeSource is only available on Java 1.5 and higher         at org.springframework.beans.factory.xml.XmlBeanDefinitionReader.doLoadBeanDefinitions(XmlBeanDefinitionReader.java:385)         at org.springframework.beans.factory.xml.XmlBeanDefinitionReader.loadBeanDefinitions(XmlBeanDefinitionReader.java:313)         at org.springframework.beans.factory.xml.XmlBeanDefinitionReader.loadBeanDefinitions(XmlBeanDefinitionReader.java:290)         at org.springframework.beans.factory.support.AbstractBeanDefinitionReader.loadBeanDefinitions(AbstractBeanDefinitionReader.java:142)         at org.springframework.beans.factory.support.AbstractBeanDefinitionReader.loadBeanDefinitions(AbstractBeanDefinitionReader.java:158)         at org.springframework.web.context.support.XmlWebApplicationContext.loadBeanDefinitions(XmlWebApplicationContext.java:124)         at org.springframework.web.context.support.XmlWebApplicationContext.loadBeanDefinitions(XmlWebApplicationContext.java:92)         at org.springframework.context.support.AbstractRefreshableApplicationContext.refreshBeanFactory(AbstractRefreshableApplicationContext.java:97)         at org.springframework.context.support.AbstractApplicationContext.obtainFreshBeanFactory(AbstractApplicationContext.java:411)         at org.springframework.context.support.AbstractApplicationContext.refresh(AbstractApplicationContext.java:338)         at org.springframework.web.context.ContextLoader.createWebApplicationContext(ContextLoader.java:251)         at org.springframework.web.context.ContextLoader.initWebApplicationContext(ContextLoader.java:190)         at org.springframework.web.context.ContextLoaderListener.contextInitialized(ContextLoaderListener.java:45)         at org.mortbay.jetty.handler.ContextHandler.startContext(ContextHandler.java:530)         at org.mortbay.jetty.servlet.Context.startContext(Context.java:135)         at org.mortbay.jetty.webapp.WebAppContext.startContext(WebAppContext.java:1218)         at org.mortbay.jetty.handler.ContextHandler.doStart(ContextHandler.java:500)         at org.mortbay.jetty.webapp.WebAppContext.doStart(WebAppContext.java:448)         at org.mortbay.component.AbstractLifeCycle.start(AbstractLifeCycle.java:40)         at org.mortbay.jetty.handler.HandlerCollection.doStart(HandlerCollection.java:147)         at org.mortbay.jetty.handler.ContextHandlerCollection.doStart(ContextHandlerCollection.java:161)         at org.mortbay.component.AbstractLifeCycle.start(AbstractLifeCycle.java:40)         at org.mortbay.jetty.handler.HandlerCollection.doStart(HandlerCollection.java:147)         at org.mortbay.component.AbstractLifeCycle.start(AbstractLifeCycle.java:40)         at org.mortbay.jetty.handler.HandlerWrapper.doStart(HandlerWrapper.java:117)         at org.mortbay.jetty.Server.doStart(Server.java:210)         at org.mortbay.component.AbstractLifeCycle.start(AbstractLifeCycle.java:40)         at org.mortbay.xml.XmlConfiguration.main(XmlConfiguration.java:929)         at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)         at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)         at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)         at java.lang.reflect.Method.invoke(Unknown Source)         at org.mortbay.start.Main.invokeMain(Main.java:183)         at org.mortbay.start.Main.start(Main.java:497)         at org.mortbay.start.Main.main(Main.java:115)     Caused by: java.lang.IllegalStateException: AnnotationTransactionAttributeSource is only available on Java 1.5 and higher         at org.springframework.transaction.config.TxNamespaceUtils.getAnnotationTransactionAttributeSourceClass(TxNamespaceUtils.java:43)         at org.springframework.transaction.config.AnnotationDrivenBeanDefinitionParser$AopAutoProxyConfigurer.configureAutoProxyCreator(AnnotationDrivenBeanDefinitionParser.java:117)         at org.springframework.transaction.config.AnnotationDrivenBeanDefinitionParser.parse(AnnotationDrivenBeanDefinitionParser.java:79)         at org.springframework.beans.factory.xml.NamespaceHandlerSupport.parse(NamespaceHandlerSupport.java:69)         at org.springframework.beans.factory.xml.BeanDefinitionParserDelegate.parseCustomElement(BeanDefinitionParserDelegate.java:1246)         at org.springframework.beans.factory.xml.BeanDefinitionParserDelegate.parseCustomElement(BeanDefinitionParserDelegate.java:1236)         at org.springframework.beans.factory.xml.DefaultBeanDefinitionDocumentReader.parseBeanDefinitions(DefaultBeanDefinitionDocumentReader.java:133)         at org.springframework.beans.factory.xml.DefaultBeanDefinitionDocumentReader.registerBeanDefinitions(DefaultBeanDefinitionDocumentReader.java:90)         at org.springframework.beans.factory.xml.XmlBeanDefinitionReader.registerBeanDefinitions(XmlBeanDefinitionReader.java:468)         at org.springframework.beans.factory.xml.XmlBeanDefinitionReader.doLoadBeanDefinitions(XmlBeanDefinitionReader.java:363)         ... 34 more


title: Python error with open layers plugin QGIS 2

text: I've been using QGIS for a while but since using QGIS 2 the OpenLayers plugin doesn't work for me on my home pc. (However, I have installed the same version at my workplace with success.) I have looked through other posts, and tried setting the crs and allowing on the fly but it doesn't work for me.   Please can anyone offer any other guidance? I have also tried uninstall/ re- install. Error code below: > An error has occured while executing Python code: > > Traceback (most recent call last): File > "C:\Users\Lee/.qgis2/python/plugins\openlayers_plugin\openlayers_plugin.py", > line 48, in addLayer self. **plugin.addLayer(self) File > "C:\Users\Lee/.qgis2/python/plugins\openlayers_plugin\openlayers_plugin.py", > line 194, in addLayer layer = OpenlayersLayer(self.iface, > self.__coordRSGoogle, self.olLayerTypeRegistry) File > "C:\Users\Lee/.qgis2/python/plugins\openlayers_plugin\openlayers_layer.py", > line 64, in __init** self.page = OLWebPage() File > "C:\Users\Lee/.qgis2/python/plugins\openlayers_plugin\openlayers_layer.py", > line 38, in **init** proxy = getProxy() File > "C:\Users\Lee/.qgis2/python/plugins\openlayers_plugin\tools_network.py", > line 42, in getProxy proxy.setPort(settings.value("/proxyPort", 10, > type=int)) TypeError: unable to convert a QVariant of type 10 to a QMetaType > of type 2 > > Python version: 2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit > (AMD64)] > > QGIS version: 2.0.1-Dufour Dufour, d94c044 > > Python path: ['C:/PROGRA~1/QGISDU~1/apps/qgis/./python/plugins\processing', > 'C:\Users\Lee\.qgis2\python\plugins\GeoCoding', > 'C:/PROGRA~1/QGISDU~1/apps/qgis/./python', 'C:\Users\Lee/.qgis2/python', > 'C:\Users\Lee/.qgis2/python/plugins', > 'C:/PROGRA~1/QGISDU~1/apps/qgis/./python/plugins', > 'C:\PROGRA~1\QGISDU~1\bin\python27.zip', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\DLLs', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\plat-win', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\lib-tk', 'C:\PROGRA~1\QGISDU~1\bin', > 'C:\PROGRA~1\QGISDU~1\apps\Python27', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site-packages', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site-packages\PIL', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site-packages\win32', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site-packages\win32\lib', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site-packages\Pythonwin', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site- > packages\Shapely-1.2.18-py2.7-win-amd64.egg', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site-packages\six-1.3.0-py2.7.egg', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site-packages\wx-2.8-msw-unicode', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site-packages\xlrd-0.9.2-py2.7.egg', > 'C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site-packages\xlwt-0.7.5-py2.7.egg', > 'C:\Users\Lee\.qgis2\python\plugins\mmqgis/forms', > 'C:\PROGRA~1\QGISDU~1\apps\qgis\python\plugins\fTools\tools']


title: How to run a shortest path analysis for DEM-derived streams and pour points?

text: I have a .shp for streams and another for pour points. Those features were extracted from a DEM by using the Hydrology package from the Spatial Analyst Tools of ArcGIS 10. I want to know the distance between every pour point along the flow path network. What procedure should I use?


title: How to correctly use Datepart( ) in ArcGIS?

text: I have two date fields dd.mm.yyyy with a start and an end date in my attribute table. Now I want only to extract the years from these fields to new fields. So I created another field, and used `year_start_date = DatePart ( "YYYY", [Start_date] )` in the field calculator. But the result is a complete date (dd.mm.yyyy) more than 100 years ago. Just to test, I created a new field and added with the function `Now()` the actual date (13.05.2014) and tried again to use `year_start_date = DatePart ( "YYYY", [Start_date] )` and I got "06.07.1905" as a result. What am I doing wrong? Is there another way to do this? I use ArcGIS 10.2.


title: Mouse right-click is not working anymore when using ArcGIS10

text: The mouse right-click option has stopped working with ArcGIS10. It happened suddenly - the mouse is ok (not broken) as it works with other programs. Do you know what happened and how to fix it? The right-click option is useful for editing as you can just right-click select from the drop down list, absolute coordinates, direction, delete sketch etc. It is also not working, when I'm in layout view as I can't change grid properties.


title: Are there any OpenSource ArcObjects Projects?

text: One of the major benefits of developing GIS applications based on opensource libraries, is that many projects spring up on top of them which are also opensource. Well designed opensource projects are great for seeing:   * which coding patterns should be used   * what coding styles to follow   * how to structure larger projects   * what unit tests to include   * autodocumentation tools and procedures   * sourcecontrol workflows The ArcObjects Code Gallery seems the best place to look, but most examples I've looked at are very small pieces of functionality or single tools, rather than full applications. The SharpGIS project is a nice .NET example to follow, but are there any ArcObjects specific projects? * * * **Updated with Examples:** http://eagl.codeplex.com/ > Easy ArcGIS Library is a set of C# .net classes that wrap the common > functionality of ArcObjects, that help ArcGIS developers do a lot of common > functions in less time and direct way. And from the same developer: http://arcgisslcontrols.codeplex.com/ > ArcGISSLControls is a set of ArcGIS Silverlight controls that is built on > top of the ArcGIS Silverlight SDK, it is developed in C# programming > language.


title: polygon map overlay algorithm

text: There are many algorithms for binary map overlay operation in vector data format which take two layers of map and produce resultant layer i.e overlaid layer as output. I am wondering whether there are any algorithms which take more than two layers say 3 layers simultaneously and produce the overlay result?


title: How do you add a WKT specified coordinate system to Arcgis?

text: I've got a WKT specified coordinate system, how do I add this to ArcGis?


title: Errors recognizing ArcHydro 2.0 and ArcGIS 10 when calculating ServiceSheds with InVEST HelperTools

text: after validating arguments, creating output folders, setting script variables, and checking DEM properties, the Create ServiceSheds tool from InVEST Helper Tools is telling me that ArcGIS10 and ArcHydro 2.0 are required. I have downloaded and installed both of these, and so I do not understand where the issue is coming from. It also says further down that the script failed and raises an exception in the last line of code. Any tips from anyone with experience with either ArcHydro 2.0 or Helper Tools or problems like these would be greatly appreciated!


title: i want to create a dam for a particular location

text: i have created a DEM for a particular location, but i am stuck on how to find the best posible location to place a dam (and how to even make a dam) of 3om in height and 300m in length. please i need every possible idea on how to come about this. thanks


title: Where can I find GIS Data for Environmental Constraints for Hawaii?

text: I am looking for all the types of environmental constraints in shapefile format for Hawaii. Is there any environmental agency website or other source where I could find it? Any kind of help would be much appreciated.


title: Longitudinal Issues with NetCDF data

text: I am seeking advice on working with world data sets where the GCS has a positive range (0 degrees to 360 degrees). The data I am working with are NetCDF oceanographic data and have positive coordinate values as mentioned. When displayed in a normal GCS WGS84 in ESRI's ArcGIS 10, it is offset from other data that exist in the range of -180 to 180. If I reproject it, gaps occur around the prime meridian (mostly due to its origins, sometimes as far off as 25 degrees W). My thought is to create a custom projection/coordinate system that has the positive values and hope that other world data sets reproject on the fly just fine. Does anyone see any issues with this or have other solutions? Hope this makes sense. I appreciate any help in advance. ** Update ** This is a screenshot of a countries data set with graticule on top to show the normal CRS of GCS. Just importing the NetCDF data as-is brings it in but the values don't start until 20.5 degrees east of the Prime Meridian. ![Normal GCS](http://grafa.co/rnd/img/GCS_Normal.png) If I choose to use the NetCDF's CRS (which is actually the same) it will reproject the world data on the fly if it is defined. Note the graticule is undefined so it does not reproject. ![Sort of reprojected GCS](http://grafa.co/rnd/img/GCS_Other.png) But then if I try to reproject everything to a normal GCS with negative values, it's like wrapping a flat map of the data around the world and then it disappears when it hits the Prime Meridian. ![Sort of reprojected GCS](http://grafa.co/rnd/img/reprojected.png) Now, I know there are no values from 0-20.5 as is stated in the metadata. But why can't the values from rest of the data display in the gap? I even tried a Shift in the Raster Tools to no avail.


title: Select and highlight features

text: I have 2 ESRI feature layers on my map; points and polylines. I can select them just fine either by rectangle or polygon. What I am missing is how to highlight (larger size and different color) the selected points and polylines for display. Do I need to create new layers with just the selected features and use a renderer or place a graphic container on the map with those features highlighted? The ultimate goal is to have the selected and highlighted features available for deletion. Using ArcEngine and ArcSDE. C# .NET. Thank you.


title: Parse Data Through ArcMap Toolbox

text: In a toolbox of mine I want to add a function to add several shapefiles by assigning the datatype to shapefile and enabling the multi-value function. I have noticed however that ArcMap parses this data as a long string deliminated by a ';' of which I have made into a list by the following.               def main(inputs,...)         shp_files = inputs.split(';')      However these are still strings in a list and therefore when I run a geoprocessing task the error I receive is that the shapefile has no table 'C'. C being my directory C:. What I need is either a more elegant approach to my current step or a way to add r' ' to surrond each of my list strings. I assume there must be a way to append each string in a list and add the r '' to indicate that is a pathname.


How can I overwrite the url of OpenLayers.Protocol.HTTP? I tried > searchformPanel.protocol.url and it works (checked with console.log), but at the end the original url is sent (in the code below: url: '/fs/') (see attached picture). ![output from firebug](http://i.stack.imgur.com/mY1pq.png) This is the code:               var searchformPanel = new GeoExt.form.FormPanel({       border: false,       width: 250,       protocol: new OpenLayers.Protocol.HTTP({         url: '/fs/',         format: new OpenLayers.Format.GeoJSON()       }),       items:[{         xtype: 'combo',         id: 'idcombo',         store: new Ext.data.SimpleStore({           fields:['fsclass','ollayer'],           data:[["Boreholes","Boreholes"],["All_layers","All layers"]]         }),         displayField: 'ollayer',         valueField: 'fsclass',         fieldLabel: 'Layer',         emptyText: 'select a layer',         submitValue: false,         selectOnFocus: true,         mode: 'local',         typeAhead: true,         editable: false,         forceSelection: true,         triggerAction: 'all'       },{         xtype: 'textfield',         id: 'idtextfield',         fieldLabel: 'Find features',         emptyText: 'enter word',         name: 'comments__like',         allowBlank: false       }],       listeners:{         actioncomplete: function(form, action){           searchShowTip(action.response.features);         }       },       buttons:[{         text: 'search',         listeners:{           click: function(){                  var comboLayer = Ext.getCmp('idcombo').getRawValue();             var keyword = Ext.getCmp('idtextfield').getRawValue();                  var newUrl = '/fs/' + comboLayer + '?format=GeoJSON&comments__ilike=' + keyword + '&queryable=comments';             console.log('1:' + newUrl);                  //this gets '/fs/' from the searchformPanel             console.log('2:' + searchformPanel.protocol.url);                  searchformPanel.protocol.url = newUrl;             console.log('3:' + searchformPanel.protocol.url);                  searchformPanel.search();           }         }       }]      }); Please any support about this is very welcome, thank you!


title: What are the Best-practices for deploying services, data, & code to an ArcGIS production server?

text: We have a production GIS viewer that runs on a production server. Testing new development/data is done on a test server and we develop on a different server also. At the moment, a lot of the deployment is manual work, with risk for errors. I would like to automate whatever is possible to reduce that risk. I was wondering what kind of best practices are used to maintain and deploy new data (Almost all data for our viewer is kept in Filegeodatabase format). My main questions are these:   * What kind of scripts/tools do you use to deploy new maps,services,data to your production server,      * do you deploy them manually? In ArcCatalog or the web manager?     * do you copy and paste the cfg files, do you change the settings of those files manually or with scripts?   * Do you keep stuff in version control? What do you keep there?   * How do you manage different settings between servers (min and max number of SOC's per service for example)   * What kind of best practices in general do you recommend. Are there good docs to read about it? Thanks!


title: Geoexplorer displays the layers I import, but not the table behind them

text: If I import layers into GeoSrver, they display in GeoExplorer, and I can use the identify tool to find more info about areas of the layer I select. However the table from which this data comes from does not show at the bottom of the screen, as it does for existing example layers. I assume there is a setting somewhere to enable this?


title: GeoServer: Link meta data to raster image?

text: I want to use the OpenGeo Suite to publish different sets of raster imagery with associated data. Doing this with vectorial data works fine:   1. Setup the data base in postGIS: import the vectorial data, create tables containing data describing the objects.   2. Setup a data store in GeoServer: make a connection to the data base.   3. Setup a layer in GeoServer: combine vectorial data and attributes with a SQL view. After that I understood that I can easily access attributes via `GetFeature` requests to a WFS. However, the same with raster data instead of vectorial data seems not to work like that. One can import raster data from postGIS (via the ImageMosaicJDBC plugin) but I see no means to link it in the same way as above. I'm new to web mapping, so: Is it my ignorance of the matter or the lack of a feature in GeoServer?


title: PyProj/Proj4 convert geographic to projected coordinates for Geostationary Data

text: I am trying to convert geographic coordinates (in degrees east/north lon lat) to projected coordinates in meters for geostationary full disc data (GVAR Goes East data from CLASS). To do this, I tried using PyProj/ Proj4 and am running into errors due to the off earth pixels in the corners of the image. This is the proj4 string given:               projection_coords = Proj("+proj=geos +h=35774290 +a= 6378137 +b= 6378137     +lon_0=-75 +units=meters +no_defs")          ll_x, ll_y = projection_coords_geos( LL_y_deg, LL_x_deg, inverse = False, errcheck=True)     'x=%9.3f y=%11.3f' % (ll_x,ll_y)          ur_x,ur_y = projection_coords_geos( UR_y_deg, UR_x_deg, inverse = False, errcheck=True)     'x=%9.3f y=%11.3f' % (ur_x,ur_y)      And the input 'LL_y_deg', 'LL_X_deg' etc are the values of the corner point extents lower left y pixel, lower left x pixel, upper right y and upper right x. For full disc geos data I am receiving incorrect (1e30) values for the extents and was wondering if anyone has ever run into the same error or have a suggestion to correct this. Here is an image to describe the problem further: ![Bounding Box Issue](http://i.stack.imgur.com/ElVRg.png) We did try drawing a limiting box inside the earth extents to get all valid data points within, and this works but cuts out data in the edges of the earth (the data in between the blue dotted box and the earth). The geographic coordinate values in the red shaded off earth areas are very large (2.14 E9) and these are the values that throw the error. Also, for remapping GOES data, should the 'sweep axis' always be included in the Proj string and set to 'y'? Any suggestions or comments would be greatly appreciated!


title: How to remove duplicate points using ArcGIS ModelBuilder?

text: I want to remove duplicate points and I have **this model**. I want to deselect all feature that have the same value of the OID. So I want to add value of the OID in the query of Select by Attribute. How do I incorporate this value in a SQL query ??? I try : "%Feature OID%" = "OID" (does not work) and a lot of other tests and nothing work Someone have an idea?


title: Making a KML layer "always on"

text: Is there a KML tag or specification by which you can force a given layer (specifically, a screen overlay) to be "always on", in that there is no way for the user viewing the file to make that layer invisible, so long as they are viewing anything in that file? I'm looking in particular for a solution that will work when Google Earth is used to view the KML file. That is, if my KML file has multiple features and layers, one of which is called "Overlay", say, then if the user loads my KML file and turns on any of the layers in it, then Overlay will automatically turn on. If the user tries to turn Overlay off while something else in the file is still visible then either (a) Overlay cannot be turned off, or (b) everything else in the file is made invisible along with Overlay. Of course I have no problem with the user removing my KML file from the GE session, or disabling it entirely, but if any of it is shown, then Overlay should be as well.


title: What is the meaning of "MULTIPOLYGON EMPTY" Message?

text: I am trying to run ST_Intersection function. SELECT ST_AsText(ST_Intersection(ST_Point(42.729123082000058,-96.619302355999992), ST_Point(42.72899715300008,-96.619410439999967))) FROM clu_08aug; I am getting Resultset with the "MULTIPOLYGON EMPTY". Can someone help me what does it means? Dharmesh


title: r.sun delievering unrealistic Values

text: I need your help. I'm trying to calculate the global irridiation for a specific area. I'm interested in the results for rooftops in special. I have a DEM, calculated slope, aspect and horizons and started `r.sun` for every day with a monthly changing Linke Factor. All calculated maps (slope, aspect, horizons) seem to have believable values. My Problem is: On the yearly sum (addition of all 365 days) I get values for total global irridiation of about 1900 KW/m². Which is 1,9times more then expected. Especially when the solar constant is only 1368 kW/m².... Could you imagine what went wrong? I'm lost.


title: Preserve polygon boundaries when vectorizing a raster according to a polygon vector

text: I have a raster which I clip by some polygons and then convert to polygon. I would like to obtain the same number of polygons in both my vector files. The problem is that some areas in the vector clipper are separated by a line (i.e. adjacent polygons) are now together (as I would expect). When I perform the vectorization of my cropped raster, I get less polygons than what I want (for the reason I just mentioned). Is there a way to preserve polygon boundaries when performing raster to vector transformation according to a given polygon vector (that is, my clipper file)? I use ArcGIS, QGIS, GRASSgis, python, MATLAB/octave. Any help is appreciated. ![enter image description here](http://i.stack.imgur.com/sWAH3.jpg)


title: High-DPI printing with ArcGIS API for JavaScript and esri.dijit.Print?

text: I've got cached tiles for all map layers at various zoom levels, and would like to print at DPI higher than the default (around 96) by using tiles from a deeper zoom level. It doesn't matter if map features become very small or unreadable. I've tried passing a higher DPI setting:               template=new PrintTemplate();     template['layout']='MAP_ONLY';     template['exportOptions']={dpi:300,width:3200,height:2400};     template['preserveScale']=false;     template['format']='PDF';          printWidget=new Print({         async:true,         map:map,         url:url,         templates:[template]     },dom.byId('printButton'));          printWidget.startup();      However pixels in the resulting PDF are still very large. The server seems to choose the zoom level of the tiles based on the map extent only, but always according to about 96 DPI. The preserveScale option does nothing (I checked that it doesn't get passed to the server at all) and the DPI setting is sent to the server but doesn't affect the resulting image. What can I do to increase the output DPI by printing tiles from a deeper zoom level than is shown on screen? Creating an A3 or A2 sized printing template and instructing users to configure their printing settings to scale it to A4 fixes the problem, but it would be preferable to produce a correct PDF.


title: How to load external geojson file into leaflet map

text: I would like to load a geoJSON (polygon) file into my leaflet map. I have seen examples where geoJSON is embedded into the javascript code but I can't find any examples showing how it is done with an external file.               <meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=no">          <link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.5/leaflet.css" />


title: Rounding off coordinates in an SpatialPolygonDataFrame in R

text: EDIT AFTER COMMENTS After a comment of Whuber (Thank you!) I realized it will probably be most easy to incorporate a round(...) statement within my code to be able to work with less decimals. Underneath I have added the functions I am using. My first function (allcoordinates) extracts all coordinates of a SPDF. My second function (transformloop) transforms in a loop. After every transformation I want to round off my coordinates. Simply using round(file@polygons, 2) or signif(file@polygons, 2) does not work. The polygons of a SPDF are not seen as numeric values apparently. Anybody an idea how to round off coordinates within a SPDF?               allcoordinates <- function(x) {       ret <- NULL       polys <- slot(x,"polygons")       for(i in 1:length(polys)) {         pp = polys[[i]]@Polygons         for (j in 1:length(pp))           ret = rbind(ret, coordinates(pp[[j]]))       }       ret      }          transformloop <- function(inputfile, count){       intermediate <- spTransform(inputfile, CRS("..."))       inputfile <- spTransform(intermediate, CRS("..."))       if (count==0){         return(inputfile)       } else           newcount <- count-1       transformloop(inputfile, newcount)           }      I have added the first function cause I am thinking that possibly the best option is to extract the coordinates, round them off, and then replace the existing coordinates with the coordinates I rounded off. ORIGINAL QUESTION (outdated): In RStudio I am transforming data from coordinate system to coordinate system and back. Then I compare the coordinates of the original file with the endfile. This shows minor differences 7 digits behind the comma. However, I want to know what happens with the numbers if a system would work with less numbers behind the comma. Within R I managed to find the code options(digits=...). This option only visualizes the numbers in a certain amount of digits. The calculation appears to do exactly the same. This is absolutely correct and usually a very good characteristic. But it is not what I want right now. I want it to be less accurate. Is there a way to manipulate R in order for it to work with a different amount of numbers behind the comma constantly?


title: OpenLayers SVG sprite externalGraphic only works in Firefox

text: I am trying to set the externalGraphic property of a vector layer's style to an SVG sprite that has various viewboxes to position the icons dynamically by creating different views within my SVG file:               <view id="bikeview" viewBox="0 0 52 52"></view>     <view id="boatview" viewBox="52 0 52 52"></view>      ...and then accessing them like so:               var defPoiStyle = {         externalGraphic: 'images/activities/activitySprite.svg#boatview',         pointRadius: 30,         cursor: "pointer"     }      The problem is that this only works in Firefox, although with Chrome and Safari I can access each view by entering the URL directly into the browser like this, just not in the map. I have also noticed that in Chrome's network log (and likely Safari's) that the `#boatview` part of my image path gets stripped out, but not so in Firefox. I assume it is some kind of XML request issue or something, but I feel like I got my declarations right:               <?xml version="1.0" encoding="UTF-8" standalone="no"?>     <!DOCTYPE svg PUBLIC "-//W3C//DTD SVG 1.1//EN"        "http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd">     <svg        xmlns:svg="http://www.w3.org/2000/svg"        xmlns="http://www.w3.org/2000/svg"        xmlns:xlink="http://www.w3.org/1999/xlink"        x="0px"        y="0px"        width="52px"        height="52px"        id="svg3753"        version="1.1">      Sorry, kinda complicated for jsFiddle, but here is my map. Give it a go in Chrome and then in Firefox, and you should see what I mean (the boat icon is what is supposed to display). Any help is appreciated!


title: ArcGIS Server 10.1 - GPU Processing Support

text: Does ArcGIS Server 10.1 or any of its extensions (e.g. 3D Analyst, Spatial Analyst) take advantage of GPU processing? Dell have a new server, the R820, that has an option to include a GPU. Hence the interest.


title: Is it possible to re-point operational layers ArcGIS Mobile projects v.3.0

text: We have an ArcGIS mobile project with two operational layers for data collection. The layers source from a GIS server on the cloud, running on top of an SDE. The database has a feature dataset, and the editable feature classes within. The editable feature classes have a lot of fields, and the mobile project has all the layers' and their fields' settings. We want to deploy projects for a number of different regions, thus've we built a separate feature dataset and its child feature classes for each region, and published the GIS services for the mobile devices to interact with. We want to replicate the mobile project and have one project for each region. Question: is it possible to make a copy of a mobile project and re-point the operational layers to a different service and still keep the settings (fields are the same as all feature classes have the same schema), to avoid having to configure the layers by manually doing all the ticks again? Or, is there a way to have the settings in a text file to read from? We're using Mobile Project Center 3.0. Thanks


title: How do I construct a file path within Python for a Feature Class within a Feature Dataset?

text: Based on my understanding of the ArcGIS path structure, a feature class inside a feature dataset inside a geodatabase can be accessed by the following:               C:\Menu\Lunch.gdb\Sandwich\Cheese      Where Sandwich is the feature dataset and Cheese is the feature class. In a python script, I want to construct this full path, so that I can check it using arcpy.Exists(x). I have the following variables (pseudocode):               WSpace = C:\Menu\Lunch.gdb\     FDSName = Sandwich     FCName = Cheese      Concatenating WSpace and FDSName gives the correct path to the feature dataset. What I need help with is tacking FCName onto the end. I would think the following code would work, but it renders a space instead of the backslash character:               fullPath = WSpace + FDSName + r"\" + FCName      Instead of r"\", I have also tried "\", and "/", but they all render a space. This seems like a simple syntax problem. Any pointers would be appreciated.


title: Why are results from Python tool not in Results window?

text: I am relatively new to ArcPy 10. and I want to return the list of feature class from a dataset. I can see them through the message, but cannot see them from the result window. Actually there are only inputs, environments and messages I can view from the results window. Am I missing something in the script? Thanks.                   import arcpy         import json         from arcpy import env         env.workspace=r"C:\sde\pan18.sde"         featureDatasets = arcpy.ListDatasets("SDE.POINT_OF_INTEREST")         fd=featureDatasets[0];         fcList=arcpy.ListFeatureClasses("*", "ALL", fd)         pf=[]         for fc in fcList:            pf.append(fc)         arcpy.AddMessage("--------Encode the list to a json string-\n")         layerlistEncode=json.dumps({'resource':pf})         print layerlistEncode         arcpy.AddMessage("--------Encode the json string to a list\n")         layerlistDecode=json.loads(layerlistEncode)         print layerlistDecode         print layerlistDecode["resource"]         arcpy.SetParameterAsText(0, layerlistEncode)              arcpy.AddMessage(layerlistEncode)         arcpy.AddMessage("finished1")          Edit / Delete Edit Post   Reply With Quote Reply With Quote   Multi-Quote This Message  Top Bottom


title: Trying to understand / troubleshoot intersection behavior of postgis

text: I'm trying to calculate the percent of intersection between two polygons using postgis. I was able to do so successfully for 90% of my polygons, but about 1500 have very unexpected values. I extracted one of these cases to examine the behavior. The orange polygon resides w/in the purple polygons so the proportion should be near 100%. ![the two polgyons](http://i.stack.imgur.com/XTPC5.png) However, after running:               UPDATE orange      SET pct_undevelopable = ST_Area(ST_Intersection(orange.geom, purple.geom)) / ST_Area(orange.geom)     FROM purple WHERE ST_Intersects(orange.geom, purple.geom);      I get the value for percent of area of intersection as: 0.000828452095134945 So I inspected the intersection part of the query and found that I got 3 rows as geometry outputs. I exported all three as individual shape files: one simply doesn't work...I tried to open it in QGIS and ArcGIS and I also tried in CSV format for QGIS, no luck. And two polygons. Polygon 1 is what I imagine PostGIS is calculating percent of intersection from: ![enter image description here](http://i.stack.imgur.com/Wn99a.png) Polygon 2 is what I imagine I would like PostGIS to calculate it from: ![enter image description here](http://i.stack.imgur.com/lzp2s.png) I would love it if anyone can 1) explain what is going on (especially with the weird non functioning intersection geometry) and 2) how to fix it. Thank you!


title: Adding time stamps to animation slides with OSS tools

text: I would like to create an animation with time lapse values printed on the PNG file. For this purpose I have created in QGIS a vector layer that stores timestamp and temperature value pairs for various points in a harbour. I interpolate the vector points layer in GRASS using v.surf.icw that creates a raster layer and used those layers with the GUI animator to create the movie. Everything looks good except that I would like to have the time lapse printed on each slide. Therefore, I have added in the vector layer a time stamp for each data point. the header looks something like: > point1, lat, long, time1, temp11, time2, temp12, ...... > > point2, lat, long, time1, temp12, time2, temp22, ...... I would like now to create a script that reads the time values and write then on the PNG image. One way I found is to use d.text in grass as documented here. However, this example shows how to create new slides not how to write on existing ones. In addition I would need to extract the time stamp information from the vector layer I assume with d.vect or something similar. I would appreciate any input on these 2 matters:   1. Write with d.text on an exiting PNG image   2. Extract an attribute value from a vector layer


title: Multi layer, multi control

text: I am new to GIS technolgies. I'm working on the final project of my career. I need to have a background map of somewhere on this map I need to create 1 .. N layers in which a particular user can add, modify or delete information. I was able to create layers, and I was able to add or delete information about these. But I can not change information on the layers. By modified I mean move a point or polygon My code is like this.               var objectListCustomLayer = jQuery.parseJSON(listCustomLayer);     for (var i=0; i < objectListCustomLayer.length; i++){         var item = objectListCustomLayer[i];         loadCustomLayers(item.nameLayerkey,item.layerkml,item.layerForm);     }      the objectListCustomLayer take info from server and have information over the layers. and here call the loadCustomLayers... function loadCustomLayers(name,customKml,customForm){ var renderer = OpenLayers.Util.getParameters(window.location.href).renderer; renderer = (renderer) ? [renderer] : OpenLayers.Layer.Vector.prototype.renderers;               customLayer = new OpenLayers.Layer.Vector(name,renderer);          var info = new infoCustomLayer(customLayer.id,customForm,customLayer);     vectorCustomLayer.push(info);          map.addLayer(customLayer);            deserialize(customKml,customLayer);          var drag = new OpenLayers.Control.DragFeature(customLayer,{ onStart: startDrag});     var select = new OpenLayers.Control.SelectFeature(customLayer, {box: true, multiple: true,onSelect: addSelected,onUnselect: clearSelected});          var drawOptionsPoint = {         callbacks : {"done": doneHandlerPoint,                     "point": pointHandlerPoint                      },         };          var drawOptionsPath = {         callbacks : {"done": doneHandlerPath},         };          var drawOptionsPolygon = {         callbacks : {"done": doneHandlerPolygon},         };          controls = {         point:  new OpenLayers.Control.DrawFeature(customLayer,OpenLayers.Handler.Point,drawOptionsPoint),         line:   new OpenLayers.Control.DrawFeature(customLayer,OpenLayers.Handler.Path,drawOptionsPath),         polygon:new OpenLayers.Control.DrawFeature(customLayer,OpenLayers.Handler.Polygon,drawOptionsPolygon),         drag:   new OpenLayers.Control.DragFeature(customLayer,{ onStart: startDrag}),         select: select         };               for(var key in controls) {                 map.addControl(controls[key]);             }      } The function loadCustomLayers is call N times. My problem is when I select OpenLayers.Control.DragFeature this works over one layer, but does not work over the layer selected by me. controls[key] is a vector that stores the name of the controls, is used by a call for a javascript event that allows me to indicate what to do (draw, select or drag) on a given layer.


title: How to display vectors in arcgis json format in openlayers?

text: For this project I am working on, I would like to display an ArcGIS Online FeatureService (JSON) in a webmap that I have built using OpenLayers. I am able to read the JSON from the service, but I can't seem to 'convert' into a vector layer. My code is the following:               var data = ($.ajax({         url: '...',         success: function (data) {         console.log(data);         }              }));             var geojson_format = new OpenLayers.Format.GeoJSON();        var vector_layer = new OpenLayers.Layer.Vector();        map.addLayer(vector_layer);        vector_layer.addFeatures(geojson_format.read(data));      This unfortunately produces the following errors:               TypeError: g is null (FF)          TypeError: 'null' is not an object (evaluating 'g.layer=this') (Safari)          Uncaught TypeError: Cannot set property 'layer' of null (Chrome)       Does anyone know how to fix this? Thnx in advance!


title: How to get from a 7-folder-deep directory structure to a 3-folder structure for TMS?

text: so I got a Tiled Map Service up and running using TileCache to store all my tiles and OpenLayers to actually show them in a browser using an Apache server. But the demo map tiles I found were stored in a 7 folder deep directory, e.g.               13/000/002/123/000/004/936.jpg       would be zoom level 13, column 2123, row 4936. But I've been given tiles formatted just 3 folders deep, e.g.               13/2123/4936.jpg.       I'm guessing the first, 7-segment is the TMS standard since it worked correctly and everything displayed correctly, but I was wondering if there was any easy fix to display the 3-segment directory correctly, or am I just going to have to change it to fit the 7-folder format? And is there any easy way to do that? Because that'd be a lot of files to have to change. Thanks.


title: kml divide map into regions

text: in google map,               how to show Douglas county MAP in regions manner. (east, west, north, south).       These regions should cover with border and show the picker when click the balloon. Here is my KML file, which shows the Douglas county, now i want to split under region. (East,West,North,South, NE,NW,SE,SW) my kml file this kml output display here.![map](http://i.stack.imgur.com/xSoNv.jpg) i posted this post in stackoverflow, but it was closed due to ambiguous. can you please some share your ideas.


title: Changing font size, when using the Easy Custom Labeling plugin?

text: I’m using the Easy Custom Labeling plugin. Useful plugin, but when I want to change the font size with the Layer Labeling Settings. Nothing changes, I can only change font size manually for each label? Also changing it in the attribute table is not possible, as only one integer value is allowed. Removing the column and adding a new one for two integers, does not work either? I am using QGIS 1.8.0 Lisboa


title: How to create hibernate layers at runtime?

text: We need to dynamically create hibernate layers from shapefiles. As far as we understand, the workflow is following:   1. Import shapefile to database(creating new table for each shapefile?) Is it good approach?   2. Dynamically create hibernate entity from the table   3. Make the Entity visible for the Hibernate mechanism   4. Create Hibernate layer in runtime Are we moving in the right way? How would you handle our problem? Is it good approach to creating new table for each shapefile? If it so, how can resolve steps 2-4.


title: Geoserver api: how to determine a layers respective workspace?

text: I am using the geoserver api to find information about layers that are available. I am using:               http://localhost/geoserver/rest/layers.json      to get a dynamic list of all layers. If I then select a layer with:               http://localhost/geoserver/rest/layers/<myLayer>.json      I get information about the layer. How can I determine the layers respective workspace?


title: Divide shapefile

text: I'm new working with QGIS 1.8.0 and wanted to know if its possible to subdivide a map (shapefile, which is round-like) into small squares but following the outer round contour of the map? I tried with the vector grid tool but the subdivision does not limit to the map and I get only a big square. Any ideas? ![Sample Image](http://i.stack.imgur.com/3ht9D.png)


title: Removing sliver polygons, ArcGIS 10

text: I unioned 2 layers into 1 layer, and it created sliver polygons. I used the elimination tool to merge them with neighboring polygons. As as result, it merges all sliver polygons to only one polygon (selected in the picture). What I need, is that, every sliver polygon to be added to its neighboring polygon. How that could be achieved? ![enter image description here](http://i.stack.imgur.com/U98oC.png)


title: Calculating area with python qgis

text: I am a little stumped with determining the area of polygons in the python console. I want the output to be in meters^2, so using: `qgis.utils.iface.mapCanvas().setMapUnits(0)` `qgis.utils.iface.mapCanvas().refresh()` `cLayer = qgis.utils.iface.mapCanvas().currentLayer()` `for f in cLayer.getFeatures(): print f["ID"], f.geometry().area()` I get the the output of: `1 0.000114895410903 2 0.00014639997019...` - which suggests to me that the map units are still in in decimal degrees and not meters squared like I expected. with some background reading, I then specified a localised CRS: `wkt = 'PROJCS["SACRS_NO_27",GEOGCS["WGS 84",DATUM["unknown",SPHEROID["WGS84",6378137,298.257223563],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",27],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]'` and applied it with: `coo = QgsCoordinateReferenceSystem(wkt)` and `cLayer.setCrs(coo)` I can check the mapUnits of cLayer with `print cLayer.crs().mapUnits()` and this gives me '0' for metres. Yet `for f in cLayer.getFeatures(): print f["ID"], f.geometry().area()` still gives me the same output as before, as if the units are still in decimal degrees. What am I missing here?


title: Source code to download OSM files in qgis 2.0

text: How to get the source code to download OSM files in qgis 2.0?


title: Architecture for serving OSM-Mapnik Tiles

text: I am currently trying to set up a Mapping-Service for serving OSM-Mapnik Tiles. In my project we want to use the Mapnik-Layer as a BaseLayer for our GI-Application. I am quite new to Mapnik (only used a public Mapnik-Service as a BaseLayer in OL so far) and i am now struggling with the decision about a good service architecture for my use case. Maybe somebody could give me a hint. Currently from my point of view i have to decide between two possible service architectures: 1.) In the first solution i serve mapnik tiles via mod_tile. The OSM data lies in a PostGIS database and is updated regularly. Mod_tile in combination with rendered does than the rendering and tiling of the OSM data and serves the tiles via tileserver interface. This data could than be integrated in an OL- Application via OpenLayer.Layer.OSM. 2.) In the second possible service architecture the OSM data again lies in a PostGIS database and is updated regularly. But the data is now served via a wms interface (are there any good solutions for this?). The wms interface would then be cached by a separate mapping cache. After that the maps could be used in the OL-Application via wms or wmts interface. From my initial point of view, the second solution would better fit in my existent architecture, where i already use a mapserver which is cached by a mapcache. Also if using wms and osm layers together in Openlayers, i always get some troublesome behavior of the map extent, where the mapnik tileset define the map extent. The problem is here, that the mapnik map extent is always quite bigger than the actual visible map extent. This also speaks for a wms / cache solution, where the behavior is from my point of view a bit more like i expect it. On the other hand i am new to the mapnik world and so i am not sure if this solution does better fit to the osm data and idea. So what is your opinion and experience? What is are benefits of the first architecture compared to the second one? What are good software tools (in a debian environment) for publishing osm / mapnik data via wms interface? I hope i make my point clear and some people could share there experience with me.


title: Set projection/CRS in map and shapefile to produce buffer

text: I simply want to create a buffer of a specific size, around a single point. When I create a buffer it uses degrees/ or radians, or whatever, and the buffer is created that covers the whole screen. I am working on a South African map, so I have tried setting the project CRS to WGS84/UTM30S. I then load my shapefile, and set the shapefile CRS to WGS84/UTM30S. Still the buffering tool creates this massive buffer. I just want to be able to create a buffer, of known radius in meters. I am using QGIS 1.8.0 Lisboa Anyone with ideas?


title: MMqgis (hub distance) error in distance

text: So I am trying to calculate the distance of several points to a coastline, using hub distance (mmqgis plugin) in Qgis Dufour 2.01. I am using the default coordinate system (WGS84) and I, already, transform the land in nodes in order to not get the centroid distance. Problem is, I am not allowed to select the units. Please I need help with this. Normally I use ArcGIS but I realy want start to use freeware. (see the print-screens) ![Coordinate system](http://i.stack.imgur.com/PhAGD.png) ![Coastline nodes](http://i.stack.imgur.com/2ZAtt.png) ![Project settings](http://i.stack.imgur.com/0qQHS.png) ![hub distance attempt \(locked units\)](http://i.stack.imgur.com/Dq4MO.png)


title: Export KML from points plotted in javascript Google Maps API?

text: I want to export a KML file from points plotted in Javascript using the Google Maps API. Is there a way to do this? The user logs in and checks Google Maps, then plot points on the map and finally export those points in a KML format. **NOTE:I am not using PHP**.


title: How to convert Raster to points in MapInfo?

text: **I wanted to convert Raster data to point data in MapInfo.Does there any tool to convert? can anyone help me?**


title: How to update file locations in QGIS projects?

text: I don't know if there is an answer to this question, aside from the obvious. I have about 30 different QGIS Project files for 30 different clients. Each is saved onto my computer and all their relevant vector and raster files are all saved locally. We are now in the process of moving office and upgrading our computers. As a result, I need to move all my files from my computer onto the network drive, thus destroying any previous file references. Aside from manually clicking on each file and looking for each ones new location, is there any other way around this. (Each project could have a few hundred links to update)


title: How to run sextante algorithms outside of qgis python console?

text: I'm kind of stuck trying to figure out the way to run sextante from a standalone python from OSGeo4W distribution. The reason I want to do this is that I got tired entering parameters in the dialog every time I want to test a model from Model Builder. So here is the python script let's call it `test.py`               # as per http://qgis.org/pyqgis-cookbook/intro.html#using-pyqgis-in-custom-application     from qgis.core import *     # supply path to where is your qgis installed     QgsApplication.setPrefixPath("C:/OSGeo4W/apps/qgis", True)     # load providers     QgsApplication.initQgis()          from sextante.core.Sextante import Sextante     Sextante.alglist()     Sextante.alghelp("saga:slopeaspectcurvature")      That I'm calling from my batch file               @echo off          set OSGEO4W_ROOT=C:\OSGeo4W     set PYTHONPATH=%OSGEO4W_ROOT%\apps\qgis\python;%OSGEO4W_ROOT%\apps\qgis\python\plugins;%HOME%/.qgis/python/plugins     set PATH=%OSGEO4W_ROOT%\bin;%OSGEO4W_ROOT%\apps\qgis\bin;%OSGEO4W_ROOT%\apps\qgis\plugins          python test.py      The problem is that it says `Algorithm not found` whereas I get meaningful output from QGIS python console. I feel like I'm missing to initialize something. But what? Is there a better way to test a Model other than via entering tons of parameters using GUI? **UPDATE 7/2/2012** I'm looking for generic pythonic solution to test with "mine" algorithms. Aforementioned algorithm is just an example showing that something probably was not initialized. **UPDATE 7/27/2012** An alternative to Script Runner is to use IPython console to debug scripts. Other than that there does not seem to be a way to do simple unit testing with sextante with nothing else running:( **UPDATE 7/30/2012** As Victor Olaya suggests, I try to initialize Sextante like in the code below.               #!/usr/bin/env python          import sys     from PyQt4.QtGui import QApplication     from sextante.core.Sextante import Sextante          def main():         """ main function or something """         # as per http://qgis.org/pyqgis-cookbook/intro.html#using-pyqgis-in-custom-application         from qgis.core import *         import qgis.utils              app = QApplication(sys.argv)         # supply path to where is your qgis installed         QgsApplication.setPrefixPath("C:/OSGeo4W/apps/qgis", True)         # load providers         QgsApplication.initQgis()         # how???         # qgis.utils.iface = QgisInterface.instance()         Sextante.initialize()         run_script(qgis.utils.iface)          def run_script(iface):         """ this shall be called from Script Runner"""         Sextante.alglist()         Sextante.alghelp("saga:slopeaspectcurvature")          if __name__=="__main__":         main()      However I get something like               Traceback (most recent call last):       File "test.py", line 29, in         main()       File "test.py", line 20, in main         Sextante.initialize()       File "C:\Documents and Settings\user\.qgis\python\plugins\sextante\core\Sextante.py", line 94, in initialize         Sextante.addProvider(GrassAlgorithmProvider())       File "C:\Documents and Settings\user\.qgis\python\plugins\sextante\grass\GrassAlgorithmProvider.py", lin     e 17, in __init__         self.actions.append(DefineGrassRegionAction())       File "C:\Documents and Settings\user\.qgis\python\plugins\sextante\grass\DefineGrassRegionAction.py", li     ne 16, in __init__         canvas = QGisLayers.iface.mapCanvas()     AttributeError: 'NoneType' object has no attribute 'mapCanvas'      Well... it all becomes a mailing list discussion alike. Perhaps it worth moving to qgis-user or qgis-developer instead of SE.


title: Draw svg-circle on a coordinate with polymaps

text: I'm making a tracker system which uses polymaps to display a satellite image. The code is copied from the bing-example on the website of polymaps. I now want to periodically call an `update` function which retrieves the current location of my trackers in geoJSON format:               {     "type":"FeatureCollection",     "features":[{                 "type":"Feature",                 "id":"T_xxxxxxx",                 "geometry":{                               "type":"Point",                               "coordinates":["5.xx","51.xx"]  // 51N, 5E; Netherlands                             },                 "properties": {                               "snelheid":"5.00",                                "richting":"0.00",                               "satAantal":"5",                               "devMessage":"",                               "devBatVolt":"4.11",                               "devBatStatus":"F",                               "devBatLaden":"0",                               "dt":"2012-06-28 07:07:48","dif":"-6667"                               }                 }]     }      The update function:               function update(){       $.getJSON('/_xml/trck/current.php', function(data) {         $.each(data.features, function(key, val) {           //update the tracker information in the list               [....]           //display an svg circle                 /*                  *  But how?                  */         }        });     };      So far so good, the information does get updated with the information from the geoJSON. **Can anybody point me in the right direction how to add an svg-circle on my polymap?**


title: Using the PgRouting driving_distance function can we see the edges?

text: I have a directed graph and am using PgRouting's driving_distance function. My application needs it (or a related PostreSQL procedure) to return a full set of EDGES emanating out of the source, and not the just the VERTICES. The Drivetime function (the essential "Dijkstra algorithm") returns the following columns:   * Vertex_id   * CostToSource (from that vertex)   * Edge_id (where edge.source = vertex_id ) Of course a vertex can have several edges associated with it and only one edge per vertex is returned in this table (for some reason). Qu: Can anybody think of a way to relate the above information with the network so that I can get what I need, which is:   * Related Pair ( Vertex_ID, Edge_ID)   * CostToSource (From vertex_id) Hint: the obvious problem with joining vertex_id to the "source" in the network table is that some of the links do not belong in the result set. Any link heading back to the source should be excluded. Thanks all.


title: using mod() with IF THEN field calculation

text: I am trying to use the field calculator to assign a "odd" or "even" value to a field of small integers. I usually do this with by selecting where:               MOD([field],2) = 0      and then just calculating it as "even", then reversing the selection and setting it to "odd" I was hoping to do something like this in the field calculator:               DIM oddeven     IF mod([Field],2) = 0        oddeven = "even"     else        oddeven = "odd"     end IF      This does not work and I assume there is an alternate use of the mod() in this setting, or perhaps you cannot use an operator in the IF/THEN condition statement. Any ideas would be appreciated! Thanks.


title: How is raw OSM data processed for openstreetmap.org

text: Can anyone provide insight into how OSM data is processed or rendered for www.openstreetmap.org? A specific example... I extracted data from a recent planet.osm PostGIS dataset for an area in Missouri. The OSM data needs a lot of cleaning before it can be rendered using the correct styles. Many water bodies are stored as line strings that don't close properly, so I have to use FME for snapping and then polygon building so that I can have blue filled rivers / lakes. If I look at the same data here the water bodies are rendered as expected. I'm having trouble identifying all the cases where snapping is required (e.g. which 'Natural' types require it and what the tolerance should be). Also I suspect there are many other data issues that I will never see as I am dealing with all of North America. Does everyone who downloads and uses OSM data go through their own cleanup process? Does anyone know how this cleanup is handled by www.openstreetmap.org? It seems like their process would be the best informed and most tested. Any insight much appreciated. **EDIT** : Here is more information on my workflow A planet.osm file is downloaded and loaded into PostGIS, using Osmosis, into the pgsql schema. I then extract OSM xml from PostGIS for lots of small areas, again using Osmosis. Each of these small xml files is then converted into Shapefiles using FME and its broad feature categories. It is this stage (OSM xml -> Shp via FME) that I am expecting to convert lines into polygons and perform other cleanup on the data. These Shapefiles are served up through GeoServer (and cached using GWC).


title: How to pan to a specific feature using arcpy?

text: I am attempting to use python and arcpy to automate the creation of some standardized maps I need to export for feature classes from various projects. I have a set of points and I need to export a single map for each point in the feature class, showing only one point at a time. I can iterate through the features and export the maps just fine. However, the maps that come out are not panned to the right area, most of the time they are totally off. I am using lyr.definitionQuery to hide all the features but one, then I need to pan my data frame so that it is centered on that feature, then export the map, and repeat for all the features in the feature class. The best I could find to recenter the data frame was the df.panToExtent function. Which works fine when I run it from the python window in ArcMap but doesn't seem to keep up when it runs inside the loop. Any help is MUCH appreciated! Here is my loop:               while isoPosition < isoCount: #iterate through the feature class         query = "[Type] = 'Isolate' AND " + "[FieldID] = " + "'" + FieldIDs[isoPosition] + "'"              #center on feature         df.panToExtent(lyr.getSelectedExtent())              #Refresh the view to load new image tiles         arcpy.RefreshActiveView()              #Set the layers definition query equal the query text         lyr.definitionQuery = query              #String for a unique file name for each map, corresponding to the current Field ID in FieldIDs         filename = FieldIDs[isoPosition]              #Export the map document:         arcpy.mapping.ExportToPNG(mxd,filename)              #Move forward in the FieldIDs array, so that next time through we export the next map         isoPosition += 1


title: QGIS: Expression WHEN property1 = 'FOO' THEN property2 = 'BAR'

text: Our Foundation is adding additional properties to points. If property1 = FOO , then add property2 = BAR                   WHEN             property1 = 'FOO' THEN property2 = 'BAR'         END      When I run this the property1 column is deleted. Can you help me with the syntax, here, please? Assumptions/Knowns: This is a spatialite layer. property1 and property2 columns both exist. I am using QGIS 2.4. I opened the attributes tables to use the Expression dialog. The property1 column is selected when the expression is executed. property1 and property2 are columns, and I am aware that " is sometimes used to enclose columns. I do need to store these values for use in other programs where they cannot be so easily derived from the data. I saw several examples that changed the current value of property1: `WHEN property1 = A THEN 1` but this only changes the value in the query column, and I am looking to insert another value property2 based on the value of property1.


title: How to find polygons with area size between 300 and 500?

text: I have bunch of polyygons, all with different area-size. I want to find polygons with areasize between 300-500. Normally I would pick out everyone above 300 (query1) and then commit a new query on query1, picking out evryone below 500. But there must be an easier way to do this?Solutions in both MI and QGIS are appreciated. Thanks.


title: Open Layers Performance issue with IE 7 & 8

text: I'm using open layers to draw the world map and below is the code. Basically there's a json file which contains the coordinates of all countries which will be read by open layers and the fill color (which will be dynamic and it can be any color) of each country is got from the DB through a separate AJAX call(single call will get the color for all countries). It works fine in firefox and chrome but in IE 7 and 8, it takes a very long time to draw the map. Please let me know how to improve the performance.. Code: Layer creation:                var vecLayer = new OpenLayers.Layer.Vector("vector", {                           rendererOptions: { zIndexing: true },                           styleMap: new OpenLayers.StyleMap({                               'default': allCountryStyle,                               "select": new OpenLayers.Style({                                   strokeWidth: 3,                                   strokeColor: "#000000",                                   graphicZIndex: 2                               })                           }                           ),                           protocol: new OpenLayers.Protocol.HTTP({                               url: url,                               format: new OpenLayers.Format.GeoJSON()                           }),                           strategies: [new OpenLayers.Strategy.Fixed()]                       });      Fill Color:               var allCountryTemplate = {                               cursor: "pointer",                               fillOpacity: 1,                               fillColor: "${getAllCountryColor}",                               pointRadius: 5,                               strokeWidth: 1,                               strokeOpacity: 0.2,                               strokeColor: "black",                               graphicName: "triangle"                       };      Please let me know if you need any further information


title: Georeference vector linework

text: I'm new to GIS and a lot of searching didn't answer my specific question so far, so I hope someone can help me here. What I want to do is creating a small web application to visualize election results on a per poll site basis for my city. The result is going to be very similar to what has been done for Berlin here. Hopefully this helps in understanding what I'm going for. The election raw data on the scale I need is published and I know how to visualize this on the grounds of a shapefile with separate polygons lining out the election districts. My problem is - as you could probably guess - that something like this doesn't already exist. What actually is published though, is a PDF file lining out the borders of the districts I'm interested in: http://muenster.de/stadt/stadtplanung/pdf/a3_stimmbezirk.pdf I'm not keen on redrawing the linework myself, so I'm looking for a way to make a shapefile out of this. Luckily, the PDf includes the linework as separate vector objects, I extracted these into a SVG file: https://www.dropbox.com/s/bn7698yrdh5tdqj/a3_stimmbezirk_ungrouped.svg I already tried converting the SVG to a DXF and importing it into qgis 2.0.1. This actually works, but (as expected) places the vectors in the atlantic ocean and I can't find a way to edit it and fix this. What can I do now to georeference these vectors with open source tools and finally produce a shapefile? Any help is highly appreciated.


title: Dissolve column based on attribute in another column

text: I have 3 columns within a polygon layer I'm working with. Column 1 = Landowner, 3 different values Column 2 = Parcel, N different values Column 3 = Management type, 5 different values I would like to dissolve column 3 based on unique values in column 2, and again within column 1.   Example   C1 = A, B, C   C2 = 1, 2, 3, 4, 5   C3 = Z, Y, X   A 1 Z   A 1 Z   A 1 Y   A 2 X   A 2 Z   A 2 Z   B 1 Z   B 1 Y      Results   A 1 Z   A 1 Y   A 2 X   A 2 Z   B 1 Z   B 1 Y   Is there a tool in any version of ArcMap that I can set this up in or will I have to create something within the model builder or python? Thanks


title: Problem in saving a file in a specific output folder

text: I have a very simple script that I'm testing. It has always been working so far, but now I'm running in an unexpected error. When I try to save my files in a folder different than the environment one, it outputs this error:                   outSnowmask.save (newpath_01+ '\\'+ 'snow_'+ MODIS)     UnicodeDecodeError: 'ascii' codec can't decode byte 0xe0 in position 47: ordinal not in range(128)      My code is:               import arcpy, os, sys, shutil          newpath_01= r'C:\\Users\\unimi\\Documents\\Umberto\\Università\\PhD\\Karakoram\\MODIS_Pakistan\\Aspect\\temporary_01'     if not os.path.exists (newpath_01):         os.makedirs (newpath_01)          arcpy.env.workspace = 'C:\\Users\\unimi\\Documents\\Umberto\\Università\\PhD\\Karakoram\\MODIS_Pakistan\\MODIS_Clip\\Sample_test'     arcpy.env.overwriteOutput= True     from arcpy.sa import *          listMODIS = arcpy.ListRasters ()     print listMODIS     arcpy.CheckOutExtension ('Spatial')          for MODIS in listMODIS:         outSnowmask = ExtractByAttributes(MODIS, "VALUE = 200")         outSnowmask.save (newpath_01+ '\\'+ 'snow_'+ MODIS)


title: How to import and patch multiple ASTER files in QGIS?

text: I'm a new user and just began with GRASS, where I learned to import 30x30 resolution ASTER files and then patch them (the six I'm working with) into one map. On QGIS, it doesn't appear that there's a similar function and all six tiles are kinda staggered regarding the color alignment. THen, how do you save the patched file and change the color schemes? Lastly, how to change the lat- long to UTM? And all this, from the same GeoTIFF file?


title: Any way to extract route in ArcGIS Network Analyst into line segment with "from-to nodes"?

text: I would like to extract every from-to nodes from every line segment from a shortest-path route (produced from ArcGIS Network Analyst), so that I can calculate the relationship between buildings along the route with the point where the users standing, by the method of cross product. Is this a reasonable method(using cross product) to extract relationship between a oriented route and a building? Any other methods which can arrive the same result? ![enter image description here](http://i.stack.imgur.com/F0r5M.jpg)


title: Intersection of polygons QGIS 2.2

text: I have a polygon layer with many circle polygons in it, which are overlapping each other. I now would like to cut of the overlapping areas at their cutting edge: ![enter image description here](http://i.stack.imgur.com/A6cMF.jpg) ![enter image description here](http://i.stack.imgur.com/zq5wg.jpg) I use QGIS 2.2 and couldn't find a tool to do it, all vector tools are only working if there are two different layers to cut (there is now way to convert it to different layers because I have over 100 polygons, which are overlapping each other different times.) Is there any chance to do this? Thanks, Jan


title: Convert Autocad blocks to ArcGIS symbols

text: Im trying to import Autocad blocks into ArcGIS. The polyline part of the symbols in these blocks i want to store as a symbol in ArcGIS. Some symbols i was more or less able to find in existing symbols or edit from 2 or more existing symbols, but its cumbersome to change color for example that way. Anyone know an easy way to convert blocks to ArcGIS symbols ? Thanx


title: Migrating from PostGIS to Oracle Spatial

text: We are now using a solution based on PostGIS extension of Postgresql database for various spatial queries. But we are now thinking about data security. This is why we have a plan to migrate existing PostGIS database to Oracle spatial. I do not know whether it is good idea or not. Can anyone suggest me what will be the best for us and what kind of challenges we have to face for the migration and queries?


title: Geoserver overload

text: we've been using Geoserver v2.2.5 for quite a while and we've started experiencing difficulties recently. Currently we have about 7 000 layers published in Geoserver served via WMS. What happens is we get "exceeded memory buffer" error when accessing the map. It seems to happen randomly, unrelated to any of the published layers. What bothers me is that the problem occurs with very small number of users using the map (actually we are able to overload Geoserver with 4 concurrent PCs accessing its layers). The request is sent to Geoserver, but it is never finished and hangs somewhere between the server and Geoserver. It all ends up with Geoserver getting stuck and unresponsive (needs to be shutdown manually). We've tried rising timeout for requests (did not help), switching from OpenJDK to Oracle Java (did not help). I've run out of ideas. Has any of you ever experienced anything like this? I'm looking for solution as well as for your user experience: how many layers is published on your Geoserver, how many people access them? Thanks for any help.


title: Suitable projection for a map of Japan

text: I'm making a map of Japan in QGIS which will span UTM zones 54 and a small part of zone 53 (central and North Honshu). I'm using Natural Earth raster and shape file datasets with some WSG84 based location points plotted on, but would like to set a projection that represents the distances between points relatively accurately across the whole map (also giving an accurate scale bar) . Would just using Japan UTM zone 54 for the project CRS be suitable, or would someone be able to direct me towards a more suitable projection and how to implement this in qgis? Thanks


title: the layer attribute table is not consistent with the pointvalue

text: The point value of my grid is right, but the output value ranges from -2147483648 to 2147483647. I haved tried to recalucate statistics, but it didn't work..


title: How to calculate F_NODE and T_NODE fields in a street shapefile

text: I have developed a Streets shapefile from papermaps and using GoogleEarth. Can someone tell me a very simple method to bring values in F_NODE and T_NODE fields in my attribute table, so that I can ultimately develop Turn Tables and use my network for route optimisation?


title: SUDS (SOAP) from the Python console in QGIS

text: I´m trying to install a SOAP API called SUDS, on a Windows 7 OS, to be able to run SOAP from the QGIS Python console and pyQGIS/Python scripts. The only Python I have installed is the one comming with QGIS 2.01 standalone. How can I install an external Python library like SUDS, to be used within QGIS Python scripts and console?


title: None stop loading blank map in Qgis web Client

text: I have a similar problem if its not the same I am using windows 7 64 bit See what I have done after installing osgeo4w I check the get capability if its working and yes it is working http://127.0.0.1./qgis/qgis_mapserv.fcgi.exe?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetCapabilities I have saved qgis web client-master under this Directory C:\OSGeo4W\apache\htdocs\QGIS-Web-Client-master\site 3.I checked if its reading it and it is reading it http ://127.0.0.1./QGIS- Web-Client-master/site/ I open http ://127.0.0.1./QGIS-Web-Client- master/site/qgiswebclient.html?map=/home/web/qgis-web- client/projects/naturalearth_110million.qgs This is where I am, The map has been loading since but nothings shows (The Frame is there) I have done all the steps (edited 06 October 2012 by underdark) I went to Global options and change line 33 and 133 (var epscode=4326 and var enableGoogleCommercialMaps=false Restart Apache I still get the blank map , the frame is there with loading map but its been blank for an hour now Thank you Please help


title: Is this pgRouting perfomance acceptable for my tiny server?

text: I'm testing pgrouting functions in a tiny server with 2GB of RAM and a dual core processor. I wrote a program that create client process and they connect to the server in a specific hour to perform a query, in fact the query is a pgrouting function call. My trouble is that after 50 process, the server crash and report this error in the postgres log: > terminating connection because of crash of another server process 2012-09-05 > 15:59:09 VET DETAIL: The postmaster has commanded this server process to > roll back the current transaction and exit, because another server process > exited abnormally and possibly corrupted shared memory. 2012-09-05 15:59:09 > VET HINT: In a moment you should be able to reconnect to the database and > repeat your command. and this error occurs for every process. The table that contains the graph information has about 300,000 rows. The graph has 250 nodes aprox. My question: is it ok that this behavior occurs with a server like this and using a table with this amount of rows?


title: ArcMap adjust z values for all verticees in polylines

text: Is there way in Arc Editor (10.0) to adjust all the vertices z values by an offset? (translate in z axis?) I have many polylines which need to be lowered (eg: subtract 1m from every vertex z-value). Thanks


title: Determine if points fall upstream or downstream from a point on a line (river)

text: I have a set of fish species occurrences (lat/long) as points and set of rivers as lines. The fish species can of course only occur in the rivers (and not on land). I would like to determine whether these points fall up or downstream of a point on the river (like a dam for example). The way I think may work is to convert the river to a polygon with a small buffer around the polygon to include points that weren't exactly in the river (but should be). Then break up the new polygon to that below the dam, and that above. Then do a spatial overlay, e.g., `over(points, line_downstream)` and `over(points, line_upstream)` to get that set of points down and upstream of the dam. Is this the correct approach. If so, any tips on the code for this? I have tried converting lines to polygons and creating a buffer, but no success yet. p.s. I did notice a similar question done in ArcGIS - I need to do this in R.


title: Is it possible to create Non Earth map in openlayers?

text: I'm working on a little hobby project for a game, it has some information about locations that are non earth coordinates with extents of 99999999999999999 to -99999999999999999 might even be larger. I have a few questions regarding such large coordinates.   1. Can openlayers handle such large coordinates    2. Would it be best to scale it down to fit lat lon coords to mimic earth coords   3. what would the maximum decimals be for lat long accuracy in openlayers if I cannot use the full non earth coordinates?