title: VSI File API support in ogr2ogr

text: It seems that support for streaming input is wishy-washy in ogr2ogr. For example, why can I do this:               curl "http://data.usgin.org/arizona/wfs?service=WFS&version=1.1.0&     request=GetFeature&typeName=azgs:activefaults&maxFeatures=10" |      ogr2ogr -f "KML" /vsistdout/ /vsistdin/      ... but I can't do this:               curl "http://data.usgin.org/arizona/wfs?service=WFS&version=1.1.0&     request=GetFeature&typeName=azgs:activefaults&maxFeatures=10&     outputformat=application/json" | ogr2ogr -f "KML" /vsistdout/ /vsistdin/      The first request uses the "WFS" driver for the source layer, and the second request uses the "GeoJSON" driver for the source layer. The first one works fine, but the second gives me:               ERROR 1: GeoJSON parsing error: unexpected end of data (at offset 6000)     ERROR 4: Failed to read GeoJSON data     FAILURE:     Unable to open datasource `/vsistdin/' with the following drivers.


title: What is the correct coordinate system for google geocoded text file layers?

text: When importing a text file containing lat/longs geocoded from something like GPSVisualizer.com into QGIS, what is the correct coordinate system to select? Most web services (GpsVisualizer included) geocode in Web Mercator, but QGIS has no Web Mercator coordinate option.


title: How many active QGIS users are there in total?

text: I'm the CEO of MangoMap which is a online GIS data management and web map publication platform, we're big QGIS fans and make extensive use of it. We would really love to integrate MangoMap more tightly with QGIS by making the transfer of data and edits completely seamless between the desktop projects in QGIS and MangoMap data/maps in the cloud. We already allow QGIS to get layers from MangoMap using WMS but want to make it much easier to send data and edits from QGIS to MangoMap. Before we get started though I would really like to know if there are any up to date numbers out there on how many active QGIS users there are at this time. I've seen an article on Spatial Galaxy from 2011 that estimated 100k, but it's my feeling that QGIS has matured and gained a lot of traction since then so the figure is unlikely still valid.


title: Merge images first, or georeference first?

text: I'm wanting to georeference and vectorise some features from scans of the OS 1896 Town Plans that I have purchased from the NLS (http://maps.nls.uk/os/london-1890s/info.html). The scans provided are not georeferenced and when I try to crop the borders they are not quite square so will need rotating first. The problem is the features lie at the intersection of 4 separate sheets and I'm uncertain of the best / most accurate approach to take to merge them: 1) Treat each image separately: Rotate and crop each image in Gimp, georeference in QGIS, then "nudge" the individual GCP points to try get them to line up properly, then merge with gdal_merge? 2) Merge first: Rotate, crop and merge the images together in Gimp, then georeference the whole image in QGIS 3) Some magic automated solution I'm not aware of :-) Thanks!


title: Intersecting Polygons

text: I have 10 shape files containing polygon features. I would like to produce a single polygon layer containing portions of the input polygons where no fewer than two of the aforementioned feature classes intersect. Any two (or more) polygons which intersect will have those portions included in the output feature class. Additionally, each of the input feature classes contains a unique field name representing itself. I would like for the output feature class to contain EACH of the unique field names from the input feature classes, as well as the value contained in the attribute, or 0 if the polygon did not intersect. So, if a polygon from feature class A and feature class B intersect, that intersecting portion will exist in the output feature class. 8 of the 10 attributes in the output will contain the value "0", and the attributes for the fields representing feature classes A and B will contain the values from those input features. I hope that makes sense. I would like to accomplish this task using ArcGIS. Thank you very much.


title: label for circle marker in leaflet

text: I am able to add label to circlemarker like this               L.circleMarker(points[i],{title: 'unselected'}).bindLabel('Destination').addTo(map);      This adds label which appears on mouse hover on circle marker. But I want to add static label which will appear regardless of mouse is on that circle marker or not. I am referring this demo http://leaflet.github.com/Leaflet.label/ for adding static label to circle marker but some how I am not able to do it. It is working fine with markers but with circle Markers static label is not working. Also is there any other method to add label on circle marker ?


title: Python in ArcGIS, CalculateField_management, updating with single value

text: I have a float field which I want to simply update all values with 1. I can easily do this within the field calculator but I have started using python today, what is the syntax to simply update a field with the value 1. I'm stuck on the expression bit of the script: > > > arcpy.CalculateField_management("ROCHDALE_TEST_NET","FRMEAS",??????


title: How to get route between Two lat/long using ArcGis Map in Android?

text: What I'm trying to do is have two edit boxes for source and destination. Based on those two locations I need to draw the route on an ArcGIS map in Android. Any ideas or suggestions would be greatly appreciated!


title: Ext JS 4.2 and OpenLayers without GeoExt

text: Is it possible to use only Ext JS 4.2 and OpenLayers 2.12 without GeoExt or anything else? (to make a basic map application) I couldn't find any example.


title: Unable to reflect JSON in grid

text: I am trying to fetch the data from JSON and display it onto the grid. The `store` has all the JSON data, but I it is not visible in the Grid. The grid comes out completely blank. This is the application, `http://128.196.142.12/geo/test/test_new.html`, select the `Bioclim Data` tool from the toolbar, on the top right and then click on the map. Here is the entire code,               action = new GeoExt.Action({             text: "BioClim Data",             control: new OpenLayers.Control.Click({                  //control.events.on({"click": function(evt) {                      trigger: function(evt) {                           var lonLat = evt.xy;                           var proj = new OpenLayers.Projection("EPSG:4326");                           var nlonLat = new OpenLayers.LonLat(panel.map.getLonLatFromViewPortPx(lonLat).lon, panel.map.getLonLatFromViewPortPx(lonLat).lat);                           nlonLat.transform(panel.map.getProjectionObject(), panel.map.displayProjection);                           if (popup) popup.close();                           var  request = Ext.Ajax.request({                           url:"/cgi-bin/proxy.cgi?url=http://128.196.142.12:9000/info?lat="+ nlonLat.lat +"&lon="+ nlonLat.lon,                           disableCaching: false,                           success: function(response) {                                   var obj = Ext.decode(response.responseText);                                   var jsonData = Ext.util.JSON.decode(response.responseText);                                   //console.dir(obj);                                   store = new GeoExt.data.FeatureStore({                                       layer: features,                                       fields: jsonData,                                       autoLoad: true                                   });                                   addToPopup(lonLat, response.responseText);                                        var grid = new Ext.grid.GridPanel({                         title: "BioClim Table",                         store: store,                         width: 600,                         height: 250,                         stripeRows: true,                         columns: [{                             header: "Annual Mean Temperature",                             width: 200,                             dataIndex: "mat"                         }, {                             header: "Mean Diurnal Range",                             width: 200,                             dataIndex: "mdr"                         }, {                             header: "Isothermality",                             width: 100,                             dataIndex: "iso"                         }, {                             header: "Temperature Seasonality",                             width: 200,                             dataIndex: "tsd"                         }, {                             header: "Max Temperature of Warmest Month",                             width: 200,                             dataIndex: "maxtwm"                         }, {                             header: "Min Temperature of Coldest Month",                             width: 200,                             dataIndex: "anntmin"                         }, {                             header: "Temperature Annual Range",                             width: 200,                             dataIndex: "tar"                         }, {                             header: "Mean Temperature of Wettest Quarter",                             width: 200,                             dataIndex: "meantwq"                         }, {                             header: "Mean Temperature of Driest Quarter",                             width: 200,                             dataIndex: "meantdq"                         }, {                             header: "Mean Temperature of Warmest Quarter",                             width: 200,                             dataIndex: "meantwaq"                         }, {                             header: "Mean Temperature of Coldest Quarter",                             width: 200,                             dataIndex: "meantcq"                         }, {                             header: "Annual Precipitation",                             width: 150,                             dataIndex: "annprec"                         }, {                             header: "Precipitation of Wettest Month",                             width: 200,                             dataIndex: "precwm"                         }, {                             header: "Precipitation of Driest Month",                             width: 200,                             dataIndex: "precdm"                         }, {                             header: "Precipitation Seasonality",                             width: 200,                             dataIndex: "precs"                         }, {                             header: "Precipitation of Wettest Quarter",                             width: 200,                             dataIndex: "precwq"                         }, {                             header: "Precipitation of Driest Quarter",                             width: 200,                             dataIndex: "precdq"                         }, {                             header: "Precipitation of Warmest Quarter",                             width: 200,                             dataIndex: "precwarmq"                         }, {                             header: "Precipitation of Coldest Quarter",                             width: 200,                             dataIndex: "preccq"                         }],                         //sm: new GeoExt.grid.FeatureSelectionModel()                         listeners: {render:function(grid){                             grid.getSelectionModel().selectFirstRow();                         }}                     });                          /*var exportCSV = new Ext.ux.Exporter.Button({                                       component: grid,                                       text : 'Download CSV'                                   });                                        grid.getTopToolbar().add(exportCSV);*/                                        var win = new Ext.Window({                                       width: 600,                                       height: 300,                                       autoScroll: true,                                       title: "BioClim Data",                                       /*tbar: [{                                         text: 'Download CSV',                                         handler: new Ext.ux.Exporter.Button({                                             //alert('Downloads the CSV');                                             component: grid                                         })                                         }],*/                                       listners:{show:function() {                                           this.loadMask = new Ext.LoadMask(this.body, {                                                msg:'Loading. Please wait...'                                           });                                        }},                                        items: grid                                         });                                 win.show();                               },                           failure: function(response) {                                   console.log('server-side failure with status code ' + response.status);                                   addToPopup(lonLat, "Server side Failure");                               }                                    //callback: myhandler                      });                  }}),             map: map,             // button options             toggleGroup: "draw",             allowDepress: false,             tooltip: "Fetch BioClim Data",             // check item options             group: "newTool"         });         actions["bioclim"] = action;         toolbarItems.push(action);


title: Reading Attributes from Attributed Relationship Class ArcObjects

text: I have an M-N relationship class with attributes. I am listening to the **`IRelationshipClassEvents.OnCreate Event`** event to see if any new relations are created. How can I read these attributes from an `IRelationship`?               private void m_editEvents_OnCreateRelation(IRelationship rel)     {         //how to get the attribute here?     }


title: Is PgRouting DARP Function still being supported?

text: I'm asking this, because I need some tool to solve a DARP problem. After researching I found PgRouting had included this functionality in one of its branches. But after installing it, getting stuck and looking for other people problems similar to mine it seems as if no one is answering DARP Releated questions. (Not just here but in general, PgRouting mailing list included) Example: How can I make Darp function of pgRouting working? And this one which has a year already hasn't been answered neither: Help making the DARP example run (I have the same doubts) So finally my question is this functionality being supported? Is anyone working with it? or it has just being deprecated?


title: Modify dataframe extent using data driven pages and arcpy.mapping

text: I'm customizing an automated mapping method using arcpy.mapping and data driven pages. Essentially I have a script that loops through each page, and modifies the location of the inset map, legend, and other map or graphic elements. I've got everything working, except for some maps, it is necessary to change the extent of the main dataframe so that the inset map, or legend does not overlap with the region of interest. As an example, the data frame may need to shift 500 metres North and East, so that there is no overlap. This is what I can't get to work. I'm using a for loop, with conditional statements to make the necessary adjustments. An example condition is:               if (pgIndex == 4):              # Set position of the Legend         elLeg = arcpy.mapping.ListLayoutElements(tempMap,"LEGEND_ELEMENT", "")[0]         elLeg.elementPositionX = 0.0186         elLeg.elementPositionY = 0.4124         # Move scalebar         scaleBar = arcpy.mapping.ListLayoutElements(tempMap, "MAPSURROUND_ELEMENT", "Alternating Scale Bar")[0]         scaleBar.elementPositionX = 0.1         scaleBar.elementPositionY = 0.1         #Move North Arrow         NorthArrow = arcpy.mapping.ListLayoutElements(tempMap, "MAPSURROUND_ELEMENT", "North Arrow")[0]         NorthArrow.elementPositionX = 0.75         NorthArrow.elementPositionY = 0.1         #Move Legend NeatLine         Neatline = arcpy.mapping.ListLayoutElements(tempMap, "GRAPHIC_ELEMENT", "")[0]         Neatline.elementPositionX = 0         Neatline.elementPositionY = 0         #Move inset map position         inset = arcpy.mapping.ListDataFrames(tempMap, "LocatorDF")[0]         inset.ElementPositionX = 0         inset.ElementPositionY = 3.25         #Set extent for main data frame (THIS IS WHAT IS NOT WORKING)         mainMap = arcpy.mapping.ListDataFrames(tempMap, "MainDF")[0]         extent_4 = arcpy.Extent(423017.450685167,4991568.14747271,438257.481165228,4976328.11699265)         mainMap.extent = extent_4


title: MapBasic Create Circle Round a Point

text: I've been tasked with creating a simple MapBasic plug-in that will loop through a table in MapInfo, zoom to the location of each point, create a circle to highlight the point, and convert the map to JPEG. I've got all the steps working except the loop (I'll post a separate question for that) and the creating a circle. This is the code I use for outputting a single JPEG;               levelOfZoom = 2     circRad = (levelOfZoom / 10)     Select * from TabName where StructureNo = sName Into Zoomtable     Set Map  Center (Zoomtable.Easting, Zoomtable.Northing)     east = Zoomtable.Easting     north = Zoomtable.Northing     Add Map layer Zoomtable     Insert Into Zoomtable (obj)      Values(CreateCircle(east, north, circRad))           Close Table Zoomtable     Close Table Selection          Set Map  Zoom levelOfZoom Units "km"     Save Window my_window_id As "U:\MapInfoOutputs\"+sName+" at "+levelOfZoom+"km.JPG" Type "JPEG" Width 5.20833 Units "in" Height 3.75 Units "in" Resolution 128      I'm getting a 'Zoomtable is read-only' error when I try and execute. If I change the `Insert Into Zoomtable (obj) Values(CreateCircle(east, north, circRad))` To; `Update Zoomtable set obj= CreateCircle(east, north, circRad)` It creates the circle, but replaces the original point. I need the original point to show with the circle round it to highlight it. Also, the circles need to be removed once the jpeg has been created, I'm not sure how to do this, but having a map with


title: How do I convert many .dxf files to one gdb file?

text: I am using ArcView and ArcCatalog 9.2 with the Data Interoperability extension and I have 159 .dxf files that are comprised of 6 files; a dessin DAO, an annotation DAO, a multipatch DAO; a point DAO, a polygon DAO, and a polyline DAO. It's all in French and I don't know the english words for some of the tools. So I apologize for any unclarity. They create a parcel layer with some building foot prints along with other features. Using a tool (Eclater un DAO; translation: Explode a DAO) in ArcCatalog that converts DAO to GDB works, but I will have to do it 159 times getting 159 different GDB when all I want is for them to be together. Will I need to combine the 159 .dxf files into one first? If so, how do I do that without CAD? Thank you for any suggestions.


title: Postgis make error on MacOS

text: I could compile it with the following command on Macbook mavericks.               ./configure --with-sfcgal=/usr/local/Cellar/sfcgal/1.0.4/bin      Then, I receive following results along with lists: > PostGIS is now configured for x86_64-apple-darwin13.3.0 make gives an error: > Can't fetch local revision (neither .svn nor .git found) Not updating > existing rev file at 12547 for s in liblwgeom libpgcommon postgis regress > raster topology loader utils extentions; do \ ... > "sfcgal_triangulated_surface_triangle_n", referenced from: _SFCGAL2LWGEOM in > liblwgeoma(lwgeom_sfcgal.o) "sfcgal_version", referenced from: > _lwgeom_sfcgal_version in liblwgeom.a(lwgeom_sfcgal.o) (may be you meant: > _lwgeom_sfcgal_version, _pg_finfo_postgis_sfcgal_version, > _postgis_sfcgal_version ) ld: symbold(s) not found for architecture x86_64 > clang: error: linker command failed with exit code 1 (use -v to see > invocation) make[1]: *** [postgis2.1.so] Error 1 make: *** [all] Error 1 Postgresql was installed using brew. I uninstalled Postgis brew installation in order to have SFCGAL functions compiled with Postgis.


title: Python version of MAP_2POINTS

text: I think I may have stumbled onto an **IDL** function which does not have a python equivalent... :'( I am working with a Whiskbroom imaging radiometer (Terra/Aqua MODIS), with a corresponding bounding box such as [ _UL, UR, LR, LL_ ] Latitude/Longitude for the final image. I also have a point located somewhere else (maybe outside the image, maybe inside) I'll call _PX_ What I want to do is:   1. Find the midpoint between _UL_ and _UR_ points, and between _LR_ and _LL_   2. I want an array of Latitude/Longitude coordinates for a line that connects those midpoints   3. For all items in the array, I want to see if the great-circle distance to the _PX_ is less than 55km Basically, I just want to see if my point (PX) fell almost at Nadir for my imager (if Nadir is within 55km of the center line)... In **IDL** I can implement this by the MAP_2POINTS function:               cpt=MAP_2POINTS(UL_Lon, UL_Lat, UR_Lon, UR_Lat, NPATH=3)     cpt2=MAP_2POINTS(LR_Lon, LR_Lat, LL_Lon, LL_Lat, NPATH=3)     nad=MAP_2POINTS(cpt(0,1), cpt(1,1), cpt2(0,1), cpt2(1,1), NPATH=2001)     ndist=fltarr(2001)     for q=0,2000 do begin        ndist(q)=(MAP_2POINTS(PX_Lon, PX_Lat, nad(0,q), nad(1,q), /METERS))/1000     endfor      Is there any simple way to do this in **Python2.7**?


title: ArcGIS JavaScript API - WFS workaround?

text: Does anybody know of an existing script or a workaround that will convert WFS returned data into something that the ArcGIS JS API can digest and display? Essentially, it would be some sort of script that either returns data that looks like it came from the REST service, or tricks the API into thinking it's an ArcGIS feature layer? I'd prefer not to reinvent the wheel if I don't have to. Thanks


title: Can't figure out wfs url filter - what am i doing wrong?

text: I'm trying to figure out how to use the filtering function on a WFS layer. Anywhere I look, there is not an example of the xml AND the query alongside it explaining how it works. In the url query below, my filter is completely ignored. (This is just an example, i have tried several modifications to this url, to no avail.) Can someone provide an example of a working url query with the filter parameter?               http://maps.co.pueblo.co.us/outside/services/pueblo_county_roads_wms_only/MapServer/WFSServer?request=GetFeature&service=WFS&typename=pueblo_county_roads_wms_only:Roads&maxfeatures=2&filter=%3CPropertyIsLike%3E%3CPropertyName%3EOBJECTID%3C/PropertyName%3E%3CLiteral%3E66580%3C/Literal%3E%3C/PropertyIsLike%3E      Thanks!


title: Error "Unable to open data source" - import .shp into GRASS

text: I'm new to GRASS and GIS in general, so please forgive my ignorance... I have problem with loading Vector Data (.shp) into GRASS I used the following steps 1.1. and 1.2: http://qgis.org/en/docs/training_manual/grass/grass_setup.html?highlight=grass In the end I have choice loaded layer (.shp) and when I clicked the "Run" button to begin the import, appear error: Unable to open data source I have all files like .shx, .shp in the same directory. I don´t know what I´m doing wrong.


title: Jsonp Callback error

text: Recently we have started migrating application from Javascript API 3.1 to 3.5 and we are finding this jsonp error coming sometimes during adding the layer.When we add the layer we also zoom to the layer .Below is the snapshot of the error I get. ![Jsonp Error](http://i.stack.imgur.com/CcXh2.jpg)


title: Performing a Python scripted tool based on whether a field is populated or not

text: My question is if I'm running the ArcMap Dissolve tool through a Python script, how can I tell it to check for attributes in one field, and if that field is empty, to perform the tool on another field. To be more specific, I would like the Dissolve tool to check for route numbers in a field called [rt_shrt_nm] (which are provided when the route name is a numerical value). But if that field is empty (which is sometimes is if the field has textual route names) to perform the tool on the [rt_long_nm] field instead. This is what I'm working with now, just to give you a better idea, even though I know its probably wildly wrong:               if [rt_shrt_nm] is null:             arcpy.Dissolve_management(OutShapesFCname, outGDB, ["rt_long_nm"], "", "", "")     elif:             arcpy.Dissolve_management(OutShapesFCname, outGDB, ["rt_shrt_nm"], "", "", "")


title: Why are towgs84 definitions missing on Ubuntu 12.04 compared to Windows?

text: Why for some of my most used CRS the 'towgs84' part of the definition of the layer spatial reference system is missing, compared to the QGIS for Windows? EPSG 3004, 23032, 23033 when used with QGIS for Ubuntu are missing the 'towgs84' definition part, that is present in QGIS for Windos. Why this? I've experimented this with QGIS 1.8 for Win and QGIS 1.8 and 1.9 Master under Ubuntu 12.04


title: Leaflet: how to change default style on a point GeoJSON Layer?

text: I need to change style on a point GeoJSON layer in a Leaflet map. I'm using the following code:               function onEachFeature(feature, layer) {                           if (feature.properties && feature.properties.popupContent) {                              layer.bindPopup(feature.properties.popupContent);                           }                          }          var myStyle = {      "color": "#ff7800",      "weight": 5,      "opacity": 0.65     };          myGeoJSONLayer = L.geoJson(myGeoJSON, {                           style: myStyle,                           onEachFeature: onEachFeature,                  });          myGeoJSONLayer.addTo(map);                               All is working but there is always the standard default blue marker on my map. Any suggestions are appreciated.


title: Is there a circle tool for line features in ArcMap 10?

text: I want to draw a roundabout in a road map in ArcMap 10. Those buggers are usually perfectly round circles. How can I do this? I am struggling to find a good way to do this. It takes ages with the Line tool and Curve tool. Can anyone help, please?


title: TIleMill multiple design feature of the label road

text: Does TileMill have possibility to design features/labels in my case based on more than one command? I have this carto CSS code               #mainroad_label[type='primary'][zoom>13], {       text-name:'[name]';       text-face-name:@sans;       text-placement:line;       text-fill:@road_text;       text-halo-fill:@road_halo;       text-halo-radius:1;       text-min-distance:60;       text-size:9.5;      But i need one type of "primary" road at the zoom level from 12! SO i have written this code:                #mainroad_label[type='primary'][zoom>13], {               text-name:'[name]';               text-face-name:@sans;               text-placement:line;               text-fill:@road_text;               text-halo-fill:@road_halo;               text-halo-radius:1;               text-min-distance:60;               text-size:9.5;         [name='Leave Road'] [zoom>=12] {           text-name: '[name]';           text-face-name:@sans;           text-placement:line;           text-fill:@road_text;           text-halo-fill:@road_halo;           text-halo-radius:1;           text-size:10; }      }      It did not show me this road at Z12 because "Leave Road" is "primary" road and obviously all "primary" roads are shown from `zoom>13`. When i have changed it to `zoom>=12` it showed me all "primary" roads but i need show only "Leave Road" primary road at this zoom level. Is there any way how can i do it? Simply i would write something like [type='primary'], [name='Leave Road'] {.....} The result will be only this one "primary" road with "Leave Road" name shown at my target zoom level. Thanks


title: Why are layers imported into a DB with Spatialite-Gui locked from editing?

text: I have a Spatialite DB which I have imported some shp files too using Spatialite-Gui. I have full read-write access to these layers via the Spatialite-Gui, however when I open these layers from QGIS I only have read access. I am able to successfully insert a new table column using the QGIS editor, but cannot modify any rows. I have read the documentation on the Sqlite page and it seems that something is holding a Reserved lock on the table. I have tried disconnecting Spatialite-Gui from the DB but the lock still remains. This problem only occurs with DBs that have been created with Spatialite-Gui. Spatialite DBs originally created in QGIS retains read-write access in both QGIS or Spatialite-Gui. I do notice that the Spatialite DBs created in Qgis are in Spatialite 3.0 could this be the problem? I am using QGIS 1.8, Spatialite 4.0, and Spatialite-Gui 1.6 TIA for your insights.


title: How to create a Winkel-Tripel map of countries with different colors?

text: I need a good quality rasterized world map, Winkel-Tripel projection, in which all countries have different (random?) colors. No borders, text, grid or anything else on the map. Is GIS the right way to go? I found a lot of map data on naturalearthdata.com. Can I create this map with QGIS? Do I need ArcGIS? I'm completely new to GIS; it looks like an exciting field. I'd appreciate any help to get started.


title: Can you script publishing a GP service?

text: I have a geoprocessing service in ArcGIS for Server 10.1. I would like to share this service with a number of external organisations that also have ArcGIS Server 10.1.   i.e I dont want them hitting my service, but for them to have the ability to easily re-publish the service at their end. Ideally, I would like to bundle this into a simple installer, which prompts end-user for some admin credentials to ArcGIS Server, and then makes use of a python script to publish the service via the new ArcGIS Server Administration API. However - To publish a GP service, the help states: > To publish a geoprocessing service, you need a result in the Results window I therefore can't think how an installer would be able to publish a GP service from a python script, as sounds like it is restricted to being published from within ArcGIS for Desktop. I also believe that it is not possible to achieve publishing a GP service from a GP package either. My question is, are my assumptions right, and **what would be the best workflow for me to share my GP service, so that my end-users can publish it to a GP service in the simplest approach possible.** i.e. I want to avoid any human error, permission issues, and minimise the time and effort it would take for the end-users to publish the service seamlessly to their ArcGIS for Server.


title: Feature comparison between QGIS and ArcGIS for Desktop

text: I am learning to use QGIS and ArcGIS for Desktop but I cannot find the differences between them. What are the differences between the capabilities of QGIS and ArcGIS for Desktop?


title: Can I access the request URL in a AGS Server Object Extension (SOE)?

text: I would like to obtain the full URL of a REST request inside of an ArcGIS Server SOE. I would like to return the full URL to an output file. But, the service may be exposed on different servers or through different urls on the same server via a proxy. I'm working with ArcGIS Server 10.1 on Windows with IIS.


title: What do I need to create web based GIS application to show GTFS data

text: I need to create several GIS based application and embed them on a web page for professional development. I need a cost effective solution and don't have big budget. Any help in right direction and to get started is appreciated (case study etc.).


title: How do I get the distance meter value between two geometries in PostGIS?

text: I have a simple question about calculating distances in PostGIS. I would like to get the distance between two geometries. I am using this sid : 4269 in meters what i am doing now is this : `ST_Distance((a.geom,b.geom)) FROM ...` but I am getting result in degrees. I think that I should work with geography but how can I cast a geom to a geography? I tried with (a.geom::geography) but i am getting an error. What can I do to get my result in meters? Thank You


title: Point size based on the scale and point color based on the attribute in SLD

text: I want to edit points size based on the scale and also I want to set up colors based on one attribute column (Quantities in ArcMap, 10 classes). I already wrote SLD rules for colors, but I don't know how to insert maxscale and minscale for the same points. For example my SLD looks like this:               <sld:Rule>               <sld:Name>0,000001 - 0,363959</sld:Name>               <sld:Title>0,000001 - 0,363959</sld:Title>               <ogc:Filter>                 <ogc:PropertyIsBetween>                   <ogc:PropertyName>attribute</ogc:PropertyName>                   <ogc:LowerBoundary>                     <ogc:Literal>1E-06</ogc:Literal>                   </ogc:LowerBoundary>                   <ogc:UpperBoundary>                     <ogc:Literal>0.363959</ogc:Literal>                   </ogc:UpperBoundary>                 </ogc:PropertyIsBetween>               </ogc:Filter>               <sld:PointSymbolizer>                 <sld:Graphic>                   <sld:Mark>                     <sld:WellKnownName>circle</sld:WellKnownName>                     <sld:Fill>                       <sld:CssParameter name="fill">#5ABA00</sld:CssParameter>                       <sld:CssParameter name="fill-opacity">1.0</sld:CssParameter>                     </sld:Fill>                   </sld:Mark>                   <sld:Size>4</sld:Size>                   <sld:Rotation>0</sld:Rotation>                 </sld:Graphic>               </sld:PointSymbolizer>             </sld:Rule>             <sld:Rule>               <sld:Name>0,363960 - 1,414790</sld:Name>               <sld:Title>0,363960 - 1,414790</sld:Title>               <ogc:Filter>                 <ogc:PropertyIsBetween>                   <ogc:PropertyName>attribute</ogc:PropertyName>                   <ogc:LowerBoundary>                     <ogc:Literal>0.36396</ogc:Literal>                   </ogc:LowerBoundary>                   <ogc:UpperBoundary>                     <ogc:Literal>1.41479</ogc:Literal>                   </ogc:UpperBoundary>                 </ogc:PropertyIsBetween>               </ogc:Filter>               <sld:PointSymbolizer>                 <sld:Graphic>                   <sld:Mark>                     <sld:WellKnownName>circle</sld:WellKnownName>                     <sld:Fill>                       <sld:CssParameter name="fill">#83CF00</sld:CssParameter>                       <sld:CssParameter name="fill-opacity">1.0</sld:CssParameter>                     </sld:Fill>                   </sld:Mark>                   <sld:Size>4</sld:Size>                   <sld:Rotation>0</sld:Rotation>                 </sld:Graphic>               </sld:PointSymbolizer>             </sld:Rule>      So, where to insert max scale and min scale rule? How to insert it? What value to set up (max zoom out must be Australia). It is pretty emercency, and thank you so much.


title: Continents and islands

text: What would you call the polygon corresponding to a coastline theme? Land mass? It corresponds to continents and islands, but I'm not sure about the best word. The opposite of the sea.


title: PostGIS 3D Extrude function

text: I selected buildings from OSM (which represents my city). I insterted a new column 'height' and I filled with a random number exactly number 2, so every building have 2metres altitude. Buildings are in 4326 and these are stored in the geom field. I would like to extrude all building with the ST_Extrude function. I tried : SELECT gid, ST_Extrude(geom, 0,0, height) as geom FROM csong_buildings; not working: ![enter image description here](http://i.stack.imgur.com/i0Few.png)


title: Prefixing graduated color labels using Python?

text: I have a series of (many) layers in an ArcGIS 10.1 map that are all symbolized using graduated colors. Each layer has the same number of breaks; however, the break values are different for each layer. For the map legend, I have to prefix the labels for each of the layers so that it looks something like this:   0 - 50%: (break values)   50 - 75%: (break values)   75 - 90%: (break values)   etc... I'm trying to find a programmatic way of inserting the prefixes into the labels. I know that arcpy gives read/write access to the labels via the classBreakLabels property. I'm a Python newbie (but learning) so I'm wondering if it's even possible to read the current labels, load them into a list (or an array?), then add the prefix values (which are static and don't change). Any suggestions? Or code snippets? Or other ways of accomplishing this? Any directions would be much appreciated.


title: Change ellipsoid reference for altitude in a DEM

text: I have DEM data stored as RGF93 (Lambert-93) for lat/lon and NGF-IGN69 for altitude. I'm trying to change the reference ellipsoid for the elevation data (ultimately to NTF, but WGS84 would be fine to understand). I achieved to reproject the file using QGis or GDAL, but I can't find a way to change the altitude values. Shouldn't gdalwarp do the trick ? Here's an example of a command I used :               gdalwarp     -s_srs "+proj=lcc +towgs84=0.0000,0.0000,0.0000 +a=6378137.0000 +rf=298.2572221010000 +lat_0=46.500000000 +lon_0=3.000000000 +lat_1=44.000000000 +lat_2=49.000000000 +x_0=700000.000 +y_0=6600000.000 +units=m +no_defs"     -t_srs "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0 +unit=m"     input.tif output      What did I miss ?


title: Style manager not working to batch style group in QGIS 2.0.1

text: I am trying to apply batch styling to a group of layers using the style manager by setting the conditions, but the symbology does not seem to be applying. This is considering the multiqml plugin is not ready for QGIS 2.0 yet. In any case, is it possible to also upon adding vector data, specify a preloaded style, i.e. qml style descriptor or sld when importing them into the layer work-space? For that matter, how do you enable batch labeling for a group?


title: Unable to get python-qgis-common on Kubuntu 14.04

text: just upgraded my Kubuntu from 13.10 to 14.04 LTS and lost all plugins. Found out that "python-qgis-common" is not installed. Tried to install it but was told that it depends on "sip-api-10.1" which cannot be installed. Searched the web and according to almost all ubuntu sites the "sip-api-10.1" is virtual package provided by "python-sip" which is installed (sip4 (4.15.5-1build1) trusty). BTW, I also had problem with libspatialite3 - had to find libspatialite3_3.1.0~rc2-2_i386 on launchpad and manually install it. Have 32bit Kubuntu, use standard QGIS repos: deb http://qgis.org/debian trusty main deb-src http://qgis.org/debian trusty main and have absolutely no idea what is wrong. I also performed autoclean after distro upgrade. Do you know any solution of this situation?


title: ArcGIS 10.1 Editor Widget

text: I am trying to use the built in editor widget in the toolkit which already contains a template picker but when i load multiple feature layers, all graphics that are in those layers are being listed in the template picker which makes it very crowded. Is there any way to manually set the templates inside the template picker? or to actually remove the template picker part? Any ideas guys? Thanks.


title: How do I completely fill polygons with text labels on a per-feature basis?

text: I would like to create a **map of regions, where each region is completely filled with its specific chunks of text**. Think a map of Wordles in the shape of regions that are assembled together into a map. This is an example (not mine) for the UK but where all UK is one region instead of divided into counties. Either approach – web map (using e.g. D3 or cartodb) or traditional offline map (made e.g. in QGIS) – would be fine. I suspect this is actually quite a hard problem (fitting text inside polygon, collision detection, automatic placement, ...).


title: Network Analyst data for the UK

text: Where can I locate data that can be used to calculate driving distances in the UK using the Network Analyst extension to ArcGIS Desktop 10? The data that came with my software (I'm in the USA) only covers North America. I see that ESRI sells StreetMap Premium which covers europe: http://www.esri.com/data/streetmap How much does this data cost? I am assuming it is expensive since they don't even list a price and I would prefer some cheaper options. It looks like you might be able to buy data from Navteq/teleatlas (referenced here http://tereshenkov.wordpress.com/2011/03/08/street-data-for-network- analyst/ ) but that they will require a lot of manual processing. Is there anything like TIGER for the UK that is provided by the government? The Ordnance Survey website didn't appear to have any streetmap style data.


title: GDAL and Python: How to "pad the extent" or otherwise handle band.ReadRaster() for area falling slightly outside the image bounds?

text: In a nutshell, I'm trying to create a bare-metal tiling utility using GDAL and Python. As you can see in the graphic below, which shows the perimeter of a raster and the OSM-structured tile index, some tiles overlap the raster's native extent. This prevents me from calling `band.ReadRaster()` on the image, as doing so throws an exception. So I'm looking for a way to handle this. I have a couple ideas, but I'm not sure how to implement them; I'll go ahead and list them: **1) Somehow "Pad" the raster extent with alpha/`NODATA` values before running my script**. Usually people want to crop collars, so I'm not real sure how to _add_ a `NODATA` collar. And I'm not fond of this approach for two big reasons. First, I would need to remember to pad every image before tiling it this way. Second, each tile index will have a different footprint, meaning each iteration would need its own padded source raster. If I choose this route, though, should I just make a huge `NODATA` VRT for the full extent of the tile index, then `gdal_merge` it with my source raster? That sounds somewhat efficient, as I could just output another VRT and avoid waiting on a big render. But how would I actually make the `NODATA` VRT?? Calling `ogrinfo` would provide the extents, but then what? **2) Somehow pad on-the-fly in the script/runtime**. This approach seems more elegant, but how can I handle `band.ReadRaster()` on an area that falls outside the image extent? Any tips/tricks? Alternatively, is there something even better I haven't thought of? Here's a graphic illustrating the problem. The tile index is in red, and you can see where the perimeter tiles extend beyond the source raster. ![enter image description here](http://i.stack.imgur.com/lI7Sb.jpg)


title: Why can't Android phones render this KML poly?

text: Our outage map works well in browsers and most mobile devices. The map can be found here: Test Outage Map However, on Android phones the perimeter polygon is rendered as black. Have tried changing the polystyle but that does not work. The outage polygon is displayed properly if the boundry is commented out. Screenshot of Android render Any ideas?


title: How to handle real time GPS data with ArcGIS Javascript API?

text: I'm developing a web app using Arcgis Server 9.3.1 and Javascript API I have a database with gps data from several vehicles and I need to display each vehicle real position in a map. I display this points in a Dynamic Layer, but I need to add an Identify action on hover, so I created a graphics layer, and every 10 seconds, I clear the graphics layer and query the database to get the new coordinates values, then, I create and add each point to the graphics layer, so I can hover over the features and ifentify them. This is already working in IE and Firefox, but in Chrome this causes a memory overflow. Does anyone have any suggestions to fix this? Thanks.


title: How to summarize the values in a column

text: I really need some help with an assignment I'm working on... In my vector field I would like to summarize all values in a column, how yo do so? If I right-click in the attribute table it is possible to choose summarize, but if I do so and select the column I want to include the values in the column is not summarized I only get a Sum_output. In the sum_output I have a new table with the selected column but no summary! In my case I have created a buffer zone and want to know the length of it. The length of the buffer zone is divided into the areas where the buffer zone is. So I have around 1000 of pieces of buffer length. How can I summarize these pieces? Thanks a lot for your help!


title: How to add Mapserver layer to Leaflet?

text: I have a layer in PostGIS in epsg:4326. Put in in MapServer:               MAP     IMAGETYPE PNG     EXTENT 60.0 56.0 61.2 57.4701     SIZE 800 600     OUTPUTFORMAT         NAME "png"         DRIVER        GD/PNG         IMAGEMODE     PC256         TRANSPARENT   ON     END         PROJECTION                "init=epsg:3857"         END     LAYER         NAME "boundary"         CONNECTIONTYPE postgis         CONNECTION "user=postgres password=postgres dbname=ekb_1 host=localhost"         DATA "the_geom from bige"         TYPE POLYGON         CLASS                 COLOR 200 120 123         END     END      END Now i add it to leaflet map:                           var map = new L.Map('map',{                     center: new L.LatLng(56.8391677245,60.7011796115)                     ,minZoom: 12                     ,maxZoom: 20                     ,zoom: 14                     ,zoomAnimation: true,                     measureControl: true,                 });                      var nexrad = L.tileLayer.wms('http://127.0.0.1/map/?map=/srv/www/localhost/public_html/map/ekb_test.map&mode=map', {                layers:'boundary',                 format: 'image/png',                 transparent: true,                 attribution: "Weather data © 2012 IEM Nexrad"             });             map.addLayer(nexrad);      And see multyple layer images on map. What can be wrong? **UPDATE** There is my /usr/lib/cgi-bin/mapserv -v output:               MapServer version 6.2.1 OUTPUT=GIF OUTPUT=PNG OUTPUT=JPEG SUPPORTS=PROJ      SUPPORTS=GD SUPPORTS=AGG SUPPORTS=FREETYPE SUPPORTS=CAIRO      SUPPORTS=SVG_SYMBOLS SUPPORTS=ICONV SUPPORTS=FRIBIDI SUPPORTS=WMS_SERVER      SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT      SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=FASTCGI SUPPORTS=THREADS      SUPPORTS=GEOS INPUT=JPEG INPUT=POSTGIS INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE


title: Does QGIS preserve Feature Class Relationships (1to1 or 1to many or Many to many relationships)?

text: Does QGIS preserve Feature Class Relationships (1to1 or 1to many or Many to many relationships) and can the related attributes be clicked and viewed through the identify feature tool?


title: How to overlay two raster layers

text: How to match two raster layers, I have one raster and a shapefile and I want to convert the shapefile into raster BUT I want it to have the same cellsize and its cells overlay (same positions) to the first raster. I already set as a 'mask' the first raster in the Environments options, but the resulting raster cells doesn't overlay with the first raster.


title: How can I automatically calculate geometry area when adding a polygon in ArcMap 10

text: Now I realize I can manually do this after the polygon feature is calculated by right-clicking the column in the attribute tables to calculate geometry or create a python script to do this but neither of these seem to allow me to automate this for new features without having to take those manual steps after the polygon is created and added to the table. I want this to happen automatically once the feature is created and shows up in the attribute table, much like the auto-generated shape_area column works. Maybe I can use those values to convert them to acres in a new column?


title: Can I create blended images with GDAL command line utilities?

text: I've written a program in Python with GDAL bindings that generates a single band georeferenced GeoTIFF as output. From there I use the GDAL commandline tool `gdaldem` to colorize that image to highlight certain classes of pixels. I'd like to also lay it on top of a basemap DEM as final output but I need to do this procedurally in my program as opposed to load both layers into a GIS tool and blend from there. Are there any GDAL command line utilities I can use to blend two GeoTIFFs together and output a pretty png? Here's a specific example with some random datasets:                 DEM                                  Pixels      ![](http://i.stack.imgur.com/Zngsv.png) ![](http://i.stack.imgur.com/SPRmr.png)                 Pixels + DEM Blended (This is what I want to do with command line tools)      ![](http://i.stack.imgur.com/Fona0.png) **Update:** I ended up using the `blend` function from PIL like this:               basemap_image = Image.open('basemap.png').convert('RGBA')     overlay_image = Image.open('overlay.png').convert('RGBA')     combined = Image.blend(basemap_image, overlay_image, 0.8)     combined.save("out.png")


title: SQL sub-queries for -say- dijkstra_sp

text: I want to combine the following 3 commands (and eventually other ones such as functions from matching.sql right after...):               routing=> select source as s from ways where gid=10213;       s        ------      8217     (1 row)          routing=> select target as t from ways where gid=64552;        t        -------      44926     (1 row)          routing=> insert into test(the_geom) select ST_UNION(the_geom) from dijkstra_sp('ways', 8217, 44726);     INSERT 0 1      So I tried many things that look like               WITH depart AS (SELECT source AS s FROM ways WHERE gid=10213),          arrivee AS (SELECT target as t FROM ways where gid=64552)      INSERT INTO test(the_geom) SELECT ST_UNION(the_geom)      FROM dijkstra_sp('ways', depart.s, arrivee.t);      but any of my experiments worked; any idea someone, please?? Edit: Thanks to dkastl, I wrote the following, certainly not optimized but still working one-line function to get a path from coordinates :               INSERT INTO test(the_geom) SELECT ST_UNION(the_geom)          FROM dijkstra_sp('ways',              (SELECT source from ways      where distance(the_geom, geometryFromText('POINT(5.697 45.27433)', 4326)) =              (select min(distance(the_geom, GeometryFromText('POINT(5.697 45.27433)', 4326))) as dist from ways where the_geom && setsrid('BOX3D(0 30, 12 60)'::box3d, 4326))),              (SELECT source                  from ways                      where distance(the_geom, GeometryFromText('POINT(5.05169 45.0482)', 4326))                      = (select min(distance(the_geom, GeometryFromText('POINT(5.05169 45.0482)', 4326))) as dist                              from ways                              where the_geom && setsrid('BOX3D(0 30, 12 60)'::box3d, 4326)                         )              )      );      Thanks!


title: Calculating driving distance from many nodes (sources) to one target using PgRouting?

text: I'm using a PostGIS / pgRouting database setup with a road network. And I want to calculate a **driving distance** and **shortest path** from **many sources to one target**. I found that **pgr_kDijkstra** can calulate the reverse thing ( Returns the costs for K shortest paths using Dijkstra algorithm, with many targets). I'm using the following query for the drinving distance :               SELECT id, id1 AS node, id2 AS edge, cost, the_geom     INTO "10min_nodes11"       FROM pgr_drivingdistance(     'SELECT gid as id, source, target,  length as cost FROM ways',     11, --source_id     10, --select #all nodes less than 10 units (minutes) of cost      false, false       )      as di      JOIN vertices_tmp pt      ON di.id1 = pt.id;      I saw this post but I didn't understand what they mean by "all_to_all".


title: Meaning of Simplify's tolerance parameter

text: The standard parameters of Douglas-Peucker's simplify algorithm are _geometry_ and _tolerance_ (e.g ST_Simplify in PostGIS). What's the meaning of the _tolerance_ parameter? I know that the bigger the value, the coarser the geometry will be. But does the number has any unit or it is just arbitrary?


title: OpenLayers 2.13 floating components resize not working

text: In my application i am using ExtJS4. After upgrade from OpenLayers 2.12 to 2.13.1 i noticed that i am unable to resize windows (floating components) that are on top of the map, like the resize action was losing focus. When i drag the window outside map div the resizing works fine. The problem does not exist in 2.12. Also when i set the window to modal, the resize is working properly. What may be causing such behaviour?


title: What is the internal format of PostGIS?

text: PostGIS stores its geometry entities as long strings:               SELECT ST_GeomFromText('POINT(32.0 34.0)',4326);                       st_geomfromtext     ----------------------------------------------------      0101000020E610000000000000000040400000000000004140      For nearby points like (32.0 34.0) and (32.01 34.01), the string has a common prefix `0101000020E6100000` which might suggest a suggest a geohash technique:               SELECT ST_GeomFromText('POINT(32.01 34.01)',4326);                       st_geomfromtext     ----------------------------------------------------      0101000020E6100000E17A14AE47014040E17A14AE47014140      Any idea what is the storage format of PostGIS geometries?


title: How to determine which tiles cover your Area of Interest?

text: I'm relatively new to using a GIS. I've been wondering how could I determine tile(s) on a selected area/point at a specified zoom level?


title: Unable to cast MapResourceLocal to Graphics.MapResource?

text: Unable to cast object of type 'ESRI.ArcGIS.ADF.Web.DataSources.ArcGISServer.MapResourceLocal' to type 'ESRI.ArcGIS.ADF.Web.DataSources.Graphics.MapResource'. MY code:               MapResourceManager MapRM = mMap.MapResourceManagerInstance;     MapResource MapRes = (MapResource)MapRM.GetResource(mGraphicsLayerName);      at the second line i am getting that cast error... and i am using Web ADF with ArcGis Server and the version is 9.2.0


title: Can custom map tiles be properly rendered via ArcObjects within my ArcGIS 10.1 Desktop AddIn?

text: Within ArcObjects for ArcGIS 10.1, we are attempting to render pre-generated map tiles as a basemap (group) layer. The tiles themselves are of the typical web mercator aux. spherical tiling scheme. Under ArcMap 10.1 the tiles do not render properly: ![incorrect rendering result](http://i.stack.imgur.com/hRXDP.jpg) and the ArcMap 10.1 UI is displaying the following drawing error: "Missing raster band wavelength properties" [Expected] rendering result under ArcMap 10.0: ![enter image description here](http://i.stack.imgur.com/VQVzt.jpg) While this is related to this (my) question here, that was a holistic approach question. Here I am specifically asking how to properly _render_ the tiles in ArcMap 10.1. We successfully accomplished the layer adding/rendering within 10.0 via the following code:               _tileLayer = new MyCustomTileLayer();    // implements ILayer among others..     _coverageBasemapLayer = new BasemapLayerClass();     _coverageGroupLayer = _coverageBasemapLayer as IGroupLayer;          // add the layer to the map     IMap map = ArcMap.Document.FocusMap;     map.AddLayer((ILayer)_coverageBasemapLayer);          ...          // wire the AfterDraw Event     IActiveViewEvents_Event avEvent = ArcMap.Document.FocusMap as IActiveViewEvents_Event;     avEvent.AfterDraw += new IActiveViewEvents_AfterDrawEventHandler(avEvent_AfterDraw);          private void avEvent_AfterDraw(IDisplay Display, esriViewDrawPhase phase) {         // force re-draw         _tileLayer.Draw(esriDrawPhase.esriDPGeography, Display, null);     }          ...          // effectively, within _tileLayer.Draw() here...     IRasterLayer rl = new RasterLayerClass();     // 'file' is the full path to the map tile PNG on the local file system, downloaded     // and/or cached     rl.CreateFromFilePath(file);          _coverageGroupLayer.Add(rl);      Presumably many things have changed within the renderer for 10.1. However can anyone provide insight into the ArcObjects changes required? Specifically how can we [properly] accomplish this within ArcObjects under 10.1? Thank you in advance.


title: gdal_retile.py tile naming scheme

text: Is there a way to define custom tile names during a gdal_retile.py process, other than the default _01_01. For instance, suppose I want to start naming my tiles at 100_100, is this possible?


title: HELP! Starting up GRASS - Error messages.

text: I am just trying to start using GRASS, but am having problems setting up. This is the error message I get:               WELCOME TO GRASS              Version 6.4.2         2012             1) Have at your side all available GRASS tutorials             2) When working on your location, the following materials           are extremely useful:           - A topo map of your area           - Current catalog of available computer maps             3) Check the GRASS webpages for feedback mailinglists and more:           http://www.grass-gis.org           http://grass.osgeo.org          Hit RETURN to continue          Starting GRASS ...     Traceback (most recent call last):       File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gis_set.py", line 35, in <module>         from gui_modules import globalvar       File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/wxpython/gui_modules/globalvar.py", line 76, in <module>         import wx       File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/python/wx/__init__.py", line 45, in <module>         from wx._core import *       File "/Applications/GRASS-6.4.app/Contents/MacOS/etc/python/wx/_core.py", line 4, in <module>         import _core_     ImportError: dlopen(/Applications/GRASS-6.4.app/Contents/MacOS/etc/python/wx/_core_.so, 2): Library not loaded: /Library/Frameworks/UnixImageIO.framework/Versions/E/UnixImageIO       Referenced from: /Applications/GRASS-6.4.app/Contents/MacOS/lib/libwx_macud-2.8.0.dylib       Reason: image not found     Error in GUI startup. If necessary, please     report this error to the GRASS developers.     Switching to text mode now.     Hit RETURN to continue...          dyld: Library not loaded: /Library/Frameworks/GDAL.framework/Versions/1.9/GDAL       Referenced from: /Applications/GRASS-6.4.app/Contents/MacOS/lib/libgrass_gis.6.4.2.dylib       Reason: image not found     /Applications/GRASS-6.4.app/Contents/MacOS/etc/Init.sh: line 635: 13465 Trace/BPT trap: 5       "$ETC/set_data"     logout          [Opération terminée]      I don't know GRASS very well. So I do not know where to begin to look to solve these problems, any help would be great!


title: Offset a Polygon defined in Lat Lng points by a given distance of miles or km

text: I am trying apply a `positive offset`/`buffer` (enlarge a shape equally from the borders) to a polygon that is made up of a series of latitude/longitude coordinates by a distance that is defined in either miles or kilometers, as apposed to a distance specified as a magnitude of the shape. My research has dug up a number of libraries that will enlarge a polygon by a given magnitude, however this magnitude is relative to the units that the polygon is specified in which isn't helpful in the case of lat/lng coordinates given the shape of the Earth and the changing proportions depending on position on Earth. I need to accomplish this server side preferably and have dabbled with arcgis's rest service which isn't working for me. Also have looked into ClipperJS and gotten close, however couldn't translate it's `delta` variable, the scale of increase into miles or km. Can anyone tell me how this might be achieved? I imagine this is not a straightforward problem but certainly one that has needed to be solved before.


title: Editing without clicking on a map - ESRI JS API

text: I'm looking to build a new application that allows staff to edit a service/layer. They will be inputting new survey benchmark data based on lat/long. All of the current samples I have seen on the ESRI site show the editing done by clicking a point on a map and then populating the fields in the popup. My goal would be to input the fields first and then have the app recognize the x,y and draw new point. This would eliminate the error of the spot being clicked not actually being the correct coordinates. Also, they will be doing the input at the office, so enabling GPS or similar in the app won't be helpful. Does anybody know of a sample or way to do this?


title: GeoServer: Shapefile or PostGIS?

text: I'm going to build a web application, which will render around ten maps using GeoServer. I will be also changing the attribute data in the application for some layers. So, what is the recommended choice for performance and deployment: A map rendered from GeoServer using Shapefile or PostGIS database ? ( or they don't really matter) Update: I don't deal with any spatial query.


title: Which conditions does the geoJSON file read from the server in?

text: I've added an alert to the function myReadFunction in order to detect function calls. As according to my observations, the function readMyFunction is called sometimes while I have been navigating on the map, not every time. The function readMyFunction has not been called when a feature is being viewable or not during map navigation. Could you please explain in which conditions the geoJSON file is read from the server?                mVectorLayer = new OpenLayers.Layer.Vector("Overlay", {             strategies: [new OpenLayers.Strategy.BBOX()],             protocol: new OpenLayers.Protocol.HTTP({                 url: 'http://localhost:56786/jlist.geojson',                 format: new OpenLayers.Format.GeoJSON({                     'read': myReadFunction,                     'internalProjection': map.baseLayer.projection,                     'externalProjection': new OpenLayers.Projection("EPSG:4326")                 })             }),             projection: new OpenLayers.Projection("EPSG:900913")         });      Thanks a lot for your helps and clarifications


title: Is there any way to reproject multiple layers in QGIS?

text: I'd like to be able to consolidate not only the layers of a QGIS project (using QConsolidate) but also the projections used by each layer. On-the-fly projection reduces performance in QGIS, and this is particularly problematic when working on a slow platform such as an Android tablet; ideally I'd like to be able to quickly recreate a project with all its layers in a single CRS and turn off on-the-fly projection. Is there a plugin that could take a selected **set** of vector or raster layers and reproject them, renaming appropriately? Is this perhaps something that could be done directly from the Python console, iterating over selected layers?


title: How to perform relate in ModelBuilder for one to many relationship

text: I'm not sure how to do this as I don't' think there are any ModelBuilder tools for this task. I'm not real savvy with Python but basically I want to relate field1 in table A to field2 in table B with a one to many relationship. So I want all records in table B that match Table A.


title: GeoWebCache OpenLayers WMS request structure

text: I'm developing a GIS app that makes use of GeoServer and GeoWebCache for serving requests, and OpenLayers for frontend. WMS requests to the GeoServer are served with no issues. But when I make a request through GeoWebCache, nothing is returned to the client. I've enabled Direct WMS Integration in GeoWebCache configuration. I'm using GeoWebCache integrated in GeoServer. My request is as shown below:               var wmslayer = new OpenLayers.Layer.WMS(        "Kenya Administrative Boundaries",        "http://localhost:8081/geoserver/gwc/service/wms?",         {layers: 'kenya_admin', format:'image/png', transparent: true},        {opacity: 1.0, isBaseLayer: false, visibility: true}     );     map.addLayer(wmsLayer);      But the layer is not loaded. Any advice on what I may be doing wrong. Thanks in advance.


title: Merge two maps with OR operation in GRASS

text: I have two vector maps. The two maps cover the same area but are split up into polygons differently (i.e. if you were to imagine one map placed above the other, the lines defining the polygons would not line up). Map1 has its polygons colored red where crime is greater than some threshold, else white. Map2 has its polygons colored red where poverty level is greater than some threshold, else white. I wish to create a single map with polygons coloured red where either crime or poverty are above their thresholds. As the polygons from each map are different, new polygons would have to be created. I'm using GRASS. Each file has an attribute table with several columns describing the polygons. v.overlay with an 'OR' operation looks like a good candidate, but I cant figure out how to specify the correct attribute column to use for the OR merge.


title: Creating Relationship classes

text: I am trying to create relationship classes in a Feature Dataset. The Relationship Class created successfully, I want when I update a record in the Origin table the Destination's would be updated. But it didn't.Any solution? OR procedure? Thank U in advance


title: how to convert qgis generated tiff images into jpg/jpeg using gdal command line utilities?

text: I want to develop a web portal that will show the images dynamically produced by qgis . but qgis provided tiff images are too big in size and it is not possible to show it on website using image viewers. That's why I have decided to convert those tiff images into jpg/jpeg to reduce in size as well as browser friendly so that I can show those images on my website easily. But I don't know how to convert tiff images into jpg/jpeg using gdal commands. Please help me in this regard. thanks in advance. Enamul


title: Creating an AttributeInspector via JavaScript API v2.7

text:
I'm trying to display an `AttributeInspector` when I click on an object from a layer. I'm following this sample, but I cannot get it to work. Here is my sample code:               dojo.require("esri.map");     dojo.require("esri.dijit.Scalebar");     dojo.require("dijit.layout.ContentPane");     dojo.require("esri.dijit.AttributeInspector-all");          var map;     var earthquakes, earthquakesFeatureLayer;     var visible;     var updateFeature;          function init() {         esri.config.defaults.io.proxyUrl = "/arcgisserver/apis/javascript/proxy/proxy.ashx";              var initExtent = new esri.geometry.Extent({             "xmin":-1089520.5238523486,             "ymin":4205238.517449941,             "xmax":5392339.474728875,             "ymax":7037689.037584678,             "spatialReference": {                 "wkid":102100             }         });         map = new esri.Map("map", { extent: initExtent });              dojo.connect(map, "onLayersAddResult", initSelectToolbar);              // National Geographic World Map         var basemap = new esri.layers.ArcGISTiledMapServiceLayer("http://services.arcgisonline.com/ArcGIS/rest/services/NatGeo_World_Map/MapServer");         map.addLayer(basemap);             // set the image parameters to toggle layer visibility         var imageParameters = new esri.layers.ImageParameters();         imageParameters.layerIds = [1];         imageParameters.layerOption = esri.layers.ImageParameters.LAYER_OPTION_SHOW;              // Earthquakes From Last Seven Days         earthquakes = new esri.layers.ArcGISDynamicMapServiceLayer(                 "http://sampleserver3.arcgisonline.com/ArcGIS/rest/services/Earthquakes/EarthquakesFromLastSevenDays/MapServer",                 {"imageParameters": imageParameters});         earthquakes.setDisableClientCaching(true);         map.addLayer(earthquakes);              earthquakesFeatureLayer = new esri.layers.FeatureLayer(                 "http://sampleserver3.arcgisonline.com/ArcGIS/rest/services/Earthquakes/EarthquakesFromLastSevenDays/MapServer/0",                 {                     mode: esri.layers.FeatureLayer.MODE_SELECTION,                                 outFields: ["objectid", "datetime", "depth", "magnitude", "region"]                 });         var selectionSymbol = new esri.symbol.SimpleFillSymbol();          earthquakesFeatureLayer.setSelectionSymbol(selectionSymbol);              dojo.connect(earthquakesFeatureLayer, "onEditsComplete", function() {             earthquakes.refresh();         });              map.addLayer(earthquakesFeatureLayer);              // add a scalebar to the map         dojo.connect(map, 'onLoad', function(map) {             var scalebar = new esri.dijit.Scalebar({                 map: map,                 scalebarUnit: 'english'             });         });     }          function initSelectToolbar(results) {         var earthquakesFL = results[2].layer;         var selectQuery = new esri.tasks.Query();         console.log(results[2].success);              dojo.connect(map, "onClick", function(evt) {             selectQuery.geometry = evt.mapPoint;             earthquakesFL.selectFeatures(selectQuery, esri.layers.FeatureLayer.SELECTION_NEW, function(features) {                 if (features.length > 0) {                     //store the current feature                     updateFeature = features[0];                     map.infoWindow.setTitle(features[0].getLayer().name);                     map.infoWindow.show(evt.screenPoint, map.getInfoWindowAnchor(evt.screenPoint));                 } else {                     map.infoWindow.hide();                 }             });         });              dojo.connect(map.infoWindow, "onHide", function() {             earthquakesFL.clearSelection();         });              var layerInfos = [{'featureLayer': earthquakesFL,             'showAttachments':false,             'isEditable': true,             'fieldInfos': [             {'fieldName': 'datetime', 'isEditable':true, 'label':'Earthquake Date:'},             {'fieldName': 'depth', 'isEditable':true, 'label':'Depth:'},             {'fieldName': 'magnitude', 'isEditable':false,'label':'Magnitude:'},             {'fieldName': 'region', 'isEditable':false, 'label':'Region:'}             ]}];                  var attInspector = new esri.dijit.AttributeInspector({                 layerInfos: layerInfos             },             dojo.create("div")             );                  //add a save button next to the delete button         var saveButton = new dijit.form.Button({label: "Save", "class": "saveButton"});         dojo.place(saveButton.domNode, attInspector.deleteBtn.domNode, "after");              dojo.connect(saveButton, "onClick", function() {             updateFeature.getLayer().applyEdits(null, [updateFeature], null);         });              dojo.connect(attInspector, "onAttributeChange", function(feature, fieldName, newFieldValue) {             //store the updates to apply when the save button is clicked              updateFeature.attributes[fieldName] = newFieldValue;         });                  dojo.connect(attInspector,"onNext", function(feature) {             updateFeature = feature;             console.log("Next " + updateFeature.attributes.objectid);         });              dojo.connect(attInspector, "onDelete", function(feature) {             feature.getLayer().applyEdits(null, null, [feature]);             map.infoWindow.hide();         });              map.infoWindow.setContent(attInspector.domNode);         map.infoWindow.resize(325,210);     }          function updateLayerVisibility() {         var inputs = dojo.query(".layer");         visible = [];         var inputsCount = inputs.length;         for (var i = 0; i < inputsCount; i++) {             if (inputs[i].checked) {                 visible.push(inputs[i].id);             }         }              if (visible.length === 0) {             visible.push(-1);         }              earthquakes.setVisibleLayers(visible);     }          dojo.addOnLoad(init);      The sample says I have to use a proxy to implement the example, but I'm not sure how to do it. Is it just a proxy issue that is getting in my way? Do you see any other problems in my code? **EDIT:** And here is the HTML:                                                 <link href="http://serverapi.arcgisonline.com/jsapi/arcgis/2.8/js/dojo/dijit/themes/claro/claro.css" rel="stylesheet" text="text/css">         <link href="styles.css" rel="stylesheet" text="text/css">                  <title>ArcGIS Project</title>
The following project is a HTML/JavaScript-based application, which uses the ESRI ArcGIS API for JavaScript v2.7
         Layers List:                     <input type="checkbox" class="layer" id="0" value="0" onclick="updateLayerVisibility();"/> Earthquakes

             **EDIT** : It seems that the `onLayersAddResult` event never fires, therefor the `initSelectToolbar` function never gets called. Why is this happening? What am I doing wrong?


title: Open source software for car navigation on iOS platform

text: I'm trying to make a car offline navigation app for iOS. Now I analyze OpenStreetMap features what to use. I would like to know if there is open source solution for the following parts of app:   * best storing of geodata - do you know Spatialite? Or use standard file system?   * rendering OSM (whatever format/own?) - vector rendering, not tile's system   * routing OSM, finding shortest/fastest path I did some experiments but if you have any experiences, i will be glad if you help me with start.


title: Different hierarchical clustering results

text: I'm running a hierarchical clustering on a sample of data using the steps below:               library(RODBC)          setwd('D:/r/cluster2')     channel <- odbcConnectExcel('cluster.xls')     data <- sqlFetch(channel, 'clust9')          y9 <- data.frame(inf=data$infest, faible=data$faible, moyen=data$moyen, fort=data$fort, lon=data$Lon, lat=data$Lat)          y9 <- na.omit(y9) # listwise deletion of missing     y9.use <- y9     y9 <- scale(y9) # standardize variables          wss <- (nrow(y9)-1)*sum(apply(y9,2,var))     for (i in 2:15) wss[i] <- sum(kmeans(y9, centers=i)$withinss)          plot(1:15, wss, type="b", xlab="Number of Clusters", ylab="Within groups sum of squares")          # K-Means Cluster Analysis     fit <- kmeans(y9, 5) # 5 cluster solution          aggregate(y9,by=list(fit$cluster),FUN=mean)          y9 <- data.frame(y9, fit$cluster)          # Ward Hierarchical Clustering          d <- dist(y9, method = "euclidean") # distance matrix          fit <- hclust(d, method="ward")           plot(fit) # display dendogram          rect.hclust(fit, k=5, border="red")      i got this results: ![enter image description here](http://i.stack.imgur.com/sDW89.png) But when i did the same steps tomorow i got different results: ![enter image description here](http://i.stack.imgur.com/FISWG.png) They are not different in everything, but there are individual that they now belong to another cluster! so i don't know why this behaviour? i'm interested in interpreting and explaning the results, so when i get different results each time, that will make my previous interpretation wrong, what can i do for now ?


title: Geoserver, Arc 10 or Windows 7 problem?

text: I wonder if you can help me resolve a problem. I have given access to a WMS served up through Geoserver to various members of staff and it works fine on Arc 10 running on Windows XP. However I get an ArcMap drawing error when the WMS is called through Arc 10 running in a Windows 7 environment. I get an ArcMAp Drawing Errors box appear with the following text within it: **One or more layers failed to draw: Geoserver Web Map Service: WMS service exceptions:Service Error Description: Rendering request would use 1328087KB, whilst the maximum memory allowed is 1000000KB** The dilemma I have is I'm not sure if this is a memory allocation problem with Windows 7, Arc or Geoserver seeing as it runs without problem on XP. If it is Geoserver, I'm presuming I would need to increase the JAI Maximum Memory; how do I do that? The web isn't being very helpful I'm afraid. Any help gratefully received, Many thanks, Red.


title: Mapnik projection for rectangular coordinates

text: I am trying to use mapnik to render indoor floorplans and am having a little trouble figuring out the SRS projection. I have geoJSON that represents the floorplan. The coordinates are just cartesian points between 0 and 1000. If I render the map without providing an SRS projection, everything looks fine _except_ the y-axis is flipped (I assume it somehow auto-detects the projection?). Now I would like to provide the parameter +axis=esu in order to flip the y-axis, but that seems to require providing an SRS projection. Any ideas?


title: ArcGIS API for Javascript - InfoTemplate on a polyline

text: I have a web map with a few points and a single polyline. I have infoTemplate(s) attached to the points and tried to replicate this to the polyline however, I can't get it to work. I looked on the ArcGIS API for JavaScript site and found an example against a polygon (hover) which is close but I need an example based on a polyline. Does anyone have a code sample that works for this request? Thanks!


title: Query for Creating a buffer around a particular point

text: I am working on a project for ambulance management using GIS. using Qgis,Postgresql,Postgis I am finding it difficult querying using St_dwithin as well as St_buffer. My aim is to find any hospitals available to a victim within a particular radius.I have the information of the victims locations in form of geometry.plz kindly help.


title: Grass 7 / ArcGIS 9.3 - tiff to raster

text: Again a Grass 7 question. I have the output of a modelling run as matlab file, which I can export as tiff and load as geotiff. OK - now can I extract the data/convert this tiff into a raster proper? I also have ArcGIS 9.3, should this be the better way (I think this would work with the raster calculator - In ArcGIS 9.3, how do I split a 3-band raster (.tiff) file and summarize the band values?, but I do not have ArcGIS at the moment to try this). Still curious if Grass 7 would be able to do this. Additional info: The data contain calculated wind speeds, but the display is as color band 0-255. So the question is how to retrieve the winds speeds? I can only see them in Matlab, but in order to correlate and spacially reference I have to export the result - as tiff... Thanks in advance, Peter


title: How does snapping in QGIS work?

text: I have read a number of posts and tutorials, even the QGIS manual on snapping. I set the snapping options under _settings>snapping options_ ; selecting the two layers that I would like snapping on, setting the mode to vertex, the tolerance of 0.00002 map units and checked avoid int. and enable topological editing. I moved the on polygon on layer one, grabbing the approximate vertex to the adjoining approximate vertex on layer two. at a scale of 1:500 it looks as though it snaps correctly, but when I zoom in to a scale of say 1:1 they are not aligned properly. Am I missing something or attempting to be to accurate? Thanks for the input.


title: Is there a way to prevent a WMX user from adding an embedded attachment?

text: I would like to allow `Linked` and `URL` attachments but **prevent** embedded attachments. ![enter image description here](http://i.stack.imgur.com/G1jww.png) The `ManageAttachments` ![enter image description here](http://i.stack.imgur.com/UviUG.png) privilege seems to be all or nothing. Does anyone else require this particular configuration? How did you accomplish it?


title: raster2pgsql: Strange error when creating a table

text: Hi GIS Stackexchange people, When I try to load a asc file into my psql setup, I get a strange error message. I have created the sql input file via               raster2pgsql -s 4236 -I -C -M gis_data.asc -F public.lightning > lightning.sql      The parameters follow exactly this link here. When I try to load the data into pgsql with               psql -d database_name -f lightning.sql      I get the following error message:               BEGIN     psql:lightning.sql:2: NOTICE:  CREATE TABLE will create implicit sequence "lightning_rid_seq" for serial column "lightning.rid"     psql:lightning.sql:2: ERROR:  type "raster" does not exist     LINE 1: ...lic"."lightning" ("rid" serial PRIMARY KEY,"rast" raster,"fi...                                                                  ^     psql:lightning.sql:3: ERROR:  current transaction is aborted, commands ignored until end of transaction block     psql:lightning.sql:4: ERROR:  current transaction is aborted, commands ignored until end of transaction block     psql:lightning.sql:5: ERROR:  current transaction is aborted, commands ignored until end of transaction block     psql:lightning.sql:6: ERROR:  current transaction is aborted, commands ignored until end of transaction block     ROLLBACK     psql:lightning.sql:8: ERROR:  relation "public.lightning" does not exist      I thought that the sql input file would instruct pgsql to create the required tables, especially since the lightning.sql file starts with               BEGIN;     CREATE TABLE "public"."lightning" ("rid" serial PRIMARY KEY,"rast" raster,"filename"        text);       **Do I need to create the tables manually? What is the required type`raster`? Where is the mistake i do to get the gis data into the pgsql database?** Thank you for your advice! PS: Interestingly, my `lightning.sql` file contains a strange bit sequence. Is that normal? I thought it would be a full ascii file with sql instructions? BEGIN; CREATE TABLE "public"."lightning" ("rid" serial PRIMARY KEY,"rast" raster,"filename"   text);               INSERT INTO "public"."lightning" ("rast","filename") VALUES                                 ('01000001001975ADBD4F55B53F1975ADBD4F55B5BF0000000000003BC0B4CA4C69FDFF514000000000000 ... ')     CREATE INDEX "lightning_rast_gist" ON "public"."lightning" USING gist                       (st_convexhull("rast"));     ANALYZE "public"."lightning";     SELECT AddRasterConstraints('public','lightning','rast',TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,,     FALSE,TRUE,TRUE,TRUE,TRUE,TRUE);     END;     VACUUM ANALYZE "public"."lightning";      **Was maybe the conversion not ok?**


How to view or import a pbf file (OSM data,binary format) in Arcgis or Qgis? Are any 3rd party applications needed?I would like to edit this file on any GIS Software.


title: How to output the ActiveView to a Jpeg image

text: I have an ArcGIS Desktop extension and I need to extract images from the map. How can I programmatically output the active view to a jpeg image?


title: Network Dataset -- can't solve route

text: I've built out numerous network datasets before, so I'm confused as to why I'm getting this error. I've built out a polyline file for the sidewalks in my community. Warning: No route from location "Graphic Pick 1" to location "Graphic Pick 2". Error: No solution found. I can solve a route on a single line, but not on any route that requires taking multiple lines. I checked all the options under 'location snap options', but that didn't solve it. Figuring connectivity might be an issue, I merged the entire file into a single feature, but I still got the same error. Does anyone have an idea of why it won't link up? Did I need to make sure the polygons only crossed each other on nodes so they intersected? I don't remember running across this problem before.


title: set values in feature class attribute table to null

text: I've created a point-data featureclass (fc) containing data from a simple .txt file using Arcpy as follows:               arcpy.CreateFeatureclass_management(out_path= outPath, out_name= 'fcTest', geometry_type= "POINT", spatial_reference= spatRef, has_m= m, has_z= z)      I've read the text file to create fields and fill each field with data:               #lstFields= data from the text file     #fNames= list of fieldnames     #fClass= list of datatypes     for i in range(len(fNames)):         arcpy.AddField_management(fcTest, fieldN, fieldT)      Let's say the .txt file contains the xy coordinates and two fields: fA and fB               with arcpy.da.InsertCursor('fcTest', ("SHAPE@XY", "fA", "fB")) as cursor:          for point in lstFields:              vals= point.split(";")             x= float(vals[0])             y= float(vals[1])             a= int(vals[2])             b= float(vals[3])             rowValue= [(x, y), a, b]             cursor.insertRow(rowValue)      So far so good. However, some values of fA and fB have value 9999, which actually means 'no data'. My question is: how do I convert the 9999 values to a 'real' NoData value? I thought               setNull(field_name)      might be it, but I can't get it to work. See ArcGIS help here. I tried something like:               with arcpy.da.UpdateCursor('fcTest', ("fA", "fB")) as cursor:         for row in cursor:             for i in range(len(row)):                 if row[i] == 9999:                     row.setNull(row[i])                     cursor.updateRow(row)     except Exception as e:         print e.message      Which gives:               'list' object has no attribute 'setNull'      I would appreciate any help. I'm new to arcpy and have based most of the above on the Programming ArcGIS10.1 with Python Cookbook by E. Pimpler. Many thanks! UPDATE! This works flawlessly:               with arcpy.da.UpdateCursor('fcTest', ("fA", "fB")) as cursor:         for row in cursor:             for i in range(len(row)):                 if row[i] == 9999:                     row[i] = None                     cursor.updateRow(row)     except Exception as e:         print e.message      Thanks @Artwork21


title: Variogram Software

text: I would like to know wich software are suitable to create Variograms and to know sill, nugget and range easily. I would prefer if this software could run in Linux. I've tried SAGA but lacks a lot of documentation, if you know some other software would be aprecciate.


title: Proximity search: ST_DWithin vs ST_Distance?

text: I have records stored in a table with latitude/longitude coordinates stored in a geometry field. I want to find all records nearby a user supplied reference point. Note "nearby" probably means less than 100km (possibly even smaller). Most examples I see use `ST_DWithin`. Is there a reason you can't use `ST_Distance`? What is the difference between using `ST_DWithin` and `ST_Distance` to do this? For example:               SELECT name, ST_AsText(coords)     FROM places     WHERE ST_DWithin(coords, ST_GeomFromText('POINT(-12.5842 24.4944)',4326), 1)      vs               SELECT name, ST_AsText(coords)     FROM places     WHERE ST_Distance(coords, ST_GeomFromText('POINT(-12.5842 24.4944)',4326)) < 1


title: Python Geoserver gsconfig create resource from Postgis store table

text: I've working Workspace and PostGIS store in GeoServer. I would like to use python and gsconfig to create new resource/new layer from existing PostGIS table, with: computing bbox,   adding abstract,   adding keywords,   setting style,   other attributes. Is it doable with gsconfig? Could you provide any examples?   The same goes with removing resources (layers).


title: Road number boxes dependent on string length?

text: I would like to have a coloured box behind each road number on my map, matching the colour of the road it relates to and scaled to the length of the road number. I understand this could be done with custom SVG symbols but it seems like it would be a time consuming process. The road numbers in question are labels attached to points. Could I then work a solution where a box is created _around_ that point as some of kind of data representation, basing the width on the string length of the road number? EDIT: I think I may have found a way to do this with Ellipse Marker symbols. If I work it out I will post an answer.


title: PostgreSQL Queries, how can I handle logic and complex mathematical functions to update a column

text: I have a Postgres Table of u and v wind velocity components as follows:               ugrid  vgrid     -2.7   -6     -2     -6.3     -1.7   -6.4      I need to calculate the speed and direction as per the JScript Psuedocode below and store the result in two new new columns in the table:               uu = U Value     vv = V Value          if vv = 0 and uu = 0 then everything is zero     if vv = 0 and uu > 0 dd = 270, ff = uu     if uu < 0 dd = 90, ff = -uu     if vv < 0 dd=Math.atan(uu/vv)*180/Math.PI     else     dd=Math.atan(uu/vv)*180/Math.PI+180;          if dd <0 dd=dd+360     ff = Math.sqrt(uu*uu+vv*vv)     returns dd as wind direction in decimal degrees     returns ff as wind speed in m/s      I have never worked with mathematical functions in Postgres. What is the best approach, can a Query handle this type of logic, or should I process the data in some other system first? If a query is possible, how would I construct it?


title: Call Identify Click Function on Dropdown Change

text: I am a web programmer with no GIS experience and I am working with a GIS specialist with no web service experience. I am using the JavaScript API 3.7 but I am not very familiar with the API. I have a select list containing unique identifiers for line segments within a layer. I can currently zoom to the location of the segment on change of the select but I would also like to highlight the selected segment on the map. I have set up a executeIdentifyTask on click of the map. I am not sure if I should call the click method and pass it the coordinates of the segment or if there is a different identify method that I should be using instead. Any help would be appreciated, Thanks


title: How to calculate correct angles for lat/lon points with the QGIS API?

text: I am trying to use qgis API to calculate the angle between two points by using function azimuth of class QgsPoint ( http://qgis.osgeo.org/api/classQgsPoint.html ), and my CRS is EPSG:4326 - WGS 84 , but I believe I am getting the angles wrong, as if QGIS were ignoring that the coordinates are Latitude and Longitude. Could it be possibly true? If yes, how can I properly calculate the angle and why would one want to make such a confusing function?   Thanks in advance


title: Drupal-Geoserver WFS error

text: > **Possible Duplicate:**   >  Access-Control-Allow-Origin error in web map. I am using Drupal and geoserver for my project. Drupal is running on localhost:80 and geoserver on localhost:8080. I want to send data from drupal to geoserver, and get it back to drupal to display it on a map. When I send the data to geoserver through the GetCapabilities URL I am not able to see the points in Layer->preview and the GML file is empty. When I view the error in firebug its giving: > XMLHttpRequest cannot load http://localhost:8080/geoserver/wfs. Origin > http://localhost is not allowed by Access-Control-Allow-Origin. can anyone tell me what this means and how to fix it, thanks.


title: Does a national(USA) interstate rest stop/rest area dataset exist?

text: And I don't care about format--worry about that later! I found some maps in Google maps, usually per state. I also have one map made my someone else, that no longer appears in the Google Maps search for user content. Maybe that user made it private, but I still can see it. The reason I ask is related to my google maps navigation question on android.se. I think the easiest way to get this info on the droid is to create or use a shared map in Google Maps. I can find Florida, but no others show up. Maybe it's a vocabulary issue? What else are are rest stops called?


title: generate a list of all themes keywords in ArcCatalog?

text: is it possible to generate a list of all themes keywords in ArcCatalog v9.3.1? similarly, i would like to do the same for place keywords as well? FYI. my office only have about 10 standalone ArcMap licenses, without any databases or fancy servers, other than a NAS server. many thanks!


title: Is toggle button's state stored in some property in Mapinfo?

text: Well can we retrieve On or OFF state of toggle button. Do we have any property of toggle button that can be read from toggle button that is ON or OFF?


title: Formatting Dates in QGIS labels

text: I would like to display a date field as a label. The data is from a shapefile. The field in question is a date field. When I use the date field for the label, it displays in the yyyy/mm/dd format. I would like to get it to the dd/mm/yyyy format. Is this possible. I know in MapInfo I would use the FormatDate function, I wondered whether something similar existed in QGIS.


title: How to change the image size (column and row number) into the same size

text: I’m using IDL program. But I have a question how to convert different row and column number (row/col) to the same size of column and row number of 36 time series NDVI images. These images grouped into three different column and row numbers. When I do the subset I used the same shapefile and Geographic Lat/long (WGS 84) projection for all of them. So I’d like to ask why col/row numbers of these images is different. Second how to correct their size into the same size (col/row number)? Thanks for your help


title: How to route over waterways?

text: There is a lot of information about (pg)routing over normal ways. Normal ways are man made and kind of ideal for routing. However i am interested in routing over waterways. In theory this should be the same as routing over normal ways, except that waterways are often not man made and so they consist out of large water bodies like lakes. Lakes are not a neat line that you can use for routing and apart from that you also have to consider water depth, bridges and locks that may or may not allow you to cross. And all this is of course again dependent on the size/depth of your boat. Is there any tooling available that tackles some of the issues i describe above?


title: Convert geodata from .rdb and .mpr files

text: I've got an COM-server and activeX client (looks like delphi app), that display some geodata. Data stores in files with .rdb and .mpr extensions. I find out that rdb is oracle database file format. There is no running oracle server on my machine, but application still works. Is there any tools which could help me convert this data to another preferred format (.shp or .tab or something else)?


title: Using Geodesic Measurements for Circle Radii?

text: I am currently developing an OpenLayers mapping site. Measurements can be made using a line tool, and an area tool. Both of these are currently set to calculate geodesic measurements as outlined in the OpenLayers API. I use geodesic measurements rather than planar measurements as during user testing people questioned the tool's measurements for distances they already knew (such as driving between towns). A new feature of the site is for a user to be able to draw a circle on the map of a set radius. OpenLayers only allows for drawing circles using planar distances, so when a user measures the circle with the geodesic measurement tool the values do not match. In the image below the circle planar radius is 10km, but the geodesic line measurement for the diameter is 12km. Clearly this will leave a user (and myself) wondering which is correct. ![alt text](http://i.stack.imgur.com/y2EaF.png) Looking at this answer it is seems most desktop GIS systems "ignore" this issue, and return planar measurements and distances. So what is the best practice in terms of user interface and accuracy to deal with planar and geodesic measurements? **Update** I found this Google example which illustrates the issue of radii and the Mercator projection: http://maps.forum.nu/gm_sensitive_circle2.html The JavaScript code to draw the circle is as follows:                   var lat1 = (PI/180)* center.lat(); // radians         var lng1 = (PI/180)* center.lng(); // radians              for (var a = 0 ; a < 361 ; a++ ) {             var tc = (PI/180)*a;             var y = asin(sin(lat1)*cos(d)+cos(lat1)*sin(d)*cos(tc));             var dlng = atan2(sin(tc)*sin(d)*cos(lat1),cos(d)-sin(lat1)*sin(y));             var x = ((lng1-dlng+PI) % (2*PI)) - PI ; // MOD function             var point = new GLatLng(parseFloat(y*(180/PI)),parseFloat(x*(180/PI)));             circlePoints.push(point);             bounds.extend(point);         }      Does this circle take the curvature of the earth into account? **Final Update** Working code posted at http://geographika.co.uk/creating-a-geodesic-circle-in- openlayers


title: Data Source: Address ranges for every block in NYC

text: I'm looking for data that defines the address ranges used on NYC blocks. On this screenshoot you can see google has noted the first addresses on each block of West End Avenue and 92nd/93rd/94th st. (658, 680 and 700 respectively). Based on that info I can extrapolate assuming that there's an even/odd side of the street and there's no unrepresented addresses to get something like this (inefficient representation used only as an example):               Street | btwn1 | btwn2 | first | last | side          -------------------------------------------------      W End | w 92  | w 93  |  658  | 678  | west/even      W End | w 92  | w 93  |  659  | 679  | east/odd      W End | w 93  | w 94  |  680  | 698  | west/even      W End | w 93  | w 94  |  681  | 699  | east/odd      Can anyone suggest a source of this type of data? Basically it's a poor mans geocode so I can convert an address (Jeryl Metz Books - 697 West End Ave) into a suggestion of how to get there (West End btwn 93/94) even without an internet connection. I'm sure this data is contained within OpenStreetMap, but I don't have a clue as to how one might extract/generate what I need. ![West End btwn 92nd and 94th st](http://i.stack.imgur.com/3JqNw.png).


title: Alternating ESRI scale bars in QGIS

text: I'm a recent convert to QGIS from arc and think it's amazing, so much more efficient and usable. It also doesn't crash constantly which helps me keep some of my sanity. I was having issues finding some features though. Mainly an alternating scale bar like in Arc, where there are smaller grading's at the start of the scale bar moving to larger near the end, (to help emphasize distance) and i think it looks cleaner.