Identifying streets of same name where split by osm

I have a large map of every street in the country. This is a vector layer in QGIS many streets are split by OSM , IE have a different OSM ID yet the street has the same name. Often these streets are split at intersecting points. Ideally I'd like to join them using Qgis and then export my shapefile My software previously looked at points in a route and determined if the street the user was on, or turning onto, continued or ended at the intersection: hence allowing me to determine road rules for the intersection. For example our main street is one road, but is split into about 10 sections. Is there an "easy" way to join streets or determine if the street is the entire street or just a section of the street ?


Red x marks along lines after adding a length field?

I added a "length" field to a shapefile of lines called "roads". QGIS puts red x's along the length of the road afterward. Why is this, and can I make these go away?


How can I secure my REST endpoint using a custom database and Basic Auth

I need to provide Google Earth users access to my REST endpoints. Currently there is a database that has the users and pwds as well as there roles. Roles determine what data in the website they can see and could be used to tell me which endpoints they have access to. I'm currently using Windows security for authentication, but this means I have to make changes to my users when ever changes in the database occur. I want to use that database for authenticating the endpoints. Google Earth only supports Basic Auth. I've gotten as far as extending IIS to invoke Basic authentication but authenticate against a non-asp.net sql database, however, I don't know how to apply what I have now to securing the endpoints. Any help/new ideas would be great


Exit script with Python

I've got a Python script for ArcGIS that I'm working on, and I'd like to have the ability to have the script quit if it doesn't have the necessary data available. I tried a straight up sys.exit() but that would give an exception in ArcMap that I'd like to avoid. I found This thread that suggests using a try block, so I made this function:               def quit_script(message):     log_msg(message) # already defined; writes a message to a file     if log_loc:         output.close() # close the file used with log_msg()     try:         sys.exit()     except SystemExit:         pass      Unfortunately, that didn't work either. Well, it doesn't make that error on ArcMap anymore, but it also doesn't, well, quit. Right now, I have the bulk of my code in an if/else statement, but that's ugly. Anybody have any other suggestions? Thanks! Brian


Using FME to convert a MapInfo TAB to GML with style attributes

I have a FME transformer that converts MapInfo TAB files to GML. What I would like is for the pen colour, pattern and width information to be brought across to the GML somehow. It really doesn't matter how it appears in the GML as the resulting file is parsed by some custom code anyway. It would probably be ideal for it to appear as an XML attribute on each feature in the file, though. I put in an AttributeExposer with the following attributes:   * mapinfo_pen_color   * mapinfo_pen_pattern   * mapingo_pen_width and I was hoping to see these in the output, but sadly not. What am I doing wrong? * * * Here's what I have currently: ![enter image description here](http://i.stack.imgur.com/aY111.png) The `Pen_Color` isn't coming through, so I added another attribute `Test` set to a constant to see whether the expression on the former was the issue. That does not appear in the output GML either.


How to split a polygon using R

I have an irregular polygon which is a spatial object. ![enter image description here](http://i.stack.imgur.com/3kQ8t.png) I want to divide it into multiple polygons (say, n=25) with equal height. To show how I want to divide it, I've plotted grids on the background using:               grid(NA,25)      If it is not possible, can anyone at least tell me how I'd be able to get the coordinates of the intersection points between the gridlines (which is not a spatial object, however) and the polygon's boundary?


Is it possible to create a geocoding and reverse geocoding service based on the shapfiles data

We have some shapfile data and we have export them to the mysql table. The table have these columns:               id name description fulladdress code lng lat ...      Now I wonder if it is possible to create a geocoding/reverse geocoding service based on these data? If yes, what jobs will to be done?


HELP with installing QGIS 1.7.3 for a mac

Could someone please help me install QGIS 1.7.3 for a Mac and possibly give step by step instructions? I have already installed QGIS 1.8 for a Mac, however I have heard that the 1.7.3 will run smoother without any problems. Thankyou


Stream Order: Strahler vs Shreve

I am writing a research paper on watershed delineation and I am looking for experiences/opinions/authoritative answers on: **What guides a hydrologist when it comes to choosing between Strahler's (1952) or Shreve's (1996) stream ordering method?** Does morphology play a role? Or the quality of the stream network extraction help decide? When do I choose either of the methods? Shreve, R. L. 1966. “Statistical law of stream numbers”. J. Geol. 74:17-37. Strahler, A. N., (1952), “Hypsometric (Area-Altitude) Analysis of Erosional Topography”. Geological Society American Bulletin, 63: 1117-1142.


Field Calculate with a function written in python within a larger python script

I'm writing a python program where I need to field calculate using a python function. The function takes two strings; one used for an if statement and the other to become an integer complete a simple addition equation. It then returns a string of that integer. I've completed the calculation on the table in ArcCatalog and it works, but in Python it doesn't and gives me the following error: ExecuteError: Failed to execute. Parameters are not valid. ERROR 000989: Python syntax error: Parsing error SyntaxError: invalid syntax (line 1) Failed to execute (CalculateField). My function:               def yearStr(age, date):         if age == 'BC':              start = str((int(date) * (-1)) + 2251)         if age == 'AD':             start = str(int(date) + 2251)         return start      Called using this:               yearStr(!Period!, !Year!)      In Python it looks like this:               arcpy.CalculateField_management(transp_table,"Start_Time","yearStr(!Period!, !Year!)","PYTHON_9.3","def yearStr(age, date):/n  if age == 'BC':/n    start = str((int(date) * (-1)) + 2251)/n  if age == 'AD':/n    start = str(int(date) + 2251)/n  return start")      I would appreciate any help! Have tried adding ' ' around field names and variables Many thanks in advance to the person who puts me out of my misery. EDIT: Managed to find a solution using search cursor will post in next couple of days


Shift between Open layer and vector layer QGIS 2.20

Can you please help me solve this issue with Open Layer Plugin. I see some shift with a vector layer of Laos water boundaries when the background is Google Physical Layer.I am using CRS EPSG: 32647 (UTM 47N). QGIS version is 2.2.0. Please see the picture attached. Thank you.![enter image description here](http://i.stack.imgur.com/nGldP.jpg)


How do I clip/mask a satellite photo to a polygon?

Given a shapefile composed of contiguous polygons, I'd like to clip (or mask) regions of satellite photo to the polygon outlines, and export each resulting raster in medium resolution for further processing. I've gotten as far as overlaying the shapefile on Google Maps imagery (in QGIS, with Openlayers plugin), but I can't find any operation to intersect vector data with raster data (presuming that the Openlayers layer is indeed treated as a raster layer). A simple solution would be great, and suggestions of other free OSX-compatible software are fine. A programmatic solution would be ideal (though not requisite) as I'd like to export many such images, each named by an attribute field in the shapefile. I hope I haven't over-explained a simple problem, but alas, I'm rather new to GIS software. Trying to move forward with this project as quickly as possible, but stymied by data acquisition... Thanks a ton!


Why does my WMS layer overlap the Google Maps v3 logo?

I'm using Google Maps v3 as a base layer and on top of that I create an overlay pointing to my WMS (streets and avenues). A few months ago all works fine, my WMS layer was not overlapping the Google Maps footer (Google logo and copyright), but today it is. Here is how it looks: ![Image Issued](http://i.stack.imgur.com/3SUrx.png) Here is the code:               var       map,       stylez = [         {           featureType: "administrative",           elementType: "geometry",           stylers: [             { visibility: "off" }           ]         },{           featureType: "landscape",           elementType: "geometry",           stylers: [             { visibility: "off" }           ]         },{           featureType: "road",           elementType: "geometry",           stylers: [             { visibility: "off" }           ]         },{           featureType: "transit",           elementType: "geometry",           stylers: [             { visibility: "off" }           ]         },{           elementType: "labels",           stylers: [             { visibility: "off" }           ]         }       ],       gmap = new OpenLayers.Layer.Google("Google Styled", {         type: 'styled'       }),       styledMapOptions = {         name: "Styled Map"       },       styledMapType = new google.maps.StyledMapType(stylez, styledMapOptions),       projection = new OpenLayers.Projection("EPSG:900913"),       displayProjection = new OpenLayers.Projection("EPSG:4326"),       bounds = new OpenLayers.Bounds(-57.81887, -25.53276, -57.24308, -25.17649).transform(         displayProjection,         projection       );          // avoid pink error tiles     OpenLayers.IMAGE_RELOAD_ATTEMPTS = 3;     OpenLayers.Util.onImageLoadErrorColor = "transparent";     OpenLayers.ImgPath = "http://js.mapbox.com/theme/dark/";          function init() {            var         options = {           controls: [             new OpenLayers.Control.Navigation(),             new OpenLayers.Control.PanZoomBar(),             new OpenLayers.Control.LayerSwitcher(),             new OpenLayers.Control.MouseDefaults(),             new OpenLayers.Control.KeyboardDefaults()           ],           projection: projection,           displayProjection: displayProjection,           units: "meters",           numZoomLevels: 12 /* 18 */,           minZoomLevel: 6,           maxResolution: 156543.033928041,           maxExtent: new OpenLayers.Bounds(-2.003750834E7,-2.003750834E7,2.0037508345578495E7,2.0037508345578495E7)         };            map = new OpenLayers.Map('map', options);            // base layer       var tiled = new OpenLayers.Layer.WMS(         "Roadmap",         "http://foo.com:80/geoserver/gwc/service/wms",         {           LAYERS: 'Foo',           STYLES: '',           format: 'image/png',           tiled: true,           transparent: "true",           tilesOrigin: map.maxExtent.left + ',' + map.maxExtent.bottom         },         {           isBaseLayer: false,           opacity: 1.0,           visibility: true         }       );            tiled.events.register("loadstart", tiled, function() {           console.debug("Foo, Load Start");       });            tiled.events.register("loadend", tiled, function() {           console.debug("Foo, Load End");       });            map.addLayers([gmap, tiled]);            map.zoomToExtent(bounds);            gmap.mapObject.mapTypes.set('styled', styledMapType);       gmap.mapObject.setMapTypeId('styled');          }


Unable to create PopUp in OpenLayers

I am unable to create Popup in Openlayer. I have followed the examples and my code is also correct but I am not sure where am I going wrong. Here is the link. Here is my entire code,               <!DOCTYPE html>                           <title>iPlant Map</title>             <link rel="stylesheet" href="openlayers/theme/default/style.css" type="text/css">                          <
iPlant Map

                 Position


What's the difference in using OpenLayers and Mapbox?

Hey guys I found this tilemill to develop WebGIS application using MapBox. I wanted to know how is it better than OpenLayers and if I have to develop which would be a better options to go for, as I will be developing a complete GIS Architecture.


how to get the boundaries of a city in KML from google maps

im trying to write a google maps application. and i wanted to only highlight one city in the world ( im using the masking method explained: https://github.com/vasile/geomask) but now im bugged with the question of how to only show the boundaries of the city as displayed in google maps. can any one suggest how to extract the boundaries of the city from google maps in KML file so that i can use the same. in my project.


How to create dynamic sub-menu in ArcMap?

I created a static menu contains 3 static items like in this link create menus, I want to add a sub-menu and get the menu's items from database like in this image: ![](http://i.stack.imgur.com/f9zhb.png) I already did the left part with static items, but i couldn't add the sub-menu like in the right side. I tried to do what in this link create dynamic menu, but it's not working. EDIT: Menu code:               public sealed class TestMenu : BaseMenu, IRootLevelMenu {      public ReportingMenu() {                 AddItem("Reporting.Desktop.UI.Commands.ReportSuitesCommand");// commands                 AddItem("Reporting.Desktop.UI.Commands.ReportsCommand");                 BeginGroup(); //Separator                 AddItem("Reporting.Desktop.UI.Commands.AboutReportingCommand");             }                  public override string Caption {                 get {                     //TODO: Replace bar caption                     return "Reporting";                 }             }             public override string Name {                 get {                     //TODO: Replace bar ID                     return "ReportingMenu";                 }             }         }      multi item code:               public class HowToMultiItem : IMultiItem, IMultiItemEx         {             private ESRI.ArcGIS.Controls.IHookHelper _mHookHelper;                  //Step 2.             public string Caption {                 get {                     return "Layer Zoom";                 }             }             public string Name {                 get {                     return "CSNETSamples_HowToMultiItem";                 }             }             public int HelpContextID {                 get {                     return 0;                 }             }             public string HelpFile {                 get {                     return "";                 }             }             public string Message {                 get {                     return "MultiItem single message";                 }             }                  //Step 3.             public int OnPopup(object hook) {                 _mHookHelper = new ESRI.ArcGIS.Controls.HookHelperClass {Hook = hook};                 return _mHookHelper.FocusMap.LayerCount;             }                  //Step 4.             public string get_ItemCaption(int index) {                 var layer = _mHookHelper.FocusMap.Layer[index];                 return "Zoom to " + layer.Name;             }             public bool get_ItemEnabled(int index) {                 //Enable zoom in only when the layer is visible.                 var layer = _mHookHelper.FocusMap.Layer[index];                 return layer.Visible;             }             public int get_ItemBitmap(int index) {                 return 0; //Not implemented.             }             public bool get_ItemChecked(int index) {                 return false;             }                  //Step 5.             public void OnItemClick(int index) {                 //Get the layer and its spatial referenced extent.                 var layer = _mHookHelper.FocusMap.Layer[index];                 var env = layer.AreaOfInterest;                 //Zoom to the extent and refresh the view.                 var activeView = _mHookHelper.FocusMap as                     ESRI.ArcGIS.Carto.IActiveView;                 if (activeView == null) return;                 activeView.Extent = env;                 activeView.Refresh();             }                  //Step 6–Optional interface implementation.                  public string get_ItemMessage(int index) {                 var layer = _mHookHelper.FocusMap.Layer[index];                 return "Zoom to " + layer.Name + "layer.";             }             public int get_ItemHelpContextID(int index) {                 return 0;             }             public string get_ItemHelpFile(int index) {                 return "";             }         }


How to Merge Raster Datasets of Different Resolutions?

I have two DEMs in .tif format and I am looking to merge the two. The main problem with merging these are that they are different resolutions (one is 10m and one is ~20m). Also, another problem with merging the two is that they overlap and I am not sure of the most effective way of dealing with that overlap. Obviously I want to keep the values of the better resolution. So what are the best methods here?   * Is there a python script available that works best with merging rasters? I would like to automate this with a script if possible.   * Which process of resampling is best for DEM data? Bipolar, Cubic... ?   * What is the best way to merge overlapping DEMs?


Slippy map format

When do you decide which to use? I've started using Geoserver and it seems like map files are downloaded to the browser (using Leaflet) in the WFS format. When I look at OSM, Mapbox etc, they use what I believe is the slippy map format: `/zoom/x/y.png` 256x256 tiles. Should Geoserver be configured to use slippy map format? And will using WMS prevent slippy map interactions like Google maps and Leaflet? What are your recommendations if I'm using Geoserver and Leaflet and want to serve some map tiles.


GeoServer Layer Preview Opens OpenLayers instead of GeoExplorer

I am using OpenGeo 3.0 and trying to get GeoExplorer to launch through the Layer Preview. I imported PostGIS tables into GeoServer and created a layer. However, when I go to view the layer and select "GeoExplorer-->Go", it opens OpenLayers. Is there something I need to configure?.


OGR Error : Silly size requested

I'm scripting OGR functionality in order to filter some mapinfo tables and convert them into sqlite files.               ogr2ogr -f "SQLITE" newtest.sqlite test.tab -sql "SELECT cols FROM test WHERE .."      where TAB file is about 2GB size (including ID,MAP,DAT & TAB extensions within) I receive error like this :               ERROR 1: CPLRealloc(-1532233231): Silly size requested.      What does it mean exactly,how to resolve this ?


How to share ArcMap shapefiles with QGIS?

I'm trying to send a ArcMap 10 shapefile for a QGIS user and something happens with attribute table during the process. Columns will be rearranged, hidden digits will appear to coordinate columns and x-coordinate will be changed to y-coordinate and vice versa. Can anyone tell is there any way to prevent this and get all the fields right in QGIS? I'm sending all the files assosiated with .shp-file by zipping those together to one file (and then as email attachment). Files in the zip-file are .dbf, .prj, .sbn, .sbx, .shp, shp.xml and .shx.


Display ESRI Shapefiles in a map in a browser. (Not in google maps)

I have a series of shapefiles that I need to render as a map on a webpage. I have an ASP.Net server, although I would prefer something that will work on Linux servers. I would like to find an easy to set up opensource or freeware solution if at all possible. This cannot be overlaid over Google maps. It is just a map of the county I live in and that is all that needs to be shown. The shapefiles consist of the streets, boundary lines, street names, etc. Each one is a separate shapefile. This map will need to be zoomable as well. Thank you in advance. P.S. I do not have a copy of ArcGIS so I cannot do any conversions with it. Although I have QuantumGIS installed. It's server is giving me issues though.


MRT (MODIS Reprojection Tool) error: 'Input file does not exist'

EDIT: Solved, see my answer. I've downloaded 8 MOD44B granules from Reverb (using the FTP script method), and am trying to use MRT to mosaic them. I'm falling at the first hurdle though - when I add the 4 files (4 from 2000, 4 from 2007) I get an error saying firstly `One of the input files does not exist or the inputs are not the same data product type` then `Input file does not exist!`. If I add files individually, I get an error saying `Error in Module: ReadHrFile, Opening Input Header File`. Has any one had this error before, or know how to fix it? I've updated Java and reinstalled MRT. Versions:   * MRT: 4.1 (March 2011)   * Windows: 7   * Java: 6 (Updated to 7, still doesn't work) I would use the online interface, but I need MOD44b which isn't available with that for some reason.


openlayers attribute substitution failing on kml file

I'm trying to read in mile markers for a marathon from a .kml file and use the MILE attribute to set the label. My layer code looks like this:               // add layer     mile_markers_layer = new OpenLayers.Layer.Vector('Mile Markers',        {         projection: new OpenLayers.Projection('EPSG:4326'),         protocol: new OpenLayers.Protocol.HTTP({             url: 'Mile_Markers.kml',             format: new OpenLayers.Format.KML({                 extractAttributes: true                          //extractStyles:  true             })         }),         strategies: [new OpenLayers.Strategy.Fixed()],                 style: {                     'fillColor': '#888888',                     'fillOpacity': 1,                     'fontColor': '#000000',                     'fontFamily': 'arial, sans-serif',                     'fontSize': '.9em',                     'fontWeight': 'bold',                     'label': '${MILE}',                     'pointRadius': 8,                     'strokeColor': '#ff0000',                     'strokeWidth': 3                 }     });      I looked in firebug and all the features of the mile_markers_layer have attribute of `MILE` with a value of "1" thru "26" (as expected). Layer is displayed with a label of `${MILE}` \- by this I mean all 26 mile markers are displayed and each has a label of `${MILE}` I was expecting the value of MILE instead. I'm a newbie, so I could be doing something extremely stupid. thanks Dennis The offending kml file is at: https://docs.google.com/open?id=0BxhE8rRos0u1X2ExNldTdWwwWkk Thanks for looking! Dennis Thanks for the suggestion, but every mile marker ended up with ${MILE.value} for a label.


MySQL GeoSpatial function for ST_Distance_Spheroid? Metric used in return Type of GLength?

>   * Is there a ST_Distance_Sphere, ST_Distance_Spheroid equivalent function > in MySQL's GeoSpatial Support? > >   * If not, any way to emulate this? > > * * * > >   * What is the metric used in return type of GLength (length of the > LineString)? > >   * The manual for GLength says that _GLength() is a nonstandard name. It > corresponds to the OpenGIS Length() function._ But I couldn't find any Docs > for OpenGIS Length(). All I found was for ST_Length are they the same? > > * * * > >   * Also for accuracy do I need to specify the SRID(4326)? How Do I do this > in MySQL? > >


How to run a simple python script for QGIS from outside (e.g. Sublime Text)?

I would like to run a few small and simple python scripts for QGIS (Mac) from "outside" of QGIS. With outside I mean in this context, either the normal os command line (terminal.app) or even better, directly out of Sublime Text (text-editor), but definitely not via the inbuilt QGIS python console. I've read through various tutorial e.g. http://www.qgis.org/pyqgis- cookbook/intro.html#python-applications and I am able to get a reference to the QGIS app, but unfortunately not to `qgis.utils.iface` or something else deeper. This little code snippet should for instance print out the name of the active layer ... here is what I have:                   import sys         sys.path.append("/Applications/QGIS.app/Contents/Resources/python")              from qgis.core import *         import qgis.utils              print "helo" # console output: helo              QgsApplication.setPrefixPath("/Applications/QGIS.app/", True)         QgsApplication.initQgis()              print QgsApplication # console output: <class 'qgis.core.QgsApplication'>              print qgis.utils.iface # = console output: none              aLayer = qgis.utils.iface.activeLayer()         print aLayer.name()              QgsApplication.exitQgis()      Please don't reply me now how to create a plugin or something else "heavy". I am just looking for a quick an easy way to shoot scripts out of a comfortable text-editor to QGIS. Many thanks!


Getting json feature id on OpenLayers

I'm projecting some features using json format on OpenLayers. I need to get the id from the feature (not the OpenLayers' id but the one in the json file). The file looks something like this.               {"type": "FeatureCollection","features": [{ "type": "Feature", "id": 7655199, "properties": {"power":"line","wires":"NA","frequency":"NA","voltage":"380000","operator":"NA","cables":"6"}, "geometry": {"type":"LineString", "coordinates":[[8.5310699,49.8531139],[8.5289255,49.8563936],[8.5268405,49.8594383]]}} etc..      Now in OpenLayers I'm trying to do something like this.               var feature_col = gjson.read(ReadFile(json_path + file[i]) , "FeatureCollection");     //Readfile is just a function for reading the file, what a surprise... :P     var power_res = feature_col[0].attributes.power;     var id_res = feature_col[0].id) //for the first feature      So the value in power_res gives line which is correct. However, the value in id_res is something like OpenLayers.Feature.Vector_362 and I need the id in this case to be 7655199. Thanks!


ArcGIS Javascript API - How to submit a defect?

I would like to know the best way to notify ESRI of a possible defect in the JSAPI and gis.stackexchange seems like the best mechanism for doing so (Derek reads it and fixes it). That being said, what is the _proper_ way to notify ESRI of a defect? The defect in this case is posted on jsfiddle and can be seen by placing the mouse over the map. Notice the onMouseMove event continuously fires. This required a workaround on my part because I was using the onMouseMove event to re-project (esri.config.defaults.geometryService.project) coordinates and if the user left the mouse over the map it was eating up bandwidth. Documentation states: > Fires any time the mouse pointer moves over the map region. A common use for > this event is to show the current x,y coordinate of the map as the user > moves the mouse pointer. Hence the claim that this is a defect. This is true in 2.8 and 2.9 (beta)


other than GrADS and NCL, any other map/chart scripting and plotting software dedicated for climate and meteorology?

NCL is quite good indeed, but I would like to know whether there are other options.


Osgeo live - Geonetwork install issues

I have a machine running OSGEOLIVE 6 Geoserver and most packages work without any issues but I can't get Geonetwork to run. The ApplicationsGeoserversGeonetworkStart seems to work fine but then I get an unable to connect on localhost:8880/geonetworks/srv/en/main.home The Geoserver is on localhost:8082 and I can login to it but can not get to the geonetwork home page. Osgeo seems to have also installed tomcat in localhost8080 -i don't need this. I tried the following as well   1. Install Geonetwork2.8 beta (same issue with login to the home page). When I try sh start-geonetwork.sh I get a rm/mv error (see below)   2. Install 2.4 using git clone instructions from geonetwork-opensource.org/download.html but I can't run the final command. See below. Any ideas? I just want to open geoserver and geonetwork (also postgis) best,               a@a-beta-osgeo:~$ cd master     a@a-beta-osgeo:~/master$ git checkout master     Already on 'master'     a@a-beta-osgeo:~/master$ git submodule update --init     Submodule 'docs' () registered for path 'docs'     Submodule 'gast' () registered for path 'gast'     Submodule 'geoserver' () registered for path 'geoserver'     Submodule 'installer' () registered for path 'installer'     Submodule 'maven_repo' () registered for path 'maven_repo'     Submodule 'release' () registered for path 'release'     Submodule 'web-itest/webspecs' () registered for path 'web-itest/webspecs'     a@a-beta-osgeo:~/master$ mvn clean install     The program 'mvn' can be found in the following packages:      * maven      * maven2     Try: sudo apt-get install <selected package>     a@a-beta-osgeo:~/master$ sudo mvn clean install     [sudo] password for atgis:      Sorry, try again.     [sudo] password for atgis:      sudo: mvn: command not found     a@a-beta-osgeo:~/master$ sudo apt-get install      Reading package lists... Done     Building dependency tree            Reading state information... Done     You might want to run 'apt-get -f install' to correct these.     The following packages have unmet dependencies:      grass-dev : Depends: grass (= 6.4.3rc2-2~precise) but 6.4.2~upstream2-3~precise3 is installed     E: Unmet dependencies. Try using -f.     a@a-beta-osgeo:~/master$ sudo apt-get -f install      Reading package lists... Done     Building dependency tree            Reading state information... Done     Correcting dependencies... Done     The following packages were automatically installed and are no longer required:       libdapserver7 libdap-dev libxerces-c2-dev linux-headers-3.2.0-33-generic-pae       libjasper-dev libgeos-dev linux-headers-3.2.0-33 libsqlite3-dev unixodbc-dev       libmysqlclient-dev libpng12-dev libgif-dev libhdf5-serial-dev libfltk1.1       libxml2-dev uuid-dev libpq-dev libgdal-dev libltdl-dev libspatialite-dev     Use 'apt-get autoremove' to remove them.     The following extra packages will be installed:       grass     Suggested packages:       gpstrans     The following packages will be upgraded:       grass     1 upgraded, 0 newly installed, 0 to remove and 29 not upgraded.     1 not fully installed or removed.     Need to get 0 B/21.8 MB of archives.     After this operation, 53.7 MB of additional disk space will be used.     Do you want to continue [Y/n]? Y     (Reading database ... 212956 files and directories currently installed.)     Preparing to replace grass 6.4.2~upstream2-3~precise3 (using .../grass_6.4.3rc2-2~precise_i386.deb) ...     Unpacking replacement grass ...     Replacing files in old package grass-core ...     dpkg: error processing /var/cache/apt/archives/grass_6.4.3rc2-2~precise_i386.deb (--unpack):      trying to overwrite '/usr/lib/grass64/bin/xganim', which is also in package grass-gui 6.4.2~upstream2-3~precise4     dpkg-deb: error: subprocess paste was killed by signal (Broken pipe)     Processing triggers for menu ...     Processing triggers for man-db ...     Errors were encountered while processing:      /var/cache/apt/archives/grass_6.4.3rc2-2~precise_i386.deb     E: Sub-process /usr/bin/dpkg returned an error code (1)     a@a-beta-osgeo:~/master$ ^C     a@a-beta-osgeo:~/master$


OpenLayers: How to add permalink long/lat anchor to links outside of map on page

I would like to extend the funtionality of OpenLayers that appends, on-the-fly as the user zooms and pans the map, the longitude/latitude parameters as a anchor to the URL in the browser's address bar, e.g.               http://mywebsite.com/#zoom=12&lat=50.90306&lon=-1.38329&layers=BT      I would like certain other links on the same page, outside of the map, to have the the anchor `#zoom=12&lat=50.90306&lon=-1.38329&layers=BT` appended to them. The reason is so that when these links are clicked on, the map will be loaded at the same pan and zoom position. These links contain other 'arguments' that determine what data appears in the overlay. The site runs on the Drupal 7 content management system, with the OpenLayers module and these other arguments control a Drupal View to control what is displayed in the overlay. So my approach is to implement some self-contained javascript/jquery that doesn't change openlayers javascript (to avoid overheads of burdening the authors with testing, commits to the source etc) but instead is loaded after openlayers javascript, this self-contained code would at run time 'register'/'bind' as necessary to be able to listen to pan and zoom events and obtain the longitude and latitude anchor data and add these to the certain other links. I would identify those certain other links with a class attribute in the html `<a />` element and use the JQuery DOM functions to search for the class and prepend the anchor. What I need is to know how the OpenLayers code is updating the location field/address bar with the # anchor in order to figure out how to write some separate code to listen for the pan zoom events. Please share your thoughts. Thanks.


OpenGeo displays layer in openlayer, but not geoexplorer

My test server updated today I think it was a php update, and since I can no longer view data in geoexplorer, only openlayers. The geoexplorer can not see my local server. any hints? Update: Firebug says:                  TypeError: e is undefined     ...&&!Ext.isEmpty(f.title)&&(e.set("title",f.title),a.setValue(e.get(d))):c.remove(...           GeoExplorer.js (line 1535)      Another error from Geoexplorer:               Non-existing source 'local' referenced in layer config.      Further information from the geoserver log org.geoserver.wms.capabilities.GetCapabilitiesTransformer$CapabilitiesTranslator.handleLayer(GetCapabilitiesTransformer.java:899) ... 87 more 2013-07-19 13:16:12,458 INFO [geoserver.wms] - Request: getServiceInfo 2013-07-19 13:16:12,459 ERROR [geoserver.ows] - org.geoserver.platform.ServiceException: Could not determine geoserver request from http request org.geoserver.platform.AdvancedDispatchFilter$AdvancedDispatchHttpRequest@1965791


QGis Cloud Real Time Data

Through Qgis 1.8 i have a layer whose data are updated every 3 seconds.But the question is whether it is possible that QGis Cloud represents dynamically these changes. Thank you very much


improve Google map loading time?

My website has a number of pages that show a Google map with a bunch of markers, here's an example. As you can see, the maps take a long time to load and I'm looking for ways to improve this. I was hoping to use GeoWebCache to cache the map tiles on the server, but I was informed that this would violate the terms of use for Google maps. The code that I use to display a map and add a marker is appended below. It's a pretty straightforward usage of the Google Maps V3 JavaScript API, so I don't think there's much scope for optimizing it. Are there any obvious steps I could take to reduce the map-loading time?


Relevant data for cost attribute to calculate route time

I have imported OSM data into PgRouting database. I would like to use driving_distance() function to calculate route time. How to add relevant time cost into database ? Thanks


Arcpy ERROR 000732 when ApplySymbologyFromLayer

I want to split a shapefile by certain attribute and then apply a symbology layer, finally output it as a KMZ file. R1 = outpath + 'R1.lyr' arcpy.MakeFeatureLayer_management(fc, R1, "WHERE CAUSE") arcpy.ApplySymbologyFromLayer_management(R1, symbologyLayer) arcpy.LayerToKML_conversion(R1,outpath+'R1.kmz') It raised error: arcgisscripting.ExecuteError: Failed to execute. Parameters are not valid. ERROR 000732: Input Layer: Dataset R1 does not exist or is not supported Failed to execute (ApplySymbologyFromLayer).


How to find out if OpenStreetMaps contains low emission zones?

Since I never really digged into the data contained in OpenStreetMap I ask myself what would be the recommended way to identify layers or tags within the overall data. In my current use case I want to find out whether there is spatial information about the boundaries of low emission zones. First of all I looked into the taginfo-master database and filtered the table `popular_keys` by `%zone%`. ![SQLite Manager / taginfo-master filtered by %zone%](http://i.stack.imgur.com/4rDQe.png) As shown in the screenshot this results in 11 keys from which none indicated what I am searching for - but still it is possible that a more general tag name such as `zone` is used. While reading another post I get to know about JOSM and installed the latest version. Since it is not recommended to download the whole dataset in order to avoid torturing OSM servers I would like to know: How can I only load OSM data tagged with `zone` for a specific city such as `Munich`? Is QGIS a better tool for this? In the end it would be nice if I can somehow download spatial information which I can convert into common formats such as ESRI ShapeFiles.


QGIS 1.8 will not load world file

My maps georeference just fine in QGIS and I save them to jpg with jpgw. The jpg image loads but the jpgw file does not. I keep getting the error message "this is an unsupported file", or words to that effect. I have tried bmp and bmpw, tif and tifw and tiff and tiffw with the exact same result. I got the same result with Global Mapper but surprisingly, OK Maps loaded both the image file and the world file perfectly. My GPS receiver tracks accurately with OK maps also. Please help any way you can. I really enjoy using QGIS. Thanks, Jim


How to create polygon that goes over +/-180 (dateline) in ArcMap?

I want to create a polygon feature that will go through the dateline (180 long.) in ArcMap. I am using WGS_1984_Web_Mercator_Auxiliary_Sphere (WKID: 3857) for my data frame and my polygonal feature class is defined with the same projected coordinate system. My feature classes are stored within enterprise geodatabase 10.2.1 (SQL Server Geometry spatial type). However, the same behaviour can be reproduced with the file geodatabase feature class. From Esri help: > Data in a geographic coordinate system is treated, in a sense, as two- > dimensional, at least to the point that there are edges, just as in a > projected coordinate system. Generally, the edges are at -180/+180 for the > east–west extents, and -/+90 for the north–south extents. These extents are > called the "horizon" of the coordinate system. > > In a geodatabase, a feature that crosses an edge will be clipped to the > horizon. If you are editing features in ArcMap, the data frame's coordinate > system will affect what happens to a feature that crosses the +/-180° line. > If the data frame is using a geographic coordinate system, that is, the +180 > and -180 lines are separated, a feature that crosses outside the horizon > will be clipped to the horizon. **If the data frame is using a projected > coordinate system, with the +/-180° line inside, a feature that crosses it > will be split into multiple parts**. I am able to create a polyline feature that crosses the 180° line. A multipart feature is created just as outlined in the Help. When I run the Multipart To Singlepart GP tool, I can verify that my line was originally multipart and is represented in the output feature class as two singlepart features. However, this doesn't happen when I am trying to create a polygon feature crossing the 180° line. I am getting the error message: > Create The Create Feature task could not be completed. The coordinates or > measures are out of bounds. Why don't polygons get splitted just lines did? Is there any workaround that would help me to achieve my goal?


How to merge two polylines in code?

We are currently working on a project where we have to merge sections (polylines) to a route, but the big question is how you do this. We have currently tried adding those two as segmentcollections but it gives some rather wicked results.


Why do people pronounce GDAL (Geospatial Data Abstraction Library) like goo-dal?

Ok, a bit of a silly question, but: I've heard people pronouncing the library GDAL (Geospatial Data Abstraction Library) like goooo-dal. My natural inclination was to pronounce it geeee-dal, as in the letter 'g'-dal. Any explanation on how goo-dal came to be popular? Was goo-gle involved?


Tabulate Area using GDAL and Numpy

How could I perform a `Tabulate Area` using GDAL and Numpy? I am trying to summarize how many pixels of each class in a Landcover raster are within each polygon of a Census Blocks shapefile. I've tried playing with `numpy.histogram2d`, but didn't have any luck.


How to get XY coordinates and field values for each row in a SearchCursor?

I am trying to enumerate each row of a points shape file and get a certain field's value (elevation) and the xy coordinates. But I see I can only use:   * arcpy.SearchCursor to get the field; OR   * arcpy.da.SearchCursor(fc, "SHAPE@XY") for the xy coordinates Can I do both at once?


Geoserver REST missing datastores

When I query Geoserver at               http://localhost:8080/geoserver/rest/workspaces/mike/datastores.xml      I get no results, but I definitely have a data store under mike/. Moreover, the stock geoserver installation came with a workspace called "nurc" with 4 datastores and none of them are listed either. What's the deal? Viewing isn't a big issue, but when i go to delete a store via rest, it can't seem to find it. Again, the datastore is clearly there per the Geoserver web GUI. Thanks, Mike


Styling layer popup within arcgis.com

while sophisticated styling of popups (infowindows when clicking on a point) are probably better in css, are there any options for styling infowindows within an arcgis.com webmap? I really just want a less-technical user to be able to style the background, header color, and basic layout. Any ideas are appreciated.


How to join USDA SSURGO soil table data to features in QGIS?

I am trying to make a soil map using QGIS. The only problem is that I cannot join table data to features in QGIS. In ArcGIS for Desktop you just open the folder and then join the tabular table with the spatial data table. Does anyone know how to join table data to features in QGIS? I had some suggestions, one is using Python, but I do not know how to use it very well. I am open to all ideas including GRASS, programing and any plugins.


Calculate Width and Height of a Google Static Maps Image

Do you know how to calculate the width and height (in degrees minutes) of a Google Static Maps image? I’ve found it extremely difficult. The only information you have is:   \- The centre of the image (as a Longitude Latitude coordinate)   \- The zoom level of the satellite (1 >= zoom <= 25)   \- The width and height of the image in pixels See this image for a visual explanation of the problem: ![enter image description here](http://i.stack.imgur.com/Bg8BQ.jpg) And from this limited information its very difficult to calculate the image dimensions in degrees minutes, mainly because the zoom level heights (in feet or metres) is not exposed by google AFAIK. AND because you cant move up on the y axis at a constant value because you have to take into account that the earth is a sphere and not a flat surface so moving up at a constant value wont work (which is what all the answers on Stack Overflow provide). It appears meridional points is a solution but the math involved is completely beyond me. Any advice would be greatly appreciated.


Display attributes of feature with clustering

I have one KML layer that has roughly 900 features and I added OpenLayer.Cluster.Strategy to help the layer display a little quicker. My KML displays the attributes when clicked on before clustering was enabled, but now the bubble displays 'undefined, undefined'. How did the clustering change the extractAttributes property?                           var state = new OpenLayers.Layer.Vector("State Probation", {                 projection: map.displayProjection,                 strategies: [new OpenLayers.Strategy.Fixed(), new OpenLayers.Strategy.Cluster()],                 protocol: new OpenLayers.Protocol.HTTP({                     url: "state.kml",                     format: new OpenLayers.Format.KML({                         extractStyles: true,                         extractAttributes: true                     })                 })             });                       select = new OpenLayers.Control.SelectFeature([sample, state],                 {                     clickout: true, toggle: false,                     multiple: true, hover: false,                     toggleKey: "ctrlKey", // ctrl key removes from selection                     multipleKey: "shiftKey" // shift key adds to selection                 });                  sample.events.on({                 "featureselected": onFeatureSelect,                 "featureunselected": onFeatureUnselect             });                  state.events.on({                 "featureselected": onFeatureSelect,                 "featureunselected": onFeatureUnselect             });                  map.addControl(select);             select.activate();                     map.addControl(new OpenLayers.Control.LayerSwitcher());                        function onPopupClose(evt) {             select.unselectAll();         }              function onFeatureSelect(event) {             var feature = event.feature;             var selectedFeature = feature;             var popup = new OpenLayers.Popup.FramedCloud("chicken",                  feature.geometry.getBounds().getCenterLonLat(),                 new OpenLayers.Size(100,100),                 "<h2>"+feature.attributes.name + "</h2>" + feature.attributes.description,                 null, true, onPopupClose             );             feature.popup = popup;             map.addPopup(popup);         }         function onFeatureUnselect(event) {             var feature = event.feature;             if(feature.popup) {                 map.removePopup(feature.popup);                 feature.popup.destroy();                 delete feature.popup;             }         }


How do services like Google Maps compute the set of instructions for a route?

The problem of computing the driving directions from A to B can be divided into two equally important parts:   * Finding the best path from A to B.   * Identifying all the relevant maneuvers along that route. While there is loads of information and research papers on the first topic, no one seems to be talking about the second one. And from my experience, computing instructions is a non-trivial problem for at least the following reasons:   * Maps are very complex, diverse and sometimes full of geometry inconsistencies.   * It is very hard to define what a "relevant maneuver" is. Some users might want to hear a "Go straight" instruction at each crossing, just for reassurance. Others would find this totally annoying.   * It is hard to compare two sets of instructions and decide which one is better (because of the previous point). This makes automatic regression testing impossible. So, are there any resources on this topic? Some literature, best practices, or actually any kind of advice would be very welcome.


Can I reset "latlng" coordinates on Google Maps?

I have made a custom map type using the Google Maps API. I have used maptiler to create the different zoomlayers and the 0 layer (really zoomed out) is a 256 x 256 tile. The image only covers 134 x 112 however. I want to be able to make the map coordinates relative to the whole map and not the default earth coordinates of google maps. For example, the top-left corner should be (0,0) while the bottom-right should be (50,50). I have received help from Heitor Chang here: http://stackoverflow.com/questions/11192281/can-i-reset-latlng-coordinates- for-a-custom-map-type The problem lies here:               dayzProjection.prototype.fromLatLngToPoint = function(latlng) {     // pixel size of tile at zoom 0     var max = 134;          // define bottom-right corner as (50, 50)     var x = max * latlng.lng() / 50;     var y = max * latlng.lat() / 50;     return new google.maps.Point(x, y);     };          dayzProjection.prototype.fromPointToLatLng = function(pixel) {     // inverse conversion          var lng = pixel.x / max * 50;     var lat = pixel.y / max * 50;     return new google.maps.LatLng(lat, lng);     };      Can anyone help me transform these coordinates. This may help: https://developers.google.com/maps/documentation/javascript/maptypes#Projections


Searching Pixel Location (idrisi) on Another GIS

I'm Searching similar command Pixel Location (like idrisi) on Another GIS, I, need create raster file with values lat, long by separately... * * * Update: I need a raster file where the value of each cell are Latitude and another raster Longitude, (similar centroid of feature point) but on raster.. So this command on IDRISI Pixel location copy rows & col similar raster e.g srtmdem, and create 2 raster with a distribution for each pixel Latitude an Longitude.


Web map tool with high interactivity

I'm looking for an ArcGIS Online alternative. It needs to be   * linked to a private MySql server   * Use shapefiles   * Display attribute data and charts upon click/hover   * Toggle shapefile to display/disappear   * Privacy options   * User can drag one layer on top or below other layers. I like ArcGIS but it's very expensive. I like TileMill but haven't seen options where you can drag layers above or below.


ArcPy error: Add Feature Class to Terrain.

I'm writing a script to convert lidar data (in .las format) into a terrain without going through all the clicking and so I can send it out to neighboring states with similar programs (I'm using ArcGIS 10.1 btw). I've gotten almost there with the script but issues arise when I try adding the Multipoint to the Terrain in the Feature Dataset. I've done this step in model builder and exported that to a python script and am using the "in_features" that's output from the model. Incidentally, the model runs just fine in arcmap, but when I try to run the exported python script I get errors.               import os     import arcpy, math, glob     from arcpy import env          arcpy.CheckOutExtension("3D")     arcpy.env.overwriteOutput = True     #make workspace in memory     arcpy.env.workspace = "in_memory"          #create a feature dataset for export     #change geodatabase to a getparameter     geodatabase = r'C:\Users\peter.a.metzger\Desktop\Projects\LiDAR\test_data\lidar_test.gdb'     outgeodatabase = 'test_lidar_output'     projection = r'C:\Users\peter.a.metzger\AppData\Roaming\ESRI\Desktop10.1\ArcMap\Coordinate Systems\NAD83_UTM_zone_14N.prj'     created_FC_tool = arcpy.CreateFeatureDataset_management(geodatabase, outgeodatabase, projection)     created_FC = created_FC_tool.getOutput(0)     #Create a Terrain          terrain_name = 'test_Terrain' #Change to a getparameter     arcpy.CreateTerrain_3d (created_FC, terrain_name, 1.24605815948, 50000, '', 'WINDOWSIZE', 'ZMIN', 'NONE', 1)          output_terrain = os.path.join (created_FC, terrain_name)     arcpy.AddTerrainPyramidLevel_3d (output_terrain, 'WINDOWSIZE', '2.5 1000; 5 2000; 10 4000; 20 8000; 40 16000')     #process the .las file and add to the feature dataset and terrain          os.chdir (r'C:\Users\peter.a.metzger\Desktop\Projects\lidar_model_test_data\R06305328')     for LASfiles in glob.glob("*.las"):         print LASfiles         pointinfo = "point_file_information"         #Get Pointfile information         arcpy.PointFileInformation_3d(LASfiles, pointinfo, "LAS", "", projection)              #get average point spacing variable         pointspacing = arcpy.SearchCursor (pointinfo)         field = "Pt_Spacing"         for row in pointspacing:             avg_space = (row.getValue(field))         print avg_space              #convert las  to multipoint              #create lable for multipoint file         listlable = list(LASfiles)         del listlable[-4]         multipointlable = ''.join(listlable)         print multipointlable              print created_FC         outputPlable = os.path.join (created_FC, multipointlable)         #class_code, change to a getparameter         class_code = 2              print outputPlable              outputPlable_parameters = outputPlable + ' Shape Mass_Points 1 0 40 true false R06305328las_embed <None> false'         print outputPlable_parameters         arcpy.AddFeatureClassToTerrain_3d (output_terrain, outputPlable_parameters)      The error I'm getting is: > "Traceback (most recent call last): File > "C:\Users\peter.a.metzger\Desktop\Projects\LiDAR\las2terrain.py", line 64, > in arcpy.AddFeatureClassToTerrain_3d (output_terrain, > outputPlable_parameters) File "C:\Program Files > (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\ddd.py", line 1843, in > AddFeatureClassToTerrain raise e ExecuteError: Failed to execute. Parameters > are not valid. ERROR 000800: The value is not a member of . Failed to > execute (AddFeatureClassToTerrain)." And it's from the parameter that I set of: "outputPlable_parameters = outputPlable + ' Shape Mass_Points 1 0 40 true false R06305328las_embed false'" These values are what the model builder is using, yet it's not working.


How to add a new road to an existing pgRouting graph?

I have an existing pgrouting network topology that I would like to add a road to. The road is a linestring from a shapefile and will intersect several roads in the existing network. Is there an easy way to add the road into 'ways' so it is segmented properly and new vertices are created at the intersections?


How to reorder KML layers in ArcGIS API for JavaScript?

How do I reorder KML layers in Esri ArcGIS javascript version?   I want to be able to reorder the order of the rendering of layers. I tried `this._map.reorderLayer(kmlLayer, 0)` This did bring the layer's ID to the first in `this._map.layerIds` but the layer still stays on the top level I also tried `this._map.addLayer(kmlLayer, 0)` and this is not putting the newly added layer to the bottom level. Is it a bug or the methods are not doing what they sound like they should be doing??


Dikes extraction from DEM

I do have DEM of 1x1 meter and original LiDAR point cloud in .las, that DEM is made of. I need to extract river dikes (highest points of the dikes) to vector feature (point, polyline). Any idea for the algorythm or an existing tool? Dikes are in light grey on the 1st image and below is point cloud image for the example area with the dikes. River centerline in blue. ![Dikes are in light grey](http://i.stack.imgur.com/gV7Du.jpg) ![Point cloud with the dikes](http://i.stack.imgur.com/R1gj0.jpg) Any tool for linear features extraction from DEM? I belive I need something like ENVI’s Spatial Feature Extraction Module (page 7), but without ENVI's part, as I can not afford it :)


Recommended approach to integrate GeoNetwork with SQL Server database

I am new to GeoNetwork. But now my organization wants to invest in it to see how we can make the best use of it for publishing, sharing and distributing metadata. We currently have all our GIS data (mostly feature classes) stored in ArcSDE databases on SQL Server 2008 R2. Only brief information about the data has been captured, therefore no standardization has been pursued. My question is if we get a GeoNetwork catalogue set up, what is the best way to integrate it with the existing SQL Server databases. Also how does each metadata record in GeoNetwork get associated with the actual dataset in SQL Server? Can we set up GeoNetwork for it to automatically extract simple metadata of our specified from the SQL Server and push it to the GeoNetwork catalogue?


Is Terrain and Topographic Correction the same for Landsat8?

Landsat 8 L1 products are "Terrain Corrected", that said it means that are they already "topographically corrected", like both words are synonyms or you need to apply correction techniques as well? (eg i.topo.corr for Grass)


Appropriate pixel size when converting vector to raster

I am converting vector polygons into raster datasets. Is there a standard calculation for converting the scale of the vector into an appropriate pixel size (e.g. 25 metres)? My vector has been derived from Aerial Photo Interpretation (API) at a scale of 1:50000.


Census Data Intersection with Unique Boundaries

Basics: I have a country wide dataset of Census blocks, I need to find the population of unrelated boundaries that do not ALWAYS follow census blocks. I would like to use Select by Location-centroid as my method for determining what census blocks should be in the different boundaries. Any tips on how to do this in a model? * * * Clarification So I originally was intersecting my Boundaries (we will call them FBoundaries) with census Blocks to find the population of each of the fboundaries. However the fboundaries do not always follow the census blocks so I was getting double counts. So what I want to do is somehow tell the model to only take the census block if the centroid of the census block is in the Fboundary. Unfortunately I have the entire country I need to work with. And i'm lost –


How to extrude a single polygon shape to a interpolated surface elevation in arcscene?

I have a new problem to solve which is: how to extrude a single polygon shape to the curvature of a defined surface but not only to a specific one value height? I have several interpolated surfaces and I want to extrude a rectangular polygon shape inclined to those surfaces so that I can make a block diagram of several geological layers.:-)


How to proxy Geoserver through Apache in Windows

i am using apache web server on localhost:80 and geoserver served from jetty servlet on localhost:8888. i want to setup a proxy to redirect the requests such that requests to **localhost/geoserver** are routed to **localhost:8888/geoserver** i am following tutorial from http://gistutor.com/geoserver/21-intermediate- geoserver-tutorials/38-configuring-geoserver-proxy-for-public-and-remote-data- access.html and i am stuck in the part where it says               copy the text below and right click to paste the ProxyPass into your httpd.conf file      right above the </VirtualHost> tag.          ProxyRequests Off     ProxyPreserveHost On     <Proxy *>     Order deny,allow     Allow from all     </Proxy>     ProxyPass /geoserver http://localhost:8080/geoserver     ProxyPassReverse /geoserver http://localhost:8080/geoserver      i cannot find where to insert above text in httpd.conf file. the virtual hosts directive in httpd.conf file points to httpd-vhosts.conf file which looks like below               #     # Use name-based virtual hosting.     #     ##NameVirtualHost *:80          #     # VirtualHost example:     # Almost any Apache directive may go into a VirtualHost container.     # The first VirtualHost section is used for all requests that do not     # match a ServerName or ServerAlias in any <VirtualHost> block.     #     ##<VirtualHost *:80>     ##ServerAdmin postmaster@dummy-host.localhost     ##DocumentRoot "C:/xampp/htdocs/dummy-host.localhost"     ##ServerName dummy-host.localhost     ##ServerAlias www.dummy-host.localhost     ##ErrorLog "logs/dummy-host.localhost-error.log"     ##CustomLog "logs/dummy-host.localhost-access.log" combined     ##</VirtualHost>          ##<VirtualHost *:80>     ##ServerAdmin postmaster@dummy-host2.localhost     ##DocumentRoot "C:/xampp/htdocs/dummy-host2.localhost"     ##ServerName dummy-host2.localhost     ##ServerAlias www.dummy-host2.localhost     ##ErrorLog "logs/dummy-host2.localhost-error.log"     ##CustomLog "logs/dummy-host2.localhost-access.log" combined     ##</VirtualHost>      i have tried and changed this httpd-vhosts.conf file with above codes, but every time i do this i am unable to restart apache server. this is the httpd.conf file as in c:\xampp\apache\conf\               #     # ServerRoot: The top of the directory tree under which the server's     # configuration, error, and log files are kept.     #     # Do not add a slash at the end of the directory path.  If you point     # ServerRoot at a non-local disk, be sure to point the LockFile directive     # at a local disk.  If you wish to share the same ServerRoot for multiple     # httpd daemons, you will need to change at least LockFile and PidFile.     #     ServerRoot "C:/xampp/apache"          #     # Listen: Allows you to bind Apache to specific IP addresses and/or     # ports, instead of the default. See also the <VirtualHost>     # directive.     #     # Change this to Listen on specific IP addresses as shown below to      # prevent Apache from glomming onto all bound IP addresses.     #     #Listen 0.0.0.0:80     #Listen [::]:80     Listen 80          #     # Dynamic Shared Object (DSO) Support     #     # To be able to use the functionality of a module which was built as a DSO you     # have to place corresponding `LoadModule' lines at this location so the     # directives contained in it are actually available _before_ they are used.     # Statically compiled modules (those listed by `httpd -l') do not need     # to be loaded here.         #     # Example:       # LoadModule foo_module modules/mod_foo.so     #     LoadModule actions_module modules/mod_actions.so     LoadModule alias_module modules/mod_alias.so     LoadModule asis_module modules/mod_asis.so     LoadModule auth_basic_module modules/mod_auth_basic.so     LoadModule auth_digest_module modules/mod_auth_digest.so     #LoadModule authn_alias_module modules/mod_authn_alias.so     #LoadModule authn_anon_module modules/mod_authn_anon.so     #LoadModule authn_dbd_module modules/mod_authn_dbd.so     #LoadModule authn_dbm_module modules/mod_authn_dbm.so     LoadModule authn_default_module modules/mod_authn_default.so     LoadModule authn_file_module modules/mod_authn_file.so     #LoadModule authnz_ldap_module modules/mod_authnz_ldap.so     #LoadModule authz_dbm_module modules/mod_authz_dbm.so     LoadModule authz_default_module modules/mod_authz_default.so     LoadModule authz_groupfile_module modules/mod_authz_groupfile.so     LoadModule authz_host_module modules/mod_authz_host.so     #LoadModule authz_owner_module modules/mod_authz_owner.so     LoadModule authz_user_module modules/mod_authz_user.so     LoadModule autoindex_module modules/mod_autoindex.so     #LoadModule bucketeer_module modules/mod_bucketeer.so     #LoadModule cache_module modules/mod_cache.so     #LoadModule case_filter_module modules/mod_case_filter.so     #LoadModule case_filter_in_module modules/mod_case_filter_in.so     #LoadModule cern_meta_module modules/mod_cern_meta.so     LoadModule cgi_module modules/mod_cgi.so     #LoadModule charset_lite_module modules/mod_charset_lite.so     #LoadModule dav_module modules/mod_dav.so     #LoadModule dav_fs_module modules/mod_dav_fs.so     LoadModule dav_lock_module modules/mod_dav_lock.so     #LoadModule dbd_module modules/mod_dbd.so     LoadModule deflate_module modules/mod_deflate.so     LoadModule dir_module modules/mod_dir.so     #LoadModule disk_cache_module modules/mod_disk_cache.so     #LoadModule dumpio_module modules/mod_dumpio.so     #LoadModule echo_module modules/mod_echo.so     LoadModule env_module modules/mod_env.so     #LoadModule example_module modules/mod_example.so     #LoadModule expires_module modules/mod_expires.so     #LoadModule ext_filter_module modules/mod_ext_filter.so     #LoadModule file_cache_module modules/mod_file_cache.so     #LoadModule filter_module modules/mod_filter.so     LoadModule headers_module modules/mod_headers.so     #LoadModule ident_module modules/mod_ident.so     #LoadModule imagemap_module modules/mod_imagemap.so     LoadModule include_module modules/mod_include.so     LoadModule info_module modules/mod_info.so     LoadModule isapi_module modules/mod_isapi.so     #LoadModule ldap_module modules/mod_ldap.so     #LoadModule logio_module modules/mod_logio.so     LoadModule log_config_module modules/mod_log_config.so     #LoadModule log_forensic_module modules/mod_log_forensic.so     #LoadModule mem_cache_module modules/mod_mem_cache.so     LoadModule mime_module modules/mod_mime.so     #LoadModule mime_magic_module modules/mod_mime_magic.so     LoadModule negotiation_module modules/mod_negotiation.so     LoadModule proxy_module modules/mod_proxy.so     LoadModule proxy_ajp_module modules/mod_proxy_ajp.so     #LoadModule proxy_balancer_module modules/mod_proxy_balancer.so     LoadModule proxy_connect_module modules/mod_proxy_connect.so     LoadModule proxy_ftp_module modules/mod_proxy_ftp.so     LoadModule proxy_http_module modules/mod_proxy_http.so     LoadModule rewrite_module modules/mod_rewrite.so     LoadModule setenvif_module modules/mod_setenvif.so     #LoadModule speling_module modules/mod_speling.so     LoadModule ssl_module modules/mod_ssl.so     LoadModule status_module modules/mod_status.so     #LoadModule substitute_module modules/mod_substitute.so     #LoadModule unique_id_module modules/mod_unique_id.so     #LoadModule userdir_module modules/mod_userdir.so     #LoadModule usertrack_module modules/mod_usertrack.so     #LoadModule version_module modules/mod_version.so     #LoadModule vhost_alias_module modules/mod_vhost_alias.so          <IfModule !mpm_netware_module>     <IfModule !mpm_winnt_module>     #     # If you wish httpd to run as a different user or group, you must run     # httpd as root initially and it will switch.       #     # User/Group: The name (or #number) of the user/group to run httpd as.     # It is usually good practice to create a dedicated user and group for     # running httpd, as with most system services.     #     User daemon     Group daemon          </IfModule>     </IfModule>          # 'Main' server configuration     #     # The directives in this section set up the values used by the 'main'     # server, which responds to any requests that aren't handled by a     # <VirtualHost> definition.  These values also provide defaults for     # any <VirtualHost> containers you may define later in the file.      #     # All of these directives may appear inside <VirtualHost> containers,     # in which case these default settings will be overridden for the     # virtual host being defined.     #          #     # ServerAdmin: Your address, where problems with the server should be     # e-mailed.  This address appears on some server-generated pages, such     # as error documents.  e.g. admin@your-domain.com     #     ServerAdmin postmaster@localhost          #     # ServerName gives the name and port that the server uses to identify itself.This can often be determined automatically, but we recommend you specify it explicitly to   prevent problems during startup.     # If your host doesn't have a registered DNS name, enter its IP address here.     #     ServerName localhost:80          #     # DocumentRoot: The directory out of which you will serve your     # documents. By default, all requests are taken from this directory, but     # symbolic links and aliases may be used to point to other locations.     #     DocumentRoot "C:/xampp/htdocs"          #     # Each directory to which Apache has access can be configured with respect     # to which services and features are allowed and/or disabled in that     # directory (and its subdirectories).      #     # First, we configure the "default" to be a very restrictive set of      # features.       #     <Directory />     Options FollowSymLinks     AllowOverride None     Order deny,allow     Deny from all     </Directory>          #     # Note that from this point forward you must specifically allow     # particular features to be enabled - so if something's not working as     # you might expect, make sure that you have specifically enabled it     # below.     #          #     # This should be changed to whatever you set DocumentRoot to.     #     <Directory "C:/xampp/htdocs">     #     # Possible values for the Options directive are "None", "All",     # or any combination of:     #   Indexes Includes FollowSymLinks SymLinksifOwnerMatch ExecCGI MultiViews     #     # Note that "MultiViews" must be named *explicitly* --- "Options All"     # doesn't give it to you.     #     # The Options directive is both complicated and important.  Please see     # http://httpd.apache.org/docs/2.2/mod/core.html#options     # for more information.     #     Options Indexes FollowSymLinks Includes ExecCGI          #     # AllowOverride controls what directives may be placed in .htaccess files.     # It can be "All", "None", or any combination of the keywords:     #   Options FileInfo AuthConfig Limit     #     AllowOverride All          #     # Controls who can get stuff from this server.     #     Order allow,deny     Allow from all          </Directory>          #     # DirectoryIndex: sets the file that Apache will serve if a directory     # is requested.     #     <IfModule dir_module>     DirectoryIndex index.php index.pl index.cgi index.asp index.shtml index.html index.htm \                    default.php default.pl default.cgi default.asp default.shtml default.html default.htm \                    home.php home.pl home.cgi home.asp home.shtml home.html home.htm     </IfModule>          #     # The following lines prevent .htaccess and .htpasswd files from being      # viewed by Web clients.      #     <FilesMatch "^\.ht">     Order allow,deny     Deny from all     Satisfy All     </FilesMatch>          #     # ErrorLog: The location of the error log file.     # If you do not specify an ErrorLog directive within a <VirtualHost>     # container, error messages relating to that virtual host will be     # logged here.  If you *do* define an error logfile for a <VirtualHost>     # container, that host's errors will be logged there and not here.     #     ErrorLog "logs/error.log"          #ScriptLog "logs/cgi.log"          #     # LogLevel: Control the number of messages logged to the error_log.     # Possible values include: debug, info, notice, warn, error, crit,     # alert, emerg.     #     LogLevel warn          <IfModule log_config_module>     #     # The following directives define some format nicknames for use with     # a CustomLog directive (see below).     #     LogFormat "%h %l %u %t \"%r\" %>s %b \"%{Referer}i\" \"%{User-Agent}i\"" combined     LogFormat "%h %l %u %t \"%r\" %>s %b" common          <IfModule logio_module>       # You need to enable mod_logio.c to use %I and %O       LogFormat "%h %l %u %t \"%r\" %>s %b \"%{Referer}i\" \"%{User-Agent}i\" %I %O" combinedio     </IfModule>          #     # The location and format of the access logfile (Common Logfile Format).     # If you do not define any access logfiles within a <VirtualHost>     # container, they will be logged here.  Contrariwise, if you *do*     # define per-<VirtualHost> access logfiles, transactions will be     # logged therein and *not* in this file.     #     #CustomLog "logs/access.log" common          #     # If you prefer a logfile with access, agent, and referer information     # (Combined Logfile Format) you can use the following directive.     #     CustomLog "logs/access.log" combined     </IfModule>          <IfModule alias_module>     #     # Redirect: Allows you to tell clients about documents that used to      # exist in your server's namespace, but do not anymore. The client      # will make a new request for the document at its new location.     # Example:     # Redirect permanent /foo http://localhost/bar          #     # Alias: Maps web paths into filesystem paths and is used to     # access content that does not live under the DocumentRoot.     # Example:     # Alias /webpath /full/filesystem/path     #     # If you include a trailing / on /webpath then the server will     # require it to be present in the URL.  You will also likely     # need to provide a <Directory> section to allow access to     # the filesystem path.          #     # ScriptAlias: This controls which directories contain server scripts.      # ScriptAliases are essentially the same as Aliases, except that     # documents in the target directory are treated as applications and     # run by the server when requested rather than as documents sent to the     # client.  The same rules about trailing "/" apply to ScriptAlias     # directives as to Alias.     #     ScriptAlias /cgi-bin/ "C:/xampp/cgi-bin/"          </IfModule>          <IfModule cgid_module>     #     # ScriptSock: On threaded servers, designate the path to the UNIX     # socket used to communicate with the CGI daemon of mod_cgid.     #     #Scriptsock "logs/cgi.sock"     </IfModule>          #     # "C:/xampp/cgi-bin" should be changed to whatever your ScriptAliased     # CGI directory exists, if you have that configured.     #     <Directory "C:/xampp/cgi-bin">         AllowOverride None         Options None         Order allow,deny         Allow from all     </Directory>          #     # DefaultType: the default MIME type the server will use for a document     # if it cannot otherwise determine one, such as from filename extensions.     # If your server contains mostly text or HTML documents, "text/plain" is     # a good value.  If most of your content is binary, such as applications     # or images, you may want to use "application/octet-stream" instead to     # keep browsers from trying to display binary files as though they are     # text.     #     DefaultType text/plain          <IfModule mime_module>         #         # TypesConfig points to the file containing the list of mappings from         # filename extension to MIME-type.         #         TypesConfig "conf/mime.types"              #         # AddType allows you to add to or override the MIME configuration         # file specified in TypesConfig for specific file types.         #         #AddType application/x-gzip .tgz         #         # AddEncoding allows you to have certain browsers uncompress         # information on the fly. Note: Not all browsers support this.         #         #AddEncoding x-compress .Z         #AddEncoding x-gzip .gz .tgz         #         # If the AddEncoding directives above are commented-out, then you         # probably should define those extensions to indicate media types:         #         AddType application/x-compress .Z         AddType application/x-gzip .gz .tgz              #         # AddHandler allows you to map certain file extensions to "handlers":         # actions unrelated to filetype. These can be either built into the server         # or added with the Action directive (see below)         #         # To use CGI scripts outside of ScriptAliased directories:         # (You will also need to add "ExecCGI" to the "Options" directive.)         #         AddHandler cgi-script .cgi .pl .asp .py              # For type maps (negotiated resources):         #AddHandler type-map var              #         # Filters allow you to process content before it is sent to the client.         #         # To parse .shtml files for server-side includes (SSI):         # (You will also need to add "Includes" to the "Options" directive.)         #         AddType text/html .shtml         AddOutputFilter INCLUDES .shtml     </IfModule>          <IfModule mime_magic_module>         #         # The mod_mime_magic module allows the server to use various hints from the         # contents of the file itself to determine its type.  The MIMEMagicFile         # directive tells the module where the hint definitions are located.         #         MIMEMagicFile "conf/magic"     </IfModule>          #     # Customizable error responses come in three flavors:     # 1) plain text 2) local redirects 3) external redirects     #     # Some examples:     #ErrorDocument 500 "The server made a boo boo."     #ErrorDocument 404 /missing.html     #ErrorDocument 404 "/cgi-bin/missing_handler.pl"     #ErrorDocument 402 http://localhost/subscription_info.html     #          #     # EnableMMAP and EnableSendfile: On systems that support it,       # memory-mapping or the sendfile syscall is used to deliver     # files.  This usually improves server performance, but must     # be turned off when serving from networked-mounted      # filesystems or if support for these functions is otherwise     # broken on your system.     #     #EnableMMAP off     #EnableSendfile off          # Supplemental configuration     #     # The configuration files in the conf/extra/ directory can be      # included to add extra features or to modify the default configuration of      # the server, or you may simply copy their contents here and change as      # necessary.          # XAMPP specific settings     Include "conf/extra/httpd-xampp.conf"          # Perl settings     Include "conf/extra/httpd-perl.conf"          # Server-pool management (MPM specific)     Include "conf/extra/httpd-mpm.conf"          # Multi-language error messages     Include "conf/extra/httpd-multilang-errordoc.conf"          # Fancy directory listings     Include "conf/extra/httpd-autoindex.conf"          # Language settings     Include "conf/extra/httpd-languages.conf"          # User home directories     Include "conf/extra/httpd-userdir.conf"          # Real-time info on requests and configuration     Include "conf/extra/httpd-info.conf"          # Virtual hosts      Include "conf/extra/httpd-vhosts.conf"          # Distributed authoring and versioning (WebDAV)     # Attention! WEB_DAV is a security risk without a new userspecific configuration for a      secure authentifcation      # Include "conf/extra/httpd-dav.conf"          # Implements a proxy/gateway for Apache.     Include "conf/extra/httpd-proxy.conf"          # Various default settings     Include "conf/extra/httpd-default.conf"          # Secure (SSL/TLS) connections     Include "conf/extra/httpd-ssl.conf"       #     # Note: The following must must be present to support     #       starting without SSL on platforms with no /dev/random equivalent     #       but a statically compiled-in mod_ssl.     #     <IfModule ssl_module>         SSLRandomSeed startup builtin         SSLRandomSeed connect builtin         SSLSessionCache "shmcb:logs/ssl.scache(512000)"         SSLSessionCacheTimeout 300     </IfModule>          # AJP13 Proxy     <IfModule mod_proxy.c>     <IfModule mod_proxy_ajp.c>     Include "conf/extra/httpd-ajp.conf"     </IfModule>     </IfModule>


WorldWind wms services show reversed coordinate in ArcGIS Desktop

I'm adding WorldWind wms services (http://data.worldwind.arc.nasa.gov/wms?) to ArcMap (ArcGIS 10.2.2) but maps are misplaced, looks like they have reversed coordinates (as explained here: http://viswaug.wordpress.com/2009/03/15/reversed-co-ordinate-axis-order-for- epsg4326-vs-crs84-when-requesting-wms-130-images/). Is there any workaround to get maps properly placed?


Snapping with polgon drawing

I am drawing a polygon with the drawing function within ArcGIS 10.1. Is there a way to snap orthographically at 90 degrees while using the polygon draw command.![enter image description here](http://i.stack.imgur.com/oox3v.jpg) I added a picture of what i would like to do but as you can see the callout is not drawn at 90 degrees.


QGis road graph plugin can't be configured

I am using QGIS 2.2.0 Valimera on a Windows 7 ultimate 64 bit version. I want to calculate the shortest path on a line layer (OSM street layer) with the road graph plugin. The settings option to calibrate the plugin doesn't work as the layer option and following features are grayed out and can't be selected. CRS of the street layer and the project layer are matching and on the fly projection is enabled. Could somebody help me please with this problem? It already appeared for other users but I didnt find appropriate answers.


Removing highlighting when selecting feature in Arcmap

Whenever I select a feature in an attribute table it gets highlighted which prevents me from seeing the symbology of that feature. I tried choosing No Color in the appearance option, but that results in it being highlighted white. This is what it looks like after a select by attribute query (how I'd like it) ![](http://i.imgur.com/Vv1rF3w.png) But this is how it appears when I double click the feature to go to it ![](http://i.imgur.com/ChH7FyA.png) Anyone know a way to make it hollow or adjust the transparency?


ArcGIS and Python Issues with Different Versions

As we know, there are few versions of Arcgis (9.x,10) and Python(2.x,3.x) are in use. Issues are facing when we run scripts/toolbox of arcgis since the script written in one version of python and the python installed with arcgis is different version. Q1) How to identify python script version so that we can directly execute and save time? Q2) How to convert the existing python script to the installed python version on the system and get result without hassle of programing? Q3) Is there any procedure for arcgis to take python version 2/3.x as default so that, we can install multiple versions of python and change default version to run according to our requirement?


Why is getting a GPS fix so slow?

I understand the basic concepts of GPS positioning system and the use of trilateration but I don't understand why it takes so much time to get a fix. On commercial GPS like TomTom it usually takes about one minute and on my phone it can take more than 5 minutes to be located. Is it because of the computation time or it needs to collect more data before having a good approximation? Thank you


Repairing disconnected layers with selection sets

Has anyone had any success of repairing disconnected layers that have selection sets that have their data sources broken in ArcGIS? Is there a way to detect this through code and then re-create that selection set? I tried:               IFeatureLayerDefinition.DefinitionSelectionSet.Count      and               IFeatureSelection.SelectionSet.Count      of which just give an "Object reference not set to an instance of an object" error. We are on version 10.0 and code in VB.NET or C#. Thanks Rob


ESRI UK Productivity Suite: Critical Error - Help Resolving

I am using the OS Conversion tool under ProductivitySuite v2.1. to convert OS MasterMap gml to a feature class. Having carefully followed the instructions I am getting the following: Error 18/05/2012 09:56:57 CRITICAL ERROR: Unable to validate OrderNumber This is the first time I have used the tool and unsure as to How would I resolve this? Apparently the data has been copied directly from the CD as it was provided by the OS. Log File below, if it helps: 18/05/2012 09:56:57 Log File created. Info 18/05/2012 09:56:57 ================================================= Info 18/05/2012 09:56:57 OS Data Converter Info 18/05/2012 09:56:57 Version: 2.1 Build: 557 Info 18/05/2012 09:56:57 Info 18/05/2012 09:56:57 Start Processing Info 18/05/2012 09:56:57 Profile name: OS MM SWW Info 18/05/2012 09:56:57 Start condition: New Info 18/05/2012 09:56:57 Info 18/05/2012 09:56:57 ================================================= Info 18/05/2012 09:56:57 Validation: Info 18/05/2012 09:56:57 - Validating data order number Info 18/05/2012 09:56:57 **Error 18/05/2012 09:56:57 CRITICAL ERROR: Unable to validate OrderNumber** Error 18/05/2012 09:56:57 CRITICAL ERROR DETAILS: Error 18/05/2012 09:56:57 ESRIUK.Products.UKPS.DC.Core.DataConversionException: Unable to validate OrderNumber ---> System.Exception: Unable to validate order number Error 18/05/2012 09:56:57 at ESRIUK.Products.UKPS.DC.Profile.OrderNumber.Validate(List`1 sourceItems) Error 18/05/2012 09:56:57 at ESRIUK.Products.UKPS.DC.Profile.OrderNumber..ctor(List`1 sourceItems) Error 18/05/2012 09:56:57 at ESRIUK.Products.UKPS.DC.Profile.DataConverterSourceList.ValidateOrderNumber() Error 18/05/2012 09:56:57 at ESRIUK.Products.UKPS.DC.Core.DC.DataConverterProcess.ValidateOrderNumber() Error 18/05/2012 09:56:57 --- End of inner exception stack trace --- Error 18/05/2012 09:56:57 at ESRIUK.Products.UKPS.DC.Core.DC.DataConverterProcess.ValidateOrderNumber() Error 18/05/2012 09:56:57 at ESRIUK.Products.UKPS.DC.Core.DC.DataConverterProcess.DoValidate(String paramName) Error 18/05/2012 09:56:57 at ESRIUK.Products.UKPS.DC.Core.DataConverterEngine.ValidateProfile() Error 18/05/2012 09:56:57 at ESRIUK.Products.UKPS.DC.Core.BaseEngine`3.ProcessProfile(T_Profile profile) Info 18/05/2012 09:56:57 Info 18/05/2012 09:56:57 OS Data Converter many thanks....


How to upload and download to gps

Do you have a step by step guide on how to upload/download points/wayponts and polylines/routes from a gps? Regards Christo


Problem with arcpy.CalculateField_management

I have a table with a column named "Kategorie", filled with values. I added a new column named "RULE". I want to fill this column with Strings depending on "KATEGORIE". I have a PY-Script and it doesn't work. What is wrong. Thank you for your help ! Expression:               "getKat(!KATEGORIE!)"      Codeblock:               """def getKat(KATEGORIE):         if KATEGORIE = "0":             myval = "Autobahn"         if KATEGORIE = "1":             myval = "Bundesstrasse"         else:             myval= 3"""     arcpy.CalculateField_management(ERG, "RULE", expression, "", codeblock)      * * * Thank you for your answers ! But both don't work. I think the problem is the line "expression" or "codeblock". I try to print the result of "myval".The result is a error message "name myval is not defined" expression:               "getKat(!KATEGORIE!)"      codeblock:               """def getKat(KATEGORIE):         if KATEGORIE == "0":             myval == "Autobahn"         elif KATEGORIE == "1":             myval == "Bundesstrasse"         else:             myval == "3"         return myval"""     print "myval", myval     arcpy.CalculateField_management(ERG, "RULE", expression, "PYTHON", codeblock)


Raster Calculator Bug at start

I use ArcGIS 10.1 on a Windows XP. ArcGIS completely shuts down when I attempt to launch the "Raster Calculator" tool. Has this happened to someone else? Any ideas on how to solve this annoying issue? Otherwise, can I calculate the sum of raster values using command-line GDal?


Python ArcGIS 10.1 Search Cursor

I have a file geodatabase with a polygon feature class named "AOI". The "AOI" feature class currently holds two records. I want to loop through each record and run a geoprocessing task, for example the ExtractByMask tool. When I run my script it only runs against the first records and stops successfully. Below is my code....               with arcpy.da.SearchCursor(AOI, ("SHAPE@")) as cursor:         for row in cursor:             Extract = arcpy.sa.ExtractByMask(costraster, row)


How to use gs:GetFullCoverage from Geoserver?

I'm asking this question because apparently the process needs a name as input, but I don't know what I should put.


Smoothing Lines in PostGIS

I want to smooth Lines in postgis, to be more specific, I want to smooth road data from OpenStreetMap. Basicly I thought of something like this scenario. ![This would be a possible road segment.](http://i.stack.imgur.com/gEQpE.png) This would be a possible road segment. ![enter image description here](http://i.stack.imgur.com/1kLAV.png) I want to use ST_LineInterpolate_Point() to select the middle Points ![enter image description here](http://i.stack.imgur.com/L4fTg.png) After this, I want to create New Points between the original ones and the darkgreen helpline. (I have no Idea how to achieve this...) ![enter image description here](http://i.stack.imgur.com/EyPJO.png) I read about circularlinestrings. I hope it will be possible to lay one of these over the darkgreen and red points. ![enter image description here](http://i.stack.imgur.com/I1zkQ.png) Finaly I want to use ST_LineInterpolate_Point() again on the circularlinestring (if this is possible) and create some new points(yellow ones), to make a new normal LineString. To sum up, I am new to PostGIS, but since I haven't found anything helpful for LineSmoothing I need a function where I can input a LineString and the number of yellow points per segment. The function should return a new LineString. First of all: Is this possible to achieve with just pl/pgsql? Second, since I am new to all of this, can somebody give me some advice where to start?


Where can I find GeoJSON resources files?

I'm looking for demographic information and various things, and most of the time data are available as CSV (my application use geoJSON format). Shall I encode the CSV to geoJSON or is there a "library" for geoJSON ? Some kind of a database of geojson files ? Any help is welcome :) Cheers


Symbology issue with feature class limited by domain

I added a field to a feature class in order to maintain symbology from the original KML file (this is a new question following up on this one). The new field is limited to three options by a **domain**. When I go to symbolize based on this field, all values are considered to be as opposed to being separated into their appropriate category. Does this have to do with being controled by a coded domain or is something else going on here? ![Undesired Symbology Result](http://i.stack.imgur.com/76ZkM.png)


Importing data from excel and maintaining the time format

I am trying to import a file from Excel which is a GPS data set of consecutive points at 15minute intervals. When I import the file and save it as a shape file the time column is not read properly and GIS converts to 00:00:00 and then further down the column shifts to 01/01/1900. I have tried to format the column in Excel as either hh:mm, hh:mm:ss, dd/mm/yyyy hh:mm or no number format and each has the same end result. The imported table shows the correct times but as soon as the data table is converted to a shape file to be displayed the time format is lost. I have never noticeably had this problem when importing data previously and I did not need to format the time. Now it seems that no formatting method will work properly. How can I format the time column to retain the time values when imported and saved as a shape file? Thanks for your help.


Fully automating a workspace to create maps with different file names but same information?

I have a lot of maps (over 600) to create. They all contain the same column titles and the layouts will be the same. I don't want to have to go through and repeat the same process everytime, I would like to just have to open the file and have it map on its own. I have created a workspace to open the other tables I need for every map, but I need to know if it is possible to write a program or code with mapbasic to ask the user to specify what table to open, then using that finish the program by joining the tables and mapping. The first and last parts I can do. I just need to know how to create a dialog box to have the user select the table, or if it is even possible.


How to install GDAL with Python on windows?

Can anyone please explain to me in step by step instructions how to install GDAL/OGR with Python on Windows. I have Windows Vista and I have tried following the information on the website and it does not seem to bind the exe files for me. So, if someone would be so kind to walk me through the process, including links to the files/folders I will need I would be happy to do a screen cast to continue to share this with the rest of the world. thanks in advance \--20100930 update-- I have now tried to run the gdal setup with minGW, but this has also failed: ![GDAL_FAIL](http://i.stack.imgur.com/pPOH1.jpg)


Installation of GDAL 1.10 in Windows 8 and Ubuntu 12.04

Is there an easy way to install GDAL/OGR 1.10 in Windows 8 and Ubuntu 12.04.


Import Error while Importing osgeo package

> **Possible Duplicate:**   >  How to install GDAL with Python on windows? I have downloaded Fwtools for python2.5 from http://fwtools.maptools.org(windows installer) after installing it, it was giving error as Unable to import osgeo. Then I checked the whole FWtools for osgeo but nothing was there .. can anyone guide me how can I get "osgeo" and "gdal/ogr" package for the Python (as I don't have any python documentation for osgeo and the command pydoc -g osgeo also doesn't work) Presently I am using all these packages from inside the Qgis Console..Can't I have the same package for my python bundle..