Insert image in attribute table(unable to save edits operation attempted on empty geometry)

I have been trying to associate/insert images in an attribute table. I have created a field called pictures and data type is Raster. ( I have used **this resource video** to do the task) Once I load the picture in the attribute table, every thing is ok. I can see the picture in the attribute table. But once I want to save the edits, I am getting this error: unable to save edits operation attempted on empty geometry Thank you in advance for your help. ArcGIS 10.1 (Advanced), Windows 8


why readOGR() is not reading projection information from .spj - r

0 down vote favorite I successfully import a polygon shapefile with readOGR, however the information in the .prj file is not being imported.               test class : SpatialPolygonsDataFrame nfeatures : 19407 extent : 35551.4, 1585917, 6318047, 9408727 (xmin, xmax, ymin, ymax) coord. ref. : NA          summary(test) Object of class SpatialPolygonsDataFrame Coordinates: min max x 35551.4 1585917 y 6318047.3 9408727 Is projected: NA proj4string : [NA]      the projection should be utm21n complex zone


Snapping doesn't work when many layers loaded!

I am using Arcmap 10, when I have only limited number of layers loaded, the snapping experience is perfect but when I have many layers loaded specially polygons, snapping process is painful. My cursor turns to "Busy" shape (the sand clock) and can't snap to the feature. I have no idea why it happens. Any solution?


List of actual GPS position

I have a question - where can I find a list of actual GPS position, for example - ships, airplanes etc. I need it for my project for College - I will create mobile/desktop/web app, which will show on map current positions of this elements, their last road and save info about that elements in database. Also, is there any website with informations about weather problems with lat/lng info? Thanks in advance, I've tried to find something about it, but without result..


ArcGIS Server REST get polygon centroid

I have a Site reference feature service data in polygons served up by ArcGIS Server. I need to get the centroid coordinate of a polygon by the Id but can’t find a method out of the box. I’m working on a GeoEvent Processor Service that will use this Site service to grab a centroid from a site Id. So far I believe my options are:   * Custom GP service, python – takes in polygon rings, returns point   * Server Object Extension – same as above   * Sql view, query layer served up as a feature service that uses STCentroid (The database uses a SQL GEOGRAPHY type) on the Site layer. Would query the layer by ID to return the geometry.  The first two options keep me on track with a service orientated architecture however it certainly seems overkill to be writing an SOE for one function. The final option could be a ‘goer’ but it means my logic is coupled to SQL Server.


repeat skobbler maps horizontally with OpenLayers?

I'm trying to get skobbler maps to repeat horizontally if you load this page http://developer.skobbler.com/demo/skobbler-tiles-openlayers.html and zoom out, the map will not repeat horizontally However if you use this link and zoom out http://developer.skobbler.com/web- tiles/embedded it works as I want (i.e. it repeats across the screen) However, I'm using OpenLayers and not an Iframe Is what I want to do possible with OpenLayers? Thanks Mark


Applying same color ramp over multiple maps with different values in ArcGIS10?

I'm doing a rainfall analysis for multiple events in ArcGIS10.1, and I want to apply the same color ramp (with a constant range) for each map. For example, all of the events I have range from zero to roughly 6.5 inches. The color ramp I created is a simple rainbow scheme, with red representing 0 and purpleish representing 6.5. My problem is is that obviously, not all of the events range from 0-6.5. Some go from 0.3-0.7, others from 0.1-1.3, etc. How do I go about applying the 0-6.5 color ramp on each map so that the values are represented with the proper color? Let me know if this question needs any clarification. Thanks! EDIT -- So perhaps I'm doing something wrong. ![](http://i.imgur.com/TfEydve.jpg) I made the color ramp as shown (saved it as a layer file and transformed it into a raster) and when I go to import a symbology to one of my interpolated layers ![](http://i.imgur.com/CuTq03g.jpg) as shown, it stretches the entire color ramp from the lowest value to the highest value, not just the parts of the color ramp it "should". Any advice on what I'm missing?


Is it real to install PostgreSQL 9 on the virtual machine Windows Server R2?

I have Virtual Box with GuestAddition. And virtual machine "Windows Server 2008 R2" When I use real computer with Windows Server 2008 R2 I can install postgresql-9.2.3-2-windows. But when I install postgresql-9.2.3-2-windows on virtual machine I have error "Problem running post-install step. Installation may not complete correctly" Please, help me.


Why are the results different for the same point using GEOS Point.transform and Postgis ST_Transform?

I'm trying to define a default point setting for OSM layers in my application and I want users to be able to define the default latitude and longitude in _EPSG:4326_. What I'm doing right now is to use GEOS `Point.transform` to turn the EPSG:4326 input into _EPSG:3857_ :               _default_point = Point(         getattr(settings, 'DEFAULT_LON', None) or 0,         getattr(settings, 'DEFAULT_LAT', None) or 0,         srid=4326     )     _default_point.transform(3857)     default_lat, default_lon = _default_point.get_coords()      Then, i.e. for `(-2.41544723477347,43.2140883725249)` I'm getting `(-268966.0382412352, 4810570.312724996)`. That, as for what I can see in the map, is wrong. Now, when trying over psql               SELECT ST_AsText(ST_Transform(ST_SetSRID(ST_Point(-2.41544723477347,43.2140883725249),4326),3857));                      st_astext                      -------------------------------------------      POINT(-268886.356213004 5344615.26048709)      Which is right. What is going on then? Why are the results different? EDIT: `select srid, proj4text from spatial_ref_sys where srid = 3857 or srid = 900913 or srid = 4326` yields:                  4326 | +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs         3857 | +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +units=m +k=1.0 +nadgrids=@null +no_defs      900913 | +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs


Open Arabic Shapefile in Qgis

I have a shape-file all it's fields are in Arabic. When I opens it in qgis 2.0 all the data is replaced by "?". I changed the encoding to Arabic or UTF-8 but with the same result. How can I Open an Arabic Shapefile in Qgis.


How to use GDAL with PyPy?

I am trying to optimize my python code using pypy. While I am able to import most libraries and run succesfully, I am not able to import GDAL libraries. The same libraries when imported in Python 2.7 run fine. I have tried adding GDAL location to PYTHONPATH, and even tried downloading GDAL Python bindings to PyPy directory, somehow pypy does not want to recognize GDAL, and comes with the following error.               Traceback (most recent call last):            File "app_main.py", line 51, in run_toplevel            File "ImportENC-P.py", line 50, in <module>              from osgeo import gdal, ogr, osr            File "B:\PYPY\site-packages\osgeo\__init__.py", line 21, in <module>              _gdal = swig_import_helper()            File "B:\PYPY\site-packages\osgeo\__init__.py", line 13, in swig_import_helper              import _gdal          ImportError: No module named _gdal


Geoserver WFS GetFeature returns no features when reprojecting PostGIS layer

I may be misunderstanding the GeoServer documents, but when accessing a GeoServer 2.2 WFS 1.1 layer, should I be able to reproject to an arbitrary SRS via the SRSNAME parameter? I'm generating a bounding box from the location of a click on a WMS map where I know there is a feature. The map is projected into UTM zone 15 (EPSG:32615). So I generate the following URL: http://server.domain.com/geoserver/MyWorkspace/wms?SERVICE=WFS&VERSION=1.1.0&REQUEST=GetFeature&TYPENAME=MyWorkspace%3AMyLayer&MAXFEATURES=9&SRSNAME=EPSG%3a32615&BBOX=382138.16774714%2c4587006.11489167%2c382138.16774714%2c4587006.11489167 This query returns 0 features, and I don't know why. It's essentially the same query I was using for another layer, also served from GeoServer, that was sourced from a shapefile in EPSG:32615, where the query worked. So here's what's different about this layer:   1. I'm using PostGIS for a backend on the new layer; the old one was a shapefile   2. The shapefile was in EPSG:32615; the PostGIS table is in EPSG:4326   3. The PostGIS table is using the "geography" data type rather than the "geometry" data type, if that's relvant. Any ideas?


I can't seem to access the plugin repository. Is this an issue with my system or is there something going on at the server end?! I get the following error as wells sometimes a window saying: "Error reading repository: QGIS Official Repository HTTP request failed" It's snowing like crazy here and that often messes with the internet but every thing else seems to be working fine.


openlayers load local gml

I'm trying to load a simple gml file with openlayers, but it doesn't load. Pretty sure it's a simple issue, but I figure out what's wrong. Thanks in advance, Pedro                           var meu = new OpenLayers.Layer.Vector("States", {                 protocol: new OpenLayers.Protocol.HTTP({                     url: "teste.GML",                     format: new OpenLayers.Format.GML()                 }),                 strategies: [new OpenLayers.Strategy.Fixed()],             });                  map.addLayer(meu)      Tried firefox and safari in mac os and firefox in linux


How to programmatically add a new plugin repository in QGIS?

How can I add the name of a repository and URL to the third-party plugin repository using python in QGIS 2.01?


For loop within SDE

I have a list of SDE connection files in a folder in my computer. I can use arcpy.listWorkspaces("*", "SDE") to find the workspaces. But I then want to dig deeper into each of the SDE files and list the feature classes. I tried:               sdeFiles=[]     listWorkSpaces=arcpy.ListWorkspaces("*", "SDE")     for sde in listWorkSpaces:         arcpy.env.workspace=sde         getFeatureClasses=arcpy.ListFeatureClasses()         getFeatureClasses.sort()         sdeFiles.extend(getFeatureClasses)     print sdeFiles      So when I print the sdeFiles list, it will list all the SDE feature classes in all the SDE files, but sometimes it will list them 6+ times, so it's duplicating it. I imported the list to excel and removed duplicates. It created 1514 duplicates and only 650 unique values. Frustrating! Can anyone tell me why?


How to display information in a div outside the map?

I tried a bit around and created a simple map with two polygons(most of the code I got from the Interactive Choropleth Map Tutorial). I added a `"L.control()"` to show information about the polygons in a `div`. Contrary to the tutorial I would like to have this `div` outside the map. To get a clue what I mean please look here: http://jsfiddle.net/Tbr9p/


postgis and osm2pgsql: authentication failed

Id like to import osm into my postgis to use it with qgis. I followed exactly these instructions: http://wiki.openstreetmap.org/wiki/PostGIS/Installation At step #Activate_PostGIS:               psql --username=gisuser --dbname=gis --command="\d"      This error occured:               psql: FATAL:  Peer authentication failed for user "gisuser"      (Same after importing my *.osm by doing this:               osm2pgsql -s -U gisuser -d gis path-to-file.osm      Error:               Error: Connection to database failed: FATAL:  Peer authentication failed for user "gisuser"      I have no idea what to do, so im asking here. Using xubuntu 12.10, postgres 9.1.11


How to use MOLUSCE plugin in QGIS?

I am trying to explore MOLUSCE ver 3.0.4 in QGIS Dufour. I am trying to simulate the land use change using this tool. I have following querries:   1. It asks for Initial and final year at the first step. Where to enter the exact year for which user is interested for prediction?   2. It gives output file of simulation. For which year this simulation is? Since we already have given Initial and Final year.


Open source/free software for generating a 3D view of mapping?

Is there any software that will generate a 3D view given a DTM/height data and raster mapping, and output an image in a printable resolution? I am looking for something similar to what you'd find in a commercial hike planning application, but free or open source preferably.


How to utilise single 'randomForest' model for several rasters in R?

I am using randomForest package in R to derive species distribution models. When I predict the model to the RasterStack containing the predictors in which the values were extracted to build my model, everything goes fine. The problem is that I would like to predict this model I generated in the present scenario to the same set of predictors in the past and future. However, R gives me an error message. Past and present predictors RasterStack are the same extension, resolution, etc. This is not the problem. Reading the help of the 'predict' function {raster package} in R, I found the following: "The names in the Raster object should exactly match those expected by the model. This will be the case if the same Raster object was used (via extract) to obtain the values to fit the model". So, according to this, I can't predict my model to another RasterStack that was not used to obtain the values for the model. So, I would like to know another way that I could predict values from a present scenario to past or future scenarios. Thanks a lot, Mariana


How to get hypsometric curve from elevation data?

I have my elevation data but cannot make an attribute table from it and so I'm not sure how to get any of that data into a table in order to make a hypsometric curve.


Missing level for cache map service

I am having this issue where caching is missing one of the levels (level 13). I am trying to find out how to resolve this issue. L9 - 1,155,581 - OK L10 - 577,790 - OK L11 - 288895 - OK L12 - 144447, OK L13 - 72223, MISSING, don’t show at all for some areas L14 - 36111, OK L15 - 18055, OK L16 - 9027, OK L17 - OK L18 - OK Thanks Jay


Disappearing Point Features in Geoserver using WMS

I have a shapefile with roughly 6,500 points worldwide in it that I am trying to serve with Geoserver 2.2.1 using WMS. All has been apparently fine with it until I implemented a filtering feature on my client application which uses leaflet. When I append a CQL_FILTER (attribute filter, not spatial) into the WMS query, I noticed missing features when zoomed out. When I zoomed back in, they would sometimes reappear, but not always. See the image below - ![Side by side comparison](http://i.stack.imgur.com/IoA2B.png) At the zoom level on the left, Atlanta is not displayed. When I zoom in, it is. However, sometimes not even the point in Tampa is displayed at the zoom level on the left. If I zoom out 3 more levels, no points at all will be displayed. I'm not sure the issue is the CQL_FILTER param, as with 6,500 points it's hard to notice a few missing points at a global scale, but the particular filter I'm showing here as an example filters down to only 3 features, and when 1 to 3 of them are missing depending on the zoom level it's particularly noticeable. I can recreate this same behavior in geoserver's own layer preview using the CQL Filter in the advanced map options, so I'm fairly certain it is not my client application. As for my own steps at debugging this issue, I have tried the default Points style/sld file to rule out my own layer styling. I have turned off all the caching that I am aware of. I double checked that my projections are correct - I created the shapefile in ArcGIS 10 using WGS_1984_Web_Mercator_Auxiliary_Sphere as the projection, and the layer is set to EPSG:3857 in geoserver, which I think is equivalent. I also upgraded from geoserver 2.2 to 2.2.1 and have had the same issue in both. I have also deleted geoserver's spatial index file (.qix) and let it get recreated since I've seen similar issues in Arc with corrupt spatial indexes, but obviously that did not work either. Here's a snapshot from Geoserver's own layer preview with the CQL Filter on and zoomed into the same area as shown above. The red circle is approximately where I should be seeing another point (Atlanta). ![Openlayers Example](http://i.stack.imgur.com/Zutqm.png) I've tried tweaking every other setting I can think of, but am having no luck. I've also looked up and down the geoserver logs and enabled verbose logging, and do not see any errors/exceptions. I also do not see any bad requests in chrome's developer tools. Thanks in advance for any help. If I'm missing any critical info, I'll provide what I can, but this is an internal/non-public application.


Multiple queries in a row

I have osm data imported into PostgreSQL with osm2po tool (wich makes the data routable). Also I have imported data relationed with nodes and node_tags with osmosis tool. I need to obtain a set of nearest nodes from a given location (maybe 30 or 40). Then I need to calculate the shortest path from the given node from each others. It means to execute shortes_path query multiple times in a row, wich takes about 2 seconds (1 query 2 seconds: 30 queries 60 -> too long).               SELECT * FROM shortest_path('        SELECT id AS id,               source::int4 AS source,               target::int4 AS target,               cost::float8 AS cost        FROM tt_2po_4pgr',     275742,     274494,     false,     false)      What would be the best way to do this in a reasonable time? Thank in advance! Pd: excuse my english **Edited:** In accordance to _dkastl_ answer, I tried with wrapper functions and effectively I obtain better results 2 or 3 times better, but isn't acceptable for me yet. With the following query:               explain analyze SELECT *         FROM dijkstra_sp_delta('dr_2po_4pgr', 278174, 289786, 0.3);      I obtain this results: explain


Left click then right click is treated as double click in Mac QGIS Layers

I am having an issue with QGIS 1.8.0 on my Mac where I left click then right click which gets treated as a double click if the timing is too close together and brings up the Layer Properties box. This usually happens when I am focused on a different program, left click to gain QGIS' focus and then right click to get a context menu for a layer. This is force of habit I guess.   1. Is there a way to disable the Layer Properties from appearing with a double click? I know that you can choose between the Attribute Table and the Properties Box.   2. Should this issue be added as a bug report to QGIS?


Enabling ArcSDE Archiving with Python and C API using ArcGIS 10.0?

I have a script that does some register as versioned and updating privileges. How would I call into the c api for ArcSDE to enable archiving with python? http://help.arcgis.com/en/geodatabase/10.0/sdk/arcsde/api/capi/versioning/archcreate.htm Any direction would be helpful.


why args.ScreenPoint.IsEmpty is true?

I use a map control and set a MapClick event for it. And i have a tool that ClientAction set to point. Event fired when i clicked on the map. But args.ScreenPoint was empty. What's wrong ? I want get x and y in map click.


How to create a line with start and end points which are in another layer in QGIS?

In my map, there are lots of points in layer **A**. I want to make a line layer **B** and I want to make start & end points of lines snapping to points in layer A manually. Can I use points of layer A as a coordinate attribute selection index for my line layer B? Thank you.


Very slow connections to ArcSDE with over 20,000 datasets

I have a problem with ArcSDE that I have never come accross before and I was wondering if anyone could shed any light on this problem. **Problem ArcSDE 9.3.1:** > Connnecting to an ArcSDE instance (Oracle Spatial) with over 20,000 datasets > (raster) causes ArcGIS/ArcObjects code to connect very slowly - minutes in > fact before a successful connection is made to the database. However if a > native SQL connection is made, the connection is almost instant. Has anyone else had this problem and if so what can be done to remedy this problem.


Command doesn't show up on ICommandBar

I'm trying to put a custom command into the Help menu on startup if it isn't already there. This is my code to do so:               private void OnCreate()     {         // other stuff in OnCreate()...         // registering my method for the right event         m_appStatusEvents = m_application as IApplicationStatusEvents_Event;         m_appStatusEvents.Initialized += new IApplicationStatusEvents_InitializedEventHandler(AddCmd);     }          private void AddCmd()     {         UID myCmdId = new UIDClass();         myCmdId.Value = myCmdIdValue;  // myCmdIdValue is my command's UID              ICommandBars commandBars = m_application.Document.CommandBars;              UID helpMenuId = new UIDClass();         helpMenuId.Value = helpMenuIdValue; // helpMenuIdValue is esriArcMapUI.MxHelpMenu         ICommandBar helpMenu = null;         helpMenu = (ICommandBar)commandBars.Find(helpMenuId, false, false);         // here, helpMenu.Count is 3         if (helpMenu != null)         {             ICommandItem myCmd = helpMenu.Find(myCmdId, false);             if (myCmd == null)             {                 Object index = Type.Missing;                 helpMenu.Add(myCmdId, ref index);                 // here helpMenu.Count is 4, so my command was added             }         }     }      I think it's pretty straight forward: I'm looking for the Help menu and add my command if it isn't already there. In the debugger, I can see that an item is added to the Help menu, as its `Count` goes from 3 to 4. Problem is: The command doesn't show up. I can add it manually, so I know it's loaded and correct, it just doesn't show up when I try to add it automatically. Edit: I just tried to refresh the Help menu with `(helpMenu as ICommandItem).Refresh()` but that didn't help either. Edit2: I just added this to debug:               for (int i = 0; i < helpMenu.Count; i++)     {         Debug.WriteLine(helpMenu[i].Caption);     }      And it shows   * &ArcGIS Desktop-Help   * &ArcGIS Resource Center   * &About ArcMap...   * myCommand Caption So it's definitely there, it just won't show!


How to change default zoom in the Google embedded map?

I am having difficulties to change the default zoom level when I open the iframe below in my site. Can anyone please help me to set zoom level when following iframe is opened?               <iframe src="http://mapsengine.google.com/map/u/0/embed?mid=zYk_BmXXTdxw.k83snDbX5SJc" width="745" height="480"></iframe>      Thank you so much everyone in advance. PS: This is the newest version of Google Maps.


How to read vector datasets using GDAL library

Up to now, I have used GDAL in order to read raster information and geotiffs and show them with the help of OpenGL. Now I have a vector dataset, and by that, I mean I have the files with types: **DBF File , PRJ File , SBN File , SBX File , SHP File , LOCK File , XML Document , SHX File** I want to read the vector information **(e.g vertexes of the geometries and geometrical and topological relationships between feature)** and then show them using OpenGL.   * What part of GDAL library should I study?   * Can you introduce me a tutorial link that has fully covered what I want to read from the vector dataset?


please Introduce me a free 3rd party library appropriate for showing Geospatial information

Up to know I have used GDAL to extract raster data from geotiffs and then opengl's used to show that raster data.   But when it comes to extract vector data with the help of ogr library, what library would you suggest to show these vector data in visual-studio picture controls opengl that is more appropriate than opengl? * * * I've found out from your comments that you have thought I'm searching for a way to read vector data. this is the subject of my another question. Here, I'm just asking suppose that we have extracted all the geometries. Which environment is the best to show these geometries?


Is there a work-around for a completely_within spatial join in ArcGIS? This is bugged (confirmation by Esri Inc)

I tried to perform a 1-to-many spatial join between polygons (target feature) and lines (join features) to join the lines that are completely covered by a polygon to that polygon. (the polygon is actually a result of buffering the lines with a dissolve option) I used the option 'completely contains', but this returned a background processing error. After sending parts of my data to Esri Inc, I got the confirmation that there is indeed a bug they need to fix (#NIM086228). Anyway, my project needs to be completed. In search of a work-around, I already tried to turn the lines into polygones (small buffer around it), but it returns the same error. When trying the contains_clementini (or actually all contains and within options) suddenly give me an error even before running it. The error is 'Relationship invalid for selected layers'. Here is what I need to do: > I have (pipe)lines that I want to 'merge' when they have common > characteristics. So I create a buffer and let it dissolve when the value of > those characteristics is equal. A polygon layer is created, but there is not > link back to the underlying lines that are the basis of the polygon. In > order to have these ID's, I wanted to perform the spatial join with match > option 'completely contains Has someone any experience with the last type of invalid error messages? Or another work-around for me? Thank you already! Mieke


Scalar relationships ArcGIS

I have point data which has a number of fields/columns, each of which is at a different scale. I'm looking to find a way of constructing a very simple index which would combine this data. Say I have 3 levels of scale (so 3 fields) - small, medium, large. Each of the points either has or hasn't got a text string at each of those levels. I'm looking to assign a number (or not) depending on whether the text string is present or not, and aggregate that data. Basically I want to see if there are patterns in the presence/absence of text strings at these different scales (possibly re-combining the value given to each of them in more than one way). I'll then figure out a way of showing this map format (Arc). I don't think I've been very clear but please ask if you need clarification. Here's a sample table with something along the lines of what I was thinking of doing (in grey). I'm not sure this will give me the outcome I need - ideally each column should have a different weighting reflecting relative size. But should give a rough idea of the sort of thing I was thinking of doing. ![enter image description here](http://i.stack.imgur.com/y01Xt.png)


How to Interpolate through ModelBuilder?

I would like to create an iterated interpolation tool through ModelBuilder. I have an excel file containing different columns and each column represents a data that will be interpolated using Spline Tool. How can I do this through ModelBuilder? I am using ArcGIS 10. ![enter image description here](http://i.stack.imgur.com/d0t0m.jpg) Image contains point data with lat/long and sample# are the columns to be interpolated.


GPS suitability

I was wondering what type of GPS/GNSS receiver should I use for recording and mapping trees in a particular farm (inventorying trees). Will handheld GPS like Garmin do the job? Or should I go for the expensive RTK system? Please help. *Remark: The area is a mountainous region.


How to package up a Python script for an end user to run

I have a Python script that I have made to some basic data management tasks. I want to package up a Python script in a nice user friendly way for an end user to run, they do not know how to run Python scipts and I need them to input some parameters in an easy way i.e. say a filename - whats the best way to do this. Cheers


How to cache online tiles on a development machine?

I am developing a slippy map web application using openlayers. Right now my tiles come straight away from OSM and google, etc. For the moment I always work on the same area for development purposes and I would like to have those tiles cached on my dev machine so I can work on the web app offline yet still see some map and to reduce network consumption. What I would like is to set up some sort of proxy which would just cache tiles I watch and if they are not on the local machine download them or show a dummy tile if I am offline. I thought tilecache was the tool for this, but I get quite confused trying to set it up. It seems more appropriate for mapserver generating tiles from map data. I know firefox does some caching, but I would like some more control and be able to clear that cache to make sure I see the last update of my work. What is the most straightforward way to achieve this?


Elevation statistics for a given Area

I am attempting to determine the best way to calculate the Minimum and Maximum elevation for a defined area. The area is defined by an `IPolygon` The elevation is defined as an `ISurface` ## **Possible Solutions:**   1. The only way I can think of is to check every cell of the surface. Is there a method that performs this?   2. In the comments, @whuber, suggested a Zonal Maximum. Expanding on this idea, I've begun looking into the `IZonalOp`. _Are there other options?_


Feature linked annotation subtype removal in migrating to ArcGIS 10

I have a question about emergent properties of deleting an annotation subtype. If we delete an annotation subtype how will this effect the annotation feature class, will there be any changes to the visible annotation the active data frame? I haven’t seen any changes with the removal of the subtype, but it seems too good to be true. Background info: In our enterprise implementation, we are in the process of upgrading from a 9.2 to 10 data format, to help achieve this we came up with a Python script, turns out there are quite a few databases that fail at the arcpy.UpgradeGDB_management() method. The exception is:"The feature class' default subtype code cannot be retrieved or is invalid” In order to resolve this issue we removed the Annotation’s subtype and the script will run just fine. Thanks in advance.


Why are ArcGIS 10.2.1 for Server services stuck in perpetual Stopping... state?

Running into a problem repeatedly with our ArcGIS Server installation where all services are in a perpetual state of "Stopping..." whenever we do an upgrade (just upgraded from 10.2 to 10.2.1). The same thing occurred when we changed the ArcGIS Server account from a local account to a domain account. The services cannot be stopped and started manually. Restarting the ArcGIS Server service does not seem to help either. Rebooting the server that ArcGIS Server is running on does not help with this issue either.


How can I change the size of the infowindow in CartoDB?

I'm working with info windows, but I'd like them to be a bit bigger than the defaults in CartoDB. To display pictures and additional information. How can I make it so they are bigger? Is there any way in the GUI to do this? Is there a more advanced way that has some template I could easily apply?


Multiple geometry in a cursor SQL

I have a cursor which I create from a table which has one geometry field only. But in my cursor I would be using two geometry fields as I need to calculate two different areas.Below is the sql statement which I am using . After that I need to insert those two geometries in the table which also has two geometry columns. But when I try to do the same it gives me an error **"inconsistent datatypes: expected NUMBER got MDSYS.SDO_GEOMETRY""** How could I fix this particular geometry issue.                 Declare          Cursor c1 is           select 0 as BG_ID,store_number ,geometry , geometry bg_geometry,0 as total_area          from   TMCS_All_Stores_TA  where CLIENT_ID= 1 and ORG_ID = 1 and  rownum <= 1;         Type C1_TAB_TYPE is table of c1%ROWTYPE;                c1_list c1_TAB_TYPE;               Cursor c3 is           select BG_ID,0 as BLK_ID ,store_number,0 as total_area_blk           from   TEMP_PROP_AREA_BGLEVEL ;          Type C3_TAB_TYPE is table of c3%ROWTYPE;                c3_list c3_TAB_TYPE;                      Begin                     For r1 in c1                       Loop                 Select  BG_ID,store_number,geometry,bg_geometry,total_area                  Bulk Collect into c1_list from(                      SELECT a.BG_ID,r1.store_number,(SDO_GEOM.SDO_INTERSECTION(a.geometry,r1.geometry, 0.005)) GEOMETRY,r1.bg_geometry bg_geometry,                      ( (SDO_GEOM.SDO_AREA(SDO_GEOM.SDO_INTERSECTION(a.geometry,r1.geometry, 0.005), 0.005, 'unit=sq_mile') /SDO_GEOM.SDO_AREA((a.geometry), 0.005, 'unit=sq_mile'))*100) Total_area                      FROM   TMCS_BG_STI_DATA a                       Where sdo_relate(a.geometry, r1.geometry, 'mask=anyinteract') = 'TRUE');                                             For j in 1..c1_list.count loop                                        DBMS_OUTPUT.PUT_LINE('block group geometry--> '||c1_list(j).store_Number);                                        insert into   TEMP_PROP_AREA_BGLEVEL                                        values c1_list(j);                                         End loop;                             For j in 1..c1_list.count loop                         DBMS_OUTPUT.PUT_LINE('area--> '||c1_list(j).total_area);                       if c1_list(j).total_area !=100 then                                             select BG_ID,BLK_ID ,store_number,total_area_blk                                 BULK COLLECT INTO c3_list from                              (                              SELECT c1_list(j).BG_ID BG_ID,a.BLK_id,c1_list(j).store_number store_number,                              ( (SDO_GEOM.SDO_AREA(SDO_GEOM.SDO_INTERSECTION(a.geometry,c1_list(j).geometry, 0.005), 0.005, 'unit=sq_mile') /SDO_GEOM.SDO_AREA((r1.bg_geometry), 0.005, 'unit=sq_mile'))*100) total_area_blk                              FROM   TMCS_BLK_STI_Q213_DATA a                               Where sdo_relate(a.geometry, c1_list(j).geometry, 'mask=anyinteract') = 'TRUE'                              And BG_ID = c1_list(j).BG_ID);                                               For i in 1..c3_list.count loop                                       --   DBMS_OUTPUT.PUT_LINE('geo_id--> '||c3_list(i).Blk_id);                                       --    DBMS_OUTPUT.PUT_LINE('area--> '||c3_list(i).total_area_blk);                                          insert into   temp_Prop_area_BG                                          values c3_list(i);                                          End Loop;                            End if;                   End Loop;                        Insert into   temp_prop_cal                   (bg_id,  store_number, total_area)                   SELECT bg_id,  store_number,  SUM(total_area)                   FROM        (                   SELECT bg_id,  store_number,  total_area as total_area          FROM   TEMP_PROP_AREA_BGLEVEL a          WHERE store_number='336655'          UNION ALL          SELECT bg_id,   store_number,  total_area_blk as total_area          FROM   temp_Prop_area_BG b          WHERE store_number='336655')        GROUP BY bg_id, store_number;              End Loop;           End;


How to populate ESRI GeoPortal programmatically?

I'm trying to find either some documentation or an example of how to manually populate the GeoPortal metadata datastore. I'm using a MySQL backend. Given a ISO XML file I can manually harvest it but this is tedious and doesn't fit with our model. I can easily write a script to populate entries in the SQL database but it appears that there is more to it than just adding entries to SQL, there is a part that goes into the Lucene engine as well. If anyone has any experience or can provide a simple example, how-to or preferably some documentation on populating the datastore it would be greatly appreciated. The goal is to have GeoPortal as an external catalog of locally registered datasets, the trouble is registering the datasets with the GeoPortal catalog in a preemptive manner rather than a reactive (having an administrator manually "harvest" or ingest resources).


Trying to write a python script to crawl a directory and list all datasources

I'm editing a scrip from ESRI to try and create a list of all datasources within MXD's. I've got the script to crawl the file tree and find all the info, but I would like it to output the data for each layer/datasource on a single line, so it is easy to sort in excel. If anyone could take a look at the code I have so far and let me know how to accomplish that, I would appreciate it.               import arcpy, datetime, os          try:          arcpy.gp.overwriteOutput = True          #Read input parameters from GP dialog     folderPath = arcpy.GetParameterAsText(0)     output = arcpy.GetParameterAsText(1)          #Create an output file     outFile = open(output, "w")          #Loop through each MXD file     count = 0     for root, dirs, files in os.walk(folderPath):         for filename in os.listdir(folderPath):           fullpath = os.path.join(folderPath, filename)           if os.path.isfile(fullpath):             if filename.lower().endswith(".mxd"):                      #Reference MXD                 mxd = arcpy.mapping.MapDocument(fullpath)                 count = 1                      #Write MXD data to file                 MapDoc = os.path.basename(mxd.filePath)                 MapDocPath = mxd.filePath                       #Reference each data frame and report data                 DFList = arcpy.mapping.ListDataFrames(mxd)                 for df in DFList:                     #Format output values                     if df.description == "": descValue = "None"                     else: descValue = df.description                               #Reference each layer in a data frame                     lyrList = arcpy.mapping.ListLayers(mxd, "", df)                     for lyr in lyrList:                         lyrName = lyr.name                          if lyr.supports("dataSource"):                             lyrDatasource = lyr.dataSource                                                      else: lyrDatasource = "N/A"                      outFile.write("Map Document, MXD Path, DataFrame List, Description, Layer name, Layer Datasource")                 str = ",";                 seq = (MapDoc, MapDocPath, DFList, descValue, lyrName, lyrDatasource);                 outFile.write(str.join(seq))                      del mxd                                outFile.close()          #Open resulting text file     os.startfile(output)          #Delete variables that reference data on disk     del folderPath, output, outFile, fullpath          except Exception, e:     import traceback     map(arcpy.AddError, traceback.format_exc().split("\n"))     arcpy.AddError(str(e))      The error I recieve when I try to run this is as follows:               \MultiMXDReport_crawl.py", line 60, in <module>     outFile.write(str.join(seq))     TypeError: sequence item 2: expected string or Unicode, list found     <type 'exceptions.TypeError'>: 'str' object is not callable     Failed to execute (MultiMXDCrawler).


Select by Location tool

I am using the Select by Location tool to select certain polygons on my shape- file layer. I start by having a given point as my Selecting Feature and that its relationship is completely_ contains. After I run the tool the first time, I get the polygon feature that my point is in (lets call it Polygon A). When I run another select by feature tool that its relationship is boundary_touches my Polygon A, I get a layer that also contains Polygon A and the features that share boundaries with it. Is there a way I can prevent this from happening? Meaning that every time I run the select by location tool it doesn't select polygons from the previous time. Or is there a tool or a way that will allow me to join multiple layers together and remove duplicate features? I have tried the Union tool as well as the Intersect tool. Any suggestions? Thanks


Using American Fact Finder to download data for all block groups in a state at at time?

I am having issues using the advanced search on the American Fact Finder to download census block data for a whole state at a time. The best option that it gives me, when trying any route, is to select geographies of all census block groups within a county within a state. This leaves me with the option then of selecting every county, but in some states (like Texas) that have more than 250 counties, I am out of luck with this method. Does anyone have a solution for this? To duplicate my error:   * Visit http://factfinder2.census.gov/faces/nav/jsf/pages/index.xhtml   * Then select Advanced Search,   * Then select on the left "Geographies",   * Then select "all geographic types" radial button,    * Then in the dropdown, select Block Group or part,   * Then select state in the next drop down, and you will start to see the issue. Also, if instead of using the list search and you use the tab that says "name" I run into the same issue. Does anyone have any experience with American Fact Finder downloads? I would appreciate any suggestions, Thanks! Also, I am not asking the same question here as I am looking to see if this is possible just by accessing AFF not their ftp. Gathering data for a large number of census block groups


topology problem

Can anyone tell me how to fix an error after running the topology rule in ArcGIS? I tried to subtract and create a new layer, but a failed warning appears every time.


Symbology of contour lines in ArcMap

I have contour lines with height values from 0 to 1030 (sorted by height 110, 115, 120, 125). Values with in-between number (like 5, 15, 115, 135...) should be grouped and colored with same color, and other values should be grouped this way 10, 60, 110, 160, 210, 260... (yellow), 20, 70, 120, 170, 220... (red), 30, 80, 130, 180, 230... (cyan) etc. I hope you get it despite my bad English :) So, I don't know how to change symbology in ArcMap to color them that way. Can somebody help me? Thank you.


ArcGIS Server 10 and Desktop 9.3 interoperability

I want to install ArcGIS Server 10 but our users only have ArcGIS Desktop 9.3.1. So does ESRI support:   * a) Access/use AGS services from 10 in 9.3.1?   * b) Manage (create services etc) AGS 10 from ArcCatalog 9.3.1?  I know you can use the web manager but our users are used to manage AGS from desktop.


How can I use fortify() to create a filtered R data frame from a shape file?

I am in the process of building a `ggplot` choropleth map of population in administrative areas in Wales. I have downloaded the Boundary-Line data from the Ordnance Survey and extracted what seems to be the right shapefile (community_ward_region.shp). Using R, I have got as far as reading in the shapefile.               require(maptools)     shape <- readShapePoly(wards)     str(shape)      Which gives me this promising output:               Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots       ..@ data       :'data.frame': 1690 obs. of  4 variables:       .. ..$ NAME      : Factor w/ 1507 levels "Abbey Cwmhir",..: 969 90 111 200 441 477 1455 249 255 305 ...       .. ..$ DESCRIPTIO: Factor w/ 4 levels "COMMUNITY","COMMUNITY WARD",..: 2 2 1 1 2 2 2 2 1 1 ...       .. ..$ COMMUNITY : Factor w/ 858 levels "Abbey Cwmhir",..: 67 67 81 128 152 152 152 152 157 190 ...       .. ..$ FILE_NAME : Factor w/ 23 levels "ABERTAWE_-_SWANSEA",..: 1 1 1 1 1 1 1 1 1 1 ...       .. ..- attr(*, "data_types")= chr [1:4] "C" "C" "C" "C"       ..@ polygons   :List of 1690       .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots       .. .. .. ..@ Polygons :List of 1       .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots       .. .. .. .. .. .. ..@ labpt  : num [1:2] 259009 188524       .. .. .. .. .. .. ..@ area   : num 1923892       .. .. .. .. .. .. ..@ hole   : logi FALSE       .. .. .. .. .. .. ..@ ringDir: int 1       .. .. .. .. .. .. ..@ coords : num [1:1629, 1:2] 259413 259420 259427 259427 259432 ...       .. .. .. ..@ plotOrder: int 1       .. .. .. ..@ labpt    : num [1:2] 259009 188524       .. .. .. ..@ ID       : chr "0"       .. .. .. ..@ area     : num 1923892      Now if I do this:               bar <- fortify(shape, region = "NAME")      I get a nice data frame called `bar` that looks pretty much as I expected:               > str(bar)     'data.frame':   4744053 obs. of  7 variables:      $ long : num  302962 302970 302974 303013 303015 ...      $ lat  : num  280066 280076 280078 280097 280105 ...      $ order: int  1 2 3 4 5 6 7 8 9 10 ...      $ hole : logi  FALSE FALSE FALSE FALSE FALSE FALSE ...      $ piece: Factor w/ 29 levels "1","2","3","4",..: 1 1 1 1 1 1 1 1 1 1 ...      $ group: Factor w/ 1762 levels "Abbey Cwmhir.1",..: 1 1 1 1 1 1 1 1 1 1 ...      $ id   : chr  "Abbey Cwmhir" "Abbey Cwmhir" "Abbey Cwmhir" "Abbey Cwmhir" ..      However, this is a large data frame and `ggplot` runs out of puff when trying to display it. In reality I only want to look at one area at a time. It looks as if the `FILE_NAME` factor in the shape object is what I want as it mostly corresponds to counties and the major conurbations.               > unique(shape@data$FILE_NAME)      [1] ABERTAWE_-_SWANSEA      [2] BLAENAU_GWENT_-_BLAENAU_GWENT      [3] BRO_MORGANNWG_-_THE_VALE_OF_GLAMORGAN      [4] CAERDYDD_-_CARDIFF      [5] CAERFFILI_-_CAERPHILLY      [6] CASNEWYDD_-_NEWPORT      [7] CASTELL-NEDD_PORT_TALBOT_-_NEATH_PORT_TALBOT      [8] CONWY_-_CONWY      [9] GWYNEDD_-_GWYNEDD     [10] MERTHYR_TUDFUL_-_MERTHYR_TYDFIL     [11] PEN-Y-BONT_AR_OGWR_-_BRIDGEND     [12] POWYS_-_POWYS     [13] RHONDDA_CYNON_TAF_-_RHONDDA_CYNON_TAFF     [14] SIR BENFRO - PEMBROKESHIRE     [15] SIR_BENFRO_-_PEMBROKESHIRE     [16] SIR_CEREDIGION_-_CEREDIGION     [17] SIR_DDINBYCH_-_DENBIGHSHIRE     [18] SIR_FYNWY_-_MONMOUTHSHIRE     [19] SIR_GAERFYRDDIN_-_CARMARTHENSHIRE     [20] SIR_YNYS_MON_-_ISLE_OF_ANGLESEY     [21] SIR_Y_FFLINT_-_FLINTSHIRE     [22] TOR-FAEN_-_TORFAEN     [23] WRECSAM_-_WREXHAM     23 Levels: ABERTAWE_-_SWANSEA ... WRECSAM_-_WREXHAM      Q. How can I select only a subset of the data from the `shape` object I extract from the shapefile? For example, only the `POWYS_-_POWYS` parts? If I can somehow include the `FILE_NAME` in the data frame that is created with `fortify` then I could easily subset the `bar` data frame but I don't know how to do that. Or is there a way to use `fortify` to extract only parts of the object?


How to pass an output stored in a single variable into another script using arcpy

I am running set of processes each with different outputs in a python script. Is there a way to set a variable as global so that just the output of that single step can be accessed from other scripts?


Python Error while installing Cartogram

I am quite new to GIS and qGIS especially. I would like to install cfarmers cartogram plugin. From what I understand, this has to be done manually - so I 've followed the instructions from RobinLoveLace (see link below). Where has the cartogram plugin for QGIS gone? When I try to click the icon for the Cartogram Creator, I get an error message (see below). Any tips on how to get it to work? > An error has occured while executing Python code: > > Traceback (most recent call last): File > "C:/PROGRA~2/QGISCH~1/apps/qgis/./python/plugins\cartogram\plugin.py", line > 136, in run d = doCartogram.Dialog(self.iface) File > "C:/PROGRA~2/QGISCH~1/apps/qgis/./python/plugins\cartogram\doCartogram.py", > line 61, in **init** for name, layer in layermap.iteritems(): File > "C:/PROGRA~2/QGISCH~1/apps/qgis/./python/plugins\cartogram\doCartogram.py", > line 61, in **init** for name, layer in layermap.iteritems(): File > "C:\PROGRA~2\QGISCH~1\apps\Python27\lib\bdb.py", line 49, in trace_dispatch > return self.dispatch_line(frame) File > "C:\PROGRA~2\QGISCH~1\apps\Python27\lib\bdb.py", line 68, in dispatch_line > if self.quitting: raise BdbQuit BdbQuit > > Python version: 2.7.4 (default, Apr 6 2013, 19:54:46) [MSC v.1500 32 bit > (Intel)] > > QGIS version: 2.4.0-Chugiak Chugiak, 8fdd08a > > Python path: ['C:/PROGRA~2/QGISCH~1/apps/qgis/./python/plugins\processing', > 'C:/PROGRA~2/QGISCH~1/apps/qgis/./python', u'C:/Users/HEC/.qgis2/python', > u'C:/Users/HEC/.qgis2/python/plugins', > 'C:/PROGRA~2/QGISCH~1/apps/qgis/./python/plugins', > 'C:\PROGRA~2\QGISCH~1\bin\python27.zip', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\DLLs', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\plat-win', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\lib-tk', 'C:\PROGRA~2\QGISCH~1\bin', > 'C:\PROGRA~2\QGISCH~1\apps\Python27', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\GDAL-1.11.0-py2.7-win32.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages\PIL', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\jinja2-2.7.2-py2.7.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\markupsafe-0.23-py2.7-win32.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\python_dateutil-2.2-py2.7.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\pytz-2014.2-py2.7.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages\win32', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages\win32\lib', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages\Pythonwin', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\Shapely-1.2.18-py2.7-win32.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages\six-1.6.1-> > py2.7.egg', 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\wx-2.8-msw-unicode', > 'C:\Users\HEC\.qgis2\python\plugins\mmqgis/forms', > 'C:\PROGRA~2\QGISCH~1\apps\qgis\python\plugins\fTools\tools'] > The plugin for Cartogram Creator for qGIS is downloaded from here: https://github.com/cfarmer/cartogram-plugin/tree/master


Cartogram plugin for leaflet or OSM

Does anyone know if there is a plugin to make cartograms for leaflet ? I have sets of geotagged Twitter datas and I'd like to distord the map according to the density of tweets of every area of the world. As cartograms are mostly changing the tiles I'm not sure if Leaflet is the most interesting plugin to do so (as it's mostly done to display datas on a static map layer like an OSM one), but as I begun to use it a lot for other purposes it would be easier to implement in my code... Another approach would be to find a plugin to transform directly the OSM layer without using Leaflet, but I don't know if that exists too. Thx in advance !


Cannot view published map content on ArcGIS Online outside of firewall

I successfully published and added a map onto ArcGIS Online using ArcGIS for Server, which is located on a secure server. I can view the map just fine within the server environment (with localhost:6080) and from my work desktop (with [servername]:6080). However, when I try to view the same map from home, it fails to find http://[servername]/arcgis/rest/services/[mapname]/MapServer. Is there any way to bypass this firewall issue without having to obtain a web server? I do not have server admin rights, and obtaining permission for installing a web server would take longer than the timeframe that I'm working under. (Would a web server even resolve this issue?)


GIS Analysis - best way to combine and display demographic data?

I am currently undertaking a demographic research project for my employer. We have the data on income, age, sex, language, transport etc and I have mapped most of it and created some thematics of said variables. Whilst the maps are perhaps adequate from a raw data representation perspective, I am looking for a technique or thematic style which displays statistical and demographic data in a more intuitive and interpretive fashion. Here is an example of what i would like to improve: ![http://i894.photobucket.com/albums/ac142/mellingham/Income_zps95372b84.png](http://i.stack.imgur.com/aXzp1.png) **Is there a specific technique or method that people are using which displays thematic data which is more easily interpreted by a non GIS user?** Similar to this **question** but perhaps with more technical descriptions on how to achieve some of the better cartographic outputs. **EDIT: Thanks Ian who suggested using rasters to generate a heat map type approach:** ![Heat Map of Non-English speaking families](http://i.stack.imgur.com/mWlnq.png)


How to validate geometry in openlayers?

How can I validate the geometry in Openlayers? I need to check " **MULTILINESTRING()** " (which is empty) is valid/not valid. In postGIS, for Geometry validation we are using " **isValid()** " function. Like this, is there any options in Openlayers?


Efficient algorithms for recovering watershed boundaries?

This is somewhat related to a previous question I asked: What Algorithm is used by ArcGIS Watershed tool? I was wondering if there is a more efficient way to calculate catchment areas if you are interested only in the polygon defining the catchment area (rather than all cells contributing to the catchment area). Doing a recursive search seems very expensive and would return a lot of information we don't want (namely each cell within the area). Assuming the same starting point - a raster grid and a pour point/outlet (I leave it open for now whether the raster is a DEM, Flow Direction or Flow Accumulation). My first thought was to simply follow the 'steepest path' upstream of the outlet point, which would quickly take us to the 'edge'. But I could not think of a way to follow the 'edge' around. Are there any other possibilities?


Linking up GPS points taken from rail network

Bit of a newbie here, trying to work out whether any PostGIS functions can help with my task. Apologies in advance if my use of GIS terminology is not up to scratch yet. I have a large database of GPS points, ie lat/long coordinate pairs (large being "millions" of points). Currently they have no relational information, so they are not linked together. However, they were collected over the UK rail network, typically every few metres. So all points map to rail tracks. I need to process these points and link them together, so they map onto real railway lines. I want to identify which GPS points are buffers (ie ends of track) and points (in the rail sense - ie where track splits) as well as those that are "just sections of track". This needs to be a process that runs automatically with no user intervention. There are many areas where multiple tracks run in parallel, or areas such as large stations with lots of points & sidings, where there are large numbers of GPS points in a small area. Being rail, there is a limit to how much the line will curve, so that must be a key feature to aid analysis. As the GPS points are close together, it should be possible to see which form lines that could be track, and dismiss those that turn too sharply (which must be points from parallel tracks). Initially I had been planning to write some java code that would pick GPS points from the DB that fall within a limited geographical area of each other, and look for GPS points within that region that continued in a straight-ish line. The add/update "previous node id" and "next node id" fields for each point. Then repeat, using the last linked node as the centre of the next geographical region to analyse. When the end of a track was reached, then I'd pick the next "untraversed" point in the DB and start again. Then repeat this cycle until every point in the DB had been traversed. Initially the data is in a mySQL database table. I'd added an extra field, in which I'd converted the lat/long pairs to a geo "Point" field. Then used some rudimentary "MBRContains" queries to return points within a certain area (eg 100m2). I was then going to write some java code that explored/traversed the points, with repeated calls to get points within small areas. So my question(s) are: Does PostGIS offers any functions that could assist with the task of identifying related points ? Or are there any better approaches out there to achieve this ? Maybe there are tools/routines to do just this kind of thing, but I don't know the terminology enough to know about it yet :) Thanks for any advice.


Problem with taverna

I want to use taverna to call OGC WFS,and the file returned from WFS is GML, but what I want is SVG,so I have to convert GML to SVG by using taverna.The Method I used ExecuteProcess_gml2svg,which comes form a WSDL,.The process and the GML datasource is as shown ![enter image description here](http://i.stack.imgur.com/pJL3z.png) The problem is when I use the datasource url in the picture (that is the 'random_points' gml datasource is http://rsg.pml.ac.uk/wps/testdata/random_points.gml ),it works well,but when I download the datasource ,and use Appache server in my pc to publish the same gml datasource(the url is like http://chenhongyan.com/random_points.gml ),the process does not work,and the errors are: Processor 'ExecuteProcess_gml2svg_2' - Port 'ProcessOutputs': An error occurred invoking the WSDL service object of type 'exceptions.IOError' has no len() and I just don't know what can I do to solve it.Can you help me?


Record Set (Recordset) for Join in ModelBuilder

I am trying to use a recordset to develop a Table View to use to join to a feature layer in model builder. I can export the table view directly with the user input included. However, when I attempt to use this table view in a join i get no features in the resulting dataset. ![My Model](http://i.stack.imgur.com/KE7kk.png) FYI. The ultimate goal is to allow users to input multiple attributes from a feature layer (Using domain in the field of the recordset) and assign a unique value (i.e cost) for each of those attributes. I would then use the join to filter the table on those results and include the user added value in the resulting feature layer. So I am also open to another approach. Thanks!


Is there a way to force a featureclass to prevent the creation of non-simple shapes when using native geometry storage?

After spending a good chunk of the morning trying to figure out why a feature was displaying a circle when I was setting it as a rectangle polygon via SQL, I found the SDE_GEOMETRY table and some documentation about how ArcGIS stores curves in a side table. So even though I was setting the shape field correctly, ArcGIS was using the original geometry that was stored in the SDE_GEOMETRY side table (for those features that had been created with curves). Is there a recommended way to prevent a this side table from being used and have ArcMap disable any tools that would create the specialized geometry. My current workaround is to delete the row in an after insert/update trigger (seems to work but I haven't really tested it much).


Creating a new time column in PostGIS / CartoDB

I received a CSV file that has year (just two digits, '70'), month (as 'JAN'), day and time (as '0310') as separate columns. Bringing in to CartoDB / PostGIS they all appear as strings. I'd like to convert them to a single column that is a real 'date' type so I can display with Torque. What's the best way to go about doing that?


Projection problem? -- OpenLayers, Drupal, and MapBox

I've installed OpenLayers in Drupal 7 and created a view that pulls geofield data for a "Geography" vocabulary (states, tribes, and U.S. jurisdictions) and shows points on a map that I'm storing on MapBox's server. The map tiles I created in TileMill. The map displays fine, the points appear in the proper locations, but I can't get the map to center properly -- it appears somewhere off the coast of western Africa. I'm guessing that this is a "projection" problem, right?   * TileMill maps are always projected to "Web Mercator" (EPSG:3857)   * So, for the layer that I've added, I've used layer type "XYZ" and set the "Projection" to "EPSG:3857"   * In the "Layers & Styles" tab for the map, I've selected "EPSG:3857" for the "Map Projection"   * I've also set the "Display Projection" to "EPSG:3857"   * Because I'd like the map to be centered over/near California, I've made the "Centerpoint" in the "Center & Bounds" tab equal to "-122, 33"   * I've also tried "Zoom to Layer" in the "Behaviors" tab -- but that doesn't seem to work either. I'm guessing there's something conceptual that I'm missing here -- that there's something about projections that I don't understand, right? Or something else? I'd appreciate any help anyone can offer... I'm so confused!


FME remove polygon holes of under size ### sqft and chop the ends off

Quetion #1   I have a wokspace that takes the latest version of my data and creates buffers and then merges those to one polygon.   When 2 buffer polygons are near each other and close to create a hole. I would like to find a reasonable area size and remove only the holes. If there are small polygons for one reason or another I need to keep them.   ![remove small hole](http://i.stack.imgur.com/tuZtM.png)   Another question not exactly related mentions the donuthole extractor. Do I need to extract these and then delete them?   **EDIT:**   I used the solution with a tester and WAS successfull. The part I was missing was I had to put all of this "behind" my first disolver.   It wasn't an actual donut until the two parts were together. I suppose I could get fancy and put some other locgic in front to make it that way.   The other problem I was having was it is using sq meters not hectares or acres. :)   On to the second problem which I have not had any luck with yet.   ![image3](http://i.stack.imgur.com/gyRWm.png) ![workspace](http://i.stack.imgur.com/BlGYw.png) Question #2   I like the rounded version not the butt of bufferer. However would like to limit the ends (after the merge) to being only a specified distance from the end of the actual line.   ![bufferchop](http://i.stack.imgur.com/D3pEg.png) ![image6](http://i.stack.imgur.com/SbNHX.png)   OK I figured out I need to use the parameter created when I set the AreaCalculator. But I still don't see proper results. I would like to chop off the end flat at maybe 20' offset from the actual end of the line. Can I get a direction or pointer on accomplishing this?


Modelbuilder Parameters - cannot select input files when model dialogue box is open

I created a model in modelbuilder and set the inputs to be parameters. When I open the model and the tool dialogue box shows up, the input files that I wish to select are not shown. These inputs are used for the "Select layer by location". Why is that and how do I fix it? I have all the input data set up in a geodatabase, and I can see it when I navigate to the .gdb in ArcCatalog.


ArcGIS: Rapid digitizing of Point Data

I have several big maps with a lot of points that I have to digitize. All these Points have different Values (depth in Meter, double) and will be stored in 1 Shapefile. Is there a way to speed up my digitization process? e.g. Create Point with my Mouse in right Hand and directly Type the Value into a NumPad in my left Hand? Right now I always have to select the Item in the Attribute Table, where i can edit the Value. (Too much clicks) Thanks!


Collector cable optimizer

I'm hoping to find some pointers for a project that i'm supposed to be working on as part of a team. We're going to create an extension within ArcMap that uses python to process information within a shapefile that the user drags into a toolbox. The program is supposed to determine as final output the best routes for laying collector cables. The model will need to account for, among other things, road allowances, properties signed onto the project, the maximum number of Wind Turbines per electrical circuit, environmental features, cost for running x number of collector cable circuits overhead, cost for running x number of collector cable circuits underground etc It's a pain because i've got to learn a new language from scratch. Anyway, i understand that such models exist already and i was wondering if someone could shoot me off in the right direction by pointing one or more out and possibly any other useful resources as well. Thanks


QGis points distance from shore in metres

I have a grid of depth points in the Irish Sea that i'm trying to add other data fields to, in this case distance from shore. I've tried both of these suggestions to no avail, but am having some success with this one. My problem is that the resulting raster ![resulting raster](http://i.stack.imgur.com/tpNgk.png) has distances as decimal degrees rather than metres. I've changed my project properties > canvas units to metres, re-done the rasterise, re-done the proximity, but it's still decimal degrees. Any ideas how I can get the results in metres?


Null Values After "toreal" Command

I am very new so help is greatly appreciated! I have joined a simple CSV file to my project and am attempting to change a string value to a real value. In editing mode, I bring up the field calculator and create a new field from the first field. I try to change the string to real by doing the command "toreal(nameofmyfield)". When I create the new field though, I get NULL for every value, even though I see the normal numbers in the column I am trying to change. Any idea what's going on?


Calculate attributes to table based on other attributes

I am going through thousands and thousands of road lines and giving them a value based on QC with imagery. So for instance for the add_data field I give it one of these four numbers 0 = not examined 1 = include 2 = consider excluding 3 = exclude 4 = DC But for each of these I also need to fill out comments and type of road data. So I was wondering if there is a way to calculate the comment and other fields based on the add_value field? So maybe it might look like this for the add data: "1" = "include" comment = looks good RD_Typ = Major Road.   "2" = "include" comment = looks good RD_Typ = Dirt Road.   "3" = "include" comment = Joined data to other data RD_Typ = Major Road.   "4" = "include" comment = Joined data to other data = Dirt Road. ..... and so on you get the point. So the question is can I calculate the other fields based on the number I calculate in the add data field. I know I could do this with a join and then calc field but that just does not seem to be as much fun. Thanks for your help.


How do I limit update/insert/delete capabilities for a WMS layer in a QGIS project?

In QGIS 2.0 I have a layer based on a query to a PostgreSQL/PostGIS database and would like to serve this as a WFS layer. Users of the layer should be able to select subsets of the data based on properties such as dates or unique IDs, but should not be able to update, insert, or delete records. The Creating a WMS/WFS from a QGIS project section of the documentation simply says, (emphasis mine) > **In the WFS capabilities area you can** select the layers that you want to > provide as WFS, and **specify if they will allow the update, insert and > delete operations**. How, exactly, do I specify this? What is the correct way to deny, and to allow, these actions? * * * Here's what I've done so far. I added the layer I want to publish to a new QGIS project. Then in the Project Properties window, under OWS server, I checked the box in the Published column for the layer. Now I don't know what to put in the free-form text input boxes in the Update, Insert, and Delete columns. I've tried putting values of "0", "NO", and "FALSE" in here, as well as more facetious entries such as "Sure" and "Why not?". No matter what I enter though, when I click Apply or OK, QGIS gives no error or feedback of any sort. When I return to the project properties later, the text boxes are all blank as they were initially. Oh, and in case it makes a difference, I'm running this under Windows (currently 7, but it will be deployed on Server 2008 or 2012 eventually).


How to create Esri Flex Geometry from WKT

I am working with Esri flex 2.3 API, I am also using Abode BlazeDS for server side binding. I am sending a geometry as WKT from server into my flex client and I am looking for a way to convert the string into a flex geometry so that i can make a Graphic object and add to Feature Layer! or do I have to write a parser myself to accomplish this?


Symbology or Styles

Anyone know where to find me some more symbology or styles that I can use for my ArcMap 9.3.1 ? I have lots of mxd files on my computer and how do I apply all of my symbology or styles turn on ? Every time I open my mxd files I have to turn on each in the Style Reference dialog box because they are not enable check.


Looking for an annotated vectorial map of the world

I'm looking for an annotated vectorial (PS/SVG/...) map of the world (+countries). Wikipedia contains many vectorial maps (e.g: http://commons.wikimedia.org/wiki/Category:Blank_SVG_maps_of_the_world )but the paths are defined like that:               <path fill="#ccc" d="M98.171,199.174c-0.505-0.076-0.938-0.322-0...      ... and I'd like to find an **identifier** so I can colorize a given country (using XSLT or another method....). Something like that:               <path id="JAPAN" fill="#ccc" d="M98.171,199.174c-0.505-0.076-0.938-0.322-0...      thanks


Where can I find the various Esri tools to design a geodatabase?

There is a confusing multitude of tools written by Esri and others (some supported, some unsupported) for designing and diagramming a geodatabase. Here is an attempt at creating a definitive list of where to find them. Related to What are convenient uml tools to create geodatabase models?


Is it possible to create color tables for floating point raster maps in GDAL?

I've got a script that dumps a floating point raster that I'd like to automatically color. Currently I load this raster in Arc or QGIS and manually create a color table every time I re-run my model. I'd like my model to automatically create a nice color table, but the only functionality I can find is to create a color table for integer only rasters. My approach so far is to output 2 rasters, one with floating point values rendered in grayscale and a second integer raster that has the colors on it. I can flip between the two to see values, but it's far from optimal. Is there a technique to procedurally create floating point color tables in GDAL that automatically load with QGIS or Arc?


ArcGIS 10: AddIn Extension loading toolbox issue

I have an extension class in an ArcMap addin project of which loads a custom toolbox into ArcToolbox when the user opens ArcMap or opens a previously created map. My toolbox gets added to ArcToolbox fine except when the user creates a new map from File > New. My toolbox just disappears (no messages). I'm running AG10 SP2 but it doesn't work in AG10 base either. Anyone else experience this? Solutions?               protected override void OnStartUp()     {        AddToolboxes();        #region Wire events        ArcMap.Events.NewDocument += delegate() { ArcMap_NewDocument(); };        #endregion     }     void ArcMap_NewDocument()     {        AddToolboxes();     }     void AddToolboxes()         {                  // Load toolbox into ArcToolbox             IArcToolboxExtension arcTbxExt;             try             {                 arcTbxExt = (IArcToolboxExtension)ArcMap.Application.FindExtensionByName("ESRI ArcToolbox");                 if (arcTbxExt == null) { throw new Exception("ArcToolbox extension could not be found."); }             }             catch (Exception ex)             {                 System.Windows.Forms.MessageBox.Show("Error finding ArcToolbox extension." + Environment.NewLine + " - " + ex.Message);                 return;             }                  try             {                 IGPToolbox2 toolbox = Maersk.Functions.Toolbox.OpenToolbox("C:\\Toolboxes\\MyToolbox.tbx");                 arcTbxExt.ArcToolbox.AddToolbox((IGPToolbox)toolbox);                 arcTbxExt.ArcToolbox.Refresh();                  }             catch (Exception ex)             {                 System.Windows.Forms.MessageBox.Show("Error launching  toolbox." + Environment.NewLine + " - " + ex.Message);                 return;             }              }


Moving a Staging ArcGIS Server to the DMZ?

I am trying to move my EDN Arcgis Server with all the mapservices and caches to the DMZ .. The move is from our intranet to the DMZ to provide access to mapservices to outside world using the internet port 80 or 8080.. Everything is fine when we try accessing applications and mapservices from within our network to this box in the DMZ!! When it comes to external access it cannot access any ?? **I have tried a similar setup using windows xp box in the DMZ! Both internal and external machines can access and works well.. This should account for Network configuration correctness, hopefully..** _I am wondering if it has to do with ArcGIS Server on Windows 7 ?? Am I missing something wrt settings with ArcGIS Server or IIS ??_


Icon missing when adding custom combobox command to Toolbar

I have a custom command:               class ATL_NO_VTABLE CSliceSelector :         public CComObjectRootEx<CComSingleThreadModel>,         public CComCoClass<CSliceSelector, &CLSID_SliceSelector>,         public ISliceSelector,         public ICommand,         public IComboBox     {        ...     }      I also have a custom Toolbar where I try to add my custom command:               class ATL_NO_VTABLE CSeismicToolbar :         public CComObjectRootEx<CComSingleThreadModel>,         public CComCoClass<CSeismicToolbar, &CLSID_SeismicToolbar>,         public ISeismicToolbar,         public IToolBarDef     {        ...        STDMETHOD(GetItemInfo)(long pos, IItemDef * itemDef) {             if (itemDef == NULL)                 return E_POINTER;                  if (pos == 0) { // CSliceSelector                  itemDef->put_Group (VARIANT_FALSE);                 HRESULT hr = itemDef->put_ID(CComBSTR(L"{96A2B223-BF0D-410C-B22B-892FED3D9A15}"));                 if (FAILED (hr)) {                     return hr;                 }             }            }     }      However when I start ArcScene; I see a red circle icon and the text "Missing": ![enter image description here](http://i.stack.imgur.com/N7W7K.png), instead of my command :-(. However if I add it manually from Customize toolbars; it works fine. I have tried both returning a null pointer and a valid bitmap:               CSliceSelector()     {          m_hBitmap = ::LoadBitmap(_AtlBaseModule.GetModuleInstance(),MAKEINTRESOURCE(IDB_SLICE));     }     STDMETHODIMP CSliceSelector::get_Bitmap(OLE_HANDLE * Bitmap)     {          if (Bitmap == NULL)             return E_POINTER;        *Bitmap = (OLE_HANDLE) m_hBitmap;          return S_OK;     }      I have other commands that are not combobox; and they work fine.


Load partial shapefile into Postgis / GeoDjango project

I have a shapefile with Canadian postal codes, but am only looking to load a small subset of the data. I can load the entire data file and use SQL or Django queries to prune the data, but the load process takes about 2 hours on the slower machines I'm using. As the data I'm actually after is about 10% of the dataset, this isn't a very efficient process. I'm following the instructions in the Geodjango tutorial, specifically the following code:               from django.contrib.gis.utils import LayerMapping     from geoapp.models import TestGeo     mapping = {'name' : 'str', # The 'name' model field maps to the 'str' layer field.                    'poly' : 'POLYGON', # For geometry fields use OGC name.                    } # The mapping is a dictionary     lm = LayerMapping(TestGeo, 'test_poly.shp', mapping)     lm.save(verbose=True) # Save the layermap, imports the data.      Is there a way to only only import data with a particular name, as in the example above? I'm limited to the Linux / OS X command line, so wouldn't be able to utilize any GUI tools.


QGIS proximity raster with extent beyond original coverage?

I need a raster coverage of distances from points outside the geographic region of interest, and am having trouble doing this with QGIS. Questions like this one (Calculating distance to points in QGIS) seem close to what I want, but I'm getting hung up on the first step. The raster that I create with the QGIS "Rasterize" tool has an extent that does not overlap with the region of interest (because all of the points I want the distance from are outside of that region). I don't see any option for changing the extent of the generated raster. I am then unable to generate a Proximity raster with the extent I need. How can I get the correct extent? In addition, when I run "Proximity", I consistently get an error that my raster cells are not square and so distance calculations may not be accurate. I don't know how to make square cells because the tool asks for how many cells I want in the x and y directions, not the x and y dimensions of each cell. Is there a way of fixing this? Ultimately, I'm a bit uncomfortable with the necessity of converting point data to raster in order to calculate proximity, and wonder whether there is an open source tool that can do what I want without the "Rasterize" step...


Set symbols for a layer in ArcMap with two classes (greater than zero and less than zero)

I need to symbolize features in a feature class with two colors thus defining what features have positive (>0) and negative (<0) values. The only two ways I found this to be possible are: To use a static field with the values calculated beforehand OR Use definition queries using two identical layers in the TOC. Since it is more time consuming to update labels for two layers, I have approached this problem in another way by using Graduated colors symbol setting with two classes (Manual classification). I have added two extra features in the feature classes with really big and really small values (to cover possible values that can come later on while the feature class will be edited) and then edited the intervals in the Range columns in the Symbology tab of the Layer Properties. Afterwards, those two features are deleted. Now features' values will fall within these two ranges and will be symbolized properly. ![enter image description here](http://i.stack.imgur.com/qkBMb.png) Considering the need to update multiple documents with these settings and having different datasets as layers, are there any other approaches (without using Python to automate the hack described above)? All other GIS.SE posts I've found so far refer to the methods outlined above.


Arcmap point symbol/marker size, How it works

I am trying to figure out how the size of markers in arcgis works. I understand that a square marker set to 1 pt is 1 pt for specific scale. but when I import emf or font, How does the change size works then? if like the symbol is rectangular or a picture of a bus? if we take the bus, I set 20 pt size does it take 20 point of the smallest part of the picture or what? please help! :) and also, how can I set a specific size to be correct for illustrator export? the size gets different, is there any script for arcmap or illustrator to set correct size of symbols?


How do I format a text box in MapInfo

I want to import a table (created and formatted in MS word?) to appear as text in the MapInfo layout window. Obviously this isn't individual data records which appear in the MapInfo browser window. So would I need to format the MS Word table, save the file as a JPG/png format then import into MapInfo as a raster file? Is there another way in MapInfo of doing this - just type the table contents as text into the layout window then draw a rectangle around the text and make this box transparent but with a border?


How to load interlis1 files directly into QGIS 2.0.1 Dufour 64Bit (W7Pro)

with earlier version of QGIS (1.8) it was possible to load data in interlis 1 format directly. If I try now it says something like unknown format. Is this because I'm using the 64Bit version of the installer? A possible workaround could be using the examples at http://giswiki.hsr.ch/HowTo_OGR2OGR "interlis to shp". It would be very nice to be able to open Interlis Files directly again. Thank you for any hints. Regards, Adrian


Why does one part of GML -> KML conversion not work on Google Maps?

Brand new to all of this, not even sure this is the right forum to ask in... I have a client who wants to visualise cancer incidence and mortality rates by local/unitary authority across Great Britain. I've used Ordnance Survey data to get the boundaries of the authorities in GML format and saved all of them as one large GML file. I used QGIS to convert the data to KML format. I have the data in Fusion Tables now and it all looks great except for one blatant gap in the coverage. When I look at the KML data for the missing local authority, it hasn't imported at all (ie the cell is empty). I've tried re-uploading the geometry (and re-downloading and re-converting) and Fusion Tables is definitely rejecting it as KML. I've read somewhere that Fusion Tables can mess around with your geometry somewhat, so I uploaded the KML file and pointed Google Maps directly at it and that didn't work either (I don't know how good a test that is). I don't know that I originally created a valid GML file (I don't know the standard) but they seemed to import ok to QGIS and the final result looks great apart from this one place. The area _is_ visible in QGIS when I load the GML file, so I don't know if QGIS is producing invalid KML, or mayble Google Maps/Fusion Tables have limitations? **Edit:** I've downloaded the OS Open data pointed to by nhopton, loaded the appropriate layer into QGIS ( _Layer_ -> _Add vector layer_ ), saved it as KML ( _Layer_ -> _Save as_ ) and uploaded it to Fusion Tables. (Note that I didn't explicitly do anything with the CRS). Generally it's worked well (and importantly the particular area that had failed with my first method worked with this, so I have the missing geometry!). However, there's a number of missing areas again. And in this case, there's much less that _I've_ done to the data myself, so it's not so easy to assume I've just messed up the data. Any ideas why my final results look like this? I'd really appreciate some guidance from anyone who knows this area a bit better:) TIA


How to find the elevation differences between adjacent pixels in ArcGIS?

Is there a way to determine what the difference in elevation is between two adjacent pixels? I am using an SRTM DEM to map faults. I am using the vertical error (>16m) as a criterion for the difference in elevation to accurately map these faults. There must be a more efficient way to do this than by using the i-tool. Cheers


SQL expression for finding strings including letters in ArcGIS 9.3.1

I'm working on ArcGIS 9.3.1. I have a shapefile with a string field containing point names like these:               10     12     12/13     t/9     6577y     1234     12t9     124y67     u9      and so on (and yes, I know it looks like a colossal mess. It's not my design). Using `WHERE MyColumn LIKE '%Search String%'` I can find (using select by attribute) the ones containing a specific letter or a backslash (/). MY QUESTION IS THIS: How can I select all values containing a letter (in this example: t/9, 6577y, 12t9, 124y67, u9), not caring what letter it is or it's position? (is there a way to specify a position?)


QGIS converting to ASCII question

I am attempting to do some analysis on a raster file and I need to convert it to ASCII. In QGIS, I used [Raster -> Conversion -> Translate] with the asc file format, and I suspect it is simply copying the file without converting anything. Are there any other suggestions on how to do the conversion?


points on google layer using openlayers

HI, I want to load some points, whose lat-long values are saved in a .txt file over google maps. The following is the code which I am using for the same. The google map layer is succesfully loaded if I keep lines 1-10 & 13-17. For loading the data points I am including 11-12 lines but I am not getting any output at all. I request you to please find the error and suggest how to do the same. "data.txt" contains the lat-long values with tab delimited present in my server.               1.     2.      3.<title>OpenLayers Example</title>      4.      5.      6.      7.

      8.      16.      17.       When I am replacing above 11-12 lines with               11.var pois = new OpenLayers.Layer.Text( "My Points",                      { location:"./data.txt",                        projection: map.displayProjection                      });      12.map.addLayer(pois);      Even then I am not able to get the output. I require some needful help from any one of you


Does anyone know if there are any sample questions available for the ESRI Web Application Developer Certification exam? Thanks, Seth _**_*** _UPDATE_ ** _*_** ESRI has released 3 sample questions for this test: http://downloads2.esri.com/campus/downloads/certification/samplequestions/EWDA10_Sample_Questions.pdf


WMS layer disabled when i grouped

i am using geoserver 2.4.0 for my GIS based web application.My Wms layer disabled when i grouped more than 3 layer in my application.Previously it worked fine.This problem occurs in real time only,in geoserver this problem does not exist.i don't know! any updation happened in geoserver? How can i fix it?kindly provide me solution for this problem?


Why is the geonames not loading?

I want to import a geonames file. I downloaded the geonames .txt file from http://download.geonames.org/export/dump/ I used: `v.in.geonames input=PH.txt output=ph` but it displayed this error: `ERROR: File 'C:UsersRhenDesktopPHPH.txt' not found` I don't know why it cannot be imported..


create temporary point geometry feature from coordinates in OGR with Python

I have two shapefiles containing two types of geometry features: one contains polygons and the other linestrings. Because I'd like to use the Within method to check whether the end(or start) of a linestring falls within a polygon located in the other shapefile, I need to access the first and last points (as geometry features). At the moment I iterate through the linestring feature and can identify the points using geometry.GetPoint(x) which returns the x and y coordinates of that point. But when I'd like to use within to check whether a given point falls in the polygon object I can't because as Paolo explained in this question I need the geometry reference (using GetGeometryRef())of a point which can be retrieved like:               ...     pointshp = driver.Open('test.shp')     pointlyr = pointshp.getLayer(0)     point = pointlyr.GetNextFeature()     point_geom = point.GetGeometryRef()     point_geom.Within(polygon_geom)      but I'm actually opening a shp that has linestring features. So after skimming through Utah State University GIS slides, I came up with this:               def within(pointcoord):     driver = ogr.GetDriverByName('test.shp')     ds = driver.CreateDataSource('test.shp')     layer = ds.CreateLayer('test',geom_type=ogr.wkbPoint)     fieldDefn = ogr.FieldDefn('id',ogr.OFTInteger)     layer.CreateFiled(fieldDefn)     featureDefn = layer.GetLayerDefn()     feature = ogr.feature(featureDefn)          point = ogr.Geometry(ogr.wkbPoint)     point.AddPoint(pointcoord[0],pointcoord[1])     feature.SetGeometry(point)     feature.SetField('id',1)     layer.CreateFeature(feature)     pointlyr = layer.GetLayer(0)     point = pointlyr.GetNextFeature()     point_geom = point.GetGeometryRef()     for x in xrange(len(polygons))             point_geom.Within(Network.polygons[x].GetGeometryRef())      Which looks really inefficient to me. So just for the sake of using a method(within) I need to create not only a layer, but a shapefile, when all I really need is a simple point feature from which I can get the geometry reference to compare it with the one of the polygons. The question is, how to create a point feature (temporary one, to be destroyed) from x and y coordinates? Or is there a simpler way of doing that than the example above? On another hand, is there any other way of getting the geometry reference of the start or end point (as a geometry feature) of a linestring feature?


How to select a raster file that is under a polygon in a vector file?

I am using: v.db.select map=PHL_adm1 where=NAME_1='Bulacan' My vector map is over the raster map. I want to select the part of the raster map under the polygon(NAME_1='Bulacan') of a vector. Is it possible to do?


calculate kernel density for polygons in ArcGIS

I have a population data set at the census block level. My goal is to use the data set to identify some population centers in the region based on population density. The simplest way to do it is to calculate the density by census blocks and use an arbitrary density cut-off value to determine those centers. One shortcoming of this method, however, is that it relies on census blocks, which have irregular size and shape. The other approach I am thinking of is to covert census blocks into points based on their centroids, and calculate kernel density based on these points. But I am not sure if this makes sense. In other words, does this approach help me overcome the problems that are associated with the first approach? Any advice is appreciated. Also, if it makes sense to do that, how the output raster size should be determined?


Water Quality Data Acquisition Sources

I was wondering what sources are available, or if anyone has water quality data for the Great Lakes and specifically Lake Ontario. We are interested in mapping Phosphorus, Organic Matter, Nitrogen and Dissolved Oxygen concentrations in Lake Ontario as they relate to algal blooms and hypoxia. Thank You


PostGIS: convert point data to raster

I have a table of points and I'm trying to visualize their distributions. Each point has a POINT geometry (WGS84) and some numeric attributes: a, b, c.               | id | geom       | a  | b | c |     |-----------------|----|---|---|     | 1  | POINT(...  | 12 | 9 | 2 |     | etc ...                      |      The first thing I tried was to project to an equal area projection (EPSG: 3175), use `ST_SnapToGrid` and then group on the resulting point geometry.               SELECT        avg(a), avg(b), avg(c),        ST_Transform(ST_SnapToGrid(ST_Transform(geom, 3175), 50), 4326) AS bin     FROM points     GROUP BY bin      This groups the points into 50x50 meter bins with their aggregate values which I can import into qgis to play with. I optimized this by adding another column to the points table called `bin_id` which was a foreign key to a `bins` table:               | id | geom      |     |----|-----------|     | 1  | POINT(... |      Then I could write queries like this:               SELECT avg_a, avg_b, avg_c, geom     FROM (       SELECT          avg(a) AS avg_a,         avg(b) AS avg_b,         avg(c) AS avg_c,         bin_id       FROM points       GROUP BY bin_id     ) AS tmp     JOIN bins ON tmp.cell_id = bins.id      This all works fine but i want to use Rasters instead. All the examples I've seen show how to import existing rasters or shapefiles into PostGIS. But I need to know how to generate one from data already in PostGIS. Ideally I want a raster with 3 bands: `avg_a`, `avg_b` & `avg_c` Can someone point me in the right direction? EDIT 1: I'm not opposed to doing it outside of the database.   EDIT 2: Seems like I can use the GDAL API to do it.


How to set custom break points based on formula?

I work on projects that each require the creation of about 40 maps. In all of these maps, I use Graduated Color symbology to represent concentrations of a different population in area Census tracts. The break points I use are not equal interval--basically the first break point is the mean, the second is the mean+10 (the "disproportionate share threshold"), and the third and fourth break points are equal intervals between the disproportionate share threshold and the maximum value in the dataset. The maximum value is the last break point. I have a workaround that I'm using. Well, I've got the "around" part--the "work" part is not quite there: I use SPSS to collapse population per tract into either the minimum value, the mean, the disproportionate share threshold, then the two break points between that threshold and the maximum value, and the maximum value itself. Once I've done that, I join the table to a Census tract shapefile, and open up symbology. When I set the classification to "Natural Jenks", the break points are ALMOST exactly what I want, except that the minimum value is in a class by itself, and one of the upper break points are omitted. So, I have to adjust it by hand anyway. Is there some way to have the break points "skip" minimum value and snap to all of the values above that? Ideally, each of the five breaks should include a range between:   1. the lowest value and the mean,   2. the mean and the disproportionate share threshold (DSP),   3. the disproportionate share threshold and ((maxvalue-DSP)/3+DSP,   4. that last one and (2((maxvalue-DSP)/3)+DSP,   5. that last one and the maximum value. Thanks in advance for your help--I hope this was clear. I'm pretty new to ArcGIS, and just getting to the point where I "know what I don't know".


how to create a map with oneway flows using FlowMapper for QGIS

my plan is to create a map that shows the top 50 (by amount of passengers) international flights departing the United States (e.g. New York -> London, Honolulu -> Tokyo...). I want to use the plugin FlowMapper for QGIS to display those flights and have the width of the flows correspond with the amount of passengers. I have read the tutorial and had a look at the examples given (all in the download file of FlowMapper). The problem is, in those examples FlowMapper is only used for back and forth flows. Thus, the data matrix in the examples looks like this, each row and column have the same: ![enter image description here](http://i.stack.imgur.com/NxBTQ.jpg) But I want to use one way flows, so row and column have different geographical locations, like this: ![enter image description here](http://i.stack.imgur.com/khsth.jpg) So my question is, how does the flow matrix- and coordinates-file have to look like in order to be correctly displayed by the FlowMapper plugin? This is my complete speadsheet that I wanna use as basis for my map: http://www.sendspace.com/file/ckaejd


ArcPy Cursor iteration update

I'm trying to get the cursor to iterate through rows and fields looking for blank or null values and change them to `999999` but it's not changing anything.               fc = "C:\Users\\bbrock\Documents\ArcGIS\Ports.shp"          # Create a search cursor      #     rows = arcpy.SearchCursor(fc)           # Create a list of string fields     fields = arcpy.ListFields(fc, "", "String")          for row in rows:         for field in fields:             if field.type != "Geometry":                 if row.getValue(field.name) == '':                     row.setValue(field.name, '999999')                                   print "%s: Value = %s" % (field.name, row.getValue(field.name))                      if row.isNull(field.name):                     row.setValue(field.name, '999999')                                   print "%s: Value = %s" % (field.name, row.getValue(field.name))


QGIS/MATLAB add geo-coordinates flag

I've created some polygons in QGIS with which I'd like to create some thematic maps. Unfortunately, if I want to import the shapefile in MATLAB with: shaperead('nameoffile, 'UseGeoCoords') i get an error which states: Expected geo-coordinates flag (T/F) to follow parameter name 'UseGeoCoords'. Below you'll find a screenshot from MATLAB: ![MATLAB Screenshot](http://imgur.com/uNQCTKu.png) As far as I see there should be an option in QGIS, with which I can say that programs can use GeoCoords and not only X,Y coordinates. Has anyone an answer for this specific problem? Thanks and best regards!


Tools for planning a building

I'm planning to build a cabin, and I'd like to combine all the GIS data available from my home country with my own measurements for placing and the cabin layout. That is, I need a tool which can:   * save data in layers,   * import third-party GIS data for elevation, infrastructure, etc., and   * import a PDF, SVG or other vector graphics format to get the cabin layout. Is there any software useful for this available for Linux (at a pinch Windows XP), preferably for free?


Publish MXD to server in ArcObjects C#

I have an MXD that I would like to publish to an ArcGIS server using C# ArcObjects (preferably as an MSD). I have been unable to find anything regarding this. Essentially I would like to implement the following ArcPy code in ArcObjects:               arcpy.mapping.ConvertToMSD()     arcpy.mapping.PublishMSDToServer()      Any help would be appreciated! Thanks!