title: SVG trouble in QGIS

text: I am having some issues using some SVGs that I have created in Inkscape. I am trying to use them in QGIS but they're not displayed properly. This: ![http://i.imgur.com/BuxMWeO.png](http://i.stack.imgur.com/xW8yZ.png) is how the SVG looks like in Inkscape and how it should look on the map. This:![http://i.imgur.com/qoIFsJm.png](http://i.stack.imgur.com/bXQYX.png) is how QGIS "sees" it. I am having similar issues with other SVGs. I would also like to mention that this issue only occurs on SVGs that contain text. To be more precise, the text is the element that is not displayed properly. There was this very simple SVG; a red, thick circle with a "V" inside of it. In Inkscape, it looked good. In QGIS, the "V" would not show up. What am I doing wrong?


title: Can't add qgis repository for openSUSE

text: Linux noob here. I keep getting "Unable to create repository" when trying to add the repo url on the qgis download page to YaST: `http://download.opensuse.org/repositories/Application:/Geo/<version>/` I've tried replacing `<version>` with a few things with no luck:   * 1.8   * 1.8.0   * 12.3   * qgis-1.8   * qgis-1.8.0 My search has not come up with any similar cases so far. Anyone know a working repo for installing qgis on openSUSE 12.3?


title: QGIS for Mac: How to restore menu bar and tool bar in print composer

text: I managed to hide both the tool bar and the menu bar in my print composer (Mac version 1.8.0). To get the tool bar back I saw you need to click on the menu bar. What then when you also miss the menu bar? ![enter image description here](http://i.stack.imgur.com/YhZqj.png) Here what I see. I right-clicked everywhere. I can't get the the check-box has described here.


title: Incorporate Google Streetview into a Leaflet map

text: The person who asked the question at the link below way back in 2011 said, "I understand that it is possible to incorporate Google StreetView into non- Google apps". Can anyone tell me how this is possible? I'd like to add street view functionality to my leaflet map. Open source alternatives to Google's StreetView engine?


title: Open Source + Offline Geo Notebook

text: I am looking for an open source offline software to record points of interest that I might visit at a later date. Let's say I am reading a book about Alaska, I'd like to take notes about interesting places, probably with the help of google maps to locate them on the map. Later I'd like to review the points, select the best ones and plan my vacation. Important are a simple and quick input of positions, for example cut and paste urls from google maps, and the possibility to add annotations. Trip planning and export to gps would be great, but if it can work with another program for planning, I'd be fine with that. I'd like to keep it offline and open source as it may take a few years till I get to the marked points, and I don't want to rely on / trust an online service either.


title: change the area in a osm-request: from lon & lat to whole country

text: I am new to PHP's SimpleXML. the script below makes a starting point for my project. Below I have added my own answer to the question, how to refine the code to ad more tags. I want to filter the data to get the nodes with special category. Here is sample of the OSM data I want to get the whole schools within an certain area. The first script below runs well - but now I want to refine the search and add more tags. I want to store all into MySQL.So we need to make some XML parsing with PHP: The following is a little OSM Overpass API example with PHP SimpleXML NOTE: i want to change the values in the boundary-box. So that we can gather the shools from a certain area - eg. the country Italy Q: is this possible? see the script:               <?php     /**      * OSM Overpass API with PHP SimpleXML / XPath      *      * PHP Version: 5.4 - Can be back-ported to 5.3 by using 5.3 Array-Syntax (not PHP 5.4's square brackets)      */          //     // 1.) Query an OSM Overpass API Endpoint     //               $query = 'node       ["addr:postcode"~"RM12"]       (51.5557914,0.2118915,51.5673083,0.2369398);       node       (around:100000)       ["amenity"~"school"];                out;';          $context = stream_context_create(['http' => [         'method'  => 'POST',         'header' => ['Content-Type: application/x-www-form-urlencoded'],         'content' => 'data=' . urlencode($query),     ]]);          # please do not stress this service, this example is for demonstration purposes only.     $endpoint = 'http://overpass-api.de/api/interpreter';     libxml_set_streams_context($context);     $start = microtime(true);          $result = simplexml_load_file($endpoint);     printf("Query returned %2\$d node(s) and took %1\$.5f seconds.\n\n", microtime(true) - $start, count($result->node));          //     // 2.) Work with the XML Result     //          # get all school nodes with xpath     $xpath = '//node[tag[@k = "amenity" and @v = "school"]]';     $schools = $result->xpath($xpath);     printf("%d School(s) found:\n", count($schools));     foreach ($schools as $index => $school)     {         # Get the name of the school (if any), again with xpath         list($name) = $school->xpath('tag[@k = "name"]/@v') + ['(unnamed)'];         printf("#%02d: ID:%' -10s  [%s,%s]  %s\n", $index, $school['id'], $school['lat'], $school['lon'], $name);     }          ?>      if we change this               $query = 'node       ["addr:postcode"~"RM12"]       (51.5557914,0.2118915,51.5673083,0.2369398);       node       (around:100000)       ["amenity"~"school"];                out;';          $context = stream_context_create(['http' => [         'method'  => 'POST',         'header' => ['Content-Type: application/x-www-form-urlencoded'],         'content' => 'data=' . urlencode($query),     ]]);      we get more than 2170 results:               #2171: ID:2837001831  [51.1783749,-0.7970866]  More House School     #2172: ID:2865091022  [50.9090614,-0.5565425]  Dorset House School     #2173: ID:2882477853  [51.6261198,-0.7349665]  Bowerdean Primary School     #2174: ID:2901434856  [51.6542477,-0.3098923]  The Haberdashers' Aske's Girls School     #2175: ID:2901434857  [51.6565707,-0.3129822]  The Haberdashers' Aske's Boys School      NOTE: i want to change the values in the boundary-box. So that we can gather the shools from a certain area - eg. the country Italy or something thelike is this possible? **update:** can i do some rework of the area-data like so:               [place] => SimpleXMLElement Object             (                 [@attributes] => Array                     (                         [place_id] => 331526                         [osm_type] => node                         [osm_id] => 96209423                         [place_rank] => 15                         [boundingbox] => 44.4361381530762,44.4361419677734,26.1027431488037,26.1027450561523                         [lat] => 44.436139                         [lon] => 26.1027436                         [display_name] => BucureČti, Sector 2, Bucuresti, RomĂ˘nia, European Union                         [class] => place                         [type] => city                         [importance] => 0.73231672860554                         [icon] => http://nominatim.openstreetmap.org/images/mapicons/poi_place_city.p.20.png                     )                  )      is this doable!`?


title: Draw area from street borders

text: I have searched online for almost a week and can not find an answer to my issue. I need to generate polygons from street borders in a given city. I have downloaded data from openstreetmaps and imported it into QGIS. Is there a way to generate polygons and use streets as borders? Thanks in advance.


title: Displaying a Selected Location on ArcGIS Map View

text: I am New to ArcGIS. I am Developing an iOS application Using ArcGIS iOS SDK. In that application i want to implement general search of an address like if I type a key word in my search field i want to display all the addresses,blocks,Highways,Avenues etc.. matching to that key word. I write the below code for implementing general search. ViewController.m               - (void)viewDidLoad     {         [super viewDidLoad];         self.mapView.layerDelegate = self;         self.mapView.callout.delegate = self;         self.searchBar.delegate = self;         AGSCredential* cred = [[AGSCredential alloc] initWithUser:@"username" password:@"password"];              AGSTiledMapServiceLayer *tiledMapServiceLayer = [AGSTiledMapServiceLayer tiledMapServiceLayerWithURL:[NSURL URLWithString:@"url"] credential:cred];         [self.mapView addMapLayer:tiledMapServiceLayer];              self.dynamicLayer = [[AGSDynamicMapServiceLayer alloc] initWithURL:[NSURL URLWithString:@"Dynamic Layer Ursl"] credential:cred];         [self.mapView addMapLayer:self.dynamicLayer withName:@"POIS"];         [self.mapView addMapLayer:self.dynamicLayer withName:@"ADDRESSES"];         [self.mapView addMapLayer:self.dynamicLayer withName:@"HIGHWAYS"];         [self.mapView addMapLayer:self.dynamicLayer withName:@"AVENUES"];         [self.mapView addMapLayer:self.dynamicLayer withName:@"ROADS_LANES"];         [self.mapView addMapLayer:self.dynamicLayer withName:@"BLOCKS"];         [self.mapView addMapLayer:self.dynamicLayer withName:@"AREAS"];         [self.mapView addMapLayer:self.dynamicLayer withName:@"GOVERNORATE"];              self.graphicsLayer = [AGSGraphicsLayer graphicsLayer];         [self.mapView addMapLayer:self.graphicsLayer withName:@"Results"];         self.graphicsLayer.calloutDelegate = self;              self.findTask = [[AGSFindTask alloc] initWithURL:[NSURL URLWithString:@"findTaskURL"] credential:cred];         self.findTask.delegate = self;              self.findParams = [[AGSFindParameters alloc]init];          }          - (void)searchBarSearchButtonClicked:(UISearchBar *)searchBar {          self.mapView.callout.hidden = YES;     self.findParams.contains = YES;     self.findParams.layerIds = [NSArray arrayWithObjects:@"0",@"1",@"2",@"3",@"4",@"5",@"6",@"7",nil];     self.findParams.outSpatialReference = self.mapView.spatialReference;     self.findParams.returnGeometry = TRUE;     self.findParams.searchFields = [NSArray arrayWithObjects:@"SearchFields",nil];     self.findParams.searchText = searchBar.text;     [self.findTask executeWithParameters:self.findParams];     [searchBar resignFirstResponder];     }      -(void)findTask:(AGSFindTask *)findTask operation:(NSOperation *)op didExecuteWithFindResults:(NSArray *)results     {          [self.graphicsLayer removeAllGraphics];          poiArray = [[NSMutableArray alloc]init];     addressArray = [[NSMutableArray alloc]init];     highwayArray = [[NSMutableArray alloc]init];     avenueArray = [[NSMutableArray alloc]init];     roadLanesArray = [[NSMutableArray alloc]init];     blocksArray = [[NSMutableArray alloc]init];     areasArray = [[NSMutableArray alloc]init];     governarateArray = [[NSMutableArray alloc]init];          poiObjectArray = [[NSMutableArray alloc]init];     addressObjectArray = [[NSMutableArray alloc]init];     highwayObjectArray = [[NSMutableArray alloc]init];     avenueObjectArray = [[NSMutableArray alloc]init];     roadLanesObjectArray = [[NSMutableArray alloc]init];     blocksObjectArray = [[NSMutableArray alloc]init];     areasObjectArray = [[NSMutableArray alloc]init];     governarateObjectArray = [[NSMutableArray alloc]init];          AGSFindResult *result = nil;          NSMutableDictionary *d = [[NSMutableDictionary alloc]init];     NSMutableString *d1 = [[NSMutableString alloc]init];     NSMutableString *d2 = [[NSMutableString alloc]init];     NSMutableString *d3 = [[NSMutableString alloc]init];          for (int i=0;i<[results count];i++)     {         result = [results objectAtIndex:i];         NSDictionary *json = [result encodeToJSON];              NSInteger i = [[json valueForKey:@"layerId"]integerValue];         NSLog(@"%ld",(long)i);         d = [json valueForKey:@"attributes"];              if(i == 0)         {             d1 = [d valueForKey:@"NAME"];             [poiArray addObject:d1];             [poiObjectArray addObject:result];             d1 = nil;         }         else if(i==1)         {             d1 = [d valueForKey:@"BUILDING NO"];             d2 = [d valueForKey:@"BLOCK NO"];             d3 = [d valueForKey:@"ROAD NO"];             NSMutableString *str = [NSMutableString stringWithFormat:@"Bul-%@,Blo-%@,Road-%@",d1,d2,d3];             [addressArray addObject:str];             [addressObjectArray addObject:result];             d1=nil;             d2=nil;             d3=nil;         }         else if(i==2)         {             d1 = [d valueForKey:@"ROAD NAME ENGLISH"];             d2 = [d valueForKey:@"ROAD NO"];             NSMutableString *str = [NSMutableString stringWithFormat:@"%@-%@",d1,d2];             [highwayArray addObject:str];             [highwayObjectArray addObject:result];             d1 = nil;             d2 = nil;         }         else if(i==3)         {             d1 = [d valueForKey:@"ROAD NAME ENGLISH"];             [avenueArray addObject:d1];             [avenueObjectArray addObject:result];             d1 = nil;         }         else if(i==4)         {             d1 = [d valueForKey:@"ROAD NAME ENGLISH"];             [roadLanesArray addObject:d1];             [roadLanesObjectArray addObject:result];             d1 = nil;         }         else if(i==5)         {             d1 = [d valueForKey:@"BLOCK NO"];             [blocksArray addObject:d1];             [blocksObjectArray addObject:result];             NSLog(@"%@",json);              }         else if(i==6)         {             d1 = [d valueForKey:@"AREA NAME ENGLISH"];             [areasArray addObject:d1];             [areasObjectArray addObject:result];             d1 = nil;         }         else if(i==7)         {             d1 = [d valueForKey:@"NAME ENGLISH"];             [governarateArray addObject:d1];             [governarateObjectArray addObject:result];         }                   if ([result.feature.geometry isKindOfClass:[AGSPoint class]]) {                 AGSSimpleMarkerSymbol *symbol = [AGSSimpleMarkerSymbol simpleMarkerSymbol];             symbol.style = AGSSimpleMarkerSymbolStyleDiamond;             result.feature.symbol = symbol;              }         else if ([result.feature.geometry isKindOfClass:[AGSPolyline class]]) {             AGSSimpleLineSymbol *symbol = [AGSSimpleLineSymbol simpleLineSymbol];             symbol.style = AGSSimpleLineSymbolStyleSolid;             symbol.color = [UIColor blueColor];             symbol.width = 2;             result.feature.symbol = symbol;              }         else if ([result.feature.geometry isKindOfClass:[AGSPolygon class]]) {                  AGSSimpleLineSymbol *outline = [AGSSimpleLineSymbol simpleLineSymbol];             outline.style = AGSSimpleLineSymbolStyleSolid;             outline.color = [UIColor redColor];             outline.width = 2;                  AGSSimpleFillSymbol *symbol = [AGSSimpleFillSymbol simpleFillSymbol];             symbol.outline = outline;                  result.feature.symbol = symbol;         }         [self.graphicsLayer addGraphic:result.feature];       }          }     - (BOOL) callout:(AGSCallout *)callout willShowForFeature:(id<AGSFeature>)feature layer:(AGSLayer<AGSHitTestable> *)layer mapPoint:(AGSPoint *)mapPoint {      return NO;               }          - (IBAction)list:(id)sender     {     tableViewController *tab =[[tableViewController alloc]init];         tab.poiArray1 = poiArray;         tab.addressArray1 = addressArray;         tab.areasArray1 = areasArray;         tab.governarateArray1 = governarateArray;         tab.blocksArray1 = blocksArray;         tab.highwayArray1 = highwayArray;         tab.avenueArray1 = avenueArray;         tab.roadLanesArray1 = roadLanesArray;              tab.poiObjectArray1 = poiObjectArray;         tab.addressObjectArray1 = addressObjectArray;         tab.areasObjectArray1 = areasObjectArray;         tab.governarateObjectArray1 = governarateObjectArray;         tab.blocksObjectArray1 = blocksObjectArray;         tab.highwayObjectArray1 = highwayObjectArray;         tab.avenueObjectArray1 = avenueObjectArray;         tab.roadLanesObjectArray1 = roadLanesObjectArray;         [self.navigationController pushViewController:tab animated:YES];      }      By Using the above code I Was able to show all the data matching with that key word in a table View section wise. Now my Problem is i want to Display Only One Location on the map whatever the user selected from that table. Can Any One Please me.


title: Merge from arcpy.ListTables() producing duplicate rows?

text: I was trying to manually merge file geodatabase tables and was getting Null values for the output file rows and decided to write a python script to merge the files. Here's my script: import arcpy               from arcpy import env     env.workspace = "G:/US county_Climate Hazard models/NWS_County/SUMCounty_NWS_new_frq.gdb"     arcpy.env.overwriteOutput = True     tableList = arcpy.ListTables()     arcpy.Merge_management(tableList, "C:/Documents/ArcGIS/Default.gdb/Merge_sum_freq")      The code takes some time to run but runs without errors. However, when I open ArcMap back again and check for the output, I get duplicate entries in the table I made. Here's how the output looks like: ![enter image description here](http://i.stack.imgur.com/QVbsM.jpg) As you see every entry is duplicated. Could someone explain why is that the case and how could I fix it? Is there something going on with the way I wrote the code? ****EDIT # 1****** I also added:               print tableList      and I got:               [u'CT_09001', u'CT_09003', u'CT_09005', u'Merge_sum_freq']      (I'm not sure why I get this "u" thing) ****EDIT #2***** The original question is now resolved thanks to Michael Miles-Stimson. I accidentally created Merge_sum_freq inside of SUMCounty_NWS_new_frq.gdb ****EDIT #3 NEW Question _*_** I have a performance related question: every geodatabase table I need to merge has 1 row and 6 fields (including OBJECTID), which you could see on the screenshot I posted. When I run the above script for 10 tables, it runs in 10sec. For 30 tables - 8sec, 264 tables - ~2min. In the end I need to merge >3000 tables (every table with 1 row, so final table will have >3000 rows). Not sure if I'm extrapolating this right, but it seems like they would complete in >20min. Is this considered a good performance? If not, is there any way I could optimize anything?


title: How to save a QGIS workspace with database link for offline use?

text: Is it possible to save current workspace of QGIS for offline usage like in mapinfo? **My workspace has link to database.** and when I reopen workspace unlink in mapinfo it can't open it.


title: Use GeoTIFF masks to optimise combine layers algorithm

text: My app is a C# program that takes in an extent from the user, then loads GeoTIFFs from a folder that overlap the extent. It layers them according to a simple number on each GeoTIFF, then exports the result based on the user- supplied extent and saves that cropped version to a new file. I use the Global Mapper SDK for this and it works great. I want to speed up the app by not loading any GeoTIFF which is already "covered" by other data. I notice that the GeoTIFF format has a mask. For example, I can see in the thumbnails for the GeoTIFF files in Windows Explorer a black background with a white foreground that is the actual data. Is there a way to read this mask into a .NET image or similar? Then I would be maintain a growing bitmap of all masks ORed together, and able to skip loading any GeoTIFF that did not increase the white area. That would greatly reduce the time taken for the GM SDK to do the export because there would be less total layers loaded.


title: How can I use MOBAC atlases in GIS Software (e.g., QGis, uDig)

text: How can I use the output from MOBAC (Mobile Atlas Creator) as background layer in open source desktop GIS software (e.g., QGis or uDig, TileMill etc.)? **My goal** : To put my question into context: I want to use my downloaded map tiles as background raster map. On top of this, I want to add GPX tracks to create a map with mountain bike routes. **My tiles** : I have a set of tiles downloaded using MOBAC. The tiles are in PNG format, each having a size of 256x256px. The projection is GlobalMercator (as in Google Maps). The naming scheme is `$(z)_$(x)_$(y).png`. I only have z (zoom level) 14. I can access the tiles through a local web server using a URL such as `http://localhost/~username/MyMaps/14_8622_5722.png` **The problem** : I can't find a way to feed output from MOBAC into GIS software. **Previous efforts** :   * I have successfully created an OpenLayers web map as follows:              var name = 'Test Layer';     var url = "/~username/MyMaps/${z}_${x}_${y}.png";     var options = {       type: 'png',       sphericalMercator: true,        isBaseLayer: false,        minZoomLevel: 14     };     var layer = new OpenLayers.Layer.XYZ(name, url, options);        * I have tried to directly add the image tiles as a layer in QGis. But the tiles are not georeferenced.    * I have exported the tiles to MbTiles format using Mobac. But - surprising to me - not even TileMill seems to be able to load the `.mbtiles` file as background layer.   * I have served the `.mbtiles` file using TileStream but could not access the TileStream server using QGis or uDig.   * I have tried adding the tiles as GeoServer store using the Image Mosaic plugin, but I could not get it to work. ...Can you provide hints and help as to which route to take? I.e., which MOBAC export format should I use so to be able to import the map data into GIS software?


title: How to enable GRASS plugin in QGIS 1.8.0 (Linux Mint 13 - 64 bit)?

text: Every time I enable GRASS plugin in QGIS (Plugins -> Manage Plugins -> GRASS), I got a message: Exception - G_getenv(): Variable GISDBASE not set.   I googled and found some information here but don't know how to "unset GISRC".   Anyone please show me how to do it?


title: V.Digit Grass Error

text: I am trying to slowly switch from ArcGIS to GRASS but the interactive V.Digit tool is giving me a headache. I am trying to add supervised classifications training samples (i.e. polygons) for use in a maximum likelihood classification but when I type,               v.digit -n map=Tsample bgcmd='d.rast map=landsat@image'      I receive the error...                'Option <map> does not accept multiple answers'       Removing the bgcmd shows me that this error is association with the bgcmd not the first map declaration. I assume that this error has something to do with my procedure as im rather new to GRASS than a bug... I am executing the function via the command console in the GRASS 6.4.3svn GUI.


title: v.surf.idw on grass return : no data points found

text: I'm a newbie with grass and I'm using grass (6.4.3) with the GUI... running Ubuntu 13.04 for what it's worth. Here is my region :               projection: 99 (Lambert Conformal Conic)     zone:       0     datum:      ** unknown (default: WGS84) **     ellipsoid:  wgs84     north:      871548     south:      871035     west:       176353     east:       176851     nsres:      1     ewres:      1     rows:       513     cols:       498     cells:      255474      Here are the five first lines of the files from which I created a vector map with `v.in.ascii` :               871207.125|176633.609375|322.54595947|1     871209.75|176633.84375|323.36795044|2     871207.4375|176623.8125|321.85391235|3     871213.75|176628.359375|322.46286011|4     871207.25|176658.4375|324.73287964|5      In table manager I have :               cat dbl_1 dbl_2 dbl_3 str_1     integer double precision double precision double precision character      I'd like to get a raster from this, so I'm trying to use `v.surf.idw`. I cannot understand what to put in the column argument of this function.               When I put dbl_3, I get :     6423 records selected from table     WARNING : No record for point (cat = 6424)     WARNING : No record for point (cat = 6425)     0 points loaded     ERROR: No data points found      When I put double precision it is not working, I've tried every ways I can think of. Thanks for reading me.


title: ArcGIS Javascript - Editor API methods?

text: I was looking at the API for the ArcGIS Editor widget (esri.dijit.editing.Editor), and I couldn't find a list of methods that are on the widget. I've seen other examples online that call destroy(), but that method isn't listed. Where can I find a complete listing of the Editor's methods?


title: Does r.in.gdal -o reproject the imported raster data to the location projection?

text: Does r.in.gdal -o actually reproject the imported raster data to the location projection, or do I need to call gdalwarp to reproject prior to importing?


title: How to merge TopoJSON features?

text: I currently have a coordinate based map of the US and it's states for topojson (https://github.com/mbostock/topojson). The coordinates are a JSON array (as topojson works), data looks like this:               [-123.233256,42.006186], [-122.378853,42.011663]...      Now what I'd like to do is combine states into regions, outside of doing this by hand (and that idea makes me feel sick), are there any tools which I can use to accomplish this? I have QGIS but have no clue where to even start (let alone how to combine states in QGIS).


title: How to setup default toolbars and menu item?

text: I need to customize the UI for a group of users at another office. I want to turn off a bunch of toolbars and menus that are not needed. I know how to save my configurations (Settings -> Customization...) to a INI file. How do I setup the other users computers to load the INI at startup? I am using Windows 7 and Qgis 2.0.1.


title: Layer name in OpenLayers 3.0

text: I'm looping through all the layers in my map using               map.getLayers().forEach(function(layer, i) {        ...     });      I'd like to get the name of each layer but I haven't been able to find anything in the examples or the documentation. Any pointers? I'm using OpenLayers 3.0.0.


title: Is it possible to pan the map in QGIS composer?

text: I am relatively new to the composer feature of QGIS, and often have difficulty getting the map centred exactly where I want it. In particular the area of interest often moves when I set the scale. Is it possible to pan or otherwise centre a plot that is easier than adjusting the extents manually. Thanks


title: Automate Select and Export to PDF using Arcpy Search Cursor

text: I'm trying to automate a process to select each feature of a feature class one by one, zoom to it and export as pdf. using the code below, the first feature is selected and the pdf is exported, and the loop stops there. How can I modify this code to loop through all records?               import arcpy          mxd = arcpy.mapping.MapDocument('CURRENT')     df = arcpy.mapping.ListDataFrames(mxd, "Layers") [0]     fc = "This_FC"           with arcpy.da.SearchCursor(fc, ("SHAPE@", "FeatureID", "OID@")) as cursor:         for row in cursor:             query_str = """ "FeatureID" = '{0}' """.format(row[1])             arcpy.SelectLayerByAttribute_management(fc, "NEW_SELECTION", query_str)              df.extent = row[0].extent             df.scale = df.scale * 5                  arcpy.RefreshActiveView()             arcpy.mapping.ExportToPDF(mxd, "C:/test/" + fc + "_" + str(row[1]) + "_" +  str(row[2]) + ".pdf")


title: How do you load a GeoJSON file in Leaflet?

text: I'm stumped as to how to add a GeoJSON layer to an OSM map. Through Google, I've found a bunch of different answers, none of which seem to work. Following the Leaflet documentation, the simplest way to do so is with this code:               var map = L.map('map').setView([51.505, -0.09], 15);          L.tileLayer('http://{s}.tile.openstreetmap.se/hydda/full/{z}/{x}/{y}.png', {     attribution: '© OpenStreetMapcontributors'     }).addTo(map);          L.geoJson(data, {             style: function (feature) {                 return { color: feature.properties.color };             },             onEachFeature: function (feature, layer) {                 layer.bindPopup(feature.properties.description);             }     }).addTo(map);      Questions: what exactly do I enter for the "data" variable? I keep a file in the path "~/Scripts/worldmap.geojson", but a direct URL string doesn't work. I'm using jQuery, and it's $.getJSON function does not seem to work.


title: Relational operation in eCognition

text: Is there anyone who has a experience with relational operation in eCognition. I'd like to ask 2 questions related to eCognition. 1\. How can I classify the pixel values between 0.23 and 0.46 as a class 1, the pixel values between 0.46 to 0.64 as a class 2 etc. How can I achieve this purpose in eCognition. 2\. Is it possible to use two variables for classification of one class. For example, can I use given ranges of temperature and precipitation to classify class 1. How should I build rule set for this purpose? Thanks a lot


title: Batch rectifying 200+ historic aerial images?

text: I have acquired historic aerial imagery for the surrounding towns, circa 1954, and the images have no georeference. I am looking for software to help me batch-rectify the 200+ images. I have current orthophotos, DEMs, and other vector information for the area.


title: opening older 9.x Esri .gdb in QGIS?

text: I want to open an Esri .gdb with QGIS. I have adjusted GDAL and it works fine with an example .gdb provided on another thread on the topic. I still cannot open the .gdb of interest to me and suspect it is because it may not be a 10.something gdb, but probably an older .gdb file version. I have currently no access to Esri software. Anybody out there knowing a solution to the problem?


title: PostGIS select sites within buffer

text: I have a WMS query (GeoServer) that returns info from my layer into a popup in Leaflet.JS.However I was wondering how would do this outside of leaflet, i.e. on a standard html/php page. I assume I would use the `ST_buffer` function in `POST GIS` when I do my query, order by distance and then limit by 1. For example I have a table with the_geom and lets say house name, How would I do a simple buffer that selects the nearest house to that given point which happens to be from the same table using the same geom column. I think if I can just get that working I should be able to extend the query myself to provide the info I need. Thanks


title: What is the cause and solution of this ArcMap tint band drawing issue?

text: I installed ArcGIS 10 SP1 last week. My graphics card meets the spec and I'm up to date on drivers. I have a layer of parcels that are tint-banded in the MXD. At some scales, straight lines that are the same color as the tint band shoot across the map. After this happens, the next, second, or third extent change crashes ArcMap (usually it's the 1st or 2nd). When I successfully change extent, I'll get a different set of straight lines. In the image below, only two of the polygons are (or should be) tint-banded: ![alt text](http://i.stack.imgur.com/HhnBx.png) Has anyone else experienced this and know the cause? And does anyone have a solution ( _besides not using tint bands_ )? * * * ## Update: I did a check/repair geometry, but the issue continues. I was able to replicate the issue with US Counties: here's a layer package of the US Counties with a green tint band. As I zoom in/out, I experience the same issue. ![alt text](http://i.stack.imgur.com/qlu58.png) Perhaps this is caused a result of the hacky way I create the tint bands ( _maybe that's good second question_ ).


title: Symbols showing sex of a species in a label in MapInfo

text: I am trying to indicate the sex of a species using the symbols ♀ and ♂ in MapInfo. I would like the symbols to display in the label of each species point, but I cannot find out how to do this. I had thought I could maybe use chr$(), but I cant find the code for these symbols.


title: SOLVED: "Open Source Routing Machine" works with browser but not with Python

text: What I have in Python:               import requests     url_osrm = 'http://router.project-osrm.org/viaroute?loc=50.12345,8.12345&loc=50.23456,8.23456'     print url_osrm     osrm = requests.get(url_osrm)     # print osrm.text      What I get in Python (Error): > requests.exceptions.ConnectionError: > HTTPConnectionPool(host='router.project-osrm.org', port=80): Max retries > exceeded with url: /viaroute?loc=50.1234 5,8.12345&loc=50.23456,8.23456 > (Caused by : [Errno 10061] Es konnte keine Verbindung hergestellt werden, da > der Zielcomputer die Verbindung verweigerte) In English: > [Errno 10061] No connection could be made because the target machine > actively refused it The same happens with urllib2. * * * If I paste the url in a browser, it works and I get a json-file back. * * * **Edit: SOLVED: My bad, I forgot the proxy:**               import requests     url_osrm = 'http://router.project-osrm.org/viaroute?loc=50.12345,8.12345&loc=50.23456,8.23456'     print url_osrm     proxies = {"http": "http://YOUR-PROXY:8080"}     osrm = requests.get(url_osrm, proxies=proxies)     # print osrm.text


title: How display outlines in a Map, only from a scale, but the colored polygon in all scales?

text: I want to display a shape with UMN-Mapserver. It has different coloured polygons and also outlines. But I want to display the outlines only, from a specific scale. I could add `MINSCALE 150000`, but in that case the polygon is completely not rendered, if that doesn't match. How I display the polygon without outline in that case?


title: Why Do I Get this Strange Serialization Exception?

text: I have the following data contract classes:               [DataContract]     public class MapServerInfo : ServerInfo     {         [DataMember]         public string mapName { get; set; }         [DataMember]         public string description { get; set; }         [DataMember]         public LayerInfo[] layers { get; set; }     }          [DataContract]     public abstract class ServerInfo     {         [DataMember]         public string currentVersion { get; set; }         [DataMember]         public string serviceDescription { get; set; }     }          [DataContract]     public class LayerInfo     {         [DataMember]         public int id { get; set; }         [DataMember]         public string name { get; set; }     }      And I try to serialize JSON data from a map service in the following class:               public class ServerInfoCollection     {         private EventHandler<AddCompleteEventArgs> _addComplete;         public event EventHandler<AddCompleteEventArgs> AddComplete         {             add             {                 if (this._addComplete == null || !this._addComplete.GetInvocationList().Contains(value))                 {                     this._addComplete += value;                 }             }             remove { this._addComplete -= value; }         }              private Dictionary<string, ServerInfo> _items = new Dictionary<string,ServerInfo>();              public int Count { get { return this._items.Count; } }              public ServerInfo this[string key]         {             get { return this._items[key]; }             set { this._items[key] = value; }         }              public void AddMapService(string serviceName, string serviceUrl)         {             string url = serviceUrl.TrimEnd("/".ToCharArray()) + "?f=json";             WebClient webClient = new WebClient();             webClient.OpenReadCompleted +=                 delegate(object sender, OpenReadCompletedEventArgs e)                 {                     if (e.Error == null)                     {                         ServerInfo info = null;                         DataContractJsonSerializer serializer = null;                              try                         {                             Stream responseStream = e.Result;                             serializer = new DataContractJsonSerializer(typeof(MapServerInfo));                             info = (MapServerInfo)serializer.ReadObject(responseStream);                         }                         catch (Exception ex) { string message = ex.Message; }                         finally                         {                             this._items.Add(serviceName, info);                             this.TriggerAddComplete(new AddCompleteEventArgs(serviceName, serviceUrl));                         }                     }                     else { string message = e.Error.Message; }                 };                   webClient.OpenReadAsync(new Uri(url));         }              public void AddGeocodeService(string name, string url)         {         }              public void AddFeatureService(string name, string url) { }              private void TriggerAddComplete(AddCompleteEventArgs args)         {             if (this._addComplete != null) { this._addComplete(this, args); }         }     }      Whenever I run the app in debug mode I get the Visual Studio Just-In-Time Debugger popping up saying, "An unhandled exception ('Unhandled Error in Silverlight Application Code: 4004 Category: ManagedRuntimeError Message: System.Runtime.Serialization.SerializationException: 'html>". The exception is thrown right after the curly brace that comes after my finally block under AddMapService(), just before the else statement. Nothing is caught in my catch statement. At this page: http://msdn.microsoft.com/en- us/library/bb412179%28v=vs.110%29.aspx it says, "The JSON serializer throws a serialization exception for data contracts that have multiple members with the same name." This is not the case. I eliminated all of the properties I didn't need and the ones that remain do not have the same name. I asked this question at StackOverflow and the one answer I got was that the presence of 'html' and 'head' in the exception indicates that the data is not JSON data. This cannot be the case since 1) I can check the string that gets returned and see that it's not (X)HTML but json and 2) my properties actually get populated as I expect them to with the correct data at info = (MapServerInfo)serializer.ReadObject(responseStream);. However, the exception still gets thrown and my catch block does not catch it. Does anyone know what might be going on with this? How can I get more detailed exception info than what shows up in the 'Just-In-Time Debugger' window?


title: "Open Source Routing Machine" querying and parsing

text: What I have in Python:               import requests     url_osrm = 'http://router.project-osrm.org/viaroute?loc=50.12345,8.12345&loc=50.23456,8.23456'     print url_osrm     proxies = {"http": "http://YOUR-PROXY:8080"}     osrm = requests.get(url_osrm, proxies=proxies)     print osrm.text      * * * Question A: How many requests can I send until I get blocked? * * * Question B: How can I extract the first lon/lat from the json-output?


title: Generate Mobile Map Cache from ArcObjects

text: As of ArcGIS 10, is there a way to generate a mobile map cache using ArcObjects? This can be done in a stand-alone Python script: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//003z00000005000000.htm Thanks, Seth


title: Error message from output workspace when using script tool

text: I am having trouble running a script tool that I just created in Arc 10.1. The script runs fine in PythonWin, but when I create a script tool with it, the tool seems to not want to accept folders or geodatabases as its output workspace. I have the output parameter (argument) set as both an output and a workspace in the script tool, but it still keeps giving me an error that the "Dataset.....already exists" I'm a new user, so I can't post a picture, but the error is 000725 Here is my script, if that helps out.               import arcpy     #allow for overwrites     arcpy.env.overwriteOutput = True          #set the workspace     inWorkspace = arcpy.GetParameterAsText (0)          #set the erase feature     eraseFeature = arcpy.GetParameterAsText (1)          #set the output workspace     outWorkspace = arcpy.GetParameterAsText (2)               #get a list of all the features in the workspace     arcpy.env.workspace = inWorkspace     featureClassList = arcpy.ListFeatureClasses()          try:              #loop through all of the features in the workspace         for featureClass in featureClassList:                  #construct the output path             outEraseFeature = outWorkspace + "\\erase_" + featureClass                  #perform erase             arcpy.Erase_analysis(featureClass, eraseFeature, outEraseFeature)                  arcpy.AddMessage("Wrote clipped file " + outEraseFeature + ". ")             print "Wrote clipped file " + outEraseFeature + ". "          except:              # Report if there was an error         arcpy.AddError("Could not erase feature classes")         print "Could not erase feature classes"         print arcpy.GetMessages()      Any help would be much appreciated. Thanks, in advance. Neil Loehlein Portland, OR


title: How can I get a road network to run a simulation on?

text: I'm new to GIS and I'm asking this a bit out of desperation, but I'm trying my best not to get overwhelmed. :) I want to download a file that contains data of a road network (of what road network is not a concern to me), visualize this road network, but then be able to run a simulation on it (so it can't be an image, i.e. a .jpeg). I'm looking for what resources/software I should use. (Is there a way I can chat with someone who knows about GIS? I feel that would be better.)


title: How are Google serving up their Styled Maps?

text: As we all know, Google changed web mapping when they introduced Google Maps. At the time of writing this, tiled maps are now the norm with web mapping. For example, with ArcGIS Server, you cache as many of your layers as possible into cached map services, and then if need be, you pull in any of your operational layers via either dynamic or feature services. Once you have a cached service, there is no way to tweak the cartography without having to build a new cache again. The v3 Google Maps API allows you to style the Google Tiled Basemap. The styling wizard lets you tinker nearly all aspects of the Google Map, to suit your needs. I can understand from an API perspective, how you can send your style request in a JSON object. **What I want to know is, how are they doing this at the back end?** (Would be good to include a reference) If I look at the Styling Wizard with Firebug, and make changes to anything, I can see that new GET requests are made to return PNG images back to the client. ![enter image description here](http://i.stack.imgur.com/gLg9K.png) There are near infinite variations you could apply, so they cannot have all the tiles pre-cooked and ready to serve out. Therefore, I can only assume they are creating these tiles on demand?


title: Debugging AGS API for Silverlight - returned featureset to graphics layers not drawing

text: I have geometries returning from the ArcGIS Server as a featureset from a gp service to the silverlight client. In debug mode in VS 2010 I see the geometries coming back to the client and I am adding them to a graphics layer that I have defined in the xaml. I figure I am doing something wrong but I can't figure out how to debug why the graphics aren't showing up in the map. foreach (Graphic graphic in gpLayer.FeatureSet.Features) {                                   ESRI.ArcGIS.Client.Geometry.MapPoint point = graphic.Geometry as ESRI.ArcGIS.Client.Geometry.MapPoint;                         point.SpatialReference = MyMap.SpatialReference;                         {                             graphic.Symbol = LayoutRoot.Resources["NewPointSymbol"] as ESRI.ArcGIS.Client.Symbols.Symbol;                             graphic.Geometry = point;                         };                              //MapPoint geographicPoint = _mercator.ToGeographic(point) as MapPoint;                         _graphicsLayer.Graphics.Add(graphic);      Could this be a scope issue or a projection issue?


title: Can one edit xls, and convert xls to csv - with a portable python or Perl module?

text: I am on a network and computer where nothing can be installed. On this computer i run portable apps through a USB, however I need also to to edit xls, and save as csv. Since i cannot install anything the usual xlrd xlwt wont work, except if they can be made portable? Any suggestion how to solve it?


title: Automatically saving ModelBuilder output in the same folder it used for the input?

text: I am trying to create a model that takes an input layer, which will always be located in a unique folder on my hard drive, and spit out a set of analyses in the same unique folder the input was referenced from. For example, I want to be able to add "X.shp" to the map from C:\Folder1, run the model, and have the outputs be saved into Folder1 automatically. If "X.shp" was taken from C:\Folder2, I would want the output to be saved automatically in Folder2. Is it possible to automate this?


title: How to format labels with leading zeroes?

text: I'm trying to display a number next to a symbol in QGIS 2.0.1 - The numbers are from an integer-field in a postgresql database. The goal is that single digit numbers get a leading 0. Examples:   * 12 --> 12   * 10 --> 10   * 9 --> 09   * 0 --> 00 In the menu                layer properties --> labelling --> formating      I switch on "number formatting" and then enter a custom expression. I tried these two codes with various combinations of ", ' and more brackets, but none of them work:   * 1              CASE                WHEN "Value" <10 THEN rpad("Value", 2, 0)          END        * 2              CASE                WHEN "Value" <10 THEN "0" || "Value"          END      EDIT: For any other person having trouble with this: **Turn of "number formatting".** It seems to be the reason why leading zeros are removed from the rpad function. Enter the following in the regular "layer labeling" box:                 rpad("Field", 2, 0)      and everything will work fine. But if you also activate number formatting it will take away the leading zeros.


title: Export additional data with ogr2ogr

text: I am trying to import my ESRI-shapefile with OGR2OGR to an existing mysql table. This works fine but unfortunately only the featureID and the geometry data is exported. How do I add the labels of the different polygons? These are stored in the .dbf-File if I am not mistaken. I am using the following command:               ogr2ogr -f MySQL MYSQL:Polygons,host=localhost,user=user,password=user,port=3306 -nln route -a_srs EPSG:4326 -append /home/Fields3.shp      if using ogrinfo on the file it tells me:               INFO: Open of `Fields3.shp'           using driver `ESRI Shapefile' successful.     1: Fields3 (Polygon)      How can I get the addition information into Mysql? The Shapefile is created using Qgis 2.0 with standard settings.


title: How to convert line simple features to topological network?

text: using Gdal/ OGR is there a way to convert a shp, kml, or PostGIS line layer into a network of nodes and segments so it can be used in packages such as networkx?


title: How to find out the Stream Order from DEM

text: I have DEM Data for the study area, I want to know how to find out the stream order using DEM. I saw Stream order tool in Hydrology toolset in ArcGIS. in Stream order tool what is input stream raster and how to create it, which method of stream ordering is mostly used ![](http://i.stack.imgur.com/B6qDb.jpg) Thanks in Advance for helping me


title: QgsVectorDataProvider.createAttributeIndex() fails in QGis 1.8

text: Adding an Attribute Index to a Shapefile fails for a Python Plugin in QGIS: For example - the following python Code in a Plugin               ly = QgsVectorLayer(data_source, layer_name,"ogr")     prov = ly.dataProvider()     prov.createAttributeIndex(i)      results in the error Message: > 'QgsVectorDataProvider' object has no attribute 'createAttributeIndex' Whereas Adding a Spatial Index works without any problem               ly = QgsVectorLayer(data_source, layer_name,"ogr")     prov = ly.dataProvider()     prov.createSpatialIndex()      Is this a Bug? Thanks in Advance


title: Ways to load OS VectorMap data

text: Us lucky fellows in the UK have been able to access public geodata from Ordance Survey under an open licence since April this year. Some of the datasets, including VectorMap, are broken into 100 files per grid square, each containing multiple layers, in this case in shape format. To cover the wider London area with VM you need to load 400 x 12 or so layers, with each layer identically named. Any tips on the best ways to do this, preferably using free or open source tools? http://www.ordnancesurvey.co.uk/oswebsite/opendata/


title: OpenLayers filter: OGC (xml) or CQL?

text: While having not used filtering that much in _OpenLayers_ , we are discussing which type of filtering we should use: The **OGC** XML filter or the more simple **CQL** filter?     * Are one of them more appropriate than the other?    * What about support from various frameworks and servers - what kind of trouble will we encounter? We are using _GeoServer 2.x_ and _OpenLayers 2.9.x_ , but we are also querying external servers for data.


title: Failing install GeoCoding 2.6

text: I'm a beginner in Qgis... trying to install GeoCoding 2.6. I have this Python error : > Impossible de charger l'extension 'GeoCoding' depuis ['C:/PROGRA~2/Quantum > GIS Lisboa/apps/qgis/./python', 'C:/Users/CÃ©line/.qgis//python', > 'C:/Users/CÃ©line/.qgis//python/plugins', 'C:/PROGRA~2/Quantum GIS > Lisboa/apps/qgis/./python/plugins', 'C:\PROGRA~2\Quantum GIS > Lisboa\bin\python27.zip', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\Python27\DLLs', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\Python27\lib', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\Python27\lib\plat-win', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\Python27\lib\lib-tk', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\qgis\bin', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\PIL', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32\lib', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\Pythonwin', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\wx-2.8-msw- > unicode', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\qgis\python\plugins\fTools\tools'] > > Traceback (most recent call last): File "C:/PROGRA~2/Quantum GIS > Lisboa/apps/qgis/./python\qgis\utils.py", line 143, in loadPlugin **import** > (packageName) File "C:/PROGRA~2/Quantum GIS > Lisboa/apps/qgis/./python\qgis\utils.py", line 309, in _import mod = > _builtin_import(name, globals, locals, fromlist, level) ImportError: No > module named GeoCoding Version de Python : 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)] Version de QGIS : 1.8.0-Lisboa Lisboa, f440b60 > Chemin vers Python : ['C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python', > 'C:/Users/C\xc3\xa9line/.qgis//python', > 'C:/Users/C\xc3\xa9line/.qgis//python/plugins', 'C:/PROGRA~2/Quantum GIS > Lisboa/apps/qgis/./python/plugins', 'C:\PROGRA~2\Quantum GIS > Lisboa\bin\python27.zip', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\Python27\DLLs', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\Python27\lib', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\Python27\lib\plat-win', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\Python27\lib\lib-tk', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\qgis\bin', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\PIL', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32\lib', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\Pythonwin', > 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\wx-2.8-msw- > unicode', 'C:\PROGRA~2\Quantum GIS > Lisboa\apps\qgis\python\plugins\fTools\tools'] Thanks for your help!


title: Googlemaps Layer on Openscales

text: I'm using openscales in a web mapping project, i would like to add a Google maps layer, i know that there's some licence problems , but i don't care for these licence issues, i'm looking to add googlemaps layer to openscales like it has been made on OpenLayers. Thanks in advance


title: Is there a way to label a layer using its file name or layer name?

text: Is there a way to label a layer using its file name or layer name (if I'd change it in Layers window TOC)? I tried some functions but I'm not good on that. I have lots of shapes (lines) and it would be easy to identify them this way. Of course, this name could be a table attribute but I'm trying a alternative way (maybe more dynamic). Thanks.


title: MMFeederSource Usage

text: I am trying to to write a custom feeder tracing routine, but facing some issues.   I am using the following ArcFM method: _IMMFeederTracer.TraceFeeder(IGeometricNetwork, IMMFeederSource)_ _MMFeederSource_ is simply initiated as _new MMFeederSourceClass();_.   There does not seem to be any QI from _IRow_ and _IFeature_ to set the source. In the current form, the Feeder Manager starts but after a while throw an HResult 0x80044020 error. Any ideas how to initialize the feeder source?   I could not find anything on the ArcFM Resource Center. Thanks in advance   Ujjwal


title: Shapfile created in QGIS does not fit in ecognition project

text: I have various shapefiles - some created by others, some created by my self. They should all have the same size of the underlying Raster file. The shapefiles have been created on top of a Raster file (in QGIS) with Coordinate System WGS84 /UTM zone 45N. I checked the Options in QGIS - Shapes and Raster use the same Coordinate System. Wenn importing the Shapefiles to Erdas Imagine, everything looks fine. In eCognition the shape files doesn't fit to the extend of the Raster file. Where can I solve this problem? Do I have to change something in QGIS or is it a eCognition problem? ![Shapefile created in QGIS - Displayed in eCognition](http://i.stack.imgur.com/9n9ax.png)


title: Importing x y coordinates from spreadsheet in ArcGIS

text: I converted degree, min, sec to decimal coordinates. I made a simple spreadsheet with X and Y columns and one set of coordinates: x= -73.533678 and y = 41.05523 I added the data (File, add data, add x,y data) and it always plots to 0,0. I'm in the projection of the base map. What am I doing wrong???im


title: Restrict featureinfo fetched with the WMSGetFeatureInfo request in geoserver

text: i followed a GeoEXT example to create popup menu that gets WMSFeatureinfo on click on the map, What i want to do is to get some specific fields not all the fields, i try the propertyName (with name1,name2 and name3) and defined the wanted fields but id doesn't work, also i want to change the title by bringing a specific value not the fid but it doesn't work, here's the code: am new to this plz how to restrict the fields fetched from WMSGetFeatureInfo request                feature Control, a "button" control         action = new GeoExt.Action({         control: new OpenLayers.Control.WMSGetFeatureInfo({     autoActivate: true,     infoFormat: "application/vnd.ogc.gml",     maxFeatures: 6,     layers: [layers1],     propertyName: "name1,name2,name3",           eventListeners: {         "getfeatureinfo": function(e) {             var items = [];             Ext.each(e.features, function(feature) {                      items.push({                    xtype: "propertygrid",                    title: feature.fid, // here i change it to feature.name1 it doesn't work                     source: feature.attributes                      });             });             new GeoExt.Popup({                 title: "infos",                 width: 250,                 height: 300,                 layout: "accordion",                 map: map,                 location: e.xy,                 items: items             }).show();           }         }        }),         map: map,         text: "infos",         tooltip: "Infos "           });         actions["infos"] = action;         toolbarItems.push(action);         toolbarItems.push("-");


title: How to get custom color palette in rasterVis plot3D in R?

text: I am trying to insert an RColorBrewer palette to visualize an elevation raster in R using the rasterVis function `plot3D`. The rasterVis reference document shows the following parameters:               plot3D(x, maxpixels=100000, zfac=1, drape=NULL, col, rev=FALSE, adjust=TRUE, ...)      However, when I use `col=` I get the following error:               Error in col(zlen) :        a matrix-like object is required as argument to 'row/col'      When I use `color=` I get:               Error in f(...) :        formal argument "color" matched by multiple actual arguments      I've also tried this:               color.palette = brewer.pal(11, "PRGn")      which renders the 3D-plot, but uses the default color palette. This article does not provide a solution. Nor have I found anything useful on GISSE or StackOverflow: Any suggestions, assistance clarifying the errors, or input is much appreciated.


title: how can i read an attribute value from qgsfeature?

text: How can i read an attribute value from qgsfeature ? i am trying to update some attribute using other and i have to read value from attribute in the qgsfeature . my source code:               provider = self.theLayerToEdit.dataProvider()     Feature = QgsFeature()     while provider.nextFeature(Feature):         dACREVal =  Feature.attributeMap()["ACRE"].toDouble()         immwaterVal =  Feature.attributeMap()["mmwater"].toInt()         dCNValue = Feature.attributeMap()["CN"].toDouble()         dCarbonFixeValue =  Feature.attributeMap()["CarbonFixe"].toDouble()         dDenitrificValue = Feature.attributeMap()["Denetrific"].toDouble()         dFloodProtVal = iWaterRetentionValue * (dACREVal * 43560) * immwaterVal *  0.00328084         dStableClim = dCarbonFixeValue * iStableClimate         dUsableWater = dDenitrificValue * iUsableWater         Feature.changeAttribute(iFloodProt,QVariant(dFloodProtVal))


title: Removing time from date field labeling in VBScript

text: I'm not at all experienced with writing VBScript for labeling purposes. Here is my code so far:               Function FindLabel ( [Date_of_Birth] )       if not isnull ( [Date_of_Birth] ) then              FindLabel = "DOB: "  &  [Date_of_Birth]                 FormatDateTime [Date_of_Birth],????       else              FindLabel = ""       end if     End Function      Basically, I want to label features that have a DOB with the DOB from the field, and if that field is null I don't want it to label. That part works great, but I noticed that some of the dates also had a time tacked on to the end. I would like to format that part out if possible. The line that reads FormatDateTime [Date_of_Birth],???? is what I found searching online, but nothing I've entered after that comma works. I'm curious if I'm on the right track or not with how I'm trying to format the date. Thanks, Chase


title: How to select vector features with partial (>33%) overlap with another vector feature?

text: I have a grid and a range map, and would like to select those grid squares which overlap at least 33% with the range map. The range map is a single vector feature. (I'm new to QGIS but have some Arc experience) Thanks!


title: Convert Graphics to Vector in QGIS

text: I am wondering if there is a way to convert a graphic produced from symbology in QGIS to a vector layer (ie: shapefile). I have a point layer with bearing direction in the attribute table. Following the answer from How to display an arrow pointing to a specific bearing in QGIS, I have been able to make a line graphic depicting the look direction from a photo collection. I would like the line graphics to be converted into a vector layer, so I can perform some spatial analysis. Any thoughts? ![graphic representations](http://i.stack.imgur.com/Bbv26.jpg) Thanks....


title: GeoWebCache - WMTS Directory Structure

text: I am trying to convert an existing tile cache (created by custom process) into a usable format/structure for a GeoWebCache WMTS layer. I cant seem to figure out the directory/naming structure required - using a layer created and seeded by GeoWebCache as a reference. The first level is easy to understand EPSG_4326_00 - that would be spatial reference and tile matrix level. The next level of directories, and the names of the tiles themselves is totally alien to me. If anyone could explain - or point me at some documentation - that would be greatly appreciated.


title: QGIS2--How long to open 3.5GB File.GDB directory vector layer? Seconds, minutes, hours?

text: I am a new GIS user, and have started with QGIS 2.01 Dufour. I am trying to open an ESRI-style File.GDB (v10) directory that has a 3.5GB gdbtable file that is of Texas state hypsometry/topography. The machine is clearly busy as it tries to open this layer because the hard disk is spinning, and one of the cpu cores shows elevated activity at ~20%. But even after waiting up to 25 minutes the layer has still not opened successfully (one small exception below). As I am new to GIS, I simply don't know what to expect--should I expect opening such a layer to take seconds, minutes, or hours? I know this is not a straightforward question given so many details are involved, but I provide more info below. I am using a fairly capable system. It is a notebook with an Intel Core i7 3632QM Quad Core cpu, 8GB ram, and Windows 8.1 (64-bit). I am running the 64-bit version of QGIS 2.01 that was installed with the standalone installer. The File.GDB directory vector layer that I am trying to open is on the local hard drive (so it is not constrained by USB/other transfer limits). I mentioned above that there was one small exception to opening this 3.5GB File.GDB directory vector layer successfully, and it opened essentially instantly (well, within probably seconds). This was the very first time I tried to open it. I was zoomed in to an area of Texas just a few kilometers wide, and I went through the process of adding this layer to my project. Woo- hoo, it opened fine within seconds, and I could see the topography in the area of my viewpoint. Then I had the "bright" idea of zooming to the extents of the 3.5GB topography layer (i.e. the whole state of Texas), and that's where the problem started. The machine tried to refresh the view, and after five minutes of a blank screen and spinning blue ball, I thought there was a problem so I force-closed QGIS via the Task Manager. Now, every time when I attempt to open this 3.5GB topography layer, the machine puts itself to work, and even after waiting up to 25 minutes, the layer has still not opened. This happens whether I try to open the layer directly, and thus, I should be at the default full extents. And it also happens when I start with a zoomed in view of several kilometers width (like the first and only time it opened successfully), and try to add it as a layer so the machine should not try to render the full extents of the layer. I have twice deleted this File.GDB directory, and unzipped a completely new, fresh copy, and this still happens. I have deleted the .qgis2 directory in my home directory, and also the main QGIS key in the Windows registry. But again, this still happens. Anyway, this gets back to my question for you guys. How long should it take to open a layer like this? Should it be seconds, minutes, or hours? Can you provide any other wisdom or advice given that I have already refreshed the File.GDB directory (twice) and I have completely reset my QGIS defaults with no success? Thanks in advance for your responses! Mark * * * Ok, I am adding some commentary to my question after seeing responses from Get Spatial and Burton449. Here it is... Wow, incredible community! Thank you so much, Get Spatial, for your comprehensive replies, and also to you, Burton449, for the helpful suggestion. So it seems like I must have been positioned and zoomed just right that when I opened this layer originally, it happened quickly. And that waiting many minutes to hours to potentially days (I hope not!) might be in order for a data set of this size. From the various suggestions that were made like chopping the layer into quadrants, simplifying the lines via point reduction, filtering on a subset of the contours, etc., it seems like I would have to open the layer first to perform any of these operations. So this would mean waiting it out at least once. Or is there a way that I could perform some of these operations before actually loading the layer, esp. for something like the contour filtering operation? In layman's terms, I'm asking if I can set up something like a filter rule first (open only every fifth contour line, for example), and then open the layer, and have this filter applied on-the-fly. I'm not averse to command line tools so I am wondering if there is something like a cookie-cutter operation (or Python script) that I could use offline to extract a subregion of this larger layer, and spit it out to a new layer. And preferably in a format that is not FileGDB since that seems to be proprietary. My whole goal is to avoid this lengthy opening process (even if it is only once) because I simply have no idea if we are talking about hours or days. Indeed, since writing my original question, I have let QGIS chomp on this layer, and after 1.5 hours it is still thinking. Anyway, thanks again for your help, and for this wonderful community! Mark


title: Is there a way to use an existing Elastic Block Store (EBS) with an FME Cloud instance?

text: I am seeing help topics for all items except `Resources`. What is the best way to make EBS data available to FME Cloud instance? Would I need to make it available via one of my own EC2 instances? Or download it locally and upload via the usual FME-Cloud upload? ![enter image description here](http://i.stack.imgur.com/wQr01.png)


title: Symbolizer with external GIF image is not rendered pixel-perfect

text: We want to pinpoint some locations on a map with some external GIF symbolizers. The issue we have is visible in the attached image*: looks like somethings happens on a sub-pixel level, so that the resulting pluses look mutated :), most of them are somehow distorted/blurred/interpolated. The original external GIF is shown in the right lower corner of the image* (simple and sharp 5x5 pixel "+" symbol). Thank you in advance for any clue or solution. Regards, skipr Original image: ![enter image description here](http://i.stack.imgur.com/cR2eD.png) Magnified image: ![enter image description here](http://i.stack.imgur.com/9CtuC.png) Update: Same thing happens with well known symbolizers (e.g. MARK_CROSS):               ...     final Mark mark = sb.createMark(StyleBuilder.MARK_CROSS, Color.red, Color.yellow, 1);     final Graphic graphics = sb.createGraphic(null, mark, null, 1, 10, 0);     final Style style2 = sb.createStyle(sb.createPointSymbolizer(graphics));     ...      ![enter image description here](http://i.stack.imgur.com/PvPQt.png)


title: How can I save a column from the attribute table of a vector layer to a variable?

text: I am trying to write a python plugin for QGIS that can convert three vector layers into an object File for an urban climate model. Unfortunately I am stuck again with a possibly easy question: How can I save a column from the attribute table of a vector layer to a variable (for later use)? For example, I have a layer with buildings (`build`) with the hight of the buildings in field 35. I want to save all the building hights into a variable (`theINPolyHField`). I tried the following:               build = canvas.layer(2)     build_prov = build.dataProvider()          theINPolyHField = build_prov.select([35], build.extent(), False, False) # but returns "None"     build.selectedFeatureCount() # returns 0      What am I doing wrong?


title: Geodjango distance between two points appears incorrect

text: I want to show you something. I am using django+postgis and I am experiencing wrong distance values between two points. This is my model:               class Location(models.Model):         name = models.CharField(max_length=255)         latitude = models.FloatField(blank=True, null=True, verbose_name='Latitude')         longitude = models.FloatField(blank=True, null=True, verbose_name='Longitude')         location = models.PointField(blank = True, null=True)              objects = models.GeoManager()              def save(self, *args, **kwargs):             self.location = Point(self.latitude, self.longitude)             super(Location, self).save(*args, **kwargs)      This is the data:               bogota = Location(name='bogota', latitude='4.61779093242', longitude='-74.0893133545')     bogota.save()     barranquilla = Location(name='barranquilla', latitude='10.9731710819', longitude='-74.8020513916')     barranquilla.save()      This is the query:               ls = Location.objects.distance(bog.location).order_by('distance')     for l in ls: l.name; l.distance.km      This is the output:               Out[29]: u'BOGOTA'     Out[29]: 0.0     Out[29]: u'BARRANQUILLA'     Out[29]: 205.27667079004999      But Barranquilla city is to ~700km from Bogota city, see here: "http://www.travelmath.com/flight- distance/from/Bogota,+Colombia/to/Barranquilla,+Colombia" So, what is wrong? the SRID are the default values.


title: OpenLayers 2.12 + click handler + get the name of the vector

text:
Basically, I add several vector layers (KML content got from jQuery GET) in this way:               var layer1Vector= new OpenLayers.Layer.Vector('activeInfo');     layer1Vector.addFeatures(featuresLayer1Vector);     mp.addLayer(layer1Vector);      It works quite well. Instead of adding the               selectControlKml = new OpenLayers.Control.SelectFeature(         [layer1, layer2, ***], {             clickout: true, toggle: false,             multiple: false, hover: false,             toggleKey: 'ctrlKey',             multipleKey: 'shiftKey'         }     );      I knew that it's possible to add a "general" click handler event               OpenLayers.Control.Click = OpenLayers.Class(OpenLayers.Control, {         defaultHandlerOptions:{             'single': true,             'double': false,             'pixelTolerance': 5,             'stopSingle': true,             'stopDouble': false,             'delay': 0         },           initialize: function(options) {             this.handlerOptions = OpenLayers.Util.extend({}, this.defaultHandlerOptions);             OpenLayers.Control.prototype.initialize.apply(this, arguments);              this.handler = new OpenLayers.Handler.Click(                 this, {'click': this.trigger},                  this.handlerOptions             );         },         trigger: function(e) {                  if (e.target._featureId) {                 //var feature = vectorLayer.getFeatureById(e.target._featureId);                 alert('click on fea' + e.target._featureId)                 //$(e.target).css('fill', '#000000');                 //$(e.target).css('cursor', 'wait');                 //$("div#info").append("You just clicked on " + feature.id + "
");             }     ***         If I click on a vector layer I have to open different custom popups (qTip2). The code:               if (e.target._featureId) {      works, but how can I get the layer name based on the e object? I analyzed all the methods/properties, but without having success. E.g. One of my layers was defined as:               var layer1Vector= new OpenLayers.Layer.Vector('activeInfo');      I need to get 'activeInfo'. Thanks a lot!


title: Technicial Features of server to hold Quantum GIS, geoserver and geonetwork

text: I work at my company as Project Assistance and we think of migrating of ArcGis to another open source program similar to ArcGIS like Quantum GIS, geoserver and geonetwork and I need to know the technical characteristics of a server to do this and buy one. These tools let me to publish maps, to process maps, I guess. I want to know which are technical specifications that a server needs to support the job of these packages (QGIS, geoserver and geonetwork) for example, recommended memory for 40 or 50 users working with maps to they don't have any problem with their jobs, about storage space of 40 terabytes or more, etc. Thanks you in advance for your help. I want for 40 or 50 users publishing maps, processing maps, sharing data, availability in disk such as mirror, storage space for maps 40 or more terabytes, memory I don't know, how many do you suggest me? of which depends on the size of the maps? I am not restricted in budget, the important is that server was high performance and it was dedicated and centralized all maps, actually everybody makes his map and it's a problem. It's for a project to national level for disasters, earthquakes, etc. I'm sorry for my writing in English. I speak Spanish. I'm new and I can't anwer my post until tomorrow so I edit my first post. * * * On http://es.wikipedia.org/wiki/OpenStreetMap says: "El conjunto de servidores consta de un servidor de base de datos de gran rendimiento, un servidor de aplicaciones para el sitio web, tres servidores de aplicaciones para la API y un servidor destinado al renderizado del mapa" Do I also need an application server and another application server for rendering map? application server (GeoServer?) application server for rendering map (Geonetwork?) application server for API ?????? for example, what application? Thanks you * * * That is I don't know. "renderizado del mapa"


title: What software programs can run a supervised classification on an image segmented into objects?

text: I am trying to run a supervised MLC classification on a multispectral remotely sensed image, however, I want the image to be classified using a shapefile of segmented image objects (as opposed to the traditional pixel-based method). I have created the segmented shapefile in eCognition (Definiens) and can export it to ERDAS, PCI or ArcMap along with the original image file I am just not sure if these programs have the capability of applying the shapefile to classify the image as objects. From my searching it seems like IDRISI has this function (SEGTRAIN/SEGCLASS) but I do not have access to this program. Does anyone know if this is possible in ERDAS, PCI or ArcMap? Any help would be appreciated! Thanks.


title: Using Python, Batch Project output won't drop in a feature dataset (ArcGIS 9.3 or 10)

text: When attempting to batch project in ArcGIS (using either arcpy or arcgisscripting), I cannot get my geoprocessing results to drop into the feature dataset that I specify. The results instead are dropped in the root of the Geodatabase (*.mdb). Here's an abbreviated version of what I have so far:               import arcgisscripting, os     gp = arcgisscripting.create(9.3)          #set date     dmy = datetime.now()     currentDate = str(dmy.month) + "-" + str(dmy.day) + "-" + str(dmy.year)          outputfolder = gp.getparameterastext(0) #where we're going to drop the gdb     citydata = gp.getparameterastext(1)          gp.workspace = outputfolder     projectedfolder = os.path.join(outputfolder, "projected")     gp.CreateFolder_management(*os.path.split(projectedfolder))     gdb = "Projected_Data_" + currentDate + ".mdb"          #set coordinate system     outputCS = gp.CreateObject("SpatialReference")     outputCS.createFromFile(os.path.join(gp.getinstallinfo()['InstallDir'], "Coordinate Systems\\Geographic Coordinate Systems\\North America\\NAD 1983.prj"))     transformation = ""          gp.AddMessage("Creating City feature dataset")     gp.CreateFeatureDataset(os.path.join(projectedfolder, gdb), "City", outputCS)          gp.AddMessage("Loading in City feature classes")     arcpy.BatchProject_management(citydata, os.path.join(projectedfolder, gdb, "City"), outputCS, "", transformation)      Obviously, the part I'm having trouble with is the Output Path in the last line. I've also hard-coded the path to the "City" feature dataset, but still drops it in the root. Any ideas? (Cross-posted at ESRI Forums here)


title: Showing Google Maps using ArcObjects?

text: Do you know how we can show Google Maps as a base map in arcObject stand alone map? For example is there any wms address from google map that we can use it to connect to it or somthing like that? I have used open layers before and it has ability to set google maps as base map..But I want to use google maps in my stand alone application without using GeoServer and .. Do you have any suggestions?


title: How can I use Python to save several MXDs within the same script?

text: I have a program (720 lines of code) that runs multiple tasks on Arcmap .mxd files. Such as update layers, update text, add layers, etc. Within my program I have multiple "mxd.save()" instances. All of which occur within for loops. For example:               for item in mxds:         mxd = arcpy.mapping.MapDocument(item)         df=arcpy.mapping.ListDataFrames(mxd,"Project Area")[0]         print "HOWDY!"                       mxd.save()        Question is can I just have at the very end of my program have one for loop that saves all the mxds I was working on, while still retaining all the changes? UPDATE: Problem I think am having is this. I have one of my for loops that adds a layer to the map. If I do not mxd.save() after that layer is added. My program later when it attempts to find the layer that was added cannot find it. I guess it's not held in memory. UPDATE2: THANKS! I was doing something similar, but wrong none the less. Each time I had a for loop I would create a new .MapDocument object which I realized (thanks to you) was unnecessary. I now define the mapdocument object once at the start of my program and work from there. Do this has also greatly increased the speed of my program.


title: How to get Leaflet interactivity from Geoserver WMS?

text: I have geoserver serving shapefiles as a WMS into Leaflet. I'd like to add interativity that shows a popup with WMS attributes in the Leaflet map when the user hovers over or clicks on a WMS feature. I've seen the wax/leaflet interactivity, but is that possible with WMS and Leaflet?


title: How to add a Google Maps layer to Openscales?

text: I'm using Openscales in a web mapping project. I would like to add a Google Maps layer, but I haven't found any snippets or examples. Any solution will be very helpful.


title: Trigger on SQL Server 2008 ArcSDE 10 and ArcGIS Server 10

text: I need to build a sql server db trigger to update a date field with the current date when a record is inserted or updated, I think. I will be using the date field as a way to recognize updates/inserts programmatically. I went to sql server mngt studio and navigated to the database in the tree view >> opened up the programmability folder >> opened up the Trigger folder and right clicked >> Create new trigger >> I then copied the below sql into the window and hit F5. It messaged back to me that the command was successful but when I tried to open the Triggers folder I didn't see it there. When I tested in arcmap inserting a record I did not get the date updated to my date field. What am I doing wrong and are there special considerations for creating Insert/Update/Delete triggers against SDE? Thanks...the sql is below.               >CREATE TRIGGER t_I_TestTable     >ON dbo.TestTable     >FOR INSERT     >AS     >     >UPDATE TestTable SET TestTable.LastUpdated = Getdate()


title: Raster Calculator

text: I have a slope dataset. How do I calculate the real surface area for this slope dataset using the raster calculator. This is the example I was given: orthogonal area of a raster cell in square meters over cos([slope raster] * 3.1415/180) How do I know what the orthogonal area of the raster cell in square meter is? Where do I check for this?


title: Convert raster_stats list output to array / csv file

text: I usually program in R, but recently I have started to use python for most of my spatial programming. My python skills are still limited, and please excuse the simplicity of the question. I am struggling to convert a list, which is an output of a spatial function, into an array / csv file. In particular, I struggle how to disaggregate the list output that i receive from the spatial function. Here is the spatial function I have used:               from rasterstats import raster_stats     import time          os.chdir("C:/Users/Documents/veg_dt/data")          start = time.time()     stats_n_w = raster_stats("counties.shp", "veg_ndv.tif",stats=['count'], categorical=True)     end = time.time()     print end - start      The function raster_stats counts the occurrence of a categorical variable in a set of polygons. The output of the function is a list, and this is the actual output:               {0: 31728, 1: 35, 2: 31, 3: 30, 4: 29, 5: 40, 6: 44, 7: 82, 8: 85, 9: 142, 10: 20, 11: 58, 12: 76, 'count': 32400, '__fid__': 39989}     {0: 32204, 1: 6, 2: 13, 3: 5, 4: 15, 5: 27, 6: 33, 7: 5, 8: 28, 9: 28, 11: 3, 12: 33, 'count': 32400, '__fid__': 39990}     {0: 31983, 1: 38, 2: 11, 3: 21, 4: 42, 5: 17, 6: 19, 7: 52, 8: 66, 9: 34, 10: 9, 11: 87, 12: 21, 'count': 32400, '__fid__': 39991}     {0: 32400, 'count': 32400, '__fid__': 39992}      Every element of the list counts the categorical variable in a given polygon. The first row shows 31,728 zeros, 35 ones, 31 twos, 30 threes, 29 fours, 40 fives, etc in polygon 39989. The second element has 32,204 zeros, six ones, etc in polygon 39990. I would like to convert the list into a data frame / array / eventually csv file that states the fid (id of the polygon cells) in the first column, then lists the occurrences of zero in the second column, lists the occurrences of ones in the third column, etc. I want to have every element of the list on one row. Also, to complicate things further, the list is around 60,000 elements and I will do the same operation approximately 50 times for different sets of polygon and tiff images. Can anyone help with this?


title: while reading a tiff file filenotfoundexception

text: I am unable to create File object while instantiate a tiff file using TestData.file function it is giving file not found exception . But the tif file is exist on the same path > file exists > > Exception in thread "main" java.io.FileNotFoundException: Can not locate > test-data for "images/no_crs.tif" at > org.geotools.test.TestData.url(TestData.java:306) at > org.geotools.test.TestData.file(TestData.java:328) at > org.geotools.gce.geotiff.GeoTiffTest.main(GeoTiffTest.java:17) here is my code               public class GeoTiffTest      {         public static void main(String[] args) throws FileNotFoundException, IOException          {         File file = new File("images/no_crs.tif");         if(file.exists())                 System.out.println("file exists");         else                 System.out.println("file not exist");             final File noCrs = TestData.file(GeoTiffTest.class, "images/no_crs.tif");         }     }


title: QGIS-server and mapproxy. Management of the transparency of a layer

text: I have problems using mapproxy with QGIS-server via openlayers (mapproxy-1.6.0, qgis-server 2 and openlayers-2.13.1). I have a layer in which i defined the background of the polygon with transparency. When I asked directly qgis-server via openlayers, I get the png wanted (png background is entirely transparent and the transparency of the background is what i need). When I pass through maproxy, the image background is transparent but not the background of the polygons (the setting of the mapproxy may be ok because I get a transparent png image). What's weird is when I set the background of a polygon with less than 50% transparency in QGIS, through mapproxy the background of the polygon is completely opaque and when I put a transparency equal to or greater than 50%, the background of polygons is completely transparent. Thank you. Exemple directly from qgis-server : ![enter image description here](http://i.stack.imgur.com/f0TBz.png) Exemple via mapproxy : ![enter image description here](http://i.stack.imgur.com/QTC1M.png)


title: WMS render geotiff as a small image

text: Hye! I am facing Strange problem,My Geoserver 2.2.4 was stop working suddenly so i made a backup of Dir folder and install Geoserver 2.3.5 and replaced workspace folder with new one. Now all shapefiles is rendering correctly but Geotiff image is not render accurately it render as a small image on zero zoom level. I am unable to understand the issue because it was working fine with previous version I am using openlayers tu call wms i.e geotiff = new OpenLayers.Layer.WMS("geotiff image",   url + "wms?", {   LAYERS: "Sindh:" + featureType + "_geotiff",   transparent: true }, { isBaseLayer: false, singleTile: true } ); map.addLayers([geotiff]); I thought it was projection issue my Geotiff image have EPSG:32642 projection and i change its declare projection into EPSG:4326 but stilll issue is not resolved Please Guide me as soon as possible


title: ArcGIS javascript API Printing: fuzzy result

text: I created a real simple web app, copied from samples actually, to try figuring out something about the Print Dijit. Here's the thing. When you click print here there are two major concers with the result:   1. The symbology, especially the rectangle in the bottom left corner, are mashed up.   2. The background layer is not that high quality as it is in the browser Any thoughts how this is possible and how to resolve both issues?


title: Overpass XML query for finding all buildings inside a polygon

text: I want to run following Overpass QL query in order to get way IDs of all buildings inside a polyon:               way(poly:"50.7 7.1 50.7 7.2 50.75 7.15")[building];     >;     out ids;      Since I'm using the osm-common library, I need to convert this script to Overpass XML. I did it using the convert form here. The result:               <osm-script>         <query into="_" type="way">             <polygon-query bounds="50.7 7.1 50.7 7.2 50.75 7.15" into="_"/>             <has-kv k="building" modv="" v=""/>         </query>         <recurse from="_" into="_" type="down"/>         <print from="_" limit="" mode="ids_only" order="id"/>     </osm-script>      But when I run it, I get node IDs, not way IDs.               <?xml version="1.0" encoding="UTF-8"?>     <osm version="0.6" generator="Overpass API">     <note>The data included in this document is from www.openstreetmap.org. The data is made available under ODbL.</note>     <meta osm_base="2013-11-06T07:51:02Z"/>            <node id="252553109"/>       <node id="252553111"/>       <node id="252553112"/>       <node id="252553114"/>       <node id="252553115"/>       <node id="252553117"/>       <node id="255302992"/>       <node id="255302993"/>       <node id="255302995"/>       <node id="255763355"/>       <node id="255763356"/>       <node id="255763357"/>       <node id="255763358"/>       <node id="255763430"/>       <node id="255763431"/>       <node id="255763432"/>     ...       <node id="2519505892"/>       <node id="2519505893"/>       <node id="2519505894"/>       <node id="2519505895"/>       <node id="2519505896"/>          </osm>      How should I modify the Overpass XML query in order to get way IDs of the buildings?


title: OpenLayers - redrawing map after container resize

text:
In my web application, I want to enable users to set the size of the map container. Everything worked fine when the container was expanded slightly (apparently this is because the tiles that were just behind the border were already loaded). However, when the container is expanded significantly (in the following example, from 300 to 1000px in width), there is blank space left. How to make the map redraw and adapt to the new size? Calling `redraw()` on all layers didn't help. Neither did zooming in and out. I tested this with the described results in Opera, Chrome and Firefox. In IE8, surprisingly, the problem didn't occur and the map automatically adapted. A simplified webpage for testing:

         <input type="button" id="test" value="resize">


title: Simulate (Mock) GPS Tracks

text: I have the following scenario: 1) I need to create GPX files for use in testing of a mobile device. 2) The tracks to be described within will span large distances over long time periods (i.e. it's not practical to just go out and record a track). 3) The people performing the testing are not very GIS knowledgeable and will struggle with things like manual file conversion, etc. I have in mind some simple application (it could easily be a QGIS python plugin) in which you set certain parameters (time, direction & speed of track "stages") and then generate a GPX file accordingly. But before I throw myself into something like that, does anyone know of a simpler way to do this? I've looked into various options (creating a tour in kml and then converting into gpx, manually cut-and-paste the xml together by hand, etc.) but it all seems like a lot of hard work for what I would have thought was quite a common use case...


title: How to install OpenLayers plugin in QGIS1.8

text: I am using QGIS 1.8.0 on Windows 7 (32-bit) but am unable to install the OpenLayers plugin. The Python plugins menu does not include OpenLayers on the plugin list. Can it be somehow added manually?


title: tools from Geospatial Modelling Environment in ArcMap model

text: How to extract raster DEM elevations to vector attribute tables in ArcGIS 10? Can GME tools be used in model for publishing a tool on web? I want to use it in integration with other tools in ArcMap. Thanks


title: Openlayers shapefile polygons hover effect

text: I have shapefile polygons. I'm using OpenLayers with OpenStreetMap tiles. Is it possible to create hover effect which will show up info window with some informations about that polygons? Is it the best to put it in some database maybe? Thank you.


title: How to Display Point from Form

text:
For the last few days I have been trying to come up with a way to display a point on a map from a coordinate pair entered into a form. I began with the sample below and then added the formToJSON Dijit. Both components seem to be working but I can't figure out how to "hook" them together. Does any one have any ideas or am I going about this all wrong? Links to the original sample, my edits and code are below. Original Sample: http://help.arcgis.com/en/webapi/javascript/arcgis/demos/map/map_topo_graphics.html dojo.formToJson http://dojotoolkit.org/reference-guide/dojo/formToJson.html My Edits: http://geoville.org/viewers/DP/               <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN" "http://www.w3.org/TR/html4/strict.dtd">                     <meta http-equiv="Content-Type" content="text/html; charset=utf-8">                  <meta name="viewport" content="initial-scale=1, maximum-scale=1,user-scalable=no"/>         <title> Display Point</title>         <link rel="stylesheet" type="text/css" href="http://serverapi.arcgisonline.com/jsapi/arcgis/2.1/js/dojo/dijit/themes/claro/claro.css">

         <form id="myform2">         <input type="text" dojoType="dijit.form.TextBox" name="X" value="-13043577">         <input type="text" dojoType="dijit.form.TextBox" name="Y" value="4034040">         <button dojoType="dijit.form.Button" onclick="drawPoint()" id="convertFormDigits"> Create Point </button>         </form>

                  Thanks for your help. -Mike


title: PostGIS: split a polygon in a given ratio

text: Is there a convenient way to split a polygon into `n` parts, more-or-less equal in size in PostGIS?


title: automatic snap line to line

text: I need to automaticly overlay/snap the red line (gps-track) exact on the lines of my roads-shape. Is that possible? ![enter image description here](http://i.stack.imgur.com/i5Nw5.png)


title: Alrign Track from GPX file with roads

text: I'm looking into possibilities to align track data from a GPX file with known roads, perhaps coming from OSM. So far I'm able to show the tack as such on a map powered by leaflet. Since the app should be as easy to use as anyhow possible I look for web based and if possible fully automated solutions. The GPX files are stored in a standard CouchDB, so is the App itself. I think about an approach in which users can upload their GPX file. Once they do the map goes into an edit mode that shows the track and may already suggest a "road only" version of the track. In a best case scenario all the user has to do is to accept the suggestion. How could something like that be done? Are there any libraries or plugins doing that?


title: Making Map Matching Points to Roads

text: I'm having some trouble while map matching these points for some streets. Could you lend me a hand? ![enter image description here](http://i.stack.imgur.com/85787.png)


title: How to implement a "Snap to Road"-functionality?

text: I am having a problem with a little application I'm trying to develop. What I do is take a phone with me and drive around the city on a bike. I collect all the Latitude/Longitude-pairs and at home I plot them on a map. What I get, however, is zig-zag lines, because the GPS location is not perfectly accurate. Here's an example plot: ![Location data visualized](http://i.stack.imgur.com/04Z4v.png) Obviously I was driving in a more or less straight line and not through the hotel. ;) This is perhaps basic for GIS-experts, but here's my question: How do I get better results, with the red line (and thus the recorded points) being on the roads? These are the things I've considered:   * Rolling-average over measurements. E.g. take 5 coordinates and calculate average coordinate. Tried it and it doesn't work. It didn't work with the data shown in the picture. (But the line was smoother..)   * Only accepting measurements where the device (iPhone) reports a high accurracy. Tried this and it helps a little bit, but still not good enough.   * Correlating of GPS data with acceleration and direction sensors of the phone. I haven't tried this yet and it sounds non-trivial.   * Using some sort of "Snap to Road"-algorithm. I think it should be possible to get data on roads for a given location. (OSM?) Then perhaps the nearest road and the nearest point on that road from the current position could be calculated and the position could be adjusted to be ON the road. Sort of like the "snapping to grid" function of graphical interfaces. If this is the right approach, what are some pointers for software and methods I can use? I assume I will have to deal with the following problems:   * Huge dataset for any given country or city (I'd have to slice tiny portions out, but also store the whole dataset on my own server?)   * Transformation of formats..   * High CPU/memory-usage(?)   * A complex algorithm(?) Looking at car navigation software it seems that this problem has long been solved. Are there any implementations that can be re-used, preferably of course for the iOS? Thank you **UPDATE:** I have implemented a simple map matching algorithm:   1. I downloaded several hundred MB of road data from OpenStreetMap (OSM) and inserted them into my DB.   2. For each measured GPS-point I fetched all OSM-nodes in a certain radius.   3. For each fetched node I fetched all roads (series of nodes) from DB and treated them like line segments in 2d-vectorspace.   4. I calculated the closest point (distance from gps-measurement) for each segment and selected the one with the shortest distance.   5. I've changed the lat-lon coordinates of the measurement to those of the matched point. Here is the actual result for the example plot above: ![example plot with matched coordinates](http://i.stack.imgur.com/H6CbC.png) It's just slightly better. I think the next step is to calculate _multiple_ points in step 4 and use each of them in a probability / plausibility step 4a, where each one is checked for things like wether the bike really took a turn onto another street. It's not real-time so it should be possible to take the future into account (e.g. to look at the next measurement that followed, or the next few.), and to calculate a probability score for each candidate point. Anyway, just wanted to post a quick update. Perhaps this will help someone.


title: Find distance and heading from a point to a polygon in lat, lon coordinates

text: I have counties mapped using polygons in Shapely with latitude/longitude for coordinates. Given a complex polygon x and single point y, how would I find the shortest distance between x and y as well as the heading needed? Using the distance function in Shapely, I can get the closest distance, but I'm unable to come up with a fast solution to get the heading. I tried determining the closest point, but this solution was unreasonably slow to run on 10 polygons.


title: Passing GPString value via Flex

text: In this example from the ArcGIS for Flex samples page: Clip and Ship I am curious how from this function how the label property of each check box that is passed as input for the Layers_to_Clip parameter of the GPService are interpreted by the GPService as actual pieces of data.               private function getLayersToClip():Array             {                 var layersToClip:Array = [];                      if (incidentPointsCheckBox.selected)                 {                     layersToClip.push(incidentPointsCheckBox.label);                 }                 if (incidentLinesCheckBox.selected)                 {                     layersToClip.push(incidentLinesCheckBox.label);                 }                 if (incidentAreaCheckBox.selected)                 {                     layersToClip.push(incidentAreaCheckBox.label);                 }                      return layersToClip;             }      I have looked at the Clip and Ship data from the ArcTutor files and don't understand how the GPString values - which in this case are "Input Incident Points", "Input Incident Polylines", and "Input Incident Polygons" - correspond to any of the map layers. Part 2) Ideally I would like to change my Flex application to take items selected from checkboxes and pass them into my model in the same fashion, except I won't be extracting them but using them in analysis. If preprocessing is done, every input would be a raster dataset regardless of the original geometry and the output would be a raster as well. How can I have those checkboxes refer to selected input rasters in a database location or within the map document/gp service? Thanks in advance, I've been banging my head against the wall with this one, Sean