title: Openlayers requests to geoserver don't produce cached tiles

text: I've created a grouped layer made from 4 MrSID rasters. I have a Openlayers map with basic Google layer. When I preview the grouped layer in geoserver , I see that tiles are being created in the GWC folder (tried both EPSG:4326 and EPSG:900919). When I view the layer in my openlayers application (on same localhost machine) , no cached tiles are being created (I've cleared the browser cache to make sure the tiles are being generated from geoserver). In my Openlayers I'm using GeoExt2 ... but still its regular WMS call. Any idea why no tiles are being generated ?!


title: How to copy a group of layers between QGIS projects?

text: In ArcGIS I am used to being able to have a couple of MXD documents open and copy-paste groups of layers from one MXD's table of contents to another. This can save a LOT of time ensuring consistent styles, but also in setting up the drawing order, pulling in the data, defining its CRS etc. for the group I would like to do the same thing in QGIS. I have two QGS documents open but there is no apparent option to save or copy a group of layers and dragging/dropping doesn't work (I get a little 'No Entry' symbol). Is this possible and, if so, how? Many thanks


title: Can I update my Goole Earth browser in real time?

text: I am about to start a project that will entail updating my Google Earth browser with data from a MySQL database package (WAMP) like is the case for web pages, but after thinking about it I really did not figure it out, so I was wondering if this is possible in real time/ I will be glad with any constructive contribution Thanks in advance. Thanks for your reply, I really realised how vague it was, thanks for pointing that out, you are awesome. For my data flow: * I will like to insert geo-coordinates on html via php to MySQL(using float INT, is there a better solution?) * Then retrieving the data and feeding it in Google Earth in My Places. So far I have done so but upon the addidtion of new data, the kml file must be loaded again before featuring in Google Earth. **My problem now is to know if an automatic browser refresh is possible, so that I won't have to load every kml file after every new data is fed into the database.** As for my SRS; I (Admin) will be using WAMP on Win 8, Php 5.3, MySQL 5.5. The users will input via HTML (frontend); and Google earth as my backend browser. The end point is I will like to be able to locate those in my network in a given perimeter with their geo-coordinates, and dynamically too. Hope I have given enough useful information, thanks in advance.


title: Hidden Features of ArcGIS

text: ArcGIS is one of the most popular software-packages available for GIS. What are the lesser-known but useful features of the ArcGIS? Which little/hidden tricks do you known that make your experience with this software better?


title: I want to edit corresponding postgresql data of a vector from openlayers- how to do it?

text: I want to edit the corresponding text/numeric data type columns of a vector through openlayres, when a user click/mouse over a vector. like, http://dev4.mapgears.com/bdga/bdgaWFS-T.html# any tutorial page for how to do it? I gone through it , it was difficult for me to understand from its js codes. -- asked in stackoverflow, but no answers. based on a comment, i posted the same question here.


title: How do you view multiple tables at the same time in ArcMap 10?

text: My original question was: > How do you undock the tables in the ArcMap attribute window? because that's how I thought it should work in the new ArcGIS 10 docking environment. Either it isn't intuitive or I am missing something obvious. I have an ArcMap user who went up to ArcGIS 10 yesterday and he wants to view two tables at the same time. If he opens two tables, they are loaded as tabs in the table window. He(we) try to drag one tab elsewhere and nothing happens. The only window control button is an 'X' which closes the table. If I right-click on the tab, **Close** is the only option.


title: Open Layer and Vector seems not working

text: I configured geoserver (2.4.4) with a layer in EPSG:4326 projection. In my client side I added the openLayer code to read and show this layer as in the following:                  map = new OpenLayers.Map ("map", {                     controls:[                         new OpenLayers.Control.Navigation(),                         new OpenLayers.Control.PanZoomBar(),                         new OpenLayers.Control.Permalink(),                         new OpenLayers.Control.ScaleLine({geodesic: true}),                         new OpenLayers.Control.Permalink('permalink'),                         new OpenLayers.Control.MousePosition(),                                             new OpenLayers.Control.Attribution()],                     maxExtent: new OpenLayers.Bounds(-20037508.34,-20037508.34,20037508.34,20037508.34),                     maxResolution: 156543.0339,                     numZoomLevels: 19,                     units: 'm',                     projection: new OpenLayers.Projection("EPSG:4326"),                     displayProjection: new OpenLayers.Projection("EPSG:4326")                 } );        .....        var kmlLayer = new OpenLayers.Layer.Vector("KML" ,{                     strategies: [new OpenLayers.Strategy.Fixed()],                     protocol: new OpenLayers.Protocol.HTTP({                         url: "http://localhost:8080/geoserver/wms/kml?layers=avl:mobiles"},                         format: new OpenLayers.Format.KML({                             extractStyles: true,                             extractAttributes: true,                             maxDepth: 2                         })                     })                 });              map.addLayer(kmlLayer);      Unfortunately I'm not able to see anything for the overlay layer. Instead if I add the folowing layer                   var avl = new OpenLayers.Layer.WMS("AVL Mobiles",                                        "http://geoserver-1:8080/geoserver/wms",                                        {                                            layers: "avl:mobiles",                                            transparent: true                                        }, {                                            opacity: 0.7,                                            singleTile: false                                        });         map.addLayer(avl);      all seems to be OK. I need to use a Vector layer insted of a WMS one, as I need to add the refresh strategy and WMS layer doesn't support it. Any idea why nothing is showed in map? Any help will be apreciate.   Thanks in advance.   I discovered that the problem could be on request OpenLayer sent to geoserver. In fact I found this exception in geoserver log:               org.geoserver.platform.ServiceException: GetMap requests must include a BBOX par     ameter. but I was not able to add BBOX parameter to my vector layer.


title: How to import PostGIS layer with geography type to QGIS?

text: I tried to import in QGIS a layer from a PostgreSQL table with a geography type column created through the extension PostGIS. I noticed that QGIS doesn't seem to like the geography column... Here the step to create the column... `ALTER TABLE traffico ADD COLUMN geog geography(POINT,4326);` then to update the column... `UPDATE traffico SET geog = ST_GeogFromText('SRID=4326;POINT(' || lon || ' ' || lat || ')');` but then when I connect QGIS to my server I see this: ![Import PostGIS layer to QGIS](http://i.stack.imgur.com/56OSw.png) There's no way I can find to tell QGIS to treat my column as a geography type instead of geometry type. Indeed if I try to upload the layer like this the point coordinates do not correspond to reality.


title: How to use ArcPy cursor to define each feature sequentially as current selection of layer?

text: I want to implement a script that runs through each feature of a feature class, zooms to the extent of the feature, sets the scale to 1:5.000 and adjacent takes a picture of the scene. The steps are in general very easy. My idea is to use the following code lines:               import arcpy.mapping     myMXD = arcpy.mapping.MapDocument("CURRENT")     df = arcpy.mapping.ListDataFrames(myMXD)[0]     layer = arcpy.mapping.ListLayers(myMXD, "", df)          # Zoom and Scale     df.zoomToSelectedFeatures()     df.scale = 5000      My Question is: How can I run through the whole list of features so that the script will zoom on each feature, set the scale back to 1:5000 and takes a picture? My first idea was to use the Search Cursor to run through the attribute table, but how can I than declare that the feature that is actually "visited" by the cursour is the selected one? Thanks for any help and suggestions. Regards


title: netCDF does not display properly in QGIS 1.8

text: I have a netCDF file of fire emissions on a regular grid of 12 km by 12 km grid cells with 200 east-west grid cells by 161 north-south grid cells centered on the continental United States. The map projection is Lambert Conformal. I am using QGIS v1.8. When I open the raster it displays as a very, very tiny raster located in mid-northern reaches of the State of Kansas. Can someone tell me what blunder I am performing?


title: DLL problem with PyQGIS

text: I am having problems with Python Shell importing the module "qgis.core". When I type "import qgis.core" the Shell ("idle.pyw") gives me this error:               Traceback (most recent call last):       File "<pyshell#0>", line 1, in <module>         import qgis.core     ImportError: DLL load failed: Impossibile trovare il modulo specificato (Trad "Impossible to find the specified module").      I have already set the environment to point the right folders following "PyQGIS cookbook" instructions. In my case, the paths are:               PYTHOPATH=C:\"QGIS_path"\apps\qgis\python;     Path=C:\"QGIS_path"\apps\qgis\bin.      Am I missing something? What are supposed to be the right libraries the Shell is unable to find? I am working with windows 7 64 bit, QGIS Lisboa (version 1.8), Python 2.7.2 (the one that is included in the QGIS package).


title: How to create a routing application with ArcGIS Server?

text: I have the task to find out what is necessary to build a routing application using ArcGIS Server. I am not a GIS beginner, but have had no experience with ArcGIS Server. As this project is at the "ideas" stage, I would like to start my research here by asking if anyone could give me some guidelines regarding what would be a sensible way to start a routing project using primarily (but not necessarily limited to) ESRI Software. Which extentions? What type/quality of data? Planning tips? Which organisations etc etc?


title: How to rename .met file with model builder which uses an iterator to rename .sid files?

text: ![iterator for renaming .sid rasters](http://i.stack.imgur.com/lvxen.png) I have a model that iterates through a folder containing .sid raster datasets, amongst other types of rasters (so I'm using a *.sid wildcard). I note, however, that ArcCatalog/Model Builder/Rename tool cannot rename the .met file of this .sid dataset, hence the renamed dataset does not have projection information. Does anyone know how to force model builder to rename the .met component of a .sid file? Possibly I'd have to save it as a Python script and add some code to that, if someone can show me. Wil


title: Good advice regarding Mobile GIS Hardware and software for outdoor tasks?

text: Could anyone offer advice regarding outdoor mobile hardware and software for GPS-photos, data capture, and textual input. We don´t want to necessarily write long text, so the ability to write into shapefile format would be optimal. Obviously the price is a big factor and we have a budget of around 2-3,000 Euro/Dollar for everything. The colleague who needs this also has an iphone...so this could also be an option. I know that this has been covered in some extent in other posts, but I also would like to know if anyone has used pure mobile GIS opensource software on Android or Apple based operating systems. The Ipad also has GPS build in and would also be a great platform for GIS mobile tasks.


title: Spatial analysis and clustering near features

text: Im working on marking behaviour in group living animals and im interested in how marking behaviour is affected by certain characteristics of the neighbouring groups. I have plotted the territories of each group from 95% density isopleths created using a-LoCoH. The marks of each individual from a group are plotted as a separate points layer. Marking behaviour does not appear to be random, and they do seem to cluster around boundaries with neighbouring groups. What I would like to investigate whether individuals mark more near some groups than near others, for example, do individuals mark less near groups that are closely related? is it possible to show that points are clustering near some features but not near others? Any help that can be provided would be greatly appreciated, Dave Humphries ##### update For demonstration, below is an image of four territories (A,B,C, and D) and the marking locations of two individuals from group A. Groups A and B are closely related. Groups C and D are unrelated to group A. ![enter image description here](http://i.stack.imgur.com/eE7TD.jpg)


title: calculating terrain elevation line

text: I found too many hits on the web to even see the right direction to go. I have to provide a small proxy server serving requests for a line of terrain elevation data between two coordinates. It is a Linux system where GDAL is installed already. There are standard raster elevation maps as inputs. The proxy needs to calculate about 1000 points/sec at runtime. I'm not that much familiar with python but what I've read so far is that GDAL/python seems the way to go. Is my feeling right to go this way? thanks for feedback Wolfgang R.


title: Styling map elements based on their properties

text: I am currently generating tiles based on OSM data using Mapnik and Python and am quite happy using the XML styling but I am at a loss as to how to vary the styling on an element by element basis. I would like to render OSM map tiles where certain elements are styled according to their attributes, for example, a roads colour is based on its length, or an area is coloured based on its size. Is this possible?


title: How to find polygon neighbors by rook and queen contiguity in a feature class?

text: I have a feature class containing several polygons. I would like to find all polygons that have a common border to a polygon of my choice. This can be done with the ISpatialFilter and esriSpatialTouches. This is known as the queen contiguity. However, the rook contiguity does not contain the polygons, which have only a point in common. It seems that this is something I can not figure out with the DE-9IM. My idea is to test all polygons which touches my selected polygon with the Intersect method from the ITopologicalOperator interface and see if the result is from the type Point. If this is the case I have to remove this polygon. Is there a better way to achieve this?


title: How can I change the size of labels and points in GeoServer?

text: I am beginner in GeoServer. I was able to do some task by my professor they want me to change different size of points in map of GeoServer example I have my database named "PLACES" then these are the data inside the database.               Name          Place     Addr:city        Baluarte                Vigan     Ayusan Norte            Vigan     Barangay 8    Village   Vigan     Chinos                  Los Baños     Casubha                 Los Baños      Then I want to enlarge the labels and points of addr:city before the small cities in the field name how can I do that is that possible?


title: How to find out the Basin Length?

text: I am doing Project on Drainage basin, to find out the Elongation ratio I need to find out Basin length, Can any one of you tell me how to find out BASIN LENGTH using ArcGIS or any other way. My data base is like below Image ![enter image description here](http://i.stack.imgur.com/8zXQa.jpg) (i just uploaded this image from internet)


title: How to create a heatmap using values?

text: How do I create a 'heatmap' in QGIS based on the values of the different points? I did do this once before with rainfall, but can't remember how!


title: Impossible to open the .qgs file I've just saved with QGIS 1.8

text: I'm using QGIS 1.8.0 on Mac OS X 10.6.8 (last stable version of QGIS I downloaded yesterday). After working yesterday and this morning on it, I'm trying to open it again, and QGIS asks me the username and password to access the PostGIS database and crashes after 10 seconds. I also have QGIS 1.7 and 1.9 installed on my computer. Here is the crash message, if this can help:               Process:         QGIS [342]     Path:            /Applications/QGIS.app/Contents/MacOS/QGIS     Identifier:      org.qgis.qgis     Version:         1.8.0 (1.8.0 ())     Code Type:       X86-64 (Native)     Parent Process:  launchd [124]          Date/Time:       2012-07-13 14:54:37.581 +0200     OS Version:      Mac OS X 10.6.8 (10K549)     Report Version:  6          Interval Since Last Report:          253852 sec     Crashes Since Last Report:           18     Per-App Interval Since Last Report:  15946 sec     Per-App Crashes Since Last Report:   7     Anonymous UUID:                      4D1290E1-470F-4017-A8E7-1EA41258D661          Exception Type:  EXC_BAD_ACCESS (SIGSEGV)     Exception Codes: KERN_INVALID_ADDRESS at 0x0000000000610163     Crashed Thread:  0  Dispatch queue: com.apple.main-thread          Thread 0 Crashed:  Dispatch queue: com.apple.main-thread     0   libdelimitedtextprovider.so     0x000000011a2edeaf QgsDelimitedTextProvider::~QgsDelimitedTextProvider() + 35     1   libdelimitedtextprovider.so     0x000000011a2ede62 QgsDelimitedTextProvider::~QgsDelimitedTextProvider() + 18     2   org.qgis.qgis_core              0x0000000103e2d5b4 QgsProviderRegistry::provider(QString const&, QString const&) + 192     3   org.qgis.qgis_core              0x0000000103e50e2b QgsVectorLayer::setDataProvider(QString const&) + 89     4   org.qgis.qgis_core              0x0000000103e5b825 QgsVectorLayer::readXml(QDomNode const&) + 463     5   org.qgis.qgis_core              0x0000000103dfc763 QgsMapLayer::readXML(QDomNode const&) + 3999     6   org.qgis.qgis_core              0x0000000103e1e815 QgsProject::addLayer(QDomElement const&, QList<QDomNode>&, QList<QPair<QgsVectorLayer*, QDomElement> >&) + 1291     7   org.qgis.qgis_core              0x0000000103e1c910 QgsProject::_getMapLayers(QDomDocument const&) + 1378     8   org.qgis.qgis_core              0x0000000103e1f908 QgsProject::read() + 3760     9   org.qgis.qgis_core              0x0000000103e1ea23 QgsProject::read(QFileInfo const&) + 77     10  org.qgis.qgis                   0x00000001000240f4 QgisApp::addProject(QString) + 140     11  org.qgis.qgis                   0x0000000100026525 QgisApp::openProject(QString const&) + 55     12  org.qgis.qgis                   0x000000010001622e QgisApp::openFile(QString const&) + 108     13  org.qgis.qgis                   0x00000001000169ac QgisApp::event(QEvent*) + 40     14  QtGui                           0x0000000104d0193d QApplicationPrivate::notify_helper(QObject*, QEvent*) + 189     15  QtGui                           0x0000000104d07dc4 QApplication::notify(QObject*, QEvent*) + 1060     16  org.qgis.qgis_core              0x0000000103da0762 QgsApplication::notify(QObject*, QEvent*) + 100     17  QtCore                          0x0000000104a4a17c QCoreApplication::notifyInternal(QObject*, QEvent*) + 124     18  org.qgis.qgis_core              0x0000000103da0ada QgsApplication::setFileOpenEventReceiver(QObject*) + 158     19  org.qgis.qgis                   0x000000010000bf52 QgisApp::QgisApp(QSplashScreen*, bool, QWidget*, QFlags<Qt::WindowType>) + 5508     20  org.qgis.qgis                   0x0000000100008f52 main + 7938     21  org.qgis.qgis                   0x0000000100006c04 start + 52          Thread 1:  Dispatch queue: com.apple.libdispatch-manager     0   libSystem.B.dylib               0x00007fff895ffc0a kevent + 10     1   libSystem.B.dylib               0x00007fff89601add _dispatch_mgr_invoke + 154     2   libSystem.B.dylib               0x00007fff896017b4 _dispatch_queue_invoke + 185     3   libSystem.B.dylib               0x00007fff896012de _dispatch_worker_thread2 + 252     4   libSystem.B.dylib               0x00007fff89600c08 _pthread_wqthread + 353     5   libSystem.B.dylib               0x00007fff89600aa5 start_wqthread + 13          Thread 0 crashed with X86 Thread State (64-bit):       rax: 0x000000011a2f34a0  rbx: 0x000000011fbe61b0  rcx: 0x00007fff710b7650  rdx: 0x00000000005e8300       rdi: 0x0000000000610163  rsi: 0x000000011fb00000  rbp: 0x00007fff5fbfdee0  rsp: 0x00007fff5fbfdec0        r8: 0x0000000000000001   r9: 0x000000011fbe67b0  r10: 0x0000000000000000  r11: 0x0000000000000015       r12: 0x000000011fbe1fd8  r13: 0x000000011fbe1f10  r14: 0x000000011fbe5180  r15: 0x000000011fbe1f50       rip: 0x000000011a2edeaf  rfl: 0x0000000000010206  cr2: 0x0000000000610163


title: Can I collect Seismic data with an array of cell phones?

text: Suppose I have an array of cell phones, each having:   * an accelerometer    * a clock   * a GPS receiver   * Windows Phone 7 (WP7) How hard would it be to write a program so that synchronizes the phones to a common clock (preferably a GPS clock)? How hard would it be to write a program that starts recording accelerometer readings when it receives text message? I would want it to collect for a duration specified in the text message, and then forward the collected data (including location) back to an address specified in the text message. How hard would it be to write a program (on a PC or a cell phone) that recieves the data, processes it to compute a "fix" (ellipse) on the source signal location, then sends that fix back to each cell phone as a text message? Alternatively, could I record audio spectrograms instead of accelerometer readings?


title: Most timely source of California Fire Map data?

text: Fire season always has me scouring the web for satellite fire map information. Unfortunately many of the locations I find are outdated or not maintained well. What would be the best site to go to for the most timely fire data?


title: Is there any problem with Service Pack 4 for FME?

text: Let me try to explain a problem. I am testing exactly the same workbench on the to computers. One has installed FME 2012 without SP, another has the same but with SP4. Tested workbenches behaviour are different. First, one without SP4 process data min. 10x faster. Other thing is difference in information about processing in log window. For example, for FME 2012 without SP:               Reading table/feature class 'name1'      Reading table/feature class 'name2'      Processed CANDIDATE 7070      Reading table/feature class 'name33'      Processed CANDIDATE 19460      Processed CANDIDATE 27470      And for FME 2012 SP4 at the same part of log:               Sorting 566 features...     Calculating inside points...     Determining uppermost containment...     Constructing donuts from 566 features...      And that takes long time. Input data are in File geodatabase API format. Any clues what should I do? Why constructing donuts and holes? I was quite satisfied with the way and processing time before. Thank you.


title: FME - process only matching name datasources

text: I'm trying to find out how to deal with the below scenario with FME 2012 on board. I have a bunch of geotiff's tiles (A) and a bunch of others geotiff's tiles (B) that I would like to mosaick together. The name of each tile from (A) match eaxactly one tile in (B). I would like the process to be run tile by tile due to the amount of the tile I have to process (roughly 13k). What whould be the scenario with FME to: - take a pair of tiles with exact names from A and B, - run a workbench, - output the result, - take another pair of tiles with exact names from A and B, ... and so on. The problem is that I actually don't know how to force FME to read a tile A then find a coresponding tile in B (not reading all tiles from B) and run a workbench for the matching pair. Thanks in advance for help with that. G.


title: How to select all parcels of a street? [ArcGIS Server]

text: I have an Arcgis Server Map Service with a `Street` layer and a `Parcel` layer, I need to get the info of each parcel of a street that is selected by the user, does anyone knows how to do that? I'm a little bit lost. I'm using Arcgis Server 9.3, Visual Studio 2008 with Web ADF for .NET, thanks.


title: Importing arcgis.rand error from python shell

text: Not sure what's going on here with this module. Maybe a pathing issue?   1. If I import the module from the python window in ArcMap it works.   2. If I import it in a script, and use that script in a scriping tool the tool works.   3. If I import it in a script and run the script, it fails. Error message says no module named arcgis.   4. If I call the script tool from a Python add-in tool it fails and gives the same error message. This worked on my computer at home but does not work on my computer at work. Any ideas? I'm using ArcGIS 10.1 and Python 2.7.2.


title: How can I get the map to repeat to fill the screen?

text: I'm programming my first ArcGIS application. I've created a simple map. When I zoom out to a world view, the map doesn't fill the width of the screen. Ideally I'd like it to repeat until the available space is filled. How can I do this? Here's the code I have so far:               var map = new esri.Map("map", {          basemap: "topo",     center: [0, 0],     zoom: 2     });      And here's the result: ![enter image description here](http://i.stack.imgur.com/exJB8.png) How can I get the map to repeat to fill the white space on either side? (Crossposted: http://forums.arcgis.com/threads/79757-How-can-I-get-the-map-to- repeat-to-fill-the-screen) **Update** : I've set the spatial reference to (I think) Web Mercator, and wrapAround180 to true, but this has no effect. New constructor:               var map = new esri.Map("map", {                 basemap: "topo",                 center: [0, 0],                 zoom: 2,                 sliderPosition: "bottom-left",                 spatialRefence: new esri.SpatialReference({ "wkid": 3857 }),                 wrapAround180: true             });      What am I doing wrong?


title: Bounded box on postgis

text: I have a few MULTIPOLYGON geometries that are very complex to work with, and I need to find POINT geometries from a ~40M rows that are inside these MULTIPOLYGON geometries.   Right now I'm working with ST_Within but things are going on extremely slowly. I don't need all of the POINT geometries, just a few of them, but they all have to be inside of the MULTIPOLYGON (I need a few points per geometry). I thought that if there was something reverse to an envelope, that is to find the largest box that is within a geometry. I couldn't find it in the docs or Google or here. Is it possible to create such a box? If so, how? I do not want to use simplify since it is:   1. Not accurate enough.   2. Still slow.


title: How to make polygon from cells of a raster?

text: I have 100 rasters, each raster has the same cell size and extent and i want to make a polygon of each cell. Can anyone help me doing this?


title: MapFish Print: does it support ArcGIS Server?

text: we are looking for a map print solution, MapFish printing module quite fits. but as far as the searches I've done, it seems MapFish doesn't support ArcGIS Server. One thing I am not sure about is, whether it supports ArcGIS Server WMS or not? wish someone in this site who knows this can give me some advice. thanks!


title: Random value calculation in QGIS

text: Is there a way to do a random value calculation within a range in QGIS? I have been able to do this within ArcMap Field Calculator with the following code:               Dim max, min     max=20     min=5     x=(Int((max-min+1)*Rnd+min))          __esri_field_calculator_splitter__     x      Thanks


title: How to store and edit data with a defined topology using open software and formats?

text: I'm looking for a way to store and edit spatial data that preserves topological information. In my concrete application, I need to represent triangular meshes, which consist of unique nodes and triangular polygons. So far, I have been working with two separate shapefiles, one for the nodes and one for the triangles. This has the drawback that the two aren't topologically connected: If I move a node, the affected polygons in the other shapefile can't be updated accordingly. Is there a standard open format for storing and editing spatial data of different types (points and polygons) that are connected topologically, similar to what I understand ESRI Coverages to be?


title: How to migrate a Geoserver setup to Tomcat?

text:
I have some problems with tomcat/geoserver. I was create project in stand- alone geoserver and all works fine. This is my code.               <%@ page language="java" contentType="text/html; charset=UTF-8"      pageEncoding="UTF-8"%>       <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"        "http://www.w3.org/TR/html4/loose.dtd">                    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">      <title> My first JSP </title>      <link href="http://localhost:8080/geoserver/openlayers/theme/default/style.css"type="text/css"/>

                But when I try to put this in tomcat nothing happens and I don't see a map. Im gonna do something with tomcat to run project which uses the Geoserver? PS im know my english is terrible.


title: Use a WMS raster in ArcGIS 10 and combine it with vector layer

text: I have connected to a WMS server which provides raster images, more specifically, LIDAR surface elevation maps. I also have a vector layer which contains the outline of buildings. I would like to measure the height of each building above the ground, and add that as an attribute to the vector layer. I have read around but can't find any information on how to work with the WMS layer beyond adding it to the map, and that you can't edit WMS layer. I have two questions:   1. Is it possible to download part/whole WMS layer so it can be edited?   2. Is it possible to extract the numerical value of the height from the raster? I'm using ArcGIS 10.1 Desktop Licence Type Advanced.


title: US Census Data - An orderly means of deducting the area-weighted population of segments of incorporated places from the population of counties?

text: I'm working in ArcMap 10.1. Currently examining the segment of states with a geographical hierarchy of   1. Counties, coextensive with state boundaries   2. Within counties, unincorporated areas and incorporated places But incorporated places are in most of these cases not restricted to one county; i.e., a given city can span multiple counties. Although these are to some extent exceptional -- most cities are entirely contained within one county -- not all are of this type. I know how to do an intersect, calculate area-weighted populations for the resulting segments of places within counties, etc. But I'm wondering if anyone could give me some feedback on how to take **an orderly approach** to deducting, from the county population numbers, the population of incorporated places lying within, either area-weighted (where the place spans more than one county), or entire: the result I'm looking for would be the population of the balance of the county (i.e., the population of unincorporated areas). By orderly approach, I mean, how should I manage and reshape the data, use database operations, etc. The process of doing what I'm talking about is pretty straightforward when I'm looking at a small set of features and can analyze them piece-by-piece, but I'm looking at a huge set of features. It _may be worth noting_ that the attribute tables for the places layers do not have a foreign key field for the counties in which they lie, which, perhaps, was obvious from the fact that they may span multiple counties. My apologies if I've set this forth in a confusing way. Thanks in advance for any help you could provide! **EDIT** : I ended up taking a different, more straightforward approach to this. See below ...


title: Is there a way to export online ESRI basemap layers to raster?

text: I have been working on polygonizing the vegetation areas in some counties from imagery. When I have the mosaic raster locally I can easily get the polygons using Image Classification Toolbar. But when it comes to using online ESRI imagery I don't know what to do. I was wondering if someone could suggest me a way to help in saving these basemap as raster locally. When I right lick on the layer it doesn't give me the Export option.


title: GeoDjango and Spatial Indexes

text: I have three questions about spatial indexes management in GeoDjango, and I cant find the answers in the doc. First, I know GeoDjango creates a spatial_index for a geometry field by default when a table is loaded. Is the spatial_index is automatically re- created when the table is updated. For instance, when new features are added to the table through the web application. Second, If a table has more than one geometry fields, is GeoDjango creating spatial index for each of them or just for the first it encounters? Third, I would like to create many spatial indexes on the same geometry table based on an id. Is it possible to do it with GeoDjango? For instance, many countries limits are stored in a single table and I would like to create a different spatial index for each country in the case each country are always queried individually. Thank you very much, any advices and information would be appreciated.


title: Copy>Paste from meter to feet coord-sys, and coordinates don't change in ArcMap 10

text: My dataframe coordinate system is in UTM with linear unit Meter. If I have a Polygon shapefile in a UTM coordinate system with a linear unit of Foot_US and I am sent a shapefile in the same UTM coordinate system with a linear unit of Meter. I reluctantly copy/paste the linear unit Meter feature into the linear unit Foot_US shapefile and the vertices' xy are identical in both the linear unit Foot_US and Meter. What am I missing. I would have figured they would be different.


title: Raster statistics along polygon edges

text: I would like to summarize raster statistics along the edge of each parcel/polygon in a shapefile. For each parcel shown below in yellow, I would like to calculate the average, min and max value of the raster along the edge bordering the polygon. There will likely be 100,000 polygons that need to be summarized. Note that the polygons will never overlay the raster, they are only adjacent and share an edge. Ideally those stats could be joined directly to the parcels. I have ArcGIS 10.2.2 with Spatial Analyst and am comfortable scripting in Python if that is a good path. ![enter image description here](http://i.stack.imgur.com/77hy5.jpg)


title: How can I apply stretch in ArcGIS 10 Image Analysis Window?

text: I need to apply stretch to ArcGIS 10.1 exported raster. The best way to accomplish this task: Image Analysis Window, because there are fine tools to equalize raster data. 1 - This is my IRS-P6 24m cell size, Landsat's ETM+ equivalent. My goal: export the raster with histogram stretch: ![Raster with ArcGIS 10 Color Balance](http://i.stack.imgur.com/geI4V.png) Full Resolution: http://goo.gl/0KlvS 2 - Now we have the exported raster. Unfortunately, the exported raster does not save equalization or color balance adjusted in the original raster data. The file was exported without histogram curves. ![Exported Raster without color balance](http://i.stack.imgur.com/9Y0QS.png) Full Resolution: http://goo.gl/PWpG0 So, the only function of the ArcGIS 10.1 Image Analysis stretch is for visualization?


title: "Window" definition in Google Maps JavaScript API v3

text: I am trying to use Google Maps Javascript API in a java application through Nashorn (basically lets you run javascript in java). Anyway when I load the API I get an error that "window" is not defined in the API. Does anyone know where I might be able to find documentation on "window" or where in the API it might be defined? I can't seem to find any documentation and any information would be greatly appreciated.


title: Change OpenLayers 3 view center

text: I'm using OpenLayers 3 to interact with some maps. I first declare my map:               map = new ol.Map({                 target: 'map',                 layers: [                   new ol.layer.Tile({                       source: new ol.source.OSM()                   })                     ],                 view: new ol.View({                     center: [0, 0],                     zoom: 2                 })             });      I've an event that triggers an action, to change my map's view center. This way, (my coordinates are in `EPSG:4326` format):               function CenterMap(lat, long) {          console.log("Lat: " + lat + " Long: " + long);          map.setView(new ol.View({                 center: ol.proj.transform([lat, long], 'EPSG:3857', 'EPSG:4326'),                 zoom: 5          }));     }      When the function runs I get this at the explorer console:               Lat: 9.0412851667 Long: -79.5658145000       But the maps goes to `[0,0]`, does anyone knows why this happen?


title: calculate least cost pathways for a region and create distance rasters to rivers

text: How do you calculate all the least cost pathways that exists for a region? I would like the output map to display all routes (and was wondering whether this can be done in one set of operations). I want to determine this by taking into account two things: slope and valley floors. I have been informed that I need a DEM raster, and distance rasters to valley floors and rivers. I have a river system shapefile (once is in a polygon and the other is in a line format) but I don't know exactly how to create these rasters. As it stands, this is what the DEM raster and the riverways shapefile look like: ![enter image description here](http://i.stack.imgur.com/RgcHG.png) I have followed the tutorials in ArcGIS help 10.1 but I just don't seem to be going anywhere. I am able to create least cost pathways between two points based solely on slope, and I have tried reclassifying via building models but have been unsuccessful in creating images that display meaningful results - I think the biggest issue is inputting or creating distance rasters to valley floors/rivers and combining that with the slope raster. Can anyone explain in layman's terms how I could create distance rasters to valley floors/rivers from a polygon or line shapefile and how this can be combined with a slope raster to create an output map displaying all existing least cost pathways for a selected study area? Please, any help, advice, suggestions or comments would be immensely appreciated! ![enter image description here](http://i.stack.imgur.com/UneSm.png) ![enter image description here](http://i.stack.imgur.com/kUA1j.jpg)


title: Time Slider and data driven pages

text: I have parcel polygon data from 1900 to 2000. Using the Time Slider tool I can animate the parcel map. I would like to parcel attribute information like grantor name as the map is populated over time. I tried data driven pages but could not get the time slider tool and the dynamic text data drven page name to sync. Can this be done? Thank you!


title: How do you organize your geodata?

text: A geodatabase doesn't provide much flexibility in organizing your data hierarchically. All feature sets live either at the top level of the database, or one level down in feature datasets. A feature dataset can't, in turn, contain another dataset. I'm interested in hearing how others organize their data, especially in environments where they must use a geodatabase and are stuck with these limitations.


title: Calculating Runoff / Streams using DEM in ArcGIS

text: I am just curious about the workflow that many of you use for runoff/stream analysis using a DEM and create a "water network" from the results in ArcGIS 10.1? I know the average rainfall for different seasons throughout the year, but that's really all I have to work with for my area. I mainly want to investigate where water will "end up" given a fairly large rain event. I am guessing you would use ArcHydro for this? Is there a different software suite or suite of tools that you would recommend for this? Thanks!


title: Polygonizer plugin halts in Windows

text: I've had some success and some problems using the Polygonizer plugin in the past. I'm using Polygonizer 2.1 in QGIS 1.8 on Windows 7. I'm aware there are sometimes problems using Polygonizer in Windows, but I was wondering if any solutions or work-arounds had come up since I researched the issue the last time. The linear layer I want to convert at present is fairly complex and consists of 18693 lines. When I use the "new method", it reaches 44% completion and just stops. I've had it running for about 3 hours now at 44%. Only rarely do I get an error message. When I use the old method, it is extremely slow, completing about 1% per hour. Last year I was able to run Polygonizer successfully on an older computer with Windows 2000, but that machine has since died! Does anyone have any suggestions? The rarely appearing error message reads: > An error has occured while executing Python code: > > Traceback (most recent call last): File > "C:/Users/George/.qgis//python/plugins\Polygonizer\PolygonizerDialog.py", > line 103, in threadFinished > > msg = QMessageBox.question(self, 'Polygonizer', 'Polygonization finished in > %03.2f seconds. \n %d polygons were crested. \n Load created layer?' % > ((self.t2 - self.t1), polyCount), QMessageBox.Yes | QMessageBox.No, > QMessageBox.Yes) > > NameError: global name 'polyCount' is not defined > > Python version: 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit > (Intel)] > > QGIS version: 1.8.0-Lisboa Lisboa, 6416f38


title: Can't perform mosaicking of rasters since upgrading to ArcGIS version 10.1

text: I just updated my ArcGIS to version 10.1. When using ArcMap, I've not been able to perform mosaicking without getting an error when I use the "Mosaic" tool in Data Management Tools -> Raster -> Raster Dataset. The error does not tell me anything. When using "Mosaic To New Raster" I don't get an error, but it also does not create the mosaic. I've checked Environments settings to see if there was a problem there, but they seem to be correct. I'm using elevation rasters for the mosaic. Is there something that I need to do since the upgrade in order to correct this problem? Thanks, Renee


title: How do I classify or spot a pixel, based on a specific reflectance values of the six Landsat Bands using Erdas?

text: I have converted the DN values to the reflectance of a Landsat image to classify the rice field by its stages. To specify the reflectance values of each stages, I have sampled the reflectance values from the converted image. However, I would like to find (spot) the rice paddies that have the same or similar reflectance values of the previous samples... How do I do this in ERDAS 9.2? I have tried the ARCGIS but seems to not have such functions in the software. In ERDAS I found 'Target Detection' in the 'Spectral Analysis' but it could only be used based on the library(correct me if I am wrong), not my own specific values. Is there anyone who have tried this kind of work before? Please Help!


title: World map shapefile including disputed areas

text: I am looking for a global shapefile that allows you do show data for all disputed areas (like Waziristan, Spratlay Islands, Kahsmir, etc.). Do you have any idea where I can find a shapefile with not only the recognized countries but all disputed territories as well? Thanks a lot for any helpful input! Stephanie


title: How do I import specific Coordinates from Namibia in QGis?

text: I have a map here with Coordinates from the Schwarzeck Map Datum Format 1903. They have the following format e.g.: y: 41 944,31 x:363 810,71 I want these point to be pictured in QGis... how is it done? Thank you in advance for your help and time! Lukas


title: Clean small artefacts from polygon shapefile

text: With some success, I've been generalising and cleaning Ordnance Survey OpenData VectorMap District woodland polygon shapefiles to make them fit for mid-scale mapping. However, the resulting shapefiles contain small 'tails' that I wish to eliminate. I've put black dots next to some of these on the attached picture. The problem is that the 'tails' are actually part of the larger polygons, so I can't get rid of them using (for example) GRASS v.clean rmarea. If anyone can help with this problem I'd be grateful for suggestions. A sample shapefile (cleaned.shp) is included in this zip archive. Nick. ![enter image description here](http://i.stack.imgur.com/u2Crg.jpg)


title: How to get the shortest distances between features on separate polygon layers?

text: Suppose I have a county-level map for the state of North Carolina, U.S. Suppose I would like to compute the shortest distance of each county to any point along the border of Georgia. How would this be done? I have looked all over this site and the internet for suggestions and cannot find much. I am very new to GIS, so I would be very grateful for any detail. Thank you very much for your help!!


title: How to keep track of changes (version control) to MXD files?

text: I am trying to use ArcMap for the first time, so please bear with me and don't down vote the question. My superior at the work place has asked me to do some research about the life cycle of an .mxd file. Is there anything called a life cycle of an mxd file, like the life cycle of servlets?


title: PostGIS Implementation of Brabyn and Skelly Road Network Travel Time Model

text: Is there a PostGIS implementation of the Brabyn and Skelly road network travel time model for the New Zealand road network? This model is outlined in: Brabyn L, Skelly C. Modeling population access to New Zealand public hospitals. Int J Health Geogr. 2002;1(1):3 and Lauder C, Skelly C, Brabyn L. Developing and validating a road travel time network for cost path analysis. The 13th Annual Colloquium of the Spatial Information Research Centre. 2001.


title: Question on ArcGIS Online credits

text: I've been browsing the new ArcGIS Online credits policy, but I don't understand very much how much money this service would cost to my company. I've read two webpages:   * http://www.esri.com/news/releases/12-2qtr/arcgis-online-will-change-how-you-think-about-mapping-and-gis.html   * http://gis-elektrika.blogspot.com.es/2012/06/arcgis-online-pricing-what-is-service.html but I haven't drawn a conclusion. In my company:   * GIS web traffic (with ArcGIS Server): 4 people, who use 0.2 Mbps of the bandwidth   * Vector information: 20 GB   * Raster information (satellite imagery): 100 GB I'm trying to guess the monthly bill, but I'm not able. Any similar experience is welcome. Thank you very much.


title: How to add EPSG:5235 to QGIS?

text: First of all I should say I am new to EPSG and QGIS world.Then I want to Georeferced Raster image to EPSG:5235 in QGIS. But in Qgis I can't found above EPSG. How I add this EPSG to Qgis ? Pls help me, Thanks


title: add description to a new polygon and save -> postgres + geoserver + openlayers -> wfs-t insert

text: actually it is possible to draw a new polygon with openlayers and save it over WFS-T-insert (geoserver wfs-t 1.1.0) in my postgres-db:               var save = new OpenLayers.Control.Button({     title: "Save Changes",     trigger: function() {         if(edit.feature) {         edit.selectControl.unselectAll();         }          var win = new Ext.Window({         title: 'Beschreibung des Fehlerreports',         html:'<br>',         html:'Bitte beschreiben Sie den aufgezeigten Fehler in der Darstellung:',              layout:'fit',         width:500,         height: 300,         modal: true,         autoScroll: true,         shadow: false,         closable: true,         buttons: [         {             scope: this,             text: 'Speichern',             handler: function() { win.hide(); saveStrategy.save(); }         },{             scope: this,             text: 'Abbrechen',             handler: function() { win.hide(); }         }]     });     win.show();         //alert('saved');     },     displayClass: "olControlSaveFeatures"     });      In the saveing-window i need 3 forms to store some information (character) together with the polygon in the DB. This is the structure of the DB and the the published layer from geoserver: DB - Postgres: http://dl.dropbox.com/u/19066498/db.jpg Geoserver - public Layer: http://dl.dropbox.com/u/19066498/geoserver.jpg Any idea what is a possible way to store these informations together with the wfs-t insert post??


title: How to create new layer within a DXF file?

text: Using OGR in Python, I'd like to create some geometries within a DXF file, but not in the "0" layer, but in new ones. I tried with:               import ogr, os, osr, sys     ds = ogr.GetDriverByName('DXF').CreateDataSource('output.dxf' )     myLayer = ds.CreateLayer('entities')     newFeature = ogr.Feature(feature_def = myLayer.GetLayerDefn())     newFeature.SetGeometryDirectly(ogr.CreateGeometryFromWkt( 'LINESTRING(0 0, 100 100)'))     myLayer.CreateFeature( newFeature )     newFeature.Destroy()     ds.Destroy()     ds = None      But I don't have any idea about how to continue and create new layers. I tried also with               layers = ds.ExecuteSQL("INSERT INTO entities (Layer) VALUES ('Layer1')")      but it does not work.   * Any suggestion?   * Is it worth using another DXF+Python library? Thank you very much.


title: Need help getting started with QgsComposition

text: I am using QgsComposition for the first time, but I fail at the first step :-(               mapRenderer = iface.mapCanvas().mapRenderer()      > line 31, in mapRenderer = iface.mapCanvas().mapRenderer() > > AttributeError: 'NoneType' object has no attribute 'mapCanvas' Thanks for help or any advice. ### PS: The Script               from PyQt4 import *     from PyQt4.QtCore import *     from PyQt4.QtGui import *     from glob import glob     from qgis.gui import *     from qgis.core import *     from qgis.utils import*          QgsApplication.setPrefixPath("C:/OSGeo4W/apps/qgis/", True)     QgsApplication.initQgis()          name="test.tiff"     fileName = name     fileInfo = QFileInfo(fileName)     baseName = fileInfo.baseName()     rlayer = QgsRasterLayer(fileName, baseName)          QgsMapLayerRegistry.instance().addMapLayer(rlayer)          mapRenderer = iface.mapCanvas().mapRenderer()     c = QgsComposition(mapRenderer)


title: Using Field Calculator to Convert Data Type In Field

text: I am trying to convert a field that is listed as "String" because it was imported into QGIS as a CSV (rural poverty population percentages by county in NC). I need to convert the percentages field from "string" to "Real" with at least one tenths decimal place. I have been using the field calculator to create a new field with a data type conversion (toreal( "Field_Name")) with the precision at 2 (for 2 decimal places I am assuming). The expression is valid and I am able to get a result as a new field. The _problem_ is that I keep getting resulting values from the adjacent column to the one that I am trying to convert. I get real numbers with 2 decimal places but from a column next to the one I am trying to convert! (This adjacent field happens to be a "string" data type (RCU codes) so I am definitely not mistaking the data in it's converted state. Has anyone had a problem with the field calculator and it selecting a different field that the one you have identified in the conversion expression?


title: Obtain virtual output directory path in geoprocessing service

text: I am publishing a geoprocessing service and I want to be able to access the virtual output directory path configured for the service. Now I know that gp services will magically return virtual paths for output parameters of data type file, but is there a way to access the virtual output directory path string directly? right now I'm doing something very hackish:               def get_results_virtual_path(resultsFilePath):         file_url = urlparse.urljoin('file:', urllib.pathname2url(resultsFilePath))         if 'directories' in file_url:             return SERVER_VIRTUAL_DIRECTORIES + file_url.split(r'directories')[1]         else:             return file_url      Its basically hard-coding the virtual path. Now I know the virtual path can be complicated by load-balancers and subdomains and such, but all I want to do is read the environment variable from the service properties.


title: Sync features from ArcGIS Online to file geodatabase

text: Using ArcGIS Desktop Advanced 10.1 SP1, ArcGIS Online for Organizations I'm new to ArcGIS Online and we just purchased an ArcGIS Online Organizational account... I'm interested in taking a feature layer (eg. fire hydrants) and publishing it on ArcGIS Online (which I'm able to do...) and allowing users/public edit the layer by adding and updating hydrants. What I'd like to do then on a regular basis is to sync or download the hydrant layer back to my file geodatabase that I have on my laptop. The only way I've been able to figure out to get data back to my geodatabase is to download the layer from ArcGIS via shapefile and copy/replace the features in my geodatabase... wondering if there is a less cumbersome way to go about this... I've referenced this post but when I download the .pkinfo as @artwork21 suggested, it only seems to reference my feature layer on ArcGIS Online, not download anything to a geodatabase as he indicated... Thanks


title: Why does my buffer take up the entire screen?

text: I am using QGIS. I am aware there is a similar question here but the solution (changing the layer CRS) is not working for me. I cannot change my CRS as I am using information for the entire Asia-Pacific region. Is there another option to create a 1000 km buffer zone from all volcanoes in the region? It is currently projected as geographic coordinates, the default WGS:84.


title: Moving Vehicle through know digital data

text: I have a set of 3D digital map data for objects within an area. I want to drive a vehicle with known dimensions (3D points for extremities) through the area without hitting anything. I would have sensors on the vehicle to work out the positions of the extremities based on vehicle movement. I know very little about GIS apart from using ArcInfo and MapInfo about 20 years ago. Are there any tools out there that I can use to work out in real time if I will hit anything? Many thanks


title: Mapnik OGCServer doesn't work if theres no active internet connection

text: Didn't find anything similar so: I have a mapnik OGCServer on an ubuntu 12.04 virtual machine (vmware player) which works fine. The problem occurs when i disconnect all network connections, when ogc server kind of hangs. On some occasion there was a fault saying something about a broken pipe. I actually noticed this on the airport where i had no connections. Any ideas what might be going on?


title: how to load gps information in my map and track way points?

text: i am using QGIS 1.8 in ubuntu 12.04 . i am working with ROBOTIC OPERATING SYSTEM(ROS) and QGIS simultaneously . My ROS system is connected to LIDAR(sensors) that will publish gps message so i could subscribe and print latitude and longitude values continuously in terminal using c++ program. My next task is to take this latitude and longitude values from terminal and plot it in qgis map. then i want to track waypoints. can anyone help to do so? i dont have any idea as i am new to QGIS. please note that everything should be done without depending to internet connection. please suggest your ideas. . . .


title: QGIS very slow I don't know what to do

text: I'm new in QGis and I already have a big Problem with the speed of this Program. Till now I have only worked with ArcGIS and have never had a problem like this. Everytime I start QGIS and add Layers it already takes a lot of time to load which I could get over, but the whorst thing about is even when I try to drag the card a little bit with the cursor it takes ages till the progam has finished loading just to move the card a little bit to the site. Is that normal? Has anybody ever had problems like this or maybe a conclusion how I can solve this problem?? As I said I have already worked a lot with ArcGIS and never had any problems like that, I was also told that QGIS is usually faster than other GIS Programs. Please help me. Thanks. Ps. I have already tried several versions of QGIS.


title: Export shapefile to KML for Google Earth - using QGIS

text: I am a real beginner with QGIS. I have shape files of my study site. I would like to have the nice google earth background. I would like to be able to see my files in Google earth, then I can zoom and see details at a fine scale. from what I understood, I need to export my files as kml files. I clicked right on the layer and clicked on "saved as", I chose a place to save and a name for the file, I chose WGS 84, and saved. When I open the kml layer in Google earth and/or QGIS it is completely wrong. In Google earth, appears some segments every where on the globe; in QGIS, I have on the right a vertical line. Can anyone explain to me what goes wrong? Thank you very much.


title: How to add non spatial sqlite table to map using python in qgis?

text: Is it possible to load pure sqlite 2.0 (without spatial extension) table into qgis using python? I would like to get the same result as in the case of using Add Vector Layer -> All files. I have already modified this example http://www.qgis.org/pyqgis-cookbook/loadlayer.html but it doesn't work in my case. Can someone help me? I use in the python console this code:               >>> db = r'c:\Program Files (x86)\Quantum GIS Lisboa\data\base.db'     >>> uri = QgsDataSourceURI()     >>> uri.setDatabase(db)     >>> schema = ''     >>> table = 'table'     >>> geom_column = ''     >>> uri.setDataSource(schema, table, geom_column)     >>> display_name = 'table'     >>> vlayer = QgsVectorLayer(uri.uri(), display_name, 'sqlite')     >>> QgsMapLayerRegistry.instance().addMapLayer(vlayer)     <qgis.core.QgsVectorLayer object at 0x0CFA08A0>      i also tried with 'spatialite' in vlayer. When I use Add vector layer i can see result table in the TOC. Using this code I can't display table in table of contents. I think the reason of my trouble is trivial, but I have no idea what I should do to see my table


title: Why do I get unexpected colon syntax error in Gemfile when I try to create dev user?

text: I run `sudo sh script/create_dev_user ${SUBDOMAIN}` and then I receive the following error:               Gemfile syntax error compile error     on line 9: syntax error, unexpected ':', expecting $end     gem 'sequel_pg',                   '1.6.3', require: 'sequel'      with a carrot below the space after "require:" indicating that the problem is there. What's the problem? I have downloaded and re-downloaded the Gemfile from the CartoDB github repo to make sure it is correct (FYI: when I run the command without `sudo`, it hangs on `\--- Creating databases`


title: ESRI javascipt agsjs TOC works in version 3.0 but not 3.5 of the API

text: Is there a way to get the TOC.js to work in version 3.5?


title: Converting shapefile to raster: Character attributes

text: I've got a shapefile that I would like to convert into a raster, but the trick is that the attribute I want a raster version of is character data. I'm familiar with geoprocessing using GDAL (in the command line) and R, so solutions using either of those two are preferred. However, I have access to QGIS if necessary, and even ArcMap. I'm using this shapefile. It contains polygons of municipalities in part of Brazil. The attribute NAME_2 is the name of the municipality, and I want a raster version of that. I know that you can't make a raster with characters, so I created a numeric version of NAME_2 in R by doing               library("rgdal")     munisCrop.shape = readOGR(dsn=getwd(), layer="munisCrop")      # Read shapefile     munisCrop.shape$NAME_2_NUM = as.numeric(munisCrop.shape$NAME_2)  # Make new attribute     writeOGR(obj = munisCrop.shape, dsn = getwd(), layer = "munisCrop", driver = "ESRI Shapefile") # Save new version of shapefile      Values in `munisCrop.shape$NAME_2_NUM` range from 1 to 787. I was hoping to then use `gdal_rasterize` on my command line (GDAL version 1.10.1 on Mac OS X 10.9.2) to make a raster version of the shapefile. Here's the code I was trying (separate lines for clarity):               gdal_rasterize -of GTiff -ot UInt32 -a NAME_2_NUM \                    -tr 0.0002777778 0.0002777778 -ot BYTE \                    -co COMPRESS=LZW -l munisCrop \                    munisCrop.shp munisCrop.tif      This completes easily enough. However, when viewing munisCrop.tif in QGIS, I see that the data range from 0 to 255. This is confirmed in R with               library("raster")     munisCrop.raster = raster("munisCrop.tif")     sort(unique(getValues(munisCrop.raster)))      0 makes sense, as NA values got converted to that. However, the maximum value of 255 is no good; remember, the maximum value of that attribute in the shapefile is 787. I think I understand why this is happening, but I'd appreciate if somebody could really explain it to me and suggest a workaround. I've tried various options in `gdal_rasterize` instead of `-ot UInt32` (for example, `BYTE`), and the result is always the same. I've also tried rasterizing it in R:               munisCrop.blankraster = raster(ext = extent(munisCrop.shp),                                     resolution = 0.0002777778,                                     crs = projection(munisCrop.shp)                                     )      munisCrop.Rraster = rasterize(x = munisCrop.shp,                                    y = munisCrop.blankraster,                                    field = "NAME_2_NUM"                                    )      But that seems like it takes forever. `gdal_rasterize` is way faster (a matter of seconds instead of... indefinite), so I'd love a solution using that. Thanks in advance!


title: Help in importing excel data

text: I have downloaded the hourly solar radiation data from NREL which is in CSV format but the problem is that the file does not have an exclusive column called latitude and longitude. There are 100s of CSV files like this and each file represents a latitude/longitude(we can say, point location).I have attached a screenshot of the file along with this post. I am not able to figure out how to import it in ArcGIS/QGIS. THe screenshot is for Bangalore city, India. The highlighted portion in RED shows the latitude longitude for the city. Can somebody please help me to import this in ArcGIS/QGIS? Screenshot of the sample file: http://tinypic.com/r/53r86a/8 Please find the same file attached here: http://www.sendspace.com/file/m9t19n This file represents hourly solar radiation data for a particular latitude longitude(In this case, Bangalore)


title: gdal polygonize in python error, blank polygon created

text: I am having trouble using the Polygonize function in python. The cookbook example for this can be found here. The relevant portion of my code is:               sourceRaster = gdal.Open('myraster.tif')     band = sourceRaster.GetRasterBand(1)     bandArray = band.ReadAsArray()     outShapefile = "polygonized"     driver = ogr.GetDriverByName("ESRI Shapefile")     if os.path.exists(outShapefile+".shp"):         driver.DeleteDataSource(outShapefile+".shp")     outDatasource = driver.CreateDataSource(outShapefile+ ".shp")     outLayer = outDatasource.CreateLayer("polygonized", srs=None)     gdal.Polygonize( band, None, outLayer, -1, [], callback=None )     outDatasource.Destroy()     sourceRaster = None      I know that the band has relevant information, here is a snippet of `bandArray`:               array([[ 4.,  4.,  3.,  3.,  3.,  2.,  2.,  2.,  2.,  3.,  3.,  3.,  3.,              3.,  3.,  3.,  3.,  4.,  4.,  4.,  4.,  4.,  4.,  4.,  4.,  4.,              4.,  4.,  4.,  4.],      When I open the attribute table in QGIS, it is empty: ![QGIS screencapture](http://i.stack.imgur.com/v00BR.png) **Edit:** The conversion works just fine in QGIS using Raster -> Conversion -> Polygonize tool Screenshot of raster to be polygonized: ![raster to polygonize](http://i.stack.imgur.com/zEsAQ.png) And screenshot of resultant conversion from QGIS tool: ![polygonized raster from QGIS tool](http://i.stack.imgur.com/ERxbA.png) I am using the Enthought distro on Windows 7, GDAL version 1.10.0-3 **The problem is that I cannot polygonize a raster in python using GDAL and the cookbook example, I can polygonize this same raster with no problem in the QGIS GUI**


title: Modelbuilder - Join - Keep matching records only - export to feature class

text: I need to join a table to a feature class in a model. Have it keep only the matching records and then export that joined feature class into a new feature class.


title: How to export from ModelBuilder as script and then run using Python. 3.4?

text: I'm currently studying for a Msc in GIS. In my assignment I have been asked to produce a model, export that model as a script file and then run that script as a stand alone in IDLE. So I have created the model in model builder and exported it as a script and have even had a go at cleaning the script. But whenever I try to run the script I get an error. I have attached a jpeg of my model and have also attached the script file. Does anyone know where I'm going wrong? Im using Arcmap 10.1 and python 3.4               # -*- coding: utf-8 -*-     # ---------------------------------------------------------------------------     # floodrisk.py     # Created on: 2014-04-21 11:34:26.00000     #   (generated by ArcGIS/ModelBuilder)     # Usage: floodrisk <River_thames> <POI> <SOA> <POI_clipped> <SOA_clipped__6_>      # Description:      # ---------------------------------------------------------------------------          # Import arcpy module     import arcpy     arcpy.env.workspace = "C:\Users\Imran\Documents\MASTER GIS\ASSIGNMENT4\Floodrisk.gdb"     arcpy.env.overwriteOutput = True          # Script arguments     River_thames = arcpy.GetParameterAsText(0)     if River_thames == '#' or not River_thames:         River_thames = "River_thames.shp" # provide a default value if unspecified          POI = arcpy.GetParameterAsText(1)     if POI == '#' or not POI:         POI = "POI.shp" # provide a default value if unspecified          SOA = arcpy.GetParameterAsText(2)     if SOA == '#' or not SOA:         SOA = "SOA.shp" # provide a default value if unspecified          POI_clipped = arcpy.GetParameterAsText(3)     if POI_clipped == '#' or not POI_clipped:         POI_clipped = "POI_clipped.shp" # provide a default value if unspecified          SOA_clipped__6_ = arcpy.GetParameterAsText(4)     if SOA_clipped__6_ == '#' or not SOA_clipped__6_:         SOA_clipped__6_ = "SOA_clipped.shp" # provide a default value if unspecified          # Local variables:     River_thames_buffer1 = "River_thames.shp"     SOA_clipped = "River_thames_buffer1.shp"     SOA_clipped__2_ = "SOA_clipped.shp"     SOA_clipped__3_ = "SOA_clipped__2_.shp"     SOA_clipped__4_ = "SOA_clipped__3_.shp"     SOA_clipped__5_ = "SOA_clipped__4_.shp"          # Process: Buffer     arcpy.Buffer_analysis(River_thames, River_thames_buffer1, "300 Meters", "FULL", "ROUND", "NONE", "")          # Process: Clip     arcpy.Clip_analysis(POI, River_thames_buffer1, POI_clipped, "")          # Process: Clip (2)     arcpy.Clip_analysis(SOA, River_thames_buffer1, SOA_clipped, "")          # Process: Add Field     arcpy.AddField_management(SOA_clipped, "AreaB", "DOUBLE", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")          # Process: Calculate Field     arcpy.CalculateField_management(SOA_clipped__2_, "AreaB", "!shape.area! ", "PYTHON", "")          # Process: Calculate Field (2)     arcpy.CalculateField_management(SOA_clipped__3_, "AreaB", "[AreaB]/1000000", "VB", "")          # Process: Add Field (2)     arcpy.AddField_management(SOA_clipped__4_, "areaP", "DOUBLE", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")          # Process: Calculate Field (3)     arcpy.CalculateField_management(SOA_clipped__5_, "areaP", "[AreaB]/ [AreaSQKM]", "VB", "")      ![enter image description here](http://i.stack.imgur.com/unACW.jpg)


I am trying to add content to the middle of a string by using `MyStr = Mid([MyField], 1, 1)` and I get an error. I've tried `MyStr = Mid([_smallx], 63, 6)` and `MyStr = Mid(["_smallx"], 63, 6)` The string is a file path totaling 76 characters and I need to add `_small` at 63 and keep the 13 characters at the end.


title: How to cite QGIS in publication?

text: What is the best way to cite QGIS in a formal publication? I mean, what text should be included in the References / Literature Cited section of a paper or report. I'm also specifically referring to the software itself, rather than the user's manual. I've looked in the user's manual and in Help->About, but I can't find an answer. Thanks in advance for the help.


title: MySQL: Query the envelope (mbr, extents) of a spatial table?

text: this question provides a solution to get the evelope (aka: mbr, aka: extents) of a spatial table in PostGIS, but I'm having trouble doing the same in MySQL. The ogr MySQL driver does this, and a solution is noted here on a Ruby forum (circa 2006), but surely this doesn't require external programming? Does anyone have a clean SQL recipe for this case? * * * **[Answer]** After getting pointed in the right direction by @GIS-Jonathan, the working query ends up looking like what follows. _CAVEAT: You may not need to change the`group_concat_max_len` value for a small table, but large tables with many features may need the larger `group_concat_max_len` buffer. It seems safe to me to just increase this value in the session as part of the solution._               # increase group_concat_max_len to be safe          SET SESSION group_concat_max_len = 104448;          # air8 is my geodata table     # shape is the geometry column          SELECT        asText(         Envelope(           GeomFromText(             Concat(Concat("geometrycollection(", Group_Concat(asText(shape))),")")           )         )       )        AS extent FROM air8;      **[Concluding Remark]** I was curious how well this would scale against larger tables and found that it starts to break down around some upper limit I don't have the time to identify. For a table with 60864 features, I had to use `SET SESSION group_concat_max_len = 10444800;`, then it worked (that's up to about 10MB). However, for a table with 242439 features, I tried 104448000 (about 100MB), then 1044480000 (about 1GB), then most rediculously 10444800000 (about 10GB). At that point I quit caring if it would eventually work because surely that violates sensible practice! **Anyway here is the take-home message: If you run this on a large dataset and it returns`null`, you may very well be hitting the upper limit.**


title: How to use the List field function for a toolbox script?

text: i´m trying to write a simple script for a toolbox, similar to the Addfield_datamanagement toolbox. The script adds a new field to a table of a feature class or a shapefile with attribute tables. So far so good               # Import system modules     import arcpy     #Workspace     arcpy.env.workspace="C:\\Users\\xxx.gdb"     # Set feature class     inputFC = arcpy.GetParameterAsText(0)     # Get List of Fieldnames     inputString =arcpy.GetParameterAsText(1)     inputString =arcpy.ValidateFieldName(inputString,arcpy.env.workspace)     fieldList =inputString.split(";")     fieldType = arcpy.GetParameterAsText(2)     for fieldName in fieldList:               arcpy.AddField_management(inputFC , fieldName, fieldType)     arcpy.AddMessage("Field generated:" + fieldName)     arcpy.AddMessage ("Script executed")      So this seems to work, but now i want to check first if the fieldname already exists before it´s created and print the existing name! I thought about the list fields together with the built-in len() function:               if len(arcpy.ListFields(fieldName(?) or inputString,?)==1:                  arcpy.AddMessage("Field already exists"+fieldName)      Not sure about the synthax of the ListFields command and how to integrate it in the code above! So, any help is welcome and appreciated!


title: QGIS Identify tool on vector Crash

text: I have built QGIS (1.7.4) on Mingw successfully and everythings seems to be alright. My own plugins are also identified and working fine. There is a crash when using identify tool/ attribute table over vector layer. On Raster identify works. GEOS version : 3.3.8 GDAL version : 1.9.2 Proj4 version : 4.8.0 I have built everything in Debug. When run with gdb ( for identify tool over vector) Program received signal SIGSEGV, Segmentation fault. 0x70fddf8f in geos::geom::Polygon::computeEnvelopeInternal (this=0x71067038) at Polygon.cpp:249 249 return Envelope::AutoPtr(new Envelope(*(shell->getEnvelopeInternal()))); (gdb) and for attribute table , it simply crashes away. Please let me know if there are any known issues before I start debugging it. Thank you. Srikanth


title: What might cause invalid syntax from ArcPy with Raster Calculator and Multiple Rasters?

text: I am having trouble for writing script in Python (ArcGIS 10.1). The task is generally divide each raster by 100. I have many of these raster files.... So what I've done is the following:               from arcpy.sa import *       arcpy.workspace="C:/Users/Vaio/.../Monthly_Original.gdb"       rasterlist=arcpy.ListRasters(p)       for k in range():          OutRaster=(Raster(rasterlist[k])/100          OutRaster.Name='p_2000_*.*'          OutTaster.Save(OutRaster.Name)      I have an Error in line 6. I am very new to any sort of programming. Thank you in advance!


title: How to use WKT in Google Map V3?

text: We have a shape with the following WKT. The problem with google map we normally extract out each of the lat and long value and put into an array and finally will do something like this. It works fine but in this WKT there are few sections in it but with openlayer I just put the whole WKT into this function var feature1 = new OpenLayers.Feature.Vector(buffer.transform(P4326,P900913), null, styleObject); it works perfectly fine without any issue. How to resolve this for google map?               polygon1 = new google.maps.Polygon({                                                 paths: lineArray,                                                 strokeColor: '#1E90FF',                                                 strokeOpacity: 0.8,                                                 strokeWeight: 2,                                                 fillColor: '#1E90FF',                                                 fillOpacity: 0.35                                           });               POLYGON((100.54157334877087 5.124795579269911,100.50169157216737 5.175633231048347,100.48893060559693 5.190782519229377,100.48883685890004 5.190920674118928,100.48041685890006 5.206620683576352,100.48035625412166 5.206766524898229,100.48032538396345 5.206921338503033,100.48032543474828 5.207079174986367,100.48035640452454 5.207233968780704,100.48041710314108 5.207379771251712,100.48050519798412 5.207510979300314,100.48061730361823 5.20762255068534,100.48074911188677 5.20771019779251,100.4808955574719 5.207770552403195,100.4810510125519 5.2078012951316595,100.48120950307525 5.2078012445564505,100.48136493833998 5.2077704026211435,100.48151134505608 5.207709954559655,100.48164309689552 5.207622223348936,100.48175513070903 5.207510580439686,100.48184314109994 5.207379316194805,100.49022348868367 5.191753262740997,100.50293939440304 5.17665749276832,100.50295708150922 5.1766357380872146,100.55158138575345 5.114653464985341,100.61784836581363 5.047494669000269,100.63328064200411 5.045760273775513,100.63344538038619 5.045724117696659,100.63359900089658 5.045654696108101,100.63373481573457 5.045555031237853,100.63384691226155 5.045429461939012,100.64478691226155 5.030239472910412,100.64485714103375 5.030122067725934,100.64490643054278 5.029994524131959,100.64493336141338 5.029860514958455,100.64619300597558 5.018766440281411,100.89601933027687 4.76550530772389,100.90397540484798 4.76251382120199,100.904125498221 4.7624388147297605,100.90425680769607 4.762334529640712,100.90436366963635 4.762205463944471,100.91498366963637 4.746255474257608,100.9150580432529 4.746115668549328,100.91510356406134 4.74596407285823,100.91511847627667 4.745806534382425,100.91510220471923 4.745649129543159,100.9150553769999 4.7454979296113295,100.91497979931266 4.745358766531834,100.9148783867677 4.745237007978542,100.91475505095255 4.745137350316068,100.91461454905823 4.745063637454637,100.91446230038925 4.7450187125847965,100.9143041773355 4.745004308511336,100.91414627886799 4.745020980816285,100.91399469529519 4.74506808642927,100.91385527335358 4.745143808431712,100.91373339069273 4.745245226138226,100.89511556261954 4.764122863341007,100.89441924808347 4.764384677730915,100.92861460308248 4.729467265814491,100.94306861105032 4.717019660578351,100.94308862499464 4.717001837872207,100.9600648087706 4.70137140183913,100.97805150444233 4.690005163726286,101.0032978733881 4.688855369346559,101.02139744198907 4.698874409253193,101.02153571050398 4.698934873115156,101.02168284847642 4.698968658695525,101.02183372382656 4.698974587578459,101.02198307411716 4.698952452968966,101.02212569010408 4.698903026905675,101.0222565974315 4.698828033332999,101.02237123013421 4.698730087971819,101.02246558989572 4.698612607085813,101.03186692696475 4.684358662732975,101.04641937919912 4.681631913397546,101.04654522151539 4.6815976608810255,101.0466638985581 4.681543682997702,101.0626338985581 4.672663691919618,101.06265873346082 4.672649300079806,101.07193873346081 4.667049305582631,101.07206478644287 4.666955391993337,101.07217019642661 4.666838917858445,101.07225096922653 4.666704296595284,101.07230404420704 4.6665566292515965,101.07232741025584 4.66640151121859,101.07232018198869 4.666244820211681,101.07080124929709 4.655762184380661,101.06942217992568 4.64481956709879,101.06939272311594 4.6446829890795245,101.0693400432311 4.64455353405661,101.06926572348885 4.644435092618359,101.06457125429455 4.638224697426946,101.06319926512629 4.632065339533069,101.0631864279256 4.632015335107842,101.05965809799591 4.619901065378956,101.05729160630031 4.611301837098225,101.05827771116607 4.608078035110399,101.05837451798986 4.608017540493808,101.058489797227 4.607909126498835,101.05858164229866 4.607780378456879,101.0586465236503 4.607636244073435,101.05868194792791 4.607482262344451,101.05868655379639 4.607324350696963,101.05863655379638 4.60660435064062,101.05861127892128 4.606452707295678,101.05855757613483 4.606308606552095,101.05847739948697 4.606177291720301,101.0583736663157 4.606063540879937,101.05825015109548 4.605971493021525,101.05811134809807 4.605904497442384,101.05796230786196 4.605864991876384,101.05780845342177 4.60585441379272,101.05765538298354 4.605873148090809,101.05750866622635 4.605920513095127,101.05737364164183 4.60599478535917,101.05725522228572 4.606093262376247,101.05715771700969 4.6062123609150225,101.05708467367768 4.6063477474016885,101.05703875007214 4.60649449560441,101.05696615166936 4.606839338106436,101.05567666227104 4.611054977816131,101.05564623136308 4.611203092992834,101.05564409128449 4.611354267067325,101.05567031736118 4.611503179055542,101.05810031736118 4.6203331764030775,101.05810357207437 4.620344668666436,101.06162634972249 4.632439867617015,101.06303073487372 4.638744658859191,101.06308581972252 4.638907454598495,101.06317427651116 4.639054911411789,101.0678446649334 4.64523344260049,101.06919782007431 4.655970431325576,101.0691998180113 4.65598518148874,101.07064329274297 4.665947058557962,101.06183355569814 4.671263283392991,101.04599114859657 4.680072345516371,101.03123062080087 4.6828380888643055,101.03107509912326 4.682883713338254,101.030931837252 4.682959345378033,101.0308065956257 4.683061943874136,101.03070441010425 4.683187383417033,101.02152658363796 4.697102476341069,101.00388255801093 4.687335579811719,101.00374644848105 4.687275831277903,101.00360166061274 4.68724194258745,101.00345309409626 4.687235060553534,100.97776309409626 4.688405061901499,100.97762450889043 4.688423433455393,100.97749117130599 4.688465325320313,100.97736706846887 4.6885294848210455,100.95913706846886 4.700049496057161,100.95902137500535 4.700138148283603,100.9420012335367 4.715809084690873,100.9275313889497 4.728270349847366,100.92748140376511 4.728317234131123,100.89084289136436 4.765729386924607,100.88857903739637 4.7665805940663155,100.86875468604576 4.764747755821863,100.86852748525929 4.7647587762658,100.8479974852593 4.768688780791709,100.84779668516134 4.768755321284555,100.83633102855727 4.774306706231912,100.82703778821273 4.77294291672017,100.82687125841765 4.772935785811107,100.7928690382914 4.774982310566,100.7619067259743 4.772427068381073,100.76173502182992 4.772431140443018,100.7615680553451 4.772471265016521,100.76141336170892 4.772545631275547,100.76127792223743 4.772650883064091,100.7611678493072 4.772782270361311,100.6687278493072 4.910582361568241,100.66865301195101 4.91072185418806,100.66860697601827 4.910873227876582,100.6685915165044 4.911030646151936,100.66864137306756 4.972952505860876,100.66009346236922 4.9913004519727435,100.64595576422197 5.007756306956835,100.64586418982786 5.007887198342606,100.64580031240318 5.008033503784978,100.64576663858664 5.0081894820473565,100.64477291047628 5.016941706842589,100.61625125341138 5.0460532463404055,100.6130793579959 5.046409728012444,100.61291638011534 5.0464453024506675,100.61276420329115 5.046513443519312,100.61262930836556 5.0466112492457205,100.59598930836556 5.0615212630263375,100.59591729333239 5.061594575864115,100.58139729333239 5.078394589533011,100.5813745376085 5.078422128797729,100.568241718084 5.095052757503959,100.5558565329245 5.107692725364091,100.55190792815084 5.111630975719314,100.55184291849078 5.11170421172008,100.55150998857583 5.1121286259010175,100.50958140376511 5.154917596322795,100.50948245300847 5.155040676578281,100.50940948200127 5.155180602156629,100.50936528873167 5.155332007765464,100.50935156773907 5.155489087927182,100.50936884513898 5.155645819584137,100.50941645844956 5.155796193046173,100.509492581994 5.155934442424661,100.50959429690371 5.156055266717915,100.5097177030391 5.1561540330704085,100.50985806853575 5.156226954412539,100.51001001124239 5.156271234669731,100.51016770509305 5.1562851759732355,100.51032510350089 5.156268243762022,100.51047617120751 5.1562210872795635,100.51061511569779 5.156145514679759,100.51073660930729 5.156044423696323,100.54157334877087 5.124795579269911),(100.89871905240886 4.762767948943637,100.90896767267513 4.752376354194918,100.90316054789211 4.76109794021234,100.89871905240886 4.762767948943637),(100.64661349655054 5.015063009244003,100.64734548944385 5.008616000844607,100.66139423577802 4.992263705850945,100.66151334893807 4.99207906958164,100.67018334893807 4.973469079183367,100.67023954646196 4.9733033059172405,100.6702584834956 4.973129353908616,100.67020868057224 4.911274134430876,100.76224763699163 4.774072078648479,100.7927932740257 4.776592927682631,100.79290874158235 4.776594210842545,100.82688510588919 4.774549247104856,100.83634221178725 4.775937080726917,100.83650305821145 4.775944533666873,100.83666219668962 4.77592007401154,100.83681331483866 4.77586467197369,100.8484073218576 4.770251155446057,100.86871970000875 4.766362818665806,100.88843722762824 4.768185776370008,100.64661349655054 5.015063009244003),(100.89174944562149 4.767110780318271,100.89241584048753 4.766860216725576,100.83228546584598 4.827825906907718,100.89174944562149 4.767110780318271),(100.61964713897322 5.045671590739228,100.64435410761305 5.020630262940911,100.64335031023455 5.029471003274368,100.6327424559458 5.044199836724733,100.61964713897322 5.045671590739228),(100.57884271020993 5.084233488822919,100.58263439073833 5.079431891466035,100.59710904740169 5.062684368787789,100.61351808034122 5.047981338205527,100.61446663282223 5.047874732959801,100.57884271020993 5.084233488822919))      Control for 2 different polygon.                               if(pl==0)                     {                           var polygon1 = new google.maps.Polygon({                                                         paths: cords,                                                         strokeColor: '#1E90FF',                                                         strokeOpacity: 0.8,                                                         strokeWeight: 2,                                                         fillColor: '#1E90FF',                                                         fillOpacity: 0.35                                                       });                         polygon1.setMap(map);                      }                      else                      {                           var polygon1 = new google.maps.Polygon({                                                         paths: cords,                                                         strokeColor: '#FFFFFF',                                                         strokeOpacity: 1,                                                         strokeWeight: 2,                                                         fillColor: '#FFFFFF',                                                         fillOpacity: 0.25                                                       });                         polygon1.setMap(map);                      }


title: Run FME workspace using python

text: ![error](http://i.stack.imgur.com/q4PCV.jpg)I'm trying to run a second workspace using python command in "Shutdown Python Script" in the first workspace using the following command. But its giving me an error. Is there anything I need to do different? I am using FME 2012 SP 2 Thanks, Sam So as suggested by Safe tech support, I created a batch file with the following code: `fme.exe c:\test2.fmw --SourceDataset_ARCGIS_LAYER C:\ScratchParcel\DevEnv.gdb\AddressPoint --DestDataset_SHAPE C:\temp` and added the following code in the shutdown python script window: `import os os.system("C:\\\temp\\\test.bat")` When I run it, although it says transaction was successful, I don't get any results i.e., the output shapefile is not created.


title: Replacement of QVariant and setAttributeMap in PyQGIS

text: I am updating a plugin for QGIS 1.8 to QGIS 2.x. I have the code below.   1. The setAttributeMap does not exist anymore. I have replaced it with fet.setAttributes and it seems to work.   2. When using fet.setAttributes({0:QVariant(fid}), I get the error: PyQt4.QtCore.QVariant represents a mapped type and cannot be instantiated. I have tried to just remove the QVariant() and used fet.setAttributes({0:fid}), but then I get the error: TypeError: QgsFeature.setAttributes(list-of-attributes): argument 1 has unexpected type 'dict' The code worked well with 1.8, so what can I do to replace QVariant or setAttributeMap? I think it could be something with changing setAttributeMap to setAttributes               vl = QgsVectorLayer("MultiLineString", layerName, "memory")     pr = vl.dataProvider()     vl.startEditing()     pr.addAttributes([QgsField("id", QVariant.Int)])     fet = QgsFeature()     addRange = northRange * math.tan(abs(math.radians(bearing)))     noPoints = ((eastRange + addRange) / projEast) + 1     for i in range(0,int(math.ceil(noPoints)), 1):         fet.setGeometry( QgsGeometry.fromMultiPolyline( [[ QgsPoint(minE + offset + (projEast * i) ,minN), QgsPoint(minE + offset + (projEast * i) - addRange,maxN) ]] ))         fet.setAttributeMap({0:QVariant(fid)})         pr.addFeatures( [ fet ] )         fid = fid + 1     vl.commitChanges()     QgsMapLayerRegistry.instance().addMapLayer(vl)      * * * I have change my code to:               vl = QgsVectorLayer("MultiLineString", layerName, "memory")     pr = vl.dataProvider()     vl.startEditing()     pr.addAttributes([QgsField("id", QVariant.Int)])     fet = QgsFeature()     fields = vl.pendingFields()     fet.setFields( fields, True )     addRange = northRange * math.tan(abs(math.radians(bearing)))     noPoints = ((eastRange + addRange) / projEast) + 1     for i in range(0,int(math.ceil(noPoints)), 1):         fet.setGeometry( QgsGeometry.fromMultiPolyline( [[ QgsPoint(minE + offset + (projEast * i) ,minN), QgsPoint(minE + offset + (projEast * i) - addRange,maxN) ]] ))         fet['id'] = fid         pr.addFeatures( [ fet ] )         fid = fid + 1     vl.commitChanges()     QgsMapLayerRegistry.instance().addMapLayer(vl)      but I still need to make an attribute named 'id' I guess before it works. If I use the _setAttributes_ function it includes a QVariant which I am not allowed to do! Or is it me who does not understand. * * * This is the final code, and it is working:               vl = QgsVectorLayer("MultiLineString", layerName, "memory")     vl.startEditing()     vl.addAttribute(QgsField("id", QVariant.Int))     fet = QgsFeature()     fields = vl.pendingFields()     fet.setFields( fields, True )     addRange = northRange * math.tan(abs(math.radians(bearing)))     noPoints = ((eastRange + addRange) / projEast) + 1     for i in range(0,int(math.ceil(noPoints)), 1):         fet.setGeometry( QgsGeometry.fromMultiPolyline( [[ QgsPoint(minE + offset + (projEast * i) ,minN), QgsPoint(minE + offset + (projEast * i) - addRange,maxN) ]] ))         fet['id'] = fid         vl.addFeatures( [ fet ] )         fid = fid + 1     vl.commitChanges()     QgsMapLayerRegistry.instance().addMapLayer(vl)


title: software architecture for creating routes?

text: I have a system that can alert staff by email if their device goes into/out of a geo-zone However, we now need to create a route, and then email if someone goes off this route Technically, how would i do this? I have no idea! Do you just check their GPS location and see if it exists within the GPS co-ordinates of a route? I'm using openlayers framework ....an afterthought Is a route a really narrow geozone?


title: Cartographic technique for detailing a 'Focus Area'

text: I'm working on updating a planning document. Whoever was in charge of cartography last time around used a great technique to identify and showcase focus areas (where to focus money and effort for maximum return on investment). I'm trying to use this same type of technique by interpreting and recreating the work-flow. ![Digital terrain](http://i.stack.imgur.com/nWxnb.jpg) As you can see, it certainly tells more of the story than a standard flat- filled polygon especially when considering the landscape from a restoration perspective. The focus area seems to be a terrain / colored hillshade -- could I recreate this with the national elevation dataset (NED) from USGS? Suggestions for a solid workflow to get me going would be great.


title: Calculate complete spatial randomness using Python

text: I currently use the csr() function from Splancs for R (http://cran.r-project.org/web/packages/splancs/splancs.pdf) to calculate complete spatial randomness (CSR) points to be used for spatial statistics. Does anyone perhaps know python code which can create CSR points? Thanks


title: Conversion of ArcGIS data to ADRG data?

text: I have maps generated using ArcGIS (.tfw, .fig, .prj, .sbn, .shp and .dbf files). I need to convert them to ADRG format. Can ArcGIS do that? If so, how?


title: How to generate thumbnails from shape file and display them on my website?

text: I tried to generate thumbnails from shape file and display them on my website(using jsp).I had been read this and this question .Very similarity questions ,but I know nothing about python... Is there any open source tools or projects to deal with it?Is there any solution or idea to implement it via ArcGIS javascript api or other solutions ?


title: Looking for CRS for NAD 83 UTM Zone 30N

text: I'm curious as to why the UTM zones are incomplete for NAD 83 in QGIS? I have a WGS84 shapefile that I would like to project to NAD83 UTM Zone 30N but can't find it in the QGIS CRS list. I see how one could create a custom coordinate system but the cut and paste I did of proj.4 code for zone 15 (simply changed it to 30) didn't work. Thanks for pointing me to proper proj.4 code or how to project this shapefile from the UK.


title: Incorrect coordinates when reprojecting NAD27 to NAD83

text: I have some historic UTM coordinates as NAD27 (556787, 7169250) which I entered into a spreadsheet and have imported into QGIS. I assign them the projection Nad27 / UTM zone 12N (EPSG: 26712) I save as... Projection Nad 83 /UTM zone 12N (EPSG: 26912) When I compare coordinates, the coordinate shift is different than what I know it should be (in particular, there is no east-west component to the transformation), and when I use a separate NTV2 transform there is a shift of ~ 30m north, 50m west. I have looked at This help article which involves projecting across UTM zones, which I am not doing. and This other help article which is solved by using a more robust version of the nad27 projection (EPSG: 26712 instead of EPSG:2027) which I am already doing Am I missing something very simple here? I am running QGIS 1.8.0 on window


title: Creating a multiple bar graph in ArcMap

text: I would like to graph my data from a shapefile. The data I wish to display is the sum of the nominal data per each year that it occurred. I have enclosed a screenshot to illustrate my idea. There is the field of data (string) I wish to graph. There are only A's and W's in the field (125 in total) and I want the W's totaled per year and the A's totaled per year as multiple bars like shown in the graph below. Then I can display those bars on a yearly basis with their respective totals. If that makes sense, please help me. Thank you. ![enter image description here](http://i.stack.imgur.com/J8F44.png)


title: Why is GDALPolygonize so much slower than ArcGIS Raster to Polygon?

text: I am attempting to polygonize a raster using GDALPolygonize() in a Python script. The script began polygonizing yesterday at 5pm and is still polygonizing now at 9:30 am. I have no clue how far along it is, but I know it is still going because when I refresh my Windows explorer I can see the file size change for the output shapefile. My raster is rather large, but I still don't expect it to take this long. My raster is 35,486 Columns and 23,682 Rows with a 1 meter cell size. It is a binary raster where a value of 1 represents data and 0 is NoData. When I polygonized in ArcGIS using Raster to Polygon in the Conversion Toolbox it took 56 seconds. The resulting shapefile is 200mb while the shapefile still being created by GDALPolygonize is still only 100mb. That makes me think GDAL is about half way done after running all night. Specs: Windows 7 64bit, 8gb RAM, GDAL 1.10 64bit, ArcGIS Desktop 10.2, 64bit Background Geoprocessing for ArcGIS Desktop, Python 2.7.3 64bit **UPDATE** Day 2 - GDALPolygonize is still running. It has gone overnight 2 nights in a row and through a whole day without completing. ArcGIS took 56 seconds.


title: gdal2tiles.py generating black tiles from .img file

text: I have a source HFA Erdas Imagine (.img) file which I'm trying to convert into tiles using gdal2tiles.py. I'm running:               gdal2tiles.py -w google -g GOOGLEKEY input_file.img      The HTML is generated correctly, and the tiles are correctly positioned. However, I'm getting black tiles instead of the actual content of the .img file: ![wrong tiles](http://i.imgur.com/cdZJULB.png) I've verified (using QGis and MapTiler) that the .img contains the correct data, i.e. the source file has raster data. Any ideas? Should I transform the .img to something else? Am I missing an option for the command?


title: IProximityOperator.ReturnDistance returns incorrect result

text: I try to implement search of nearest features in ArcObjects.NET (please, don't suggest existing tools like Near) and I use IProximityOperator.ReturnDistance to get distances between features and to select nearest one. But I've found that sometimes this method returns incorrect results. Here is test map with one point layer (in general type of geometry is arbitrary): Test Data. Look at the feature with OID = 18. You can see that nearest point is feature with OID = 40 but result of IProximityOperator.ReturnDistance for this feature is greater than for feature with OID = 22. What am I doing wrong? Test code:               var f18 = featureClass.GetFeature(18);     var f22 = featureClass.GetFeature(22);     var f40 = featureClass.GetFeature(40);          var f18AsProximityOp = (IProximityOperator)f18.Shape;          var d18_22 = f18AsProximityOp.ReturnDistance(f22.Shape);     var d18_40 = f18AsProximityOp.ReturnDistance(f40.Shape);      But d18_22 < d18_40. I use ArcGIS 10.2, .NET 3.5, Windows 7. Thank you. P.S. Sorry for my English.


title: Creating square buffer around point feature using ArcGIS for Desktop?

text: I would like to create a square buffer from a point feature but I do not understand the code that goes into it. Similar questions have been asked on the forums.esri website but that was over 10 years ago, and it did not work when I tried the code. How do I create a square buffer from a point feature?


title: QGis 2.01 Fedora 19 vs Windows SQL-Server Access

text: i just installed QGis 2.01 on Ubuntu / Fedora 19 / windows 8, configured ODBC- support for Fedora 19 and created ODBC-DSNs on the windows machine and the Fedora machine to the same SQL-Server 2012 instance where i stored some layers. When i load these layers in the windows version of QGis 2.01 everything works as expected, but in Fedora i get a strange behaviour: My point-Layers are loaded and displayed but every polygon layer i open, i get nothing or just the smallest polygons in that layer. No error message at all. The connection works perfectly with iODBC, unixodbc and freetds. Since only some small polygons appear: Is this some kind of ressource problem? Thanks in advance. Erik