title: Openlayers. How to remove last sketch if invalid

text: I can currently identify if a polygon is invalid through wps as per below code. What I want to know is how to remove the last sketch polygon if it returns invalid? To give you an idea, here is the code that is run on "sketchcomplete".                   farmpaddocks.events.on({         sketchcomplete: function (e) {         var poly = e.feature.geometry;//var of sketch geom              //check if geom valid through wps         wpsClient.execute({         server: 'local', //from server address         process: 'JTS:isValid', //wps function         inputs: {         geom: poly  //geom input for testing if valid         },         success: function(outputs) {          //on completion of testing check if false = invalid   geom         if(outputs.result==false){         $('#noti').append("Invalid Polygon").show(1000);         $('#noti').css("color", "red");         $('#noti').fadeOut(5000);         draw.deactivate();         }         }                        });


title: How to Automatically Update ArcGIS.com Map Layer Data

text: We'd like to refresh layers in ArcGIS Online every five minutes with location data modifications occurring in our externally hosted SQL database. We'd like to avoid purchasing ArcGIS for Server. Any pointers for using REST to make the updates? Can we update multiple locations simultaneously? Can we push a JSON feed containing changes and deletions to ArcGIS Online every five minutes?


title: I cannot use Gdal Python but Gdal binary, how to fix?

text: **UPDATE: I can use GDAL with Python now by copy file from Gdal 1.11 compiled folder to Python package folder. But the problem with GDAL binary and GDAL library that they cannot read the HDF5 file. How can I fix this?**               ERROR 4: `GAERO-VAOOO_npp_d20130201_t0549417_e0555221_b06552_c20140121145913139281_noaa_ops.h5' not recognised as a supported file format.          gdalinfo failed - unable to open 'GAERO-VAOOO_npp_d20130201_t0549417_e0555221_b06552_c20140121145913139281_noaa_ops.h5'.      I've installed GDAL binary in Linux Ubuntu 13.10. Every GDAL binary has been run correctly (ex: gdalinfo, gdal_translate,...) > gdalinfo --version GDAL 1.11.0, released 2014/04/16 But the problem is that I cannot include GDAL library to Python script. It's error like this:                   >>> import gdal         Traceback (most recent call last):         File "<stdin>", line 1, in <module>         File "/usr/local/lib/python2.7/dist-packages/gdal.py", line 2, in <module>         from osgeo.gdal import deprecation_warn         File "/usr/local/lib/python2.7/dist-packages/osgeo/gdal.py", line 83, in <module>         from gdalconst import *         File "/usr/local/lib/python2.7/dist-packages/osgeo/gdalconst.py", line 108, in <module>         GRA_Average = _gdalconst.GRA_Average         AttributeError: 'module' object has no attribute 'GRA_Average'      How can fix this problem, manythanks for help me to install GDAL correctly or something to give me a chance to import GDAL to Python script. Thanks! NOTE: I installed GDAL by using wget and sudo make, sudo make install from source code. **Update** It look like that I compile too much GDAL version (1.9, 1.10 and newest 1.11 then now Python confused to load library). Please help:               locate _gdalconst.so     /home/rasdaman/TEMP_FOLDER/gdal-1.10.1/swig/python/build/lib.linux-i686-2.7      /osgeo/_gdalconst.so     /home/rasdaman/TEMP_FOLDER/gdal-1.9.0/swig/python/build/lib.linux-i686-2.7/osgeo/_gdalconst.so     /home/rasdaman/gdal-1.11.0/swig/python/build/lib.linux-i686-2.7/osgeo   /_gdalconst.so     /home/rasdaman/gdal-1.9.2/swig/python/build/lib.linux-i686-2.7/osgeo/_gdalconst.so     /usr/local/lib/python2.7/dist-packages/osgeo/_gdalconst.so     /usr/local/src/gdal-1.9.0/swig/python/build/lib.linux-i686-2.7/osgeo/_gdalconst.so


title: IDW: Default Maximum Distance and Coincident Points

text: The ArcGIS Help file says the Default Maximum Distance, when using the Variable Search Radius in IDW, is the extent's diagonal, however, no value is automatically filled in the corresponding box after I define the extent, so is ArcMap really using this default value if the box is empty? Regarding coincident points (different measurements at the same location) when running IDW, it is possible to set the Environment "Geostatistical Analysis" to Min, Max, Mean, etc., but is this selection really honored when running IDW using Spatial Analyst, or does it only apply to tools from the Geostatistical Analyst extension? I am confused because the ArcGIS help says: "Tools that honor the Coincident Points environment define how coincident data is treated in Geostatistical Analyst". If the Coicident Point selection is not honored when using Spatial Analyst, how do I know what IDW is doing about them? Thanks a lot for your help!


title: Polyline ZM: automatic vertex height errors correction

text: Paolo I've got many rivers (polylines ZM) with a lot of height errors in Polylines 'vertices. Vertices' Z errors are not all equal in height (i.e. they are not all at zero height). I must create a TIN using these rivers as breaklines. How can I automate this kind of correction? CAn I use an intermediate DTM (without the contribution of my river/breaklines) to solve my problem? How? I was thinking to use "Add Z information" and "feature to point" to transform each polyline's vertex in a point with it's Z value... then.. I would have to compare all the points with their neighboring and identify those with very large Z difference... and then, assign a medium Z value. Could it work? Some suggestions? Thanks in advance PAolo


title: Load a feature from geoserver using WFS protocol on localhost

text: I am a student enginner in GIS , I want to load my features from geoserver using WFS protocol. I'm using this code but it's not appearing in my browser                   OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";     var map,wfs;               function init() {     //-11593508, 5509847, -11505759, 5557774         var extent = new OpenLayers.Bounds(             474813,  211897, 491713, 227436.5         );     map = new OpenLayers.Map('map', {         controls: [             new OpenLayers.Control.PanZoom(),             new OpenLayers.Control.Navigation(),             new OpenLayers.Control.LayerSwitcher({'ascending':false}),         ]     });       var wms = new OpenLayers.Layer.WMS(                         "test_postgres:territoire", "http://mylocalhost:8080/geoserver/test_postgres/wms",                         {                             LAYERS: 'test_postgres:territoire',                             STYLES: '',                 format: 'image/svg',                         },                         {                             isBaseLayer: true,                         }                      );       var wfsLayer = new OpenLayers.Layer.Vector("WFS", {             strategies: [new OpenLayers.Strategy.BBOX()],             styleMap: new OpenLayers.StyleMap({                 strokeWidth: 3,                 strokeColor: "#333333"             })          });                       map.addLayer(wfsLayer);               map.zoomToMaxExtent();


title: Arcgis, Convert analyse result in numeric value

text: I have a problem that seems really simple but blocking me... in Arcgis i have script python where i want to extract the resolution of a raster and use this result like numeric value for use in a mathematical function. But when i try to use it, here i have my problem because it is not a numeric value. So i want to convert the result of my request (Get raster Properties) in numeric value. Some body have a solution? > > > a = arcpy.GetRasterProperties_management(r"Topographie\Raster > d'élévation \Lidar\Elevation LastReturns.tif","CELLSIZEX") Result: Cellsize > in x direction = 1.000000 >>> >>> type(a) Result: class 'arcpy.arcobjects.arcobjects.Result'               and if i try :               File "<string>", line 1, in <module>               TypeError: float() argument must be a string or a number


title: Difficulty converting standalone script to ArcGIS tool

text: As my title says, I am currently trying to make my (working) standalone script into an ArcGIS tool so I'm not the only person in my office that can use it. It essentially selects all the data within a specified location and copies it to a new geodatabase. Since replacing my old variables with sys.argv[#], new errors keep popping up. The following section of code is giving me problems right now:               out_folder_path_gdb = sys.argv[1]     out_name_gdb = sys.argv[2]     out_name_fc = "buffer_point"     geometry_type = "POINT"     has_m = "DISABLED"     has_z = "DISABLED"     sr = sys.argv[3]     fc = out_folder_path_gdb + "\\" + out_name_gdb + "\\" + out_name_fc     lyrlst = []          ######################     # Create geodatabase #     ######################          try:         arcpy.CreateFileGDB_management(out_folder_path_gdb, out_name_gdb)     except:         print "Geodatabase already exists."         arcpy.AddIDMessage("Error", 12, out_name_gdb)         sys.exit(0)          ##########################     # Create point to buffer #     ##########################          arcpy.env.overwriteOutput = True     arcpy.CreateFeatureclass_management(out_folder_path_gdb + "\\" + out_name_gdb, out_name_fc, geometry_type, "", has_m, has_z, sr)          ############################################     # Sets given coordinates as a point object #     ############################################          rowInserter = arcpy.InsertCursor(fc)     x = sys.argv[4]     y = sys.argv[5]     pointGeometry = arcpy.Point(x,y)      When running the tool within Arc, I get the following error: : Failed to execute. Parameters are not valid. ERROR 000732: Feature Class Location: Dataset C:\Documents and Settings\benoyn\Desktop\DaS GIS work\data_extraction_test\pleasepleaseplease does not exist or is not supported Failed to execute (CreateFeatureclass). The only reason I can think of for this error is the spaces in "Documents and Settings". Thanks everyone!


title: Free Imagery WMS Feeds in QGIS

text: I was wondering if anyone knows of a free WMS of imagery that I could use within QGIS. The imagery does not have to be high resolution (Blue Marble even). I would like to have a base imagery feed from a WMS. Any help would be appreciated. Thanks...


title: Where to download USGS DOQQ boundaries?

text: I downloaded a shapefile of USGS DOQQ boundaries a few years ago and cannot seem to find one now. The attached image shows an example of the shapefile I'm looking for (green polygons showing the quarter quad boundaries). Where can I download a shapefile covering the contiguous USA that shows USGS DOQQ boundaries with relevant attributes? ![enter image description here](http://i.stack.imgur.com/asIjl.jpg)


title: ESRI ArcGIS JSAPI 3.4 Defect when setting LODS on a map. How to fix?

text: I tend to get better help here then the ESRI forum so... Using http://developers.arcgis.com/en/javascript/sandbox/sandbox.html I was able to create a simple application which demonstrates that if you set a non- trivial lods on a map the basemap will not behave correctly (try zooming beyond level 15). Any workarounds?               <!DOCTYPE html>                 <title>Create a Map</title>       <meta http-equiv="Content-Type" content="text/html; charset=utf-8">       <meta http-equiv="X-UA-Compatible" content="IE=7, IE=9, IE=10">       <link rel="stylesheet" href="http://serverapi.arcgisonline.com/jsapi/arcgis/3.4/js/dojo/dijit/themes/claro/claro.css">       <link rel="stylesheet" href="http://serverapi.arcgisonline.com/jsapi/arcgis/3.4/js/esri/css/esri.css">


title: CartoDB toggle map view - can I share the map as an iFrame?

text: I followed the tutorial - toggle map view and have successfully made the map. I'm now wondering what is the best way to share this map? I'd like to be able to give a client an iFrame, can I do this after I've edited the code?


title: QGIS Bathymetry Display, Bathymetry along a path

text: This post has two related questions. I am new to QGIS or any GIS software so please forgive me if the answer is obvious. I wish to view the Bathymetry point for any given coordinates (e.g. as I move the cursor, coordinates are displayed in the status bar) I would like to see next to coordinates in that status bar the bathymetry data. I created a layer using WMS to connect to Global Multi-Resolution Topography @ http://www.marine- geo.org/services/wms_2.0_merc?request=GetCapabilities&service=WMS&version=1.0.0. Another application, GeoMapApp, provides the feature describd above but I need QGIS to draw shapefiles (e.g. paths) that meet certain bathymetric requirements. Finally, I wish to obtain a bathymetry profile a long the selected path (Google Earth calls this feature "elevation profile". Can QGIS provide this? Many thanks


title: BBox missing when converting TIGER to GeoJSON

text: I can't get the bounding boxes of polygons if TIGER/Line shape files when converting them to GeoJSON using FWTools. The other default information is present and there is no difference between the GeoJSON files created with or without -lco WRITE_BBOX=YES. Here is the command line I used: ogr2ogr -f GeoJSON myF.json myF.shp -lco WRITE_BBOX=YES


title: Why are z-values replaced with zero after in_memory use in Python?

text: My script takes separate 3D feature classes, calculates a field, appends to a target feature class, and then deletes the originals. This is a huge data processing task of over 4700 mass point feature classes appending to a single one. I know it was a big mistake to put so many feature classes into SDE, but given they're in there now, I've got to deal with them there (yes, the sheer fc count is bringing SDE performance to its knees--hence my urgency). Early on it worked great, in that the merged result showed valid z-values. I was able to build a terrain and it looked great. I'd been running this in the background for weeks, and have just noticed z-values of 0 in the output, where the originals were not 0! Oh no!!! I'm certainly glad I still have the file- based originals! Pseudocode: > An in_memory feature class is created from a z-aware template > > For each fc in the feature dataset: > >> It is appended to the empty in_memory fc >> >> The TILE_ID field is calculated using a portion of the name of the source fc >> >> It is then appended to the target SDE fc >> >> The in_memory feature class is emptied > > Each processed fc is then deleted EDIT: When I initially posted, I was curious if it was my use of in_memory that caused this trouble. I had not initially been using it, and only switched to it because performance was far too slow otherwise, and WAY faster in_memory. To find where trouble started, I ran the Add XY Coordinates tool, which proved the trouble started with in_memory use. At that time I had also switched from CalculateField_management to da.UpdateCursor, so that's likely the bigger reason for the performance enhancement. I will now change my code to avoid in_memory. Thank you to eric and PolyGeo.


title: Imported GPS data does not overlap with map

text: I´m working on QGIS and I geocoded 2 addresses of a city in Germany and tried to input the text-file as a layer. It worked, but the 2 points must be in the map of the city and if I zoom to the layer-file they are far away. I tried EPSG 3857 (like the other layers) and ESPG 4326. With the last one it was much nearer, but still not in the city-map. On-the-fly-projection is activated. Does anybody know, how I get the points into the right place in my map? Am I using the wrong projection? Thank you.


title: How to force ArcGIS iOS to use GPS data from External GPS device?

text: I created an app that will display the map and the location of the external GPS device. The scenario is using a wifly module to broadcast the GPS data of the GPS receiver in wifi. The app is now receiving the data from the wifly module. My question is: How do i force ArcGIS to use the data from the wifly module? Im using the AGSLocationDisplay of ArcGIS iOS but its using the data from the CLLocation.


title: styles for QGIS

text: Could someone point me to a source of line and point styles for QGIS (*.qml) - railways and transmission lines are the first two I've run up against, but I know there will be more.


title: GPS units compatible with Galileo

text: It seems like after lots of initial problems Galileo project is taking off (ESA, BBC). I was wondering what the state of preparation for that is on hardware side. Are there any consumer grade GPS units that are ready to take advantage of Galileo?


title: Modifying Fields in QGIS

text: I'm pretty new to Quantum GIS, and I was looking to add zeros to a field I was having some trouble with. I learned through browsing a few posts here that existing fields can be modified by defining actions. My question is, how does this process work? How do I go about defining actions and what is the language (for example here) used to define these actions? Example: `substr( ( '000' || tostring ("myfield")) , ( length ( ( '000' || tostring("myfield")) ) -2) , 3 )`


title: Google Map - display specific markers popup using a list of markers

text:
I have the following code:                     What I want to do is dynamically add the list item to the page in divs then allow the user to click each div the relevent infoWindow will display on the map. I have done this with openlayers but now need to do it with Google maps. I have got as far as adding the the list to page with the following:               // add items                   var  available;                     var  propertyimg;                     var  location = "http://cherrybird.fourstudio.co.uk/" ;                 var listItem= document.createElement('div')                     listItem.className = "propertyhot";                 listItem.innerHTML = ''  +'
' + datanew.property_name +'
' + datanew.property_number +' ' + datanew.street_name +'
' + datanew.county +'
' + datanew.postcode +'
' + datanew.rooms +' Rooms

 Available from:' + datanew.availability + '
';                 list.appendChild(listItem);      The next step is to get it to display the infoWindow for the relevent marker when it is clicked in the list. Can anyone help with this.


title: How to programatically modify the placement properties for labels

text: Does anyone know how to modify the location of a label programatically? I'd like to make the labels sit directly below the symbol!


title: remove duplicate points based on a specified distance

text: i have a table of points with each point acquired on a specified date, by a specified source, and having a measured value. i want to consider points as duplicates if they have the same date and if they are within a certain distance of each other (say <10 m), even if they came from different sources. here's a brief example using 8 points:               create table points (gid serial primary key, date date, source varchar, value integer);          select addGeometryColumn('points', 'geom', 32633, 'point', 2);          insert into points (date, source, value, geom) values       ('1985-06-19', 'B', 22, st_geomFromText('POINT(593835 3910494)', 32633)),       ('1985-06-19', 'A', 48, st_geomFromText('POINT(602920 3910591)', 32633)),       ('1987-04-17', 'B', 12, st_geomFromText('POINT(592317 3910479)', 32633)),       ('1993-07-22', 'A', 16, st_geomFromText('POINT(592323 3910476)', 32633)),       ('1993-07-22', 'B', 30, st_geomFromText('POINT(592317 3910479)', 32633)),       ('2000-07-09', 'A', 20, st_geomFromText('POINT(592108 3909121)', 32633)),       ('2006-07-14', 'A', 22, st_geomFromText('POINT(592933 3909842)', 32633)),       ('2006-07-14', 'B', 18, st_geomFromText('POINT(592931 3909843)', 32633));      i'm trying to create a query that should result in 6 points, since points 4/5 and 7/8 would be considered duplicates because they have the same date and <10 m from each other. i also want the query to select the duplicate point that has the maximum value, so in my example the result set would be ('geom' field not shown):                gid |    date    | source | value      -----+------------+--------+-------        1 | 1985-06-19 | B      |    22        2 | 1985-06-19 | A      |    48        3 | 1987-04-17 | B      |    12        5 | 1993-07-22 | B      |    30        6 | 2000-07-09 | A      |    20        7 | 2006-07-14 | A      |    22     (6 rows)      been trying various queries using self joining and the st_dwithin and st_distance functions, but can't seem to nail it down. i'm using postgis 2.0 and postgresql 9.1. thanks! update: based on the link whuber suggested, i explored using the st_snapToGrid function and this query seems to produce the result i'm looking for:               select distinct on (date, grid_location) grouped.*     from     (select *, st_snapToGrid(geom, 10, 10) as grid_location     from points order by date, grid_location, value desc) as grouped;      update 2: the above query eliminated most of the duplicates in my actual dataset, but i found several still making it through since sometimes points would be snapped to different grids if they happened to be bisected by a grid boundary. i've posted another solution as an answer to my question.


title: Transformation data austria - slovenia

text: can anyone please help with fitting together these 2 shape-files? many many thanx https://www.dropbox.com/s/saa9yp7ecg1qxh1/Grenzen.rar


title: How to remove popup window from KML layer in ArcGIS API for JavaScript?

text: I am using the ArcGIS API for JavaScript for creating a map. I have added some CSV and KML layers to the map. Now the KML is colored base layer(different colored for country) on the map and CSV is some green dot layer on map but top on the KML layer.When I click on the green dot(CSV) the green dot is selected and a popup window is shown on the for green dot(CSV layer).But on clicking the green dot the base KML layer is also selected and also give popup window for the KML layer on goto the arrow sing on the header of popup. Here is a link to my demo. You can check here click on the green dot and click on the arrow sing on the right of popup black header. Can anyone suggest how to remove the popup from the KML layer?


title: How to calculate exact area of a loaded jpg map in to arcmap?

text: I calculated the polygon area through arcmap.but i'm not sure that the area is not actual,may be scale is not right. so how to calculate the exact area of a jpg image?


title: Removing "the_geom" field in the .csv output of Geoserver

text: Good day sires. My team is currently working to develop a web-based GIS for a government client and we're using Geoserver 2.1.1 as our server application. We're trying to play around with Geoserver's GetFeature function to generate a .csv file of the attribute fields of features saved in a PostGIS database. The function below does exactly that: _function generateFeatureTable(feature){ var featureTableURL = "http://localhost:8080/geoserver/prism/wfs?REQUEST=GetFeature &SERVICE=WFS&VERSION=1.1.0&TYPENAME=prism:" + feature + "&OUTPUTFORMAT=csv&SRSNAME=epsg:4326"; window.open(featureTableURL); }_ The problem is that the .csv file is also including "the_geom" as one of its fields. If we make use of the PROPERTYNAME argument to specify the fields we like to include in the table, we would have to specify each and every field name excluding that of "the_geom". Any ideas on how to specify the ones you would NOT want to include in the output file? We will highly appreciate the insights you can provide us.


title: pgRouting: unable to execute routing_core.sql

text: I'm unable to execute the routing_core.sql, receiving following error. I'm using pgAdmin III, Windows 7, PostgreSQL 9.1, PostGIS 2.0. ERROR: could not load library "C:/Program Files/PostgreSQL/9.1/lib/librouting.dll": %1 is not a valid Win32 application. _**_*** _**Error **_** _*_** ERROR: could not load library "C:/Program Files/PostgreSQL/9.1/lib/librouting.dll": %1 is not a valid Win32 application. SQL state: 58P01 Regards, Nas


title: Spatial Database Connection String

text: From MS Access, I'd like to create a link to the attribute table of an SDE feature class in a spatial database. For example, I have the following connection string which works for an _OLE DB connection_ :               ODBC;DRIVER=SQL Server;SERVER=WCARCSDE00;Trusted_Connection=Yes;DATABASE=MyDb      I want the equivalent connection string for a _Spatial Database Connection_. Something like the following (which does not actually do what I want because the SERVICE parameter is ignored):               ODBC;DRIVER=SQL Server;SERVER=WCARCSDE00;SERVICE=sde:sqlserver:wcarcsde00;Trusted_Connection=Yes;DATABASE=MyDb      Is this possible? Do I have some fundamental misunderstanding about how this all works?


title: Public PostGIS Server available?

text: can anyone point me to a (if possible at all) list of public PostGIS Servers, where one can set up it's own database/tables ? It's just that I would need a PostGIS 24h on, while dont having these resources by myown, for a or in fact multiple humanitaerian projects in the south pacific.


title: Automate end-snap of lines mid-way between nodes?

text: I have a polyline feature class of field-collected (Trimble GeoXT with ArcPad10) stream segments. For whatever reason (still troubleshooting why), my employee returned from the field and the majority of the endpoints were not snapped (despite occurring within distances of the snapping distance I'd specified on the XT). To automate generation of a continuous series of stream segments, I already ran the "Snap" tool, but it only moved one endpoint and left the other one fixed. In this scenario, I cannot assume that one end-point is more accurate than the other, so I would like to average them. Is there a command (or maybe an environment setting within "Snap") I'm missing that would move endpoints for adjacent lines and snap them half-way along a straightline distance in between? Alternatively, a command that generates a point midway between ends that I could subsequently snap lines to would probably work as well. I have an ArcInfo license v10 SP4. Thanks, Will


title: Averaging satellite data with missing values (using Python)

text: I am using sea surface temperature, chlorophyll, and other data parameters as covariants in a model, and I would like to average them over a particular time period. I am currently doing this using a simple average of the values across the days in my time period, but this is flawed due to missing data from cloud-cover, etc. (and perhaps other reasons). `filledsst = np.average((datafile.variables["FilledSST"][idx1:idx2,:,:]), axis=0)` What technique is typically used for accounting for such missing data?


title: OpenLayers generate random polygon

text: I have a world map. I need randomly to generate polygons and render them on map. I have no problem with rendering, but can anybody suggest how to generate polygons with random geometry? (I need `n` polygons, `n` \- random number). The main problem is how to generate geometry to them.


title: USGS Topographic Map Online API

text: I am creating a website for a non-profit organization and would like to plot points on a map similar to this: http://www.arcgis.com/home/item.html?id=931d892ac7a843d7ba29d085e0433465 I am currently plotting points on google maps through its API, but it isn't a great solution and I have heard that Google has rights to that point. (Correct me if I'm wrong) Anyway, I was wondering if there was a good way to use the USGS Topo layer on the USGS website to plot points dynamically through javascript with an API or some other means... Is there anyway I can set this up through mapnik, similar to this (http://www.toposm.com/us/index.html)? On another note I have to be "licensced" in order to access ArcGIS javascript API, don't I? Or is this incorrect? Aaron


title: Problem compiling QGIS against trunk - georeferencer

text: I got problems compiling QGIS (trunk) this morning. The georeferencer seems to be the source of the problem:               [ 81%] Building CXX object src/plugins/georeferencer/CMakeFiles/georefplugin.dir/qgsgeorefplugingui.cpp.o     /usr/local/src/Quantum-GIS/src/plugins/georeferencer/qgsgeorefplugingui.cpp: In member function ‘void QgsGeorefPluginGui::createActions()’:     /usr/local/src/Quantum-GIS/src/plugins/georeferencer/qgsgeorefplugingui.cpp:864:3: error: ‘mActionLocalHistogramStretch’ was not declared in this scope     /usr/local/src/Quantum-GIS/src/plugins/georeferencer/qgsgeorefplugingui.cpp:868:3: error: ‘mActionFullHistogramStretch’ was not declared in this scope     /usr/local/src/Quantum-GIS/src/plugins/georeferencer/qgsgeorefplugingui.cpp: In member function ‘void QgsGeorefPluginGui::addRaster(QString)’:     /usr/local/src/Quantum-GIS/src/plugins/georeferencer/qgsgeorefplugingui.cpp:1124:3: error: ‘mActionLocalHistogramStretch’ was not declared in this scope     /usr/local/src/Quantum-GIS/src/plugins/georeferencer/qgsgeorefplugingui.cpp:1125:3: error: ‘mActionFullHistogramStretch’ was not declared in this scope     make[2]: *** [src/plugins/georeferencer/CMakeFiles/georefplugin.dir/qgsgeorefplugingui.cpp.o] Error 1     make[1]: *** [src/plugins/georeferencer/CMakeFiles/georefplugin.dir/all] Error 2     make: *** [all] Error 2      Bit weird as I've been successfully compiled QGIS trunk on this machine running Ubuntu 12.04. Any help would be most appreciated.


title: problem creating and editing features at extent edges

text: I'm in the process of creating a polyline feature layer of tectonic plates. I'm working with WGS-84 reference system, which means the extent is +-180 in lon, and +-90 in lat. In order to create the features in a correct way, I imported a CSV file with the 4 "corners" and displayed them as points I then snap my features to when I create them. The problems arises when I want to snap the line to the edge point. It simply doesn't want to create the line all the way, cutting it shy of the actual point. Here's what I see: ![enter image description here](http://i.stack.imgur.com/JZjsq.png) The coordinate is that of my cursor. The green point is (-180,-90). The green background is another layer I used (continents) to see if there's some kind of error or a general problem, but as you can see, it sits perfectly in place. I have tried different editor tools (line, polygon, autocomplete), but the problem persists. I have also tried to extend that segment the the point, as well as moving a node to the point, but it doesn't edit it, and jumps back to it's previous position. Does anyone know why this happens?


title: Copy, not export, a table from Arcmap

text: How can I **copy** a table, or selection of rows or columns in a table, from Arcmap _without_ going through the pointy-clicky-typey exercise of Export table? The destination of said copy could be anything: notepad, excel, word, whatever. The existing context menu only allows copying the value for a single cell and using the keyboard [Ctrl-C] also just copies a single cell value: ![context menu expanded to showing copy cell value](http://i.stack.imgur.com/VnkAX.png) I'm using Arcmap 10.1, but this has been annoying for years, so don't feel the need to restrict your answer to a particular version.


title: What are cloud-based alternatives to Google Fusion Tables?

text: I'm a beginner in web-mapping. I was playing around with Fusion Tables, but I realized that I had too much data to display. So here's my question, what are some alternatives to Fusion tables? I have 9 layers of data (one layer is parcel data so there's over 500000 records) and I have to be able to query the data (i.e. how many parcels are within 150m of a school etc.) Thanks for your help!


title: Generating very simple slippy map tiles

text: I have an app that uses slippy maps. The tiles I'm currently using contain a lot of information that is superfluous to my users. I've been investigating generating replacement map tiles that contain only very basic features such as country/province boundaries, rivers, lakes, EEZ lines (the application focuses more on water than land) and maybe railways and major motorways when at high zoom levels. First I tried loading `planet.osm.bz2` (a 21GB file) into Maperitive. This appeared to be a bad idea as the process memory usage rose quickly beyond that which my machine could tolerate. I've tried again on smaller files (such as Antarctica alone) and have more success, but I really need to generate tiles for the entire Earth. Most of the features in these large OSM data files are not of interest to me. Can someone please recommend an approach that should run on my x64 Win7 PC with 4GB RAM?


title: Most accurate GPS hardware for less than 1 mile radius

text: I apologize if this is a duplicate but I can't find a good answer or area to know where to find the correct answer. I am looking for the technology or hardware around wearables or similar technology that can detect movement or at least log location within about a 2 foot area. I have looked around at it looks like the most common devices use triangulation but aren't more accurate than around 400 feet. Is there any tech that is more specific than a two foot radius? I am not trying to waste anyone's time but I can't find a better place to get an answer....


title: What GPS unit can I use to track myself all day for a 10 day trip?

text: > **Possible Duplicate:**   >  Logging GPS and Accelerometer Data - Device Suggestion? I will be visiting China this October as part of a field studies course and I want to be able to take the opportunity to try out something unusual (maybe usual?). I want to be able to track my movement throughout China using the GPS unit attached to my bag/clothing. I want something that will be reliable/accurate and where the battery lasts a decent amount of time (might go without being able to charge 8 hours at a time). Bonus, this GPS will stay with me for the years to come so I will think of it as an investment, so my budget is fairly loose at the moment. I've drawn some inspiration from these maps:   * http://www.youtube.com/watch?v=iC56P-Ro4PY   * http://www.flickr.com/photos/walkingsf/sets/72157629014750905/   * http://www.flickr.com/photos/walkingsf/sets/72157628738161697/ Anyways, just thought I'd poke around in gis.se because there should be a fair amount of geographers in here with GPS usage experience!


title: PostGIS and invalid SRIDs

text: I'm having difficulty with a PostGIS query. I am trying to query my database (`world_countries`) using a bounding box drawn on a google web map. The objective of this query is to find any countries that lie within a box drawn by the user, and to return them in GeoJSON format. The query is as follows:               SELECT ST_AsGeoJSON(the_geom)      FROM      "world_countries"     WHERE      ((the_geom && 'POLYGON((     -1572563.9181513 7593219.6338481,      -1709539.0728192 6507202.3361235,      -413167.07328311 6477850.5172661,      -271299.948805577549191.905562,      -1572563.9181513 7593219.6338481))'));      and the error message is:               ERROR: Operation on two geometries with different SRIDs     SQL state: XX000      The SRID of this data is 900913 (`srid(the_geom) = 900913`). When I look at this data in ArcGIS and pan around, the values from the box above match up with ArcGIS. Is there a reason why a polygon geometry requires a specific SRID? I assumed that my values were within the range. * * * A quick note about my projection. I originally had a map of countries of the world in geometric form using WGS_84. I then used ArcGIS to re-project this into 'WGS 1984 Web Mercator.prj' based on this thread. As I am using the following version of PostGIS `"POSTGIS="1.3.6" GEOS="3.1.1-CAPI-1.6.0" PROJ="Rel. 4.7.1, 23 September 2009" USE_STATS"` I ran the following `SQL` to include the google projection info:               INSERT into spatial_ref_sys (srid, auth_name, auth_srid, srtext, proj4text) values (900913 ,'EPSG',900913,'GEOGCS["WGS 84", DATUM["World Geodetic System      1984", SPHEROID["WGS 84", 6378137.0, 298.257223563,AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], NIT["degree",0.017453292519943295], AXIS["Longitude", EAST], AXIS["Latitude", NORTH],AUTHORITY["EPSG","4326"]], PROJECTION["Mercator_1SP"],PARAMETER["semi_minor", 6378137.0],      PARAMETER["latitude_of_origin",0.0], PARAMETER["central_meridian", 0.0], PARAMETER["scale_factor",1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0],UNIT["m", 1.0], AXIS["x", EAST], AXIS["y", NORTH],AUTHORITY["EPSG","900913"]] |','+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m      +nadgrids=@null +no_defs');


title: How to identify grid shift file used by osgeo in Python?

text: I wrote a Python script reprojecting points from a source coordinate reference system (SCRS) to a target one (TCRS) both identified by an EPSG code. I would like my python script to retrieve details about the grid shift file it uses when doing reprojections from NAD83 to NAD27. Such information is given by CS2CS in a terminal window when the variable PROJ_DEBUG is set to ON. Can my python script retrieve the same information? Any ideas? Thanks My script is:               from osgeo import osr          #Input coordinates     xin = 436898.0     yin = 5577554.0          # Transformation setup from NAD83 UTM Z15 to NAD27 UTM Z15     srs = osr.SpatialReference()     srs.ImportFromEPSG(26915)     trs = osr.SpatialReference()     trs.ImportFromEPSG(26715)     rsTrans = osr.CoordinateTransformation(srs,trs)          # Doing transformation     (xout, yout, zout) = rsTrans.TransformPoint(xin,yin)


title: can I export service area (network analysis) to shapefile?

text: Is it possible to export service area (network analysis) to shapefile? if so, how can I do it? I am using ArcGIS 10.


title: How to perform a table look-up using Python in ArcGIS 10.0?

text: Using ArcGIS 10.1 the arcpy module has a TableToNumPyArray command, which converts a table into an array. Unfortunately I don't have access to 10.1 on my work PC. Is there a workaround which will accomplish this using ArcGIS 10.0? **Background** I'm iterating through features in a featureclass, and I need to look-up a price depending on an attribute value. I need to find a price for multiple fields per feature, so I can't use a table join. I can't create a relationship class as I only have an ArcView license. A messy workaround could be to load the table into a searchCursor and iterate through each record until I find the desired value. But this would be intensive as I'd need to work through the searchCursor for every feature in the featureclass. Thanks for any advice, Steve


title: Mercator Didn't work on Geoserver

text: Got some problem with Mercator in Geoserver, I tried to display Flight Information Region that has been stored in my Postgre. **The result is more likely drawing using Cartesian projection:** ![enter image description here](http://i.stack.imgur.com/ZYQVl.png) **I tried redrawing the FIR based on the same data using WorldWind and got the result:** ![enter image description here](http://i.stack.imgur.com/NBRjv.png) Here is the configuration of the Coordinate Reference System (EPSG:4326): ![enter image description here](http://i.stack.imgur.com/k6njX.png) And I tried with EPSG:900913 (the result same like the first picture): ![enter image description here](http://i.stack.imgur.com/gsPmF.png) Also I tried to change the Native SRS to 900913 but got error:               <ServiceExceptionReport version="1.3.0" xsi:schemaLocation="http://www.opengis.net/ogc http://127.0.0.1:8080/geoserver/schemas/wms/1.3.0/exceptions_1_3_0.xsd"><ServiceException code="internalError">           Rendering process failed     java.io.IOException     null     ERROR: Operation on two geometries with different SRIDs     </ServiceException></ServiceExceptionReport>      **Is there anything wrong with the configuration?**


title: How can I control what is passed to the OnCreate() method when a command or tool is initialized?

text: In an ArcEngine application can I control what is passed to the OnCreate() method when a command or tool is initialized? (I am using a synchronization object because I have both a map and a pagelayout in my form. Rather than just passing the active control I would prefer to just pass the synchronizer because it contains all the controls.)


title: Sort feature type order in FME KMZ writer

text: I have 4 feature types I am writing to a single kmz file using FME. Currently, these feature types are written as folders in the kmz file and appear in alphabetical order. Is there any way to sort these folders so a particular folder is always first regardless of the first character in the folder name? I know there is a way to sort the features within each folder, but I cannot determine how to sort the feature types themselves within the document. Thanks in advance for your help.


title: How to contour irregular point data for overlay on Google maps?

text: I am having some serious trouble finding a solution to fit my current needs. I am attempting to create contour maps from irregularly spaced point-data on the fly to overlay on a Google map. I have looked at almost everything I can find, but I haven't come across anything that seemed to do this. I cannot generate the contours before use, as there are literally thousands of possible maps, and they change regularly. If anyone knows of anything that might help me accomplish this, it would be extremely helpful. EDIT: I am fairly new to all of this stuff, so there is a good chance I have missed something basic. Please let me know if that is the case.


title: Manage representions by updating rules with original data

text: In case your not familiar with representations.   There is new since 9.2 a cartographic representation available.   It allows for rules based on data fields to determine symbology. Admittedly it is very useful and powerful for single entity manipulation to show it how and where you want.   One problem with it is that there is only one way to get a representation from existing symbology. Using "convert symbology to representation" it creates a new feature class (db only) and creates all the rules and symbology for the representation. These symbols are not the normal esri symbols so there are some issues and restrictions for maintining that.   However at that point you have a (personal, file or sde) gdb with symbology in it.   Whenever you add that data to a new map the symbology is already there and if you have more than one rep you simply click between the reps you have defined. As I said very powerful.   The REALLY big drawback I just found is that when you edit the original value that created the rule the rule is not updated. An easy enough fix if you are symbolizing with only one field.   As luck would have we are not.   Now the question. **Has anyone used a script to update the rule_id field from more than one originating value?**


title: Create a definition query for groups of layers: Is it possible?

text: I have a large number of groups created, with each group having many layers. I have looked for a way to set a definition query for all the layers in each group but haven't found anything. Is it even possible, maybe with model builder or python? I have ArcMap 9.3 and 10.


title: How to visualize shortest path that is calculated using Networkx?

text: I have a graph G (made from a road network shapefile), and a source node S, and a destination node D. I have calculated the length of shortest path using single_source_dijkstra_path_length, but now I need to save the actual path in a shapefile. This will not help, because it returns the graph representation of my network, which looks nothing like my original network. Is there a way to save this path separately as a shapefile?


title: Iterating shapefiles using ModelBuilder or Python?

text: What tool can I use to run many shapefiles to intersect with a shapefile as a output? It looks like that Iterate tools is use for feature classes or datasets and I don't see any for the shapefile. I can use the Batch but I want to use either ModelBuilder or ArcPy because I might need to add more tools later.


title: what is the output coordinate system when I use Near3D tool in ArcScene?

text: I added a 3D rectangular surface (rectangle formed by 3D points) and a 3D point as two different shapefiles and then I tried to use Near3D tool. Then it gives an error "Datum conflict between input and output". I checked my coordinate systems and all the three (Scene, point and polygon) are in same geographic coordinate system which is "WGS_1984". What am I missing?


title: How can I make Leaflet tiles appear in a page?

text: I am having trouble with leaflet working, but not displaying any of the background tiles in the page. This one is using the leaflet-rails gem in a Rails app. ![enter image description here](http://i.stack.imgur.com/Dv1lh.png) This is just in an html sample page: ![enter image description here](http://i.stack.imgur.com/x1sgR.png) I have tried every browser I have and I am not sure where to start --since apart from the backgrounds the map seems to be working fine.


title: How do I use code (preferably c#) to clip all points within a shapefile that are not contained within a polygon?

text: So the task is simple: I want to take a point shapefile as input and remove all of the points that are not contained within a polygon. I'm using DotSpatial to do this, but the catch is that there are millions of points, and it's just... hilariously slow. I'm using some code that looks like this:               for (int i = 0; i < points.ShapeIndices.Count - 1; i--)     {         if (polygons.ShapeIndices[activeRow].Intersects(points.ShapeIndices[i]))         {             allPoints.Add(points.Features[i]);         }     }      This takes minutes to do 80,000 points. There's obviously something janky with their implementation but I don't have time to run it down. So my question is, if you had to get this done and introduce minimal dependencies, how would you do it? I could install PostGIS and load the shapefiles into PostGIS to do the contains operation, but I am loath to introduce a database into this project.


title: Getting the outer boundary of a group of polygons

text: ![enter image description here](http://i.stack.imgur.com/98Tjg.png)I have an issue where I have feature classes of groups of zip codes in different sectors of the US that I have created. I want to be able to define the outer boundary of these groups so I can create a single shape file (one polygon). I was going to use this new polygon to clip the original file. Is there a way to take many polygons and create a file that only shows the outer extent of the group? And how would I go about such an endeavor?


title: Moving graphics in Esri Flex Feature Layer

text: I have an ESRI Flex 2.3 FeatureLayer with some graphics in it. I am able to access the graphic`s geometry and parse the geometry as Esri`s MapPoint class for Flex Api. I am calling the update(x,y,srs) function..but the change doesnt seem to reflect on the fly unless i zoom in or zoom out. I am trying to make an animation in which the point graphics would move to given locations all at once, but the update() function is not effecting on the run, I called featurelayer.refresh() but that did not help at all. Is there a way to invoke feature layer`s redraw and not refreshing it so that it doesnt request data from service again? my code..               //movementfeatureLayer is the layer which has point     var features:ArrayCollection=movementfeatureLayer.graphicProvider as ArrayCollection;     for(var i:int=0;i<features.length;i++)     {     var g:Graphic=features.getItemAt(i) as Graphic;     var geomPoint:MapPoint=g.geometry as MapPoint;     var x:Number=geomPoint.x;     var y:Number=geomPoint.y;     (((movementfeatureLayer.graphicProvider as ArrayCollection).getItemAt(i) as Graphic).geometry as MapPoint).update(x+500,y-500,map.spatialReference);                         }


title: can Geocat Bridge enable me to connect to existing PostGIS database and make edits through ArcGIS?

text: We are planning to set up a Webgis to share our GIS data on the web. Currently, we are on a desktop GIS and using ESRI's File geodatabase format. We are considering ArcGIS Server or using Opengeo suite. We want to retain using our ArcGIS Desktop but we plan to move our data to a Postgres/PostGIS database to enable multi-user viewing and editing when served on the web. To be able to connect to a Postgres/PostGIS database without ArcSDE license, there are a number of solutions available, ST-links, GISquirrel, Geocat Bridge... Geocat Bridge seems appealing cause we can directly publish mxd to Geoserver. But, it's not mentioned on their site whether we can **connect to existing PostGIS data**. There is a mention that we **can load TO PostGIS** but not a thing mentioned on whether we can load data on ArcGIS Desktop **from PostGIS.** Anyone who has tried Geocat Bridge? Is this possible in Geocat Bridge?


title: How to change feature color?

text: I have feature (`polygon`) was gotten by `spatial query` from the map. How to change its color? I used this code but it didn't work               IMxDocument mxDocument = application.Document as IMxDocument;     IActiveView activeView = mxDocument.ActiveView;     IRgbColor rgbColor = new RgbColorClass();     ISimpleFillSymbol fillSymbol = new SimpleFillSymbolClass();     IElement element = new PolygonElementClass();     rgbColor.Red = 255;     rgbColor.Green = 40;     rgbColor.Blue = 50;     fillSymbol.Color = rgbColor;     element.Geometry = parcel.Shape;     IFillShapeElement pFillElement = (IFillShapeElement)element;     pFillElement.Symbol = fillSymbol;     activeView.Refresh();


title: Convertfeaturestographics Tool or Program

text: Problem : With over 19 feature classes in our dataset, we need to convert all these featureclasses to graphic . But to do it programmatically, we searched the resources but could not find a Python script for converting features to graphics. Though below VBA code was available , I hope it wont help us in 10.1 http://edndoc.esri.com/arcobjects/9.0/samples/arcmap/convertfeaturestographics.htm Mightbe we are missing some trick here ?


title: How to cut short a shapefile?

text: I have a large world map in a shapefile from which I need to crop a smaller region. I am planning to use in in java based simulation. How to do that? I am completely new to GIS world so can anyone recommend a decent tutorial on basic manipulations to help me understand what is actually going on?


title: Connect points along bathymetry contours

text: I have a table with start and end point data in each row and need to convert them into lines, but I would like something other than a straight segment.   Ideally, I want to connect the points via a segment that is determined by the bathymetry of the area (so along the same bathymetry as the start and end points). I have both raster bathymetry layers and contour shapefiles.   The only way I can think to do this is create a raster layer with the difference in bathymetry from the mean bathymetry of start and end points and create a least-cost path with large differences being very costly.   I wanted to float the idea out there though to see if there are any automated or different ways to do this. Here is a diagram of what I have/want:   ![enter image description here](http://i.stack.imgur.com/FyjAC.png)   The black lines are bathymetry contours, and the red points are examples of my start and end points. The blue line is the straight line segment that is how ArcMap tools would connect the two points.   What I need is a line connecting them that follows the bathymetry contours, i.e. the path with the least amount of change in bathymetry from the start and end points. After exploring the idea, it is not feasible to simply snap the points to the bathymetry contour shapefile and split the lines, that would change the location of the points too much.


title: Using Alter Field tool of ArcGIS for Desktop on shapefile?

text: According the ArcGIS 10.2.2 help this should work:               inZoneData = sys.argv[1] # This is a shapefile     fieldList = arcpy.ListFields(inZoneData)          for field in fieldList:             if field.name.upper() == 'STATUS':  # rename status                 arcpy.AlterField_management(inZoneData, field, 'status')      but it does not. What I get is invalid parameter. If I try this:               for field in arcpy.ListFields(inZoneData):         if field.name() == 'STATUS':  # rename status             arcpy.AlterField_management(inZoneData, 'STATUS', 'status')      I get invalid unicode error. I understand that error at least. Can anyone shine some light on how AlterField actually works? I've tried several other things but cannot get this tool working. Is it a shapefile bug?


title: What is a .ssf file?

text: What exactly is a .ssf file? I am refering to the .ssf files used with Trimble GPS units. I am aware that it is some sort of post processing correction file, but I do not know more than that.


title: How to get Geodatabase into Oracle database/ Tables

text: Simple question possibly. But how can I import a Geodatabase into a oracle database. I like to put data into tables but just wondering what is the best way or tools to use? Will it be through FME or any other process?


title: GP-Tool: Detect if running in background or foreground ? (ArcObjects)

text: I got a geoprocessing tool which I want to prevent from running in background, due to some weird crashes that happens. Is it possible to programmatically detect if background processing is enabled ? (tried to look into environments with no luck). I know there is a mechanism for this in 10.1, but I am stuck with 10.0 SP 5.


title: How to create points based on the distance and bearing from a survey point?

text: I am plotting witness tree data to calculate pre-settlement tree density. I have a shapefile with points along section lines (square miles) at each half mile. I am making another table that contains tree data - distance, bearing, and diameter - related by survey point to the existing shapefile attribute table. Is there a way to create points on the map for the trees based on the distance and bearing from each survey point? I'm not necessarily asking for a detailed set of instructions (although if anyone wishes to, I'd be grateful) but rather just if it's possible and a couple key words to continue searching how to do it. Thanks for your time, Tim


title: How does one replace characters in a layer name using Python ("_" with " ")?

text: I can rename layers easily enough, but I really need to replace any "_" with a " ". This code just renames layers, but I would like to do a find and replace on each layer name               mxd = arcpy.mapping.MapDocument("current")     layers = arcpy.mapping.ListLayers(mxd)          for lyr in layers:         if lyr.name == "X":             lyr.name = "Y"          arcpy.RefreshTOC()


title: How to attach text from HTML file to objects via WMS GetFeature Geoserver

text: Hi I've been working with a new OpenLayersapplication and I'm using WMSGetFeatureInfo to get information about layers from Geoserver. I'm using polygon layers containing buildings. I can add images with content.ftl for every building. My problem now is that I want text information attached aswell from html- files. Is there a way I can call these html-files to get the information displayed when I click on certain buildings?


title: Why does OpenLayers render these two AGSCache layers differently?

text: I am using OpenLayers to try to display aerial maps for two cities (McHenry, IL and Cupertino, CA). By following the tutorial for OpenLayers I have gotten my McHenry example to render correctly, as shown in this jsfiddle: http://jsfiddle.net/je45Y/ I should be able to use the exact same code and change only the URL and the zoomToExtent() arguments at the end (such that they match the Initial Extent reported on the ESRI info page). However, when I do this, all the map tile requests generate 404 errors as shown in this jsfiddle: http://jsfiddle.net/9y90mvLo/ Eventually, if I zoom in enough, the tiles will display. This is pretty confusing to me. Presumably, there are different values for zoomToExtent() that I can provide to make this work properly right from the start, but I have no idea how I'd find out what these values are. Is their ESRI server misconfigured somewhow? Or does OpenLayers not work properly with this newer version of ESRI Server? It's confusing that the same code works for the first server but not the second. Any assistance would be appreciated.


title: Shifting problem with mapfish

text: I'm using OpenLayers 2.13 with GeoExt and mapfish print module for GeoServer. I've a map composed of a WMS layer and a vector layer. On my web browser everything is fine but when I print this map into a PDF file with mapfish print, layers are shifted and I can't figure out why. Any ideas? How the map is rendered by OpenLayers/GeoExt ![How the map is rendered by OpenLayers/GeoExt](http://i.stack.imgur.com/taa4g.png) How the map is rendered in PDF file generated by Mapfish print ![enter image description here](http://i.stack.imgur.com/Nd6oh.png) Part of my code : Map creation part:               map = new OpenLayers.Map("map", options);     var layer_dpt = new OpenLayers.Layer.WMS("dpt", "linkToGeoserver",         {             layers:'base_layer'         },         {             singleTile: 'true'         }     );     map.addLayer(layer_dpt);     layer_comm = new OpenLayers.Layer.Vector( 'communes',{     strategies: [new OpenLayers.Strategy.Fixed()],     protocol: new OpenLayers.Protocol.HTTP({     url: 'GeoServer WFS URL',     format: new OpenLayers.Format.GeoJSON()     })       });     map.addLayer(layer_comm);      Printing part:               var printProvider = new GeoExt.data.PrintProvider({         method: "POST", // "POST" recommended for production use         capabilities: printCapabilities // from the info.json script in the html     });     printProvider.setDpi(printProvider.dpis.getAt(2));     // Our print page. Tells the PrintProvider about the scale and center of     // our page.     printPage = new GeoExt.data.PrintPage({         printProvider: printProvider,         customParams: {             mapTitle: "",             comment: ""         }     });     mapPanel = new GeoExt.MapPanel({         region: "center",         map: map,         center: [map.getCenter().lon, map.getCenter().lat],         zoom: map.zoom     });     new Ext.Panel({         renderTo: "map",         layout: "border",         width: 1200,         height: 950,         items: [mapPanel],         bbar: ["->", {             text: "Print",             id: "imprimer",             handler: function() {                 // convenient way to fit the print page to the visible map area                 var opt = {"mode":"closest"};                 printPage.fit(mapPanel, opt);                 // print the page, optionally including the legend                 printProvider.print(mapPanel, printPage, false);             }         }]     });      Thank you in advance.


title: overpass-api on opensuse 13.1

text: i want to install overpass.api on opensuse 13.1 what is needed to do that. are there things i have to get via YAST - or do i need to get all via Terminal!? Any and all help will be greatly appreciated greetings


title: Is there a way to add a "current user" label to the print composer?

text: Is there a way to add a current user label to the print composer in QGIS. I thought there might be something similar to the $CURRENT_DATE format available. This would be in a windows environment.


title: Programmatically block Toolbar choosing in ArcMap 10.1

text: With an extension I can block ArcMap's Customize menu using following code line:               ArcMap.Application.LockCustomization("password");      But, it makes no sense since by right-click on Toolbar docking area, user still can add any toolbar s/he wants. Is there a way to block that right click menu too?


title: ArcGIS Online and Sign In options are greyed out in ArcMap 10.2.2

text: I have just installed ArcGIS for Desktop 10.2.2 on my Windows 7 machine. When I got ArcMap opened up I found the "ArcGIS Online ..." and "Sign in..." options were greyed out. Meanwhile in ArcCatalog there was no GIS Server -> ArcGIS Online option. Just wondering if there is any options/settings I should tick on to get them working? Thanks. Cheers,


title: geojson object with different representations and label/popoup conten

text: I have a quite complex geojson in terms of vertexes. It holds all the information which is needed in the webmap. it serves as the input for many layers like this (example for two):               var represenation1 = new L.geoJson(GEOJSON,{                     onEachFeature: popup_function_rep1,                     style: function (feature) {                         return {fillColor: feature.properties.color1,                                 weight: 1,                                 color: feature.properties.color1,                                 opacity: 1,                                 fillOpacity: 0.6};                         }                     });                 representation2 = L.geoJson(GEOJSON,{                     onEachFeature: popup_function_rep2,                     style: function (feature) {                         return {fillColor: feature.properties.color2,                             weight: 1,                             color: feature.properties.color2,                                 opacity: 1,                                 fillOpacity: 0.6};                         }                     });      Doing so for two layers might be fine but mine has eight and crashes the browser when it comes to the style I know I could do sth like:               representation2 = representation1.setStyle(     function (feature) {                         return {fillColor: feature.properties.color2,                                 weight: 1,                                 color: feature.properties.color2,                                 opacity: 1,                                 fillOpacity: 0.6};                     });      But I am struggling with the onEachFeature section. Both layers will be added to the L.control:               L.control.layers({"Rep1": representation1, "rep2": representation2}).addTo(map);      Are there any possibilities to keep ressources small and just alter the style and function onEachFeature according to the current active layer?


title: polygon to raster - not all polygons rasterized

text: I have a building footprint shapefile (polygons) that I want to convert to a raster. I use the Polygon to Raster tool, and specify a cell size by importing from another raster (which is 24m x 30m). The results appear to show that only a couple of the polygons have been rasterized. Why is it that some polygons are and some aren't? I've attached a picture below showing my results. Pink = building footprint; red = raster pixels after using tool Polygon to Raster. I would like to have a raster file that covers ALL of the building footprints (pink). ![enter image description here](http://i.stack.imgur.com/37x71.png)


title: Java, GIS knowledge...What next?

text: I know programming in Java and I am familiar with GIS concepts as well. I do GIS programming using javascript and PHP. I want to now do it using Java. How should I go about it? PLease let me know where I should begin with? Links to tutorials or simple Java based GIS project works will really be appreciated! Regards, Shreerang Patwardhan.


title: Where to create ArcGIS 10 Map Templates?

text: Sorry, I could not find the location: C:\Users\AppData\Roaming\ESRI\Desktop10.0\ArcMap\Templates Do I create it or there is a specific funtionality for the location? I would like to create different templates for the areas we mostly use the maps.


title: Drawing arcs in AutoCAD

text: My question is extremely basic but has been causing me a headache. I wish to draw a shape that looks like the following shape: ![A mock up](http://i58.tinypic.com/5cf7nr.png) I have very exact widths of each of the arcs in mind, and a constant spacing between the arcs. How can I draw this in AutoCAD? The basic arc tools have been a headache. Cheers for any tips.


title: Problems using v.split in SEXTANTE toolbox to split lines

text: I got this error： ![enter image description here](http://i.stack.imgur.com/mPxMh.png) Or don't click "open output file after running algorithm", but there is no new shapefile created on my desktop. How to solve this problem and split lines at a specified distance？ PS. I use QGIS 1.80 under WIN7-64 bits


title: Measuring the distance between lines and points in QGIS

text: Using QGIS - I have a vector of fault lines. I have point data in the form of a csv which I have also loaded into QGIS. I need to know how far each point is from the nearest fault line then enter this distance in a new column in the csv file. Any ideas?


title: Autocomplete against geocode.arcgis.com?

text: I have built a tool that allows my users to perform geocode requests against the World Geocoding Service provided from Esri. I would like to extend my tool to have an autocomplete feature, similar to how this can be achieved with the Autocomplete class from the Google Geocode API. Although there is not a similar class available at the Esri service end, I was thinking that you could still hit the geocode service with a GET request, after listening for a slight delay after each keystroke from the end-user in something like a combo-box, and just return the top 5 results based on score. I am looking for an answer that provides some code (can be pseudocode, or any language) on how to accomplish this, or if this workflow is even achievable. Something similar to this, but against the Esri Geocode service as opposed to map services.


title: Determine if point is within an irregular polygon using Python

text: I would like to create a mask of a grid where points are evaluated as either `1`'s or `0`'s depending on if they lay outside of an irregular (i.e. not rectangular) polygon. I have a set of high-latitude coordinates in lat/lon that define the boundary of an area I would like to make a polygon. From this polygon I would then like to evaluate another set of lat/lon points as to whether they lay within the polygon. Is there a library best suited for creating such a polygon (given a list of lats/lons that define it)? `shapely` appears to be geared towards 2D/projected points. What would be the best method of evaluating other points as to whether they lie within or out of this polygon? Thanks!


title: Comparison of vector file formats

text: Does anyone know of a chart or table that would explain the benefits / shortcomings of the different vector file formats, eg. .shp, .gml, .kml, .gdb. It seems that ESRI shape is still quite prolific, but since we are constantly "trading" data with other places, pieces of the shape data tend to be forgotten (.shx, etc.). We are thinking of switching to .kmz files but what would we lose by doing so? I don't believe all of the attribute data comes over with the .kml/kmz files. Is this comparison possible, or am I trying to compare apples and oranges? Thanks.


title: Qgis 2.0: Several symbols in one line behind a label?

text: Is that possible in qgis 2.0: ![enter image description here](http://i.stack.imgur.com/iAyE4.png) If a city have several POIs which I can not draw at their origin location because of the scale of the map or because there are just too much symbols, can I draw them in one line behind, or below a label?


title: How to make WFS-T layer grid panel with Geoext

text: I am using geoext,openlayers and geoserver. I am using WFSCapabilitiesStore to make the grid panel. I want to this grid panel WFS-T layer so that on each layer I can draw and edit feature after adding these layer to tree panel. My store object for for grid panel is this:               store = new GeoExt.data.WFSCapabilitiesStore({              url: "http://mindcrewgis.com/geoserver/wfs?service=wfs&version=1.0.0&request=GetCapabilities",              protocolOptions: {version: "1.0.0"},               layerOptions: {                              visibility: false,                               displayInLayerSwitcher: false,                                strategies: [new OpenLayers.Strategy.BBOX({ratio: 1}),saveStrategy ],                                projection: new OpenLayers.Projection("EPSG:4326"),                               geometryName: "geom",                               srs: 'EPSG:4326',               }              });          // load the store with records derived from the doc at the above url            store.load();      The above store object is displaying properly the grid as well as all layers, but after adding the layer to the tree panel and I draw feature on layer. The feature is not saved. Can anybody help me to that? how can i make a grid panel and layer to editable feaure? Thanks in advance.


title: Double version of PostGIS installed on my server

text: For some reason I messed up an installation of Postgis on my PostgreSQL 9.2.8 server. I had version 2.0.3 installed and when I tried to upgrade to 2.1 it did so. The problem comes when I dump the database to restore it on another machine. In the dump there is the create extension statement for version 2.1 of postgis but I also get the create statements for all the postgis associated functions. If I manually remove all of the create function statements from the dump and import the dump in another DB it will work perfectly. What I need to know is how to solve this problem on my current machine. I have posted on Database Administrators but without success. This is the link to the question. What I have tried as per suggestions is the following: `create extension postgis from unpackaged` but this will throw an error stating that postgis already exists. `alter extension postgis update to '2.1.0'` will complete but when i dump the database the create function commands relative to version 2.0.3 are still there. Dumps of the `select * from pg_available_extensions` are on the DB stack exchange question and they clearly show that default_version is 2.0.3 but installed version is 2.1.0 The database is working properly so for it's just a matter of cleaning up this mess that I have created. Suggestions are highly appreciated!


title: Python.dll could not be loaded with PyScripter

text: I'm trying to configure **PyScripter** with QGIS. I've installed it using OSGeo4W and I'm following instructions to setup a `pyscripter.bat` file by copying content of `qgis.bat` file. I have the following code:               @echo off     call "%~dp0\o4w_env.bat"     call "%OSGEO4W_ROOT%"\apps\grass\grass-6.4.3\etc\env.bat     @echo off     path %OSGEO4W_ROOT%\apps\qgis\bin;%OSGEO4W_ROOT%\apps\grass\grass-6.4.3\lib;%PATH%     set QGIS_PREFIX_PATH=%OSGEO4W_ROOT:\=/%/apps/qgis     set GDAL_FILENAME_IS_UTF8=YES     rem Set VSI cache to be used as buffer, see #6448     set VSI_CACHE=TRUE     set VSI_CACHE_SIZE=1000000     Start "PyScripter" /B "C:\Program Files\PyScripter\PyScripter.exe" --python27 --pythondllpath=C:\OSGeo4W\bin      However when I run this file I get the following error and PyScripter exits after that. ![enter image description here](http://i.stack.imgur.com/qdjE1.png) Where am I wrong?


title: Count raster value in a bounding box (spatial extent)

text: I have simple question: given a bounding box, how can I count the number of specific values (say, I am interested in the number of `value == 1` ) in a raster in 1) arcpy, or 2) raster package in R?


title: How to Merge/Dissolve Polygons By Common Attribute Field in QGIS?

text: I have joined corresponding Local Enterprise Partnerships table data into a polygon shape file containing UK Local Authorities. ![enter image description here](http://i.stack.imgur.com/lO5ou.png) How do I Dissolve or Merge the the Local Authorities polygons by LEP_NAME as shown in the screenshot? Ideally, I would like to retain the Local Authority attributes with the resultant merged polygons. I have zero scripting knowledge so am looking for a tool based solution. Any help or pointers much appreciated. ![enter image description here](http://i.stack.imgur.com/knj5b.png) Additional Note: When attempting to use the dissolve tool the LEP_NAME field is not available as shown above. The LEP attribute data was created in Excel, saved as CSV format and joined to the Local Authority polygon shape file. There are duplicate IDENTIFI0 fields as shown in the other screenshot - could this be the root of the issue?


title: Openlayers with allOverlays = true. How to create a dummy (base) layer?

text: I have an OpenLayers map with allOverlays set to true. Im adding 3 layers to the map. First, my own rendered version of OSM data (but only for levels 8-16), secondly an orthophoto layer (for levels 13-19) and 3rd layer is a vector, defined on all levels. Order is important, because of overlaying. Also, note that resolutions array i pass into the map object has no effect. If i initialize the map in this way:                   map = new OpenLayers.Map('map_canvas',          {         projection: p3857,         displayProjection: p4326,              resolutions: [156543.0339, 78271.51695, 39135.75848, 19567.87924, 9783.939619, 4891.969809, 2445.984905, 1222.992452, 611.4962262, 305.7481131, 152.8740565, 76.43702827, 38.21851414, 19.10925707, 9.554628534, 4.777314267, 2.388657133, 1.194328567, 0.597164283, 0.298582142, 0.149291071, 0.074645535],         units:'m',         allOverlays: true,         controls: [new OpenLayers.Control.Navigation(),                     new OpenLayers.Control.LayerSwitcher({'ascending':false})]         });      ...and add the layers in the order i specified, i have a problem, because of allOverlays: true it assumes the OSM layer as the fictional base layer (because its the first one added to the map) from which the map calculates its properties. These are incorrect because that layer isnt defined for all resolutions. I then simply copied the initialization code for my vector layer in front of XYZ layer (with OSM data) initialization and then it works ok, but the vector layer is rendered twice of course. I want to get rid of this dummy first layer which is needed to correctly initialize the map object. How can i do that? I suppose i could use some world coastline shapefile and render it, but if its possible i'd rather just remove it somehow.


title: What open source software should I use to geocode the addresses?

text: Need to geocode the student addresses using the regular 5-digit zip codes. Geocoder will give us an accurate point on the map, so in a Latitude and Longitude environment that should suffice for our specific data needs with respect to answering the question of what students belong to what congressional district. ArcGIS for Desktop software will geocode the student addresses and get points on the map. Then we need to create two other layers one for the representative districts and the other for the senate districts. Then we have to use the Spatial Join tool in the toolbox to combine the tables and give us one attribute table that includes all of our geocode located points and what congressional district they belong to. So ArcGIS for desktop will be able to solve my puzzle, but is their any open GIS tool available which will solve this issue?


title: How do I use a layer of Bing / Google a "WebApplication" in "ArcGIS Server Manager"?

text: I'm using "ArcGIS Server." I want to publish a "WebApplication". I created a map of my city. I want to have as the background world map with satellite images of "BingMaps" or "GoogleMaps" or "ArcGIS World Map free". See the image below. I tried it. But it bring to me this error: "Unable to add layer. Request failed with HTTP status 417: Expectation Failed." https://sites.google.com/site/christiancrca/erro%20BING%20maps%20Arcgis%20Christian%20Rogerio.jpg Any idea about how to use as background a service of image satellite maps with my map? ![Error Status HTTP 417](http://i.stack.imgur.com/gXzup.jpg)


title: Error in shapefile?

text: I have been provided a set of shapefiles (a couple of thousand). I don't own any GIS software but I would like to convert the shapes to KML and Google's encoded polyline. I have put the shapefile and associated files online here in a zip file http://www.goo.gl/T5yu6r However, there appears to be a problem with the original .shp file and I wonder if anyone can confirm it or provide a solution. When I view it using shp2kml (http://www.zonums.com/shp2kml.html), I can see a GPS trace, but there are no Lat or Lon values (which I do see if I load up a file known to work). Consequently, the resultant KML file does not work on Google Maps (and confirmed with KML verifiers to contain errors) I have tried to both view and convert it at Geocommons also, but it does not work. So, I suspect that the original .shp file is missing some data. How can I resolve this problem?


title: Modis Reprojection Tool mirror or download alternative?

text: I need to convert and reproject a complete archive of MODIS *.hdf data. NASA's Modis Reprojection Tool is literally designed for this, and I even have an old script from the last time I performed this task on this very archive. However, I can't download the MRT right now because the US Government Shutdown has shut down the download site: > https://lpdaac.usgs.gov/tools/modis_reprojection_tool > > The requested service is temporarily unavailable. It is either overloaded or > under maintenance. Please try later. Does anyone know an alternate repository where I can download the Linux64 (or 32..) source files or .deb for the MRT?


title: Create buffer around point in LeafletJS?

text: I am using `LeafletJS`. I want to create buffer around point. after Selecting a point buffer radius will be specifed by user in feet. After searching of awhile I found this link of JSTS. does JSTS buffer works with leafletJS ?? I am not able to find any sample for this. Any Other solution to archive same will be appreciated. **Update 1** After look into the code of JSTS and searching for a while, I can to a conclusion there is no conversion between JSTS geometry and Leaflet JS geometry. so it is not possible to use it without implementing conversion.


title: Handling Mapquest different zoom levels for lower 48 USA and rest of the world

text: The Mapquest Open data site states: > Please also note that global coverage is provided at zoom levels 0-11. Zoom > Levels 12+ are provided only in the United States (lower 48) I want to be able to use the Mapquest 12+ Zoom levels in the USA lower 48 and provide my own tiles for outside the USA lower 48. I am using Open Layers and with Geoserver for my own tilesets for outside the USA What is the best way to determine if my viewport or map center is contained in or out of the lower 48 so I can switch layers?


title: qgis plugin: "WindowsError: [Error 6] The handle is invalid"

text: OS: Windows 7 64bit GIS: QGIS 1.8 (qgis.org/downloads/QGIS-OSGeo4W-1.8.0-2-Setup.exe) path variable:               GISBASE=C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2     LD_LIBRARY_PATH=C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\lib          PATH=C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\etc;     C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\etc\python;     C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\bin;     C:\Program Files (x86)\Quantum GIS Lisboa\apps\msys\bin;     C:\Program Files (x86)\Quantum GIS Lisboa\apps\Python27;          PYTHONLIB=C:\Program Files (x86)\Quantum GIS Lisboa\apps\Python27     PYTHONPATH=C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\etc\python     GRASS_SH=C:\Program Files (x86)\Quantum GIS Lisboa\apps\msys\bin\sh.exe      * * * * * * I have created a plugin for QGIS, with the QGIS plugin builder. Now i would like to import grass.script into this script               # Import the PyQt and QGIS libraries     from PyQt4.QtCore import *     from PyQt4.QtGui import *     from qgis.core import *          # Initialize Qt resources from file resources.py     import resources_rc          # Import the code for the dialog     from testdialog import testDialog          # Import grass.script     import os, sys     import grass.script as grass      If i run this Script then i get the following error: WindowsError: [Error 6] The handle is invalid The whole Error massage is:               Couldn't load plugin test due an error when calling its classFactory() method          Traceback (most recent call last):       File "C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python\qgis\utils.py", line 164, in startPlugin         plugins[packageName] = package.classFactory(iface)       File "C:/Users/admin/.qgis//python/plugins\test\__init__.py", line 52, in classFactory         from test import test       File "C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python\qgis\utils.py", line 309, in _import         mod = _builtin_import(name, globals, locals, fromlist, level)       File "C:/Users/admin/.qgis//python/plugins\test\test.py", line 31, in          import grass.script as grass       File "C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python\qgis\utils.py", line 309, in _import         mod = _builtin_import(name, globals, locals, fromlist, level)       File "C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\etc\python\grass\script\__init__.py", line 1, in          from core  import *       File "C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python\qgis\utils.py", line 309, in _import         mod = _builtin_import(name, globals, locals, fromlist, level)       File "C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\etc\python\grass\script\core.py", line 1201, in          debug_level = int(gisenv().get('DEBUG', 0))       File "C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\etc\python\grass\script\core.py", line 558, in gisenv         s = read_command("g.gisenv", flags='n')       File "C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\etc\python\grass\script\core.py", line 238, in read_command         ps = pipe_command(*args, **kwargs)       File "C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\etc\python\grass\script\core.py", line 215, in pipe_command         return start_command(*args, **kwargs)       File "C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\etc\python\grass\script\core.py", line 177, in start_command         return Popen(args, **popts)       File "C:\Program Files (x86)\Quantum GIS Lisboa\apps\grass\grass-6.4.3RC2\etc\python\grass\script\core.py", line 56, in __init__         startupinfo, creationflags)       File "C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\subprocess.py", line 672, in __init__         errread, errwrite) = self._get_handles(stdin, stdout, stderr)       File "C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\subprocess.py", line 784, in _get_handles         p2cread = self._make_inheritable(p2cread)       File "C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\subprocess.py", line 823, in _make_inheritable         _subprocess.DUPLICATE_SAME_ACCESS)     WindowsError: [Error 6] The handle is invalid          Python version:     2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)]               QGIS version:     1.8.0-Lisboa Lisboa, f440b60          Python path: ['C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python', 'C:/Users/admin/.qgis//python', 'C:/Users/admin/.qgis//python/plugins', 'C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python/plugins', 'C:\\Program Files (x86)\\Quantum GIS Lisboa\\apps\\grass\\grass-6.4.3RC2\\etc\\python', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\bin\\python27.zip', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\DLLs', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\plat-win', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\lib-tk', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\qgis\\bin', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages\\PIL', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages\\win32', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages\\win32\\lib', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages\\Pythonwin', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode']      Has someone an idea ? regards Marcus


title: How to get the pixels from a geotiff file in gdal python for a given bbox?

text: I have a GeoTIFF file in WSG 84 Projection reference (it's a GeoTIFF of SRTM). I can use python & gdal to read the whole lot of data into an array (with `ReadAsArray()`), but I want to only have a subset of the data, everything inside one bounding box (defined by a top latitude, bottom latitude, right longitude, left longitude). But I'm very new to gdal, and I can't figure out how to do it. Here is some python code to confirm it is in WSG 84, etc.               >>> from osgeo import gdal     >>> gg = gdal.Open('srtm_34_02/srtm_34_02.tif')     >>> gg.GetGeoTransform()     (-15.000416327675339, 0.0008333333333333334, 0.0, 55.00041712672993, 0.0, -0.0008333333333333334)     >>> gg.GetProjectionRef()     'GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]]'      It doesn't matter if the resultant array of data includes points that are just outside the bbox, due to rounding. Just so long as it doesn't exclude any points that are in the bbox, and it's as close enough as is practical.


title: Exporting a QGIS layout with Addition on layers

text: I've been playing around with an approach in QGIS that uses Addition between the layers I'm mapping to build a tri-variate composite colour map - an example is here in the QGIS Flickr pool: http://www.flickr.com/photos/97267218@N06/12626847074/ To achieve an output I thought I had to set the background of the image to black - in order to draw out the colours. However I've been asked (demanded!!!) to set the background to white. Thinking about it a little, one approach that would seem to work would be to set the background of the Composer to Alpha 0. This is a screen grab of an map produced in a similar manner. In the composer I set the background with an Alpha of 0 - meaning that it doesn't influence the layers. Great! ![Screen Grab of composer.](http://i.stack.imgur.com/qU6kC.jpg) However when I export the final image in any format (including Pdf and SVG) I only get a white page! There has been a similar question asked about this in ArcGIS (Graphic with transparency in ArcGIS 10.1 layout does not export properly in PDF Export) - and the suggestion there was to change the graphic to a feature and insert that .... something I _don't think is possible with QGIS_ ... So - has anyone either any suggestions how this might be overcome? How can I export an image which is using Addition and Alpha _without_ needing a black background? More fundamentally ... can anyone explain (in relative layman's terms!) why the export isn't working in the first place?


title: Min. weighted average distance from multiple points

text: I am a complete beginner in mapping-related topics. Still, I´ve been wondering if anyone know this problem. It might be well known, I just don´t know its name. Initial situation: nx : x geographical locations Wx: weight per location Dx: distance between the location nx and the optimal location O Problem: Where is the geographical position O, where the Σ (Dx*Wx) is at a minimum? Numerical example: n1: Rheinfelden; Wx: 20% n2: Zurich; Wx: 40% n3: Berne; Wx: 30% n4: Lucerne; Wx: 10% Result: somewhere around Aarau in Switzerland Does anyone know the common name for this spatial problem? I have found a tool that can do this on a straight line ("as the bird flies", http://www.geomidpoint.com/) but I would like to be able to calculate this on a street routing level. Is anyone aware of such a tool?


title: Change contour line stroke weight based on elevation?

text: I have shapefiles of contour lines that I have imported in Qgis. At the moment the contour lines are all the same stroke weight. I'm a beginner to this program, but wanted to know if its possible to change the stroke weights of the contour lines based on the elevation?


title: Why might Spatial Join appear to produce incorrect values?

text: I'm trying to do a spatial join of a point file (created from a raster) to a polygon (target), where I am joining the GRID_CODE value from the point to the polygon. I've tried a variety of permutations, one-to-many, many-to-one, changing the order of the input files, putting in a mode and mean condition for field_map, checking the field properties (both Double), trying a multi- part to single-part conversion (data management) for both inputs, but the resulting joined field is always populated with incorrect data (as verified by looking at the point values). For example, here is a screenshot of the correct values from the point file, highlighting the 2 (green)![enter image description here](http://i.stack.imgur.com/pJOoR.jpg) Here is a screenshot of the spatial join output highlighting that same 2 value, which I would expect to be in the same location.![enter image description here](http://i.stack.imgur.com/4mPqJ.jpg) What is additionally bizarre, and might be a clue as to what is going on, is that the same field (GRID_CODE) has a different number of values (4 in the original, and 5 in the output![enter image description here](http://i.stack.imgur.com/5IeaS.jpg) Any suggestions would be much appreciated!


title: Get list of all tables in CartoDB account with SQL statement?

text: How can I get a list of all available tables in my own authenticated CartoDB account with a SQL statement? This returns `null`:               SELECT table_name as t FROM information_schema.tables WHERE table_schema = 'public';      Alternatively, could it be done with CartoDB.js?


title: What is the easiest one click install method for PostGIS and it's entire dependencies?

text: For someone new to database servers and their dependencies, what is the easiest assured way of installing PostGIS on a machine with no other software beyond what comes with its operating system? I'd like to know for Ubuntu, Windows and if possible a common Linux method. I'm understanding that you need to install PostGreSQL first, but before that perhaps you need to install Apache Server, and I'm not sure if that has dependencies. Perhaps there a single binary that loads EVERYTHING for Windows, and maybe there a single program from the Ubuntu GIS repositories that installs everything. If not, what's the install order?


title: DotSpatial Plugins for Measure,attributes,query

text: I want the end product user to open the attributes,measure task and query the shapefile data using the DotSpatial library in c# . Thank you in advance.I have Dll with me but when i create the folder in solution and add dll it doesn't show any change in the output.Thank you in advance.


title: Symbols in map server

text: I was just trying to render a point layer with classes using a symbol circle.Below is the part of the code.               DATA 'E:\SBI\GIS\Growth\Growth11-12.shp'     METADATA       'ows_title' 'Growth11-12'     END     STATUS on     TRANSPARENCY 100     PROJECTION     'proj=longlat'     'datum=WGS84'     'no_defs'     END     CLASSITEM 'GROWTH SCA'     CLASS       NAME "GROWTH SCA = 0_POOR"        EXPRESSION "0_POOR"         STYLE          SYMBOL "circle"           SIZE 5.25           OUTLINECOLOR 0 0 0          COLOR 255 0 0        END     END     CLASS       NAME "GROWTH SCA = 1_BAD"        EXPRESSION "1_BAD"         STYLE          SYMBOL "circle"           SIZE 5.25           OUTLINECOLOR 0 0 0          COLOR 212 42 0        END     END      The other polygon layer renders fine. But this particular point layer doesnt render and throws the following error message: msDrawMap(): Image handling error. Failed to draw layer named 'Growth11-12'. msLoadMSRasterBufferFromFile(): General error message. unable to open file /ms4w/apache/htdocs/xyz/map/circle for reading. Could anyone suggest what am I missing here. Thank you.


title: Create a single connected polygon from multiple polygons

text: I am trying to create a single continuous (neighborhood) polygon from many tax block parcels. The Multipart to Singlepart GP tool does not accomplish the task since I need the polygon to be complete; i.e., without the spatial gaps between the parcels (see image). Any ideas of a quick and dirty way to accomplish this? ![Parcels](http://i.stack.imgur.com/3NABd.jpg)


title: Apply new style to selected features in QGis

text: Let's assume that I have selected some features in QGis. Now all the features in that vector layer already has some style applied to it. Now, I want to apply some new style to the selected features. Is this possible in QGis ?


title: Find longitude for point between two coordinates with pyproj?

text: I've written a function to get the longitude of a point of known latitude between two known points (might not be working correctly yet), but I feel like there must be a better way of using pyproj to do this. There isn't a built in method for doing this, but doing this binary search method I'm using seems possibly unnecessary as opposed to just solving for the other half of the coordinate directly. Any suggestions? Ideally I'd like to generalize this to solve for latitude when the longitude is known, too. My function thing (maybe not working, will update once I'm sure):               import pyproj     import random          def unknownlon(lon1,lat1,lon2,lat2,lat_known):         ''' Binary search by distance to find unknown longitude between points                  latitude and azimuth are known'''              # Get positive for north movement and negative for south         vert_change = (lat2 - lat1) / abs(lat2 - lat1)         horz_change = (lon2 - lon1) / abs(lon2 - lon1)              # Calculate angle and distance between known points         fwd_azimuth,back_azimuth,dist = g.inv(lon1,lat1,lon2,lon2)              # Random distance between known points         dist_result = random.uniform(0.0, dist)              # Loop while difference between desired lat and result is high         lat_diff = 1.0         while lat_diff > 0.00000001:             lon_result,lat_result,back_azimuth = g.fwd(lon1,lat1,                                                        fwd_azimuth,dist_result)                  # Calc difference between desired lat and result             lat_diff = abs(lat_known - lat_result)                  # Determine whether to use a distance higher or lower than last             if (lat_result * vert_change) > (lat_known * vert_change):                 dist_result = random.uniform(0.0, dist_result)             else:                 dist_result = random.uniform(0.0, dist - dist_result)                 lon1 = lon_result                 lat1 = lat_result              return lon_result, dist_result, lat_diff          if __name__ == '__main__':              lat_known = 65.0              lat1 = 68.0         lon1 = 3.0         lat2 = 63.0         lon2 = -12.0              g = pyproj.Geod(ellps='WGS84') # Use WGS84 ellipsoid               # Calculate an unknown longitude         lon_result, dist_result, lat_diff = unknownlon(lon1,lat1,                                                        lon2,lat2,                                                        lat_known)


title: STIntersects and STContains fails on parameter

text: I am new to GIS and using GIS in SQLServer, so apologies if this is a rudimentary question, but I cannot find the answer so far. I am trying to use STIntersects() or STContains to determine what state a given lat/long combination is in. My table named SpatialTable has polygons for all the states and my example below is coded for a point inside my polygon named AL02. I need to understand why this code (cloned from example code on the Microsoft TechNet pages) works fine:               DECLARE @g geography;     DECLARE @h geography;     select @g=GeogCol2 from SpatialTable where name = 'AL02';     SET @h = geography::Parse('POINT(-88.09 31.24)');     SELECT @g.STIntersects(@h);      and this code does not:               -- this will not work as it sees my parameter as nvarchar(max)     --error message is: Cannot call methods on nvarchar(max).     SELECT a.name     FROM SpatialTable a     WHERE a.GeogCol2.STIntersects(geography::Parse('POINT(-88.09 31.24)')) = 1     and a.name = 'AL02';      I also tried STGeomFromText instead of Parse, but still get the nvarchar error.               SELECT a.name     FROM SpatialTable a     WHERE a.GeogCol2.STIntersects(geography::STGeomFromText('POINT(-88.09 31.24)',4326)) = 1     and a.name = 'AL02';      Even if I pass @h into the SQL query with @h declared as geography, the error shows that the database is treating the parameter as nvarchar. How can I dynamically load my lat and long without passing nvarchar? Any help would be greatly appreciated. Thanks.


title: Changing georef backgound colour

text: Working with a greycale scanned image: when I georef that image it rotates (expected) the blank corners are then filled with black. When the new image overlays I need to make those corners transparent. Obviously I cannot select black as a transparent Colour as the image is greyscale. Is it possible to persuade the GeoRef app to create those blank corners in, say, red - or some other colour so that I can make it transparent? Allan