title: What purpose do service definition packages serve?

text: I'm a software developer. In my world, when you bundle up a set of scripts or binaries to be run on another machine, your goal is usually to publish it to multiple environments. Sometimes the intention is that someone else will use your scripts or binaries inside of an application (libraries), and sometimes you're actually trying to publish an application. In the former case, you generally do all you can to make your tool as environment agnostic as possible, allowing users to pass in information about what resources to use (files, databases, web services). In the latter, you give the user some ability to configure the application so that it can connect to the correct resources. Making an application configurable like this still has uses even when you're only giving it to one client. In particular, it makes it simpler to have multiple testing copies of the application before you actually push to production. (E.g., one for developers to screw up without bothering anybody, one for internal QA, and of course production.) When I first heard about service definition files (.sd files) for geoprocessing services, that's what I thought they were supposed to do: let you bundle your code up so you could deploy it to different environments. However, I have found that when you use them, there are zero options for configuring anything at publish time. (See here and here.) This means that a service definition is locked to a single environment, and you have to create a new one for each environment. As a result, I cannot see any way in which service definitions are useful for geoprocessing services since you must recreate them in each environment you want to deploy to. So I pose the question: what is their intended purpose? What problem do they solve? When should I be using them?


title: Modifying properties of Attributes Table in QGIS?

text: In my QGIS project I'm using vector layers that require a lot of data inputting/editing. The problem is that some of the data I'm inputting varies greatly in matter of length and structure (various number of lines of text for example). When it comes to edit those lines in many objects or to compare the data in various objects it is painful to do it by selecting each object(line or point) separately. Instead of doing all this work by selecting each object and edit its properties I'm using Attributes Table. I would like to know is there a way so Attributes Table remembered height and width of each column? Is it possible to modify those properties in python? And last but not least, do you think there is reasonable effort in writing plugin that extends Attributes Table for functionalities mentioned above.


title: How to add point with current mouse position while mouse over on line feature?

text: I have a line feature in a vector layer "A", while i mouse over on that line. A point should be added to vector layer "B" with current mouse position. OR how can get lat lon of mouse position while mouse over on that line? Thanks in advance... :)


title: Store "template" mxds inside model

text: Since the ability to store mxds inside a file gdb is still just an idea, I've been mulling over my options for using a "template" mxd within a model. My model uses 4 different mxds, but currently it accepts one at a time (i.e. I run the model 4 times for each feature class, using the different mxds). The model does various things such as saving a copy of the "template" as a new mxd, replacing data sources etc. Up to now I have "Input Template mxd" as a model parameter. This is fine, if I'm the one running the model, so I know which mxd is the correct one. However, if someone else were to run the model, they could input any mxd they wished, which would inevitably lead to the model failing. The options I thought of were:               1. Trust that the user will enter the correct template mxd containing         the correct layers: if only.     2. Force the user to enter a template mxd which has the correct name:         won't work in this multiple template mxds scenario, as I would then have to         duplicate the model 4 times to check for the 4 different names. Alternatively,        if I wrote an if statement inside the model to check which template I'm receiving,        it still would not potentially stop the user from renaming Template1.mxd to         Template2.mxd, causing the wrong processes to run.     3. Store the templates inside the same file gdb with the model: as stated,         cannot be done.      My current scenario has my model and template feature classes inside a gdb, inside a folder along with the four template mxds. The model automatically accesses each template mxd as needed, so the branch for each template runs consecutively everytime the model is executed. However, the user could still rename the mxds, or simply remove them from the folder... My question is, is there a way to somehow store the template mxds within the model itself or in some other way so that the user cannot change it?


title: TileMill shows black lines around raster layers

text: I've taken a large number of raster images and added them to TileMill. Everything has worked great. The problem that I am seeing is that TileMill adds a small black line around each layer. I know that the layers overlap so there shouldn't be any gaps in between the layers. Is there an easy way to get rid of the black lines? These lines do not show when viewing the raster layers in QGIS. ![enter image description here](http://i.stack.imgur.com/ocscU.png)


title: Create temporary files as intermediate output to be used as input in concatenation

text: Q_PART_I I need to merge many (MANY) tables in one unique output table by means of the arcpy.Merge_management. Anyway, I'd like to avoid the creation of tons of intermediate files as my script running, having only my final big table rather than a number of "mini-table" for each file I have. I'm converting raster attribute table to table views (by arcpy.TableToTable_conversion), and then using them to run the Merge. How can I avoid having a table for each raster in the intermediate process? Hope someone could help me. I've been reading something about the creation of an InMemoryWorkspace, but I couln'd manage it so far. Thanks in advance. Q_PART_II Ok, I tried both Your suggestions, but so far I couldn't make it. It seems the Merge function cannot recognize the temporary files as input. Infact, I receive this error: ERROR 000735: Input Datasets: Value is required Failed to execute (Merge) Here is my script (the important part is at the end, I can always avoiding this problem by simply having a lot of files, but I'd prefer to be more "clean"):               #import relevant modules, create geoprocessing dispatch object     import win32com.client, sys, string, os, arcpy          gp = win32com.client.Dispatch("esriGeoprocessing.gpDispatch.1")          # Remember to change this to wherever your files are stored     gp.workspace = 'G:\\Chile_2012\\MODIS10A2.V005\\provvisoria\\ATTRIBUTE_TABLE'     # Remember to change the path to point where you want to store the final tables     out_dir = 'G:\\Chile_2012\\MODIS10A2.V005\\provvisoria\\ATTRIBUTE_TABLE'          try:         rasters = gp.ListRasters("*", "all")         rasters.reset()         rst = rasters.Next()              while rst:             # Create the new field             gp.AddField_management (rst, "FILENAME", "text", "", "", "50")                  # Apply the filename to all entries                    gp.CalculateField_management (rst, "FILENAME", '"' + rst + '"')             rst = rasters.Next()          except:         print gp.GetMessages ()          listRaster = arcpy.ListRasters ()          # Create a table view from each raster to be used as input for the next Merge function:     for raster in listRaster:         arcpy.TableToTable_conversion (raster, "in_memory", 'Tab_'+ raster)          listTable = arcpy.ListTables ()          # Change the name of the output table as needed!!!     arcpy.Merge_management (listTable, 'Cl_Rio_Aconcagua.dbf')      Hope this could help You helping me in turn. Many thanks in advance!


title: Python script for GRASS - executing from PHP

text: I wrote the following python script:               #!/usr/bin/env python     import os     import sys     gisbase = os.environ['GISBASE'] = "/usr/lib/grass64"     gisdbase=os.path.join("PATH")     location=LOCATION NAME     mapset=MAPSETNAME     sys.path.append(os.path.join(os.environ['GISBASE'], "etc", "python"))     import grass.script as grass     import grass.script.setup as gsetup     import grass.script as grass          gsetup.init(gisbase,                 gisdbase, location, mapset)     print "one"      grass.run_command('v.in.ascii','-o', input='/opt/lampp/htdocs/project/basic_statistics/sites_grass.csv', output='sites',fs=",",skip='1', x='9', y='10')     print "two"     grass.run_command('v.db.addcol', map='sites', columns="height INTEGER")     grass.run_command('g.region', rast="dem21")     grass.run_command('v.what.rast',vector='sites', raster='dem21', column='height')     print "three"     grass.run_command('v.out.ascii','-o', input="sites",output="sites_geography_stats.csv", columns='height')     print "four"      It runs without any problems when I simply execute it from UNIX terminal (/.FILENAME.py). However, my goal is to execute it from a php code:               <?php     ...     echo system("./FILENAME.py");     ...     ?>      When I run the php code the following output is printed to the website: one two three four Sadly, the sites_geography_stats.csv file is not created. It looks that the python script is executed from php but its GRASS-related part is ignored. Does anybody have any ideas how to fix it?


title: PopUps for all the Feature layers

text: I am trying to create GeoExt action button to fetch the feature info from all the layers. I am still getting the info for my point features but I am getting nothing for the usa:states, that appears in the geoserver by default. When I looked into the firebug, the html for states does not have anything in it's , http://dpaste.de/NLjzW/. But the species data that I am adding has the contents in the , http://dpaste.de/XR3Ou/. The contents in the body are the result of passing `evt.text` param, which is assigned to the `html: evt.text` in my `addToPopup()` function. Why is this happening? **Update:** I am sorry i forgot to add the addToPopup() function before. I thought that maybe there was no data in the usa:states layer, but it does have the data. If you click on the polygon in this app, `http://128.196.142.12:8080/geoserver/usa/wms?service=WMS&version=1.1.0&request=GetMap&layers=usa:states&styles=&bbox=-124.731,24.956,-66.97,49.372&width=780&height=330&srs=EPSG:4326&format=application/openlayers` You can see the data in the polygon.


title: Cannot merge dbf tables (problem with lock files???)

text: I am new to Python and to this community as well, hope someone could help me. I am running script to extract and merge toghether the tables of some Raster. Each Attribute Table presents only one row and the same field names. My script works fine until the very end of it, when I try to merge all the tables. It seems the tables are locked (I can see an "sr.lock" file for each Table I want to merge in my folder (that desappear only when arcgis programs are closed) and I cannot access them. Strange thing: I do access them if I simply use a sample (for instance) of two of these so created tables in the Python window within arcMap and type their names instead of using a FOR LOOPING... Here is my script: #import relevant modules, create geoprocessing dispatch object import win32com.client, sys, string, os, arcpy               gp = win32com.client.Dispatch("esriGeoprocessing.gpDispatch.1")          # Remember to change this to wherever your files are stored     gp.workspace = "G:\\Chile_2012\\MODIS10A2.V005\\provvisoria\\ATTRIBUTE_TABLE"     out_dir = "G:\\Chile_2012\\MODIS10A2.V005\\provvisoria\\ATTRIBUTE_TABLE"          try:         rasters = gp.ListRasters("*", "all")         rasters.reset()         rst = rasters.Next()              while rst:             # Create the new field             gp.AddField_management (rst, "FILENAME", "text", "", "", "50")                  # Apply the filename to all entries                    gp.CalculateField_management (rst, "FILENAME", '"' + rst + '"')             rst = rasters.Next()          except:         print gp.GetMessages ()          listRaster = arcpy.ListRasters ()          for raster in listRaster:         arcpy.TableToTable_conversion (raster, out_dir, raster+ '.dbf')          listTable = arcpy.ListTables ()          for table in listTable:         arcpy.Merge_management (table, 'Cl_Rio_Aconcagua.dbf')      The code is taken from others, I simply adapted it for my needs. Anyway, the error I am receiving so far doesn't have to do explicitly with LOCK topics, but since I have those files in my folder the problem might be due to that. Even if I don't use the FOR LOOP (which seems not necessary), I always receive this error:               ERROR 000732: Input Datasets: Dataset 'xxx;xxx...' does not exist or is not supported     Failed to execute (Merge)


title: Export .dbf maintaining column data type

text: I require a shapefile's .dbf to contain two types of text columns VARCHAR(n) and CHAR(n) to meet specs (attribute format) of gov receiving body. Have to send them the .shp, .dbf and .prj. Within arcmap or qgis add field only allows text and length for new fields. In postgres i can create the required format for example filed blah data type is character (10) and field blah2 data type is character varying(80), however if i export this as .dbf and then reimport it (to check data types are maintained) using the postgis plugin it changes the character(10) to character varying (10). Does anyone know a way to create a .dbf and export it maintaining set data types, i.e character (10) and character varying (80) for example. Thanks


title: Most optimized way of extracting data from planet.osm

text: I want to cut extracts from planet.osm.bz2 . I have found this command               bzcat downloaded.osm.bz2 | osmosis  --read-xml enableDateParsing=no file=-  --bounding-box top=49.5138 left=10.9351 bottom=49.3866 right=11.201 --write-xml file=- | bzip2 > extracted.osm.bz2      Only concern I have the disc space. planet.osm uncompressed is 400GB! How much space would use this command from WIKI?


title: When to use multiple geodatabases?

text: I'm working on a project involving one database , but the customer's needs have evolved, so he needs to manage multiple projects. I tried to work with adding an id for each project, which gave me more complex queries all over the application, and I need to perform a replicate per project , but the the replicates would be heavy for a one database containing multiple projects. In this application I'm using multiple arcpy geoprocessing scripts, so if I will opt for multiple databases the SDE database would be a parameter. But multiple databases would be easier to backup and to archive. I tried to make a picture about the context of the question: **Question:** Is there any risk in adopting multiple geodatabases in one application? Any advice or experiences are welcome.


title: scale dependent rendering of base layers(gmap,ghyb) causes warping/flickering

text: Hey guys I am using geoserver to host my WMS layers and I am using gmaps base layers in Openlayers. Below I have attempted to use scale dependent rendering of the base layer without any luck. What I am trying to do is once the user zooms past a certain level, Google's hybrid map will turn off and Google's aerial will turn on.....With the code I have now it will turn on and off but it tends to warp the base map or flicker. Not sure whats going on.                           /*map.events.on({ "zoomend": function (e) {                 if (this.getZoom() > 16) {                   gmap.setVisibility(false);                   ghyb.setVisibility(true);                 }                 else {                   gmap.setVisibility(false);                   ghyb.setVisibility(true);                 }               }             });*/


title: Matching points in table to polygons

text: this is my first post to this forum. I hope to start participating actively. I have a table with several records (aprox 200.000), each with its own coordinates. And I need to identify the "Polygon" where each point belongs (I got another table with the vertex for earch polygon) . I am working with Excel and Access at this point. Thanks for the help! (the polygons doesn’t have holes and doesn’t overlap with other polygon) I know there is an algorithm that checks if the number of times a line cross another line from the desired point to another point is even, the point is inside polygon, if the number is odd, the point is outside the polygon. But I don’t know how to implement this. Here are some reference links with the algorithm: http://alienryderflex.com/polygon/ Here is a sample: I got a table similar to this:   * Point_ID / Coordinates    * 1 / -70.559349,-33.411913   * 2 / -70.579004,-33.407542   * etc. (aprox. 200.000 points) And I got the following table defining the boundaries of each polygon:   * Polygon_ID / Polygon_Sequence / Coordinates   * 1 / 1 / -70.567331,-33.408803   * 1 / 2 / -70.584755,-33.414249   * 1 / 3 / -70.586987,-33.409692    * 1 / 4 / -70.585527,-33.404175   * 1 / 5 / -70.571537,-33.398084   * 1 / 6 / -70.567331,-33.408803   * 2 / 1 / -70.567245,-33.409047   * 2 / 2 / -70.564756,-33.416570   * 2 / 3 / -70.552998,-33.418576   * 2 / 4 / -70.551796,-33.408474   * 2 / 5 / -70.560122,-33.407542   * 2 / 6 / -70.562525,-33.407327   * 2 / 7 / -70.567245,-33.409047   * etc... So, this is the output I need (in **bold** ):   * Point_ID / Coordinates / **Polygon_ID**   * 1 / -70.559349,-33.411913 / **2**   * 2 / -70.579004,-33.407542 / **1**   * etc... Thanks a lot! Alejandro


title: Object does not exist in the namespace "http://schemas.esri.com/arcgis/client/2009"

text: I am trying to move all of my GIS code over the Visual Studio 2012. I have added the following references: C#:               using ESRI.ArcGIS.Client;     using ESRI.ArcGIS.Client.Geometry;     using ESRI.ArcGIS.Client.Tasks;      XAML:               xmlns:esri="http://schemas.esri.com/arcgis/client/2009"      The problem is that many of the objects are not being recognized because they do not exist in the namespace... Some of these include the InfoWindow and Legend. However it seems to be OK with objects such as the PictureMarkerSymbol and the Map and Layers. I don't know if something has changed since I last looked at it, if so I think ESRI need to update their SDK samples. Any help would be appreciated. Thank you in advance.


title: using ZonalStatisticsAsTable in python script

text: I use a python script in which the function ZonalStatisticsAsTable is called. In ArcGIS9.3 it all worked fine. Our company just went to ArcGIS 10.1 and now it doesn't work anymore. When looking inside the resulting temporary tables it seemed the name of column number 2 was changed from "VALUE" to "NR". I changed it in the script and it all seemed to work fine. But, the results are very different from these in 9.3. As input I use a rastermap with cellsize of 1x1km and a Polygon feature. I guess the polygon is converted to a raster first before the comparison is made. Whenever I first use "polygon to raster" to created a raster and use "Zonal Statistics as Table" with these two rasters as input the results are fine.   What goes wrong using then function ZonalStatisticsAsTable in the Pythonscript? Regards, Wilco.


title: How to write Process Dialog log information to text file in ArcGIS?

text: I have created a custom model. When this model is run, it shows a process dialog which has some log information. I want to write that log information into the text file for further analysis. Does anyone know how can I do this?


title: Autocorrelated DEM with QGIS or Grass

text: I was wondering if there is a way to make an Autocorrelated DEM within QGIS or Grass. I am looking to do this with stereo aerial photography. I am able to do this with a number of COTS tools, but I am looking for an open source solution. If you know of other tools I am open to suggestions.


title: Require a 3rd Party GIS service to embed in my application

text: I am looking for a service (Free or Premium) that allows offline access to maps and also gives the functionality to search map using an address or Latitude / Longitude. I want to embed it into a Web Application. I am particularly interested in Maps. Please tell me if any one know of any such service. Thanks. Note: Offline means the map's data will not reside on the Internet but on a private (intranet) server.


title: ArcGIS 10 tool (add x,y) is not working

text: I have a problem with go to X, Y tool in arcgis 10. When I want to use this tool the program crashes. Is there any service pack or something to repair that problem?


title: Making 3D landscape area appear in Vectorworks 2011

text: I am using Vectorworks 2011 to build a site plan. I would like to include a 3D view using a "Landscape Area" to display a block of vegetation. I have drawn the shape of the area, set the 3D view to "3D Plants" and added the plant information. When I turn on OpenGL mode and switch to a 3D view the plants appear only as plain grey cross-shaped 3d objects (not sure of the right word!). I can't get them to look like plants either as 3D shapes or 2D images.


title: Coordinate system in Map&Guide

text: One of my customers is evaluating transport planning software Map&Guide (http://www.mapandguide.com/) and what we require is to get the transport route data from their file to Garmin navigation device. The problem is that the M&G file contains coordinates in some system that I do not recognize and I would like to ask your help to sort it out. Slice from the M&G *.BCR file:                    [COORDINATES]     STATION1=2751856,8265742     STATION2=4506233,9477610     STATION3=4183918,7501674     STATION4=3393742,6515209     STATION5=5539499,4914891     STATION6=4983102,5109208     STATION7=3655214,4850900     STATION8=1503212,6863049     STATION9=-49653,6695422     [DESCRIPTION]     STATION1=EST 10148,At Tallinn,Pärnu Maantee,0,,0,0,0.000,0.000,0.000,0.000,0.000,0.000,0,,0,0.0000000000,0.00,0.00,,0,0,0,0     STATION2=RUS,At Archangel'sk,,0,,0,0,0.000,0.000,0.000,0.000,0.000,0.000,0,,0,0.0000000000,0.00,0.00,,0,0,0,0     STATION3=RUS ,At Moskva / Kitaj-Gorod,Bol'shoj Cherkasskij pereulok,0,,0,0,0.000,0.000,0.000,0.000,0.000,0.000,0,,0,0.0000000000,0.00,0.00,,0,0,0,0     STATION4=UA 01001,At Kyjiv,Esplanadna Vulytsia,0,,0,0,0.000,0.000,0.000,0.000,0.000,0.000,0,,0,0.0000000000,0.00,0.00,,0,0,0,0     STATION5=AZ ,At Baki,Nanman Narimanov Avenue,0,,0,0,0.000,0.000,0.000,0.000,0.000,0.000,0,,0,0.0000000000,0.00,0.00,,0,0,0,0     STATION6=GE ,At T'bilis? / Avlabari,Ketevan Camebulis Kuca,0,,0,0,0.000,0.000,0.000,0.000,0.000,0.000,0,,0,0.0000000000,0.00,0.00,,0,0,0,0     STATION7=TR 31600,Ankara / Dörtyol,,0,,0,0,0.000,0.000,0.000,0.000,0.000,0.000,0,,0,0.0000000000,0.00,0.00,,0,0,0,0     STATION8=D 12529,Berlin-Schönefeld-Terminal a,Flughafen,0,,0,0,0.000,0.000,0.000,0.000,0.000,0.000,0,,0,0.0000000000,0.00,0.00,,0,0,0,0     STATION9=GB TW6 1,London Heathrow-Terminal 4 Arrival,Salisbury Road,0,,0,0,0.000,0.000,0.000,0.000,0.000,0.000,0,,0,0.0000000000,0.00,0.00,,0,0,0,0           The coordinates should match the locations in second section. Being from Estonia, I do not recognize 2751856,8265742 pair as anywhere near me. Any help on the matter or any hints where to dig on are much appreciated


title: Why does ManageR and rpy2 not work in QGIS (1.8.0-Lisboa)?

text: I have been trying to get R scripts to function and ManageR to work within QGIS (1.8.0-Lisboa). (Windows 7). Python2.7.5. Two problems indicated that R was not working in QGis. First - when I run an R script there is not output in the pop-up window. Second - when I click on ManageR (bottom right corner of Qgis) it tells me I should either install rpy2 or check my installation of rpy2. So, installing (checking rpy2). This in itself has proved very problematic. I found the following posts and 'quick-fix' software: http://stackoverflow.com/questions/11928277/rpy2-and-python-installation-road- block https://bitbucket.org/breisfeld/rpy2_w32_fix/issue/1/binary-installer- for-win32 followed the instructions but still to no avail. Sadly I am no computer genius. I can follow simple instructions but am not up to fixing anything myself. So, guys. This is as much a request as a question. (a) how to fix, (b) is anybody working to fix this? many thanks BW


title: Symbology for Open Street Map [OSM] Data in MapInfo [tab] format

text: I have shapefiles for the different SOM layers which I have translated into TABs but I am now trying to symbolise them in a similar style to OSM online. Is there a way to easily apply the OSM symbology to MapInfo TAB files?


title: heat map of density ratio of 2 layers of points

text: I have 2 vector layers of points:   * one with a point for each house in the city (thanks to osm.org), roughly equivalent to population   * one with a point for each case of cholera in the city I'm trying to define in which locations I should site new water supply points based on a local attack rate (ie: number of cholera cases / local population). I've built separate heat maps for each layer (using QGIS heat map plugin), but it's hard to visually figure out hot spots. So I'm looking for a way to calculate this local attack rate, ideally something like (number of cases within X meters radius / number of houses within same radius) and spit that out as a raster layer. Basically just like a heat map, but using a ratio rather than a simple count. I'm guessing there's some magic keyword I don't know that will open a swath of solutions. Can anyone point me in the right direction? I'm using QGIS.


title: Using the Processing Framework in QGIS 2.0

text: I am writing a plugin for QGIS 2.0, and I want to use the processing framework to call some functions. When I open the pyhton console inside QGIS, I am able to call it like this:                   >> import processing      And then, for instance:                   >>> processing.alglist()      However, if I try to import it in my Python Plugin, it fails to find the processing include:                   from qgis.core import *         from qgis.gui import *         import processing      I googled the problem and found this directives:                   from processing.core.Processing import Processing         from processing.core.ProcessingConfig import Setting, ProcessingConfig      But unfortunately it also does not work for me :-/ Does anybody know how to include the qgis processing framework in a plugin, and how to call the algorithm (lets say for instance "processing.runalg(...) I really appreciate your answers, as I found many explanations on how-to do it in the Python console, but not outside QGIS :-/ Thanks again for your time and have a nice day, Joana


title: Efficient computation of path profile with Numpy and GDAL

text: I have numpy array with SRTM elevation data and I need calculate path profiles - vectors with elevation data from one point (transmitter) to every points of my map. Here's my current solution:               def get_elev_vector(self, coord1, coord2, step=90):         if coord1 == coord2:             return 0         lat1, long1 = coord1         lat2, long2 = coord2         dist = self.get_dist(coord1, coord2)         npt = dist / step         vect = np.zeros(npt + 3, dtype='float32')         vect[0] = npt         vect[1] = dist         vect[2] = self.get_elevation_direct(coord1)         vect[npt] = self.get_elevation_direct(coord2)         g = pyproj.Geod(ellps='WGS84')         coordvect = g.npts(long1, lat1, long2, lat2, npt)         i = 3         for lon, lat in coordvect:             vect[i] = self.get_elevation((lat, lon))             i += 1         return vect               #(tx, ty) transmission coordinates     while lat > latmax:              long = longmin              while long<longmax:                    vect =  self.carta.get_elev_vector((tx, ty), (lat, lon))                    arr[xarr, yarr] = 100                    yarr +=1                    lon += dlon              xarr +=1              lat +=dx      It works, but it take about 3 hours for 20x20 km area. Is there any way to make it work faster? Add                   def get_elevation(self, coord):         """TODO"""         lat, lon = self.xyfromlatlon(coord)         mlat = int(lat)         mlon = int(lon)         vicin = self.tile['numpy_array'][mlat: mlat + 2, mlon: mlon + 2]         temp_bool = np.equal(vicin, -32768)         np.putmask(vicin, temp_bool, 0)         elev = billin_interpol(vicin[0][0], vicin[0][1], vicin[1][0],                                             vicin[1][1], lat - mlat, lon - mlon)         return int(elev + 0.5)               def billin_interpol(tl, tr, bl, br, x, y):     """Bilinear interpolation.        Return:         float, interpolated value     """     b1 = tl     b2 = bl - tl     b3 = tr - tl     b4 = tl - bl - tr + br     return b1 + b2 * x + b3 * y + b4 * x * y


title: Calculating Z-scores in Arc

text: I am trying to take the natural log and then calculate z-scores on cell values in several composite band GRID rasters in Arc. I am not fluent in python and do not see a way to automate this within Arc - I am in need of some coding help or a suggestion on how to automate this within model-builder. Thank you for the help!


title: How to make a Python script that does a task in QGIS?

text: When I run this code in QGIS Python console it works fine, but I want want to run it from an external Python script. What I want the script to do is: 1- Open QGIS 2- Run this code 3- Close QGIS               from qgis.core import *     from qgis.gui import *     from qgis.utils import *     from PyQt4.QtCore import *     from PyQt4.QtGui import *     from PyQt4.QtXml import *          mapinstance = QgsMapLayerRegistry.instance()          uri =  "file:///C:/Users/Home/Desktop//sample.csv?delimiter=%s&xField=%s&yField=%s" % (",", "Long", "Lat")           vlayer=QgsVectorLayer(uri,"sample","delimitedtext")     mapinstance.addMapLayer(vlayer)          uri="C:/Users/Home/Desktop/Style.qml"     vlayer.loadNamedStyle(uri)               vlayer.selectAll()     qgis.utils.iface.actionZoomToSelected().trigger()     vlayer.invertSelection ()          qgis.utils.iface.mapCanvas().saveAsImage('/Users/HOME/Desktop/capture.png')


title: Why does python reads only the no data value from GeoTiff file?

text: I had a number of tiff files with a country's map. They had different rows, columns or different projections etc. I needed to make them with same projection, spatial resolution and same number of rows and columns. So I used `QGIS Raster->Projections-> Warp(Reproject)` method and a vector boundary shape file to get `tiff` files with the intended parameters. After doing this each file contained **two bands**. One band with the original data values and another band with values from 0 to 255. (I think it is for visualization, is it?) In some of these files I had specified no data values and in others, I didn't. But when I'm reading these files in `numpy` arrays in Python using `gdal`, I'm only getting an array of intended size `(53,54)` **but with no data values!** I used the following code:               from osgeo import gdal     from gdalconst import *     from scipy import stats     import numpy as np     import os, glob     import statsmodels.api as sm     import matplotlib.pyplot as plt          os.getcwd()     os.chdir(r"H:\...\Processed")     #load tiff files in numpy arrays     ds1 = gdal.Open("cases2001.tiff")     can01 = np.array(ds1.GetRasterBand(1).ReadAsArray()) #data type is float64     print can01.shape     ds2 = gdal.Open("cases2002.tiff")     can02 = np.array(ds2.GetRasterBand(1).ReadAsArray()) #data type is float 64     print can02.shape     #Load the RASTER CLIPPED REPROJECTED IMAGE AM I DOING SOME MISTAKE HERE?     ds3 = gdal.Open(r"H:\...\Annual 53-54\AOD2001.tiff")     AOD01 = np.array(ds3.GetRasterBand(1).ReadAsArray())     AOD01 =AOD01.astype('float64')  #Convertinf to float data type     print AOD01.shape     print AOD01      **Why am I not getting the required values? And instead ending up with the "no data values"?** For the files `cases2001` and `cases2002`, which are the original files according to which I re-projected my other files are having the perfect, intended values. But for the third file here i.e. `AOD2001` I'm getting an array of zeroes (I didn't change the default no data value). I've tried for both bands. But in `QGIS` "Identify Features" I'm getting the correct values for the 1st band. But not while reading ion Python. I've already tried masking and a few other suggestions for similar problems (such as`GetStatistics` not refreshing etc) but I haven't found exact problem anywhere yet. On a side note I need to plot these 2-D arrays against each other and perform regression. Which I'm obviously getting incorrect same values for each row as the `AOD2001` file contains only zeroes. Another question, **Can we compare (In my case perform regression on) two tiff files with same dimension, same spatial resolution and same end coordinates but different origins?** It's the first time I'm asking a question here, so I may have made mistakes in formatting etc, please point it out, I'd edit and take care in future. Also, please tell me if I should provide some other information or data. Thank you.


title: Zonal Statistics using Gdal and Python

text: I recently came across a loop_zonal_stats function by @ustroetz found here: Issue Trying to create Zonal Statistics using Gdal and Python.. I am using a (single) polygon shapefile and a raster file to compute the zonal statistics. The polygon shapefile (http://www.epa.gov/wed/pages/ecoregions/na_eco.htm#Downloads) is a list of North American ecoregions. I first opened the file in QGIS and saved it in the EPSG:4326 - WGS 84 coordinate reference system. The raster file (http://worldclim.org/current) contains precipitation data and is an ESRI grid file with 30 arc-second resolution. One issue I noticed with the code is that NoData values are included in the calculations (a "bug" that I noticed in QGIS as well). In order to fix this, I added the following lines of code to the zonal_stats function:               nodata = banddataraster.GetNoDataValue()     if nodata is not None:         zoneraster = np.ma.masked_equal(zoneraster, nodata)      After making this change, I am successfully able to exclude the NoData values from the calculation. However, I noticed something strange that I don't understand. Many of the polygons in my shapefile have the same attribute (and none of them overlap). I used QGIS to convert the shapefile of polygons to a shapefile of multipolygons (QGIS SinglePart to MultiPart) and then applied the loop_zonal_stats function to compute the sum and count for this multipolygon shapefile. However, the results from the original polygon shapefile are inconsistent with those from the multipolygon shapefile. For example, for a particular region with "NA_L3CODE" = 10.1.2, I get a count of 352 + 6873 + 317400 = 324625 for the original polygon shapefile and a count of 369288 for the multipolygon shapefile. Might anyone know what is going on?


title: Where to start rtk and Rinex?

text: I'm a total newbie at this filed and I want to use Rinex data to add more precision to my GPS receiver. I don't know where to start and even I spend a lot of time googling, I just got more confused.I would appreciate if anyone can help me with a starting point. **Edit 1:** I want to determine the exact 3D position of a flying quadrotor while flying (So i need the position realtime). The main project is multi-agent UAV flying, so their position has to be accurate enough to prevent crashing during flight. Also for accurate landing in a desired position. In the first phase I need to know position relative to a ground station in a not so far distance (I mean about 5-7 meter from robot). And this station (the positioning reference) doesn't move.


title: OSM digitizing road map by shapefile

text: I have shapefile with roads. I want to digitize some of the objects in shapefile to OSM, because I allways use OSM shapefiles and I have something to give to. What software should I use? I have JOSM installed. Is it enough?


title: How to use replaceDataSource to replace data source for multiple layers?

text: I appear to be having difficulty with using 'replaceDataSource'. ESRI has atrocious documentation on updating and fixing data sources, and almost everything I have read on the subject is almost what I am looking for, but not quite. What I am looking for is a way to repair multiple map documents without having to specify the name of each layer that has a broken data source. The code I have written so far **_should_** loop through the map documents in a directory, find which layers have a broken data source and determine if it is a shapefile, personal geodatabase feature class, file geodatabase feature class or a raster. The code determines the name of the feature class or shapefile, and uses the 'replaceDataSource' to fix the broken data source. When the script is run in the Python window in ArcMap, it points to their being an error on the lines that contain 'lyr.replaceDataSource' bit of code along with **Info: Layer: Unexpected error** When the script is run in PyScripter, the first datasource is a shapefile and it gives an error message but keeps running: **Info: Layer: Unexpected error** and nothing else. However, as soon as it runs into a feature class in a personal geodatabase, python crashes and gives this message: **exceptions.EOFError: [Errno 10054] An existing connection was forcibly closed by the remote host** I suspect the primary reason the script is not working may be the manner is which the broken data sources are determined, but really have no idea. Currently using ArcGIS 10.2, python version 2.7.3 and PyScripter. Here is an edited, pared down version of the code, with the full code following:               import arcpy, os, sys, traceback          workspace = arcpy.env.workspace = r'c:\users\athom\Documents\GIS\GIS_computer_06_20_2013\GIS_LOCAL_MASTER\MAPS\Map_Projects\TEST'     mxdList = arcpy.ListFiles("*.mxd")     for Map in mxdList:         filePath = os.path.join(workspace, Map)         mxd = arcpy.mapping.MapDocument(filePath)         brkLayerList = arcpy.mapping.ListBrokenDataSources(mxd)         for lyr in brkLayerList:             if lyr.supports("DATASOURCE") and lyr.supports("WORKSPACEPATH"):                 brokenDataSource = lyr.dataSource                 if brokenDataSource.startswith('n'):                     newDataSource = brokenDataSource.replace(r'n:', r'c:\users\athom\Documents\GIS\GIS_computer_06_20_2013\GIS_LOCAL_MASTER')                     layerName = lyr.datasetName                     ChangedDataSet = newDataSource.replace(layerName, '')                     try:                         if "shp" in newDataSource:                             ChangedDataSet =ChangedDataSet.replace(' ','')[:-4]                             layerName = layerName+'.shp'                             lyr.replaceDataSource(ChangedDataSet, "SHAPEFILE_WORKSPACE", layerName, True)                             lyr.save()                         elif "mdb" in newDataSource:                             lyr.replaceDataSource(ChangedDataSet, "ACCESS_WORKSPACE", layerName, True)                             lyr.save()                         elif "gdb" in newDataSource:                             lyr.replaceDataSource(ChangedDataSet, "FILEGDB_WORKSPACE", layerName, True)                             lyr.save()                         else:                             lyr.replaceDataSource(ChangedDataSet, "RASTER_WORKSPACE", layerName, True)                             lyr.save()                     except Exception as e:                         print arcpy.AddError(e)                         print "N FAILED"                 else:                     print 'Data Source Lost'              mxd.save()         del mxd      Full code with error checking.               import arcpy, os, sys, traceback          workspace = arcpy.env.workspace = r'c:\users\athom\Documents\GIS\GIS_computer_06_20_2013\GIS_LOCAL_MASTER\MAPS\Map_Projects\TEST'     mxdList = arcpy.ListFiles("*.mxd")     for Map in mxdList:         filePath = os.path.join(workspace, Map)         mxd = arcpy.mapping.MapDocument(filePath)         brkLayerList = arcpy.mapping.ListBrokenDataSources(mxd)         for lyr in brkLayerList:             if lyr.supports("DATASOURCE") and lyr.supports("WORKSPACEPATH"):                 brokenDataSource = lyr.dataSource                 print brokenDataSource                 if brokenDataSource.startswith('N'):                     newDataSource = brokenDataSource.replace(r'N:', r'c:\users\athom\Documents\GIS\GIS_computer_06_20_2013\GIS_LOCAL_MASTER')                     layerName = lyr.datasetName                     ChangedDataSet = newDataSource.replace(layerName, '')                     print newDataSource                     print ChangedDataSet                     print layerName                     try:                         if "shp" in newDataSource:                             try:                                 ChangedDataSet =ChangedDataSet.replace(' ','')[:-4]                                 print ChangedDataSet                                 lyr.replaceDataSource(ChangedDataSet, "SHAPEFILE_WORKSPACE", layerName, True)                                 lyr.save()                             except:                                 tb = sys.exc_info()[2]                                 tbinfo = traceback.format_tb(tb)[0]                                 pymsg = "PYTHON ERRORS:\nTraceback info:\n" + tbinfo + "\nError Info:\n" + str(sys.exc_info()[1])                                 msgs = "ArcPy ERRORS:\n" + arcpy.GetMessages(2) + "\n"                                 arcpy.AddError(pymsg)                                 arcpy.AddError(msgs)                                 print pymsg + "\n"                                 print msgs                              elif "mdb" in newDataSource:                             try:                                 lyr.replaceDataSource(ChangedDataSet, "ACCESS_WORKSPACE", layerName, True)                             except:                                 tb = sys.exc_info()[2]                                 tbinfo = traceback.format_tb(tb)[0]                                 pymsg = "PYTHON ERRORS:\nTraceback info:\n" + tbinfo + "\nError Info:\n" + str(sys.exc_info()[1])                                 msgs = "ArcPy ERRORS:\n" + arcpy.GetMessages(2) + "\n"                                 arcpy.AddError(pymsg)                                 arcpy.AddError(msgs)                                 print pymsg + "\n"                                 print msgs                              elif "gdb" in newDataSource:                             try:                                 lyr.replaceDataSource(ChangedDataSet, "FILEGDB_WORKSPACE", layerName, True)                             except:                                 tb = sys.exc_info()[2]                                 tbinfo = traceback.format_tb(tb)[0]                                 pymsg = "PYTHON ERRORS:\nTraceback info:\n" + tbinfo + "\nError Info:\n" + str(sys.exc_info()[1])                                 msgs = "ArcPy ERRORS:\n" + arcpy.GetMessages(2) + "\n"                                 arcpy.AddError(pymsg)                                 arcpy.AddError(msgs)                                 print pymsg + "\n"                                 print msgs                         else:                             try:                                 lyr.replaceDataSource(ChangedDataSet, "RASTER_WORKSPACE", layerName, True)                             except:                                 tb = sys.exc_info()[2]                                 tbinfo = traceback.format_tb(tb)[0]                                 pymsg = "PYTHON ERRORS:\nTraceback info:\n" + tbinfo + "\nError Info:\n" + str(sys.exc_info()[1])                                 msgs = "ArcPy ERRORS:\n" + arcpy.GetMessages(2) + "\n"                                 arcpy.AddError(pymsg)                                 arcpy.AddError(msgs)                                 print pymsg + "\n"                                 print msgs                     except Exception as e:                         print arcpy.AddError(e)                         print "N FAILED"                      elif brokenDataSource.startswith ('n'):                     newDataSource = brokenDataSource.replace(r'n:', r'c:\users\athom\Documents\GIS\GIS_computer_06_20_2013\GIS_LOCAL_MASTER')                     layerName = lyr.datasetName                     ChangedDataSet = newDataSource.replace(layerName, '')                     print newDataSource                     print ChangedDataSet                     print layerName                     try:                         if "shp" in newDataSource:                             try:                                 ChangedDataSet =ChangedDataSet.replace(' ','')[:-4]                                 print ChangedDataSet                                 lyr.replaceDataSource(ChangedDataSet, "SHAPEFILE_WORKSPACE", layerName, True)                             except:                                 tb = sys.exc_info()[2]                                 tbinfo = traceback.format_tb(tb)[0]                                 pymsg = "PYTHON ERRORS:\nTraceback info:\n" + tbinfo + "\nError Info:\n" + str(sys.exc_info()[1])                                 msgs = "ArcPy ERRORS:\n" + arcpy.GetMessages(2) + "\n"                                 arcpy.AddError(pymsg)                                 arcpy.AddError(msgs)                                 print pymsg + "\n"                                 print msgs                         elif "mdb" in newDataSource:                             try:                                 lyr.replaceDataSource(ChangedDataSet, "ACCESS_WORKSPACE", layerName, True)                             except:                                 tb = sys.exc_info()[2]                                 tbinfo = traceback.format_tb(tb)[0]                                 pymsg = "PYTHON ERRORS:\nTraceback info:\n" + tbinfo + "\nError Info:\n" + str(sys.exc_info()[1])                                 msgs = "ArcPy ERRORS:\n" + arcpy.GetMessages(2) + "\n"                                 arcpy.AddError(pymsg)                                 arcpy.AddError(msgs)                                 print pymsg + "\n"                                 print msgs                              elif "gdb" in newDataSource:                             try:                                 lyr.replaceDataSource(ChangedDataSet, "FILEGDB_WORKSPACE", layerName, True)                             except:                                 tb = sys.exc_info()[2]                                 tbinfo = traceback.format_tb(tb)[0]                                 pymsg = "PYTHON ERRORS:\nTraceback info:\n" + tbinfo + "\nError Info:\n" + str(sys.exc_info()[1])                                 msgs = "ArcPy ERRORS:\n" + arcpy.GetMessages(2) + "\n"                                 arcpy.AddError(pymsg)                                 arcpy.AddError(msgs)                                 print pymsg + "\n"                                 print msgs                         else:                             try:                                 lyr.replaceDataSource(ChangedDataSet, "RASTER_WORKSPACE", layerName, True)                             except:                                 tb = sys.exc_info()[2]                                 tbinfo = traceback.format_tb(tb)[0]                                 pymsg = "PYTHON ERRORS:\nTraceback info:\n" + tbinfo + "\nError Info:\n" + str(sys.exc_info()[1])                                 msgs = "ArcPy ERRORS:\n" + arcpy.GetMessages(2) + "\n"                                 arcpy.AddError(pymsg)                                 arcpy.AddError(msgs)                                 print pymsg + "\n"                                 print msgs                     except arcpy.ExecuteError:                         arcpy.AddError(GetMessages(2))                         print "n FAILED"                      else:                     print 'Data Source Lost'              mxd.save()         del mxd


title: How to resample polygon vertices so that they are evenly spaced in Arc 10

text: I have digitized the outline of many river channels (by both manually picking points and using streaming mode), but now I need the vertices to be evenly spaced (for input to another program). I tried to use the "densify" tool which was partially successful, but this adds additional vertices at the specified distance whereas I wanted to also simplify the sections that were streamed and have too many vertices. Has anyone any suggestions for a tool I might try? I have no experience in writing scripts, but may be able to modify another. Many thanks in advance.


title: Is it possible to generate a heatmap on the server-side using ArcGIS Server?

text: Using ArcGIS Server's JS API it is possible to generate a heatmap at runtime, as in the example at http://tmappsevents.esri.com/website/heatmap/index.html This requires having client-side access to the individual points in the heatmap, for example by downloading a feature layer to the browser. The problem is that the browser becomes unusably slow when working with more than a few thousand points. Has anyone been able to generate the heatmap on the server-side, and send it down as a cached or dynamic layer? I need the ability to generate the heatmap at runtime, so that the points could be filtered with the change immediately visible in the heatmap.


title: How to randomly select from existing points within each polygon?

text: I have a polygon layer (where each feature has a class attribute), and a point layer, and there are at least 10 points within each feature of the polygon. Is there a way to randomly select only "k" points of these 10 points for each feature (class) in QGIS? ![select only two of ten points that are within each polygon](http://i.stack.imgur.com/kRWI9.jpg) Thanks!


title: export postgis database

text: I'm totally new to the awesome PostGIS world, so sorry for the maybe silly question. After a while I was able to install postgres and postgis on my computer (debian testing) and I was also able to create some databases. Through the DB Manager of QGIS I have imported some shapefile into the database and that's ok. But now, if I want to export the entire postgis db as a single file (like SpatiaLite does) how can I do that? I mean, Spatialite creates a _real_ file into a directory while postgis no (or at least I haven't found it). Thanks guys and sorry for the simple question


title: qgis's georeferencer renders my white map with a translucent blue tint, can I prevent this?

text: I'm running QGIS 2.0.1-Dufour on MacOSX and when loading a PNG or TIFF into qgis as a raster its always getting rendered with a blue tint. Here's what it looks like before being loaded into QGis's Geoferencer. ![](http://i.stack.imgur.com/aTWnC.png) and after: ![enter image description here](http://i.stack.imgur.com/wNiTv.png) Is there any way I could save this file in photoshop or gimp that would allow it to maintain its original colorization? GDALINFO (png)                   davodavo$ gdalinfo color\ key.png         Warning 1: libpng: iCCP: Not recognizing known sRGB profile that has been edited         Driver: PNG/Portable Network Graphics         Files: color key.png         Size is 308, 588         Coordinate System is `'         Metadata:           XML_com.adobe.xmp=<x:xmpmeta xmlns:x="adobe:ns:meta/" x:xmptk="XMP Core 5.4.0">            <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">               <rdf:Description rdf:about=""                     xmlns:exif="http://ns.adobe.com/exif/1.0/">                  <exif:PixelXDimension>308</exif:PixelXDimension>                  <exif:PixelYDimension>588</exif:PixelYDimension>               </rdf:Description>            </rdf:RDF>         </x:xmpmeta>              Image Structure Metadata:           INTERLEAVE=PIXEL         Corner Coordinates:         Upper Left  (    0.0,    0.0)         Lower Left  (    0.0,  588.0)         Upper Right (  308.0,    0.0)         Lower Right (  308.0,  588.0)         Center      (  154.0,  294.0)         Band 1 Block=308x1 Type=Byte, ColorInterp=Red           Mask Flags: PER_DATASET ALPHA          Band 2 Block=308x1 Type=Byte, ColorInterp=Green           Mask Flags: PER_DATASET ALPHA          Band 3 Block=308x1 Type=Byte, ColorInterp=Blue           Mask Flags: PER_DATASET ALPHA          Band 4 Block=308x1 Type=Byte, ColorInterp=Alpha      gdalinfo (tiff) -- same blue tint                   Driver: GTiff/GeoTIFF         Files: colorkey.tif         Size is 308, 588         Coordinate System is `'         Metadata:           TIFFTAG_DATETIME=2014:07:09 17:27:31           TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)           TIFFTAG_SOFTWARE=Adobe Photoshop CS6 (Macintosh)           TIFFTAG_XRESOLUTION=143.9926           TIFFTAG_YRESOLUTION=143.9926         Image Structure Metadata:           INTERLEAVE=PIXEL         Corner Coordinates:         Upper Left  (    0.0,    0.0)         Lower Left  (    0.0,  588.0)         Upper Right (  308.0,    0.0)         Lower Right (  308.0,  588.0)         Center      (  154.0,  294.0)         Band 1 Block=308x8 Type=Byte, ColorInterp=Red         Band 2 Block=308x8 Type=Byte, ColorInterp=Green         Band 3 Block=308x8 Type=Byte, ColorInterp=Blue


title: Markers not displayed in xhtml in OpenLayers

text: I just started working with OpenLayers because I want to include a map on my JSF page (.xhtml files). I want to use the same feature that is demonstrated in this example: http://openlayers.org/dev/examples/styles-context.html I played with the code a little bit in the html example files that come with the .zip download. When I copied the code into my JSF .xhtml page, the map appears but the markers are not displayed. Here is the example code:               <script type="text/javascript">     //<![CDATA[     var map;     function init(){      map = new OpenLayers.Map('map', {maxResolution:'auto'});     var wms = new OpenLayers.Layer.WMS.InlineXhtml( "OpenLayers WMS",      "http://vmap0.tiles.osgeo.org/wms/vmap0", {layers: 'basic', format: 'image/svg+xml'} );      map.addLayer(wms);          // create 50 random features in the northern hemisphere      // give them a "type" attribute that will be used to style     // them by size     var features = new Array(50);                 for (var i=0; i<features.length; i++) {      features[i] = new OpenLayers.Feature.Vector(     new OpenLayers.Geometry.Point(                             (360 * Math.random()) - 180, 90 * Math.random()     ), {      type: 5 + parseInt(5 * Math.random())     }     );      }           // allow testing of specific renderers via "?renderer=Canvas", etc      var renderer = OpenLayers.Util.getParameters(window.location.href).renderer;      renderer = (renderer) ? [renderer] : OpenLayers.Layer.Vector.prototype.renderers;                      // create the layer styleMap with a simple symbolizer template     yer1 = new OpenLayers.Layer.Vector('Points', {     styleMap: new OpenLayers.StyleMap({     pointRadius: "${type}", // based on feature.attributes.type      fillColor: "#000000"      }),      renderers: renderer     });     layer1.addFeatures(features);           map.addLayers([layer1]);     map.zoomToExtent(layer1.getDataExtent());      }     //]]>      </script>      I added:               map.zoomToExtent(layer1.getDataExtent());      and changed the markers from 50 random to 2 markers of the coordinates of London and Berlin. And the map was centered exactly to show them, but the markers are still not displayed. But that proves that the markers are read and loaded to the map correctly, just not displayed. If I do the same back in the example .html page the map is again centered properly and the markers are displayed. I have uploaded the example .html file to my webapp and it works perfectly. So this seems to be the .xhtml issue... Is it the problem of xhtml? How to fix this? EDIT: Using the google console tool I analysed the code that was generated and it turnes out that in the .xhtml page the radious is zero:               <circle id="OpenLayers.Geometry.Point_41" cx="483.55555555555566" cy="150.7555555555556" r="" fill="#000000" fill-opacity="1" stroke="#000000" stroke-opacity="1" stroke-width="2" stroke-linecap="round" stroke-linejoin="round"></circle>      Please notice that it's always r="" but why is that? The script code is IDENTICAL in both cases - one (.html) produces the right radios, the other(.xhtml) one doesn't...


title: Minimum system requirements for ArcGIS on Mac

text: Just looking to buy a new-ish computer for study for my PhD. I will be using ArcGIS. I know the ideally it would be best to buy all the memory and ram I can afford but I am self-funding on a low budget. I am wondering what's the minimum spec Apple Mac I could use to run this software practically?


title: Qgis OpenStreetMap plugin misbehaviour

text: I used the plugin successfully. But yesterday, after changing several buildings in JOSM. Import under qgis with the plugin does not work anymore. The buildings are in the layer point! Do you have any idea? thank you! I can not post image because it requires 10 reputation.


title: Which remote sensing satellite is best used for vegetation extent mapping

text: I'm currently writing an assignment for my remote sensing class and while I think I have correct answer was wondering what people who do this stuff for a living think. The question is: which satilite out the **Landsat TM (Thematic Mapper)** and **SPOT 5** would you use to map general extent of a 300km x 300km study area. My answer was that you would use the **SPOT 5** as it is higher resolution and would let you get a finer extent vs the 30m resolution of the Landsat TM. However the SPOT 5 has a small swath area so you have to use more images. I also thought about the 2200km VMI swath on the SPOT 5 but the resolution is 1km. Thoughts?


title: Error in Select Layer by Location python script?

text: I keep on getting an error on Select by location that the file cannot be found and the inputs are invalid, but I can't figure out what I'm doing incorrectly. I realize some of my variables unnecessarilly repeat themselves, which I plan to clean up, but it still doesn't explain why SelectByLocation isn't working.               # Import arcpy module     import arcpy          #Set Geoprocessing Environments     arcpy.env.overwriteOutput = True     arcpy.env.scratchWorkspace = "C:/PATH HERE/Testing"     arcpy.env.workspace = "C:/PATH HERE/ArcGIS/gdb.mdb"               # Parameters.     PRISM = arcpy.GetParameterAsText(0)     PRISM = "C:/LOCATION HERE"          #Unchanging Variables and Workspaces     Map = "C:/LOCATION HERE"     workspace = "C:/PATH HERE/gdb.mdb"     scratchWorkspace = "C:/PATH HERE/Testing"     PRISM_pts = scratchWorkspace + "/PRISM_pts.shp"     Pts_Lyr = scratchWorkspace + "/Pts_Lyr.lyr"     Map_Layer = scratchWorkspace + "/Map_Layer.lyr"     Prism_Lyr = scratchWorkspace + "/Prism_Lyr.lyr"          # Process: Raster to Point     arcpy.RasterToPoint_conversion(PRISM, PRISM_pts, "Value")          #Process: Make a Feature Layer for the Rasters     arcpy.MakeFeatureLayer_management(PRISM_pts, Pts_Lyr)          #Process: Make a Feature Layer for the Map     arcpy.MakeFeatureLayer_management(Map, Map_Layer)               # Process: Select Layer By Location     arcpy.SelectLayerByLocation_management(Map_Layer, "INTERSECT", Pts_Lyr, "", "")          #====================================     # Runtime error  Traceback (most recent call last):     # File "<string>", line 1, in <module>     # File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\management.py",     # line 6618, in SelectLayerByLocation     # raise e ExecuteError: Failed to execute. Parameters are not valid.     # ERROR 000368: Invalid input data.     # ERROR 000732: Selecting Features:     # Dataset C:/PATH HERE/Testing/Pts_Lyr.lyr     # does not exist or is not supported     # Failed to execute (SelectLayerByLocation).     #     # =====================================               # Process: Spatial Join     arcpy.SpatialJoin_analysis(Map, Prism_Lyr, Spatial_Join_shp, "JOIN_ONE_TO_ONE", "KEEP_ALL", "", "INTERSECT", "", "")      I don't quite understand why ArcGIS can't find the layer by location file, since it just created it. Am I missing a step?


title: Geoservers communicating with other Geoservers

text: I am able to connect to multiple Geoservers using the Geonetwork but is it possible for two Geoservers to communicate with each other without the Geonetwork?


title: How to shift the Map a little bit?

text: I'm having trouble with two maps, both UTM33N/WGS84. They should be exactly on top of each other but there is a slight shift between them of probably a few hundred meters or so. ![enter image description here](http://i.stack.imgur.com/fMx1T.png) I Tried to adjust this with the Vector Affine Transformation Plugin but didn't succeed. I don't know how to use it.


title: WMTS GetFeatureInfo in android

text: I want to create something in android like this: http://openlayers.org/dev/examples/wmts-getfeatureinfo.html I'm able to load the costum tiles with WMTS but now i need to get data to show in some kind of a popup window when clicked(touched). I found that you can get extra info with the service getFeatureInfo, for this request i need some extra parameters like "I" and "J". Whitch are the number of pixels from the bottom-left point of the tile to the clicked position on the tile? But I have no idea how I can get these parameters in android. Is this the right method to do this or is there a better solution? Any advice would be appreciated.


title: mulitple line labeling arcgis

text: I try to label a field according to his length with VBScripting in Argis 10.0 The idea is : If myfield is more than a certain length, then please add a new line and split the field correctly. My field can be Null. Cant get the code below work (Verify retuns nothing) : `Function FindLabel ( [JOIN_NUMBER] ) If IsNull([JOIN_NUMBER] ) Then FindLabel = [JOIN_NUMBER] Else FindLabel = Left([JOIN_NUMBER] ,8) & vbnewline & Right([JOIN_NUMBER] ,Len([JOIN_NUMBER] ,-8)) End If End Function`


title: How can I Destripe Landsat Images?

text: I am working with Landsat 7 ETM+ SRC-off imagery. I'm doing some land-cover classification for a place in Kenya. My problem is the no-data stripes in the image. Is there a way in which I can interpolate data for such 'empty' areas using the neighbouring pixels?


title: how to keep relationship between feature classes after intersection one of feature classes?

text: I have geodatabase with 3 feature classes. one of them is circuit polyline feature class and two polygon feature class. substation and division in geodatabase there is relationship between circuit and substation. I am going to intersect between division and substation. the result of relation will be feature class with all substation and in attribute table all data about substations and division are available. but my problem is that there isnt the relation between circuit and substation in result feature class from intersect? how can I make intersect feature class between sub and div that relation between substation and other classes maintain?


title: ArcPy nested loop problem

text: I have a question connected with ArcPy. I was trying to make a nested loop with 2 cursors - 1st with 2 elements and 2nd with 5 elements. The statement should return 10 connections, e.g. 1-1, 1-2, ... 2-5. Here is part of my code:               import arcpy     import os.path          folder = "Z:\WOLGA"         na="na"         baza = "testy.mdb"         dataset = "network_analyst"         arcpy.env.workspace = os.path.join(folder, na, baza, dataset)         shapePocz = os.path.join(folder, na, baza, dataset, "dom")          listaPktPocz = arcpy.SearchCursor(shapePocz, "", "", "", "")         shapeKonc = os.path.join(folder, na, baza, dataset, "punkty_docelowe")         listaPktKonc = arcpy.SearchCursor(shapeKonc, "", "", "", "")         shapeKoncName = arcpy.Describe(shapeKonc).shapeFieldName         shapePoczName = arcpy.Describe(shapePocz).shapeFieldName          for pktPocz in listaPktPocz:         for pktKonc in listaPktKonc:             print str(pktPocz.numer) + " - " + str(pktKonc.numer)      For now I got such connections as: 1-1, 1-2, 1-3, 1-4, 1-5 only. The outer loop is not being executed properly. Pls help!


title: How do I specify the position, style and units of the ScaleBar in pyqgis?

text: I would like to add a scale bar to my composition using the python console and then export it to a raster image.   I am following the MapRendering tutorial with the "output using map composer " approach however,the Scale Bar item example do not have many features. I would like   * a scale bar with 2 ticks, one at the beginning, one at the end    * units in km   * eventually round it to the nearest km   * place it on the upper right of my figure.    * set a transparent background  Just like in the desktop QGIS version. What are the options ? So far I have found               item = QgsComposerScaleBar(c)     item.setStyle("Line Ticks Up") # optionally modify the style     item.setComposerMap(composerMap)     c.addItem(item)      Will the scale bar be "at scale" ? Thank you for your help I cannnot find anything from that page ScalBar class


title: QGIS crashes when trying to open ecw files

text: I added Ubuntugis unstable repo and installed Qgis, libgdal-ecw-src, libecwj2-3.3 but still does not work. When I attempt to open an ecw file with Qgis it crashes and tells me: 'std::length_error': basic_string::_S_create I am using ubuntu 13.04 and when I do: gdalinfo --formats|grep -i ecw it seems to work. ECW (rw): ERDAS Compressed Wavelets (SDK 3.x) JP2ECW (rw+v): ERDAS JPEG2000 (SDK 3.x) What is the problem?


title: Scripting (with) WPS?

text: As far as I understand WPS (web processing service) can be used to execute geospatial operations on the server side. Is it possible to create scripts with WPS? For example something like: 'check in a directory if there is a new file' - if so , 'process this file', 'store the result in the database' I assume that I'm expecting too much, but what is the best solution for an script when using qgis server or geoserver? D3


title: Any example of creating a time-aware layer (acrgis js api + ruby rails) using data from the server

text: I have a data set of points. Each point has latitude, longitude, date and description sent from the server I already have a graphic layer which is used for represening the points in a static way with some popups like this                         ...               var graphic = new esri.Graphic(point, symbol);                    graphic.setInfoTemplate(new esri.InfoTemplate("", "Number: " + data.number));                    gLayer.add(graphic);               ...      But Now I need to add a time-aware layer (a feature layer I believe) that allows me to show the data according to datetime information using a time slider, like here oil Filed. The only difference, I need to use the data sent from the server as json and I need to add popups to each point I believe I need to add a featureLayer (MODE_SNAPSHOT) with a featureCollection as params that contains the data sent from server but I coulnd't figgure out how to do this (I'm a newbie in ArcGis, several hours only). To be exact I'm not sure how to create that feature collection with latitude/longitude for each point and the most important thing: how to connect it with corresponding time information I.e. I have sets of {lat: xxx, long: yyy, datetime: yyyy-mm-dd, description: sometext } sent from the server Now I need to figure out how to pass the data to a featureLayer so it can be shown on a map with a time slider Any example/help will be appreciated thanks a lot


title: What can go wrong in UTM to WGS1984 conversion (shifted data)

text: I try to convert .shp files to OSM data with shp2osm and ogr2osm but the output data is shifted. _Semi-Cross-Posting Warning: I have already asked a questions on theOSM help site, but this questions is more theoretical._ The .prj file I have contains the following (line-breaks added for readability) and the source data was saved using ArcGIS:               PROJCS["Adindan_UTM_Zone_37N",     GEOGCS["GCS_Adindan",     DATUM["D_Adindan",     SPHEROID["Clarke_1880_RGS",6378249.145,293.465]],     PRIMEM["Greenwich",0.0],     UNIT["Degree",0.0174532925199433]],     PROJECTION["Transverse_Mercator"],     PARAMETER["False_Easting",500000.0],     PARAMETER["False_Northing",0.0],     PARAMETER["Central_Meridian",39.0],     PARAMETER["Scale_Factor",0.9996],     PARAMETER["Latitude_Of_Origin",0.0],     UNIT["Meter",1.0]]      I am using ogr2ogr for the coordinate transformation:               ogr2ogr -t_srs WGS84 -s_srs "ESRI::subcities.prj" "subcities_converted.shp" "subcities.shp"      When I finally import the data in the OpenStreetMap editor, the data is shifted by around 176m, 222°. Does this mean the source data and the .prj file is plain wrong? Or I am doing something wrong on the way?


title: ArcGIS Server geoprocessing service examples for shapefile data export

text: I have a web application client built using ArcGIS Server Manager which is currently showing a single map composed by three data layers. The map data are provided as a map service by my ArcGIS Server 10.1 instance. My task is to provide the map's users a way to download the full map locally as a shapefile (e.g: zipped) via her/his web browser. I was hinted to read here and to build a geoprocessing service to accomplish the task. I decided to write a specific toolbox myself using Python, then publish it as a geoprocessing service on my ArcGIS Server instance and finally include it into my web application. Can you point me towards any third-party (meaning: non-ESRI) code sample to start me up? I can find no examples other than the ones in ESRI's support website! Thanks in advance


title: Cannot add Plugins from Python: "Plugin is Broken"

text: I'm using the latest version of QGIS on a Windows 7 (32 bits), and when I try to add a plugin from the option "Fetch Python Plugins", it simply cancels the installation and I can see the following error: `Couldn't load plugin 'xytools' from ['C:/PROGRA~1/Quantum GIS Lisboa/apps/qgis/./python', 'C:/Users/ProprietÃ¡rio/.qgis//python', 'C:/Users/ProprietÃ¡rio/.qgis//python/plugins', 'C:/PROGRA~1/Quantum GIS Lisboa/apps/qgis/./python/plugins', 'C:\Program Files\ArcGIS\bin', 'C:\PROGRA~1\Quantum GIS Lisboa\bin\python27.zip', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\DLLs', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\plat-win', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\lib-tk', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\qgis\bin', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\PIL', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32\lib', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\Pythonwin', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\wx-2.8-msw- unicode', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\qgis\python\plugins\fTools\tools']` The plugin that was not successfully installed appears in red colour saying "Invalid" - and in the description field: The plugin is Broken. It happens with all plugins from the QGIS Plugin Installer. Does anyone knows how can I fix this problem? Thank you in advance for any help you can provide.


title: Sending and receiving gis data to map

text: I have a spatial database in Sql server and I am working in asp.net mvc. My project must show stored shape in web page that uses openlayers js and user must see and edit shapes and send the result to server. I want to know how can I send data to my map and receive new data from it and save data to database. Is there any library that work for me? Is is not important that it must be open source thanks a lot


title: Create DEM from stereo pairs?

text: Does ArcGIS Desktop or any of its native extensions allow me to create DEMs from stereo pairs (acquired by satellite, GeoEye 1)? If the answer is a clear no, brownie points will go to an answer that proposes a viable alternative using OpenSource.


title: Bug in gdalwarp's simple nearest neighbor? Or more likely I need some help

text: I'm trying to understand how to use the GDAL utilities a little better. I think I've defined something incorrectly or I've found a bug in gdalwarp. My bet is that I've defined something incorrectly, I would love some help in understanding this better. I'm getting answers for a simple nearest neighbor reprojection that I don't understand. I apologize in advance for the length of this post, but it's hard to describe in any fewer words. If someone who's a gdal master could follow it along and show me where I fell of the path, I'd really appreciate it. I start with a 5x5 raw 2-byte image, with values 1 through 25 ( [[1, 2, 3, 4, 5][6,...,10][...]]) file:ease_subset_324_376_5x5.bin I'm letting this represent specific grid cells in a 25km EASE2 projection.   _ISPRS Int. J. Geo-Inf. 2012, 1, 32-45; doi:10.3390/ijgi1010032_   This projection has a proj4 string: `'+proj=laea +lat_0=90. +lon_0=0. +x_0=0. +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m'` In my case the upper left corner of the upper left grid cell is located at:   `-900000.00, -400000.00` Which means the lower right corner of the lower right grid cell is located at:   `-900000.00 + 5 * 25000., -400000.00 - (5 * 25000.)` or `-775000., -525000.` I want to convert this to a format that contains the projection and location information so I create a VRT file to desibe the raw data:               <VRTDataset rasterXSize="5" rasterYSize="5">       <VRTRasterBand dataType="Int16" band="1" subClass="VRTRawRasterBand">         <SourceFilename relativetoVRT="1">ease_subset_324_376_5x5.bin</SourceFilename>         <ImageOffset>0</ImageOffset>         <PixelOffset>2</PixelOffset>         <LineOffset>10</LineOffset>         <ByteOrder>LSB</ByteOrder>       </VRTRasterBand>     </VRTDataset>      And use gdal_translate to create a GeoTiff that describes my data:               gdal_translate subset.vrt subset_EASE2_N25km_gdal.tif \         -a_srs '+proj=laea +lat_0=90. +lon_0=0. +x_0=0. +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m' \         -a_ullr -900000.00 -400000.00 -775000.00 -525000.00      Next I want to reproject this image to a 114x105 2km polar stereographic grid. The upper left corner of the upper left gridcell is located: `-366000.00 -782000.00` and the lower right corner of the lower right gridcell is: `-366000 + (114 * 2000.), -782000.00 - (105 * 2000.)` or: `-138000. -992000.` I use gdal warp to project my image into my new stereogrphic projection.               gdalwarp -t_srs "+proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs" \          -te -366000.00 -992000.00  -138000.00 -782000.00 \          -tr 2000 2000 subset_EASE2_N25km_gdal.tif subset_Moge2km_gdal.tif      I now have an image 114x105 with values 0-25 like I expect, but there are a few places, where I don't understand the values in the reprojected data. **For example: The value in gridcell[50, 26] => 5.** I would presume nearest neighbor sampling, would take the center of that gridcell, project it into the original image and select the closest value from there, but that's not what I see. the location of the center of polarstereo gridcell[50, 26] is (2km cells)   `-366000.00 + 2000./2 + 50 * 2000., -782000.00 - 2000./2 - 26 * 2000.`   or: `-265000.00, -835000.00` Using gdaltransform I can find the location in my original image:                gdaltransform           -s_srs "+proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"           -t_srs "+proj=laea +lat_0=90. +lon_0=0. +x_0=0. +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m"          -265000.00, -835000.00     -800001.11880625 -414546.034290511 0      Converting this x/y onto my original EASEgrid2 grid cell. Center of upper left gridcell = `-900000 + 25000./2, -400000 - 25000./2.` or: `-887500., -412500.`               xmeters = x - x0 = -800001.11880625 - (-887500.) = 87498.88119374996     ymeters = y0 - y = -412500 - (-414546.034290511) =  2046.034290510986      Convert to gridcells by dividing by resolution:               x_gridcell = 87498.88119374996 / 25000. =  3.4999552     y_gridcell = 2046.034290510986 / 25000. =  0.081841372      This to me would say that I would want the value at 3,0, but I was getting the value at 4,0. grid[3,0] = 4, grid[4,0] = 5 Also if you just look at the meters:               back_projected_point = -800001.11880625, -414546.034290511     grid at 3,0 = g3 => -812500.00      -412500.00     grid at 4,0 = g4 => -787500.00      -412500.00          distance(back_projected_point, g3) = 12665.239     distance(back_projected_point, g4) = 12667.447      Again, this to me looks like I should be getting the value at gridcell[3,0], but that's not what I see. Following this logic: All of these points in my polar stereo projected image are incorrect.               points =  [ [ 47, 88 ], $                 [ 59, 87 ], $                 [ 60, 88 ], $                 [ 61, 89 ], $                 [ 62, 90 ], $                 [ 50, 26 ], $                 [ 49, 25 ], $                 [ 48, 24 ], $                 [ 47, 23 ], $                 [ 46, 22 ], $                 [ 45, 21 ] ]      Again, all I can think of is that I'm not properly specifying the gdal parameters, or there's something wrong with gdal.


title: How to open a file geodatabase with tables and relationships in QGIS?

text: How can I open or connect to an ESRI file geodatabase in QGIS? I don't want to only see the features (shapes), I want to make sure that tables and relationships are opened correctly, too. Thanks in advance


title: How to select several disseminate features with ArcGIS JS API?

text: I have a cadastral map similar to this one from an ArcGIS Javascript API example: Cadastral map What I would like to do is selecting several parcels clicking on them and show the information related in a sidebar or floating panel. For example, I would like to select that big house on the right side and then few more spread on the map. Every time I click on one, I would like the application to show the information associated **AND** not losing that selection when clicking on a new one. In the end, I should have an array of selected features (or similar) that maybe I can re-use later. I have tried Feature Layer with Selection example, but it does not fully fit my requirements, because depending on the distance between features, I can't select two of them without selecting many more that I do not need. Is this possible? Can anyone explain the process? Thanks in advance!


title: ArcGIS arcpy reprojecting raster issue

text: I have a raster which I want to reproject from GCS_WGS_1984 to GCS_NorthAmerican_1983. I am using the following code:                   try:               spatial_ref = arcpy.CreateSpatialReference_management("GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]];-400 -400 1000000000;-100000 10000;-100000 10000;8.98315284119521E-09;0.001;0.001;IsHighPrecision", "", "", "", "", "", "0")         except:             print arcpy.GetMessages()              try:             arcpy.Project_management(o_dir+int_ras,o_dir+geo_int_ras,spatial_ref)         except:             print arcpy.GetMessages()      It gives the error: ERROR 000732: Input Dataset or Feature Class: Dataset does not exist or not supporte However, the input datasets most defn exist. Is there some issue with how I am reprojecting?


title: How to fix spatial reference error on CopyFeatures ArcGIS tool?

text: When I run the CopyFeatures tool to copy features from an ArcSDE geodatabase to another empty ArcSDE geodatabase (no features, no tables), I get the following result: ![enter image description here](http://i.stack.imgur.com/tLOyi.png) The error message is about the spatial reference of the output feature class. The output is created by the tool itself. Why wouldn't it set the correct spatial reference? Do you have any hint on how fixing this?


title: how to compute tile row from latitude in mercator projection

text: I work on google tile images. I want to compute a tile row from its latitude. How can I do that. I know that tile system cover all of the earth from -85 latitude to 85 latitude. 4 tiles exists at the first level. Google uses mercator projection.


title: What do I need to create 3D model of real terrain?

text: I'd like to create 3D model of some real terrain (to be exact it's about some real roads and only as much as 80m to 100m of roadside - 40m to 50m of each side). As I don't have access to highly detailed local maps, my only source of topographic information of given terrain will be Internet. Now my question - what would be the best approach to create such 3D model and make it as detailed as possible? I know I won't be able to get very detailed elevation data (AFAIK, DEM for Europe has only 90m resolution) but on the other hand it's still better than nothing. But also I don't know what would be the best source of detailed topographic data? Should I use Google Maps and work on its KML files or OpenStreetMap gives better results? Or maybe there are some other resources for getting the information I will need? BTW, would be some in-car GPS system of any help in this case? I'd like to add that I want to model only the surface - the vegetation, buildings and other objects are not needed. Also I won't need any geographic co-ordinates of terrain as I need to have only right measurements in meters. I'm totally new to this subject so any kind of help will be very much appreciated. Thank you in advance for reply. EDIT: Thank you everyone for your comments. I'll try to add more information on this project. Source of data - AFAIK, the most detailed DEMs for Europe (Western and South Europe and also North Africa are my areas of interest) can be obtained from ASTER and Landsat. Both have 30m resolution so it's as accurate as I need it to be. 3D model of terrain - I'd like to reconstruct some real roads that used to be parts of famous rallies or races in the 50's and 60's (Coupe des Alpes, Liege- Rome-Liege, Rallye del Sestriere or Tulpenrallye). Of course not whole roads as it would take ages for me, only some sections that can be about 7km long (that's game engine limitation). So I need a 3D model of surface of some specific road so I could edit it in 3DS Max 5 program. Or some other 3D modeling program (FreeCAD or Pov-Ray) as I need 3DS Max 5 only for texturing, adding some objects like buildings, trees, people, sponsor banners and finally exporting max file to bmf file (it's a format that Rally Trophy uses). And the game looks like in the link below: Rally Trophy - Game Intro Now, my biggest concern is the right method of generating such 3D model of terrain - what kind of programs do I need? And also how to mark a road on it? I'm sorry I cannot explain it any better. But please ask as many questions as you need to.


title: How to handle PostgreSql 1:m relationship data by using “Data-Driven Input Mask” QGIS plugin?

text: I have 1:m relationship tables in my PostgreSql. I like to show the main table attributes as well as the relational table attributes (like grid rows) in QGIS attribute-edit-form. What are the entries i have to add to dd_* tables to get this?


title: Finding overlapped areas and identifying their frequencies in a polygon layer

text: I have a fire patch polygon layer from 1998 to 2012 and I want to find the areas which experienced the fire several times! How I can find those overlapped areas and the times they experienced it (frequency)?


title: Convert pixel value histogram to table

text: I am working in QGIS 2.1 with the National Land Cover Database 2006 land cover layer. It is a raster layer with pixel values corresponding to different land cover types. I have clipped this raster to a buffer layer around my study site. Now I want to analyze the proportion of each land cover type within the site radius. I can view the frequency of each pixel value in the histogram tab of the layer properties window, but I would like to see the pixel frequency data in a table so that I can perform calculations on it. Is there a way to do this in QGIS?


title: why Semi-Automatic Classification plugin installation failled repeatedly?

text: I successfully downloaded Semi-Automatic Classification plugin installer OSGEO from http://download.osgeo.org/osgeo4w/osgeo4w-setup.exe. But it respond that it couldn't the installation. Again I download manually from http://plugins.qgis.org/plugins/plugins.xml and when I extract it in plugins folder still QGIS window respond that I commit error. how can I instal it?


title: How to integrate Google Drive with web mapping capabilities?

text: I use a shared Google Drive spreadsheet to coordinate a large GIS mapping project. The work involves processing thousands of DOQQ images. Essentially, work is divided into study areas and then further into DOQQ boundaries that make up the study areas. I am envisioning a dynamically updated web map that show which DOQQs have been processed within study areas. In other words, as I update the shared spreadsheet, the web map should also update. The attached image is a rough example of the type of data I wish to present. Is this feat possible with any web mapping services (e.g. CartoDB)? Which technologies should I pursue to accomplish this type of web mapping? ![enter image description here](http://i.stack.imgur.com/Fl7xi.jpg)


title: Am I correctly projecting a geometry using a WKID in ArcObjects?

text: I would like to take in any valid WKID for a spatial reference and then project a single geometry to it. The following code seems to be working, but (1.) can someone verify that this is the correct way to do this. (2) Can I just put the ProjectionHint parameter in Project5() to 0/false -- thanks               private IGeometry5 ProjectGeometry(IGeometry5 geom, int outSR)     {         IGeometryServer geomServer = new GeometryServerClass();         ISpatialReference spatialReference = geomServer.FindSRByWKID("EPSG", outSR, -1, true, true);         geom.Project5(spatialReference, 0);     }      http://help.arcgis.com/en/sdk/10.0/arcobjects_net/componenthelp/index.html#/Project5_Method/002m000001rz000000/ "The projectionHint parameter contains information about whether a geometry will cross the coordinate system horizons. If it doesn't, parts of the code that check this and intersect the geometry with the horizons can be omitted. This can dramatically improve performance. Use ISpatialReference3::ProjectionHint to initialize this parameter."               C#     public void Project5 (         ISpatialReference newSpatialReference,         int ProjectionHint     );


title: Where can I find European City boundary data?

text: Does anyone know any sources for free European city boundary data? Thanks!


title: OpenLayers - is it possible for a feature to have both a single click and double click event?

text: I'm have an OpenLayer map with a single vector layer containing many features. I've been using SelectFeature to pop up some info when a feature is clicked, and that's working fine. However, I now want to add some functionality when the feature is double-clicked - specifically I just want to zoom in, and not have the popup appear. I gave this a shot: http://blog.slashpoundbang.com/post/2342888467/openlayers- zoom-when-a-feature-in-a-vector-layer-is But it removed my single-click functionality. Essentially this is the same problem as SelectFeature with Double-Click Any ideas? UPDATE: So I'm trying a workaround using SelectFeature to handle the single click and registering an event listener on the vector layer to handle the double click. All I want to do is zoom in on the double clicked feature. This is what I have:                   VEClayer.events.register('dblclick', VEClayer, function(e){                         map.zoomToExtent(e.features.geometry.getBounds(), closest=true);                             }                             );      However now I have another problem, and I'm not sure if it's being caused by having both SelectFeature and the listener active or what. The new issue is that try as I might I can't access any of the event attributes. So when I try to get e.features.geometry.getBounds(), Firebug gives me the error "e.features is undefined". I've tried to get this multiple other ways (e.xy, e.longlat, getMousePosition) and it's always the same result.


title: Why can't I export XY data to a shapefile?

text: I have added a table from Access to ArcMap, and then used 'Display XY data' to plot points. When I then try to export the Events to a shapefile I get an error message telling me: "Error exporting data. Create output feature class failed" Followed by: "Warning: Some records may not have been exported correctly." and then: "There was an error exporting the layer." Can anyone tell me why this is happening, and how I can export my data? I was able to export an almost identical file this morning with no problems. I'm saving straight to my C drive, so no unusual characters in the files path. All the fields are either text or Long types, and none have long names/unusual characters in their names. I've tried restarting ArcMap and my machine, and exporting to several different locations. I'm using version 10, with Access 2007, and there are 7332 points in the table. Thanks, Flo


title: Change Icon of Single Vector Feature in Layer

text: Currently I am trying to change the icon of a particular feature of a vector layer that the user is focusing on. I add each feature to the map like so:               var point = new OpenLayers.Geometry.Point(pt.lon, pt.lat);     var markerStyle = OpenLayers.Util.extend(OpenLayers.Feature.Vector.style['default'],  {         externalGraphic: iconURL     });     var marker = new OpenLayers.Feature.Vector(point, attributes, markerStyle);      Later I do the following to update the feature's icon:               var marker = this.findSelectedMarker();     if (marker) {         marker.style.externalGraphic = newIconUrl;         this.layer.redraw();                  }      But when the layer redraws, all features in my layer use the `newIconUrl`, not simply the selected marker that I am trying to update. How can I change the icon of the one selected feature of my layer? Thanks.


title: How to export/save a layer composition to GeoTiff?

text: I want to export/save a QGIS composition of raster/vector layers to a GEOTIFF or ECW raster with high resolution (20 000 * 20 000 Pixel).


title: Hide labels for features without data

text: I have a FeatureLayer that uses a ClassBreaksRenderer. The renderer loads data from a collection like so:               var renderer = new esri.ClassBreaksRenderer(                 null,                 function (graphic) {                     return data[graphic.attributes.ID];                 });      Using a LabelLayer, is there a way to only render labels when the feature has a value (as determined by the ClassBreaksRenderer)? I am able to "hide" the labels in this instance, by setting the font size to 0, but the DOM elements are still created. On a layer with thousands of features where only a sub-set are displayed, that kills performance.


title: Are there any Hosting companies, that offer Geoserver?

text: I am working with a non-profit, and we are planning to have a webmap for them. Most of the hosting providers just offer a simple LAMP setup on their Linux option, and if you want to install any software on your own, you need to go for a dedicated server (either real or virtual) which is very expensive for the organisation I am working with. I have seen that there is a hosting provider for Mapserver based maps. Is there any similar provider for Geoserver based maps?


title: GeoTiff size - what to do?

text: I have a whole series of GeoTiff tiles, each one about 250mb in size, in total just over 90gb. I want to be able to display these as a wms into QGIS to 20+ GIS users and in the end make them into a tile cache. My problem is that I don't have enough server space to stick the GeoTiff tiles on to in their current state, so they are living on an external hard-drive. My server contains GeoServer and I've also got access to MapServer on the same server but I'm less fluent its configuration than I am with GeoServer. Is there a way that I could reduce the size per GeoTiff to make them more manageable or should I convert them into another file tile, png/jpeg2000 for example? Has anyone else any experience of such large GeoTiff sizes and what did you do with them? Oh and I'd like to do this entirely with OpenSource (if possible). Med


title: QGIS - "Bad Allocation" Exception - on identify

text: I have QGIS (v1.8) set up to identify in Top down mode. When I use the identify tool on the map, I frequently get a Bad Allocation exception message and I was wondering what gives rise to this error. It seems to occur, after all layers have been scanned for intersects with the click location (i.e. the status displays as 'identifying done') - it seems to appear when the popup window is being generated. Any thoughts on what I can try to do to resolve this? THanks


title: IS it possible to run GRASS commands in the QGIS graphical modeler?

text: I have a geoalgorithm set up in the graphical modeler and everything runs except for the GRASS commands. It is really frustrating because when I execute these commands on their own they work just fine but as soon as I run them as part of a larger process they simply do not produce the desired results. I am using the mapcalc command to extract raster values below 10, for example.


title: Web-based map authoring with ArcGIS Server?

text: I'm embarking on a project in which users will author, save, and share their own maps over the web. We will provide them with a large number of feature classes, but users will effectively author their own maps, map symbologies, etc. Furthermore, they will create and edit their own feature classes, which they can both map and share with other users. The model for AGS map services seems to be: author a map in ArcMap, save an MXD/MSD, publish. I'm struggling to understand how this can help us build a dynamic web mapping platform as described above. Can anyone offer some tips on how to go about it?


title: Slope, aspect and beam in R.Sun (GRASS)

text: I'm starting to mess about with the solar radiation packages of GRASS and ArcGIS and have become rather stuck on the GRASS side of matters. I have a DSM area, and that is feeding into GRASS well enough, however when processing on the default values it falls over, not outputting anything. I'm using the command...               g.region rast=elevation.dem     r.sun -s elev=elevation.dem day=180     r.colors col=grey     d.his i_map=beam_map h_map=elevation.dem      In the R.Sun man page, the following command is recommended:               g.region rast=elevation.dem     r.sun -s elev=elevation.dem slop=slope asp=aspect beam=beam_map day=180     r.colors beam_map col=grey     d.his i_map=beam_map h_map=elevation.dem      However here, where am I supposed to get the slope, aspect and beam from?


title: How to solve loading performance issues when handling multiple KMLs?

text: I am using over 10 KML in Google Maps. When I increase the number of KMLs they load too late. I use checkboxes to hide/show KML layers. With up to 5 KML layers, it worked fine. When I increase the size of the KMLs it's not working normally. When I check "show", it shows but "hide" does not work ... How can I solve this?


title: What to call these types of grids?

text: In a GIS application I'm developing there's an option to show two types of grids. First is a standard lat/lon grid. Second is a crosshair with marks and distances in four directions away from the center. ![attached image](http://i.stack.imgur.com/QpOib.png). Since English isn't my native language I wonder how to better call these two types in menus, etc. to distinguish between them?


title: Overlay analysis of topographic data

text: I have a topographic map with Buildings, Roads, Hydro and land use as the prominent layers in polygon format. Also i have newly collected Building polygons of the same area in shape file format. Now i need to have an overlay analysis of Buildings in these two resources. I can rasterize and geocode the map to orient with Digital data. Also these two can be overlaid and compared manually. But i need to do this analysis automatically with possibly using a shorter method. One possibility is to vectorize Buildings in the map and perform vector overlay. However I am looking at a more efficient method possibly filtering Buildings in raster and comparing it with Buildings converted to rater in shape file. Any idea on this is appreciated.


title: Can ESRI's File Geodatabase API be Redistributed?

text: ESRI's file geodatabase API can now be used with GDAL/OGR - http://www.gdal.org/ogr/drv_filegdb.html however each individual seems to have to login to ESRI and download the API files individually. For example this is referenced in a blog post here and in the GDAL build notes. The license agreement before downloading the API shown by ESRI is at: http://resources.arcgis.com/node/agreement/3193 Clause 47 states: >   1. Licensee may develop and distribute software or web applications that > use the Esri File Geodatabase API to Licensee's end users. > The best discussion about the license is at Paul Ramsey's blog post from January 2011. This finishes with a final update that: > Update 4: Received email from ESRI confirming that the final license will be > reviewed to ensure there are no ambiguities and that it reflects their > intent that the API be usable by any application in any application category > and the derived product freely redistributable and royalty free. So to the > extent that the current license has any ambiguity it shouldn't be considered > a red flag that the final one will. So has or when will the ESRI API be available as part of GDAL downloads or installations, and when can we distribute application which include the API?