Importing OGR from OSGEO Issues

I am trying to run:               from osgeo import ogr      and I get: Traceback (most recent call last): File "", line 1, in from osgeo import ogr File "C:\Python27\ArcGIS10.2\lib\site-packages\osgeo__init__.py", line 21, in _gdal = swig_import_helper() File "C:\Python27\ArcGIS10.2\lib\site- packages\osgeo__init__.py", line 13, in swig_import_helper import _gdal ImportError: No module named _gdal So obviously, it's finding the osgeo mod, but what the heck is the _gdal mod? A simple search within my Python paths, shows that I have a module called gdal, but not one with the underscore. Seems like a small issue?


Publishing GeoTiffs with Geoserver - ColorPalette with Alpha

I want to publish rasterdata (Geotiffs) with GeoServer. The resulting image from WMS is shown in Grayscale, but actually it is an RGB-Image. gdalinfo:               P:\OSGeo4W\bin>P:\OSGeo4W\bin\gdalinfo  Q:\data\tif\pic.tif     Driver: GTiff/GeoTIFF     Files: Q:\data\tif\pic.tif     Size is 4523, 4522     Coordinate System is:     LOCAL_CS["DHDN / 3-degree Gauss-Kruger zone 3",         GEOGCS["DHDN",             DATUM["unknown",                 SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],             PRIMEM["Greenwich",0],             UNIT["degree",0.0174532925199433]],         AUTHORITY["EPSG","31467"],         UNIT["metre",1]]     Origin = (3532474.929505631800000,5367000.280343170300000)     Pixel Size = (0.333485342299112,-0.333485342299112)     Metadata:       AREA_OR_POINT=Area       TIFFTAG_DATETIME=2012/05/08 11:38:04       TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)       TIFFTAG_SOFTWARE=WGEO 5.0 WGEO BASIS       TIFFTAG_XRESOLUTION=200       TIFFTAG_YRESOLUTION=200     Image Structure Metadata:       COMPRESSION=DEFLATE       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  ( 3532474.930, 5367000.280)     Lower Left  ( 3532474.930, 5365492.260)     Upper Right ( 3533983.284, 5367000.280)     Lower Right ( 3533983.284, 5365492.260)     Center      ( 3533229.107, 5366246.270)     Band 1 Block=256x256 Type=Byte, ColorInterp=Palette       Overviews: 2262x2261, 1131x1131, 566x566       Mask Flags: PER_DATASET ALPHA       Overviews of mask band: 2262x2261, 1131x1131, 566x566       Color Table (RGB with 256 entries)         0: 0,0,0,255         1: 128,0,0,255         2: 0,128,0,255         3: 128,128,0,255         4: 0,0,128,255        ...       253: 255,0,255,255       254: 0,255,255,255       255: 255,255,255,255     Band 2 Block=256x256 Type=Byte, ColorInterp=Alpha       Overviews: 2262x2261, 1131x1131, 566x566      I think the problem is the Alpha-Channel. Without this Channel, the resulting image/colors are correct. I'm using the default style "raster" from Geoserver.


Using Torque with Cartodb

Is torque supposed to be part of Cartodb visualizations by default? It doesn't show up for me. All I get is simple, choropeth, category, bubble. This post and video (http://blog.cartodb.com/post/66687861735/torque-is-live-try-it-on-your- cartodb-maps-today) make it seem like there's nothing to install or setup so I'm wondering why it's not working for me. I'm just using a free account, but the website above says Cartodb has released it to all users so I don't think that should matter.


search a point of interest using GeoServer

I develop an application of webmapping using GeoServer, GeoExt and ExtJS. I want to search my points of interests (layer that exist in GeoServer) so i create a two combox where i have to put the categorie of my points of interest and the name of my point of interest(optional). My data base is in PostgreSQL. Any idea guys of how i should proced ?


Having trouble geocoding ZIP 4 only using 10.0 NA address locator

I'm having a rough time getting my list of 100 ZIP plus 4's to geocode in ArcMap 10.1. I'm using the 10.0 North America geocode service which has separate drop downs for ZIP and ZIP plus 4. I do not have a full address - only a ZIP4 which (contrary to popular belief that zips aren't tied to geography) in many urban areas tells us what side of the street that range of addresses lies. All 100 records go unmatched - so the question is - is it even feasible to geocode zip plus 4 using the free ArcMap locators? I can go to Mapquest and using my developer key am able to geocode ZIP4 right on the money every time in urban areas with exact decimal degree coordinates. Google maps don't recognize ZIP plus 4 and I am beginning to believe that it doesn't work in Arcmap either.


GeoJSON formatting for Leaflet

I am using http://GeoJSON.io as the source for my data, which is then to be displayed via Leaflet.js. I found that leaflet neads the file to include "var name = {" but GeoJSON.io doesn't like it, and doesn't load the file if it starts with this. To edit the file, I need to remove this part, edit, then add it back in. I am a beginner, so could anyone help me to avoid doing this? Can I declare the vairable some other way? This is my HTML:               <!DOCTYPE html>                        <title>Vietas</title>         <meta charset="utf-8" />         <meta name="viewport" content="width=device-width, initial-scale=1.0">                  <link rel="stylesheet" href="dist/leaflet.css" />                  <link rel="points" type="application/json" href="dati1.geojson">

                                                     And here is the current GeoJSON file: http://sapnis.com/vietas2/dati1.geojson As you see, it starts with "var" which is unacceptable per JSON syntax **I marked the question as answered, but will add the solution for future visitors of this page. Solved like this:**               <!DOCTYPE html>                   <title>Vietas</title>         <meta charset="utf-8" />         <meta name="viewport" content="width=device-width, initial-scale=1.0">                  <link rel="stylesheet" href="dist/leaflet.css" />                  <link rel="points" type="application/json" href="dati1.geojson">


How to project Google Earth KML files in ArcGIS so that they display properly?

I am sure my question has been answered in the past but I would really appreciate it if someone could help me with this... I have downloaded a KML point file from Google earth and I imported it into ArcMap. I cannot figure out what to do in order to project it correctly. My dataset is set to WGS 1984 UTM zone 34N and the KML file is in WGS 1984. The difference of the file's display in ArcMap and Google earth is not that dramatic but some points fall into the sea which is bad for the anaysis I am trying to perform. I tried to project the Google earth point file in WGS 1984 Web Mercator (Auxiliary Sphere), but it still didn't display properly... Can somebody please explain how it works with Google earth and ArcGIS projections? Thank you in advance for any response :)


Error #010067 when running spline interpolation

I am attempting to run a spline interpolation on a point dataset of 585 features in ArcGIS 9.3, using the spatial analyst toolbar. The dataset holds only latitude, longitude and abundance attribute data (all as 'double' field type), in addition to the standard 'FID' and 'Shape' fields. I have successfully run IDW and natural neighbour interpolations on the data, but whenever I attempt to run a spline interpolation, Arc returns the error '#010067: Error in executing grid expression'. The aim of this interpolation is to model fish abundance data across the North Sea and eastern English Channel, whilst producing as smooth an output as possible (hence dissatisfaction from my bosses with the natural neighbour and IDW interpolations). I have queried the error and common responses on the ESRI forums suggest either a) illegal characters are used somewhere in the filename/attributes/file path; b) the output raster name is too long or c) the shapefile has coincident points. I have tested these theories by firstly making a copy of the shapefile and saving it to 'C:\', with no illegal characters. This produced the same result. I also had a raster name of three characters - 10 under the 13 limit. The shapefile also has zero coincident points. Some are close, however, so I ran the spline on a 'thinned out' dataset, with only disparate points. Again, I received the error '#010067'. Does anyone have any suggestions? I have run a spline interpolation on a test dataset that I created, just to be sure the tool wasn't in some way erroneous, but that ran without issue.


Is it possible to route shapefiles using python and without ArcGIS, QGIS, or Pgrouting?

The question is pretty much in the title. Is there a way to route through a road line shapefile using python and some library without having to resort to converting the data into PostGIS/pgrouting or going through desktop software API such as ArcGIS and QGIS?


How to iterate through only select attributes of a single field using ModelBuilder in ArcMap 10.1

The feature class that I am iterating with ModelBuilder is a polygon FC called "Alignment" in my model. My current solution to iterating through only certain attributes of a single field does not seem to be working. My initial thought when creating this model was that I would use "Iterate Feature Selection" to choose the field that I wanted to iterate and then use the "Null Value" option to specify the attributes in the field that I wanted to skip. This might still work with a syntactic correction but I do not like that this current model is based on exclusion rather than inclusion, below is a screen shot of how I have the tool set up. ![enter image description here](http://i.stack.imgur.com/mnEXX.jpg) I would prefer to have more of an inclusion set up (specify which attributes I _do_ want included) rather than an exclusion set up (specify which attributes I _dont_ want included). Does anyone know how I can alter my model to accomplish this? ![enter image description here](http://i.stack.imgur.com/gna2b.jpg) Thanks, kkaszas


Image not loading properly with OpenLayers Map in Chrome

I have an application using OpenLayers 2.13 that load serveral maps. The first image loads and displays fine. The second dosen't load properly. See image below: ![enter image description here](http://i.stack.imgur.com/vs56N.png) I am new to OpenLayers. The application uses the same code to load and build all the maps but only the first one seem to load properly. It seems to work just fine in IE and FireFox but in Chrome it loads and displays like the above image. Any help would be greatly appreciated.


FME Server API - Retrieve Output file as a URL

I have a FME Server web workspace converting a file from SHP to KML, once it runs it outputs a kml file, I need to be able to take this file and drop into my web API google Map. Is there a way to retrieve this temporary file and display in my web map?


Road Graph Settings

I am unable to configure road graph setting. in "Plug in Menu" there is no road graph settings but it is already installed. Please tell me how to solve this problem.


Inserting new Text Element into MXD using ArcPy?

I'm using Desktop 10.2 and I tried with this script to add new text element "road 90" to a lot of mxd's, but Something wrong with the code:               import arcpy     from arcpy import env                env.workspace = r"F:\Projects\Nadav\road90\gis"     for mxdname in arcpy.ListFiles("*.mxd"):         print mxdname         mxd = arcpy.mapping.MapDocument(r"F:\Projects\Nadav\road90\gis\\" + mxdname)         for TextElement in arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT"):             TextElementName = "road 90"             if TextElementName == "TextElementName":                 true                 print 'text'         mxd.save()     del mxd      Any help would be much appreciated.


OpenLayers GeoJSON: limit significant digits to 6 when using write()?

Is there a way that OpenLayers.Format.JSON write() method to serialize coordinates using 6 decimals only ? Example: instead of http://pastie.org/4200997 to have http://pastie.org/4200994


Are there no italic fonts in Print Composer?

QGIS 1.7 doesn't seem to have an 'italic' choice when annotating a map in Print Composer, only 'regular' or 'bold' fonts, a rather surprising omission if correct. Is there some other way to create an italic font in Print Composer?


How to georeference PNG images using QGIS?

How can I georeference png images using QGIS, running on OSX? I've looked at `ExifTool`, and the plugins `GeoTag and Import Photos` + `Photo2Shape`. These tools (and other answers) as far as I can tell all use georeferenced images or images with EXIF data. I would like to load the png somewhere, and drag it around on a basemap from OpenStreeMap to position it where it should be. I know the general location, and just need to position it. (I think I know the projection, but will find out if that is the case. The images are of an area ~500m x 200m so the projection variation will hopefully not matter hugely if I am incorrect).


Calculating distance from a central point in three dimensions and visualising the result

We will shortly be collecting data from our remote systems that track ADS-B signals from Aircraft. This data will provide the lat/lon in epsg:4326 (wgs84) and a height in meters. The station location is known in EPSG:4326 as a lat/lon. We want to be able to visualize the effective range of the signals for planning additional stations. The 3D is required to visualize LOS (line of Sight) shadows due to mountains, man made obstacles, etc. How do we do the following using PostGis:   * Calculate the bearing from True North from the station to the aircraft position   * Calculate the distance over the ground between the two points   * Calculate the distance to the point factoring in the height as well. Assuming that we then store in a PostGis Database the furthermost distance for each radial, how do we do the following using PostGis:   * Create a cool looking polygon that connects the points with curves in 2D for use in OpenLayers   * Create a cool looking polygon that connects the points with curves taking in the height (z) for output to a kml file for use in Google Earth to visualize the range three dimensionally


Where is the official website of the JTS project?

I'm using the JTS API, published at http://www.vividsolutions.com/jts/JTSHome.htm However, there is another JTS API, published at http://sourceforge.net/projects/jts-topo-suite/ The former is version 1.8, from 2006. The latter is version 1.13, from 2013. I don't understand. What is the right version of JTS to use? If the right version is at www.vividsolutions.com, then is this project alive? There are no new features there since 2006.


Is there a way to Query what's "NOT THERE"?

If you have land tenure mapping, and you're only interested in the areas for which there are **no polygons** in a region, is there a way to query these areas out and create a shapefile from them. A tangible example using Mining Claim Mapping: Various forms of land tenure are represented for an area. These areas are captured with polygon geometry. Areas include Mining Claims, and various forms of Patented Mining Lands (ie - Surface rights only, surface & mineral rights, 11 year lease, etc, etc). In this example, If I'm **ONLY** interested in the areas that are **not** captured (which would be _Crown Land_ ), how do I go about extracting that data (data which is "NOT THERE") without manually digitizing? Is there a way to do this with a _select by location_ query? The only thing I can come up with is to use the _bpoly_ command in AutoCAD, but this would be a bit time consuming and would also lend itself to the possibility of oversights/errors on the part of the user.


Difference between ArcObjects SDK for Java and ArcGIS Runtime SDK for Java and ArcGIS Engine

What is the difference between ArcObjects SDK for Java and ArcGIS Runtime SDK for Java and ArcGIS Engine. I have ArcGIS Desktop 10.1 on my laptop. What do I need to create ArcGIS desktop applications (arcmap and standalone) in Java? Do I need to purchase a license for creating applications? License is $1500, I can't afford it. I am confused with the terminology.


Placemark name not shown after converting from .shp to kml

I've converted a bunch of .shp files to .kml using ogr2ogr and I've noticed that the placemark names in Google Earth are not reflected in the newly made kml files. I have also used ogr2gui and found that the placemark names do appear when made with them. With that said, ogr2ogr seems to be the better tool for automatic batch conversion of the files, and hence my desire to stick with it. Edit: I have been told by Inventis (developers of Ogr2Gui) that by default the shp's column headers are replace by NameField and Description. Then, the driver fills those two tags with the attributes of the DBF. In this case, could this be a driver issue? should I be building LIBKML and using that instead of the default KML? I cannot use the 'NameField=' and "Description=" since the header names vary with each shp file. Can someone please help? Cheers!


QGIS-Processing GRASS-tools don't work

When I try to use the processing grass tools like r.series I always get the following messege. I deinstalled QGis and reinstalled the program totally fresh without success. ![enter image description here](http://i.stack.imgur.com/bKIzK.jpg)


Geoserver download for mac osx

I've been trying to set up geoserver on my mac (osx mavericks). However, at the stable GEoserver downloads page (http://geoserver.org/display/GEOS/Stable), the download link for the mac is disabled. In addition, if I download a previous release for a mac and install it; when I launch Geoserver, I get an alert that the application is damaged and should be deleted. Does anyone know of a stable build for a mac that they can point me to? Thanks, Tyler


QGIS locked into fullscreen mode on Ubuntu 12.04

No crash in this case, but I made the mistake of toggling to fullscreen mode on my netbook running Ubuntu 12.04. It's now impossible to escape this "No Menus" view, and it seems to have disabled the popout unity launcher as well. Steps I've already taken:   * Alt-Tab to move to a different program, then back - QGIS remains in fullscreen   * Alt-V, Ctrl-F and F11 get no response   * I've unistalled QGIS, deleted the /.qgis/ profile data folder, and reinstalled. Still launches fullscreen To add to the weirdness, it now launches - still in fullscreen - into a different workspace than the active one. I realize this has the sound of a Unity problem, but it's just this program that has a malfunctioning UI. Thanks to all the crew; this in no way changes my opinion of QGIS as an essential tool.


How to georeference an areal photograph given known coordinates of the four corners and center of the image

First, let me say that I am new to this field. I have an existing image/video viewing application. My primary goal is to implement an efficient method to determine the latitude and longitude of any pixel in an image or video frame. I know the latitude and longitude of the four corners and the center of the image/frame. I also know the projection system that the coordinates are based in. The images in question are not maps. They are areal photographs. I cannot make assumptions about the viewing angle or anything like that. I'm not looking to use any existing application or tool. I need to implement the necessary functionality in my own application. I am open to using a third party library if that's the best way, but my application is proprietary. So licensing is a concern. I've been looking into this for a while now. I've read up on world files and using them for rectification. But I don't know how to create one for my image. It seems like the tutorials that I have seen make assumptions for simplicity sake that I cannot make. I'm basically looking for a description of the necessary math to perform the required analysis and/or a good library I could use that is compatible with proprietary software. My existing application is written in c++ and c#. I apologize if any of terminology is incorrect. Any advice would be appreciated.


GeoJSON generator with OpenLayers support

Is there a service that allows me to draw on an OpenLayers map, and see the resulting GeoJSON?


Converting a .jpg map into a shapefile

I have a .jpg map of a set of U.S. counties that were affected by a tornado. The map shades where the tornado passed through, resulting in some counties only partially experiencing the tornado. I would like to calculate the share of each county's boundary that was exposed to the tornado. Naturally, this is straight forward for the counties that were either entirely exposed or not at all exposed. However, I am not sure what to do for those that were partially exposed. I came across this thread: How to convert image of map into vector format? I cannot understand it, however. For starters, does Paint count as an image editor. Second, what's the deal with GRASS? Is that the best way to go? Is that in lieu of ArcGIS 10.1? In any event, I do have access to ArcMap in ArcGIS 10.1. I am trying to follow the link: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//009t000000mq000000 However, I am very lost starting around step 6. I do not know if I'm adding the .jpg correctly or assigning the points correctly. I am just extremely lost. In the end, I would like an excel spreadsheet of the share of each county exposed to the tornado. Thank you for any help!! To put the calculation that I'd like to do into context, here's an exmaple state-level map: http://www.marc.org/emergency/images/tornado_risk_map.jpg In this example, suppose I'd like to estimate the share of each county (or state) within the red.


Is there a way to use the file dialog in Model Builder?

Is there a way to have the file dialog come up and let a user choose a file name and location for output? I have the table to Excel script as an geoprocessing tool in the Flex Viewer and would like to have the user select the file name and location rather than use hard coded ones. Any ideas are appreciated. thanks


How to write binary ArcGRIDs outside of ArcGIS in linux?

I've been looking for a way to create\write binary formatted ArcGRID images in Linux without much luck. It seems the GDAL project only supports reading these files, but there is not write support. My major problem with this is that I have raster images being created in custom C++ programs and I'm using ArcGIS Server to display them. Most operations in ArcGIS Server regarding raster processing via geoprocessing services seems to require the data to be in the binary ArcGRID format (or automatically converts to this format). Converting from another format to the ArcGRID format on the ArcGIS Server is just too slow for web apps.


Loaded shapefile into SDE without spatial proj creates problem

Via a Python script - various shapefiles were loaded into SDE with no spatial reference set. This causes a strange problem - when you view the featureclass in ArcCatalog (geography view) - it shows the following dialog: ![enter image description here](http://i.stack.imgur.com/6mvC8.png) When you look at the Resolution tab of the featureclass it has an invalid resolution value: ![enter image description here](http://i.stack.imgur.com/EUjbX.png) It should normally show a small value - like as 0.000000001. The resolution is effecively a mulitplier or divider that ESRI GDB storage uses to convert real numbers of the vertex coordinates into integers for storage, so in this context, ‘0’ is invalid as you can’t divide by zero. Is there any way to change this value?, as it is certainly not editable in the dialog and Calcualate Default XY Tolerance in ArcToolbox does not work. Can it be altered by ArcObjects - or could that break things.


Why is every SELECT via a geometry function null?

I'm using Spatialite (2.3.1, binary from the website) on Windows XP via Java ("Xerial Driver"). SQLite loads fine, also Spatialite seems to be loading successfully. I can even do SELECT-statements. But I can't do SELECT's on geometries. E.g.               SELECT link_id, SRID(GEOMETRY), AsText(GEOMETRY) FROM streets1      returns only null columns except for link_id. How can that be? By the way, I imported this data originally from MapInfo files using ogr2ogr:               ogr2ogr -clipsrc 8.6688 50.1018 8.6794 50.1082 -f "SQLite" streets1.SQLite ~/Desktop/G2AM10101EG2000MAADT/Streets1.TAB      Also I tried:               ogr2ogr -clipsrc 8.6688 50.1018 8.6794 50.1082 -f "SQLite" streets1.SQLite Streets1.TAB nlt GEOMETRY      This also gives the same result... Thanks, Philip


Time estimates on building overviews for mosaic datasets

I'm trying to create some time estimates for imagery caching in ArcGIS 10. Does anyone know of any forumlas or ratios to calculate the time it takes to generate overviews for a mosaic dataset? Does it have anything to do with the area, resolution, file size etc of the actual imagery, or is there even a way to calculate a time estimate for overview generation? Thanks **UPDATE:** Does anyone have experience in building overviews on massive mosaic dataset (upwards of 9,000 images) in arcGIS 10? I'm wondering if there is a way to speed up the process ... right now it looks like arcGIS is only using 1 core to build the overviews. I have a 12 core machine. Obviously things would go much faster if we can utilize more cores. My original estimates are blown out of the water ... I'm looking at over 3 weeks to complete this process. Any ideas?


ESRI Silverlight API - KmlLayer with a Schema?

I have a KML file I'd like to load with the KmlLayer in my Silverlight application. The problem is that the markup seems more complicated than the class wants to handle. I have a element at the top of the Document that adds custom elements () and for some reason changes the element name to a custom one. I've been searching for a couple days but no one seems to have covered this issue. Has anyone here come across a option or somesuch you can set so that the KmlLayer notices and adheres to the elements customizations?


Importing legacy shapefiles and attribute information to geodatabase

I am trying to take older sewer/water system shapefile data with limited, but important attribute information and importing it in to a feature class in a geodatabase, while keeping the attribute information. These older files do not have similar field names as the geodatabase in which to import the attribute information. Any advice in which to complete this task is greatly appreciated. The solution does not neccesarily have to be solved programmatically.


How to generate a Shape File in DotSpatial

I have a .csv file of coordinates such as the following:               -90.19, 38.62, St Louis     -122.42, 37.82, San Francisco     -74.04, 40.68, New York      I want to create a shape file using DotSpatial that has these points, so that I can overlay them on a shape file of the United States. Which functions in DotSpatial would I provide these points to, to generate a .shp file?


QGIS- "On the fly" projection ON by default?

In QGIS 1.6, can users set "on the fly projection" to be checked ON by default (at application launch)? I know how to turn it on, but some users find they use it more often than not, and would like it to be enabled by default. I have read this thread and it seems like there was discussion to add an "ON by default" option, but I could not find it in the Options settings. Is there something I am missing or has this function not been implemented?


Intersection Type Data

Would anyone have any input as to where I could find a table/shapefile/feature class that describes the intersection type of roads? I've looked into TIGER roads, StreetMap premium, OpenStreetMap tables. Specifically looking for "street-light", "4 way stop", etc. If such a database even exists.


How to iterate field values and use output in ArcGIS ModelBuilder?

I am using ArcGIS 10.2 Model Builder with the Field Values iterator to select unique field values. I then want to forward these values to the Add Field tool, to create a new field for each unique value. The problem I have is that the iterator's output value appears to have a temporary starter value of "1" by default, which is then replaced by the actual values once the model runs. However, the issue is that the Add Field tool throws an error before the model even starts, because it cannot use numbers as inputs for field names, so it remains out of action once the "real" values become available. Ideas for how to flush-out or replace this starting value of "1" would be much appreciated. FYI - the Iterate Field Values tool is correctly set to "String" data type, yet the temporary starter / filler value remains "1". I've attached a picture of my model and submodel. Any help would be appreciated. UPDATE: I was unable to get it working in Model Builder, it is presumably a bug. But I was successful in getting it working via a python script along the lines of the answer listed below. ![Sub-model called "sifter"](http://i.stack.imgur.com/NeEmo.png) ![Model](http://i.stack.imgur.com/Pfj9u.png)


Non-admin GeoServer users with advanced privileges

I have user, who has full access to only their workspace. How can I allow him to:   * change their password   * set password protection for reading/writing layers from his workspace   * create own service rules for own workspace/layers


How to add columns with x/y data to Spatialite Layer?

I´m new here and I am also quite new at working with QGIS. My Spatialite Database looks like this (arch_locations.sqlite)               Fields:     pkuid     short_ID     name     geometry      I´d also like to have two additional columns in that table in QGIS (in the Spatialite Database) where it displays the X and Y coordinates. In the "geometry" field I only see a very long combination of numbers and letters. While looking for a way to adding those columns I found this short guide: http://pvanb.wordpress.com/2012/02/21/extract-x-and-y-coordinates-from-a- spatial-table-in-spatialite/ Using the Plugin "QSpatialite" and modifying the following code passage I managed to create another Table in QGIS with two columns containg the x and y value. I adapted this code:               CREATE TABLE 'Test' AS          SELECT PKUID, X(Geometry) as X, Y(Geometry) as Y          FROM Points      I did not even know where to enter the code until I installed the plugin ;) Though, what I really want is to add two additional columns in the original table and to display the x/y values in those. How do I do this?


Mapnik, TileCache, Apache2 - Setup&Debug

I'm trying to setup TileCache with Mapnik working on an Apache Server. Right now I don't see any tiles in my OpenLayers Framework in the webbrowser. I really ask myself where the problem is hidden. a) Does anyone have a good debug strategie for working with these 3 components? b) Anyone working with that and can give a little working sample setup. Maybe my confs are wrong or incomplete. Thanks!


problems about setting pgrouting with QGIS under windows

I followed this website：A Beginner’s Guide to pgRouting and now facing some problems. I finished steps before "Creating a routable road network", and open QGIS 1.8, click the PostGIS Manager plugin, here's the figure：![enter image description here](http://i.stack.imgur.com/fE9We.png), the first question is that how to enable PostGIS support？ And if I click the Data button trying to load data from shapefile, here's the figure：![enter image description here](http://i.stack.imgur.com/B2GhX.png) And when I try either way：back to QGIS menu->click the Database->shp to postgis->key in my password->click connect, then it shows：![enter image description here](http://i.stack.imgur.com/FiN1V.png) How to solve these problem and correctly use pgrouting？ ps.I work under windows7 64-bit, and have QGIS1.8, postgres8.4, pgadmin3, postgis2.0 installed.


how to connect points to create a polygon boundary in QGIS

I have 31 points in a shapefile, which delimit a field (it's an irregular shape). I want to connect those points so they create the boundary of a polygon. Does anybody know a kind of "connect the dots" tool?


How to split line layer at polygons' boundary of a polygon layer?

I have a line layer and a polygon layer in PostgreSQL 9.2, PostGIS 2.0. How can I split every line each time it crosses a polygon's boundary? ![enter image description here](http://i.stack.imgur.com/ksl9B.png)


Problem with compilation

I have tried to compile Quantum GIS by it's source code on Windows 7 (32 bit).i have follow it's installation guide.i have download and install all it's dependencies . i have run it on C-Make but it gives error that QT_QMAKE_EXECUTABLE not found. so kindly tell me how to compile it successfully.


projection/pyproj puzzle, and understanding SRS format

I have a shapefile that has a .prj file, with the following wkt contents:               PROJCS["NAD_1927_StatePlane_New_York_West_FIPS_3103",GEOGCS["GCS_North_American_1927",DATUM["D_North_American_1927",SPHEROID["Clarke_1866",6378206.4,294.9786982]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-78.58333333333333],PARAMETER["Scale_Factor",0.9999375],PARAMETER["Latitude_Of_Origin",40.0],UNIT["Foot_US",0.3048006096012192]]      I am using the following python code to convert that to a pyproj projection:               prjText = open( prjPath, 'r').read()     srs = osr.SpatialReference()     if ( srs.ImportFromWkt( prjText ) ):         return ( False, "error importing .prj information from " + prjPath )     inProjection = pyproj.Proj( srs.ExportToProj4() )      I then create another projection:               outProjection = pyproj.Proj( "+proj=latlong +datum=WGS84" )      and then go through each point `p` from the shapefile and convert it:               print "before transform: " + str( p[ 0 ] ) + ", " + str( p[ 1 ] )     p = pyproj.transform( inProjection, outProjection, p[ 0 ], p[ 1 ] )     print "    -> after transform: " + str( p[ 0 ] ) + ", " + str( p[ 1 ] )      when I print the srs.ExportToProj4() text, it gives:               +proj=tmerc +lat_0=40 +lon_0=-78.58333333333333 +k=0.9999375 +x_0=152400.3048006096 +y_0=0 +ellps=clrk66 +to_meter=0.3048006096012192 +no_defs      Furthermore, using another program that I have, I import the same shapefile without using the .prj file. To do this, I specify that the projection is "State Plane", and specifically that it is "New York West NAD27 (ftUS)". This other program then determines that that means the same as epsg:32017, and indeed when I look up epsg:32017 I see that its srs string matches the one above. Thus, it seems like my inProjection variable (above) should be set up correctly. However, when I run the code above to transform each of the points, I get, for the first point:               before transform: 759232.003438, 1149854.52147         -> after transform: -70.1071453809, 50.036825617      **The puzzle:** That lon/lat, approximately (-70, 50), is not in the right place. By importing the file into either the above-mentioned program, or into ArcGIS Explorer Desktop, I see the point (correctly) in northern New York, just south of Lake Ontario. But (-70, 50) is up in central Quebec. So it's off by like 1000km to the north and east. So the questions are: (1) Does it seems likely something is going wrong in pyproj (I'm using version 1.8.9)? (since according to these other programs my srs string is correct, so everything is good up to the point I use pyproj) (2) Based on the wkt and srs string shown above, what lon/lat should I expect the the shapefile point (759232.003438, 1149854.52147) to correspond to? Is that point the x/y in meters relative to the lon_0 and lat_0 of the srs string? If that were the case, I would indeed expect a point somewhere up in Quebec, and not in northern New York. However, since both of these other programs are correctly mapping the point to northern New York, that must not be the correct interpretation of the file coordinates with respect to the file's projection? Thanks for any help.


Add row and column name to square grid feature in ArcGIS

I have a square grid feature class of 250 x 250 tiles. This will serve as a tile index for a a project. They currently have no identifying attributes other than a unique ID. I need to add a field to rename them by their row and column, e.g. "R001_C001" or something. Is there a way to do this easily? I'm okay with python.


Check if layer is a QGSVectorLayer from Python Plugin within QGIS

I am trying to extend some concepts in the QGIS Workshop tutorial. Basically, I would like to check if the layer is of type "QGSVectorLayer". In the QGIS Python console, I am able to type               canvas = qgis.utils.iface.mapCanvas()     cLayer = canvas.currentLayer()     type(cLayer).__name__      with no problem. However, when attempting the same from the plugin, I get the following error message: > if type(self.cLayer). **name** == "QgsVectorLayer": NameError: global name > 'cLayer' is not defined The code is as follows:               def handleLayerChange(self, layer):             self.cLayer = self.canvas.currentLayer()             if self.cLayer:                 if type(self.cLayer).__name__ == "QgsVectorLayer":                     self.provider = self.cLayer.dataProvider()


How to get a map surround as an IElement in ArcObjects?

In an MXD file, I have a page layout and a map with a north arrow, and I need to move it in the page layout programmatically using ArcObjects Java SDK. I managed to get the north arrow map surround object:               for(int i = 0; i < dataFrame.getMapSurroundCount(); i++) {             IMapSurround mapSurround = dataFrame.getMapSurround(i);             if(mapSurround instanceof INorthArrow) {                 log.info("North");                 // ...             }             if(mapSurround instanceof IScaleBar) {                 log.info("scale bar");                  }         }      Now, I need to get my north arrow as an IElement, in order to apply an Envelope as a geometry, that will cause the north arrow to be moved within the page layout. INorthArrow cannot be cast to IElement (ever tried). Any idea to get an IElement?


How to save python ogr2ogr shapefile to memory, and stream to web client as zip

I'm trying to accomplish something similar to the site http://ogre.adc4gis.com/ I want to be able to manipulate a geojson file, and send it to the client as a shapefile, in a single zip file. I'm writing the program in python, and I see that the site mentioned above is using Javascript. I've already been able to save the json as a shapefile with the following command:               ogr2ogr.main(["", "-f", "ESRI Shapefile", "fishnet", "fishnet.json"])      This saves our geojson file as 4 separate files in the folder `fishnet`. I'm wondering how I can avoid saving files to the harddisk and just serve the zip file to the client in memory. ## Update So far I've zipped it. Still yet to stream to the client, but I think it should be easy enough to find some resources on `StringIO`. Basically, what I've done so far is to go ahead and save the files to disk, get the data from them, convert them, and then delete them (after they've been streamed):               import ogr2ogr, shutil, uuid     name_of_files = 'custom_name_here'     # random directory name     # http://stackoverflow.com/questions/10501247/best-way-to-generate-random-file-names-in-python     directory_for_files = str(uuid.uuid4())     # -nln gives a custom name to files in folder     # string without .shp gets saved to directory as four files     ogr2ogr.main(['','-f','ESRI Shapefile','-nln', name_of_files, directory_for_files, 'test.json'])     # create zip from directory with same name     shutil.make_archive(directory_for_files, "zip", directory_for_files)     # remove directory     shutil.rmtree(directory_for_files)


Producing NDVI maps in qGIS

I'm trying to produce a simple NDVI map in qGIS using a Macbook pro (2011) from landsat raster images. I've tried using the 'Measuring Vegetation' plugin, but for some reason that goes absolutely haywire as soon as I load it up, and I have to uninstall the plugin to use the system at all afterwards. So I've tried doing it in the Raster Calculator plugin using exactly the same foruma as works perfectly in Grass GIS, namely Float(band4 - band3)/ Float (band4 + band3) I've tried it without float, with and without spaces between commands, with quotation marks around raster files.   All that is ever produced is a map that is uniform in tone (ie contains no information). What am I doing wrong? Cheers, Daniel


How do I set the working directory in QGIS?

I'm learning QGIS and lovin' it. But one thing that annooys me a little is is that each time I want to save data, the default directory is the QGIS folder under Program Files. How can I set this, globally, for all my data?


Does QGIS support internationalization?

Does QGIS support multiple languages? I have some shapefile with Arabic names in its attributes, and they appear as: `????` If I want to see the names in Arabic, do I need to purchase multi-language supporting GIS software or is there any other option? I am using QGIS 1.8.0.


Python Stand Alone ArcGIS toolbox output feature class becomes input feature class within same tool?

Developing a script in toolbox.   1. I have a table that I use to make X Y layer (input 1)   2. I am saving that created layer as a point feature class. (output 1)   3. I then want to run kriging on that newly created point feature class, based on one of the columns. So my output feature class becomes input feature class.   4. As you may see below at the "obtained from" part the field that will be used for kriging is not appearing yet, since the tool didn't created yet the point feature class.


Error Opening QGIS Lisboa Quantum Msys

![Error Image attached.](http://i.stack.imgur.com/Ho2fe.png) I am deploying QGIS in our environment. I found the silent switch (/S) to be working, it installs the 4 default components, Grass, MSys, GIS Desktop & GIS browser. The trouble begins when a user attempts to launch MSys & Grass GIS. Msys seems to have issues with naming space. The icon itself calls a batch file (MSYS.bat) located under C:\Program Files (x86)\Quantum GIS Lisboa\apps\msys but throws an error (refer to attached pic). Running MSYS.bat within the folder throws the same error. If I move MSYS folder to desktop then I can run the batch file. I think that the problem is with spaces in location address. If I eliminate the spaces it runs for example, C:\Users\Desktop\Quantum GIS Lisboa\apps\msys\msys.bat runs without problems. Is the programming team behind this software aware of the issue? If yes then what's the fix. Something similar goes for GRASS GIS 6.4.3RC2. Running the desktop/start menu shortcut fails with the attached error (the image on right). What seems to be the problem? Why will the program not run/function straight out of box? We have Windows 7 SP1 x64 bit operating systems in our environment. Any help will be highly appreciated.


Antipodal distance (or polygon fetch or polygon diameter) for concave polygons

I am working on some class examples in Python implemented within ArcMap to calculate the antipodal distance within a polygon. This is fairly routine for convex polygons, however, for concave polygons, I wish to exclude solutions (formed by a ray connecting the boundary points), that aren't completely within the polygon and not on the polygon boundary or intersecting it. Have I interpreted the definition wrong or is this beast go by another name. Consider these two polygons pnts = [ [0,0], [0,1], [1,4],[3,5],[5,4],[4,1],[0,0] ] # a closed-loop convex pnts = [ [0,0], [2,1], [1,4],[3,5],[5,4],[4,1],[0,0] ] # a closed-loop concave polygon In my interpretation, point 0,0 should not have an antipodal distance associated with it since the vector connecting it with the other points is either self intersecting the polygon or is on the polygon boundary. If anyone has any clarification on the definition or potential solutions, I would appreciate it. A visual of the convex polygon and the desired lines (shown in red) is enclosed (example vectors from point 0 are only shown). ![Interior Antipodal Distance example](http://i.stack.imgur.com/jqNb0.png) In the convex example, the first point has no antipodal vectors, however, the second point does. ![Concave Antipodal Example](http://i.stack.imgur.com/ujTGs.png) EDIT I have had some success searching using "polygon fetch" or "polygon diameter" on the web, I suspect that this is what I am after.


QGIS2.0 on Ubuntu 12.04 cannot see the hints when hover over

I am running QGIS2.0 on Ubuntu server 12.04. After installing some libraries for Grass and Orfeo toolbox (I don't remember exactly which ones), I cannot see the hints when i hover over a menu bar or a layer. I can see only the colored frame of the hint without the text. Has anyone experienced something similar? Thanks!


How to implement feature search in OpenLayers with Geoserver?

I'm looking to build an application similar to this ArcGIS Server JS example which does a find task against a layer in the map. However i want to do this with OpenLayers, Geoserver and PostGIS. Some pointers would be great. Regards Med


Coordinate System correctly renders both decimal Long/Lat and some number in the billions

I am new to ArcObjects and I have been creating an application which will (among other things) be reading data in from a sister program and plotting it on the map.I have found out that I can enter the coordinates as both decimal long/lat, or by using a pair of numbers close to a billion. Both of these methods will draw a point somewhere on the map. My question is: why? How is it possible that this works, and which one is correct?


Combine shape in folder for on/off

Is it possible to organize some shapes into folders to turn on/off all the shape togheter? Is it useful for the shape that you can punt at the top/bottom of a project.


How do you determine if feature is editable?

I am using the FocusMap.FeatureSelection to return all the currently selected features. I want to process only the currently editable features and leave the remainder alone. How can I quickly programmatically tell which features are editable? I know of the IEditLayers.IsEditable method, but think that might be a little heavy handed to try and map from a feature back to the featurelayer and then test that for all selected features. Is there another way, or is this all I have available to me?


Adding points with defined intervals in QGIS

I'm trying to make a grid with points of defined intervals (e.g. 5 m). How can I achieve this with QGIS?


R program and shape files

I have an R program that sketched a country. This is done by reading a shape file that contains the coordinate of a country (latitude and longitude).   1. How do I export the data from shape file into Excel? to understand the structure of the data (it is a bit hard to visually see the data structure in R)   2. How do I get R program to plot different colors for each area (postcode)? Thanks for the help.


Permission Issues on ArcGIS Server with SOC machines

I am trying to figure out some permission issues for ArcGIS Server. I have 2 SOC machines that are associated with SOM. I noticed this problem when I added these SOC machines to SOM. This service that I have published is a GP service. Appreciate any ideas. Thanks Jay First Error: ~~~~~~~~ WARNING Can not open directory c:\arcgisserver\arcgisoutput. The system cannot find the path specified. Probable cause: The network path is inaccessible to ArcGIS Server Object Container account. 2011-11-08 11:14:15 SDM.ServerDirectoryManager cwf414 Second Error: ~~~~~~~~~~~ ERROR Can not open file D:\Projects\GPServices.mxd. The system cannot find the path specified. Probable cause: The file is inaccessible to Server. 2011-11-08 11:07:31 PRJ/GPServices.MapServer cwf415 ERROR Server Object instance creation failed on machine cwf415. Can not open file D:\Projects\GPServices.mxd. The system cannot find the path specified. Probable cause: The file is inaccessible to Server. 2011-11-08 11:07:31 PRJ/GPServices.MapServer cwf415 ERROR Server Object instance creation failed on machine cwf415. Can not open file D:\Projects\GPServices.mxd. The system cannot find the path specified. Probable cause: The file is inaccessible to Server. 2011-11-08 11:07:31 PRJ/GPServices.MapServer cwf415


How to improve creating contours for super-high unnatural mountains?

in a current project I am working on we are not mapping the actual height of a landscape rather mapping the interest which specific stakeholders have in parts of the region. So, what we did was the following: We invented an algorithm to measure interest. In the next step, we mapped these interests to areas, buildings etc. Then, we converted this information to a layer of hexagons which, in the following, were converted to points ("of interest"). These point were converted to a raster to be finally converted to contours. You can follow this process in these screenshots: ![interest as polygons](http://i.stack.imgur.com/BVNWh.png) ![interest as hexagons](http://i.stack.imgur.com/OsB0X.png) This was accomplished with the help of http://anitagraser.com/2012/03/04/mapping-density-with-hexagonal-grids/. ![interest in points](http://i.stack.imgur.com/Ybqmb.png) ![interest as raster](http://i.stack.imgur.com/BgG21.png) Here, we used the interpolation extension of QGIS. ![interest as contours](http://i.stack.imgur.com/Zj3dX.png) This was done by using GRASS with r.contour.step . I tried to achieve this output with the contour plugin but the result was not satisfactory and sometimes the plugin finished with an error. I guess we just had too much points (round about 140 000 -- hexagons have a width of 50 meters...) When areas of interest overlap, they just add up. So the highest peek you can see has got a value of 74 000 and the lowest hast got 1 200. Nevertheless, the result is not bad at all but because of the extreme difference of "height" the peeks are not processed very well (as one can see). Does someone have an idea of maybe "normalizing" the heights without "cutting" the "mountains". I already played around with the settings of the raster layer changing the Min Max values -- the result above is the best I could get, yet. Thanks for your help, Niklaas


How do you export to georeferenced PDF from QGIS 2.0?

When I use the print composer in QGIS, there is no option to export to georeferenced PDF using GDAL. Do I need to export from somewhere else or how do I go about exporting to georeferenced PDF? I am using the newest version of QGIS. * * * Related questions for older versions: Is it possible to export GeoPDF from QGIS?


How to control the launch page of mobile project?

I have an ArcGIS mobile Window project to upgrade from 10.0 to 10.1. there are minor code changes. The project has 3 custom tasks in addition to view map and one capability. The project used to open at the page 'Task List' in 10.0 but now in 10.1, it is opening at 'Choose Project to open' even if there is only one project. I checked in project initialization in Extension, the code is almost same as that of 10.0 and I tried changing too but have not had any success so far. Noticed- default app behaves the same way. The launch is run by ArcGIS.lnk (a file different than ArcGIS Mobile 10.0 and not sure if this behavior is controlled there or in MPC) So right now it is opening here   ![ So right now it is opening here](http://i.stack.imgur.com/s5dS0.png) . but I want it to start it to next page (task list) If anyone has any clue about this? Thanks


Join intersecting lines with PostGIS

This is something that I'm almost ashamed to ask but I can't seem to get it to work for the life of me. I have a road layer with segments, each segment has a `Road ID` and a segment `type`. I would like to join all the segments together, for each `Road ID` into one linestring but only when they are the same `type` and are **touching** (all the lines are snaped together). ![enter image description here](http://i.stack.imgur.com/vymc7.png)               Road ID - Type        1       L        1       L        1       T        1       L


using GDAL/python to stack georeferenced images of different sizes

I am in the process of porting a code I wrote in IDL (interactive data language) to python but am running into a bit of a problem that I am hoping someone can help me with. The code goes like this:   1. take individual classified Landsat geotiffs (say there are N individual 1-band files per scene, each representing a different day) and further reduce these images to three binary-themed 1-band images (water and not water, land and not land, water/land and not water/land). This will be done by reading the rasters as matrices and replacing values. ** I don't actually need to have these images, so I can save them as memory or just keep them as numpy ndarrays to move to the next step   2. stack these images/arrays to produce 3 different (1 for each 'element') N-band stacks (or a 3-dimensional array-- (samples, lines, N)) for each scene   3. total the stacks to get a total number of water/land/water&land observations per pixel (produces one 1-band total image for each scene)    4. other stuff The problem I am running into is when I get to the stacking, as the individual images for each scene vary in size, although they mostly overlap with each other. I originally used an ENVI layer-stacking routine that takes the N different-sized 1-band images for each scene and stacks them into an N-band image with an extent that encompasses all of the images' extents, and then reading the resulting rasters in as 3-d arrays to do the totals. I would like to do something similar with gdal/python but am not sure how to go about doing so. I was thinking I would implement gdal capabilities of geotiffs by using the geotransform info of the images to somehow find the inclusive extent, possibly padding the edges of the images with 0's so they are all the same size, stacking these images/3-d arrays so that they are correctly aligned, then computing the totals. Hopefully there is an easier way, as I'm not sure how to pull that off. Does anyone have any suggestions or ideas as to what would be the most efficient way (or, any way really), to do what I need to do? I'm open to anything. Thanks so much, Maggie


What format is this GPS tracking string?

I am receiving tracking data from a device and I am not sure what is the format/standard of this string in order to properly process the variables, this is a sample string received: `#355488020131816##1#0000#AUT#01#2260012b9cbb47#2112.703000,E,4545.754120,N,0.07,0.00#031213#184625.000##` I am not sure about the variables 2,3,4,5, (`##1#0000#AUT#01#`) the others seem obvious, position, speed, bearing, date,time. The device appears to be `GT02A`


What projection is this and how to use it in a GIS?

I have a request to do a map in non-standard projection: ![non-standar projection](http://i.stack.imgur.com/C1bQP.jpg)   (source: http://www4.agr.gc.ca/AAFC-AAC/display- afficher.do?id=1263246127905&lang=eng)   I couldn't manage to find it my materials. As well I wasn't able to apply such projection in ArcGIS.      Do you know what projection is it? And how to apply it in ArcGIS or Manifold or any open source GIS?


Percent of change?

I'm a GIS student and my task is to map, at the county level, the growth in elderly population between 2000 and 2010 (as a proportion of total population) for Louisiana and to compare the pattern of growth to growth in the elderly population at the state level across the nation. So, two maps—one of the entire U.S., state level, and the other for Louisiana, county level. The variable I'm trying to map is: Change in Elderly = Elderly2010/Total_Pop2010 - Elderly2000/Total_Pop2000 I'm trying to map the percent of change. How can I go about doing this? Any help would be great! Thank you!


R - Occasional error when using alphahull package areaahull() function

I've recently started using the fantastic **alphahull** R package, which works beautifully most of the time but occasionally I am getting the following error when running the **areaahull()** function:               Error in row:rownew : NA/NaN argument      For most of my data the function works as expected, but with some data and for certain values of alpha, the above fatal error occurs. It does not seem to depend on the complexity of the data, sometimes the error occurs with a few dozen data points, while it will run fine with a dataset of a few thousand points. The code is very straightforward, simply by changing the value of alpha for a given dataset I can generate the error.               library(alphahull)          lonlat <- read.csv("areaahull_error.csv")          alpha <- 0.2          alphahull <- ahull(lonlat, alpha = alpha)     areaahull(alphahull)      I'd be interested to know if anyone else has seen this error and if there's anything to be done about it, other than trying to catch it, or tweak alpha to avoid it when it occurs.


shapefile terminology for key-value pair

In an .osm file, there are `<TAG ... />` elements like `<way visible="true" id="-4409" timestamp="1970-01-01T00:00:00Z" version="1"> <nd ref="-21845"/> <tag k="fcc" v="A40"/> <tag k="len_miles" v="0.037850000000000002"/> <tag k="et_id" v="4407"/>` The "k" and "v" in these refer to key - value pairs. What are the equivalent things called in shapefiles? Are they called (attribute, value) pairs? I see they're visible in the Attribute Table in QGIS and the k's correspond to columns and the v's are what is in the cells of that table.


Cartodb error: ETHON: performed EASY url= response_code=0 return_code=couldnt_connect total_time=0.000113

I have installed cartodb (master) on local ubuntu and I am able to import data / add new tables etc. When I click on VISUALIZE link and fill in the details, press ok, the window waits there for ever. I see an error in logs `(400 bad request)`               21:37:47 web.1     | (0.000947s) INSERT INTO "visualizations" ("id", "name", "map_id", "active_layer_id", "type", "privacy", "tags", "description", "created_at", "updated_at", "encrypted_password", "password_salt") VALUES ('b6b0302a-afae-11e3-9d31-5452009f0b27', 'test10', '48b2b00e-2363-4af1-8d32-69628d730e2d', NULL, 'derived', 'public', (NULL), NULL, '2014-03-19 21:37:47.544059+0000', '2014-03-19 21:37:47.544082+0000', NULL, NULL) RETURNING "id"     21:37:47 web.1     | (0.000279s) SELECT * FROM "user_tables" WHERE ("map_id" = '48b2b00e-2363-4af1-8d32-69628d730e2d') LIMIT 1     21:37:47 web.1     | (0.000268s) SELECT * FROM "maps" WHERE ("id" = '48b2b00e-2363-4af1-8d32-69628d730e2d') LIMIT 1     21:37:47 web.1     | (0.000359s) SELECT * FROM "users" WHERE ("users"."id" = '6fcfa361-9b41-453a-b453-1d12b76421ec') LIMIT 1     21:37:47 web.1     | (0.000286s) SELECT * FROM "maps" WHERE ("id" = '48b2b00e-2363-4af1-8d32-69628d730e2d') LIMIT 1     21:37:47 web.1     | (0.000336s) SELECT * FROM "users" WHERE ("users"."id" = '6fcfa361-9b41-453a-b453-1d12b76421ec') LIMIT 1     21:37:47 web.1     | (0.000256s) SELECT * FROM "maps" WHERE ("id" = '48b2b00e-2363-4af1-8d32-69628d730e2d') LIMIT 1     21:37:47 web.1     | (0.000334s) SELECT * FROM "users" WHERE ("users"."id" = '6fcfa361-9b41-453a-b453-1d12b76421ec') LIMIT 1     21:37:47 web.1     | (0.000255s) SELECT * FROM "maps" WHERE ("id" = '48b2b00e-2363-4af1-8d32-69628d730e2d') LIMIT 1     21:37:47 web.1     | (0.000350s) SELECT * FROM "users" WHERE ("users"."id" = '6fcfa361-9b41-453a-b453-1d12b76421ec') LIMIT 1     21:37:47 web.1     | ETHON: Libcurl initialized     21:37:47 web.1     | ETHON: performed EASY url= response_code=0 return_code=couldnt_connect total_time=0.000113     21:37:47 web.1     | Completed 400 Bad Request in 327ms (Views: 0.4ms | Models: 0.0ms)     21:37:47 web.1     |     21:37:47 web.1     |     21:37:47 web.1     | Started GET "/assets/layout/loader.gif" for 10.254.184.166 at 2014-03-19 21:37:47 +0000     21:37:47 web.1     | Served asset /layout/loader.gif - 304 Not Modified (0ms)


how to transform DB stored EPSG data to Geography data type in sql server

I have a full, 3 GB, database in SQL Server stored with certain points stored in pure float values. The points are stored in EPSG 3765. I would prefer to have them saved in the Geograpy datatype. I can not store them directly, because the points are stored like this 1734994.3283045683 5740529.6114465743 ( x, y ) Geograpy requires values between 90 and -90 as well as 180 and -180. How can I do this? I am googling for the last 4 days and cant find a solution. I would not mind writing any code (C#, python, .... sql, ...) I have no shapefile or anything else. just the data in the DB.


How can i store shape file in postgress database as a table using php?

I want to create a table from a shape file in Postgres. To do that, I can run this command in cmd:               shp2pgsql -s 4326 -W LATIN1 E:\Any\THA_adm public.test123 | psql -h localhost -d DATA -U postgres -w 123;      This works, but I want to run it from php, How can I do that, I tried php functions:   1. shell_exec()    2. exec()    3. system() but did not work. Is there any other way to create a table from shape in php or run this command?


How to tell when a new MXD is opened

I am trying to force a tab from a tab container to open to a specific tab once a new MXD is opened. Currently it will stay on a previously selected tab when I use the catalog(docked) to open a new MXD. My question is how do I tell if a new MXD is loaded? Thanks -Chris


I would like to change all the values in a attribute field to the same text

So I recieved a great answer to this question, How to change all values in a field in one operation?, although I need some further information. I would like to change the values in the entire field to the font name that I would like it to be. When I type in the font name, it claims that it is an invalid expression. Another example of what I would like to do is change all values of another column to half, related to the vertical alignment.


How do I split a vector layer into separate grid layers?

I need to split a vector layer into 1km x 1km squares so I can carry out some calculations for these individual areas. Is there a feature to do this?


How to migrate from ArcGIS for Desktop to QGIS?

Our lab has about 10 projects in ArcGIS 9 (a couple are made in ArcGIS 10) that includes the following in each:   1. A georeferenced (WGS1984 geographical projection) raster seafloor image, representing areas of about 20m x 20m.   2. A layer of points across the image, each of which is hyperlinked to documents residing on the same computer.   3. A few layers of points (feature class) representing features in the raster image.   4. A few layers of polygons (feature class) representing features in the raster image. With the above, we have been able to:   1. Analyse the areas outlined with the polygons and added buffers to some of them.   2. Answer questions like: Does this kind of fish (represented by a layer of points) always occur near a certain rock type (represented by polygons or points).   3. Add new seafloor images (as new layers) representing the same area, but taken at a later time, to see if the distribution of features in it have changed over time (i.e. there's 15% more fish, etc.). In the future, we would like to add bathymetric data so that we can make a 3D view of seafloor represented by our images. Finally the question: Can the above be successfully migrated to QGIS and **HOW** do I do it? (so we can save licensing fees and also use open source software, etc.) Based on my initial look a QGIS, most of what I outlined above can be done, but can anyone confirm that? If we do this, it needs to be done within the next couple of months, all by ourselves with no programming ability. (we don't have funds to hire someone to do the migration for us, or to add feature to QGIS).


Visualize text file with location info and intensity

I have an ascii text file containing location data (column 9-lat and 10-long) and intensity(column 20)                 200010  207  020311    40658.5  406593 52 344927.31 7100203.50  -26.2078720  127.4491855 345060.64 7100369.14   26.4  650.3  628.0 55471.293    20.168 55648.817 55637.523  -146.062      the text file has many lines 10k+ I am trying to export this ascii text into a GeoTIFF using GDAL, but not sure how to proceed. Ideas?


ESRI fc is called database.DBO.fcname and version is called SDE.fcname

ArcSDE 9.3.1 and MS SQL 2008 R2 In ArcCatalog the feature class is called database.DBO.fcname and when I look at the connection properties in ArcCatalog the default is called DBO.Default and the other versions are called SDE.versionname1 and SDE.versionname2. We just migrated to MS SQL 2008. I did not do the install of sde on the new server. The databases were not created new on the new server instead they were backed up and then restored to the new server from a MS SQL 2000 db and then upgraded to 2008 databases. Is the sde on SDE.versionname1 referring to the schema or the owner? Thanks.


What is a .mem file?

I’m running ArcServer 9.3.1 on a Windows Server 2008 virtual machine. I noticed the other day that my C:\ drive was almost full. I tracked the issue down to this location “C:\Users\ArcGISSOC\AppData\Local\Temp” on the server. There are some .mem files that are written everytime a service is accessed. For the most part they are small, but every now and then one gets written that’s 2GB. That is causing issues. I’ve looked on the ESRI site and can’t find a resource to explain these files and what is being written to them. Can anyone point me to a resource on these files that explains what they are, or have an explanation as to why I’m getting these large files? Thanks.


about spatial joining and getting area of attribute from each subgrid

I have grids of 3 x 3 sq km, i.e., 9 sq kms which I want to divide into 4 subgrids. I have already done so. I neach grid of 9 sq kms, I have digitised land cover types from Google Earth and so in the attribute table I have polygons having area of each land cover type. Now as I have split the grids into subgrids, I want the area of each land cover type from these individual subgrids. I have tried spatial joining However, there is one polygon having a contiguous expanse across all subgrids. After spatial joining, that agricultural land area is returning a common area value for all subgrids. It is not getting divided according to the subgrids. Please help me with this. I am using QGIS.


Pansharpening in QGIS

I was wondering if anyone has successfully been able to perform an image pansharpening within QGIS? I have utilized OpenEV to do this, but the results have been very pixelated. Any info would be great. Thanks...


Where do I download the ArcGIS 10.1 SP1 Web Adaptor?

I've been looking all over the place and I can't find the download! Where do I find the ArcGIS Server Web Adaptor?


SVN Version of GRASS and QGIS

I have GRASS6.5 SVN installed and wanted to go for SVN installation of Quantum GIS in Ubuntu 11.04. So as instructed in QGIS Wiki, I added the qgis repository to my `sources.list` file. Now I am able to see **qgis1.7.0svn** version in Synaptic but when I try to install it, **GRASS** is a dependency. How can I install only Quantum GIS and not GRASS from here.


Fastest way to convert shapefile to postgis in java

I have been using geotools to read shapefiles and import them into postgis. Unfortunately, when large files are being imported, this leads to a massive lag time --8000% of shp2pgsql piped into pgsql-- which I assume is caused mostly by individual inserts of data. Assuming that the amount of time it takes to write the features to the dump is not excessive, this should massively decimate the exorbitant amount of time required to import the shapefile. My plan is to just convert all of the feature attributes into rows of strings and then just use the copy functionality in the jdbc to bulk load the data --but I do not know of any functions in the jdbc or elsewhere that can convert java objects into postgres dump comptable strings. Below is a sample output from shp2pgsql --which is essentially what I want to replicate using java.               SET CLIENT_ENCODING TO UTF8;     SET STANDARD_CONFORMING_STRINGS TO ON;     BEGIN;     CREATE TABLE "raw_data"."mud" (gid serial,     "objectid" numeric(10,0),     "tcmud_name" varchar(50),     "shape_area" numeric,     "shape_len" numeric);     ALTER TABLE "raw_data"."mud" ADD PRIMARY KEY (gid);     SELECT AddGeometryColumn('raw_data','mud','geom','2277','MULTIPOLYGON',2);     COPY "raw_data"."mud" ("objectid","tcmud_name","shape_area","shape_len",geom) FROM stdin;     14  LAKEWAY MUD 9.60083055833e+004  1.24005179036e+003  0106000020E50800000100000001030000000100000005000000E03113F38C20474188DA1ADB9142634140726A6B94204741B07B4EE669426341007FCF91FE1F4741D0BE4D2668426341A03E7819F71F474190631C1B90426341E03113F38C20474188DA1ADB91426341     16  LAKEWAY MUD 1.99805337165e+004  5.90416771226e+002  0106000020E5080000010000000103000000010000000900000040CA100AE92F474110A13411B3466341E09D3C7DE42F474198B6CE82AA466341A0BA19BD912F4741F08A1B80AF466341940F2FE9912F4741E9D4FBCEB24663410005BDDD902F4741B0CC491BB646634120440ECA8C2F4741403B9018BC466341E0865D22ED2F4741D06DCE02BA466341E0AAB0DAEA2F474170988569B646634140CA100AE92F474110A13411B3466341     17  LAKEWAY MUD 2.55495765553e+004  6.75755882634e+002  0106000020E5080000010000000103000000010000000800000080183EE7062F4741A03A533557466341208801F3F52E4741E836A58E55466341805FC108D32E4741685B85F96446634140DE1C01CF2E4741D8E7CCEA6A466341A063D3DFC92E47415008A8CE71466341E0CBB53AF42E4741A848726374466341C0F32A05292F4741588B715F5A46634180183EE7062F4741A03A533557466341     20  LAKEWAY MUD 1.54327655151e+004  5.08037611343e+002  0106000020E5080000010000000103000000010000000500000000A3626F4F284741402A23586044634140FFC13E36284741A8F42ECF4D44634100E74A4A06284741604256135244634140943F531F284741409C39B66344634100A3626F4F284741402A235860446341     27  WTCMUD NO. 7    6.92786595456e+005  3.61805221358e+003  0106000020E50800000100000001030000000100000018000000C0063953CF464741107CA79ED538634100190A95D246474110AE1AA1DA386341C0240A3EDE474741D866062A0439634180FC86331F484741B8A81D80EE386341E03C4D782848474190155C7FC2386341203E3FAAED4747418056C1F598386341C0BE4FB8ED474741B0F4E3B28A386341A08A5238C947474180DBCB047E386341E036F588A6474741389250786438634140FEF57170474741205693EF523863416028CABB3F4747414897B23D503863410081AAA732474741989B78844F386341601450391A474741903F8F535138634140948698EB46474170D55FC754386341A0068BD6D9464741703CF3B85C386341001E6513CC464741E0D4F5E06238634140A9C648CC464741C8048B7274386341A06B007ACC4647415043AEA684386341E07EA77ECC464741D8E2993086386341A00B156BE0464741B845446E91386341008BC1BB15474741C004DF82AF386341C05D1BB7EC464741401E7AF3C0386341C0EF4FB5CB464741B0DF1CFCCE386341C0063953CF464741107CA79ED5386341     ...     \.     CREATE INDEX "mud_geom_gist" ON "raw_data"."mud" USING GIST ("geom");     COMMIT;      To sumarize, I am looking for the fastest/easist/robust way to bulk upload any shapefile to postgis in java. I do not care whether solutions use geotools, as I have found geotools very cumbersome for my very simple needs -- changing the name of the geometry column from "the_geom" to "geom", for example, requires generating a brand new schema from scratch with the new geometry name, iterating through all the features of the shape datastore, generating a new feature in the new schema, and then individually mapping each attribute of the old schema to the corresponding attribute in the new schema!


How do you make Amazon Cloud GIS Server accessible outside local environment?

I have services in the cloud and want to access them from applications outside the cloud. How do I set permissions to consume services outside of the cloud? Also, what is the format or syntax used for referencing Amazon Cloud GIS Server? So, when I try to add GIS Server what would type in?


Unknown Old Coordinate System of the Phillipines

I am reading a report made in 1979 and the coordinates show N 46,400 and E 47, 600. This is in the Philippines. Would anybody know what the coordinate system was used in this numbers? Thanks.


Append Point Shapefile Location by GPS'd text file

In my industry it is common to create preliminary design maps for point locations that have a unique identifier. Upon the actual construction on the design, points change location based on construction. Is there any way in ArcGIS 10.0 Editor license to append the new XY information (update the geometry) for those points? Python is preferred.


How to calculate average distance from points to points?

I have a point layer of 400,000 voter residencies and a point layer of 500 polling places. I am trying to find the average distance it takes voters to get to their polling location. Distance by road would be optimal, but I would settle for as the crow flies as well. I have ArcGIS Basic with Network Analyst. I just got NA last week, so I have no experience with it. I have a related field that can match individual voters to a polling location. I appreciate any ideas. Thanks!


Color balance and publish large amount of high resolution imagery in geoserver

I have a large amount of RapidEye images (+300 Gb) that I want to use as a basemap for a online application. So far I've managed to serve them as an ImagePyramid in geoserver following the steps shown in a presentation called "Geoserver on steroids". The problem is that I haven't been able to do a proper color balance so when looking at the entire dataset, the mosaic looks very poor but as I zoom in, it gets better (pictures bellow). My process is:   1. Convert all images to 8bits, epsg 4326, cubic interpolation using gdal   2. Generate a gdal virtual raster with all the images (gdalbuildvrt)   3. Generate the image pyramid (gdal_retile with compression and geotiff tilling) and publish on geoserver (using histogram stretch in the layer style) Any tips on how to improve the process and get better color balance? ![enter image description here](http://i.stack.imgur.com/aOVvI.png) ![enter image description here](http://i.stack.imgur.com/vBJME.png) ![enter image description here](http://i.stack.imgur.com/LHYRv.png)


Delete multiple features easily in QGIS?

I have a vector layer being rendered in QGis. I want to delete various features in the layer. I'm currently deleting each node by selecting the node through node tool and then pressing the delete button. But I want to delete very large portions in the map. Using this method will cause a lot of time. Is there any other method where I can draw a polygon around the feature which needs to get deleted and delete them easily ?


Raster image outside Grass screen, How to get it inside from start?

HI ALL: I'm new in Grass and need some help! I have created a LOCATION based on the coordinate system in which I want my data to be projected. So far I'm being able to Import & display the raster image (Tiff Image). However, to see this image have to click severals times in zoom-in/out, I understand this is because the data is somewhere out of the "screen" but still inside the LOCATION. Is there to get this image inside de screen from the start? **Location created:**   proj=UTM datum=wgs84 zone=30 towgs84=0,0,0,0,0 no_defs=defined **The Geo-TIFF image (MOD13Q1-NDVI and MOD13Q1-EVI ): was reprojected using MRT-Modis:** -UTM -Resampling: nearest Neighbor -utm zone: 30 -datum=wgs84 -UL 39.9999 / -13.05407 -LR 29.9999/ 0.0 Despite creating few locations to get this data "inside the screen", couldnt make it. Don't know where my mistake is. MANY THANKS IN ADVANCE!


how to overlay a polygon over SpatialPointsDataFrame and preserving the SPDF data?

I have a `SpatialPointsDataFrame` with some additional data. I would like to extract those points inside a polygon and at the same time, preserve `SPDF` object and its corresponding data. So far I've had little luck and resorted to matching and merging through a common ID, but this works only becaues I have gridded data with individual IDS. Here's a quick example, I'm looking for points inside the red square.               library(sp)     set.seed(357)     pts <- data.frame(x = rnorm(100), y = rnorm(100), var1 = runif(100), var2 = sample(letters, 100, replace = TRUE))     coordinates(pts) <- ~ x + y     class(pts)     plot(pts)     axis(1); axis(2)          ply <- matrix(c(-1,-1, 1,-1, 1,1, -1,1, -1,-1), ncol = 2, byrow = TRUE)     ply <- SpatialPolygons(list(Polygons(list(Polygon(ply)), ID = 1)))     ply <- SpatialPolygonsDataFrame(Sr = ply, data = data.frame(polyvar = 357))     plot(ply, add = TRUE, border = "red")      The most obvious approach would be to use `over`, but this returns the data from the polygon.               > over(pts, ply)         polyvar     1        NA     2       357     3       357     4        NA     5       357     6       357


How to Make a File Geodatabase Unlockable

I have several replica File Geodatabase's on my server. I publish maps to be viewed by other departments in ArcReader. These ArcReader pmf's use the replica Geodatabase's on the server. I need to update them from time to time, but the Geodatabase's are always locked during work hours. Usually I just update at night, but when someone needs something done right away, it's a pain to ask everyone to shut it down so I can update it. Is there anyway to prevent the users, who are only reading the data in the Geodatabase's, from creating lock files in the Geodatabase's? Can this be done by setting user permissions on the folder? Can it be done with some sort of ESRI security measures?


Stroked and filled vector lines in Openlayers

I am drawing vector lines, using GeoJSON, in my web application that uses Openlayers. I want the line to be displayed in yellow color with black border. I am using the following for the line styling:               my_layer = new OpenLayers.Layer.Vector("my_layer", {                 styleMap: new OpenLayers.StyleMap(new OpenLayers.Style({                         fillColor: "#ffff00",                         strokeWidth: 6,                         strokeColor: "#000000"                             }))             });      But 'fillcolor' does not affect line symbology. How can I achieve that effect?


Adding barriers to a Network Analyst layer programmatically

Right now, I have code that allows the user to indicate on a map, which contains a Network Dataset layer, stops that are needed to create a route. I need to add functionality to the interface that allows the user to also create barriers as well as stops. The way this is normally done in ArcMap is as follows:   1. Create the route layer by clicking on **Network Analyst** and **New Route**   2. Make the Network Analyst window visible   3. Highlight the **Restriction** layer inside the **Line Barriers** layer   4. Click on the button for **Create Network Location Tool**   5. Draw the barrier on the map. I already have implemented the ability to create stops by using `ControlsNetworkAnalystCreateLocationToolClass()`. My question is, how can I also implement the capability of creating barriers as well as stops into my code?


Zonal Statistics for Polygons with R

I'd like to do zonal statistics / terrain analysis with polygon overlays - I achieved to get min / max like so:               library(raster)     library(rgdal)          alt <- getData('alt', country = "AT")     gadm <- getData('GADM', country = "AT", level = 2)     gadm_sub <- gadm[1:3, ]          plot(alt)     plot(gadm_sub, add=T)          extract(alt, gadm_sub, fun = max, na.rm = T, small = T, df = T)      But for slope / aspect I get this error which I don't grasp how to prevent:               extract(alt, gadm_sub, fun = function(x) terrain(x, opt = "aspect"), na.rm = T, small = T, df = T)     # > Fehler in FUN(X[[1L]], ...) : unbenutzte(s) Argument(e) (na.rm = TRUE)


integrating spatialite and php (sqlite3) on windows platform

anyone can give me an explanation about how to integrate spatialite with php(sqlite3)-windows platforms. I have tried to refer from this link http://www.gaia- gis.it/spatialite-2.4.0-4/spatialite-cookbook/html/php.html But I do not understand is why the error message appears as follows.               "Warning SQLite3::loadextension(), Not supported in multithreaded Web Servers"      i use libspatialite-2.dll


How to convert Layout View to DXF or DWG?

How can I convert Layout View to a DXF or DWG format? I have tried the PDF way i.e.   * Export to PDF    * Convert from PDF to DWG using PDF2cad  It was not a great experience because the DWG does not seem good and is not georeferenced. Can anyone offer help or experiences to help solve this issue?


select by location

When performing an overlay or intersect I would like to know if there is a way to ignore selection when intersection or overlay is less than a certain percentage of area or less than a certain number of crossing map units of the perimeter of a polygon.


How can I display with OpenLayers the geometry returned from a ST_Buffer in POSTGIS?

Im using PostGIS, GeoServer and OpenLayers, and the thing is that i want to display the geometry result from a query like this:               SELECT ST_Buffer(ST_GeomFromText('POINT(100 100)'), 500);      To perform queries I'm using Java, the line I use to perform a query is the following:               resulset = statement.executeQuery(query);      where resulset is a ResulSet type object, when I convert it to String I get the following:               org.postgresql.jdbc4.Jdbc4ResultSet@163b102      My question is how do I interpret that to display the POLYGON I suppose is in there in a layer? or am I doing something wrong? Hope someone can help me with this question, thanks.


Limits when using the SQL API

i use the SQL-API as described in here: And i find the following issues undocumented:   1. How many requests per second can i send to the SQL-API, before the user gets slowed down/blocked? I would like to throttle myself and also make calculations, how long those requests will take.   2. How many statements can be put in a batch? Or is this purely a limit on the POST request? Whatever fits in the POST request, will be executed? Does anybody know or can clarify?


How to refresh a combo-box in a QGIS 2.2 form without closing and opening it again?

I'm setting up some embedded forms on QGIS 2.2 and I would need like a refresh button for forms to update choices in a combo-box (value relation field). Indeed, I currently have to close my form and open it again so that QGIS can offer me to choose my new data (in a combo-box). For example, I want to attach to someone a job in a form with several tabs (pictures 1&2). The 1st tab (1st picture) allows me to choose a person who's already in my database thanks to a combo-box (value relation field). ![1st form, tab 1](http://i.stack.imgur.com/lxI7k.jpg) In my second tab (2nd picture), I can add a new person in my database thanks to a relation reference to a table with all my people. ![1st form, tab 2](http://i.stack.imgur.com/VPMEu.png) A second form (3rd picture) appears in order to create the new person (after clicking on the updating pencil and the green plus). ![2nd form](http://i.stack.imgur.com/OZFck.png) After clicking on “OK”, back to the 1st form, I click again on the pencil to save my new person. Coming back to the combo-box on the 1st tab to assign the job of my new person, the person isn't present in the combo-box. To have it updated like in the 4th picture, I have to close the 1st form and open it again. ![1st form, updated](http://i.stack.imgur.com/YWuFG.png) That's why a refresh button for the forms would be handy to avoid useless clicks. Two sorts of actions may be possible:   * updating the whole form, but it could be very slow refreshing everything   * or just a combo-box **Thank you in advance!**


How to prevent transparent borders (collars) of a tif wipe out overlapping regions of adjecent tif while using gdalbuildvrt

I have a set of Geotifs that have got transparent borders (collars) of irregular shapes & dimensions around them. First I created a VRT over this set of Geotifs through gdalbuildvrt. Then I generated tiles through gdal2tiles.py using this VRT. I noticed that the transparent borders (collars) of some Geotifs wiped out the overlapping regions of the adjacent Geotifs creating holes in the output. Is there any setting that I should use while creating VRT through gdalbuildvrt or while generating tiles through gdal2tiles.py and the above created VRT, so that I can get rid of the holes in the output caused by the transparent borders (collars) of some Geotifs wiping out the overlapping regions of adjacent Geotifs.


postgis union geometries in consecutive order

`Database is postgres 9.2, POSTGIS="2.0.2 r10789" GEOS="3.3.6-CAPI-1.7.6" PROJ="Rel. 4.7.1, 23 September 2009" LIBXML="2.7.8"` i have geometries in consecutive order. Here is query:               SELECT (ST_Union(geom)) FROM roads WHERE osm_id IN (140122771,140122775,9620011,9620010,174011196,9620025,174011193,174011195,9620024,140119144,4377642,183432458,193166932,140119146,183432465,140119145,183811230,25585759,25585692,183811266,25585671)      osm_id is in order that i need them. after i run second query:               SELECT (         ST_Line_Interpolate_Point(             ST_LineMerge(                 ST_Line_Substring(                     (SELECT (ST_UNION(geom)) FROM roads WHERE osm_id IN (140122771,140122775,9620011,9620010,174011196,9620025,174011193,174011195,9620024,140119144,4377642,183432458,193166932,140119146,183432465,140119145,183811230,25585759,25585692,183811266,25585671)  ),                      0.1, 0.101                 )              ),0.1         )     ) as punkts       it produces point on line. ( i rand 9 queries from 0.1 to 0.9. 0.1 is first, 0.2 is second and so on), but points not in consecutive order.   ![points on line , must be from bottom to top](http://i.stack.imgur.com/4C3dy.png) You can see, that 4 point is at end of line .. Is there way, that i can union geometries in consecutive order? or maybe geometries can concatenate ?


Why does it take so long to run zonal statistics?

QGIS 2.0 Windows 7 I have pulled agricultural data from NASS CDL in tiff format. I want to track where a certain crop has grown over the years. My intention was to use a county shapefile overlaid against the CDL data, create buffers on the counties, and run zonal statistics to show the mean, count, etc. for each county containing that specific crop in each year. This should show where that crop's growth has moved to. However, when I try to run the zonal statistics, my computer bogs down and qgis claims "Not Responding." I have tried to replicate this process on much smaller pieces of data, but I am running into the same problem. Does anyone have any ideas? I realize that this is somewhat vague. I will answer any questions to help search for the answer.


import qgis.core problem with QgsFeature and QgsGeometry

I've tried several previous answers to similar questions - no luck. I run Win7. I installed QGIS 2.2.0 using the OSGeo4W 32-bit installer. Problem 1) from qgis.core import * failed. Solution 1) After everything related to setting Path and PYTHONPATH variables failed, I had to copy C:\OSGeo4W\apps\qgis\python\qgis to C:\Python27\Lib\site-packages (and yes, including C:\OSGeo4W\apps\qgis\python\qgis in both Path and PYTHONPATH did NOT work). Problem 2) The above approach got the qgis.core import started, but resulted in the following error message related to QgsFeature and QgsGeometry:               Traceback (most recent call last):       File "D:\Research\inhaca\OBIA\model_script.py", line 1, in <module>         from qgis.core import *       File "C:\Python27\lib\site-packages\qgis\__init__.py", line 35, in <module>         from qgis.core import QgsFeature, QgsGeometry     ImportError: DLL load failed: The specified module could not be found.      Now I'm lost. I tried copying all files from C:\OSGeo4W\apps\qgis\include and C:\OSGeo4W\apps\qgis\lib to the respective (include and Lib) folders in C:\Python27 - still didn't work. I'm fumbling in the dark here, any pointers?


Geocoding Levels - zip code versus block face

I'm looking at some massive data, most of which has been geocoded to the block face (front of the building) level, while some of which has been geocoded to the zip code level (less accurate), giving me a lat/long pair for every sample. When we run a test on the accuracy of the zip code level geocoding, there's a largish error compared to address level coding (test was possible since some limited sample addresses were available even for samples coded at zip code level) - a median difference of 0.63 miles. However, on plotting the lat/long pair on google maps, I see a very specific location.   1. What then does the error tell us?    2. Is the lat/long point plot (from zip geocoding) an erroneous one?    3. Could anyone explain the difference between zip and block face geocoding, in the sense of loss of accuracy for the former?   4. Finally, if I am to convert lat/long and zip code information to census tract, does the fact that the sample has been coded to zip level make a precision difference in this exercise? Many thanks!


Has QGIS 1.7 been released as a stable version yet? I see on the QGIS site the windows download on the is still referring to 1.6 and the OSGEO4W install has 1.6.0-7 as the stable release. I had seen some links here and here among others which referred to March/April for a release. I'm not unhappy with the timeline, I'm just curious that I might have missed something. Thanks Ando


geocoding for website directory

I want to create a Worldwide directory. I only have addresses in excel sheet. I am looking for a good geocoding tool. Could you advise me because i hesitate to use geotheme or geocoding for a wordpress website made from scratch. thank you very much, Vincent


Problem with plotting map in R

While plotting the `gadm` map of Russia in `R`, I have a problem at the 180 degree meridian, where part of the Russian region is plotted separately. You can see it on the plot using the following code (it's the edge left and part is on the right). Can anyone suggest how to avoid this problem? Thanks in advance. Here is a short version of the code for my question:               require(sp)          rus<-url("http://www.gadm.org/data/rda/RUS_adm1.RData")     print(load(rus))     gadm$regions = as.factor(1:88)     spplot(gadm,"regions")      ![Map](http://i.stack.imgur.com/SSgHj.png)


OpenLayers: Problems adding WMS layer to OSM

I need to create a map with 2 layers. The first layer (basemap) is the OSM layer, and the second layer is my custom geoserver WMS layer. This layer is in EPSG:4326 projection. Is this possible? I have read many posts but I didn't solve my problem. I don't find this geoserver layer in my map. Thanks in advance, Luca               map = new OpenLayers.Map('map'{projection:"EPSG:900913",                                    maxResolution:2319.89824519781,                                    units:'m' });           osm = new OpenLayers.Layer.OSM( "Simple OSM Map");      wms = new OpenLayers.Layer.WMS( "OpenLayers WMS", "95.172.21.2:8080/geoserver/wms";,                                      {layers: 'asl.caserta:distretti_4326',trasparent: true},                                      {isBaseLayer: false});      map.addLayers([osm, wms]);      map.setCenter( new OpenLayers.LonLat(14.22231, 41.26700).transform( new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject() ), 12 ); –


Is there a way to keep toll_booth nodes in osm2po's output?

I am writing a code that performs a custom routing on the sql relation generated by osm2po. Amongst several other things, I need to detect if the selected route (list of ids) passes through one or more toll booths, and then display them on the map. I have a table that contains the list of nodes that are toll booths with their corresponding osm id, and I would like to go through the list of ways in the generated route while matching the start and end nodes. Unfortunately, this doesn't work in all cases since some toll nodes are removed (joined) by osm2po. Therefore, is there a way to prevent osm2po from joining ways, at least in the case of toll nodes?


Using Easy Custom Labeling: Where's my missing label?

I’m using the Easy Custom Labeling plugin (with QGIS 1.8). I have one feature label which refused to show itself. It’s just not visible. When I look at the attributes table, the feature name appears in the same columns as the names of all the other features. I’ve tried moving nearby labels away from it--no difference. I’ve tried moving the feature (just to see what would happen)—nothing. If I turn on Labels in the corresponding feature layer, the label appears. I’m out of ideas. Where’s my label? Thnx TVZ


Error using ST_AddPoint for adding a new point to a polyline in PostGIS

I have two PostGIS tables, one where I store points, named " _nodes_ " and another with where I store multilinestrings, named " _lines_ ". I'm trying to append a new vertex to a line using one of the nodes executing a query similar to the following SQL sentence:               UPDATE lines SET      the_geom = ST_AddPoint(the_geom, (select the_geom from nodes where id_node='XXXXX')  )     WHERE id_line='YYYY'      But I get:               ERROR: lwline_deserialize: attempt to deserialize a line which is really a Invalid type      Any suggestions?


Learning about ArcMap's Tracking Analyst tool

I was examining the ArcMap toolboxes to learn more about their functions and found the Tracking Analyst tool. Does anybody have examples of how can this function is used and when best to use it?


How to zoom to specific lat/long point in QGIS

Does QGIS have the ability to Zoom to a specific lat/long point upon the user entering one. I can see it shows the coordinates in the bottom task bar but I'm not seeing a way to manually enter lat/long point and zoom to it. Possibly even harder, I would like to zoom to specific lat/long point even though I'm in a UTM projection.


Maximum number of point features in an OpenLayers vector layer

In your experience, how many point features can be added to an OpenLayers vector layer (new OpenLayers.Layer.Vector("Point Layer")) before it goes unusably slow? My use case is to display points from a database table. The user can decide which time frame to visualize. Therefore the result can be from very few to potentially 100,000s of points. I'd like to introduce a reasonable limit and warn the user if his query would return more features.


Open the IGxDialog in MS Access 2010 using VBA

I hope this is the right forum for this question. I have a program in MS Access that allows users to input varius aspects related to land treatments including the spatial information (point, line, or polygon). As part of the program the user can select a feature (shapefile) and the feature is automatically copied to the correct location, attributes attached, and loaded to a geodatabase. Currently I use the standard file dialog `Set fDialog = Application.FileDialog(msoFileDialogFilePicker)` to accomplish this but as ESRI moves more towards geodatabases this is no longer an option as you can only select shapefiles (.shp) with the standard file dialog. Selecting geodatabase (.gdb) does not produce the features within the geodatabase. I want to switch over to ArcMaps file dialog (IGxDialog). The code I found (see below) opens the file dialog but then immediately crashes and shuts down Access before I make a selection. Am I missing something, do I need to adjust a setting or work with something like ArcObjects in Python, or why is it crashing my Access program? Help I'm just an ecologist and in way over my head.               Private Sub Test_Click()         'Set variables         Dim pGxDialog As IGxDialog         Dim pFilterCol As IGxObjectFilterCollection         Dim pCurrentFilter As IGxObjectFilter         Dim pEnumGx As IEnumGxObject         Dim GetFile As String              '  Select Case aFilter         Set pCurrentFilter = New GxFilterFeatureClasses         aTitle = "Select Shapefile or Feature"              Set pGxDialog = New GxDialog         Set pFilterCol = pGxDialog         With pFilterCol             .AddFilter pCurrentFilter, True         End With         With pGxDialog             .Title = aTitle             .ButtonCaption = "Select"         End With              If Not pGxDialog.DoModalOpen(0, pEnumGx) Then             MsgBox "No selection : Exit", vbCritical         End If              GetFile = pEnumGx.Next.FullName     End Sub


Looking for advice regarding solar potential applications

I have the task of looking into whether it is possible to build a solar potential application which would inform the user whether a roof is suitable for solar panels and a brief report about the roof´s potential output. I have seen quite a few of these applications and they are pretty user friendly and easy to use. Being a GIS Analyst with experience in vector/raster data analysis and presentation using opensource desktop/webgis tools ( and ESRI products of course!) I am not a newcomer to this kind of project...but I must admit, the consequences of getting the methodology wrong would be quite embarassing. Does anyone have any "real" experience in this area and can report or supply advice regarding the methodology and processes involved. Is this something that can be realistically done by someone without any experience in the solar world...or should this be left to the universities? I look forward to hearing from any comments.


Failing to build PostGIS on Mac OS X 10.6.8 (Snow Leopard)

I'm getting the same problem expressed here Failing to build PostGIS on OS X 10.7 Lion But in Mac OS X Snow Leopard. Things I have checked include:   * PostgreSQL compiled to several architectures. I used the binary from http://www.postgresql.org/download/macosx/, and tested the dynlib objects with the tool file   * PROJ and GEOS compiled for i386 and x86_64. I used the frameworks from http://www.kyngchaos.com   * GDAL compiled from sources following KyngChaos instructions. I compiled it instead of using the provided framework because I'm a GDAL developer, and I'm going to work with it.   * Trying to compile PostGIS, because I'm a PostGIS developer too, and I don't want to use the binaries. It has been configured correctly. I used this: https://gist.github.com/2002495 But I'm still getting the error: `ld: symbol(s) not found for architecture i386` Other things I've tried:   * Compile ONLY for arch i386. Not working. PostgreSQL expects PostGIS compiled for both architectures. I get the error ld: symbol(s) not found for architecture x86_64.   * Compile ONLY for arch x86_64. Not working, for the same reason. I get the error ld: symbol(s) not found for architecture i386 I could provide configure and makefile logs, if needed, but the problem is this. It's a problem known in PostGIS lists, but I still haven't found a solution... Thanks in advance, and best regards


Getting text to center in a label expression

Is there an expression that can be entered in the label expression that will force text to center? Such as the ones used to change case/replace/etc. Thanks!