title: With osm2pgrouting, how to handle relations?

text: I'm trying to create a topology of the public transport in Paris. I use a OSM export, osm2pgrouting with the correspondent mapconfig.xml and Postgis. Everything works fine, except one thing. In Paris, there are two kinds of public transportation on railways: the metro and the "RER", which is basically a metro that go further in the suburbs. I'm am able to import both networks (metro and RER), but I'm not able to connect them properly, i think because on OSM the shared stations are mapped as relations of two nodes (one station for the metro and another one for the RER while it is the same). One example can be seen there: http://www.openstreetmap.org/relation/2154932 According to : http://workshop.pgrouting.org/chapters/osm2pgrouting.html It seems that osm2pgrouting doesn't handle relations. Is it true that osm2pgrouting doesn't handle relations, or am i over- interpreting ? If it true, how can i manage to overcome this problem ? Thank you !


In what coordinate system and unit measures the coordinates returned by the Google Maps API are? I ask this because I have a page that returns the geocode lat / long for me to find in my postgis basis. But the query does not return anything, I need to know which projection should I work with my data.


title: How do you build a MongoDB spatial index for GeoJSON geometry?

text: I'm attempting to build a spatial index on a collection of features that I created from: `ogr2ogr -f GeoJSON out.json ...` and then imported into MongoDB (2.6.1). When I attempt to build a spatial index(`{"features" : "2d"}`), I receive this error:               Error creating index     WriteConcern detected an error 'geo values must be 'legacy coordinate pairs' for 2d indexes      The MongoDB glossary has this for legacy coordinate pairs: > The format used for geospatial data prior to MongoDB version 2.4. This > format stores geospatial data as points on a planar coordinate system (e.g. > [ x, y ]). I don't see how the GeoJSON doesn't meet this spec, since the coordinate pairs are also [x,y]. Is there something else in the legacy coordinate that a Mongo spatial index is expecting? And if that is the case, does this mean that you can't build a spatial index directly against GeoJSON ( _meaning that the geometry JSON has to be modified in some way to make it ingestible by MongoDB_ ). Ref: What is difference between GeoJSON & Legacy coordinate pairs in terms of mongoDb?


title: creating many shapefiles from one multi-polygonal shapefile

text: Does anyone know of a way to create a standalone shapefile from every entry (row) in the attribute table of a large set of polygons (4909 in total)? preferably using python in a way that the results could be stored as a list of values which could be assigned to a variable in an equation later.


title: Get lat/lon in degrees using PostGIS

text: I downloaded the shape file https://dati.comune.milano.it/media/rdp/comunemilano/dati/MM_FERMATE.7z, it represents the Milan's subway stations. I used ogrinfo to output the shapefile's summary and guess the SRID:               $ ogrinfo -summary MM_FERMATE.shp MM_FERMATE          INFO: Open of `MM_FERMATE.shp'       using driver `ESRI Shapefile' successful.          Layer name: MM_FERMATE     Geometry: Point     Feature Count: 100     Extent: (1506722.964600, 5027236.460200) - (1534186.769400, 5043728.155000)     Layer SRS WKT:     PROJCS["Monte_Mario_Italy_zone_1",     GEOGCS["GCS_Monte Mario",         DATUM["Monte_Mario",         SPHEROID["International_1924",6378388,297]],         PRIMEM["Greenwich",0],         UNIT["Degree",0.017453292519943295]],     PROJECTION["Transverse_Mercator"],     PARAMETER["latitude_of_origin",0],     PARAMETER["central_meridian",9],     PARAMETER["scale_factor",0.9996],     PARAMETER["false_easting",1500000],     PARAMETER["false_northing",0],     UNIT["Meter",1]]     ID: Integer (10.0)     FERMATA: String (100.0)     N_ARCHI: Integer (10.0)      Then I loaded the data in PostgreSQL:               shp2pgsql -s 3003 MM_FERMATE.shp | psql -d gis_sample      Finally I tried to convert "the_geom" field in degrees but I always got the same results. Output using st_astext function:               gis_sample=# select st_astext(the_geom) from mm_fermate where n_archi = 4;         st_astext                  ----------------------------------     POINT(1515086.3214 5037713.221)     POINT(1516846.3194 5036826.4702)     POINT(1515864.6252 5036806.3372)     POINT(1514766.6816 5034612.3064)     POINT(1513860.5034 5034993.1312)     (5 rows)      Output using st_asewkt:               gis_sample=# select st_asewkt(the_geom) from mm_fermate where n_archi = 4;             st_asewkt                       --------------------------------------------     SRID=3003;POINT(1515086.3214 5037713.221)     SRID=3003;POINT(1516846.3194 5036826.4702)     SRID=3003;POINT(1515864.6252 5036806.3372)     SRID=3003;POINT(1514766.6816 5034612.3064)     SRID=3003;POINT(1513860.5034 5034993.1312)     (5 rows)      How can I get the lat/lon in degrees? Below there's my PostGIS version:               gis_sample=# select PostGIS_Full_Version ();                         postgis_full_version                                               -------------------------------------------------------------------------------------------------------     POSTGIS="1.5.3" GEOS="3.3.3-CAPI-1.7.4" PROJ="Rel. 4.7.1, 23 September 2009" LIBXML="2.8.0" USE_STATS     (1 row)


title: QGIS hides few graduated points after reopening project

text: I have a strange problem in QGIS 2.2 (problem is also present in QGIS 2.1). When I add a point (shapefile) layer, apply graduated style and change point size, everything looks fine. Then I save project, reopen it, and some of the points are missing. When I open attribute table all points are there, but some of them are not shown on the map. Usually those point are in the middle of the map. I noticed that number of mising points depends on selected "on the fly" CRS (before save). When I remove graduated style then all points show up on the map... Zooming, refreshing view and changing projection doesn't help. ![Before save](http://i.stack.imgur.com/qAqgs.png) ![After save](http://i.stack.imgur.com/e9AKT.png)


title: Distributed & Cloud GIS Resources

text: Are there any definitive research papers, web resources, or (text)books on distributed GIS? Topics I'm interested in are:   * Parallelism of sequential spatial algorithms   * Frameworks for implementing GIS in either the cloud/cluster environment   * FOSS and commercial software for the above I have seen this question, but I want something a bit more definitive. CW status would be beneficial, as well as recommendations for tags. EDIT: I have entered a number of these terms on Google Scholar, and while I am a graduate student and deal with research papers on a regular basis many of the results are not substantive or encompassing.


title: Using QGIS, features are suddenly placed 7.5 thousand kms off correct location

text: I've been counting an muskrat dens on lakes with aerial photos by creating a feature in a separate shapefile. I checked my progress after doing about 80 small lakes and suddenly, half of the lakes were off of their correct location (7.5 km south) and half of them are in the right spot. I can't seem to remember, but I think that I wanted to confirm that the projection was correct and I confirmed the right projection of the layer. Is it possible that this confirmation changed the first half of my layers and when I continued with the rest, they were placed correctly? Is there anyway to reproject some of the features in a shapefile? Thank you


title: Can I have the "store relative pathnames to data sources" box checked as the default?

text: Is there a way to have the "store relative pathnames to data sources" box checked as the default rather than always having to go into the document properties and check the box? ArcMap 10.0.


title: What are standard scales, or zoom levels, for map applications?

text: I am working on an ArcEngine program that is restricted to a set of scales. It currently has:               1:187.5     1:375     1:750     1:1500     1:3000     1:6000     1:12000     1:24000     1:50000     1:100000      I don't think there's anything wrong with my scales, but I thought it would be interesting to know if there are any "standard" sets of scales used in web mapping or any other mapping applications.


title: How to create path based on waypoints?

text: I am trying to map flight paths and compare historic paths with current ones (so I later need to perform some calculations as distances between them). For this task I am using QGIS on OSX 10.9.1 As an input, I have the waypoints of every path that I want to connect linear to a line (see picture). ![enter image description here](http://i.stack.imgur.com/2U7iw.png) When using the “add delimited text layer” function, I can map the waypoint. However here I am stuck. I already tied the points to path plugin, but I can´t see any line in my map. I read that it might be possible by transforming the points into a multilinestring in WKT format. Unfortunately, I never worked with WKT files and do not even know how to create them. Also, there is a really long list of waypoints in the actual sample. All the work should not be done in a “dirty” way (just draw the line), since there are some calculations about distances (e.g. shortest distance between paths, etc) to be performed at the end. Thank you for every little piece of information! Regards


title: Local Copy of OpenStreetMap of my City in PostGIS Database

text: I would like to use OpenStreetMap vector data for land cover analysis in my city. Here is my goal:   1. I set up a PostGIS database with the city area extent (around 20x20 km) on my server.   2. I will download the polygon data from OpenStreetMap. Maybe I could use a command-line tool like osm2pgsql   3. The local copy will be "read-only".   4. On a daily basis or more often, I want to "synchronize" my local PostGIS copy of the OpenStreetMap data of my city with the latest online edits to ensure that my data is up-to-date. Is this possible without re-downloading all of the OpenStreetMap data from my region every time? (I would like to download the 'latest changes' only and somehow merge them with my local PostGIS OpenStreetMap copy)


title: Creating a Clickable Map with OpenLayers?

text:
My girlfriend is trying to create a webmap using OpenLayers that is similar to this one, which was done in ArcGIS Server. I thought I'd be able to help her as I have GIS and some coding experience but this is over my head. She's wanting to use a WFS-capable template and already has a very basic page set up. She'd like users to be able to edit the data. Do you have any advice, know of any tutorials, or of any published maps that are similar to the Gender Map that we could take a look at? Neither of us have any real experience with webmaps and are panicking! Thank you so much!! Her current code is as follows:               <?xml version="1.0" encoding="UTF-8"?>     <!DOCTYPE HTML PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "DTD/xhtml1-      strict.dtd">                   <title>OpenLayers map preview</title>                  <link rel="stylesheet" type="text/css"       href="http://southway.ddns.net:8888/geoserver/openlayers/theme/default/style.css"/>

location

             Click on the map to get feature info


title: How accurate is the Scale Line in OpenLayers

text: Is there any source describing the accuracy of OpenLayer's scale line? You can get an scale line by using the constructor:               map.addControl(new OpenLayers.Control.ScaleLine());      However, it seems to be very inaccurate in relation to my geodesic measurements (0.5 km - 1.5 km difference depending on zoom). Is there any info I can cite to explain this different? I'm using OpenLayers "Release 2.13.1" and OSM as base layer. I'm calculating the distance between point using WGS84 projection with 3 decimals and geopy distance module. As an example, the distance between the points -2.416,55.703 and -2.434,55.698 is 2066.599 m according to my calculations and 2.5 - 3 km according to the scale line depending on zoom. Ideally, the info I'm looking for is something like "the scale is accurate in +- 1 km because it calculates the distance using X method". The official doc (http://dev.openlayers.org/docs/files/OpenLayers/Control/ScaleLine-js.html) doesn't mention accuracy.


title: How to improve OSM layer resolution in OpenLayers plugin while exporting?

text: I'm having trouble with the OpenLayers Plugin exporting a map to a png or pdf file: The output resolution of the exported osm-layer is very, very poor. Zooming to 100% in the png file shows big pixels. The rest of the map - so the output from shapes and fonts - is in the usual high resolution e.g. 300 or 600 dpi. Is there any way to improve the osm-background-resolution? Best regards desputin


title: How to create infrastructure coverage using pgRouting?

text: I'm trying to reproduce what Underdark did here.   I succesfully installed pgrouting and the two required function (nearest neighbour and driving_distance).   cdv_osm_v contains all the openstreetmap points of my region, from which I select only transportation infrastructure.   taa_2po_4pgr and taa_2po_v contain network and nodes extracted by osm2po for my region.   When I try to run this query in pgAdmin3 I got a syntax error at or near "FOR".   The first select works well, and also driving_distance outputs results fine if called separately.   I'm a noob at postgresql/postgis and I don't know how to pass a geometry to nn() to test it separately (and maybe x,y and factor parameters are not so good).   What am I doing wrong?               CREATE TEMPORARY TABLE nodes_costs (vertex_id integer,edge_id integer,cost double precision);          FOR i IN       SELECT the_geom       FROM cdv_osm_v       WHERE highway like 'bus_stop' OR railway like 'station' or amenity like 'bus_station'     LOOP       INSERT INTO nodes_cost       driving_distance(         'SELECT id, source, target, km AS cost         FROM taa_2po_4pgr',         mm(i.the_geom ,0.0001,2,2,'taa_2po_v','id','geom_vertex'), 0.4, false, false);     END LOOP;          SELECT *      FROM taa_2po_v      JOIN        SELECT vertex_id, min(cost) AS costo       FROM nodes_costs     ON     taa_2po_v.id = nodes_costs.vertex_id;      I forgot: posgresql 9.1, postgis 2.0.1, pgrouting 1.0.5 from kyngchaos.


title: Why are JPG files invisible in Georeferencer?

text: I want to georeference a jpg file (map) with QGIS 1.8.0 on a Windows 7 (64 bit) machine. The georeference window opens, but the map is not visible. When converting the jpg file into a tif file (e.g. with gimp) everything works. I tried it with several files. Always the same: The tif files are visible and can easily be georeferenced, while the jpg-files are not visible in the georeference window. Any ideas? Thanks for answers, Ulli


title: Produce a list of polygons that a line crosses

text: I have a polyline feature class and a polygon feature class. What I want to do is add a field to the lines that lists all polygons that touch the line. This is effectively performing a select-by-location finding the polygons that touch each individual line then add the names of each of these to an attribute. Does anyone know how I can do this in ArcGIS Desktop or FME? My Python is not great but I can give it a go if that's the only way.


title: How to upload a shapefile to geoserver using cURL?

text: I need to automatically (using python) upload a shapefile into my geoserver. I'm trying to do as described in here but i'm getting this error:               java.io.IOException: Wrong magic number, expected 9994, got 1347093252* Connection #0 to host localhost left intact      I'm sending this:               curl -u admin:geoserver -v -XPUT -H 'Content-type: application/zip' --data-binary @C:/x/test_data/localfile.zip http://localhost:8081/geoserver/rest/workspaces/test_VN/datastores/maisum/file.shp


title: Using Torque to show rises and falls?

text: I'm looking at a handful of different mapping projects that I'd like to do something Torque-like with, but I'd like maps to stay on the map for a limited period of time. Some fake data:               Address                   Start Year         End Year     7706 Heather Limits        2002               2006     5938 Noble Shadow Inlet    2002               2003     9128 Round Pioneer Link    2003               2005     7176 Crystal Via           2004               2005     4475 Harvest Square        2004               2007      Think of a map of buildings participating in a rent subsidy program. One building enters the program in 1988 and exits in 2008. A dozen enter the program in 2002, and will exit in 2032. It seems like I can add points to the map based on a single date, but there's not an easy way to take them off on a specific day, too. I'd have to turn the "cumulative" option off and then create basically a row for each site for each year that it should be on. Something like               TorqueYear   Address                   Start Year         End Year     2002      7706 Heather Limits        2002               2006     2002      5938 Noble Shadow Inlet    2002               2003     2003      7706 Heather Limits        2002               2006     2003      5938 Noble Shadow Inlet    2002               2003     2003      9128 Round Pioneer Link    2003               2005     2004      7706 Heather Limits        2002               2006     2004      9128 Round Pioneer Link    2003               2005     2004      7176 Crystal Via           2004               2005     2004      4475 Harvest Square        2004               2007     2005      7706 Heather Limits        2002               2006     2005      9128 Round Pioneer Link    2003               2005     2005      7176 Crystal Via           2004               2005     2005      4475 Harvest Square        2004               2007     2006      7706 Heather Limits        2002               2006     2006      4475 Harvest Square        2004               2007     2007 ...      Am I right about that or am I missing something?


title: Producing building shadows?

text: How would you generate building shadows from building footprints?   I have attribute values for the heights of buildings.   So would preferably like to have the shadows' sizes vary based on that.   The best suggestion I've found using Arc is drop shadows: http://blogs.esri.com/Support/blogs/mappingcenter/archive/2008/10/21/3d-effect- for-a-building-footprints.aspx


title: How to work with shortest_path_shooting_star() in pgRouting?

text: It's been a while we are using PostGIS to store our Geometry information in our Postgres database. We have developed vertex-to-vertex routing ourselves using iGraph library. Yesterday I read about pgRouting and I am porting our current application to pgRouting. However, I face some issues in using shooting start path finder. We stored edges and vertices information in tables and for finding a graph route with turn restriction, we are calling following SQL script:               SELECT  shortest_path_shooting_star(         'SELECT e.id,                     e.first as source,                     e.second as target,                     ST_Length(e.geom) as cost,                      st_x(st_pointn(e.geom, 1)) AS x1, st_y(st_pointn(e.geom, 1)) AS y1,                     st_x(st_pointn(e.geom, 2)) AS x2, st_y(st_pointn(e.geom, 2)) AS y2,                     (not connected)::boolean::int * 1000000 AS to_cost,                     c.first AS rule             FROM edges e LEFT OUTER JOIN edge_connections c ON e.id = c.second             WHERE e.group_id = 0             ORDER BY e.id',         1209560, 1209653, true, false);      My inner subqery woud have following results:                 id    | source  | target  |   cost |  x1   |  y1   |  x2   |  y2   | to_cost |  rule       --------+---------+---------+----+------------------+------------------+---------+---------     1087134 |  926686 |  926687| 2.3299 |51.4675|35.728 |51.4675|35.7280|         |         1209706 | 1039731 | 1039870 | 4.005 |51.4082|35.7239|51.4082|35.7239| 1000000 | 1209564     ...      In the results, when a constraint exists, rule and to_cost would not be NULL. When a turn restriction exists, the results would be exactly like what specified in the docs. However, when I call the `shortest_path_shooting_star` function, results would not differ if a `rule`-and-`to_cost` exists or not. I have changed directed and `has_reverse_cost` parameters but nothing happened. I was wondering if you guys could help me in this matter. Any advice or even a clue would be very much appreciated. Thanks, * Roozbeh


title: MapServer: apply a symbol to certain data?

text: I have a shp file which contains some points. Each point has a "POINTTYPE" attribute which contains text. Is it possible to apply a different symbol basing on this attribute? For example, if pointtype is "streetlamp" use a square icon, if it is "tree" use another one, etc...? Thanks EDIT: Here is the final mapfile               LAYER         NAME bla         TYPE POINT         GROUP bla         DATA "bla.shp"         STATUS OFF         CLASS           EXPRESSION ('[pointtype]' = 'streetlamp')           STYLE              SYMBOL 'circle'             SIZE 8             COLOR 255 0 170             OUTLINECOLOR 0 0 0             ANTIALIAS TRUE           END         END         # Other classes can go here, they will cumulate     END


title: WMS weather data for Europe

text: I am looking for a source of weather data for Europe. I use http://mesonet.agron.iastate.edu/cgi-bin/wms/goes/conus_vis.cgi? and http://mesonet.agron.iastate.edu/cgi-bin/wms/nexrad/n0r.cgi? for the US, but have been unable to find a source for Europe. I saw a similar question here: List of available online WMS services (Weather, Land Data, Place Names...), but it does not answer my question. Thanks


title: How to launch a tool when my button has been pressed

text: I use ArcGIS 10.0 and ArcGIS .net SDK. I created AddIn project which contains Tool:               <Tool id="AddSearchLocationPoint"       class="AddSearchLocationPoint"       category="Add-In Controls"       image="Images\insert.png" />      and DockableWindow:               <DockableWindow id="SearchDockableWindow"                     class="SearchDockableWindow"                     image="Images\SearchDockableWindow.png">       <InitialPlacement height="360" width="570" state="pinned" position="right" />     </DockableWindow>      The dockableWindow contains WinForm’s HostElement which displays my WPF UserControl. Now I want to add a button to the WPF UserControl which will launch my tool. Do you know a way to raise the event which launches the tool?


title: ArcGIS Desktop 10 - MyBase.UpdateBitmap()

text: I am changing one icon on my extension toggle toolbar in runtime using: MyBase.UpdateBitmap(New Bitmap(Me.GetType(), bitmapResourceName)) This works fine but I get a border/highlight around my image. What can I do to prevent that? I am using ArcGIS Desktop version 10.


title: Secure online GIS for sensitive data?

text: My company I started working for is looking to move from ArcInfo and MapPoint to an online GIS for easier data access/availability. Obviously the GIS capabilities needed are rather basic, but security is a must. What sort of options are out there? In order to ensure security of data, would we need a hosted platform? Any help would be appreciated!


title: Speeding up my Python calculated field

text: I am new to Python and have started to create scripts for ArcGIS workflows. I am wondering how I can speed up my code to generate an "Hours" double numeric field from a timestamp field. I start with a track point log (breadcrumb trail) shapefile generated by DNR Garmin, with an LTIME timestamp field (a text field, length 20) for when each trackpoint record was taken. The script calculates the difference in Hours between each successive timestamp ("LTIME"), and puts that into a new field ("Hours"). That way I can go back and sum up how much time I spent in a particular area/polygon. The main part is after the `print "Executing getnextLTIME.py script..." ` Here is the code:               # ---------------------------------------------------------------------------     #      # Created on: Sept 9, 2010     # Created by: The Nature Conservancy     # Calculates delta time (hours) between successive rows based on timestamp field     #     # Credit should go to Richard Crissup, ESRI DTC, Washington DC for his     # 6-27-2008 date_diff.py posted as an ArcScript     '''         This script assumes the format "month/day/year hours:minutes:seconds".         The hour needs to be in military time.          If you are using another format please alter the script accordingly.          I do a little checking to see if the input string is in the format         "month/day/year hours:minutes:seconds" as this is a common date time         format. Also the hours:minute:seconds is included, otherwise we could          be off by almost a day.              I am not sure if the time functions do any conversion to GMT,          so if the times passed in are in another time zone than the computer         running the script, you will need to pad the time given back in          seconds by the difference in time from where the computer is in relation         to where they were collected.          '''     # ---------------------------------------------------------------------------     #       FUNCTIONS     #----------------------------------------------------------------------------             import arcgisscripting, sys, os, re     import time, calendar, string, decimal     def func_check_format(time_string):         if time_string.find("/") == -1:             print "Error: time string doesn't contain any '/' expected format \                 is month/day/year hour:minutes:seconds"         elif time_string.find(":") == -1:             print "Error: time string doesn't contain any ':' expected format \                 is month/day/year hour:minutes:seconds"                  list = time_string.split()             if (len(list)) <> 2:                 print "Error time string doesn't contain and date and time separated \                     by a space. Expected format is 'month/day/year hour:minutes:seconds'"               def func_parse_time(time_string):     '''         take the time value and make it into a tuple with 9 values         example = "2004/03/01 23:50:00". If the date values don't look like this         then the script will fail.      '''         year=0;month=0;day=0;hour=0;minute=0;sec=0;         time_string = str(time_string)         l=time_string.split()         if not len(l) == 2:             gp.AddError("Error: func_parse_time, expected 2 items in list l got" + \                 str(len(l)) + "time field value = " + time_string)             raise Exception          cal=l[0];cal=cal.split("/")         if not len(cal) == 3:             gp.AddError("Error: func_parse_time, expected 3 items in list cal got " + \                 str(len(cal)) + "time field value = " + time_string)             raise Exception         ti=l[1];ti=ti.split(":")         if not len(ti) == 3:             gp.AddError("Error: func_parse_time, expected 3 items in list ti got " + \                 str(len(ti)) + "time field value = " + time_string)             raise Exception         if int(len(cal[0]))== 4:             year=int(cal[0])             month=int(cal[1])             day=int(cal[2])         else:             year=int(cal[2])             month=int(cal[0])             day=int(cal[1])                hour=int(ti[0])         minute=int(ti[1])         sec=int(ti[2])         # formated tuple to match input for time functions         result=(year,month,day,hour,minute,sec,0,0,0)         return result               #----------------------------------------------------------------------------          def func_time_diff(start_t,end_t):         '''         Take the two numbers that represent seconds         since Jan 1 1970 and return the difference of         those two numbers in hours. There are 3600 seconds         in an hour. 60 secs * 60 min   '''              start_secs = calendar.timegm(start_t)         end_secs = calendar.timegm(end_t)              x=abs(end_secs - start_secs)         #diff = number hours difference         #as ((x/60)/60)         diff = float(x)/float(3600)            return diff          #----------------------------------------------------------------------------          print "Executing getnextLTIME.py script..."          try:         gp = arcgisscripting.create(9.3)              # set parameter to what user drags in         fcdrag = gp.GetParameterAsText(0)         psplit = os.path.split(fcdrag)              folder = str(psplit[0]) #containing folder         fc = str(psplit[1]) #feature class         fullpath = str(fcdrag)              gp.Workspace = folder              fldA = gp.GetParameterAsText(1) # Timestamp field         fldDiff = gp.GetParameterAsText(2) # Hours field              # set the toolbox for adding the field to data managment         gp.Toolbox = "management"         # add the user named hours field to the feature class         gp.addfield (fc,fldDiff,"double")         #gp.addindex(fc,fldA,"indA","NON_UNIQUE", "ASCENDING")              desc = gp.describe(fullpath)         updateCursor = gp.UpdateCursor(fullpath, "", desc.SpatialReference, \             fldA+"; "+ fldDiff, fldA)         row = updateCursor.Next()         count = 0         oldtime = str(row.GetValue(fldA))         #check datetime to see if parseable         func_check_format(oldtime)         gp.addmessage("Calculating " + fldDiff + " field...")              while row <> None:             if count == 0:                 row.SetValue(fldDiff, 0)             else:                 start_t = func_parse_time(oldtime)                 b = str(row.GetValue(fldA))                 end_t = func_parse_time(b)                 diff_hrs = func_time_diff(start_t, end_t)                 row.SetValue(fldDiff, diff_hrs)                 oldtime = b                  count += 1             updateCursor.UpdateRow(row)             row = updateCursor.Next()              gp.addmessage("Updated " +str(count+1)+ " rows.")         #gp.removeindex(fc,"indA")         del updateCursor         del row          except Exception, ErrDesc:         import traceback;traceback.print_exc()          print "Script complete."


title: Extracting Lat/Lng from Shapefile using OGR2OGR/GDAL

text: I'm trying to extract lattitude/longitude values for the streets in this shapefile: http://gispub02.sfgov.org/website/sfshare/catalog/stclines_streets.html I've seen it reommended in many places that I use this command to do so: _(Below does work! Spelling error fixed.)_ > ogr2ogr -f CSV output.csv STCLINES_STREETS.shp -lco GEOMETRY=AS_XYZ However, when I do so, I get a CSV containing a lot of useful information about the streets -- names, neighborhoods, etc -- but no actual lat/long values. What am I doing wrong? I'm a GDAL noob so am sure I'm missing something obvious.


title: How to access an ArcGIS 9.3 geodatabase from 9.2?

text: If I'm using ArcGIS 9.2 is it possible to access a feature class from a geodatabase that was created by ArcGIS 9.3? If not are there any other ways to access the data?


title: Scale in printed map (mxd) is different in code to manual printing

text: _This question is related in part to aprevious question I have asked._ I have tried printing an A3 mxd with the Print Active View code as published by ESRI and have changed only one line of code and that is, to change it to A3 size. docPaper.FormID = 8; This forces it to a size of A3 instead of the default 0 Letter size (esriPageFormLetter). The map prints out fine(ish) but it is not to scale i.e. the map is supposed to be 1:10000, but the space between the graticule grids is 9.7cm, not 10cm as it should be. If I print the EXACT same mxd out, but via the manual method i.e. in ArcMap itself, File -> Print. The scale is correct and it measures 10cm. Important to note that I have not changed any of the Page and Print Setup settings between the two prints. Why would it do this? See pic below to help understand whats going on. ![enter image description here](http://img819.imageshack.us/img819/1959/scalemapcomparison.jpg)


title: Some objects missing after importing a .shp into PostGIS

text: I have a .shp file which I am trying to import into a PostGIS database using the "PostGis SHP and DBF Loader". This works smoothly but some (maybe 5%) of the records have an empty geometry field. (I have double checked the original files and there should be geometry on every row). I can't see anything special about these rows - it appears random but presumably isn't. Can anyone tell me why this might happen?


title: vb.net ArcObjects is there a function the same as "Select Features by Polygon" of ArcMap's tool

text: I want to know whether there is a function provided by ArcObjects (vb.net) to realize the "Select Features by Polygon" on ArcMap (shown below). I am using VS 2012 (vb.net) and ArcObjects 10.2.2, thanks. ![enter image description here](http://i.stack.imgur.com/IECzx.jpg)


title: Census Shapefiles Fit to Land

text: I'm looking for an oceans or coastlines shapefile that's a more detailed than what Natural Earth has to offer (1:10m). I'm working with 2000 Census county shapefiles, but it doesn't cut out lakes or oceans. So the county boundaries are a bit blobby. If I clip it to Natural Earth's oceans, a sliver of New York County gets left inside King's County. (I'm a perfectionist using QGIS and Tilemill). ![County shapefiles example](http://i.stack.imgur.com/sAp7l.png)


title: How to set waste water gravitational flow direction using DEM automatically?

text: I have a Geometric Network of Waste Water and I need to set line digitize direction according to the gravitational flow direction using a DEM. I'm trying to avoid flipping manually 1000 Km of pipes... I'm well acquainted with ModelBuilder but a beginner with python so any hint will help!


title: How to get a centroid of a closed line in MYSQL

text: I have a table that contains only closed lines (lines where endpoint is the same as startpoint). I would need to find a centroid for those lines. How would I be able to achieve it? It does not need, to be inner centroid, it can also be just an average of all the line nodes. Thank you.


title: How to address Geodatabase must be enabled as geodatabase admin user error?

text: I am trying to use "enable Enterprise Geodatabase" but I am getting the following error:               "Geodatabase must be enabled as geodatabase admin user.Failed to execute "      I am using postgresql 9.1.3 database. I have created the database named geo using this command: `createdb geo` and sed user using: `CREATE USER sde WITH SUPERUSER LOGIN PASSWORD 'sdepassword';` and then I connected to database as sde user. How can I solve this problem?


title: projections Basemap and pyproj

text: I am struggling to understand the projections using **pyproj**. for now my question is to understand the results of projection operations I have following coordinates that I project on x,y               lon  = [3.383789, 5.822754]     lat = [48.920575, 53.72185]          # with Basemap      from mpl_toolkits.basemap import Basemap     M = Basemap(projection='merc',ellps = 'WGS84')     q1, q2 = M(lon, lat)     print zip([x for x in q1], [x for x in q2])     # this gives me      #[(20414190.011221122, 65799915.85232639), (20685694.35308374, 66653928.94762396)]          # with pyproj     from pyproj import Proj     p = Proj(proj='merc', ellps='WGS84',errcheck = True)     print p(3.383789, 48.920575), p(5.822754, 53.72185)     #this gives me      #(376681.6684318804, 6229168.979819128) (648186.0102944968, 7083182.075116195)      why are the results different while I use the same projections parameters as a newbie in geospatial data processing I apologize in advance for a question which may be trivial P.S: I edited the questions following comments


title: layer disappear on Openlayer after specific zoom in

text: i have added some layer(s) and baselayer using openlayer v 2.11, layer does not appear when it first opened the page. But After zoom out it does appear. Why is behaving in this way, any params for always visible?               var map = new OpenLayers.Map("map");                      var wms = new OpenLayers.Layer.Google("Google Streets", {                     numZoomLevels: 12                 });                           map.addLayer(wms);                 var wms2 = new OpenLayers.Layer.WMS("Basic1", "http://localhost:8080/geoserver/wms",                {                    layers: 'SomeLayers',                    transparent: true,                    numZoomLevels: 12                }, {                    isBaseLayer: false                });                 map.addLayer(wms2);                 var ls = new OpenLayers.Control.LayerSwitcher();                 map.addControl(ls);                 map.addControl(new OpenLayers.Control.ZoomBox());                 ls.maximizeControl();                      map.setCenter(lonlat, 6);


title: Layer relations via Python

text: Is it possible to, through Python, to create the relations that it's possible to create in Layers/Properties/Fields/Relations, and if so how? I first thought these were saved in the styles, but that was not the case.


title: Trying to extract arc information from a Suburb shapefile

text: I **downloaded some shapefiles** from the **Australian Bureau of Statistics** ( **ABS** ) Its basically showing a bunch of suburbs when I view the shapefile through a GIS tool. The GIS tool that I am using at present is the **GRASS 6.4.2** tool. So basically the shapefile has the below statistics               Number of points:       0               Number of areas:      15231           Number of lines:        0               Number of islands:    7003            Number of boundaries:   31362           Number of faces:      0               Number of centroids:    15220           Number of kernels:    0      I'm trying to extract from this shape file the arcs. That is the polylines of the boundaries. This is because I need to find for each arc its accompanying adjacent boundary. For example arc1 may have boundary with Cat id - 135 on its left and Cat ID = 125 on its right. Any clues as to how to do this would be appreciated.


title: Land Cover Feature Extraction from Satellite Imagery

text: I'm interested in a low-cost or open source solution for creating land cover GIS layers that utilize both spectral and textural extraction algorithms. I have used PCI Geomatica, ENVI, and Feature Analyst VLS in the past; however these solutions are a bit beyond my price range, any software recommendations?


title: CartoDB: custom sorting infowindow entries

text: I am trying to create a visualization at runtime like explained in the Cartodb documentation:   Then I add infowindows for each entry with this command:               cdb.vis.Vis.addInfowindow(map, layer.getSubLayer(0),        ['center_name','phone','age_range']);      But the entries are displayed in alphabetical order, not the order I want it to be. Is there any way to force it to keep my order?


title: Merge single band and multi band color images

text: I do have 2 color images (orthophotomaps). One is only 1-band with color map and second is with 3-bands. How can I merge them with any FOSS tools or ArcGIS? Thanks a million


title: Can't project ArcGIS Tiled Service for Esri-Leaflet

text: I'm following this example from the esri-leaflet docs. I have an ArcGIS Tiled MapService in `ESRI:102696`. So I'm trying to use Proj4Leaflet to project the service and use it in my application. I found some proj4defs here and here. They do differ slightly, but neither seems to work. This is the Map Service I am trying to use. I put together a JSBIN example with the code below.                 // The two proj4defs I found       //var crs = new L.Proj.CRS('ESRI:102696', '+proj=tmerc +lat_0=35.83333333333334 +lon_0=-90.5 +k=0.999933 +x_0=250000 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192  no_defs', {       var crs = new L.Proj.CRS('ESRI:102696', '+proj=tmerc +lat_0=35.83333333333334 +lon_0=-90.5 +k=0.9999333333333333 +x_0=250000 +y_0=0 +datum=NAD83 +units=us-ft +no_defs', {         transformation: new L.Transformation(1, 17627100, -1, 19798600),         resolutions: [           303.8194444444444,           217.01388888888889,           86.80555555555556,           65.10416666666666,           43.40277777777778,           20.833333333333332,           8.680555555555555,           4.166666666666666,           2.083333333333333,           1.0416666666666665,           0.5208333333333333         ]       });       // pass crs into map constructor       var map = new L.Map('map', {         crs: crs       }).setView([38.5432816667, -90.2688116667], 5);  L.esri.tiledMapLayer('http://maps.stlouisco.com/ArcGIS/rest/services/Maps/PublicBase/MapServer', {         maxZoom: 10,         minZoom: 0,         continuousWorld: true       }).addTo(map);      I'm not sure if the transformation is off or what. Any help would be appreciated.


title: Need shapefile for worldwide bathymetry

text: I need a (free, no cost) shapefile for bathymetry of entire world, specifically from depths of zero to 200 metres, ideally in 50 or 100 m contour intervals. Has to be a shapefile. Help please!


title: flex/AIR based SLD generator?

text: Styling layers in GeoServer can be tedious to say the least. While Styler can relieve some of that tedium, it doesn't yet work in chrome. It can wreck existing styles if not carefully used and still involves a fair amount of manual labor. I was wondering if there were any flex/AIR SLD generators floating around?


title: Does ArcGIS for Desktop support SpatiaLite?

text: SpatiaLite is playing an increasing role in my organization. As it stood in 2012, using SpatiaLite alongside ArcGIS Desktop required an amount of importing and exporting of data that was not always trivial. Outside of products such as the Data Interoperability Extension or FME, what other options are there for SpatiaLite support in ArcGIS for Desktop?


title: How to keep layer style/symbology when saving as KML?

text: When saving as a kml, is it possible to keep the style of the QGIS layer? I have tried a couple plugins that have been mentioned on here: qgis2kml ('format not yet supported') and qgis2google2 (crashed QGIS upon activating). Is there anything else?


title: Turn-by-turn direction using PgRouting

text: How can I generate turn-by-turn directions using PgRouting? At the moment, I have a web-page where a user can enter and start/end location and a route is generated on a map (using PgRouting and OpenStreeMap). I know how to list all of the street names, but I don't know what the best approach for stating which direction the user should take at intersections. I have read through previous questions on this site: Does pgRouting support turn-by-turn routing? and driving directions in pgRouting but both of these questions are over two years old, and several of the answers link to sites that no longer exist.


title: Does pgRouting support turn-by-turn routing?

text: I have route list as result of pgRouting with street names. How to get messages "Turn left", "Go along" etc ... Can pgRouting do this or I need some other source and which?


title: Best way to join and handle multiple attributes related to the same feature

text: this is not going to be a question, I'm just trying to find some advice. Suppose I have a point shape file (for example just 2 point, actually I have many more points) representing water control stations. The attribute table looks like:               ID    X        Y          Cl       Na     A    3480240   5887204   125.70   60.30     B    3480450   5887280   230.30   120.80      Now, how can I join another table that contains repeated measures of `Cl` and `Na` for each `ID`? This could be the example of the other table (measures in different years) I want to join:               ID    Cl       Na      Year     A    100.30   210.20   2010     A    120.40   200.30   2011     A    130.40   198.20   2012     B    101.00   199.20   2009     B    130.40   210.40   2010     B    140.60   210.00   2011     B    166.90   240.20   2012      I have tried both the QGIS standard join option and the mmqgis plugin (join from csv file) but the results are not as expected. QGIS join option adds some columns to each `ID`, but many fields have "NULL" attributes while mmqgis increased the number of rows for each `ID` and each field contains the same value (e.g. `Cl` values are always 120.20 for each year). So I would just like to know if anybody has worked with this kind of data in order to understand the best way to handle multi-attribute field related to one single feature. Thanks


title: Active Directory AND integrated user store on ArcGIS Server 10.2

text: We have one License for ArcGIS Server 10.2. We need to publish Services for internal use with Active Directory authentication and for external use (internet) with username/password authentication. What would be the best way to achieve that? Is it possible to create two sites with one AGS License on one server so that I can use different user-stores?


title: Filter Vector Layer Features

text: I'm using OpenLayers - 2.11,GeoServer 2.13, ExtJS 3.4, and Java for developing GIS based web applications. I have added some point features in one vector layer. It has id, name, address, etc. Can I filter some points with 'name' from that vector layer and only show matched features in vector layer? I have to filter some points and later I have to show all hidden data. eg: ![enter image description here](http://i.stack.imgur.com/vcj6P.png)


title: Are Date columns in feature classes stored as python datetime objects?

text: I have a Date column called readings_measurement_tstamp stored in my feature class. I would like to add a String column to this feature class and Field Calculate it to be be just the %H:%M:%S component of readings_measurement_tstamp. After adding the new column I calculated it according to the following expression: `datetime.datetime.strftime( !readings_measurement_tstamp! ,'%H:%M:%S')` I was under the impression that Date columns can be referenced like python datetime objects. Therefore, I was very confused when running the field calculator with this expression gave the following error: `Descriptor 'strftime' requires a 'datetime.date' object but received an 'str'` The error implies to me that readings_measurement_tstamp is being stored as text, not as a Date. Any ideas what I need to do in order to be able to access readings_measurement_tstamp as a datetime object? Thanks!


title: Can I "undo" buffering of a point layer in QGIS?

text: I have point layer, where each point is one tree, about 1300 trees in total. I made buffers around each tree, the buffer being the height of each tree in QGIS 2.0.1. My next step was to use Spatial Query to see if the buffers (the height of each tree) reached a road. The resault was that 60 buffers/trees reached the road. I now want to show the location of each tree on a map, and I believe that the map would easier to read if it contained points instead of buffers. Is there a way to go back and make points again from the layer with the selected buffers in QGIS?


title: How to iterate RGB raster layers in ModelBuilder?

text: I have an RGB raster that I'm calculating color-based indices from. When I do the calculations manually, I can select a certain band of the RGB raster to calculate. I need to do the same operations on a folder of rasters, so I use ModelBuilder to iterate the folder, but I don't can't select the individual R, G, B layers from each raster and do calculations on them I've found this question, and have defined the Raster Calculator expression as such, but I receive the following error: ![enter image description here](http://i.stack.imgur.com/GFeLi.png)


title: Is it possible to use Arcpy in ArcGIS 10 to zoom to a selected feature?

text: The answer to this may be a simple yes or no, but I'm hoping that someone is aware of a way to zoom to the extent of a selected feature using the ArcPy modules in ArcGIS 10.


title: What GIS application can open and use .osm files?

text: A quick question that I don't see an answer to elsewhere: if I download a .osm dataset from Planet.osm (for example) which GIS applications could consume/display this data without translation? For example, this page says ArcGIS 10 "can use OpenStreetMap as a basemap" - but I'm not sure if that means it will open .osm files or read the data from a service? So, are there any GIS apps (proprietary or open source) which will read .osm files? Or perhaps taking a step back: is OSM data meant to be consumed like this, or is OSM meant to be used mostly for dynamic display of mapping on web pages?


title: Scriptable alternative to ArcScan (particularly vectorising lines within tolerances)

text: Sadly the answer to this question suggests that it is not possible to script ArcScan from within a python script. Therefore, does anyone know of any alternatives to ArcScan (either within ArcGIS, or within another GIS, or anything really...) that are scriptable? The functionality that I am in need of is the ability to convert from raster to vector by 'intelligently' vectorising lines. For example, ArcScan allows you to set a distance tolerance so that any lines that are vectorised and are less than that distance apart will be joined together. It also does various other sensible things with raster lines (for example, if a line is 3 pixels wide it vectorises it to a single line placed in the middle of the raster line). Any ideas?


title: Setting color of a feature "street" with a FeatureRenderer

text: I have a layer with streets and I'm trying to change the color of the streets. I know I asked a similar question, How do I change the colors of my layer file?, and the answer worked for points. This is my attempt at getting it to work with a street layer.               IFeatureRenderer featureRenderer = geolayer.Renderer;     int streetFeatureID = geolayer.FeatureClass.FindField("STREET");     IFeature feature = geolayer.FeatureClass.GetFeature(streetFeatureID);     ISymbol symbol = featureRenderer.get_SymbolByFeature(feature);     ILineSymbol linesymbol = symbol as ILineSymbol;     linesymbol.Color = SymbolCreator.getColor(255, 0, 0); /*my method which creates an RgbColor*/      This seems like it should work - there are no nulls or errors - but I'm missing a way to set the symbol. I'd like something like this:               featureRenderer.set_SymbolByFeature(feature, linesymbol);      Since I can't find that method I think I must be doing something completely wrong. Thanks in advance. Using ArcEngine 10 C# VS2010


title: GDB file and USA Base Map.mxd not in the same directory ArcObjects

text: I am trying to open a address locator in ArcObjects by following this tutorial: http://help.arcgis.com/en/sdk/10.0/arcobjects_net/conceptualhelp/index.html#/Creating_an_address_locator/0001000004pm000000/ The code sample has this line: String workspacePath = @"C:\UnitedStates.gdb"; I talked to an analyst and they told me to find the sample in this folder "C:\Program Files\ArcGIS\DeveloperKit10.0\Samples\data\Usa" And that has a file called USA.gdb However my "USA Base Map.mxd" file is located here, Z:\ESRI\DataMaps10\usa It seems like ArcEngine should use the same map / data to do both geocoding and projections... right? Am I doing something obviously wrong?


title: How can I extrude a polygon "below" a base height?

text: I am creating some animations and short movies in ArcScene showing how a number of disparate datasets relate to one another. In one of these movies my "big finish" is to extrude a number of polygons "negatively". In other words, I have a number of data layers stacked on one another in space (in some cases naturally as a function of z-coordinates, in some cases as a function of using a layer offset, and both where I need an exaggerated space between layers) and I want to show relationships between the different layers by using extrude to drop polygons on "upper level" layers to lower layers. I can't get polygons to extrude below their base heights... It seems weird to me as it's possible to drop layers below their base heights using negative numbers as layer offsets. When I put negative numbers in the Extrusion Value Or Expression box in the Extrude tab in Layer Properties nothing happens. Am I missing something? Is there a specific expression I need to use? Also, if negative extrusions are possible, I'd like to be able to cleanly clip the bottom of each extruded polygon to it's final intersecting surface. In other words, when the intersecting surface has texture, my extruded polygons (when I've tried extrusion "going up") the entire extent of my tops either don't sit flush with the surface or go through to the other side. I envision a "3D Extract by Mask" kind of thing. Any ideas on either of these fronts would be greatly appreciated!! I've seen a few non-ESRI tools that may do this stuff, but I'm limited to ArcGIS (in that I'd prefer not to have to learn a whole new software to deal with my other layers. Jesse


title: Help choosing a suitable routing engine

text: I am building a route planning system, but i have still to decide what underlying routing engine i will use. So far i have found pgrouting and neo4j. I have my route network in a postgresql/postgis database (imported from a shapefile). I have made queries to extract nodes (endpoints of ways where you have to make a decision what direction to go or dead ends) and to extract edges (often made up by several consecutive ways). All my edges are bidirectional. My main goal is to calculate routes on this network using an A-star algorithm where distance =cost. My feeling tells me that a graph database like neo4j is the way to go (as it seems to be made for just this purpose), but they don't seem to support A-star by default and also there is no real sense of geometry. It seems better suited for social networks instead of maps.   * Would pgrouting fulfill my needs?   * Is it fast enough for on the fly queries (+-2000 nodes, +-4000 edges)? Normally this would be a few ms for A-star, but i am unsure about this implemenatation in sql.   * Does pgrouting A-star give me a list of nodes and edges?   * In most examples i see about pgrouting i notice there is usually a list of commands after calculation of the route (like "At X turn left, etc"). Does pgrouting produce this or is this from another system? Hopefully someone can give me some information about what system to choose. Neo4j, pgrouting, or some other system.


title: QGIS Init script: How to get active layer and GID of edited object?

text: I would like to ask, how to get the GID (from Postgresql+Postgis) and active layer of edited object in QGIS Init script (Python)? Code Layer=qgis.utils.iface.activeLayer() works in QGIS console, but it doesn't in Init script, which runs when object is edited (in form made by PyQT designer). Thanks for reply! Seb * * * To get layer name we can use layerid.name() from def open_form (dialog,layerid,featureid): To get GID of edited obcjet we cen read it from dialog form GID = dialog.findChild(QLineEdit,"gid")) One more thing: How to check if active layer is enable to edit?


title: How can I get the values returned by an OpenLayers bounds object in Longitude / Latitude format?

text: I've created a map and added a layer using the Open Layers library:               map = new OpenLayers.Map('map');     var gphy = new OpenLayers.Layer.Google("Google Physical", {type: G_PHYSICAL_MAP});     map.addLayer(gphy);      The map loads and I can center my map:               map.setCenter(new OpenLayers.LonLat(-114.922485, 56.400224), 4);      I've then created a box handler, but the bounds returned aren't in Longitude / Latitude format. When drawing a box around the center of my map, my bound object returns the following:               left-bottom=(278,245) right-top=(304,211)      How can I convert these values into Longitude / Latitude? The resolution thanks to unicoletti:               map.getLonLatFromPixel(new OpenLayers.Pixel(bounds.left, bounds.bottom))


title: Determining ESRI File Geodatabase release version?

text: Is there any method in Python, ArcObjects, or the File GDB API for determining the release version (9.3, 10.0, etc.) of a file geodatabase? There is a need to find all of our shared 10.0 version file geodatabases and convert them back to 9.3 since not everyone as at 10.0 yet. I know that you can create 9.3 version file geodatabases in 10.0 so that part should be easy enough.


title: Is it possible to compile a grass addon into a commandline app?

text: As the title says, I'm wondering if it's possible to take the source of a grass add-on and compile it into a commandline app. Or, alternatively, use the add-on in a python program. Basically, I don't wish to use grass but I do want the power of an add-on for use in a separate program. thanks!


title: OpenLayers With StreetView Overlay

text: Im trying to add StreetView functionality to my OpenLayers web app, and would like to add the overlay showing StreetView availability to the map itself. I have seen an example of it being done successfully here, but i cannot get the code to work in my own solution. Specifically :               var layer = new OpenLayers.Layer.Google("Google Street", {sphericalMercator: true});          layer.mapObject.addOverlay(new GStreetviewOverlay());      results in the error:               layer.mapObject is null      Has anyone got any idea what i might be doing wrong, and if so, how can i fix this? Thanks in advance! C


title: How to always enable an ArcGIS toolbar?

text: On ArcMap 9.3.1, is there a way to ensure that a toolbar is always enabled? In my case, I want the Standard toolbar (the one with new/open/save/print icons) to always be visible. I'd like the toolbar to be enabled on ArcGIS startup, and to prevent the user from disabling it. Is there a callback that fires when a toolbar is disabled where I can write code to prevent it from disabling, as well as API calls I can make during ArcMap startup to always enable the toolbar? Thanks for any help.


title: Converting ESRI Web Client JSON to Standard Feature Class With Python Geoprocessor

text: Is there a function exposed by ESRI's python geoprocessor that converts the JSON format used in the various web clients to a valid feature class. Clearly, the function exists because something is converting these JSON strings into a format usable by the geoprocessor. I can't find it anywhere.


title: Upgrade guidelines for ArcGIS 9 to 10?

text: We have a copy of ArcGIS 10 that we are preparing to use to upgrade our existing ArcGIS 9.3 installation. Should I uninstall ArcGIS 9.3 first? Or do I install ArcGIS 10 over it? Please let me know if there any best practises or guidelines I should be aware of, tahnks!


title: Load Geoexplorer into Tomcat

text: I have installed Tomcat on a Windows server and have successfully deployed Geoserver. I am trying to deploy geoexplorer by building the war file from github and loading it into tomcat to deploy. It loads, but when I go to localhost:8080/geoexplorer, I get a 404 error. How do I load it into tomcat if not that way? EDIT: I did load Geoexplorer into the webapps folder of Tomcat and restarted Tomcat. Here's what is in the log about Geoexplorer:               Jul 07, 2014 8:34:34 AM org.apache.catalina.startup.HostConfig deployWAR     INFO: Deploying web application archive C:\Program Files\Apache Software Foundation\Tomcat 7.0\webapps\geoexplorer.war     Jul 07, 2014 8:34:40 AM org.apache.catalina.util.SessionIdGenerator createSecureRandom     INFO: Creation of SecureRandom instance for session ID generation using [SHA1PRNG] took [240] milliseconds.     Jul 07, 2014 8:34:40 AM org.apache.catalina.startup.HostConfig deployWAR     INFO: Deployment of web application archive C:\Program Files\Apache Software Foundation\Tomcat 7.0\webapps\geoexplorer.war has finished in 5,844 ms


title: Accuracy of planimetric geometry vs. ellipsoidal model

text: I'm examining the two geospatial options for IBM's DB2 database. In particular I am interested in the accuracy between these two options and whether it's worth it (the more accurate representation isn't free). The thing I am most concerned about is someone doing a query (e.g. looking for a point or other geometry within a polygon) and the database incorrectly returning an empty resultset. My use case is a catalogue indexed by several means (temporal, geospatial and XML) and entries can be anywhere on the globe. IBM states (PDF link) that: > Spatial Extender treats the Earth as a flat map. It uses planimetric (flat- > plane) geometry, which means that it approximates the round surface of the > Earth by projecting it onto a flat plane. This projection causes > distortions, which can vary across the extent of the data, but the > distortions generally increase toward the edges of the projected > region…Spatial Extender is best used for local and regional data sets that > are well represented in projected coordinates, and for applications where > location accuracy is not important. The other option is: > Geodetic Extender treats the Earth as a globe. It uses a latitude and > longitude coordinate system on an ellipsoidal Earth model. Geometric > operations are precise, regardless of location…Geodetic Extender is best > used for global data sets and applications that cover large areas on the > Earth, where a single map projection cannot provide the accuracy required by > the application. If I was using something like WGS84 in the Spatial Extender geospatial implementation, can anyone give me a ballpark figure of the worst-case accuracy errors I might encounter when compared to using the Geodetic Extender? Is there a formula or process for working it out for anywhere on the globe for a given map projection?


title: binary map file editing

text: in Mapsforge binary .map file exist a zoom table.defining of this table is: > A table indicating the number of POIs and ways in this tile for the > different zoom levels covered by the enclosing sub-file. Let Z be the number > of zoom levels supported by the enclosing sub-file (e.g. 6 for a sub-file > that covers levels 12-17). Then the table has Z rows and 2 columns (first > column: POIs, second column: ways). Each cell in the table represents the > number of POIs or ways on the specific zoom level. The table is written row- > wise and values are encoded as VBE-U INT. I want add POIs & Ways to the table.


title: Using ArcPy to fetch Information from a .GDB

text: I am new to GIS and I would appreciate if you could guide me a little. I have a .gdb folder with data and I need to process that data to find educational institutions which are within say X Km of the Roads.....And within Y Km of any other educational institution.... I need to generate a shape file with the resulting data. Could you please tell me how shall i go about this? Do i need to install ArcGIS to do this? Or is it possible to do this in ArcPy alone? I am pretty ok with Python and I would appreciate if someone could guide me with a Python Snippet.....I am not looking for full code but just an idea to get me started. Appreciate your help on this one. Regards, Guddu


title: How to route on a Shapefile?

text: We use shp.file in our VB6 to run tracking system. Now we want to convert it into the navigation system, which requires doing the route analysis on shp.file. The route is drawn in shp.file Usually, we see the route analysis is on the grid when programmed it using Visual Studio 6. How to make it to run on shp.file? Thanks.


title: Extract by mask run time error in VBA

text: I have the following script language code shown below within a main code just to run the tool of spatial analyst "Extract by mask". Unfortunately, when I run it I get a Run Time Error. I checked file names etc. but I cannot find the problem. Could anybody help me please? Thnaks.               Dim GP As Object     Set GP = CreateObject("esriGeoprocessing.GpDispatch.1")          GP.AddToolbox "C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst  Tools.tbx"               InRaster = "C:\LACONISS\LandSpaCES\valuationmapr"     InMask = "C:\LACONISS\LandSpaCES\Parcel165"     OutRaster = "C:\LACONISS\LandSpaCES\ParcelR165"          GP.CheckOutExtension ("Spatial")          GP.ExtractByMask_sa InRaster, InMask, OutRaster


title: Is there a way to see the coded value domain list of possible values from arcmap 9.3.1?

text: Is there a way to see the possible values of a specific coded value domain in arcmap 9.3.1 (I have no access to the SDE which holds all the data in arcCatalog)?


title: Openlayers FramedCloud - Add Offset

text: I'm working with the example here: http://openlayers.org/dev/examples/light- basic.html When you hover over a circle, the popup appears. However, the behaviour is 'flakey' because once you hover over the tip of the popup, it disappears. Obviously this is happening because we've moused out of the circle itself. While functional this is user unfriendly. I've been trying to add an offset so that the tip is never on the circle itself. However, since it's not known before you hover over the circle which quadrant (top right, top left, bottom right, bottom left) the popup will be in, it is impossible to calculate the offset correctly. How would I check the position of the tooltip, then apply the offset? Fixing the popup to a specific position is not acceptable as a solution, because on smaller screens and with restricted extents some popups will appear almost entirely off screen.


title: Create a Polygon

text: I'm working with ArcGis to create a habitat mapping, so I have digitized the most part of my area with different polygons. I would like to create another one but it's really difficult because it is like a background polygon. I have to snap this polygon with all the others that I have just created. Do you know if is there a good way to do this? I wouldn't do this polygon manually.   Thanks!!!


title: Where can I find a tutorial for installing an addon to GRASS GIS that a six-year-old could understand?

text: I want to use the r.stream.distance addon in GRASS GIS 6.4 on ubuntu 10.04 (64 bit). I've installed grass from the repositories and that works fine but I can't seem to figure out how to install this addon (or any other for that matter). I've tried to run g.extension but I get a "command not found" error. I've tried to install the addon from svn but the instructions are less than clear to me. I've read through the documentation on the grass wiki and it seems to assume a lot of knowledge that I seem to lack. So, does anyone know of a tutorial where I can walk through the steps (baby steps) of installing an addon in grass? Thanks.


title: Getting XY screen coords from ArcScene 3D point using ArcObjects?

text: I already know how to convert screen point (X,Y in pixels) to 3D point:               ISceneGraph pSceneGraph = m_pSceneHookHelper.SceneGraph;     IPoint pPoint;     pSceneGraph.Locate(pSceneGraph.ActiveViewer, X, Y, esriScenePickMode.esriScenePickFastLocation, true, out pPoint, out pOwner, out pObject);      But how may I get X,Y screen point from 3D point? I just use snapping to change 3D point coordinates, so I need to know new screen coordinates to use Locate function (because I need to find feature object).


title: QGIS: (Shapefile) Convert rings to fake-rings?

text: I hope I can make my question I posted yesterday clearer... sorry! Basically, I want to obtain the opposite of this: How to convert fake donut holes in polygons to real ones? Such polygons will work perfect for me, as the program I use to import shapefiles for my use cannot recognize rings properly. Instead all is loaded as polygons and rings do not appear as holes...


Is ESRI metadata standard compatible with GeoNetwork metadata standard? if yes how do i update the Geonetwork metadata using Arc-Catalogue 10 Metadata xml?


title: GDAL (Python) Aggregate Raster Into Lower Resolution

text: I have a (continent-scale) raster with 30m resolution in an Albers equal area projection. Each 30m cell has a value 1-34. I would like to aggregate information from this raster into a 1km grid such that each 1km grid cell contains a floating-point value which is the number of 30m cells with a value of 2 divided by the total number of 30m cells which fit into the 1km cell. Is there a way of doing this with GDAL? Or an efficient way to hook GDAL and Python together to perform the calculation (given the large size of the raster, not all of it fits into RAM at once)?


title: Building Custom Programs with ArcObjects

text: I found another thread at ArcObjects Resources that asked for general ArcObjects resources but most of the resources in that thread are obviously about building add-ins for ArcMap. I need help getting started with building a completely custom application. In other words, if an ESRI logo or ArcGIS button shows up on my program, something's gone really wrong. Where can I look for resources on this topic? So far, http://www.amazon.com/dp/1118442547/?tag=stackoverfl08-20#reader_1118442547 looks possibly relevant but it's hard to tell because I can't see if the book talks about custom development from the first few pages that Amazon.com provides. To be super specific, I'm trying to write a program that takes coordinates and a layer that's in an Enterprise Geodatabase (aka ArcSDE) and figures out whether or not each point is inside the layer. I've already found a helpful function in the ArcObjects SDK.


title: Call ArcPy script via CGI on IIS

text: I have an ArcPy script that does a whole bunch of processes with ArcGIS 10.1 to write the output data to a SDE geodatabase. What I want to achieve is to enable an online user to call the script and then do the job on the backend automatically. I figured out CGI and made a mapping to py on IIS. However, everything I called the script from the web. It threw out the following error: HTTP Error 502.2 - Bad Gateway The specified CGI application misbehaved by not returning a complete set of HTTP headers. The headers it did return are "Creating catid files ESRI\Server10.1\Configuration\CATID\arcgis-response- streamer.catid.ecfg Create JARID File... outfile=ESRI\Server10.1\Configuration\JARID\arcgis-response- streamer.jarid.ecfg Creating catid files ESRI\Server10.1\Configuration\CATID\arcgis_indexer.catid.ecfg Create JARID File... outfile=ESRI\Server10.1\Configuration\JARID\arcgis_indexer.jarid.ecfg Content-Type: text/plain;charset=utf-8 Hello World! starts Hello World! ends ". I am just wondering if ArcPy script can be exposed to CGI like a normal python script? Thanks Cheers, Alex


title: Is it possible to get Elevation from a Satellite Image?

text: I would like to get elevation from a satellite image. I would be needing resolution less than or equal to 5 m. Is this possible? If yes, how can I procure these images and where?


title: Multiple entries in attribute table for multiple buffers using ArcGIS for Desktop?

text: I am using ArcGIS 10.2.2 to draw 500ft buffers around a group of 133 bars. The numerous buffer tools all produce a single entry in the attribute table for the collected group of buffers instead of one entry for each buffer when I use the dissolve setting. Is there a way to produce a shapefile of dissolved buffers whose attribute table contains 133 rows, that it, one for each buffer, instead of the default that produces one row for all? ![Print screen jpeg of numerous buffers and open buffer attribute table](http://i.stack.imgur.com/lTkcC.jpg)


title: Qt Designer and QGIS -

text: I am just trying understand how qgis and designer work together. From Nathan W's excellent posts I know how to create a user form and how to hook it to a qgis layer and also how to add validation code. My questions here are pretty naive but the stuff that everyone knows is precisely the stuff that one can't find in print:   1. Is the copy of Qt Designer that comes with QGIS integrated wtih the QGIS GUI in anyway (like the forms designer within ms-access is the the rest of ms-access)?   2. Can I start Qt Designer from within QGIS or do I have to have to fish around in the QGIS installation with Windows Explorer and find the Qt Designer exe to double click on?   3. Am I better off installing a separate version of Qt Designer/Qt Creator or pyQt than using the integrated Qt Designer?   4. Does the pre-packaged windows installer for new users come with Qt Designer configured in a way that can be used for form design?   5. If I am looking for a way to create forms with python validation code, is PyQt a good solution?   6. Many data collection tasks are non-spatial, for example adding in a price list for road signs. The GIS can record which signs need replacing but the price list is needed to provide a cost estimate. Is there anyway of triggering the display of forms for non-spatial data tables?  That is a lot to ask I am sorry. Thanks.


title: Data providers equivalent to USDA/NRCS/FSA Common Land Unit (CLU)?

text: For one of my agricultural projects, I really need reasonably accurate field boundaries for fields in the US, specifically Iowa right now but later, other similar states. The USDA/NRSC/FSA Common Land Unit (CLU) data is exactly what I need, but apparently since 2008 it cannot be accessed by the public. Does anyone know of a similar source, commercial or otherwise, of this type of data? I don't need much for attributes, just field boundaries that will let me aggregate analytic data within a field, accurately.


title: Adding an ILayer to IMap in Arcgis Server 10

text: I wanted to add ILayer to IMap, i think it has been added but problem is that it does not display it on the website. I am running arcgis server 10, using visual web developer 2010 express edition for web mapping application, vb.net (web adf). Here is the code:                Public Sub ServerAction(ByVal args As ToolEventArgs) _         Implements ESRI.ArcGIS.ADF.Web.UI.WebControls.Tools.IMapServerToolAction.ServerAction              If TypeOf args.Control Is ESRI.ArcGIS.ADF.Web.UI.WebControls.Map Then                  Dim map1 As ESRI.ArcGIS.ADF.Web.UI.WebControls.Map             map1 = CType(args.Control, ESRI.ArcGIS.ADF.Web.UI.WebControls.Map)                       Try                 Dim mapFunc As ESRI.ArcGIS.ADF.Web.DataSources.IMapFunctionality = _                     CType(map1.GetFunctionality(0), ESRI.ArcGIS.ADF.Web.DataSources.ArcGISServer.MapFunctionality)                 Dim mrl As ESRI.ArcGIS.ADF.Web.DataSources.ArcGISServer.MapResourceLocal = mapFunc.Resource                 Dim sc1 As ESRI.ArcGIS.Server.IServerContext = mrl.ServerContextInfo.ServerContext                 Dim mapS As ESRI.ArcGIS.Carto.IMapServer = mrl.MapServer                 Dim mapSObj As IMapServerObjects = mapS                 Dim pMap As ESRI.ArcGIS.Carto.IMap = mapSObj.Map(mapS.DefaultMapName)                      MsgBox(pMap.Name & ": " & pMap.LayerCount)                      Dim f As ILayerFile = New LayerFile()                 f.Open("C:\data\a.lyr")                 Dim pLayer As ILayer = f.Layer                 pMap.AddLayer(pLayer)                      MsgBox(pMap.Name & ": " & pMap.LayerCount)                      mrl.RefreshServerObjects()             Finally                 map1.Refresh()             End Try              End If     End Sub      Please can you explain or show why the added layer is not being displayed? thanks!


title: Are there any Objective-C frameworks for shapefile drawing?

text: Hello can someone point me to library for this? Which can read .shp files and render them? Thanx. Mirek


title: Trying to place a Point geometry on OSM map, unsuccessfully

text:
I have a piece of code that takes a map layer and places a vector layer on top of that. Currently, I'm just trying to place a few points on top of the map, unsuccessfully. I managed to place 1 point but it would only stay in the middle of the map. I took all the html code and simplified it. Here is the code for the whole page, it's standalone, so you can just copy paste into a new html page and see the problem.

                                            To reiterate, how do you place the point at some location on the map? (I have lon/lat) Thanks! **Update** Just to make it easier for readers in the future. The answer which is marked as the solution below only hints at the full solution. The 'full solution' is to modify the point creation code as such:               // Make Point     coordinates.push(new OpenLayers.Geometry.Point(33, 33).transform(new OpenLayers.Projection("EPSG:4326"), map                             .getProjectionObject()));


title: How to join data from a CSV to a shapefile using QGIS

text: I'm new to GIS and new to StackedExchange. I need help joining a **CSV** of Accessibility data to a LSOA shapefile in **QGIS**. The process explained here does not work for me and I don't know why. The **mmqgis plugin** and the **Joins** feature in QGIS both provide me with either an empty Shapefile or an empty CSV file. I'm trying to map statistics related to accessibility to hospitals from the Department for Transport. The data can be downloaded from here and is in **CSV** format and is georeferenced by Lower Super Output Area (LSOA), with one column of LSOA codes and the rest are the separate statistics (time taken to get to nearest hospital, for example). I also have a LSOA shapefile from here, which obviously has the **exact same list of LSOA codes**. How can I join the two datasets so that I can map the accessibility stats for each LSOA with a range of colours? I've tried using the 'Join' feature but I can't seem to get it to work. Any help would be much appreciated. * * * The first image shows the Accessibility data, note the LSOA data in the first column which matches the LSOA data in the Shapefile and the last column is what I want to map/Join. The second image shows the data associated with the LSOA shapefile, note the LSOA data matches that in the Accessibility data. As far as I know the data type in each file matches. What I want to do is map the data in, say, column 5 of the Accessibility data onto the LSOA shapefile. ![enter image description here](http://i.stack.imgur.com/QCJL7.png) ![enter image description here](http://i.stack.imgur.com/VUC4E.png)


title: How to Embed KML files (Google Earth) into a website without the google gadget?

text: I am trying to embed a KML file onto our webpage and the way it was previously done, using the KML embedding tool from Google no longer works, it seems that google has taken the gadget down or something. Is there another way of doing this that isn't overly complex? Any Suggestions and or advice is greatly appreciated


title: Need to convert points to polygons of a specific size, using the point as the center of the polygon

text: I have a large collection of aerial photos with no georeferencing other than the Lat/Long location of the camera. Georeferencing them invididually is not an option (too much time, too many photos, some are oblique, some vertical). I need to automatically generate some kind of approximate polygonal footprint of each photo. I was thinking there must be some way to auto-generate a rectangle from a point layer using the point as the center and using numbers in the attribute fields as the dimensions. Or, barring that, just generating a circle at each point using an attribute as the diameter. I'm using Arcview 9.3.1.


title: Symbolizing xyz data

text: I have hydrographic data in the form of an xyz file. I would like to display it in ArcMap. I need to symbolize the points based on their depth, shallow points one colour and deep points another. With a small subset of data I was able to import the symbols, symbolize based on elevation and turn on labels at an appropriate scale. ![xyz sample points](http://i.stack.imgur.com/li8hF.jpg) I then tried creating a multi-point file with the full dataset and set the symbology the same however I end up with unnatural bands. ![xyz multi point](http://i.stack.imgur.com/xFiyr.jpg) Is there a way to get the multi-point file to display the same symbology as the sample set? I want the user to be able to quickly zoom in and out of the data to be able to check the depth of any given point.


title: Problem with custom projection definition in Qgis

text: I want to define a custom crs valid for Mexico. The context is that older cartography in Mexico is based on the NAD27 datum. I wanted to test the following proj4 string `+proj=longlat +ellps=clrk66 +towgs84=-12,130,190,0,0,0,0 +no_defs` Nevertheless , when I test my string with 19N 99W (19,-99) I get the following values ![enter image description here](http://i.stack.imgur.com/AHTR2.png) testing the proj 4 definition included in Qgis gives the same results. These are wrong since the resulting north should be around 19 and the resulting east around -99. I wonder if anyone can tell me what I am doing wrong. Thanks in advance


title: How to run MapServer on an existing XAMPP installation?

text: How can I use MapServer on existing Apache installation? I don't want to install Apache server bundled with MapServer.


title: Turn Popups On/Off in ArcGIS Online Webmap via Javascript API Web Application

text:
I am writing my first ArcGIS Javascript API web map application and consuming a web map authored in ArcGIS Online. The ArcGIS Online web map apparently does most of the work for me in terms of symbology, editing templates, and even popups. For my application, I am trying to allow the user to edit features, but not all the time. I am able to toggle the editor on and off with the click of a button. The problem is that I have not been able to toggle the default popups on and off so that the user can edit attributes. Even when editing is enabled, the default non-editable popup appears when you click on a feature. The only way to avoid this is to use the `ignorePopups` option when calling `esri.arcgis.utils.createMap`. However, this means that when editing is disabled nothing happens when you click on a feature. Here is my code. This has all been adapted from several samples provided in the ArcGIS Javascript API documentation.

                             Edit

                         This pane could contain tools or additional content

           I realize that I could recreate a custom `InfoWindow` and `AttributeInspector`, but since this is already done automatically by the API, it would be nice if I didn't have to. I tried to get the `clickEventHandle` and `clickEventListener` from the response given by `mapDeffered.then()` and then `connect` and `disconnect` when the editor is disabled and enabled, but the default popups remained even after I called `dojo.disconnect(clickEventHandle)`. Any help would be greatly appreciated.


title: Oracle Spatial - How to calculate the point of intersection of 2 line geometries

text: How to calculate the point of intersection of 2 line geometries.


title: Where can I find the data that QGIS used to custom its menus as well as it's toolbars?

text: I have noted that QGIS uses some data to restore its menuitems as well as its toolbar items everytime when it is started.And I have being trying to hunt them out in my hard disks but sadly failed...The whole registry has been searched with the key word "OSgeo" in the "find" dialog,but nothing found to be related with what i want,neither under the whole directory of "OSGeo4W". File types of " _.ini","_.xml","*.cfg" have been used for search. That has made me sad,you all know. So,any kindness to help ?


title: QGIS 2.4 (and 2.0.1): How to uninstall/reinstall a plugin where uninstall is disabled?

text: Plugin Raster->TerrainAnalysis->slope which **had** initially worked, started showing what another questioner referred to as a "crazy histogram": qgis2 slope analysis from NED data gives crazy histogram Before finding his question and the answer, I endeavored to uninstall and reinstall the plugin. I found/still find (QGIS 2.0.1, now upgraded to QGIS 2.4, on Ubuntu 14.04), the plugin "uninstall" button grayed out/disabled. It's working, but I am worried about having a plugin I installed that does not offer an option to be uninstalled. Failing anything else--and following my penchant for "trying something"--I would be tempted to use the path information to rename the .so library and see if it disappeared, but better advice or insight would be most appreciated, thank you. ![uninstall grayed-out](http://i.stack.imgur.com/GvzOB.png)


title: Is there a limitation on the numbers of layers that can be displayed in QGIS?

text: I have been looking for an answer to this question for the last several weeks, so hope that someone here has encountered the same issue and has a viable solution to my problem. I am currently using the following configuration: • Windows 7 Pro • QGIS 2.2 (also tested 2.4 with same issue on a separate system) The issue I am seeing is that that once I reach 254 or 256 layers added to a project, as soon as I add any additional layers, all layers disappear from my map and I am unable to see any of my data. I work with large amounts of data. Is this a limitation in QGIS or is it a setting somewhere? I have also tested MapInfo, which I am thinking of buying if QGIS is unable to handle the load; MapInfo does not have this limitation. Any ideas? I can provide more information if necessary. Thanks


title: How to save a google map direction into a postgresql/postgis database?

text: I have taken a direction(route) from google map. ![enter image description here](http://i.stack.imgur.com/ddX9U.png) Can I save this resultant direction/path/route to my postgresql/postgis database? How?


title: shifting/seperating overlapping polygons into several smaller polygons within the same sapce (QGIS)

text: I'm using QGIS and have shape files from the City for parcels of land. There are a few apartments and condos in which the parcels have not been separated - i.e., instead of a dozen small rectangles right next to each other, each apartment's shape is the entire parcel of land with all dozen polygons stacked right on top of each other. Is there a way to programmatically take those layered polygons and split them into a dozen smaller rectangles (or whatever) which in totality would take up the same space? or any other suggestions that would better symbolically show the dozen individual apartments, rather than one big layered blob? Right now my only option appears to be editing each shape one by one by hand, which clearly can't be the best solution. Thanks,


title: How to deploy ArcGIS Engine (Standalone) with License Agreement?

text: I've completed a standalone app in ArcGIS Engine and in the process of deploying it, but I have a license agreement and supposed to add a licenseagreement.rtf file format in the properties menu while setting up in visual studio 2010. Am I to create this file myself or am I able to find it?


title: pyqgis: AttributeError: 'QgsRasterLayer' object has no attribute 'drawingStyle'

text: qgis gurus. I have a question about pyqgis and hopefully I can get some advise here. I am following http://www.qgis.org/en/docs/pyqgis_developer_cookbook/intro.html. In python console after a raster layer is loaded to a variable "rlayer"               >>>rlayer       <qgis.core.QgsRasterLayer object at 0x113faadf8>          >>>rlayer.drawingStyle()       Traceback (most recent call last):              File "<input>", line 1, in <module>              AttributeError: 'QgsRasterLayer' object has no attribute 'drawingStyle'          >>>rlayer.setDrawingStyle(QgsRasterLayer.SingleBandPseudoColor)       Traceback (most recent call last):              File "<input>", line 1, in <module>              AttributeError: type object 'QgsRasterLayer' has no attribute 'SingleBandPseudoColor'               >>>rlayer.setColorShadingAlgorithm(QgsRasterLayer.PseudoColorShader)        Traceback (most recent call last):              File "<input>", line 1, in <module>              AttributeError: 'QgsRasterLayer' object has no attribute 'setColorShadingAlgorithm'      what did I do wrong? thank you


title: Determine overlapping points using a spatial filter in ArcObject

text: I have a point layer with thousands of points. Using a spatial filter, I have to find out those points which are intersecting each other. Then the objectid's of the intersecting points should be displayed in a list box. I'm using ArcObjects and C# _and need a solution in this environment._ How should I approach this problem?


title: How to fix "import qgis.core ImportError: DLL load failed"

text: I am trying to develop a GIS application in QGIS using the Python language. Unfortunately, I encountered a problem from the outset is that of import and qgis.core qgis.gui. the Problem in the python console is as follows:               import qgis.core          Traceback (most recent call last):       File "<pyshell#0>", line 1, in <module>         import qgis.core     ImportError: DLL load failed: Le module spécifié est introuvable.      I am using:   * OS : Windows 7 , 32 bits   * QGIS : 1.8 ( Lisbona )    * Python : python 2.7


title: WMS getFeature control doesn't take 'viewparams' into account

text: I use Geoserver to get my WMS data. Everything works great for showing data in openlayers. Geoserver has the ability to create sql views and use parameters to return certain geometry. These parameters are called viewparams. This could be an example of a sql view:               select *     from myTable     where age = %ageParam%      When using viewparams you are able to pass the value for the 'ageParam' parameter. If a parameter is not set, a predefined default value will be used. More information about this topic can be found here. In openlayers however, "getFeature" control does not take the 'viewparams' parameter into account. This means it always returns features with the default parameter values. Is this a bug, a missing feature or something I just missed. **Edit1** Here is a screencast I made to visualise the problem. http://screencast.com/t/RgrQvaoyoqz As you can see, the layer is redrawn with new points but the selectFeature control still returns the old features. **Edit2**               <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.1.0" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">     <wfs:Query typeName="test:point_query" xmlns:test="http://test.boedy.p:8080/wms">     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">     <ogc:BBOX>     <ogc:PropertyName>geometry</ogc:PropertyName>     <gml:Envelope xmlns:gml="http://www.opengis.net/gml">     <gml:lowerCorner>490488.82657346 6800062.1058854</gml:lowerCorner>     <gml:upperCorner>491071.65891403 6800556.5579121</gml:upperCorner>     </gml:Envelope>     </ogc:BBOX>     </ogc:Filter>     </wfs:Query>     </wfs:GetFeature>      **edit3** This is the javascipt code for the layer and controller               test = new OpenLayers.Layer.WMS("test",         "http://test.boedy.p:8080/wms",          {             'layers': 'test:test_query',             viewparams: 'x:'+myX+';y:'+myY,             maxfeatures: 2000,             ratio: 1,             transparent: true,             isBaseLayer: false}         );          selectCtrl = new OpenLayers.Control.GetFeature(             {                 protocol: OpenLayers.Protocol.WFS.fromWMSLayer(test),                 clickout: true, toggle: false,                 multiple: false, hover: false,                 box: true,                 toggleKey: "shiftKey", // ctrl key removes from selection                 multipleKey: "shiftKey",                 eventListeners: {                     featureselected: onFeatureSelect,                     featureunselected: onFeatureUnselect                 }             }         );


title: ESRI File GeoDatabase: Unicode/Codepage problem - Convert UTF-8 to ANSI

text: I have a Codepage issue, and for now I am a bit stumpped My original source data contains multiple country data delivered in SHP (UTF-8). The data is then merged into a **File GeoDB** as the data is too large for the 2gb SHP/DBF limit. Previously I was able to use the *.CPG method with a smaller SHP dataset to correctly display the data that contains caracters that matches ANSI (ê, ë, ô, etc.). Now that the data is in a File GeoDB I am not able to use this method. Strangely enough, I changed the codepage of a smaller SHP-UTF8 data set to SHP-ANSI using a FME process, then imported the data into a File GeoDB and the data then displayed correctly in ArcMap. Now transforming all the source SHP-UTF8 data to SHP-ANSI is quite a timely process (many files that makes up a big dataset), so my question is if it is possible to **set or convert the CodePage of the data in a File GeoDB** so that the merged dataset would display correct? The UTF8 based data looks like this : KoÃ¶perasie Street Where it should look like (ANSI): Koöperasie Street Cheers. Pieter


title: projection tool as standalone script

text: How would I automate the define projection tool in arcmap so I can just enter my layer each time? I have it as a script in toolbox, but I am unsure how to set it so I can just input my layer as getparameterastext Here is the ESRI script - i want to replace the "forest.shp" with a layer of my own choosing each time                # set local variables      inDataset = "forest.shp"      coordinateSystem ="GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"          arcpy.DefineProjection_management(inDataset, coordinateSystem)


title: Bounding box for app engine?

text: I'm using python 2.7 and app engine. To do bounding box we used geomodel but it executes more then 5 seconds in table with 350 records. It's to slow.   Can you please suggest better solution ?   Is there online solution to store data and then execute bounding box search ?


title: How to configure mapserver correctly on cherokee?

text: I am attempting to configure mapserver using fastcgi on cherokee and here's what I have got at the moment:- ![enter image description here](http://i.stack.imgur.com/GM3Pg.png) ![enter image description here](http://i.stack.imgur.com/Cgg6R.png) ![enter image description here](http://i.stack.imgur.com/UbOpp.png) ![enter image description here](http://i.stack.imgur.com/CNHJz.png) Unfortunately, that doesn't seem to be enough as I somehow need to have mapserv running as a process. But I don't see such an option for mapserver in its docs. Any suggestion how I may configure mapserver correctly on cherokee?


title: Search Toolbar Register issues

text: When i have to use search tool, it will come warning dialogue box. How can I solve this issue. I have attached an example image to show the warning dialogue box. Thanks to advance![enter image description here](http://i.stack.imgur.com/OXEx9.png)


title: Trajectory Indexing: How could I store and index the large trajectory dataset generated by cars' GPS?

text: I have a large trajectory dataset generated by cars' GPS. Since I will do some data mining work on them, I need an efficient way to store and index them. I read some papers about **_Trajectory Indexing_** ( ** _Computing with Spatial Trajectories, by Yu Zheng and Xiaofeng Zhou, Microsoft research Asia_** ), 3D R-tree, TB-tree(Trajectory-Bundle tree), SETI(Scalable and Efficient Trajectory Index) were proposed to solve this problem. But it is quite difficult to implement them. Is there any implementation on this problem?   I am not quite familiar with POSTGIS. It has R-tree, but trajectories have a "time" dimension. Could I expand postGis to index trajectories? How?   Truely grateful for all anyone's help! 1,2000 taxis' trajectories for one month. About 50GB. latitude,longtitude,time. I need to search for trajectories of regions during anytime, and so on.


title: eTrex 10 UTC or GMT time zone

text: I have an eTrex 10. I'd like my GPS to remain in UTC (or GMT) time zone, but it does not offer this option as far as I can tell. I think "Western European" and "Ignore DST" is equivalent to GMT. Is this the best or only way to have data on this device recorded as close to UTC as possible?


title: Why aren't my basemap tiles showing up via ArcGIS Server

text: I have a leaflet js script referencing the ArcGIS server imagery api. The code worked fine up until about two days ago. Now when I click on the layer there's a gray background. I am using my local machine to deploy the map. My cloudmade maptiles are appearing without error, and there in no error in the console regarding truncated images or anything of the sort for the ArcGIS imagery api. I am using this api. http://www.arcgis.com/home/item.html?id=10df2279f9684e4a9f6a7f08febac2a9 ![enter image description here](http://i.stack.imgur.com/b1wQc.jpg) ![enter image description here](http://i.stack.imgur.com/qcsSa.jpg)               var aLayerTwo = L.tileLayer('http://services.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer',{                          attribution: 'Map data © ESRI, Created @ CSUN Geography',                     });