title: Why does copying ArcSDE datasets and features in Python give ERROR 000732 & 000735?

text: Go easy, I'm new to python I keep getting a kick back that says this: ERROR 000732: Input Features: Dataset 'Database Connections\Wyoming.sde\Supporting_MT_DBO_Geology' does not exist or is not supported Failed to execute (FeatureClassToGeodatabase). _**_** _*_ ** _*_ ** _*_** _UPDATED_ ** _*_ ** _*_ ** _*_ ** _*_ *** _****_** _*_ ***** _HERE'S THE ENTIRE SCRIPT_ ** _*_ ** _*_ ** _*_** >     # Import modules     import sys, string, os, datetime, arcpy          # Sets the start time     nowstart = datetime.datetime.now()          # Set location of SDE instance - in this case a special SDE connection just for this script     SDEname = r"C:\Users\jjones\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Wyoming.sde"     SDEnameShort = r"'Database Connections\Wyoming.sde"               # Set up the log file     LogLoc = r"P:\Database_Backups\SDE\BackupLog_"     LgFile = LogLoc + nowstart.strftime("%Y_%m_%d_%H%M") + ".log"     logfile_out = open(LgFile, 'w')     # Write the start time to the log file     logfile_out.write("Backup process started on " + nowstart.strftime("%Y_%m_%d %H:%M") + "\n")          # Deletes temp timestamp     del nowstart          # Add the Conversion and Data Management Toolboxes to the geoprocessor               #Sets work environment to SDEname     arcpy.env.workspace = SDEname          #Sets environment to overwrite     arcpy.env.overwriteOutput = True          # Resets the start time     nowstart = datetime.datetime.now()          # Creates a variable for naming the backup GDB     YearMonthDayTime = nowstart.strftime("%Y_%m_%d_%H%M")          # Variables and parameters for setting up the backup GDB     outgdb = "PetroPAM_Backup_" + YearMonthDayTime + ".gdb"     outpath = r"P:\Database_Backups"     output = outpath + os.sep + outgdb     print output               #Create the backup GDB     arcpy.CreateFileGDB_management(outpath, outgdb)     print "Created BackupGDB"                         # Find and export feature classes in datasets to a backup GDB     logfile_out.write("\nCreating Feature Dataset List...\n")     fc = arcpy.ListFeatureClasses()          dsList = arcpy.ListDatasets("*", "Feature")     print dsList     print "The Datasets"     for dsGEO in dsList:         outDS = arcpy.ValidateTableName(dsGEO)         arcpy.CreateFeatureDataset_management(output, outDS)         fcListGEO = arcpy.ListFeatureClasses("*","",dsGEO)         print fcListGEO         print "Feature Classes in Dataset"         for fcGEO in fcListGEO:             print fcGEO             print SDEnameShort + os.sep + outDS + os.sep + fcGEO             arcpy.FeatureClassToGeodatabase_conversion(SDEnameShort + os.sep + outDS + os.sep + fcGEO, output + outDS)      Any ideas? I'm using ArcMap 10.1


title: Extracting flow accumulation to stream shapefile

text: I would like to create a shapefile that displays the flow accumulation for a given point in a stream channel. Is this possible? My client wants a map that shows the accumulated flow/catchment area for each point in a stream channel--as catchment area increases the color of the stream shapefile transitions from red to blue, kind of like a heat map but just for a stream channel. I am able to approximate this using a flow accumulation raster in ArcGIS 10.1, but the color (and by proxy the flow accumulation data) are lost in the final map product due to the small size of the raster cells relative to the other features. I'm wondering if there is a way to extract the data from the flow accumulation raster to a stream shapefile? I expect there would have to be some sort of interpolation necessary but I am at a loss on how to proceed. My potential work around is to query the flow accumulation raster by a series of flow accumulation thresholds, convert these to polylines and layer them on top of each other in the map product. However, there may be up to 20 different polyline shapefiles, so it may be extremely tedious. Update with workaround solution: I used Model Builder to put together an iterative model that queried the flow accumulation raster for cells that fell within a specific range of values using conditional inequalities and common tools to convert stream rasters to stream shapefile. The ultimate output was a set of stream shapefiles that corresponded to specific ranges of stream discharge. These shapefiles, though, did not have any discharge data in their attribute tables. The lack of discharge data was a problem as my goal was to create a single shapefile where the symbology/color would represent the increase in accumulation/streamflow. The workaround was that I modified my model to appended a new field to each attribute table to describe discharge, and the discharge value was linked to the upper limit of the interval from the earlier conditional statement. The resulting shapefiles were then merged into a single shapefile, and the discharge value was used to drive the symbology.


I am new to QGIS. When I use OSM or Google Maps as a basemap via Open Layer Plugin, the map should be projected as EPSG: 900913. However, the Coordinates shown in the status bar are useless for me, what I need to see is lat/lon or UTM. Is it possible to have another CRS in the status bar than used for the map? I did find the Coordinate Capture Plugin, wich can be used to follow the mouse cursor in another CRS, but this is not handy for me and conflicts with the use of other tools. Thanks in advance, Redoute


title: Mapserver openlayer postgis spatial data are in postgis but label I want to display from webservice

text: My database of spatial data wkb is in postgis and my textual data is in sqlserver now is it possible to show labels through webservice. I have common id in both the database through which I can join textual data and spatial data Any help would be highly appreciated


title: QGIS Connection PostGIS and shapefile

text: I'm working on QGIS with a database created in PostGIS and a shapefile. Inside the database, I have geometry type data about location of some points, instead on the shapefile there some selected areas. With these 2 kind of data, I would like to know if there are some points inside or close to the shapefile area. For example: is there a point inside a selected area? Can I select a part of the total area to see points also near to it? Is there a query tool to use? How can I connect these data created in different layers? Thank you so much


title: How to reposition a polygon Shapefile using ground control points and v.rectify?

text: I'm new to Grass. I would like to reposition a polygon shapefile using ground control points and v.rectify seemed just the perfect tool to do that. The raster from which the shapefile comes from was georeferenced using wrong control points for the KKJ1 projected coordinate system (EPSG:2391), now I obtained the correct ones and I am trying to find a way to reposition the vector without having to digitize everything again. After running the command, the vector does move to the desired position, yet the original conformation is not retained and in particular some of the areas "lose" the centroid and become just empty boundaries. Moreover I noticed that the areas whose centroids had not gone missing in the process, presented wrong attributes when inspected. I had a look at the output and 3 errors show up (my output is in Italian so the translation might not be exactly the same): WARNING: Number of incorrect boundaries: 1752   WARNING: Number of centroids out of the area: 139   WARNING: Number of duplicated centroids: 25   Any suggestion on how to solve this would be much appreciated.


title: Error when dividing a field with a number

text: Using arcpy CalculateField_management i want to divide a field with a spacific number:               arcpy.CalculateField_management(EZG_shp, "IDX_AT_PK", "(!EURO_SUM!/!Anz_Pers!)/1668,4","PYTHON")      This works in ArcGIS Desktop 10.2.1 fine, but the result in arcpy is an error (No NULL-Values a allowed)... Any suggestions what I am doing wrong?


title: Converting NAD27 geographic to grid

text: I have a series of lat/long coordinates in NAD27 that i need to convert to NAD27 zone 17n eastings/northings. I have ARCGIS 9.3 but no direct and reliable means for the transformation. Any ideas?


title: Unable to install QGIS 1.8 on Ubuntu

text: Can someone please walk me through installing qgis on Ubuntu 13.04. I keep getting 1.7 no matter what I do and I would like 1.8 if possible. I am a new Linux user. This was addressed in a previous post but it still did not fix my problem. I was still getting the same errors even when I followed this post: Cannot fix QGIS after upgrade to Ubuntu 13.04 Raring * * * I was finally able to figure it out: On the QGIS website I installed this http://hub.qgis.org/projects/quantum- gis/wiki/Download#262-With-updated-dependencies section 2.6.2 I actually had to download these packages from the site:               sudo apt-get install python-software-properties     sudo add-apt-repository ppa:ubuntugis/ubuntugis-unstable     sudo apt-get update sudo apt-get install qgis


title: No fetch Python plugins in qgis 1.8 + ubuntu 13.04

text: I tried to install qgis 1.8 to ubuntu 13.04 but when installation finished I couldn't find "Fetch python plugin" button in plugin menu. I didn't find any soloution for this problem. Any idea what's the problem?


title: Extracting info on polygons from underlaying raster layers

text: In ArcGIS 10.1, I have two layers: 1) raster of Topographic Wetness Index 2) vector layer of polygons representing soil types I am wondering how to accomplish what I am after, that is to get information for polygons from the underling raster. In other words, I would like to extract the mean TWI for each polygon in order to calculate what is (on average) the TWI for each soil type.


title: Couldn't load GLMapCore: findLibrary returned null

text: Im getting this Error when loading an ESRI map inside the Android Environment. This didn't happen with ESRI's .1 Alpha SDK, but it does happen with their 1.0.0 SDK. I know it has something to do with the OpenGL lib, but can't seem to find a solution. I am running the Motorola Xoom Tablet with the latest Android environment. Thanks for your help all! Also Its using Fragments, so that might have something to do with it. Here is my stack dump. > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): FATAL EXCEPTION: main 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): java.lang.RuntimeException: Unable > to start activity ComponentInfo{.map.EsriMapWidget}: > android.view.InflateException: Binary XML file line #13: Error inflating > class com.esri.android.map.MapView 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): at > android.app.ActivityThread.performLaunchActivity(ActivityThread.java:1818) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.app.ActivityThread.handleLaunchActivity(ActivityThread.java:1834) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.app.ActivityThread.access$500(ActivityThread.java:122) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.app.ActivityThread$H.handleMessage(ActivityThread.java:1027) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.os.Handler.dispatchMessage(Handler.java:99) 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): at android.os.Looper.loop(Looper.java:132) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.app.ActivityThread.main(ActivityThread.java:4126) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > java.lang.reflect.Method.invokeNative(Native Method) 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): at > java.lang.reflect.Method.invoke(Method.java:491) 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): at > com.android.internal.os.ZygoteInit$MethodAndArgsCaller.run(ZygoteInit.java:844) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > com.android.internal.os.ZygoteInit.main(ZygoteInit.java:602) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > dalvik.system.NativeStart.main(Native Method) 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): Caused by: android.view.InflateException: Binary > XML file line #13: Error inflating class com.esri.android.map.MapView 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.createView(LayoutInflater.java:596) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.createViewFromTag(LayoutInflater.java:671) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.rInflate(LayoutInflater.java:724) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.inflate(LayoutInflater.java:479) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.inflate(LayoutInflater.java:391) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > .map.EsriMapFragment.onCreateView(EsriMapFragment.java:21) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.support.v4.app.FragmentManagerImpl.moveToState(FragmentManager.java:714) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.support.v4.app.FragmentManagerImpl.moveToState(FragmentManager.java:912) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.support.v4.app.FragmentManagerImpl.addFragment(FragmentManager.java:990) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.support.v4.app.FragmentActivity.onCreateView(FragmentActivity.java:257) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.createViewFromTag(LayoutInflater.java:660) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.rInflate(LayoutInflater.java:724) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.rInflate(LayoutInflater.java:727) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.inflate(LayoutInflater.java:479) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.inflate(LayoutInflater.java:391) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.inflate(LayoutInflater.java:347) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > com.android.internal.policy.impl.PhoneWindow.setContentView(PhoneWindow.java:245) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.app.Activity.setContentView(Activity.java:1786) 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): at > .EsriMapWidget.onCreate(EsriMapWidget.java:155) 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): at > android.app.Instrumentation.callActivityOnCreate(Instrumentation.java:1050) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.app.ActivityThread.performLaunchActivity(ActivityThread.java:1782) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): ... 11 more 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): Caused by: > java.lang.reflect.InvocationTargetException 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): at > java.lang.reflect.Constructor.constructNative(Native Method) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > java.lang.reflect.Constructor.newInstance(Constructor.java:416) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): at > android.view.LayoutInflater.createView(LayoutInflater.java:576) 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): ... 31 more 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): Caused by: java.lang.ExceptionInInitializerError > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > com.esri.android.map.MapView.a(Unknown Source) 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): at com.esri.android.map.MapView.(Unknown Source) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): ... 34 more 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): Caused by: > java.lang.ExceptionInInitializerError 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): at com.esri.android.map.MapSurface.(Unknown > Source) 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): ... 36 more 01-24 > 14:20:31.290: ERROR/AndroidRuntime(4370): Caused by: > java.lang.UnsatisfiedLinkError: Couldn't load GLMapCore: findLibrary > returned null 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > java.lang.Runtime.loadLibrary(Runtime.java:425) 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): at java.lang.System.loadLibrary(System.java:554) > 01-24 14:20:31.290: ERROR/AndroidRuntime(4370): at > com.esri.android.map.Layer.(Unknown Source) 01-24 14:20:31.290: > ERROR/AndroidRuntime(4370): ... 37 more


title: Why does ArcGIS copy & paste layer in a TOC sometimes change the Data Source of the layer?

text: I run into this situation all the time; I often just copy layers and groups of layers in the TOC of one Map Document and paste then into the TOC of another Map Document. For the most part this works fine but often enough the Data Source of layer changes when pasted. As a result the pasted layer has the red exclamation mark beside and no data is displayed until I manually set the Data Source. I confirmed that the Data source of the layer I copy on one document changes when I paste the layer into another document. This is something that has been happening since the early 8.x version and still an issue in version 10 SP3. It's also not computer specific as i run into this issue on several computer. Any ideas? Layer after being pasted: ![enter image description here](http://i.stack.imgur.com/z35pK.png) Source "layer being copied" Data Source: ![enter image description here](http://i.stack.imgur.com/BMDVi.png) Destination "after pasted in another document" Data Source: ![enter image description here](http://i.stack.imgur.com/hu8pG.png)


title: Finding the exact coordinates of circle's extreme east and west points on a sphere

text: I'm currently implementing the Bentley-Ottman algorithm applied to spheres, involving intersections of many spheres on one sphere, circles in most cases of course. This requires locating "event" points on the sphere (ie: points where something "special" happened, such as extreme points of intersection circles or intersections between these). My question is very similar to this one, mainly because I'm asking for the exact same thing : the eastmost and westmost points of a circle on a sphere, in my case named the _start_ and _end_ points of the circle. These points have the particularity of being the tangency points between the circle and a special longitude. The problem resides in finding these as **exact** coordinates. I'm working with the CGAL library, which provides a great deal of exact intersection primitives, and these are in my opinion probably the most appropriate way to find the exact start-end points. So far I've tried to reduce the problem to a simpler one : instead of finding the longitudes of these points, I'm trying to find the height of the horizontal plane which contains these two points. Using the following notations :   * S : start point   * E : end point   * C : center of the circle (not on the sphere of course)   * r : radius of the circle   * O : center of the sphere   * R : radius of the sphere   * h : height of the horizontal plane containing the start-end points   * H : intersection between the pole axis and the horizontal plane containing the start-end points   * phi : the latitude containing C   * phi' : the latitude containing the start-end points. Starting out from the formula based on the Law of Sines, from the last post :               sin(phi') = sin(phi) cos(r/R)      We notably have :               h = R sin(phi')     sin(phi) = OC / R      Giving finally :               h = OC cos(r / R)      This is as far as I've gotten, and I'm stuck on finding an algebraic form of this cosine.


title: Converting KML to layer in ArcGIS 10 - nothing happens

text: I have converted many KML files to layer in ArcGIS many times before without any problem. I just did a conversion (Arc Toolbox>Conversion tools> KML to layer) and the file is showing up in the table of contents but there is no layer. Any ideas or suggestions?? ![enter image description here](http://i.stack.imgur.com/y7W7g.jpg)


title: ArcGIS 10: Adjust custom grid intervals based on map scale?

text: I am preparing a pmf which includes a custom grid overlay. To clarify: a custom grid overlay is marginalia, not a feature class with geometry (so a scale dependecy inside a layer won't work). The final user, in ArcReader, will interact with this map at any scale, zooming in and out as they please. So far, the grid intervals stay the same, regardless of scale, which results in too many grid values at small scales and too few grid values at large scales. Any suggestions? Thank you!


title: How to create one-sided buffers or parallel lines in PostGIS?

text: I'm looking for a way to create so-called one-sided buffers or parallel lines in PostGIS. Example: 1. buffer, 2. one-side, 3. both-sides ![alt text](http://trac.osgeo.org/geos/raw- attachment/ticket/215/singlesidedbuffer.jpg) I've found some discussion on the mailing list from 2009 and information that it has been implemented in GEOS, but nothing about the current PostGIS status. In case the function is not implemented yet, do you know any workarounds? Is it possible to cut off one side of a normal buffer?


title: Serving MBTiles file via MapProxy

text: I'm trying to serve an existent MBTiles file using MapProxy. My YAML file looks like this:               services:       wms:         srs: ['EPSG:3785', 'EPSG:4326', 'EPSG:23028']         image_formats: ['image/jpeg', 'image/png']          layers:       - name: my_mbtiles_layer         title: My MBTiles layer         sources: [my_mbtiles_cache]          cache:       my_mbtiles_cache:         sources: []         grids: [my_grid]         cache:           type: mbtiles           filename: /path/to/myfile.mbtiles          grids:       my_grid:         base: GLOBAL_GEODETIC         res_factor: 'sqrt2'         bbox_srs: 'EPSG:4326'         bbox: [-6.54,38.35,-3,40.65]      So, I defined my own grid, based on GLOBAL_GEODETIC, and used it for my MBTiles cache, that is the only source for my layer. And the queries look like these: http://myserver/service?LAYERS=my_mbtiles_layer&FORMAT=image%2Fpng&SRS=EPSG%3A4326&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&BBOX=-6.54,38.35,-3,40.65&WIDTH=256&HEIGHT=256 All my tiles form a coverage inside the BBOX (in epsg:4326): [-6.54,38.35,-3,40.65] But the problem is I get a blank image of 256x256 px, and a HTTP code 200. So, the response is ok, but nothing is rendered in it. The MBTiles coverage is ok. I can see it in a desktop viewer. Here, an screenshot: http://dl.dropbox.com/u/6599273/capturas/mbtiles_coverage.png So, where can be the problem? Many thanks in advance


title: Illustrator to QGIS workflow for vectors

text: I vectorized a map in Illustrator and would like to import it to QGIS now. I know I must use DXF for that purpose but when I convert the DXF to SHP with the QGIS plugin it is not georeferenced which is actually not surprising. But the QGIS georeferencing tool works for raster images only. So what to do? The base TIFF image I used as the bottom layer in Illu is georeferenced, so maybe I can use the settings from there but I didn't figure out how. I heard you can use the .prj file from another layer but my TIFF doesn't have this, it has a .points only.


title: .NET Entity Framework with PostGIS?

text: Is it possible to use the .NET Entity Framework with PostGIS? I've reviewed a few options, including hookup it up with Npgsql using http://npgsql.com/index.php/2009/08/how-to-set-up-entity-framework-npgsql- part-1/, but it appears that the data types in PostGIS are unsupported in Entity Framework. The closest I can find to a comprehensive statement of what's supported is with Devart's dotConnect for PostgreSQL EF data type mapping at http://www.devart.com/dotconnect/postgresql/docs/DataTypeMapping.html. It appears to expose only base PostgreSQL types and neither PostGIS or presumably common addons like ltree.


title: Downloading vector latlong boundary data of India's Telecom circles from Open street map?

text: I recently came across a couple of really useful sites for downloading vector latlong boundary data for free converts entered address into lat/long coordinates (has a batch mode): http://stevemorse.org/jcal/latlonbatch.html?direction=forward Tool that lets you enter a lat/long coordinate, and then get vector data for the administrative regions that point is in....has an API http://global.mapit.mysociety.org/ The source of the data in the last link is Open Street Map, and these tools together have been great for getting lat/long boundary data of the districts within India's states. However, I'm wondering if the same boundary data exists for India's telecom circles, similar to a map of US area codes. any help would be much appreciated!


title: ArcGIS Select Analysis using Unique Values from Field

text: I am trying to create separate shapefile using Select analysis with ArcPy by doing an SQL query to create them based on a unique value from the field. I want to create a list of the unique values in the table from a specific field since there are some duplicates of the same value in the field. I want to do a search of the table and find all the unique values and then append them to the list and then do the Select analysis to create individual shapefiles and eventually individual KML files. Thanks in advance.


title: Industry-based examples of using Python in geoprocessing and arcpy

text: I was asked recently by some students what types of operations GIS Analyst/Developer are usually trying to automate by using Python with geoprocessing in ArcGIS and ArcPy site package. I guess it might be useful to know when searching for some exercises to do to be sure that what you create is relevant for the industry and might be re-used later already at the workplace. The easiest answer would be "read the Esri help and go through the samples", but I was looking for more specific scenarios which are most common to implement. Thus, sharing workflows sush as "we are getting a .zip file with shapefiles, we use Python to unpack it, project them all to the X coordinate system, load into an SDE geodatabase, grant users access to these data" are most welcome. Please feel comfortable to provide a concise description of the workflow, no extreme details are required. Internet search does not reveal anything on that, so I thought to ask here. The question seems to comply with the FAQ, since it is a question where a specific answer on a certain issue is expected.


title: Why is "Add feature" of the GPS Information tool inactive for some of my PostGIS layers?

text: Sorry If I rise it again on this forum but I couldn't find any info on limitations/requirements of "GPS Information" tool for QGIS (2.2.0) against data layers. For some layers I'm able to collect data for others button "Add feature" is inactive (despite layer is read/write and in an edit mode). My goal is to collect features directly to existing PostGIS layers. Any advise/suggestion welcome.


title: What is the path for font in mapfish print module for windows in config.yaml?

text: The following path is for linux:               fonts:       - /usr/share/fonts/truetype/ttf-dejavu/DejaVuSans      What is the path in order to call font type from windows?


title: How to track ajax WMS requests in OpenLayers interface?

text: I'd like to track WMS requests on my openlayers web interface. How is it possible to log on my server or in Google Analytics every WMS requests (GetMap and GetFeatureInfo would be enough)? Then I'd like to make statistics based on query string in WMS requests. This way I'll know what layers get the most interest from my users.


title: gdal commands cannot import gdal -- possibly ArcGis conflict with separate Python?

text: I installed ArcGIS 10.1 yesterday (Windows 7 64-bit), but have a separate Python 2.7 installation. Now my gdal does not work correctly anymore. It looks like all commands without the ending "*.py" work fine -- e.g. gdal_rasterize, but all gdal commands with the ending "'.py" give me the following error:               Convert  C:\Users\max\icechart_minimum198403_201303.tif  to shapefile.     Traceback (most recent call last):     File "C:\Program Files\GDAL\gdal_polygonize.py", line 34, in <module>      import gdal, ogr, osr     ImportError: No module named gdal      I checked all the variables, even the Registry entries as described here and the Environment path has C:\Program Files\GDAL at the first position and GDAL_DATA is C:\Program Files\GDAL\gdal-data ... By the way, within Python I type "import gdal" and it is imported without any problem ... This all seems right ... any advice on what could be missing / wrong, someone having run into the same problem ?


title: Openlayer 3 - Layer Drawing order

text: I have a map with various vector layers, and I would like to be able to decide which layer goes on top of which. I thought that inserting the layers in the map layer group using: map.getLayers().insertAt(zIndex, layer); would force the layers to be drawn in that specific order, but that does not work. Any suggestions?


title: FME: Using WorkspaceRunner to batch 2 readers

text: I need a little help using the WorkspaceRunner in FME. Essentially I want to batch run polygon shapefiles in one folder (30 of them) against another load of files, in the form of a tiled grid (26 tiles). Unfortunately you can't just put two merged feature type readers into 1 workspace and press Run, the results are all wrong due to the way that FME is reading in the files and processing them. I need to send one shapefile at a time to process against all 26 tiles. Enter the WorkspaceRunner! This is the main processing workspace, 'Workspace 1'; ![Workspace 1](http://i.stack.imgur.com/bpxlK.png) It's not a complicated workflow, just bringing in some Grid tiles and some polygons and doing some joining, taking some attributes and writing out to a .dbf. The writer is fanned out by a unique combo of tile name and Year attribute from the polygon shapefile. This is the WorkspaceRunner, 'Workspace2'; ![workspace 2](http://i.stack.imgur.com/pGw5x.png) Using the Directory and File Pathnames reader (PATH), i am reading the individual file names in the folder where the 30 shapefiles for analysis are kept. In the WorkspaceRunner itself, the .fmw for Workspace 1 is identified, the location of the 26 tiles is auto-inputted into the top Source parameter and the parameter for the location to be written to is also auto-inputted (which should contain 30 shapefiles against 26 grids = 780 .dbf files). Now, the parameter field for the 2nd lot of shapefiles, i.e. the 30 polygon shapefiles, is causing problems. The results will not go beyond the 1st shapefile, i am only getting 26 results. I have used path names, the $(SourceDataset_PATH) user parameter, the path_filename attribute and i get a misture of log error messages or incorrect output (ie only 26 files). The regular log error message is confusing; i am told no .shp exist in the nominated folders, when they clearly do, but when using the User Parameter as shown in the picture, that error goes away. However, being shapefiles, there are pathnames for .prj, .dbx files etc etc and it's just not working! does anyone have any advice at all!?! thanks (FME 2014, 32-bit, Build 14230) (I've read the FMEpedia and FME blog posts about this to no avail)


title: What is a GIS Professional?

text: What is a **GIS Professional**? Is it the same as GIS Developer or Analyst? It's just that the term has been used a lot and I'm not really sure what it means. * * * **Update:** Voted to close the question as it does not conform to the Q&A format. I'm still a bit confused by how people just use the term as if everybody knows what it means when people can't agree on what it is though.


title: How to manually change from GPS (WGS84) to web based maps projection systems?

text: I have collected location point data and wish to display this on a web based map. It was collected using a Garmin Etrex GPS using the WGS84 datum. I wish to display the points using Ushahidi and have done so successfully. The site is here. http://www.flamingocensus.com/ushahidi/ Unfortunately the points do not correlate to what I know they should be, I believe it is a projection error. The points are supposed to be on the edge of the lake. Any suggestions or pointers on how to translate the points in the current projection to my new one.


title: Is there any way I can use a Key-Value store for geospatial data?

text: I have used many relational databases in the past, but I have also read about all NoSQL databases, and the Key-Value stores looks interresting. When I store geometric object I mostly use five indexed columns ID, MIN_X, MAX_X, MIN_Y and MAX_Y (where X and Y are in a map projection). I don't need an index on my other data. I need the X and Y values to lookup objects in a specified place (map rectangle), and I need the ID value if I want to update an specified object. Is there any way that I can use a Key-Value store for this?


title: IDE Problem ArcGIS Runtime SDK for Android

text: I'm trying to develop an app with ArcGIS Runtime SDK for Android. Honestly, I'm not much experienced. Just looking for the best IDE choise. Eclipse works just fine with the Runtime SDK. But has no library documentation in the IDE. ![How Runtime SDK Docs is seen in Eclipse](http://i.stack.imgur.com/dFwGd.png) IntelliJ IDEA 12+ is a good option for Android Development. But I couldn't use the Runtime SDK (v10.1.1) in it. I want to import the documentation to the one of these IDEs. Such as : ![enter image description here](http://i.stack.imgur.com/zC1LO.png)


title: Can draw points but not polylines in WPF Runtime code

text: The following C# code works when drawing in Point mode, but not Polyline . . .                 // "Draw" is ESRI.ArcGIS.Client.Draw object            _draw = new Draw(_map);            _draw.DrawMode = DrawMode.Polyline;            _draw.DrawComplete += (senderObject, drawEventArgs) =>             {                 _barrierLine = new Graphic();                 _barrierLine.Symbol = this.Resources["LineSymbol"] as SimpleLineSymbol;                 _barrierLine.Geometry = drawEventArgs.Geometry;                 _lineGraphicsLayer.Graphics.Add(_barrierLine);             };      This is what the XML code looks like (I tried several different methods):                  <Window.Resources>             <esri:SimpleMarkerSymbol x:Key="PointSymbol" Color="DarkViolet" Size="15" Style="Diamond"></esri:SimpleMarkerSymbol>             <esri:SimpleLineSymbol x:Key="LineSymbol" Color="Aqua" Width="5" Style="Solid"></esri:SimpleLineSymbol>         </Window.Resources>         <Grid>             <Grid.Resources>                 <esri:SimpleMarkerSymbol x:Key="RedMarkerSymbol" Color="Red"/>                 <esri:SpatialReference x:Key="theSpatialReference" WKID="4326"/>                 <esri:SimpleLineSymbol x:Key="LineSymbol" Color="Aqua" Width="3" Style="Solid"/>                   </Grid.Resources>                 <!--Graphics Layer-->                 <esri:GraphicsLayer ID="graphicsLayer" x:Name="_graphicsLayer">                     <esri:Graphic>                                             <esri:Graphic.Symbol>                                                     <esri:SimpleMarkerSymbol Style="Circle" Color="Red"/>                         </esri:Graphic.Symbol>                         <esri:Graphic.Geometry>                             <esri:MapPoint  X="0" Y="0" >                                 <esri:MapPoint.SpatialReference>                                     <esri:SpatialReference WKID="4326"/>                                 </esri:MapPoint.SpatialReference>                             </esri:MapPoint>                                                 </esri:Graphic.Geometry>                                       </esri:Graphic>                  </esri:GraphicsLayer>                 <esri:GraphicsLayer ID="lineGraphicsLayer" x:Name="_lineGraphicsLayer">                     <esri:Graphic>                         <esri:Graphic.Symbol>                             <esri:SimpleLineSymbol Style="Solid" Color="Aqua" Width="3" />                         </esri:Graphic.Symbol>                         <esri:Graphic.Geometry>                             <esri:Polyline>                                 <esri:Polyline.SpatialReference>                                     <esri:SpatialReference WKID="4326"/>                                 </esri:Polyline.SpatialReference>                             </esri:Polyline>                         </esri:Graphic.Geometry>                     </esri:Graphic>                      </esri:GraphicsLayer>             </esri:Map>             <StackPanel Orientation="Vertical" HorizontalAlignment="Left">                 <Button Content="Stops" VerticalAlignment="Top" Name="StopsButton" Margin="5,5,5,0" Click="StopsButton_Click" />                 <Button Content="Point Barriers" Height="23" Name="PointBarriersButton" Width="76" Click="PointBarriersButton_Click" />                 <Button Content="Solve" Name="SolveButton" Margin="5,5,5,0" Click="SolveButton_Click" />                 <Button Content="Line Barriers" Height="23" Name="LineBarriersButton" Width="76" Click="LineBarriersButton_Click" />             </StackPanel>         </Grid>


title: Possible solution for dangling nodes?

text: I was just asked in a interview what a dangling nodes was. I had no idea and asked a former co-worker. I was informed that a significant amount of time was spent fixing these. I was confused how this could become a problem and came up with few quick solutions.   1. Go through the system and identify all endpoints.   2. Reduce the list of endpoints to endpoints that are terminals. (optional just less work)   3. Eliminate all endpoints that aren't within a range of another object (within the same feature type)   4. Create a report   5. Have a system to mark the entry as resolved, if it is a valid entry (so it is not reproduced in the report). For example two roads that are really close but a 1 meter ditch separates them)   6. Repeat regularly. Solution 2:   1. Check for system percolation   2. Mark all entries that do not percolate within a small area   3. Steps 4 through 6 This would produce a larger amount of false positives, but it's running time is in order (N log N) We were just going through the map block by block looking for them. I must be missing something since I can think this problem can be solved relatively simply. I realize that the initial round would be quite a task, but with a system in place maintenance and monitoring would be minimal. What am I missing that makes dangling nodes a large issue?


title: Does anyone know of Electric Company Boundary Shape files that are out there for Free?

text: Ohio has a great website for their Electric Service Commission and I am looking for a shape file or data that would get me the service boundaries of electric companies for other states. EIA has data down to the county. Does anyone know of anyone having it down to the zip code? Or Shapefiles.


title: When converting legacy addins to VS 2010 (ArcGIS 10.1), why are references to certain ArcObjects updated while others are not?

text: I am converting a number of my addins from VS 2008 Express, ArcGIS 10.0 to VS 2010, ArcGIS 10.1 and I need to update some but not all references. Why? For example, in the current addin I just converted to VS 2010, Why would references to _ESRI.ArcGIS.System,ESRI.ArcGIS.Framework, ESRI.ArcGIS.Desktop.AddIns, ESRI.ArcGIS.ArcMapUI_ be updated automatically while references to _ESRI.ArcGIS.Carto, ESRI.ArcGIS.Geodatabase, ESRI.ArcGIS.Geometry_ were not. Until I remove references to ArcObjects version 10.0 and add the appropriate reference to version 10.1 I am getting errors such as: _The project currently contains references to more than one version of ESRI.ArcGIS.Carto, a direct reference to version 10.0.0.0 and an indirect reference (through 'ESRI.ArcGIS.ArcMapUI.IMxDocument.SelectedLayer') to version 10.1.0.0. Change the direct reference to use version 10.1.0.0 (or higher) of ESRI.ArcGIS.Carto._ Would removing ArcObjects 10.0 (DeveloperKit10.0 directory) from the machine correct this issue? Should these even be on the machine after using ArcGIS 10 program remover prior to ArcGIS 10.1 installation? If I can remove these, what is the correct procedure?


title: Where can I find annual temperature data for every city and town in the US and Australia?

text: Where can I find yearly annual averaged temperature data for every city and town in the US and Australia for a heatmap? links are good


title: Is it possible to do spatial interpolation using points in a field?

text: I would like to use Quantum instead of ArcGIS to make field map but before, I need to know if its possible to do spatial interpolation with definited points in a field (ex. pH points). Is there a manual to learn about it?


title: WGS84 CRS without elipsoid

text: What's the difference between using CRS defined as +proj=longlat +datum=WGS84 +no_defs and a CRS defined as +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs Will using one or the other pse any problems. In the first case, how can WGIS project properly without knowing which ellipsoid the data are referred to? Cheers, Andrés


title: Can't get simple map to show up on my test site using CartoDB

text: So I'm trying to get a simple test-site setup but it's not working. Whenever I load the site I can very quickly see the controls flash on the screen then go away. The when not doing local testing, the test-site is located at www.parkru.com/mapTest.html               <!DOCTYPE html>     <html>         <head>             <title>TODO supply a title</title>             <meta charset="UTF-8">             <meta name="viewport" content="width=device-width">                       <link rel="stylesheet" href="http://libs.cartocdn.com/cartodb.js/v3/themes/css/cartodb.css" />             <!--[if lte IE 8]>               <link rel="stylesheet" href="http://libs.cartocdn.com/cartodb.js/v3/themes/css/cartodb.ie.css" />             <![endif]-->             <script type = "text/javascript" src="static/jquery.min.js"></script>                  <script src="https://maps.googleapis.com/maps/api/js?v=3.exp&sensor=false"></script>             <script src="http://libs.cartocdn.com/cartodb.js/v3/cartodb.js"></script>                  <script>                 $(document).ready(function() {                     cartodb.createVis('map', 'http://willkara.cartodb.com/api/v2/viz/2dfc4142-674c-11e3-a9f5-5fd5763ac197/viz.json');                 });                  </script>              </head>         <body>             <div id="map">                  </div>             </body>     </html>


title: Get Descriptive Statistic by Geometry from ArcGIS Server Image Service

text: Image Services have a method called Compute Histogram which is akin to Tabulate Area in ArcGIS Desktop. Is there a method or raster function template method an Image Service can use to get descriptive statistics like MIN, MAX, RANGE, MEAN, etc. for a given geometry (i.e. Zonal Statistics)?


title: Is there any way to symbolize a user-supplied input parameter in ModelBuilder?

text: I have a model built in ESRI ModelBuilder that is run from its tool dialog, where the user supplies the input feature layer that then gets processed in the model. I need the input feature layer to get a certain symbology when the model is run. I've tried symbolizing it by right-clicking the input layer in ModelBuilder, going to its Properties, and setting the layer symbology using a layer file. I then save & run the model, and the input feature layer is still not symbolized - it just stays the same. Am I doing it wrong? I can get output layers in the model to symbolize using layer files just fine with no issues. Is there no way to symbolize inputs... only outputs? I'm also open to using Python to symbolize it, but not sure if this is possible?


title: In PostGIS, how do I merge a multilinestring into a single linestring, keeping the pieces in proper order?

text: I have a multilinestring made up of several pieces, some of which are facing the wrong direction. (Accessible here: http://pastebin.com/89Zcixv5 ) I think what I'm looking for is some utilization of ST_MakeLine, but when I try to use it, the pieces of the trail get put together wrong. (See: http://i5.minus.com/iJ5QVbaA9rm2g.png ) ![enter image description here](http://i.stack.imgur.com/BakkJ.jpg) I have a ton of trails in this format that I need to convert from MultiLineStrings to LineStrings containing the whole trail. If anybody has any ideas, I'd greatly appreciate all the help I can get! Thanks!


title: How to install Mapserver on Ubuntu?

text: I try to install MapServer on Ubuntu by follow this tutorial. I was installed libraries and apach. But have a problem with last command:               libtool: install: /usr/bin/install -c /home/user/mapserver /usr/local/bin/mapserv      I have libtool library in system. But its says that directory `/home/user/mapserver` not exist. When i try:                   libtool install /usr/bin/install -c /home/user/mapserver /usr/local/bin/mapserv      I get:               /usr/bin/install: omitting directory '/home/user/mapserver'      Whats wrong?   PS: first time see Linux OS.


title: PostGIS spatial reference ID (SRID) for regular cartesian coordinate system?

text: I would like to store a polygon as a set of 2D points represented as x, y coordinates of a Cartesian coordinate system. Which spatial reference can I use? Looking at my spatial_ref_sys table, all the spatial references seem to be geography related. The points of the polygon would represent satellite measurements which I would invoke `ST_ConvexHull` on to get the shape of the satellite footprint.


title: QGIS 1.8 python refresh maps

text: i have developed this plugin .. all works fine but the map wont get refreshed . here is code                   vectorCount =0     mpointCount=0     mlineCount=0     mpolygonCount=0     mrasterCount=0     mlegnd = self.iface.legendInterface()     # add the groups and save their index in variable      # self.iface.mapCanvas().refresh()     legnd.addGroup('0_Points')     legnd.addGroup('1_Lines')     legnd.addGroup('2_Polygons')     legnd.addGroup('3_Rasters')     mw = self.iface.mainWindow()     lgd = mw.findChild(QTreeWidget, "theMapLegend")  # get ref to object by type/objectName     lgd.sortItems(0, Qt.AscendingOrder)  # sort first column (Qt.DescendingOrder to reverse)     for layer in legnd.layers():         if layer.type() == QgsMapLayer.VectorLayer:             # layer.setCacheImage( None )             # layer.triggerRepaint()             vectorCount=vectorCount+1             QRect=layer.extent().toString()             f = QgsFeature()             f.setGeometry( QgsGeometry.fromRect( layer.extent()) )             area=f.geometry().area()                 #QMessageBox.information( self.iface.mainWindow(),"Info", " Layer extent area %s "%(area) )             if layer.geometryType() == QGis.Point:                 pointCount=pointCount+1                 legnd.moveLayer(layer, 0)             elif layer.geometryType() == QGis.Line:                 lineCount=lineCount+1                 legnd.moveLayer(layer, 1)             elif layer.geometryType() == QGis.Polygon:                 polygonCount=polygonCount+1                 legnd.moveLayer(layer, 2)             else:                 QMessageBox.information( self.iface.mainWindow(),"Info", " Sorry, This version of plugin supports limited data formats " )         elif layer.type() == QgsMapLayer.RasterLayer:             rasterCount=rasterCount+1             legnd.moveLayer(layer, 3)         else:                 QMessageBox.information( self.iface.mainWindow(),"Info", " Sorry, This version of plugin supports limited data formats " )        i tried all this but it won't work ... any suggestions? what i am doing wrong ?                   self.iface.setVisible(True);         self.iface.refresh()          self.iface.mapCanvas().refresh()     layer.setCacheImage( None )     layer.triggerRepaint()


title: WMSGetFeatureInfo returns an empty popup

text: I am building a web map application using Openlayers and Geoserver. Both the application and the geoserver are hosted on amazon. I would like to include a feature info tool and am using WMSGetFeatureInfo to perform WMS query. However, when I click on the map, it returns an empty popup. Here is how I set up the proxy host and the feature info control in the html file:                // Set up Proxy Host     OpenLayers.ProxyHost = "proxy.php?url=";               // Add feature info control                 var featureInfo = new OpenLayers.Control.WMSGetFeatureInfo({                     url: 'http://my_site:8080/geoserver/opengeo/wms',                      layerUrls: ["http://my_site:8080/geoserver/gwc/service/wms"],                     title: 'Identify features by clicking',                     queryVisible: true,                     INFO_FORMAT:'text/html',                     eventListeners: {                         "getfeatureinfo": function(event) {                             // Show Popup                             map.addPopup(new OpenLayers.Popup.FramedCloud(                             "chicken",                              map.getLonLatFromPixel(event.xy),                             null,                             event.text,                             null,                             true                             ));                         }                     }                 });                 map.addControl(featureInfo);                 featureInfo.activate();      when i try to load http://my_site.com/proxy.php?url=http://www.google.com on my browser. I got the following message: _This page contains the following errors: error on line 2 at column 1: Extra content at the end of the document Below is a rendering of the page up to the first error._ The following is my proxy script (proxy.php):                <?php      /*     License: LGPL as per: http://www.gnu.org/copyleft/lesser.html     $Id: proxy.php 3650 2007-11-28 00:26:06Z rdewit $     $Name$     */          ////////////////////////////////////////////////////////////////////////////////     // Description:     // Script to redirect the request http://host/proxy.php?url=http://someUrl     // to http://someUrl .     //     // This script can be used to circumvent javascript's security requirements     // which prevent a URL from an external web site being called.     //     // Author: Nedjo Rogers     ////////////////////////////////////////////////////////////////////////////////          // read in the variables          if(array_key_exists('HTTP_SERVERURL', $_SERVER)){        $onlineresource=$_SERVER['HTTP_SERVERURL'];      }else{        $onlineresource=$_REQUEST['url'];      }     $parsed = parse_url($onlineresource);     $host = @$parsed["host"];     $path = @$parsed["path"] . "?" . @$parsed["query"];     if(empty($host)) {       $host = "localhost";     }     $port = @$parsed['port'];     if(empty($port)){       $port="80";     }     $contenttype = @$_REQUEST['contenttype'];     if(empty($contenttype)) {       $contenttype = "text/xml";     }     $data = @$GLOBALS["HTTP_RAW_POST_DATA"];     // define content type     header("Content-type: " . $contenttype);          if(empty($data)) {       $result = send_request();     }     else {       // post XML       $posting = new HTTP_Client($host, $port, $data);       $posting->set_path($path);       $result = $posting->send_request();     }          // strip leading text from result and output result     $len=strlen($result);     $pos = strpos($result, "<");     if($pos > 1) {       $result = substr($result, $pos, $len);     }         //$result = str_replace("xlink:","",$result);     echo $result;          // define class with functions to open socket and post XML     // from http://www.phpbuilder.com/annotate/message.php3?id=1013274 by Richard Hundt          class HTTP_Client {        var $host;        var $path;       var $port;        var $data;        var $socket;        var $errno;        var $errstr;        var $timeout;        var $buf;        var $result;        var $agent_name = "MyAgent";        //Constructor, timeout 30s        function HTTP_Client($host, $port, $data, $timeout = 30) {          $this->host = $host;          $this->port = $port;         $this->data = $data;         $this->timeout = $timeout;        }             //Opens a connection        function connect() {          $this->socket = fsockopen($this->host,            $this->port,            $this->errno,            $this->errstr,            $this->timeout            );          if(!$this->socket)            return false;          else            return true;        }             //Set the path        function set_path($path) {          $this->path = $path;        }             //Send request and clean up        function send_request() {          if(!$this->connect()) {            return false;          }          else {            $this->result = $this->request($this->data);           return $this->result;          }        }             function request($data) {          $this->buf = "";          fwrite($this->socket,            "POST $this->path HTTP/1.0\r\n".            "Host:$this->host\r\n".            "User-Agent: $this->agent_name\r\n".            "Content-Type: application/xml\r\n".            "Content-Length: ".strlen($data).            "\r\n".            "\r\n".$data.            "\r\n"          );               while(!feof($this->socket))            $this->buf .= fgets($this->socket, 2048);            $this->close();            return $this->buf;        }                  function close() {         fclose($this->socket);       }      }                     function send_request() {       global $onlineresource;       $ch = curl_init();       $timeout = 5; // set to zero for no timeout            // fix to allow HTTPS connections with incorrect certificates       curl_setopt($ch, CURLOPT_SSL_VERIFYPEER, FALSE);       curl_setopt($ch, CURLOPT_SSL_VERIFYHOST, 2);            curl_setopt ($ch, CURLOPT_URL,$onlineresource);       curl_setopt ($ch, CURLOPT_RETURNTRANSFER, 1);       curl_setopt ($ch, CURLOPT_CONNECTTIMEOUT, $timeout);       curl_setopt ($ch, CURLOPT_ENCODING , "gzip, deflate");            $file_contents = curl_exec($ch);       curl_close($ch);       $lines = array();       $lines = explode("\n", $file_contents);       if(!($response = $lines)) {         echo "Unable to retrieve file '$service_request'";       }       $response = implode("",$response);       return $response;     }     ?>


title: Why is my linestring not valid ? was recorded from GPX file

text: I have a linestring that was built from points from GPX file. I've insert it to sql as geometry and when I try to run funcs on it (like STStartPoint ) I get error that line isn't valid and I need to run MakeValid on it (which change the GPS Points). Here is the line :               LINESTRING (34.7873 31.3182, 34.7884 31.3366, 34.7941 31.3552, 34.7942 31.3712, 34.793 31.3756, 34.7892 31.3906, 34.7883 31.4086, 34.7864 31.423, 34.7863 31.4219, 34.7863 31.4221, 34.7863 31.422, 34.7862 31.4216, 34.7862 31.4219, 34.7849 31.4241, 34.7734 31.435, 34.7656 31.4466, 34.7734 31.4668, 34.7798 31.491, 34.7826 31.5068, 34.7927 31.5311, 34.8005 31.5514, 34.8065 31.5749, 34.8132 31.5973, 34.8147 31.6196, 34.8095 31.657, 34.806399999999996 31.6748, 34.8157 31.6988, 34.8212 31.7214, 34.8227 31.7407, 34.8264 31.7604, 34.8328 31.7959, 34.8401 31.8035, 34.8576 31.8097, 34.8766 31.809, 34.8998 31.8109, 34.9156 31.8146, 34.9208 31.8193, 34.9487 31.826, 34.9687 31.8311, 34.9774 31.8349, 34.9881 31.8345, 34.993 31.8268, 35.0003 31.8261, 35.0056 31.8259, 35.0177 31.8215, 35.0216 31.8183, 35.0269 31.8133, 35.0314 31.8098, 35.0382 31.8062, 35.0454 31.8047, 35.0532 31.805, 35.0649 31.8041, 35.0758 31.8003, 35.0884 31.8003, 35.1064 31.8015, 35.1261 31.8007, 35.147 31.7972, 35.1622 31.7933, 35.1716 31.7977, 35.1877 31.7979, 35.1884 31.7965, 35.1887 31.7961, 35.1905 31.795, 35.1926 31.7946, 35.1927 31.7936, 35.1936 31.7931, 35.1947 31.7934, 35.1961 31.7916, 35.197 31.7908, 35.1979 31.7902, 35.1987 31.7899, 35.1993 31.7872, 35.1985 31.7854, 35.1994 31.7847, 35.1992 31.7823, 35.1984 31.7814, 35.2012 31.7761, 35.2008 31.7747, 35.2045 31.7747, 35.209 31.7762, 35.2086 31.782, 35.2098 31.7825, 35.2114 31.7817, 35.2135 31.7808, 35.2148 31.7817, 35.2147 31.7809, 35.2126 31.7812, 35.2086 31.7831, 35.2083 31.783, 35.2032 31.7822, 35.1985 31.7863, 35.1998 31.7875, 35.2001 31.788, 35.1973 31.7905, 35.1911 31.7951, 35.1816 31.8017, 35.1687 31.7957, 35.1592 31.7958, 35.1194 31.7988, 35.0804 31.8009, 35.0632 31.8051, 35.0481 31.804, 35.0329 31.8085, 35.0209 31.8188, 34.9901 31.8293, 34.988 31.839, 34.978 31.8351, 34.9499 31.826, 34.9274 31.8231, 34.917 31.8156, 34.9037 31.8119, 34.883 31.8105, 34.86 31.8101, 34.84 31.8036, 34.8302 31.7943, 34.8325 31.7811, 34.8271 31.7626, 34.8228 31.7437, 34.8205 31.7209, 34.8142 31.6955, 34.8063 31.6752, 34.81 31.653199999999998, 34.8148 31.6239, 34.8133 31.6027, 34.808 31.5828, 34.803 31.5607, 34.7972 31.5415, 34.7868 31.5232, 34.783 31.5045, 34.7761 31.484, 34.7734 31.4673, 34.7665 31.4504, 34.7659 31.4422, 34.7853 31.4222, 34.7891 31.4016, 34.7913 31.3816, 34.7957 31.3653, 34.7921 31.3504, 34.7865 31.3298, 34.7884 31.3094, 34.7908 31.2899, 34.7976 31.2879, 34.8037 31.2868)      and this is that I get after MakeValid()               MULTILINESTRING ((34.991771945003556 31.828729956620826, 34.989424975190637 31.83241805204112, 34.990100026130676 31.829299926757813, 34.991771945003556 31.828729956620826), (34.954345241759384 31.827439520434716, 34.968700051307678 31.831099987030029, 34.977399945259094 31.83489990234375, 34.988100051879883 31.834500074386597, 34.989424975190637 31.83241805204112, 34.98800003528595 31.83899998664856, 34.978000044822693 31.835099935531616, 34.954345241759384 31.827439520434716), (34.999167677687957 31.8262086938732, 34.991771945003556 31.828729956620826, 34.993000030517578 31.826800107955933, 34.999167677687957 31.8262086938732), (34.947309795820473 31.825666136443175, 34.949900031089783 31.825999975204468, 34.954345241759384 31.827439520434716, 34.948699951171875 31.825999975204468, 34.947309795820473 31.825666136443175), (34.922794889315966 31.81977900590601, 34.947309795820473 31.825666136443175, 34.92739999294281 31.823100090026855, 34.922794889315966 31.81977900590601), (35.022785521481737 31.817181583739565, 35.021600008010864 31.818300008773804, 35.017699956893921 31.821500062942505, 35.005599975585938 31.825900077819824, 35.0003000497818 31.82610011100769, 34.999167677687957 31.8262086938732, 35.020900011062622 31.81879997253418, 35.022785521481737 31.817181583739565), (34.9165756262933 31.815481844098283, 34.91700005531311 31.815599918365479, 34.922794889315966 31.81977900590601, 34.920799970626831 31.819299936294556, 34.9165756262933 31.815481844098283), (35.031246560569912 31.809919226177239, 35.022785521481737 31.817181583739565, 35.026900053024292 31.813299894332886, 35.031246560569912 31.809919226177239), (35.035067600420085 31.807858299827611, 35.031399965286255 31.809799909591675, 35.031246560569912 31.809919226177239, 35.032899975776672 31.808500051498413, 35.035067600420085 31.807858299827611), (34.844628452901595 31.805104311964445, 34.858650011538806 31.809661326879048, 34.876600027084351 31.809000015258789, 34.899799942970276 31.810899972915649, 34.915599942207336 31.814599990844727, 34.9165756262933 31.815481844098283, 34.903699994087219 31.8118999004364, 34.883000016212463 31.810499906539917, 34.860000014305115 31.810100078582764, 34.858650011538806 31.809661326879048, 34.857599973678589 31.809700012207031, 34.844628452901595 31.805104311964445), (35.188035494269172 31.797229039308949, 35.187700033187866 31.797899961471558, 35.187080863729726 31.797892268672587, 35.181599974632263 31.801700115203857, 35.173038050284156 31.797717795369007, 35.187080863729726 31.797892268672587, 35.188035494269172 31.797229039308949), (35.14965103796618 31.796519744421513, 35.146999955177307 31.797199964523315, 35.126099944114685 31.800699949264526, 35.106400012969971 31.801500082015991, 35.08980497387266 31.800393562901903, 35.080399990081787 31.800899982452393, 35.063199996948242 31.805099964141846, 35.057397041666896 31.804677208028728, 35.053200006484985 31.805000066757202, 35.045697030633413 31.804711327578612, 35.035067600420085 31.807858299827611, 35.0382000207901 31.80620002746582, 35.045400023460388 31.804699897766113, 35.045697030633413 31.804711327578612, 35.048099994659424 31.803999900817871, 35.057397041666896 31.804677208028728, 35.064900040626526 31.804100036621094, 35.0757999420166 31.800299882888794, 35.08840000629425 31.800299882888794, 35.08980497387266 31.800393562901903, 35.119400024414062 31.798799991607666, 35.14965103796618 31.796519744421513), (35.191488436823214 31.794811774547405, 35.191100001335144 31.795099973678589, 35.188035494269172 31.797229039308949, 35.188400030136108 31.796499967575073, 35.188699960708618 31.796099901199341, 35.190500020980835 31.795000076293945, 35.191488436823214 31.794811774547405), (35.192666237926815 31.793937906774794, 35.192600011825562 31.794600009918213, 35.191488436823214 31.794811774547405, 35.192666237926815 31.793937906774794), (35.162199974060059 31.793299913406372, 35.1673577298707 31.795714187650482, 35.168699979782104 31.795700073242187, 35.173038050284156 31.797717795369007, 35.171599984169006 31.797699928283691, 35.1673577298707 31.795714187650482, 35.159199953079224 31.795799970626831, 35.14965103796618 31.796519744421513, 35.162199974060059 31.793299913406372), (35.198880189308213 31.789089100267343, 35.198699951171875 31.789900064468384, 35.197900056838989 31.790199995040894, 35.197000026702881 31.790800094604492, 35.196099996566772 31.791599988937378, 35.194700002670288 31.793400049209595, 35.193742910444314 31.79313907120585, 35.192666237926815 31.793937906774794, 35.192700028419495 31.793600082397461, 35.1936000585556 31.793100118637085, 35.193742910444314 31.79313907120585, 35.197299957275391 31.790499925613403, 35.198880189308213 31.789089100267343), (34.831336083280775 31.787779747905049, 34.832800030708313 31.795900106430054, 34.840100049972534 31.803499937057495, 34.844628452901595 31.805104311964445, 34.8400000333786 31.803600072860718, 34.830199956893921 31.794300079345703, 34.831336083280775 31.787779747905049), (35.199178294521381 31.786926018042656, 35.19980001449585 31.787499904632568, 35.200099945068359 31.788000106811523, 35.198880189308213 31.789089100267343, 35.199300050735474 31.787199974060059, 35.199178294521381 31.786926018042656), (35.198788290813596 31.78604849517162, 35.199178294521381 31.786926018042656, 35.1985000371933 31.78629994392395, 35.198788290813596 31.78604849517162), (35.2007999420166 31.774699926376343, 35.204499959945679 31.774699926376343, 35.208999991416931 31.77620005607605, 35.208600044250488 31.782000064849854, 35.209800004959106 31.782500028610229, 35.211400032043457 31.781699895858765, 35.213500022888184 31.780800104141235, 35.213824916937241 31.781024982663222, 35.2146999835968 31.780900001525879, 35.214800000190735 31.781699895858765, 35.213824916937241 31.781024982663222, 35.212599992752075 31.781199932098389, 35.208600044250488 31.783099889755249, 35.208299994468689 31.782999992370605, 35.203199982643127 31.78220009803772, 35.198788290813596 31.78604849517162, 35.1985000371933 31.785399913787842, 35.199399948120117 31.7846999168396, 35.199200034141541 31.782299995422363, 35.198400020599365 31.781399965286255, 35.201200008392334 31.776099920272827, 35.2007999420166 31.774699926376343), (34.825456178748283 31.755374636510627, 34.827100038528442 31.762599945068359, 34.832499980926514 31.781100034713745, 34.831336083280775 31.787779747905049, 34.8264000415802 31.760400056838989, 34.825456178748283 31.755374636510627), (34.810187663125 31.652054055238416, 34.809499979019165 31.657000064849854, 34.8063999414444 31.674799919128418, 34.8157000541687 31.698800086975098, 34.821200013160706 31.721400022506714, 34.822700023651123 31.740700006484985, 34.825456178748283 31.755374636510627, 34.822800040245056 31.74370002746582, 34.820500016212463 31.720900058746338, 34.814200043678284 31.695499897003174, 34.806300044059753 31.67519998550415, 34.809999942779541 31.653199911117554, 34.810187663125 31.652054055238416), (34.814564661037444 31.62057350543299, 34.814800024032593 31.623899936676025, 34.810187663125 31.652054055238416, 34.814564661037444 31.62057350543299), (34.803976636563448 31.565016789967814, 34.80649995803833 31.574899911880493, 34.813199996948242 31.597300052642822, 34.81470000743866 31.619600057601929, 34.814564661037444 31.62057350543299, 34.813300013542175 31.602699995040894, 34.807999968528748 31.582799911499023, 34.803976636563448 31.565016789967814), (34.802189379692969 31.558016620660805, 34.802999973297119 31.560699939727783, 34.803976636563448 31.565016789967814, 34.802189379692969 31.558016620660805), (34.799052621581232 31.547632967487864, 34.80050003528595 31.551399946212769, 34.802189379692969 31.558016620660805, 34.799052621581232 31.547632967487864), (34.795644399658784 31.538762870825277, 34.797199964523315 31.541500091552734, 34.799052621581232 31.547632967487864, 34.795644399658784 31.538762870825277), (34.78429707291923 31.510882897847225, 34.792700052261353 31.531100034713745, 34.795644399658784 31.538762870825277, 34.786800026893616 31.523200035095215, 34.78429707291923 31.510882897847225), (34.781294481037484 31.4994327924487, 34.782999992370605 31.504499912261963, 34.78429707291923 31.510882897847225, 34.782600045204163 31.506799936294556, 34.781294481037484 31.4994327924487), (34.766367359702791 31.448587318836147, 34.773399949073792 31.466799974441528, 34.779800057411194 31.490999937057495, 34.781294481037484 31.4994327924487, 34.776100039482117 31.483999967575073, 34.773399949073792 31.467299938201904, 34.7664999961853 31.450400114059448, 34.766367359702791 31.448587318836147), (34.766160912887024 31.445765728619, 34.766367359702791 31.448587318836147, 34.765599966049194 31.446599960327148, 34.766160912887024 31.445765728619), (34.78629994392395 31.422100067138672, 34.78629994392395 31.421999931335449), (34.788482692723164 31.404946598016188, 34.788300037384033 31.408600091934204, 34.786399960517883 31.42300009727478, 34.78629994392395 31.421900033950806, 34.78629994392395 31.421999931335449, 34.786200046539307 31.421600103378296, 34.786200046539307 31.421900033950806, 34.784899950027466 31.424099922180176, 34.773399949073792 31.434999942779541, 34.766160912887024 31.445765728619, 34.765900015830994 31.442199945449829, 34.7853000164032 31.422199964523315, 34.788482692723164 31.404946598016188), (34.791161962913016 31.382855310863285, 34.789100050926208 31.401599884033203, 34.788482692723164 31.404946598016188, 34.789199948310852 31.390599966049194, 34.791161962913016 31.382855310863285), (34.794197849383437 31.370864674378371, 34.79419994354248 31.371200084686279, 34.792999982833862 31.3756000995636, 34.791161962913016 31.382855310863285, 34.791300058364868 31.381599903106689, 34.794197849383437 31.370864674378371), (34.794122337747623 31.35877037755764, 34.7956999540329 31.365299940109253, 34.794197849383437 31.370864674378371, 34.794122337747623 31.35877037755764), (34.787409615770848 31.320033730460121, 34.788400053977966 31.336600065231323, 34.794100046157837 31.355200052261353, 34.794122337747623 31.35877037755764, 34.792099952697754 31.3503999710083, 34.786499977111816 31.329799890518188, 34.787409615770848 31.320033730460121), (34.803699970245361 31.286799907684326, 34.797600030899048 31.2878999710083, 34.7907999753952 31.289900064468384, 34.788400053977966 31.309400081634521, 34.787409615770848 31.320033730460121, 34.787299990653992 31.31820011138916))      I'm trying to show the line and auto put flag on start and end of it.


title: Which is the best way to create a polygon with the extents of a raster file in QGIS?

text: I have a tiff image file and want to get a polygon with the extent of the tiff. Is there a direct way of doing this? Right now i have to convert the image to vector (with gdal_polygonize) and then dissolve the polygons...


title: What are best practices to compare densities and locations between two sets of points?

text: I have two sets of point shapefiles. They represent two different customer groups for a retail store. Shapefile A represents online customers, Shapefile B represents in-store customers. **What I would like to do is create some sort of heat map that displays where the concentrations of these customers are.** That way I can hopefully show the differences in where these two different sets of customers are located in comparison to our actual store locations. I even have how much money each individual customer has spent attributed in both shapefiles if I wanted to organize this by amount spent. What do you think would be the best methodology to show this? I have the Spatial Analyst and Business Analyst extensions installed.


title: What's the use of the table extended_geometry_columns in PostGIS

text: Sometimes my PostGIS database creates a extended_geometry_columns table. After some search i've not found what is the use of this table. I'm using PostgreSQL 8.4 and PostGIS 1.3.3 This is the definition of the table               CREATE TABLE extended_geometry_columns (         f_table_schema character varying(265) NOT NULL,         f_table_name character varying(256) NOT NULL,         x_min double precision,         y_min double precision,         x_max double precision,         y_max double precision,         CONSTRAINT extended_geometry_columns_pkey PRIMARY KEY (f_table_schema, f_table_name)      ); AFAIK the only tools that are connecting to the database are pgadmin3/psql, gvSIG, and some custom python scripts (that i'm sure not create the table). Thanks


title: Problems configuring MapProxy to serve WMS layer

text: I'm having problems making MapProxy render a Mapserver WMS layer. My configuration file looks like this:               services:      wms:        srs: ['EPSG:3785', 'EPSG:4326', 'EPSG:23028']        image_formats: ['image/jpg', 'image/png']          layers:      name: mylayer      title: My layer      sources: [my_wms_source]          sources:      my_wms_source:        type:wms        req:          url: http://mydomain/cgi-bin/mapserv?          layers: mylayer_in_mapserver          map: /path/to/my/mapfile.map      Direct queries to mapserver work fine. So, a query like this works: http://mydomain/cgi- bin/mapserv?styles=&map=/path/to/my/mapfile.map&format=image/png&width=800&height=600&bbox=-18.4267,27.3264,-13.0047,29.5457&layers=mylayer_in_mapserver&service=WMS&request=GetMap&srs=EPSG:4326&version=1.1.1 The key, I think, is in the bounding box. If I deploy MapProxy with debug activated and try to see my WMS source through the demo app, queries like these are raised: http://mydomain/cgi- bin/mapserv?styles=&map=map=/path/to/my/mapfile.map&format=image%2Fpng&width=800&height=600&bbox=-506.77734375,-266.8359375,370.01953125,190.1953125&layers=mylayer_in_mapserver&service=WMS&request=GetMap&srs=EPSG:4326&version=1.1.1 And that query doesn't work, because the bounding box is in a different coordinate system (I guess epsg:3785). How could I control that? The bbox is used for WMS queries. I've tried the 'srs' parameter in the wms option of the 'services' section, but it doesn't work. Any clues? Thanks in advance, and best regards


title: Carto DB - How to manipulate the drawing order of point data?

text: I’m doing a bubble/choropleth combo map with point data (cities) using CSS in CartoDB, and I want my biggest bubbles (that correspond with population size) to draw on the map first so that the map is easier to read and the infowindow feature is easier to use. I’ve assigned IDs to my data according the population size so that the city with the largest population has carto_id 1 and so forth (that’s how I would do it with ArcMap, because it draws ID 1 first). I’ve also tried to write my CSS so that the largest bubbles would be drawn first (I have little CSS experience). The CSS is below. I can’t figure out what order it’s drawing the bubbles in. How can I fix this?               #tm_db001_merge2{       marker-fill: #F11810;       marker-line-color: #FFF;       marker-line-width: 0.5;       marker-line-opacity: 0.6;       marker-opacity: 0.6;       marker-comp-op: screen;       marker-placement: point;       marker-type: ellipse;       marker-allow-overlap: true;       marker-clip: false;       marker-multi-policy: largest;     }     #tm_db001_merge2 [ population <= 2880546] {        marker-width: 62.0;       }     #tm_db001_merge2 [ population <= 2000000] {        marker-width: 52.0;     }     #tm_db001_merge2 [ population <= 1000000] {        marker-width: 43.0;     }     #tm_db001_merge2 [ population <= 750000] {        marker-width: 35.0;     }     #tm_db001_merge2 [ population <= 500000] {        marker-width: 28.0;     }     #tm_db001_merge2 [ population <= 250000] {        marker-width: 22.0;     }     #tm_db001_merge2 [ population <= 100000] {        marker-width: 17.0;     }     #tm_db001_merge2 [ population <= 75000] {        marker-width: 13.0;     }     #tm_db001_merge2 [ population <= 50000] {        marker-width: 10.0;     }     #tm_db001_merge2 [ population <= 25000] {        marker-width: 8.0;     }               #tm_db001_merge2 [ percent_inner_city <= 100] {        marker-fill: #e31a1c;     }     #tm_db001_merge2 [ percent_inner_city < 75] {        marker-fill: #fd8d3c;     }     #tm_db001_merge2 [ percent_inner_city < 50] {        marker-fill: #fecc5c;     }     #tm_db001_merge2 [ percent_inner_city < 25] {        marker-fill: #ffffb2;     }


title: Two r.mapcalc calculations return a much bigger raster cell size than the input raster

text: I tried the following calculations with a GeoTiff sample, having a cell size of 25m:   1. Calculating the togographic exposure (TOPEX) and   2. estimating the mean topographic curvature Both, the resulting raster data from calculation 1. and 2. have a raster cell size of approximately 765m. Even when I choose a 1m GeoTiff, the resulting raster has a cell size of 765m. Obviously it doesn't depend on the calculation. Could be the problem that the input file is not type of `.dem`, `.adf` or `.asc`??


title: How can I create a merged buffer in QGIS then count points that fall within it

text: I have a QGIS plan which shows ponds (as points in a shape file with an OS tile as a background) within a survey area. I need to calculate the total number of ponds within 1 km of each of the on-site ponds individually. I have created one layer with 1 km buffers for each pond as a shapefile. As there are many ponds, this is messy to look at and I would like to create a layer which has the outline of all of the 1 km buffers combined. I have tried selecting the layer with the buffers and vector-geoprocessing tools - dissolve, but the new layer had outlines for each individual buffer. **Question part 1 - how do you create combined buffers, or convert existing buffers into one?** I will then use this outline to manually check the OS plan and add points representing ponds as a new layer "off site ponds". **Question part 2 - How would I then find the number of ponds (onsite and offsite, not counting the pond being queried) within 1 km of each on-site pond (individually).** _Sub-question - if one of the buffers dissects a symbol for a pond will it count the pond, or only if the centre point falls within the buffer? (i.e. could the size of the symbol affect the result?)_


title: Operation on mixed SRID geometries

text: [I'm new to the site and my English is poor, so it is difficult for me to explain my problem. But there is no other person I can ask.] My question is related to this exercise http://workshops.boundlessgeo.com/postgis-intro/joins.html, but it gives an error when I run it.               SELECT     subways.name AS subway_name,     neighborhoods.name AS neighborhood_name,     neighborhoods.boroname AS borough     FROM nyc_neighborhoods AS neighborhoods     JOIN nyc_subway_stations AS subways     ON ST_Contains(neighborhoods.geom, subways.geom)     WHERE subways.name = 'Broad St';      The error is `Operation on mixed SRID geometries`


title: Finding the center of a quadkey

text: trying to get a center point given a quadkey like "120023231323". Have some data i'd like to display on a map for that quadkey but can only pass a lat/lon to the data visualization tool i'm playing with.


title: Is it possible to use polygons to create heatmaps?

text: I hope I can explain my question clearly, so here it goes: I've been fiddling around with heatmaps on QGIS for a few days, but came across a small obstacle on my way. I have a dataset comprised of point data, corresponding to objects inside a 1 sq meter grid. Each of those objects has several data associated, namely their weight (in grams). So, I have already created a simple heatmap based on the number of points and another one using the grams value of each point as the weight parameter of the heatmap. This resulted in two different heatmaps, both of them interesting to my analysis. Nevertheless, I have a different dataset that comes from the same grid - this dataset corresponds to objects that I know came from a specific square on the grid, but to which I don't have a specific point location (as I said, I only know they came from a specific square, but I don't have XY coordinates for them inside that square). I have all the other data for these objects, namely their weight and I wanted to use them to create a new heatmap, where I would attempt to represent the distribution of objects through the grid, based upon their weight. I was thinking of creating a new attribute on my grid shapefile, where I would insert the total object weight (both the ones to which I have XY coordinates and the ones to which I don't have that information) for each square. I've been looking at the heatmap plugin on QGIS, but it only allows me to use point shapefiles as the input layer... Is there a way of creating a heatmap using polygons as the input layer? Remember, the polygons do not overlap, so I cannot use the approaches sugested elsewhere, regarding overlapping polygons. Solutions based upon QGIS, GRASS or SAGA would be preferable, since I don't have access to proprietary software at the time being. I have already created a "points per polygon" shapefile and I guess I could also create a "total weight por polygon" shapefile, but I'd rather have a heatmap, since it would give out a less "boxy" visualization. All help will be greatly appreciated! Many thanks in advance!


title: How to count the number of data points in each area in QGIS?

text: I have a shapefile (.shp file) which defines a large number (several thousand) separate areas. I also have an Excel spreadsheet which contains a large number of data points. Each data point consists of a latitude and longitude, and a "weight". I want to determine, for each area in the shape file, the total weight of all the data points whose lat/long is in that area. I would prefer to have this information in a format that can be downloaded back into the Excel spreadsheet. What is the best way to do this in QGIS?


title: Define Workspace for SDE Connection in Python

text: How do we define SDE Connection for workspace in Python Scripting?


title: mod09a1 STATE FLAG

text: I am using MODIS MOD09A1 surface reflectance data.For each 8 day composite ,one state flg is given which is to be used for extracting a cloudy pixels. MODIS state flaga are stored in 16 bit integer format. To extract the information form state flags its necessary to convert it to binary format and select a bit word ahich indicates cloudy pixel as per given in MODIS Manual.How to convert this 16 bit to binary format in envi or erdas to build a cloud mask .Is thier any IDL programme?


title: Where are descriptions of attributes in datafiles on naturalearthdata.com?

text: I have downloaded some datafiles from naturalearthdata.com site. For example, `ne_50m_admin_0_countries` feature source contains dozenz of attributes with names like `scalerank`, `level`, `mapcolor7`, `mapcolor8`, `mapcolor13` and so on. It is possible to guess the meaning of some of them but others are mysterious. Is there any description for these attribute names?


Insert Cursor fails to write results to table. The Search Cursor is returning the results, but errors when trying to insert into field.               # Insert Cursor for GDB table to GDB table       # Import modules, environment settings     import arcpy, traceback, string     arcpy.env.overwriteOutput = True          # Database connection variables on SPGISPROC1     OutFGDB = r"P:\batch\APO\District_overlay.gdb\pa_fp_dat1"     InFGDB = r"P:\batch\APO\District_overlay.gdb\pa_fp_rel1"          #Create cursor for input table items     InCursor = arcpy.SearchCursor(InFGDB)          # Create a variable that stores the value for column in a given row of the input table          for InRow in InCursor:              parcel      = InRow.getValue('pafp_frq_PARCEL')           fp_zone1    = InRow.getValue('pafp_frq_ZONE2')               # Copy the stored value from the input table into the output table         gdbCursor = arcpy.InsertCursor(OutFGDB)         row = gdbCursor.newRow()              if fp_zone1:             fp_zone1 = fp_zone1.upper              #Fill out fields with variables from input table                 row.pafp_frq_PARCEL = parcel         row.pafp_frq_ZONE2  = fp_zone1                  #Insert new row into output file          gdbCursor.insertRow(row)          #Delete cursor and row objects to remove locks on the data      del row     del cursor      Error Message:               Traceback (most recent call last):       File "C:\Python27\ArcGISx6410.1\Lib\site-packages\pythonwin\pywin\framework\scriptutils.py", line 325, in RunScript         exec codeObject in __main__.__dict__       File "C:\gisdata\Python\Insert_table_Cursor.py", line 28, in <module>         row.pafp_frq_PARCEL = parcel       File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\arcobjects\_base.py", line 35, in __setattr__         return setattr(self._arc_object, attr, ao)     RuntimeError: ERROR 999999: Error executing function.


title: How to fix crashes after update to Ubuntu 12.04?

text: Some weeks ago I had updated Ubuntu at the 11.10 version, but I had started having problems with QGIS 1.8: I installed the program and launched it, but it didn't open and crashed continuously. I decided to update the operating system again, to the 12.04 version, but this hasn't solved the problems. So I've removed QGIS and all its related packages, in order to re-install it with the PPA from the website (for Precise):               deb http://qgis.org/debian precise main     deb-src http://qgis.org/debian precise main      To add the the qgis.org repository public key to your apt keyring, type:               gpg --keyserver keyserver.ubuntu.com --recv 997D3880     gpg --export --armor 997D3880 | sudo apt-key add -      So the program is installed and it starts, but it creates an unstable window and often crashes: where is the problem? I have two options: A) my laptop is over and unable to cope with GIS requirements B) there's something disturbing the GIS that lasts across the OS updatings.


title: How do different GIS systems determine a Polygon's interior?

text: I'm interested in understanding how real-world GIS systems and their data encode Polygons. Specifically, how do they resolve the ambiguity of a Polygon's interior on a sphere? Background: in 2D, it's trivial to pick the side of the boundary that has a finite area, since the 2D plane is infinite. However, a sphere is finite, so it's impossible to know which side is inside without making additional assumptions. Possible approaches that I know of:   1. **Right-hand-rule** : exterior boundaries are always specified in clockwise order, and holes are specified in anticlockwise order. (There is of course the left-hand-rule, too).   2. **Smallest area** : for any given ring, always pick the side with the smallest area. I'm not sure how you would specify a large-range Polygon: perhaps an empty exterior ring followed by holes?   3. **Equirectangular** : Simply consider the equirectangular projection on an infinite 2D plane. However, this assumes features are chopped at the antemeridian, otherwise a fallback would be required to one of the two methods above. My personal preference is the first approach, but I'm interested in understanding whether this is common in standard GIS systems.


title: arcgis 10, python update raster of mxd file

text: i wanted to update the raster file in a mxd file using python. for example, i created a.mxd which contains b.tif file in it. now, i obtained c.tif and want to update the a.mxd, i.e. want to delete b.tif and add c.tif, or overwrite b.tif with c.tif               import arcpy          mxd = arcpy.mapping.MapDocument(r"C:\data\a.mxd")     df = arcpy.mapping.ListDataFrames(mxd)[0]          updateLyr = arcpy.mapping.ListLayers(mxd)[0]     sourceLyr = arcpy.mapping.Layer(arcpy.Raster(r"C:\data\c.tif"))          arcpy.mapping.UpdateLayer(df, updateLyr, sourceLyr, True)          del mxd      i started as above, but i think that i am not getting most of the things correctly... plz can you help me? thanks a lot!


title: Where does MapInfo store its Coordinate Systems?

text: When choosing a projection in MapInfo, there is a large list of available coordinate systems. Where is the file that stores these and is it easily editable?


title: FME ArcSDE.... the scale is very small

text: I have a problem with ArcSde and FME. I try to import some line features to ArcSDE Oracle. When I start the script in FME, I have some errors: impossible to insert linestring in ArcSde because the scale is very small. The length of linestring is equal at 100 meters. But I have the others geometry with the same length and there is no problem to insert theses features in ArcSDE. So, i have no idea on this problem. Is it possible to clean geometry with FME ? or a correction of geometry ? (I try Geometry Filter, OGC validator, but nothing). Thanks Error: ***Line Feature not allowed in table 'CHEMN'. Migration attempted Error while generating simple line shape for Table 'CHEMIN'. Scale Factor may be too small.:SDE_ERROR_CODE(-155) Linestring or poly boundary is self-intersecting Error while setting the shape to 'sde30_nil' for feature in table 'CHEMIN'. Shape Type Not Allowed - Change Table Def.:SDE_ERROR_CODE(-35) Shapes of this entity type not allowed in this layer***


title: What tool would you use to "clip" a vector line file to a vector polygon file?

text: I am trying to select the use of a large vector line file within only the parameter of a vector polygon file to do spatial analysis. It would be the same thing as using the extract> clip tool under the analysis tools that is for two polygon files. What tool do I use that works to clip a vector line layer to a polygon layer?


title: Hypothetical Cost Saving Question: Proprietary v Open Source

text: I am writing a paper for my GIS MSc and trying to justify a move for the EPA from proprietary GIS software to open source (QGIS). Can anyone hazard an educated guess as to how much ArcGIS would cost an organisation the size of the EPA per year? Are we talking 100,000's or millions? Cheers


title: How to enable GRASS-modules in QGIS?

text: I have the latest QuantumGIS-1.9.0-master from Github with python-modules support and full GRASS-integration (GRASS-4.6.3-RC2) on Linux. When opening the GRASS-tools plugin, it only displays six tools. But there are more like 300 installed. ![Only six tools available?](http://i.stack.imgur.com/AwVyh.png) The QGIS-Wiki tells me to add _new tools to the GRASS toolbox_ menu using some xml file in:               qgis/src/plugins/grass/config/default.qgc      But where to find this file and how to update it? Isn't the `src`-directory only used prior compiling QGIS? And why are my GRASS modules not included in QGIS by default? **Update 1** : Ok, at least I found the file:               ~ $ find / -name default.qgc 2> /dev/null     /usr/local/share/qgis/grass/config/default.qgc      But the contents look good. How to trigger QGIS to load the tools? What's the point I'm missing here? **Update 2** : Okay, using Sextante-plugin seems the way to go. But this is pretty much broken. I keep getting Error-Messages like this one. I tried reinstalling, but this didn't solve it. ![Sextante: broken all over.](http://i.stack.imgur.com/3f5Dl.jpg)


title: ArcObjects IDataset.Copy causes program to crash

text: My implementation of IDataset.Copy continues to crash and I cannot figure out why. The program attempts to copy a shapefile from a temporary folder in one directory to a final folder in another directory. Here is the relevant code, where strOutVoronoiPath is a filepath containing the destination folder and g_temporary is a filepath containing the temporary folder:                       //create a shapefile workspace factory             ESRI.ArcGIS.Geodatabase.IWorkspaceFactory pSWSF = new ShapefileWorkspaceFactoryClass();                  //create a temporary IWorkspace from Shapefile IWorkspaceFactory object             ESRI.ArcGIS.Geodatabase.IWorkspace pTempWSP;             pTempWSP = pSWSF.OpenFromFile(g_tempFolder, 0);                              //create a final IWorkspace from Shapefile IWorkspaceFactory object             ESRI.ArcGIS.Geodatabase.IWorkspace pFinalWSP;             pFinalWSP = pSWSF.OpenFromFile(strOutVoronoiPath, 0);                  //create a temporary IFeatureWorkspace from Shapefile IWorkspaceFactory object             ESRI.ArcGIS.Geodatabase.IFeatureWorkspace pTempFWS;             pTempFWS = (IFeatureWorkspace)pSWSF.OpenFromFile(g_tempFolder, 0);                  //create a final IFeatureWorkspace from Shapefile IWorkspaceFactory object             ESRI.ArcGIS.Geodatabase.IFeatureWorkspace pFinalFWS;             pFinalFWS = (IFeatureWorkspace)pSWSF.OpenFromFile(strOutVoronoiPath, 0);                  //Some code goes here where I do some casting and operations             //to produce a shapefile inside the temporary folder. I now have             //a shapefile I want in temporary folder (at g_tempfolder, which is a string              //filepath) and I want to copy it over to a final folder at strOutVoronoiPath             //I have taken this shapefile and stored the data from it inside an IFeatureClass.             //called pFeatureClass. The following code attempts to take the feature class data and              //cast it into IDataset, from which I call IDataset.Copy. The program crashes              //on the final line, where I invoke the pDs.Copy() method.                   pFeatureClass.FeatureClass = pTempFWS.OpenFeatureClass(strOutVoronoiFile);              pDs = (IDataset)pFeatureClass; //cast the feature class into a dataset              pDs.Copy(strOutVoronoiFile, (IWorkspace)pFinalFWS);      I'm thinking I might have cast something wrong in the final three lines of code?


title: how to search for a Point of MultiLineString geometry in a json file using OpenLayers

text: I have a json file with large geometries(MultiLineString), here I got a Point among the json layer(baseLayer). I would like to search for a Point's LineString in the layer, from which MultiLineString that Point been selected. And how to read a json file in OpenLayers/JavaScript I can project the json file in OpenLayers, but when I click the map I am getting the Point, from that Point I would like to search for its base geometry(LineString/MultiLineString) Thank you


title: Trying to map numerous attributes on a map in MapInfo

text: Any help would be much appreciated for the following problem: Based on grouping postcode sectors together I have created 30 sales areas in the country. These sales areas have been made visible on the map using normal thematic mapping (separate colour for each area to make them easy to differentiate). For each sales area I want to show the following attributes in the centre of that sales area:   1. Product Market Share - an absolute value   2. Sales Volume - a graphical representation (e.g. size of circle)   3. Target Area - Yes or No   4. Main Competitor - Choice of 3   5. Contract Renewal Date - Choice of 3 (less than 12 months, 12-24 months and more than 24 months) Does anyone know if this is possible? My thinking so far, although untested is: If I can create a point in the centre of each sales area (not sure how to do this though!?), I could use dispersed groups as the font for this point, with 3 squares which are colour co-ordinated to show attributes 3, 4 and 5 above. If I can create another point next to the first one, I could create a pie chart showing the market share, where the size of the pie chart indicates the volume of sales. This would cover off attributes 1 and 2 above. Many thanks for any help that can be given! All the best Kev


title: Land use land cover change post- classification in ArcGIS for cross- tabulation

text: Please some one help me first I have finished the image classification process in another software package but many literary recommend Arc GIS package to process post-classification which means from which one land scape type has been changed to what. So please how can I process post-classification in Arc GIS 9.3 to look at the matrix of change or "From-to" change? Thank you in advance Alemayehu


title: KML from EPSG:4326 to EPSG:900913

text: im trying to create KML files manually (XML format) by extracting data from a postGIS database with "geography" values in a column. But i belive i need to transform the projection of each row in the database from EPSG:4326 to EPSG: 900913.. How can i do that? Code:                map = new OpenLayers.Map({     maxExtent: new OpenLayers.Bounds(-20037508.34,-20037508.34,20037508.34,20037508.34),     numZoomLevels:16,     zoom: 0,     maxResolution:156543.0339,     units:'m',     projection: "EPSG:900913",     displayProjection: new OpenLayers.Projection("EPSG:4326"),     div: "map_canvas",     controls: [                     new OpenLayers.Control.Navigation(),                     new OpenLayers.Control.PanZoomBar(),                     new OpenLayers.Control.MousePosition(),                     new OpenLayers.Control.LayerSwitcher()             ],     });               var KML = new OpenLayers.Layer.Vector("KML", {             strategies: [new OpenLayers.Strategy.Fixed()],             protocol: new OpenLayers.Protocol.HTTP({                 url: "assets/text.kml",                 format: new OpenLayers.Format.KML({                     extractStyles: false,                      extractAttributes: true                      })             })         });        * Coordinate system in database: (one "Area") / Polygon. -179.999995473772 -89.9999977368861,-89.9999954737723 -89.9999977368861,4.52622771263123e- 06 -89.9999977368861,90.0000045262277 -89.9999977368861,179.999995473772 -89.9999977368861,179.999995473772 89.9999977368861,90.0000045262277 89.9999977368861,4.52622771263123e-06 89.9999977368861,-89.9999954737723 89.9999977368861,-179.999995473772 89.9999977368861,-179.999995473772 -89.9999977368861 Current error: Uncaught Bad LineString point coordinates: -179.999995473772


title: Change text color of PopupMobile

text: I am facing the problem that with some mobile Browsers like DolphinHD and Android Browser that the browser makes the text color white, and the background is also white. So i want to change the textcolor - or the background color but there's no such thing defined in the documentation. http://help.arcgis.com/en/webapi/jav...opupmobile.htm Does anyone know how to do this?


title: Problem creating TIN from multiple contour shapefiles with ArcGIS 9.3 Spatial Analyst

text: I have about 20 shapefiles that contain contours lines, and I'm trying to convert them to one TIN using Spatial Analyst in ArcGIS 9.3. I included all of the SHP's as input to generate a TIN, but it didn't work. I have tried with hard line, soft line, and mass points options, but nothing is working. I think if I first created a single contour shapefile from the 20 then I would be able to generate a TIN. My question is: how can I create 1 contour layer from the original 20, or how can I create one TIN from multiple shapefiles?


title: ROAD COORDINATES Data (CARDIFF)

text: Does anyone know how to get road coordinates data from openstreetmaps? I know its possible I just don't know how. Thank you!!


title: QGis - Plugin libraries

text: I need to develop a couple of qgis-plugins. I made some promising first steps and completed a few easy plugins. Now its time to move to more complex plugins which need to import several tool-classes I developed. As this is quite a lot of stuff and most of my plugins need to use these libraries I wonder wether it is possible to copy them somewhere where all plugins could import them from. This needs to be a folder which is accessable for the current user. I thought about using something like "/home/username/.qgis/python/plugins/includes" but I couldn't get this working. Can anybody give me a tip? Thanks and regards Ok - hello again, on writing it occured to me - I can use "/home/username/.qgis/python" to put my classes in - I guess that's why its called like that! Thanks again


title: web map with arcGIS for server and javascript?

text: How would I go about creating a feature layer from Arcmap, and saving it as a service to arcGIS for server, which I could call from javascript API, so the web map will provide user interaction. For example, selecting polgyons/states of Australia, which will display information in a list box based on the user selected area among the web map. Being said, is it possible to just use a shapefile as a web map with the use of arcGIS for server and javascript api? I found a few sites, however, most people tend to call a shapefile upon a base layer map, such as Google maps or ESRI. I guess at this stage I am trying to understand a few things, in order to proceed further. Any help/links would be greatly appreciated. Thanks in advance. I have recieved a few suggestions/help already, but trying to complete a project with the use of ArcGIS for server and javascript api only.


title: Maki Icons to Qgis

text: Does anyone know how to import Maki icons to qgis?   Maki Icons