title: Lines as Great circles on a map

text: How can I get line features to appear as great circles on a map, instead of straight lines? Thanks,


title: openlayers multiple maps which is used

text: My page contains multiple OpenLayers map objects loaded in the JS var domWidgetMaps . Each of them with a specific navigation.               var domWidgetMaps = $('div[id^="' + this.id + '"]');     for (var i = domWidgetMaps.length; i > 0; i--) {         var options = {             projection: mercator, //EPSG:3785/900913             displayProjection: geographic, //EPSG:4326             theme: null,             /*                   Lines removed from OL 2.12                 maxExtent: new OpenLayers.Bounds(-20037508.34, -20037508.34, 20037508.34, 20037508.34),                 maxResolution: 156543.0339, //http://wiki.osgeo.org/wiki/WMS_Tiling_Client_Recommendation             */             units: 'm',             controls: [                 new OpenLayers.Control.Navigation(),                 new OpenLayers.Control.PanZoom(),                 new OpenLayers.Control.ArgParser(),                 new OpenLayers.Control.Attribution()                    ],             numZoomLevels: MAX_ZOOM_LEVEL,             allOverlays: true         };             widgetMaps[i] = new OpenLayers.Map(this.id + i, options);         var osm = new OpenLayers.Layer.OSM('osmBackground' + i, '', {visibility: true}, {buffer: 0});         widgetMaps[i].addLayers([osm]);         widgetMaps[i].zoomToMaxExtent();             }      If the user wants to compare several maps, features like: "zoom-in, zoom-out or pan" should be "shared". The easiest solution (I think) is to identify which map I'm working on (How?) and to launch               var bboxSelected = widgetMaps[<SELECTED>].getExtent();     widgetMaps[2].zoomToExtent(bboxSelected);     widgetMaps[5].zoomToExtent(bboxSelected);     widgetMaps[6].zoomToExtent(bboxSelected);     widgetMaps[9].zoomToExtent(bboxSelected);      I tried               $('.map').click(function (e) {     });      It works if I click on the map, but it doesn't work if I click on zoom-in, zoom-out or pan buttons. OpenLayers intercepts the events. How can I get the selected map based on zoom-in/out, pan? Thanks!


title: How to produce multiple maps with same layout for different layers' extents?

text: I have an area which is divided into several zones the borders of which I have defined as separate layers. The layers are defined by sql queries on a postgis table of line segments joined with a table that says which segments are part of which zone boundaries. What I want to do is produce a set of maps where each map is a zoom to layer extent for each zone and which has standard header, legend and scale (which will preferably vary with the zone but could be fixed across all maps). My first attempt was to get a new composer for each map -- that worked but was rather tedious. Next I tried producing one map and printing it and then deleting the map frame and going into qgis to zoom to the next layer extent and make the boundary visible then hopping back into the composer to load it in. This worked more or less. Map was fine, legend was OK because it was the same for all the maps and the zone name was easy to edit but the scale needed to be deleted and set up again for each map. Next I looked at Atlas Generation. The manual seems to assume that you know how it works and all you need is a reminder about the details (like so many software manuals ;) It seems to assume that you have something that you are iterating over but I can't see how to make it work -- there is nothing to iterate over to get the layers. Any suggestions welcome!! :) Russell


title: How to get a composition name in Qgis python console?

text: In Qgis python console: I want to populate a list with compositions names, like in Composer Manager, but I can't figure out how to get the names (titles). So far:               composerlist = self.iface.activeComposers()       And I get a list of `QgsComposerView` objects. Now I want something like that:               for item in composerlist:         print item.title()        ![enter image description here](http://i.stack.imgur.com/c2qwi.jpg) Thank you.


title: qgis how to split a field into field1 and field2 - comma is separator

text: How to split a field into field1 and field2 - comma is separator in field? Tryed to use field calculator /Jan


title: Can you set a QGIS print composer to render the map as default?

text: I have set up a number of templates for the print composer in QGIS. In the Print Composer, there is the option to have the map image displayed in 3 ways:   * as a blank rectangle with no updates until it is output    * as a manual refresh based on cache    * as a rendered refresh When my templates are opened, the map image loads as a blank rectangle by defualt, is there a way either in QGIS settings or in the template itself to force the map image to load as "Cache" or "Render" without having to manually change this?


title: What's wrong with my expression to control label colors?

text: I am struggling to alter the label colours according to field values (as strings) in the attributes table. In my mind what I want is: If the string in column X = 'Proposed' Then label the point using the string in column Y, with color A Else Label the point using the string in column Y, with color B End I have tried a bit, but struggling with the syntax               case when  "ownerenc" = 'Proposed'      then  (color_rgb(256,256,256))     else  color_rgb(56,56,0)        end      But I realise this is quite wrong! Many thanks for your suggestions.


title: How to use inbuilt tools in QGIS Modeler?

text: Is it possible to incorporate inbuilt tools like, for eg: vector grid of Research tools in QGIS model builder ?


title: Restrict the scale of the search widget in ArcGIS Flex API

text: I wonder if I could restrict the scale scope only for the searching widget (and not the extent of the map)? We want to narrow down the scale for searching in order to avoid too many results which may cause intranet traffic jam. I am using Adobe Flash Builder and ArcGIS Server 10.2.2.


title: Raster export for regressions analysis

text: I have the following issue: I am trying to do a regression analysis between the (z) values of two rasters. What I would ideally like to have is the raster output not as grid but rather a list of x1y1z1 x2y2z2, x3y3z3,... Is there a way to do this easily?


title: Coordinate System (Malaysia RSO Projection) not lining up in QGIS and ArcGIS for Desktop?

text: I am using QGIS 1.9.0 on Linux Mint; and I am new to (Q)GIS I have three layers (2xshapefiles, 1xopenlayer[Bing Maps]) that I am trying to line up located in peninsular Malaysia. **Shapefile_1** is a subdistrict border (WGS84; EPSG:4326) **Shapefile_2** is river system within the subdistrict (GDM2000 Peninsula RSO; EPSG:3375) **Openlayer** (set as the ProjectLayer CRC; WGS84 Pseudo Mercator; EPSG:3857) When importing the layers in QGIS, Shapefile_1 is close to the borders of the subdistrict shown in the Openlayer -- but clearly not matched. Shapefile_2 is incorrectly located hundreds of kilometers north east off the coast. Whe the same layers, with the (apparently) same CRC information are displayed in ARCGIS, they line up correctly. Any insights would be greatly appreciated.


title: Zonal Statistics-Weighting

text: I have a raster file on precipitation and I want to create zonal statistics as a table using a shape file wiht boundaries of countries across the world. However it appears that I need to take into account that in different latituted the corresponding area is different and thus I need to weight my statistic (mean etc) for area. Any suggestions on how I can do that? Thanks!


title: Loading geotiffs into PostGIS using raster2pgsql

text: I have an issue with loading some Ordnance Survey raster data into PostGIS. I have been using the following raster2pgsql command to load a tif that I have placed in my C:\Program Files\PostgreSQL\9.1\bin folder into PostGIS:               raster2pgsql -s 27700 -t 5000x5000 -I -C -M  *.tif  > os_strview.sql      I then execute the follwing psql command:               psql -h localhost -U postgres -d raster_test -f os_strview.sql      This works fine but when I tried to view the data in QGIS I could only see the raster_columns table in the selection box when trying to add a PostGIS layer. I tried using the plugin Load raster to PostGIS instead. This works but my raster image had no colour. I then saw this post. I used the Raster > Convert > PCT to RGB technique to convert a single tif to a geotif. Unfortunately though when running the psql command above I get this error message: C:\Program Files\PostgreSQL\9.1\bin>psql -h localhost -U postgres -d raster_test -f os_strview.sql BEGIN psql:os_strview.sql:2: NOTICE: CREATE TABLE will create implicit sequence "sx03 nw-geo_rid_seq" for serial column "sx03nw-geo.rid" psql:os_strview.sql:2: NOTICE: CREATE TABLE / PRIMARY KEY will create implicit index "sx03nw-geo_pkey" for table "sx03nw-geo" CREATE TABLE psql:os_strview.sql:3: ERROR: could not load library "C:/Program Files/PostgreS QL/9.1/lib/rtpostgis-2.0.dll": Invalid access to memory location. LINE 1: INSERT INTO "sx03nw-geo" ("rast") VALUES ('01000003000000000... ^ ROLLBACK I take it this has something to do with the size of the geotif I'm trying to load into PostGIS - which is 71 MB! Is there another method I should be trying? What is the limit on the size of a single image you can load into PostGIS?


title: How to create a grid constrained within a polygon in QGIS?

text: In QGIS I have created an irregularly shaped polygon. I can create a grid based on the maximum extents of the polygon via `Vector>Research Tools>Vector Grid`. However, what if I need the grid to only be displayed inside the polygon's boundaries? How can I create a grid within a polygon in QGIS?


title: How to create a polygon (fixed width buffer) around a polygon in QGIS?

text: I have created a polygon in QGIS (in this case it describes a country's maritime baseline). I need to create a second polygon whose boundary is a fixed distance from the original polygon. (in this case I need to create the territorial sea boundary, which lies 12 nautical miles outside the baseline). I think geoprocessing tools>buffer is correct tool for this job? But I can't figure out what units QGIS uses for the `buffer distance`. What should I enter in `buffer distance` to generate a 12nm buffer?


title: GDAL: C++ API to Convert Text file of Coordinate points to Raster Pixel line?

text: I am using GDAL to "Grow" areas form starting seed points. I have a text file listng the coordinates of my seed locations. The X,Y coordinates are in the coordinate system of a raster image. I want to burn these point locations into the raster image. I don't kow how to convert my point coordinates to the proper raster image pixel/line locations. I have the geotransformation of the image. I see GDAL classes in the API tutorial that return transformed coordinates from image row/col locations but not the reverse. I am famialiar with the utility that does this, but I don't see a C++ programmers API to do this. Does anyone know how to do this? This capability is fundamental to what I'm trying to do. Thanks in advance for any help you can offer.


title: Get Web Map Layer Field Value for If/else Statement

text: I've been trying to read through the ArcGIS JavaScript API but can't find what I want. I have a web map on my page that is based on my my map stored on ArcGIS online. I want to create an if/else statement that does different things based on the value of a field in the layer. But I'm having a hard time accessing the field of the layer I want. Here is part of the JavaScript in my HTML code.               var exLayer = actualmap.getLayer(actualmap.layerIds[6]).layerInfos[0];         console.log(exLayer);      With the code above, I am able to access the layer I want in my web map. The web service I want is number 6 on the map, and the layer I want on this web service is number 0. I'm trying to access the field in layer 0 called "TID", and I want to get the rows with "TID" > 0 and rows "TID" = 0 for my if/else statement. I tried the below, but it didn't work. It just gave me the layer name.               var test1 = Field(exLayer).name;         console.log(test1);      Any suggestions?


title: Why is QGIS not launching?

text: I've got a problem using qGIS. Always when I try to launch the programme with the desktop icon, nothing happens. I also checked the Task Manager and it doesn't show any process running regarding qGIS. * * * @yalestar: The target field shows this: "C:\Program Files\Quantum GIS Lisboa\bin\nircmd.exe" exec hide "C:\Program Files\Quantum GIS Lisboa\bin\qgis.bat" @Sven: I don't think that's the case, as I don't get any error message, but I'll give it try. * * * @Underdark: I had Microsoft C++ 2008 already installed, so I launched the installer and selected to repair it. However, qGIS still isn't launching. @Sven: I've tried your soulution, but qGIS still isn't launching. @mapoholic: Nope, I'm using Vista 32-bit. * * * @Andre Joost: I've tried to launch it manually, but then I get an error message telling me, that qgis.bat couldn't be found. I don't have any .qgs files on my PC, so I can't test it. * * * @M P: Thank you for the help. :) However, I can't find qgis.bat in the bin folder. But there's a file there called 'qgis.bat.tmpl'.


title: Historical Hurricane Tracks

text: I have found a great source for historical hurricane data, but I can't seem to find a way to download the track information. It appears that you can only view the data on the website. I was wondering if anyone else has had this issue, and if anyone knows a way to get access to the data (either as a .shp or .kml). Here's the website: http://www.csc.noaa.gov/hurricanes/# Thanks... * * * Yes I did see this, I was hoping there was some data already made. Wishful thinking. Guess I will make the files. Thanks...


title: Create point between two points with given distance or length

text: What I have:   * Coordinates of Point A and Point C and therefore vector and azimuth from A to C What I want:   * Coordinates of Point B (between A and C) with a distance from A to B of 35.2 meters.   * Suppose to work in every coordinate system Some code:               import shapely.geometry as g     line = g.LineString([(0, 0), (0, 100)])       print line.interpolate(35.2).wkt      ..works in UTM 32N but not in geographic coordinates?!


title: Consuming secured ArcGIS for Server services in Google Earth?

text: I am able to easily consume unsecured ArcGIS for Server (AGS) services in Google Earth (GE), but not secured services. I have created a KML Link in AGS Manager with no luck, and added the service manually as a link in GE with the same result. The service is added to GE but a red circle is displayed when the service refreshes. Has anyone been able to do this successfully?


title: Openlayers - Many broken tiles with vmap WMS in EPSG:900913

text: I would like to use the vmap WNS (http://vmap0.tiles.osgeo.org/wms/vmap0) in an OpenLayers map. I encounter many broken tiles (pink) when using the srs: EPSG:90013 contrary to EPSG:4326 and that at every zoom level. ![enter image description here](http://i.stack.imgur.com/GXlUM.png) I would like to know what could be the cause(s), and if it's a problem on the server side or on the client side.                 map = new OpenLayers.Map('map',                 { projection: new OpenLayers.Projection("EPSG:900913") });            openL = new OpenLayers.Layer.WMS( "OpenLayers WMS",                     "http://vmap0.tiles.osgeo.org/wms/vmap0", {layers: ['basic','stateboundary','coastline_01'],                                         srs: 'EPSG:900913',                                       isBaseLayer: true                                       });       map.size = new OpenLayers.Size(1000,800);       map.addLayer(openL);       var bounds = new OpenLayers.Bounds.fromArray([-180.0, 90.0, 180.0, -90.0]);       map.zoomToExtent(bounds.transform(source,dest));


title: What does property FeatureTypeName do in Styled Layer Descriptor (SLD) specification?

text: I am studying the Styled Layer Descriptor (SLD) specification. I got a general grasp of how I can use it to style a map. I still have one question about the FeatureTypeStyle though. When looking at the FeatureTypeStyle I see the property FeatureTypeName. What does this do? Say for example I got a layer called 'Rooms' and every room has a door (line) and figure (polygon) could I use the FeatureTypeName to say I only want to draw the doors?


title: When is an address in a new subdivision available for mail delivery?

text: I have customers that are complaining that some addresses validate on the USPS validation site (https://tools.usps.com/go/ZipLookupAction!input.action) but my geocoding vendor can't geocode the address. Google Maps shows recent construction in the area, so I'm guessing that the underlying data for the geocoding vendor is out of date. However, I'm wondering why the postal address is valid. What is the process by which the USPS validates addresses in new subdivisions?


title: using raster calculator

text: I have a combination method to combine two raster( raster x and y) . the method is : x+y / 1- min(|x|,|y|) my problem is about "," character .whats the "," means in the method and how can i use this method in raster calculator ?


title: ArcGIS Engine, How to change IMarkerSymbol size in runtime?

text: I'm using ArcGIS Engine 10 with C# Winforms to visualize some elements on a globe control. All elements are of type `IMarkerElement` which have symbol of type `IMarkerSymbol`. Now I want to add two buttons to increase/decrease element's size. I tried this two methods with no success 1.               ((element as IMarkerElement).Symbol as IMarkerSymbol).Size = newSize;      2.               ((element as IMarkerElement).Symbol as IMarker3DPlacement).Size = newSize;      I noticed that I can't change `IMarkerSymbol`'s properties values after creation (like size, offset, depth, ..etc.)? How to change this values?


title: How can I mosaic raster data which have different number of bands ( 3 and 4)?

text: I would like to mosaic 100 raster data. several raster data's band are 3, though last data have 4 bands. How can I mosaic raster data which have different number of bands?? If somebody knows how to solve this issue, let me know. OS:Windows XP GIS: ArcInfo 10.1 or QGIS1.8.0


title: How to interpret this pgRouting shortest path result?

text: I am trying to use pgrouting's shortestpath function. It returns a set of edges for the shortest path. However, the start node of the next edge is not the end node of the previous edge. I mean they are not connected. So while displaying I am joining the end node of the first edge to the start node of the next edge. However, I am seeing that the coordinates are repeating. I am not sure why, For eg plotting just the coordinates of the start and end node of each edge, I have something like this:               start end     A     B      C     D     A     E      As you can see the coordinates are repeating like A,A in the previous case. As a result there are lots of lines crossing over. Any suggestions or clarifications?


title: Display a georeferenced DEM surface in 3D matplotlib

text: I want to use a DEM file to generate a simulated terrain surface using matplotlib. But I do not know how to georeference the raster coordinates to a given CRS. Nor do I know how to express the georeferenced raster in a format suitable for use in a 3D matplotlib plot, for example as a numpy array. Here is my python code so far:               import osgeo.gdal     dataset = osgeo.gdal.Open("MergedDEM")     gt = dataset.GetGeoTransform()


title: Which is more efficient process in ModelBuilder - Merge/Dissolve or Union/Diissolve?

text: I'm creating a model which combines four very large files and dissolves them. At the moment, my workflow uses the Merge tool and then the Dissolve tool. However, I was wondering if using the Union tool and then dissolving would be more or less efficient. I don't need to maintain attributes for my purposes. I just need to create a layer which contains the geometry. The feature classes may overlap each other, and they may also contain overlaps within themselves. As the end result is the same, does anyone know which process would be more time efficient? Below are screenshots of my workflows: Merge then Dissolve: ![enter image description here](http://i.stack.imgur.com/JvI3o.jpg) Union then Dissolve: ![enter image description here](http://i.stack.imgur.com/02r0A.jpg)


title: Level of data available on FactFinder

text: I am trying to find data at the block level or even block group level for the SF1 files and ACS 2010 or 2009. Just recently, two months ago, I was able to get this data but for some reason I am now only able to get it at the census tract level??? I have followed all the ways in which I am suppose to extract this data. Did this change?


title: Comparison of metadata styles available in ArcGIS

text: I'm looking for an overview of the differences between the editable metadata styles available in ArcCatalog. Were the diferent styles developed for use by specific groups of users? Are there any advantages/disadvantages of using one style over another? What are the specific differences in the available fields? Complete answers to these questions might be too lengthy for this forum. A few sentances with links to pertinant resorces would be a great help. The ESRI help is severely lacking in this area (unless I've missed something). I'm using ArcGIS 10, which seems to have fewer options than earlier versions had, so I'm mainly interested in the four options built into ArcCAtalog, plus the FGDC editor add-in. However, more is always better. Thanks.


title: How to debug "msProcessProjection(): Projection library error. no system list, errno: 2" error in mapserver?

text: I think I have a misconfiguration in my mapserver that is preventing it from serving up imagery. In my map file, I have this section:                 PROJECTION         "init=epsg:4326"       END      If I leave that section in, I get this error when calling `GetMap` or even `GetCapabilities`:               msProcessProjection(): Projection library error. no system list, errno: 2      If I remove that section (not desirable), `GetCapabilities` works, but `GetMap` gives this error:               <ServiceExceptionReport version="1.1.1">         <ServiceException code="InvalidSRS">msWMSLoadGetMapParams(): WMS server error. Cannot set new SRS on a map that doesn't have any projection set. Please make sure your mapfile has a projection defined at the top level.</ServiceException>     </ServiceExceptionReport>      Which is expected because I just removed the section. So, my question is how do I debug the error: `msProcessProjection(): Projection library error. no system list, errno: 2`? Searching online, I see comments related to permissions, but apache (user running httpd) can access the following directories defined in my httpd.conf:               SetEnv LD_LIBRARY_PATH /some/place/lib/     SetEnv PROJ_LIB /some/place/share/proj/     SetEnv GDAL_DATA /some/place/share/gdal/     SetEnv LC_NUMBERIC C     SetEnv GEOTIFF_CSV /some/place/share/gdal/     SetEnv LANG=C      In `/some/place/share/proj`, there is a file `epsg` that is as I expect. I get no errors from apache.


title: How to reproject base layers in QGIS once open layers are removed?

text: I recently discovered the wonders of the osm maps. however now all my projects have all gone curved. is there a way to go back to how they were set originally. opposite to this question: How to reproject base layers in OpenLayers? So basically how to reproject base layers once open layers are removed.


title: ArcGIS constantly shut down by itself

text: I am trying to open Census tract data, downloaded from http://www.census.gov/geo/maps-data/data/tiger-data.html Every time, I open the file in ArcGIS, the software will return an "serious application error...unable to continue" and close by itself after a few minutes. Does anyone have the same problem, or know how to deal with it? I have tried different computers, different state's file, and redownlaoded files... Thank you!


title: Is it possible to link polyline features on a web map to a Google Chart?

text: I would like to display a line chart and (leaflet) map beside each other on a web page. When I turn off a series on the chart the associated polyline feature (river) should turn off on the map. I am having trouble finding examples of this functionality. Is it possible?


title: How to use pgRouting (pgr_drivingDistance) with GTFS?

text: I used gtfs2pgrouting to import gtfs data in postgresql, then I have these tables               List of relations      Schema |    Name     | Type  |  Owner        --------+-------------+-------+----------      gtfs   | agency      | table | postgres      gtfs   | calendar    | table | postgres      gtfs   | frequencies | table | postgres      gtfs   | routes      | table | postgres      gtfs   | shapes      | table | postgres      gtfs   | stop_times  | table | postgres      gtfs   | stops       | table | postgres      gtfs   | trips       | table | postgres     (8 rows)      (several tables as "shapes" are empty) But is it possible to use pgr_drivingDistance on this or is it another way to do it? My goal is to create isochrones Thanks


title: How to store both fixed and variable attributes as a feature class?

text: I have a live data feeds as an XML file from the Internet. I need to use Python and ArcPy lib to store each feed as a point feature in a feature class of an ArcSDE datbase. The xml looks like:               <array-list>         <spider-pos-wSVO brigade-no="36145" heading="280" altitude="302" longitude="146.09046833333332" latitude="-37.054605" personnel-no="0" speed="13">             <attrA>12</attrA>             <attrB>0</attrB>         </spider-pos-wSVO>         <spider-pos-wSVO brigade-no="36145" heading="280" altitude="302" longitude="146.09046833333332" latitude="-37.054605" personnel-no="0" speed="13">             <attrA>12</attrA>             <attrC>0</attrC>         </spider-pos-wSVO>         <spider-pos-wSVO brigade-no="36145" heading="280" altitude="302" longitude="146.09046833333332" latitude="-37.054605" personnel-no="0" speed="13">             <attrA>12</attrA>             <attrD>0</attrD>         </spider-pos-wSVO>         ....     </array-list>      Obviously each element has a fixed set of attributes (heading, altitude, longitude and latitude, etc etc) and a various set of child-elements (attrA, attrB, attrC and attrD, etc). I am able to store the fixed attributes as attribute fields in the SDE feature class but have no idea how and where to store the variable ones. Has anyone got a best practise for it? Thanks


title: Are US Census tracts GEOID the same across different years?

text: I want to know if it is appropriate to join tables from different years to the same census tracts TIGER file. Do all census tracts' ID stay the same? How would they change--adding new tracts, modification, merging, etc. ? Thank you! S.


title: Manipulating Azimuthal Equidistant Projections in QGIS

text: I have a shapefile of world countries that is projected as Azimuthal Equidistant with Chicago, USA as the center. I would like to change the projection such that Edinburgh, UK is the center. One would think that this should be as simple as changing 2 numbers in the coordinate system parameters... The problem that I seem to be experiencing is that QGIS makes it difficult if not impossible to view the complete proj4 code of a given shapefile's coordinate system parameters. This is what I have tried so far:   * Added shapefile of world countries projected as Azimuthal Equidistant with Chicago as the center   * Opened properties, metadata tab, copy proj4 parameter code   * Went to settings, custom projection, paste in proj4 code, changed the lat/lon from that of Chicago to that of Edinburgh   * Saved shapefile of world countries, specified my new custom projection, added new file to map   * My new file looks like a map of Pangaea, with the continents all smushed together. I think there must be other parameters aside from those shown in the proj4 code of the metadata tab. The only thing I'm changing in the proj4 code is the center of the map, but the change I'm seeing is much more than that. My observation suggests that more than one independent variable is changing. I would be happy to send my shapefile to anyone who would like to play with it. Thanks, Daniel Wolf Environmental/Geospatial Enthusiast


title: create the postGIS table error No such driver: PostgreSQL

text: I'm trying to create the postGIS table with spatial features I can do it by:               library(RPostgreSQL)     library(spacetime)     library(rgdal)     dbname = "geodatabase"     user = "postgres"     password = "***"     #password = ""          drv <- dbDriver("PostgreSQL")     con <- dbConnect(drv, dbname=dbname, user=user, password=password,host='localhost', port='5432')          data(air)     rural = as(rural, "STSDF")     p = rural@sp     sp = SpatialPointsDataFrame(p, data.frame(geom_id=1:length(p)))          OGRstring = paste("PG:dbname=", dbname, " user=", user," password=", password, " host=localhost"," port=5432",sep = "")          print(OGRstring)     writeOGR(sp, OGRstring, "rural_space", driver = "PostgreSQL")      I look for is how the error:                > source('~/R/postgres/myScript.R')     [1] "PG:dbname=geodatabase user=postgres password=zina host=localhost port=5432"     Error in writeOGR(sp, OGRstring, "rural_space", driver = "PostgreSQL") :        No such driver: PostgreSQL      please help me EDIT:                > ogrDrivers()                  name write     1      AeronavFAA FALSE     2          ARCGEN FALSE     3          AVCBin FALSE     4          AVCE00 FALSE     5             BNA  TRUE     6             CSV  TRUE     7             DGN  TRUE     8             DXF  TRUE     9          EDIGEO FALSE     10 ESRI Shapefile  TRUE     11     Geoconcept  TRUE     12        GeoJSON  TRUE     13       Geomedia FALSE     14         GeoRSS  TRUE     15            GML  TRUE     16            GMT  TRUE     17       GPSBabel  TRUE     18  GPSTrackMaker  TRUE     19            GPX  TRUE     20            HTF FALSE     21         Idrisi FALSE     22            KML  TRUE     23   MapInfo File  TRUE     24         Memory  TRUE     25   MSSQLSpatial  TRUE     26           ODBC  TRUE     27            ODS  TRUE     28        OpenAir FALSE     29    OpenFileGDB FALSE     30         PCIDSK  TRUE     31            PDF  TRUE     32            PDS FALSE     33         PGDump  TRUE     34           PGeo FALSE     35            REC FALSE     36            S57  TRUE     37           SDTS FALSE     38       SEGUKOOA FALSE     39           SEGY FALSE     40            SUA FALSE     41            SVG FALSE     42            SXF FALSE     43          TIGER  TRUE     44        UK .NTF FALSE     45            VRT FALSE     46           Walk FALSE     47           WAsP  TRUE     48           XLSX  TRUE     49         XPlane FALSE     > sessionInfo()     R version 3.1.0 (2014-04-10)     Platform: x86_64-w64-mingw32/x64 (64-bit)          locale:     [1] LC_COLLATE=French_France.1252  LC_CTYPE=French_France.1252        [3] LC_MONETARY=French_France.1252 LC_NUMERIC=C                       [5] LC_TIME=French_France.1252              attached base packages:     [1] stats     graphics  grDevices utils     datasets  methods   base               other attached packages:     [1] spacetime_1.1-0 sp_1.0-15       RPostgreSQL_0.4 DBI_0.2-7                loaded via a namespace (and not attached):     [1] grid_3.1.0       intervals_0.14.0 lattice_0.20-29  rgeos_0.3-6      tools_3.1.0          [6] xts_0.9-7        zoo_1.7-11       Edit: code R:               writeOGR(sp, "C:/file.sql", "rural_space", driver = "PGDump")      code psql via cmd.exe               Microsoft Windows [version 6.1.7600]     Copyright (c) 2009 Microsoft Corporation. Tous droits réservés.          C:\Users\zina>cd C:\Program Files\PostgreSQL\9.3\bin          C:\Program Files\PostgreSQL\9.3\bin>psql -h localhost -d geodatabase -U postgres      -f "C:\file.sql"     SET     psql:C:/file.sql:2: ERREUR:  la table « rural_space » n'existe pas     DELETE 0     BEGIN     CREATE TABLE                           addgeometrycolumn     --------------------------------------------------------------      public.rural_space.wkb_geometry SRID:4326 TYPE:POINT DIMS:2     (1 ligne)               CREATE INDEX     ALTER TABLE     psql:C:/file.sql:9: ERREUR:  la colonne « OGC_FID » de la relation « rural_space      » n'existe pas     LIGNE 1 : ...ERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" ...                                                                    ^     psql:C:/file.sql:10: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:11: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:12: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:13: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:14: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:15: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:16: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:17: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:18: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:19: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:20: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:21: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:22: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:23: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:24: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:25: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:26: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:27: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:28: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:29: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:30: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:31: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:32: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:33: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:34: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:35: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:36: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:37: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:38: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:39: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:40: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:41: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:42: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:43: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:44: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:45: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:46: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:47: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:48: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:49: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:50: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:51: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:52: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:53: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:54: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:55: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:56: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:57: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:58: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:59: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:60: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:61: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:62: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:63: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:64: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:65: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:66: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:67: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:68: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:69: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:70: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:71: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:72: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:73: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:74: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:75: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:76: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:77: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     psql:C:/file.sql:78: ERREUR:  la transaction est annulée, les commandes sont ign     orées jusqu'à la fin du bloc     de la transaction     ROLLBACK      file.sql :               SET standard_conforming_strings = OFF;     DROP TABLE "public"."rural_space" CASCADE;     DELETE FROM geometry_columns WHERE f_table_name = 'rural_space' AND f_table_schema = 'public';     BEGIN;     CREATE TABLE "public"."rural_space" ( OGC_FID SERIAL, CONSTRAINT "rural_space_pk" PRIMARY KEY (OGC_FID) );     SELECT AddGeometryColumn('public','rural_space','wkb_geometry',4326,'POINT',2);     CREATE INDEX "rural_space_wkb_geometry_geom_idx" ON "public"."rural_space" USING GIST ("wkb_geometry");     ALTER TABLE "public"."rural_space" ADD COLUMN "geom_id" INTEGER;     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000BABBCE86FC2B2340B1524145D5D54A40', 0 , 1);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000C98E8D40BC5E2340D09B8A5418C34A40', 1 , 2);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000E012807F4A952340ADDBA0F65B094B40', 2 , 3);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000772CB649454B2B405A643BDF4F394A40', 3 , 4);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000009A249694BB972A402235ED629A534A40', 4 , 5);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000F435CB65A3732A401684F23E8E3C4A40', 5 , 6);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000F530B43A39DB2140E67805A227B74940', 6 , 7);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000828B1535986621407978CF81E56C4A40', 7 , 8);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000897956D28A7F1B407593180456EE4940', 8 , 9);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000B05417F0327B244038F4160FEF0B4B40', 9 , 10);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000006762BA10AB6F2A40B9347EE195684840', 10 , 11);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000037894160E5382940CD052E8F35374940', 11 , 12);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000077D66EBBD0042A4063B83A00E2A64940', 12 , 13);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000DC9F8B868C371C4050E3DEFC86E14840', 13 , 14);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000006971C63027C02440363B527DE7474940', 14 , 15);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000060CD0182399A2140670FB40243EA4A40', 15 , 16);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000BF29AC5450E92C409F1D705D31154A40', 16 , 17);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000DF4F8D976EA22140520ABABDA4D34840', 17 , 18);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000006F2F698CD6B92640AAF06778B3A64A40', 18 , 19);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000042942F6821212840423EE8D9ACE84A40', 19 , 20);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000EF54C03DCF072C408AABCABE2B484A40', 20 , 21);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000BF4692205C811F4007EBFF1CE6174B40', 21 , 22);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000F75965A6B4E62140DD4598A25CEA4940', 22 , 23);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000006A12BC218D522740C37E4FAC53654940', 23 , 24);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000078EFA83121C61C407E18213CDADE4840', 24 , 25);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000009A7D1EA33C0B214064744012F6DB4840', 25 , 26);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000042226DE34F2C2B40B97020240B6A4940', 26 , 27);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000D3FA5B02F09F24408DD2A57F49B24940', 27 , 28);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000E085ADD9CA832C40E5805D4D9EC24A40', 28 , 29);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000005131CEDF84922D40787B1002F2694940', 29 , 30);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000003B1BF2CF0CDA2940FEB7921D1BC34940', 30 , 31);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000009A42E73576712740C9586DFE5F294940', 31 , 32);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000001EFB592C45622040CC46E7FC14774940', 32 , 33);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000081EDFDE35102A40A165DD3F16924A40', 33 , 34);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000007C7DAD4B8D002240F302ECA353354940', 34 , 35);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000046EC134031021C40D2AC6C1FF2D04840', 35 , 36);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000C939B187F6812140B9E00CFE7EBB4840', 36 , 37);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000009FEB7921DDB1C40166A4DF38EDB4A40', 37 , 38);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000363B527DE7772940A7E7DD5850584A40', 38 , 39);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000030D461855B4E1F40D0285DFA97A24840', 39 , 40);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000E2E6543200BC254012691B7FA2AA4940', 40 , 41);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000043A9336553F2A40156F641EF93D4A40', 41 , 42);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000E1455F419A592B409241EE224C8F4940', 42 , 43);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000020B24813EF382640581AF8510DF34740', 43 , 44);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000008C67D0D01F2840753FA7203F274940', 44 , 45);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000008257CB9D99082140A3772AE09EE94840', 45 , 46);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000D55C6E30D4A11F40BA675DA3E5F44740', 46 , 47);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000B262B83A008A254078B6476FB8534940', 47 , 48);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000D4449F8F324A2B402B685A62657C4A40', 48 , 49);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000005C3E92921EEE2A4024F25D4A5D544940', 49 , 50);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000058CA32C4B10E1F40386A85E97BE74740', 50 , 51);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000029AF95D05D6A22408BFCFA21362C4840', 51 , 52);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000005C1C959BA8CD2540F1D58EE21CEF4A40', 52 , 53);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000008C67D0D01F1940F148BC3C9D534940', 53 , 54);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000049A0C1A6CE4B204053910A630BC94940', 54 , 55);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000F8E12021CA2723406891ED7C3FE14940', 55 , 56);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000006556EF703B242C40AF7B2B1213DE4A40', 56 , 57);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000A29BFD817283254018CFA0A17F664A40', 57 , 58);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000FFEA71DF6A1D2C4068AED3484B454A40', 58 , 59);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000006823D74D29DF234003EE79FEB43F4940', 59 , 60);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000632B685A622523405B5EB9DE36714840', 60 , 61);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000F697DD938755264018CFA0A17F7A4A40', 61 , 62);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000006D8D08C6C1851940D7A6B1BD16224940', 62 , 63);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000002CBB6070CD9D20402B685A6265764B40', 63 , 64);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000008D0DDDEC0FE41F40A01A2FDD24624940', 64 , 65);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E610000045A0FA07912C1D40F1D8CF6229B64840', 65 , 66);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000007DE882FA9639254053910A630BE14940', 66 , 67);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E6100000B37A87DBA1712940D9EBDD1FEF374B40', 67 , 68);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000008D28ED0DBE802B40CBD6FA22A15D4940', 68 , 69);     INSERT INTO "public"."rural_space" ("wkb_geometry" , "OGC_FID" , "geom_id") VALUES ('0101000020E61000001DCC26C0B0E422401FF30181CE9E4840', 69 , 70);     COMMIT;


title: Does osm2po find routes of alternating modes e.g. car, ferry, car

text: I loaded europe-latest.pbf into osm2po using the default osm2po.config and java -Xmx8g -jar ../osm2po-4.6.9/osm2po-core-4.6.9-signed.jar prefix=eu tileSize=x,c ../europe-latest.osm.pbf and looked at car routes that must use ferries, e.g. Dover to Calais or between Athens and the greek Islands etc. osm2po doesn't find any route. Does any else observe this, too? How should I proceed to make osm2po also find these kinds of routes? Thanks.


title: Can't view PostGIS rasters after changing schemas

text: I have loaded a few rasters into the Public schema of a PostGIS 2.0 database. Using the db manager in QGIS 1.8.0 I can add the rasters to the map canvas. After moving the rasters to another schema in the same PostGIS database, I am unable to view the rasters in QGIS. However, moving the raster back to the public schema I am again able to view the rasters in QGIS.   1. Why is this the case? and   2. is it possible to fix? I have tried to load the raster directly into another schema using the someschema.sometable syntax in raster2psql but this returns an error stating that my target schema does not exist, however it does exist in the database. Therefore, I have left out the someschema.sometable syntax resulting in the rasters being stored in the public schema and manually moving them to another schema. **EDIT:** Some additional information. When loading a raster to PostGIS a rid and -i column are created. Only after changing the -i column to rast does QGIS recognise the raster layer (layer is still in public schema). However, the QGIS error "not a valid layer" returns after moving the raster table to another schema. **EDIT2:** I have figured out why the raster column heading was named -i. I used the -f parameter without specifying the colummn name. The -I constrained directly followed the -f constrianed therefore resulting in the raster column name being -i.


title: Populate a layout text element with the description from a coded value domain

text: ArcGIS 10.1+ lets you easily put dynamically populated text boxes in your print layouts. Doing so puts some text in the text element that looks something like this:               <dyn type="page" property="attribute" field="sampleTable.sampleAttribute"domainlookup="true"/>      It works fine for me when I have it just pulls from the index feature class. My problem is that I have data in other feature classes in the same database that need to be joined to the index layer so I can pull relevant attributes out of there and display them in other text elements. As soon as I join one of those feature classes to the index feature class the domain lookup no longer works regardless of what it is set to. That is to say I only get the coded values in my text elements. ESRI has logged this bug (#NIM104451), but I need a workaround sooner than when (if) they get around to fixing it. My first thought was to use some python to take care of this, but it's getting more unwieldy the more I think about it. My general thought process is this:   * Make sure the text elements I want to use share the name of the domain they will pull from    * List layout elements   List domains in the database   * For every layout element, step through the domain list until the domain equals the element name    * Take the current element text (the coded value) and replace it with the equivalent description from the domain I got this far in my code before thinking there has to be a more efficient way:               import arcpy     mxd = samplemxd          domains = arcpy.da.ListDomains("Database Connections\\database.sde")          Elements = arcpy.listlayoutelements(mxd)          For elem in Elements:         for dmn in domains:             if elem.name == dmn.name:                 coded_values = domain.codedValues                 for val, desc in coded_values.iteritems():                     if val == !layer.attribute!                         elem.text = desc      Any ideas? I'm open to non-python ideas as well, in case there's something database side I can do.


title: Importing Feature Class to Oracle SDE Geodatabase Simplifies Features

text: I recently imported a feature class from a File Geodatabase to an Oracle SDE Geodatabase. The features of the two feature classes now look drastically different despite not having any edits made. ![enter image description here](http://i.stack.imgur.com/6JsH8.jpg) The green polygon is the original feature as it looked in the file geodatabase. The purple transparent polygon on top is the feature after being imported into the Oracle SDE Geodatabase. It's as if the act of importing the feature class simplified the features. Any ideas on why that would happen? I am new to Oracle and ArcSDE, maybe I set something up incorrectly? Here is another view of the same area, with different symbology: ![enter image description here](http://i.stack.imgur.com/JO900.jpg) The storage format is ST_GEOMETRY, and there was no conversion performed, as far as I know, I simply imported the feature class into the Oracle SDE geodatabase from the file geodatabase.


title: FME Features to read parameter

text: I would like to ask if there is a way in FME to set the Feature types to read parameter, from values in microsoft access table. I have a model that translates .dgn files and then appends them to a SQL table to use in ArcMap, but the model at this moment is hard-coded to specific level numbers, and it needs to be able to change these levels according to access tables that are provided by the engineers.


title: Scaling features based on an attribute

text: I would like to scale polygon features by a specific factor I have calculated for each feature in a feature class. I know of the Scale tool, however it only allows for manual entry of a scale factor and I have hundreds of features, so I would like to automate this as much as possible. The end product will have all of my features resized so that they have the same area, but still have their own shapes and dimensions. Does anyone know how I can accomplish this, preferably via python? I can also utilize qGIS if a better function exists.


title: Russian GIS Data

text: Which companies are the best suppliers of vector, geocoding, imagery, and POI data in Russia? Which companies provide a complete service (which includes the above data types and also routing, rendering, etc.)?


title: Most efficient option to create a raster basemap for a web application?

text: I'm sure this is probably a common problem and many have already found the right answer. I have a large mosaic raster file (a city's ortophoto) in .ecw fromat and we want to use it as a basemap in a web application. Due to its size (2GB) rendering is very slow, therefore increasing performance is vital. 1) My first option was importing the raster in our ArcSDE geodatabase (on postgreSQL)-> build pyramids-> then publishing it as a map service with a 5 scale cache on our ArcGIS Server. Failed for unknown reasons (so far). 2) The second option would be importing the raster in the ArcSDE GDB then->create pyramids-> publish it as an image service to ArcGIS Server. 3) Third option would be to keep the raster in the folder-->copy the raster to the server --> publish it as an image service. 4) Fourth option would be to save the raster as a layer file-> then publish it in a map service. Which of these four options would be the best when talking about rendering speed? Thank you in advance I use ArcGIS 10.1 + POstgreSQL DB + ArcGIS Server Enteprise 10.1 Update: as brenth was saying you cannot publish an .ecw. I do have all the tiff files that created the .ecw mosaic, so I will try to create a mosaic from them and then publish using the fourth method. Does it sound ok?


title: importing data (weighted average) from a shapefile from another shapefile

text: I am new to GIS and hope this is not too basic. I am using QGIS. Say I have two shapefiles: SHP1, with attributes allocated to municipalities, and another shapefile I built from a wkt file, SHP2, which is an orthogonal 60 by 60 grid. I wish to do the following: add an attribute column to SHP which reflects the weighted average of an attribute from the underlying SHP1. thanks for any help Alex


title: How to calculate the azimuth of a line of a multilinestring in PostGIS?

text: I have two spatially enabled tables in PostgreSQL with a geometry in each: First table is with POINT type, second is with MULTILINESTRING. I need to calculate the azimuth of the LINESTRING part of a MULTILINESTRING, which is closest to the POINT. So the question is: how do I extract the LINESTRING in question, and compute the azimuth?


title: Does ogr2ogr Convert from PostGIS to File Geodatabase (FGDB_BULK_LOAD YES) work?

text: I am trying to export a lot of data from a PostGIS database to a File Geodatabase. The command I used to use was               ogr2ogr -progress -append -f "FileGDB" -gt 66536 -sql "SELECT * FROM line" C:\Temp\data.gdb  PG:"host=localhost user=postgres password=xxx dbname=xxx" -nln Line -nlt LINESTRING -a_srs "EPSG:4326"       This took 34 minutes. In the 1.9.2 release of GDAL which I am using OSGeo4W there is a new switch to bulk load data into a File Geodatabase. So I changed my code to be               ogr2ogr --config FGDB_BULK_LOAD YES -progress -append -f "FileGDB" -gt 66536 -sql "SELECT * FROM line" C:\Temp\data.gdb  PG:"host=localhost user=postgres password=xxx dbname=xxx" -nln Line -nlt LINESTRING -a_srs "EPSG:4326"      And this has no performance improvement at all. I have run it several times with and without the other switches and still nothing seems to improve. Has anyone got the Bulk Load working?


title: Spatial ETL Tools won't run "AttributeError: 'module' object has no attribute, etc."

text: ![AttributeError: 'module' object has no attribute](http://i.stack.imgur.com/bSbb1.png) I have a python script that runs 5 Spatial ETL tools that I wish to run as a scheduled task. Each of the tools will run successfully from ArcCatalog and also from the ArcCatalog Command line. If I try to run the tools from the DOS window, as would be done via python script for a scheduled task, arcpy cannot find the tools, although it does load the toolbox. In the attached screenshot:   1. I open a cmd.exe window   2. start python   3. import arcpy   4. Import my toolbox, which is verified as loaded on the next line   5. Use arcpy.ListTools to see what tools are loaded. the result is "[]", meaning no tools are recognized from this tool box   6. I use the same ListTools to show the Spatial Analyst tools that were loaded by import arcpy   7. I attempt to execute one of the tools from my toolbox and forget a " (my mistake!)   8. I attempt to execute one of the tools from my toolbox and receive the message "AttributeError: 'module' object has no attribute 'CopyFeatureClass_newgpetl'". What is preventing my tools from being recognized from the DOS CMD.EXE window while ArcCatalog and ArcCatalog Command Line can run them?


title: Python literature relevant to QGIS

text: Python for ArcGIS is a book that I have been looking at. Does anyone know if this would be applicable in QGIS? I am not looking to write application/plugins, simply get better at understanding expression etc. for standard workflow. e.g. queries and data-defined settings.


title: How to merge dbf files(zonal statistics as a table) in GIS environment.?

text: ![enter image description here](http://i.stack.imgur.com/4WKm6.jpg) I am not very savvy at using Python. Actually I have around 1000 output files from using the tool zonal statistics as table in ArcGIS. The names of the files are jan_60, feb_60......dec_13. These all lie in one folder. Each of these files contain an ID and Name field which are my identifiers. And each file has got a column of interest called Mean. So what I want is that all the files be merged, using ID and name field. In the final table, I want one column added having the mean values from each of the individual tables. Note that I want the name of the new columns as the original file name. I am trying to give you a visualization below: **original files I have:** ID Name Mean **output file needed:** ID Name jan_60 feb_60 march_60 Can someone help me out, if any of you have had a similar task. It would be great if I can have a script or something. Thanks a lot. output file needed ![output file needed](http://i.stack.imgur.com/m4G3c.jpg) folder containing individual files ![folder containing individual files](http://i.stack.imgur.com/OW7oN.jpg) one of the individual ![one of the individual](http://i.stack.imgur.com/BARWm.jpg)


title: How to speed up spatialite-query using spatialindex?

text: I have got 2 tables, one of type LineLayer called strlayer, the other one is a PolygonLayer called blayer. The LineLayer contains two rows called fromBland and toBland that should be filled with an sql-query. They should be filled for each line according to in which polygon (from the Polygon-layer) the start- and endpoint of the line is. The names are in the Table PolygonLayer in a row called 'name'. With the following query it would work basically:               update strlayer set fromBland =(select name from blayer as n, strlayer     as s where strlayer.strid = s.strid and n.blevel = 4 and     within(startpoint(s.geometry),n.geometry));      blevel=4 only selects certain polygons, as there are also polygons of another level in this table, that aren´t needed in this query. Now I'd like to speed up the query as the database is very large and the operation takes several hours. As I already read it should be much faster with a spatialindex. So I´d like to know how the query above would look like with this index. I already know how to create the index - but not how to use it. Has anybody an idea how to write that query with the spatialindex and so speed it up?


title: Esri ArcGIS server: Errors rendering annotation features

text: We have a ArcGIS Server 10 site based in Map Services that has some annotation feature classes. Most of the time the annotation text displays correctly. However, sometimes the rendering of annotation features 'goes crazy'. ![Bad Annotations][1] Our annotation features are plain text, usually numerical or alphanumerical codes with 1 to 4 characters, using font Arial. When the rendering 'goes crazy' the annotations are rendering instead with: • A different font (ex. Wingdings), or • A different text angle, or • A different size, or • with underline or strikeout, or • with the characters overlapping each other's These errors occur randomly and can affect all the features displayed from a feature class or not. The way each feature is displayed just keep changing as we zoom and pan. In attachment I'm sending two images that illustrates the phenomena. Have you ever experienced the same problem? Do you have an idea what could be the source of the problem? The client side of our site is developed with the Flex API 2.4. The Map Services are not cached as the data used is subject to updates every night. Our data is an ArcSDE 9.3.1 geodatabase, based on Oracle 10. Your feedback will be welcome. Regards, Paula Lopes Altran.com Belgium


title: How to make raster cells with a certain value stand out in QGIS?

text: I would like to make my raster cells glow in QGIS. If I have a dataset where only a few pixels has a value, they can be pretty difficult to spot (unless I know where to look). Stretching helps a bit, but one pixel still occupies one pixel on the screen (or none at all if downsampling removed it), and I would like it to stand out. Is there a method or a plugin that can do this in QGIS?


title: Postgis / QGIS Snap Point to Line

text: I'm trying to snap points to lines in postgis for use in QGIS, but cannot get it to work. I have one line in my line table (although it was created by merging multiple lines together) and about 600 points in my point table so would rather not do this manually. Any help would be much appreciated! Many thanks!


title: Convert planar distance to ellipsoid (real) distance

text: I have the following linestring which approximately marks an airport runway:               LINESTRING(-381578.819887 6574044.360475, -378513.534443 6574799.475787)  # EPSG 3857     LINESTRING(-3.427781 50.732117, -3.400245 50.73641)                       # EPSG 4326      If I use the GLength function, I get a length of               SELECT GLength(GeomFromText('LINESTRING(-381578.819887 6574044.360475, -378513.534443 6574799.475787)'))     3156.925m   # EPSG 3857     SELECT GLength(GeomFromText('LINESTRING(-3.427781 50.732117, -3.400245 50.73641)'))     0.027869°   # EPSG 4326 ~ 0.027869 x 111319.9 = 3102m      The actual length of this runway is 2076m. QGIS will display this _correct_ distance if I check the ellipsoidal box using the measure tool, so I assume this performs a projection/trigonometry conversion. Yet, I can't seem to find the correct calculation to replicate this conversion and go from the (planar?) distance returned by the GLength function to the real distance. Can anyone direct me to a tutorial for this? This question is a continuation of my pervious question How can I get the ellipsoidal values for length and area using spatialite? which got side tracked into correcting a different issue.


title: Intersection algorithm that correctly handles the 180° meridian and poles?

text: Is there a well known intersection algorithm that correctly handles 180 meridian and poles? For example, let's say we have a list of latitude and longitude values that represent Antarctica. Let's say we also have a simple polygon that represents an airplane. We want to know if the airplane is over Antarctica. Generic 2D polygon intersection does not work in this case if you simply use latitude for y and longitude for x, because the flat coordinate system will have edges at -180 and 180 longitude, and -90 and 90 latitude. Antarctica will go off the page on three of the edges. How is this problem typically solved?


title: Statistics on matching point layers ignoring nodata values in QGIS

text: I have a large area and a street crossing this area. Both are vector point layers in QGIS. The large area contains -99 as NoData value. I already created the intersection of both as a new layer. I want to calculate the basic statistics from this third layer. But I want to ignore the -99 points, because otherwise my sum is not correct. How do I define the NoData values for point layers (shapefile) or ignore them during the statistical calculation using QGIS? _Question solved by OP, answer below_ Unfortunately I cannot answer my question officially, so here is my answer: I edited my intersection layer with the attribute table. I simply chose all -99 values, removed them and saved the layer. Afterwards I calculated my statistics for the modified intersection layer.


title: How to print .dwg to .pdf with teigha true view

text: I don't have Autocad software and I can only find this free tool which can't read .dwg files. Now, I need to convert .dwg to .pdf. But the `teigha true view` can't export .dwg to .pdf. But I find it offers a plot command. Unfortunately, I aslo don't have the DWF6 ePlot.pc3 printer. So, does anyone know how to handle that or is that an impossible task? Do you know any other software can convert .dwg to .pdf?


title: Is there a free way to convert DGN v8i to... any other format?

text: (Update: We have postgis and ogr2ogr (etc.) but not any ESRI or Microstation software.) I just got DGN v8i files and I can't find any **free** way to convert them to... anything else. FME and AutoDWG are not free. OpenDGNViewer (aka DGNDirect) appears to be able to view but not export. I really would be happy with any other format:   * DGN v7   * SHP   * DWG   * ? Crazy how difficult it seems to be to work with what is apparently a standard in many municipalities.


title: What are .cld and .own files? (supposedly DGN v8i)

text: What are these files?               Twp_0608.own     Twp_0609.own     Twp_0708.own     Twp_0709.own     Twp_0710.own     Twp_0809.own     Twp_0810.own     twp_0608.cld     twp_0609.cld     twp_0708.cld     twp_0709.cld     twp_0710.cld     twp_0809.cld     twp_0810.cld      Supposedly they're DGN v8i files... but I haven't found a reference mentioning these extensions.


title: Open mail client when clicking feature in Attribute Table

text: One of the attributes/features of my layer contains an email address. I'd like to be able to open the default mail client by clicking on the mail address from the Attribute Table. After some initial investigation I've found, what seems to be, two different possibilities. 1) I could make an action 'Layer Properties/Actions', or 2) perhaps I could change the fields 'Edit widget' in Layer 'Properties/Fields' to a 'Web view', and that would take care of it(?) On option 1) I'm a bit unsure how to create an action programmatically using PyQGIS. I've found that a `QgsVectorLayer` has a function `actions()` that returns a `QgsAttributeAction` object. The `QgsAttributeAction` has an `addAction` function, so I should be able to use it like this...               layer = iface.activeLayer()     aa = layer.actions()     aa.addAction(...)      but I'm not a 100% on the arguments. From the documentation:               addAction(QgsAction::ActionType type, QString name, QString action, bool capture=false)      Anyone done this before that could give a few hints?


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

text: Dark tiles are generated when I run gdal2tiles.py on a .IMG file(a scanned and georeferenced topographical map).   This question advised to apply a color relief map using gdaldem before generating tiles.   Does gdaldem work if the .IMG has no elevation data ? I tried it using the color ramp below and the output file was all green.               0 110 220 110     900 240 250 160     1300 230 220 170     1900 220 220 220     2500 250 250 250      This is because the height for all raster cells are the same, 0. Which color ramp should i use with gdaldem to prevent gdal2tiles.py from generating dark pngs ?   gdalinfo output is as below                 Driver: HFA/Erdas Imagine Images (.img)     Files: ALOI.IMG            ALOI.IMG.aux.xml     Size is 3830, 2552     Coordinate System is:     GEOGCS["GCS_WGS_1984",         DATUM["WGS_1984",             SPHEROID["WGS_1984",6378137.0,298.257223563]],         PRIMEM["Greenwich",0.0],         UNIT["Degree",0.017453292519943295]]     Origin = (32.999437318839924,3.000470692247897)     Pixel Size = (0.000392066994127,-0.000392066994127)     Metadata:       DataType=Generic     Corner Coordinates:     Upper Left  (  32.9994373,   3.0004707) ( 32d59'57.97"E,  3d 0'1.69"N)     Lower Left  (  32.9994373,   1.9999157) ( 32d59'57.97"E,  1d59'59.70"N)     Upper Right (  34.5010539,   3.0004707) ( 34d30'3.79"E,  3d 0'1.69"N)     Lower Right (  34.5010539,   1.9999157) ( 34d30'3.79"E,  1d59'59.70"N)     Center      (  33.7502456,   2.5001932) ( 33d45'0.88"E,  2d30'0.70"N)     Band 1 Block=64x64 Type=UInt16, ColorInterp=Undefined       Description = Layer_1       Min=0.000 Max=255.000        Minimum=0.000, Maximum=255.000, Mean=225.572, StdDev=31.338       NoData Value=256       Overviews: 1915x1276, 958x638, 479x319, 240x160       Metadata:         STATISTICS_COVARIANCES=982.046346716112,845.521289084631,710.510646703791         STATISTICS_MINIMUM=0         STATISTICS_MAXIMUM=255         STATISTICS_MEAN=225.57242379015         STATISTICS_MEDIAN=0         STATISTICS_MODE=0         STATISTICS_STDDEV=31.337618714831         STATISTICS_SKIPFACTORX=1         STATISTICS_SKIPFACTORY=1         STATISTICS_EXCLUDEDVALUES=         LAYER_TYPE=athematic       Image Structure Metadata:         COMPRESSION=RLE     Band 2 Block=64x64 Type=UInt16, ColorInterp=Undefined       Description = Layer_2       Min=4.000 Max=255.000        Minimum=4.000, Maximum=255.000, Mean=226.079, StdDev=29.739       NoData Value=256       Overviews: 1915x1276, 958x638, 479x319, 240x160       Metadata:         STATISTICS_COVARIANCES=845.521289084631,884.412989496947,765.062067015664         STATISTICS_MINIMUM=4         STATISTICS_MAXIMUM=255         STATISTICS_MEAN=226.07921779459         STATISTICS_MEDIAN=0         STATISTICS_MODE=0         STATISTICS_STDDEV=29.739081853631         STATISTICS_SKIPFACTORX=1         STATISTICS_SKIPFACTORY=1         STATISTICS_EXCLUDEDVALUES=         LAYER_TYPE=athematic       Image Structure Metadata:         COMPRESSION=RLE     Band 3 Block=64x64 Type=UInt16, ColorInterp=Undefined       Description = Layer_3       Min=0.000 Max=255.000        Minimum=0.000, Maximum=255.000, Mean=172.929, StdDev=27.135       NoData Value=256       Overviews: 1915x1276, 958x638, 479x319, 240x160       Metadata:         STATISTICS_COVARIANCES=710.510646703791,765.062067015664,736.335138826829         STATISTICS_MINIMUM=0         STATISTICS_MAXIMUM=255         STATISTICS_MEAN=172.92884721066         STATISTICS_MEDIAN=0         STATISTICS_MODE=0         STATISTICS_STDDEV=27.135495920046         STATISTICS_SKIPFACTORX=1         STATISTICS_SKIPFACTORY=1         STATISTICS_EXCLUDEDVALUES=         LAYER_TYPE=athematic       Image Structure Metadata:         COMPRESSION=RLE      Is there another hack ? Thanks,   Allan.


title: ArcPy: Feature class "doesn't exist" after MakeFeatureLayer

text: I was working on a Python script to analyze the geometries of a feature class and write something out. It worked fine the first time, but gave an IOError each subsequent time. I narrowed the problem down to the following: (not my actual code, but shows the weird behavior)               arcpy.Exists("MyFeatureClass")  # returns true     arcpy.MakeFeatureLayer_management("MyFeatureClass", "MyFeatureLayer")     arcpy.Delete_management("MyFeatureLayer")  # returns true     arcpy.Exists("MyFeatureClass")  # returns FALSE      Can anyone think of why the feature class would become inaccessible after making a feature layer out of it? Something to do with locks?


title: Custom Editor Extension - how to maintain original and the new feature

text: I am editing file geodatabase using edit toolbar. I've writing an extension to listen to edit events. It gives me the new feature with OnChangeFeature. However I am trying to obtain the original feature as well as a log of edits made. I've tried OnStartOperation but IWorkspaceEdit::EditDataChanges some how comes back null for file gdb. I could try to listen to OnSelectionChanged event but sometimes features attributes can be edited without selecting the row. Ideas?


title: Is it possible to change spatial reference (and reproject) feature classes in-place, using ArcMap/GeoProcessing?

text: We have an Esri based system in production; a complicated data model with a large number of feature classes. We need to migrate the data to another spatial reference. The feature classes have class extensions and relationship classes defined on them, so the behaviour of creating new feature classes for the reprojection result that the standard projection geoprocessing tool exhibits is very cumbersome (class extensions and relationships are lost). If necessary, we'll code up a solution but would prefer not to. Is it somehow possible, using standard or 3rd party extensions, to do a bulk, _in-place_ reprojection and setting of new spatial reference? Thanks!


title: Restart IIS Admin Service or Recycle Default App Pool?

text: A weekly maintenance process requires stopping ArcGIS Server and ArcSDE services. After restarting all the services, is it better to restart IIS Admin Service or recycle Default App Pool? I currently run a python script to do the later, with inconsistent results. Sometimes it takes 2 or 4 attempts, or stopping and starting AGS and ArcSDE services again, to get the web services working. I'd like a clear procedure that will work each time.


title: How to prepare domain on multifeature using Arcmap?

text: I am using in Personal Geodatabase. I want to prepare more than one field domain setting.


title: Can we request new OpenLayers WMS image for an extent without moving the map?

text: Can we request new OpenLayers WMS image for an extent without moving the map center ? Reading OpenLayers code, OpenLayers move map center to new extent, and THEN do HTTP request to get the image. Is it possible to instruct the map to HTTP load certain extent and when 'tileloaded' is fired, THEN load new WMS image? Edit ----------------------   1. The problem is I need to "tell the map to zoomToExtent" so it gets a new bound.   2. OpenLayers Map sets its center and zoom first, and then request new map through HTTP.   3. The map response will be available in 2 seconds, until that time, map view is blank (if new bound is outside current bound). if( result.zoom ) {                                     var bounds = result.data.map.mapBounds;                             var x1 = parseFloat( bounds.x1 );                             var x2 = parseFloat( bounds.x2 );                             var y1 = parseFloat( bounds.y1 );                             var y2 = parseFloat( bounds.y2 );                                  map.zoomToExtent([x1, y1, x2, y2]);                     }        4. What i think i need is fetchWmsBound() below.                   if( result.zoom ) {                            var bounds = result.data.map.mapBounds;                                  var x1 = parseFloat( bounds.x1 );                             var x2 = parseFloat( bounds.x2 );                             var y1 = parseFloat( bounds.y1 );                             var y2 = parseFloat( bounds.y2 );                                  function wmsLoaded() {                                  map.redraw(true);                             }                                  //                             // @param { <Array.number> } bounds,                             // @param { Function } callback to exec when wms is fetched                             // fetchWmsBound might trigger a 'fetched' event also                                  map.wmslayer.fetchWmsBound( [x1, y1, x2, y2], wmsLoaded );                                                          }      Filed at https://github.com/openlayers/openlayers/issues/573


title: How to display 33k placemarks with custom symbols/labels?

text: I want to overlay address numbers on our embedded Google Map (V3 api). Even though the KML is under 2 MB, it exceeds the 1,000 feature limit; there are over 33,000 placemarks. Fusion tables solve the size restrictions, but do not support custom symbols. Clicking on a marker to display the address number in a bubble is not acceptable. Did I miss something obvious that supports creating labels? I can see the address numbers in the Name column... Are custom image overlays the answer? http://oregonstate.edu/campusmap/ http://www.mail-archive.com/google-maps-js- api-v3@googlegroups.com/msg09037.html Please share tips on how to create the png image from shapefile vector points. I have access to ArcGIS 10 if that's an option. Thanks!


title: Unable to edit PostgreSQL/PostGIS layers in ArcMap 10.2.2

text: I am connected to the database, have added a layer, started edit session, but I am unable to edit the attributeeo table of said layer. ArcMap tells me that it is not an editable layer. When attempting the same procedure in QGIS, I am able to make any edits I need to. Does the database need to be modified to allow specific users read/write ability?


title: What GPS device do I need for vehicle tracking?

text: I am a programmer with no connection to Geographic information, except that I work on application which uses google maps. So here is my problem. I need some device, which I will be able to set in a vehicle in order to keep track of it. Every second (it would be better if this parameter can be configurable) a device should send information to the server:               GET/POST request to http://mysite/getLocation      with the following information:               lat: 23.1232     lng: 34.1232     deviceId: 321      Few important things:   1. information should be sent in real time (just storing info and accessing it in the end of the day does not work)    2. vehicle will be operating in the city   3. device should work at least 24 hours till recharging My question is: what device should I buy: If possible please provide some links to the device or where I can read about it.


title: How can I draw an ellipse filling the dataframe with ArcPy?

text: I need to modify one of my template scripts to draw a circle from the centre of the (inset) data frame to the edge of the data frame. This will then be used to clip the data frame to shape. I'm using ArcGIS 10.0.


title: ArcObjects/VBA to python Network edge to junction lookup IEdgeFeature equivalent in arcpy

text: In an effort to migrate a script that transfers attributes from a geometric network to a feature class (on a given network edge segment, what are the up and downstream junctions?) I would like to determine what the equivalent function in python for the following VBA/.Net IEdgeFeature Interface:   * FromJunctionFeature    * ToJunctionFeature Thank you!


title: QGIS-QT designer. How to create a form that allows editing elements linked by one to many relationship

text: I know this question has been already asked on this forum, but I really have not found any answer, not here, not anywhere. My question is simple, I have a Spatialite DB having one spatial layer (Pole) connected to a table (equipment) via a 1-N relationship. I created a custom form using QT designer which allows me to edit the elements of the Pole layer. How can create a sub-form or how to do to allow edit the elements of the connected table in the original form. Is it possible? Thanks in advance, Yuliana


title: Can I access a raster attribute table in QGIS?

text: I have some Python code that makes Raster Attribute Tables for GeoTiffs using GDAL. In ArcGIS I can see these attribute tables fine and they behave how I like (as classifications). But in QGIS the whole concept of an attribute table seems to be missing. Is there any way to access a RAT in QGIS?


title: Can I get gml:id from GetFeature request

text: I want to make a call to geoserver like `http://someurl:8080/geoserver/wfs?service=wfs&request=GetFeature --&someParameters...--` and get from geoserver just items I wrote in the parameters. For example I would like to have only gml:id of all features. Is it possible? Thank you in advance.


title: Where are QGIS application setting file(s) stored?

text: I have some application settings (Windows/QGIS v2.0) that I would like to pass to other users (such as those settings defined under Settings menu>Options>Rendering tab). In ArcGIS they are stored in a **Normal.mxt** file located: > C:\Users\userName\AppData\Roaming\ESRI\Desktop10.2\ArcMap\Templates I suspect in QGIS they may be stored in **qgis.db** located here: > C:\Users\userName.qgis2 But I do not know for sure. Where are these application settings file(s) stored? thank you


title: differentiating stone houses from barren rock in ArcGIS 9.3

text: I'm working in ArcGIS 9.3 trying to differentiate stone houses (often red brick roofs) from barren rock in a typical Mediterranean karst region via remote sensing. I have been using Landsat imagery. I have been trying out various band combinations, multiplying and filtering bands in unsupervised classification. So far no good. Does anyone have some helpful suggestion? Thanx in advance :)


title: bulk change of certain nodes in OSM file

text: How to I change a OSM file in multiple nodes and just change the value i.e. aminity=townhall into aminity=school? Please suggest.


title: How to convert a polygon to a grid of polygons using QGIS?

text: For example, the city of Rotterdam represented as a hexagonal grid. ![Rotterdam hexagon urban identity](http://i.stack.imgur.com/60uM6.gif)


title: Which program that merges excel data onto shapefiles is easiest to use?

text: some options for merging looks like fusion tables, vizmaps, and these 2   http://mapbox.com/tilemill/docs/guides/joining-data/   http://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=About%20joining%20and%20relating%20tables.. which is easiest? the shapefiles are probably going to be TIGER census ones unless I find better ones.


title: Can you explain this WKT multipolygon format?

text: I am confused by the output I'm getting from ogr2ogr when converting a shapefile to WKT. What I'm expecting to see is something in the form below (from Wikipedia):               MULTIPOLYGON (((40 40, 20 45, 45 30, 40 40)),     ((20 35, 45 20, 30 5, 10 10, 10 30, 20 35),     (30 20, 20 25, 20 15, 30 20)))      What I'm getting out of `ogr2ogr -f CSV out.wkt source.shp -lco GEOMETRY=AS_WKT`, in very abbreviated form, is this:               "MULTIPOLYGON (((-120.6 50.8, ...)", 33.9, 137.1, 2, 1,1,1      I do not understand how to interpret the final six numbers, which are not enclosed in brackets and fall outside the quotation marks. Can someone explain this or point me towards the appropriate documentation?


title: How can a HTML table be stored within a database field as a BLOB, text or other value?

text: My Spatialite database has several thousand spatial polygon records In QGIS, each with a couple dozen field values. I also have many thousands more of additional tabular data relevant to those polygons, all stored within Spatialite, that the polygons need, a one to many relationship (a table as a value) I'd like to create a QGIS ACTION of a single Field value that pops up a whole table stored within a single value in the spatial record. My initial strategy is to use a BLOB or Text to store HTML coded tables I create, so that the QGIS ACTION would open it up in a web browser. Opening up an equivalent table in Excel or some other common application would be satisfactory. My worries are that I can't store html as a BLOB since its not binary. (but maybe that's not relevant). Also, I figure I can't store the html table as text since it would contain more characters of data than spatialite or shapefile can handle. I can create an html coded page for each field value without a problem, but...   1. How would I store the html page as a single field value in Spatialite or Shapefile?   2. Or, what would be a good alternative to doing that, such as using another format besides html or BLOB data.   3. Would it be possible to have a SQL query stored as a field value, and assign an Action that pops open the Results? I'm understanding I could use PostGIS to create a database relationship if all else fails. It may be a day before I get back with any feedback. Thanks in advance for everyone's suggestions.


title: Is it possible to change the Mouse Wheel zoom direction in QGIS?

text: I'm just starting to learn QGIS. I'm starting with the basic simple things. Coming from an ArcGIS background, I am used to the mouse wheel action zooming in and out the default ArcGIS way. QGIS default is the opposite. The only option I seem to find about changing mouse action options is to enable or disable the mouse wheel action, but nothing about the altering the action. Can I do this?


title: How to write a CSV from Python console in QGIS?

text: I am trying to write out a CSV file from a vector layer made up of points using the python console in QGIS. I've tried the following command:               QgsVectorFileWriter.writeAsVectorFormat(layer_Map, Folder + "my_shapes.csv", "CP1250", None, "delimitedtext")      but it does not like the `delimitedtext` parameter (if I use `ESRI Shapefile` it works fine). I have not been able to determine what other options can be used instead of `ESRI Shapefile`. I find the Quantum GIS API documentation to be quite confusing, but I hope I will get used to it soon. Thanks in advance for any help.


title: How to view error feedback in qgis when testing connection to postgres database?

text: Using qgis 1.8 on Ubuntu 10.04 with Postgres 8.4.13 I am creating a new (and my first) Postgresql connection in qgis. I can connect to the user and the db using the pgsql client.               jlam@lusk:~/programmation/Quantum-GIS/build-master/output/bin$ psql -U qgis -W qgis -p 5432     Mot de passe pour l'utilisateur qgis :      psql (8.4.13)     Saisissez « help » pour l'aide.          qgis=> \q     jlam@lusk:~/programmation/Quantum-GIS/build-master/output/bin$      But when testing the same parameters within the PostGIS connection creation dialog box, I get connection failure:               "La connexion échouée - Réglez la connexion et réessayez"          ("Connection failed - Fix the connection and try again")      My pg_hba.conf file includes the line:               # "local" is for Unix domain socket connections only     local   all         qgis                          md5       Is there a way I can get qgis to send me back the error code and error string it gets from postgres, if its connecting at all, or to increase debug output on the console?


title: 32bit to 64bit topology errors

text: Can generating a shapefile from a pgdb cause topology errors due to floating point conversion? I am running into gaps in a line layer of really small distances (order of 10^-4; CRS EPSG:2227 usft).


title: How can I get the geometry length in meters?

text: The ST_length() function returns the length of a linestring but not in meters. Is there a function like ST_Length() but in meters?. Thanks!


title: How to run all rulesets to get a final image which include all the results of classification together

text: How to run a all rule-sets to get final image which include all result of classification together. Once I execute on 'process tree window', it gives me a result of individual classes. I need an image with 10 classes together on it. I can't find any tool to process it from drop-down list on main screen. Thanks a lot.


title: How to fix ImportError: cannot import name MAXREPEAT on Ubuntu 12.04

text: I have problem with installing qgis 2.0 on ubuntu 12.04. If I don't install qgis-python plugin then qgis starts normally but I can't access python console in plugins. If I install that plugin then I can't even start qgis because it gives me error ImportError: cannot import name MAXREPEAT I believe someone have already had this problem so I'd appreciate if you could tell me how to fix it. Thanks


title: Sourcing GIS data from First World War (WW1)?

text: Simple question, but does anyone know where I can get some detailed (eg 1:10m scale) vector data from the First World War? I've been through the archives here, but it seems the best sites for it have already closed.


title: How to convert decision tree to ruleset in eCognition?

text: I have created a decision tree which created using decision tree tool in ENVI to classify 10 land cover types. Now I want to convert it to a ruleset in eCognition. How can I convert it? Thank you very much.


title: WMS GetCapabilities and variable substitution with MapServer

text: I have created a layer with a variable substitution to request wether or not a bird species has been observed in a town. It works ok. But if I make a WMS GetCapabilities I only see one layer as expected. I have more than hundred species so if I want to have the GetCapabilities request to show one layer for each species do I need to generate a huge mapfile? Does anyone can suggest some tips/trick to make this work easier?


title: How to create GoogleMaps/OSM hyperlink for given bbox coordinates

text: I have extracted the bbox coordinates from a shapefile and want to create a hyperlink for google maps respectively OSM. When you click on the hyperlink the map should be zoomed to the bbox extend and most perfectly the bounding rectangle should be rendered. Do you know how to do it? Thanks!